Fractional-Order Ultrasonic Sensing for Monitoring Microstructural Evolution in Cementitious Materials

TL;DR

A new ultrasonic sensing method using fractional-order models tracks cement hardening, offering better accuracy than traditional methods.

Contribution

A novel fractional-order ultrasonic framework is introduced to monitor cement hydration with physics-based accuracy.

Findings

The fractional-order model outperforms classical and empirical models in tracking hydration.

The parameter β reflects the viscoelastic transition during cement hardening.

The framework enables real-time, physics-informed monitoring of cementitious materials.

Abstract

Monitoring the early-age evolution of cementitious materials is essential for ensuring the quality and reliability of concrete structures. However, most ultrasonic approaches rely on empirical correlations and lack a physics-based mechanism to describe the continuous viscoelastic transition during hydration. This study proposes a fractional-order ultrasonic sensing framework that couples a fractional Zener viscoelastic model with ultrasonic attenuation theory to quantitatively link microstructural evolution and measured acoustic responses. A custom ultrasonic measurement system was developed to capture real-time attenuation during hydration under different water-cement ratios. Results show that the fractional-order model achieves higher accuracy and robustness than classical integer-order and empirical models. The fractional parameter β serves as a physically interpretable indicator that reflects the transition from viscous-dominated to elastic-dominated behavior and aligns with known hydration development. The proposed framework provides a compact, physics-informed sensing strategy for early-age characterization of cementitious materials and offers potential for intelligent construction and high-end structural monitoring.

Figures

Full text

1. Introduction

The rapid development of high-end infrastructure and equipment—such as shield tunnels, nuclear engineering structures, offshore wind foundations, and large-scale prefabricated components—has created an urgent demand for advanced sensing technologies capable of evaluating material quality and structural performance throughout their entire life cycle [1,2,3]. Cementitious materials, which constitute the core load-bearing components of these systems, undergo complex microstructural evolution during early-age hydration. This evolution directly governs their rheological behavior, stiffness development, defect formation, and ultimately the reliability and service safety of high-end structural equipment. Therefore, real-time and physically interpretable monitoring of the early-age evolution of cementitious materials has become a key enabling technology for intelligent construction and structural health management [4,5,6].

Ultrasonic sensing offers a unique advantage in this context due to its non-destructive nature, high sensitivity to microstructural transitions, and compatibility with embedded and automated monitoring systems [7,8,9,10]. Traditional ultrasonic methods generally rely on empirical correlations between wave velocity, amplitude, or attenuation and the hydration state [1,2,6]. However, such empirical sensing approaches often lack a rigorous physical foundation and provide limited interpretability, especially for the transition from viscous-dominated behavior to viscoelastic–solid behavior during structural formation. The need for model-based, physics-informed ultrasonic sensing has become increasingly evident in high-end engineering applications [3,11].

Recent advances in fractional calculus have demonstrated that fractional-order operators provide a compact yet powerful framework for describing non-local and memory-dependent behavior in complex materials, especially in viscoelasticity and wave propagation [12,13]. In contrast to classical Kelvin–Voigt or Maxwell models, which approximate the relaxation spectrum using a finite number of rheological elements, general fractional viscoelastic models can reproduce broad-band relaxation and power-law-type creep or relaxation with significantly fewer parameters, while maintaining strong physical interpretability [14,15,16,17]. For concrete and other cement-based materials, several recent studies have introduced fractional viscoelastic formulations to capture long-term creep, temperature-dependent deformation and nonlinear damage evolution, showing that fractional models not only achieve higher fitting accuracy but also better reflect the continuous transition of mechanical behavior under varying stress and environmental conditions [16,17]. These works collectively indicate that fractional viscoelasticity is particularly well suited to describing the multi-scale, microstructure-driven rheology of cementitious materials [12,13,14,15].

In parallel, ultrasonic techniques have been extensively employed to monitor the hydration process and early-age property development of cement pastes and concretes, typically by tracking the evolution of ultrasonic velocity, attenuation and backscattered signals [6,7,11,18,19]. Recent developments include omnidirectional or embedded PZT-based transducers and active sensing schemes for identifying key hydration stages and predicting strength development, further highlighting the potential of ultrasound for real-time, non-destructive monitoring in practical structures [19,20]. However, most existing ultrasonic monitoring approaches rely on empirical correlations or simplified linear viscoelastic or poroelastic assumptions for the cement paste matrix, while several studies have applied fractional calculus to model viscoelastic or creep behavior in cementitious materials, the explicit integration of a rigorously derived fractional viscoelastic constitutive law into ultrasonic attenuation and dispersion modeling remains limited [14,21,22,23]. Even in studies where guided waves and fractional viscoelastic parameters are jointly considered—for instance, in layered cement-grouted anchorage systems—the focus remains on parameter identification instead of hydration-induced microstructural evolution [23,24]. Consequently, there is still a clear gap between (i) the rapidly maturing fractional viscoelastic modeling of cementitious materials and (ii) physics-informed ultrasonic sensing models that can exploit fractional parameters as real-time indicators of early-age microstructural evolution, particularly in the context of high-end structural and equipment applications.

To address these limitations, this study introduces a fractional-order ultrasonic sensing framework that couples a fractional Zener viscoelastic model with an ultrasonic attenuation model to quantitatively link hydration-induced microstructural evolution with measurable acoustic attenuation characteristics. A fractional parameter $[eqn]$ is introduced to describe the continuous transition from viscous to elastic dominance, providing a physically interpretable sensing indicator. An embedded ultrasonic monitoring system is developed using PZT transducers and a custom excitation–acquisition unit, enabling continuous measurement of attenuation during the early-age period.

In recent years, structural health monitoring (SHM) research for cementitious and reinforced concrete materials has expanded significantly, driven by advances in sensing technologies, non-destructive evaluation (NDE), and data-driven analysis. At the material level, studies on novel cementitious systems—such as self-leveling sand concrete—have investigated the relationships between mixture design, mechanical performance, durability, and non-destructive indicators, demonstrating how material composition and processing influence measurable signals relevant to structural assessment [25]. These efforts highlight the growing role of NDE techniques in linking material properties to performance-oriented monitoring.

At the structural and component levels, machine learning (ML) and deep learning (DL) approaches have been increasingly applied to SHM of concrete and reinforced concrete systems. Recent studies and reviews report that data-driven methods, including neural networks, support vector machines, and ensemble learning algorithms, can significantly enhance damage detection, strain prediction, and condition assessment accuracy in beams, slabs, and bridge structures compared with conventional approaches [26]. In parallel, the integration of piezoelectric sensor networks with deep learning has enabled real-time monitoring of concrete strength and structural state, illustrating the synergy between physics-based sensing and AI-driven interpretation [27]. In addition, research on self-sensing cementitious composites further demonstrates the potential of intrinsically sensing materials for SHM applications [28]. Collectively, these studies establish a broad SHM landscape spanning advanced materials, sensing modalities, and intelligent data analytics, within which the present work positions a physics-informed, fractional-order ultrasonic framework focused on early-age cementitious materials.

Although the present study primarily focuses on the early-age behavior of cement paste at the material scale, this scope is intrinsically aligned with the broader framework of structural health monitoring (SHM). In SHM practice, reliable condition assessment and damage evolution analysis critically depend on the accurate establishment of initial material states and baseline mechanical properties. Early-age hydration governs the formation of load-bearing microstructures, residual stress development, and initial stiffness, which together determine the long-term durability, damage susceptibility, and service performance of concrete structures. Consequently, material-scale sensing during early-age stages provides indispensable baseline information for subsequent component- and structure-level SHM.

From a multi-scale SHM perspective, the proposed fractional-order ultrasonic sensing framework can be naturally extended from laboratory-scale cement paste to structural components and full-scale concrete structures. The physically interpretable fractional parameter $[eqn]$ characterizes the intrinsic viscoelastic state of cementitious materials and can serve as a state variable within hierarchical SHM systems, linking microstructural evolution, meso-scale material properties, and macro-scale structural response. When embedded or distributed ultrasonic transducers are deployed in real structures, the evolution of $[eqn]$ can be continuously tracked and integrated into digital twin or lifecycle monitoring frameworks, enabling early detection of abnormal material evolution, construction defects, or long-term degradation trends.

The major contributions of this work are summarized as follows:

1. Most existing ultrasonic monitoring approaches for cementitious materials rely on empirical correlations or integer-order viscoelastic/poroelastic assumptions, which provide limited physical interpretability of attenuation mechanisms. The proposed framework fills this gap by embedding a fractional viscoelastic constitutive law directly into the ultrasonic attenuation model, enabling a physics-informed description of material-state evolution;
2. Early-age cement hydration is characterized by a continuous transition from viscous-dominated to elastic-dominated behavior and by distributed relaxation processes across multiple time scales. The fractional Zener model naturally captures this evolving viscoelastic behavior with a compact parameter set, whereas integer-order models with discrete relaxation times are insufficient to represent such continuous transitions;
3. The fractional parameter $[eqn]$ serves as a physically interpretable state variable that quantifies the degree of viscoelastic transition during hydration. Unlike conventional ultrasonic indicators (e.g., wave velocity, amplitude, or empirical attenuation coefficients) or isolated rheological parameters, $[eqn]$ captures the continuous evolution of intrinsic energy dissipation mechanisms in the cement paste;
4. Compared with previous studies that combine ultrasonic measurements with viscoelastic descriptions in a phenomenological manner, the present work advances the state of the art by directly coupling fractional constitutive modeling with ultrasonic attenuation analysis. This integration improves interpretability and provides a clearer pathway toward multi-scale SHM and digital twin applications.

It should be explicitly noted that the experimental investigations presented in this study are conducted on cement paste only. This choice is intentional, as cement paste constitutes the load-bearing matrix of concrete and governs the intrinsic viscoelastic response during early-age hydration. By excluding aggregates, fibers, and reinforcement, the proposed framework isolates the fundamental coupling between microstructural evolution and ultrasonic attenuation, thereby enabling clear physical interpretation of the fractional parameter $[eqn]$ . While this focus inevitably limits the direct generality of the experimental results, it provides a necessary and well-controlled foundation for subsequent extension to more heterogeneous cementitious composites.

The proposed sensing framework provides a foundation for developing next-generation, model-based, embedded ultrasonic systems for quality assessment in high-end structural equipment. It also offers new opportunities for intelligent construction, defect prevention, and digital twins of cement-based infrastructure.

2. Theoretical Model

This section delineates the development of a fractional-order theoretical model to precisely characterize the ultrasonic wave propagation in cementitious materials during the early-age hydration process. The limitations of traditional integer-order models are first analyzed to justify the introduction of fractional calculus. Subsequently, the core framework of the fractional-order Zener constitutive model, based on the Caputo derivative, is elaborated, including its frequency-domain expression and the derivation of the ultrasonic attenuation coefficient. The universality of the proposed model is further verified through asymptotic analysis under classical limits.

2.1. Limitations of Traditional Models and the Introduction of Fractional-Order Theory

In the domain of viscoelastic energy monitoring for cement-based materials, traditional integer-order rheological models (e.g., the Bingham model $[eqn]$ ) are limited because they cannot capture hereditary effects arising from structural evolution and time-dependent properties. Since integer-order derivatives only reflect instantaneous states, they fail to describe the long-range temporal correlations established by the continuously developing hydration network. Therefore, this study introduces fractional calculus theory and constructs a Caputo-type fractional-order constitutive model for viscoelastic characterization. Its derivative is defined as:

$$[eqn]$$

The Caputo fractional derivative incorporates a kernel $[eqn]$ with inherent memory properties, enabling accurate depiction of the evolving internal structure of cementitious materials—such as the relaxation process after micro-crack closure or cutting damage. Compared with the Riemann–Liouville definition, a major advantage of the Caputo derivative is that its initial conditions are expressed in the same form as classical integer-order models, avoiding ambiguity in specifying initial parameters. The fractional order $[eqn]$ introduces dynamic flexibility to the model, supporting its use as a diagnostic parameter for monitoring material evolution. Physically, the fractional kernel describes the cumulative damage and relaxation behavior associated with microstructural evolution.

2.2. Fractional-Order Zener Constitutive Model and Its Frequency-Domain Representation

To establish a quantifiable rheological constitutive relationship, the Caputo fractional derivative must be embedded within a physical framework capable of simultaneously characterizing the elastic and viscous behaviors of the material. Accordingly, this study adopts the classical Zener model as the foundational structure for fractional extension. The Zener model, composed of a spring and dashpot in series, can already describe stress relaxation—a key viscoelastic phenomenon of cement-based materials—in its integer-order form. By replacing the integer-order time derivative with the Caputo fractional derivative, the model naturally incorporates power-law memory effects into the stress–strain relationship, thereby yielding an enhanced constitutive equation that captures instantaneous elastic response, long-term relaxation behavior, and the inherent dependence on loading history.

$$[eqn]$$

Here, $[eqn]$ is the fractional relaxation time, and $[eqn]$ and $[eqn]$ denote the instantaneous and long-term elastic moduli, respectively. The parameter $[eqn]$ serves as a microstructural diagnostic index: as $[eqn]$ the system approaches a viscous-like fluid, whereas as $[eqn]$ the system tends toward an elastic solid.

To facilitate the analysis of the model in the frequency domain, a Fourier transform is applied to the fractional constitutive Equation (2). Using the frequency-domain property of fractional derivatives, $[eqn]$ , one obtains:

$$[eqn]$$

The complex modulus $[eqn]$ is defined as the ratio of stress to strain in the frequency domain, i.e., $[eqn]$ . Dividing both sides of Equation (3) by $[eqn]$ ( $[eqn]$ ) and substituting the definition of the complex modulus, one directly obtains:

$$[eqn]$$

After rearrangement, the explicit expression of the complex modulus is obtained as:

$$[eqn]$$

The real part (storage modulus $[eqn]$ ) and the imaginary part (loss modulus $[eqn]$ ) of the complex modulus, respectively, characterize the elastic and viscous responses of the material:

$$[eqn]$$

The frequency-domain representation of the complex modulus $[eqn]$ depicts the viscoelastic response of cement paste under acoustic excitation. Figure 1 illustrates the variation in the normalized complex modulus magnitude with frequency for different fractional orders $[eqn]$ . As $[eqn]$ increases from 0.2 to 1.0, the overall amplitude of the complex modulus gradually stabilizes at a higher level, indicating that the fractional derivative order governs the transition of the material from a viscous-dominated state to an elastic solid-like state. When $[eqn]$ is relatively small, the system exhibits pronounced frequency dependence, reflecting stronger memory effects and a dominant viscous contribution. As $[eqn]$ approaches 1, the modulus becomes nearly frequency-independent, and the system tends toward an ideal elastic solid. This result confirms that the fractional Zener constitutive model formulated in Equation (5) can effectively capture the frequency-dependent viscoelastic behavior of cementitious materials.

To further reveal how such viscoelastic characteristics affect ultrasonic wave propagation, it is necessary to incorporate the fractional constitutive relation into the wave equation, thereby establishing the linkage between acoustic waves and material parameters.

2.3. Derivation of the Ultrasonic Wave Equation and Attenuation Coefficient

This subsection establishes the fractional viscoelastic constitutive description of early-age cement paste. From a physical perspective, the material response during hydration continuously evolves from viscous-dominated to elastic-dominated behavior, accompanied by broad relaxation spectra. The fractional Zener model is adopted here to capture this continuous transition in a compact and physically interpretable manner, forming the basis for subsequent wave propagation and attenuation analysis.

The following assumptions are adopted throughout the theoretical formulation:

• (i)The cement paste is treated as an effective homogeneous and isotropic viscoelastic medium;
• (ii)Wave propagation is approximated as one-dimensional plane-wave motion;
• (iii)Small-strain and linear viscoelastic behavior are assumed;
• (iv)Frequency range corresponds to low-frequency ultrasonic excitation relevant to early-age monitoring.

The overall modeling framework and the associated transformation from the fractional constitutive relation to the ultrasonic attenuation coefficient are schematically summarized in Figure 2.

To establish a direct connection between the fractional constitutive model and ultrasonic testing, it is essential to define the propagation scenario of acoustic waves in cement paste as a viscoelastic medium. As shown in Figure 3, the temporal evolution of the flocculated network during resting—characterized by the fractional order ( $[eqn]$ )—is a key factor influencing ultrasonic propagation. A plane wave generated by the transmitting transducer (T) undergoes attenuation while traveling through this evolving medium because its energy is scattered by the flocculated structures. This leads to a reduced signal amplitude at the receiving transducer (R), and such attenuation can be quantified by the attenuation coefficient ( $[eqn]$ ). The following wave equation is formulated and solved to theoretically describe this physical process.

Assuming the material is unbounded and homogeneous, boundary effects are neglected. Although practical cement paste exhibits particle-level heterogeneity, it can be regarded as macroscopically homogeneous. In such an unbounded homogeneous viscoelastic medium, consider a plane ultrasonic wave propagating along the (x)-direction. Its dynamics are governed by the classical wave equation, whose frequency-domain form is given by:

$$[eqn]$$

Here, $[eqn]$ is the density of the paste, and $[eqn]$ is the displacement field in the frequency domain. The wave equation simultaneously involves the stress $[eqn]$ and the displacement $[eqn]$ , which must be coupled through the constitutive relation of the medium. The infinitesimal strain $[eqn]$ induced by the acoustic wave is defined by the displacement gradient as:

$$[eqn]$$

Combining the constitutive relation $[eqn]$ with the strain definition (8), the stress expression can be obtained as:

$$[eqn]$$

Taking the spatial derivative of Equation (9) and substituting it into the wave Equation (7) yields:

$$[eqn]$$

Combining Equations (7) and (10) and eliminating $[eqn]$ , the expression for the complex wavenumber (k) is obtained as:

$$[eqn]$$

Therefore,

$$[eqn]$$

Since $[eqn]$ is complex, $[eqn]$ is also complex, and its imaginary part corresponds to the ultrasonic attenuation coefficient $[eqn]$ , i.e., $[eqn]$ .

To obtain an explicit solution, the numerator and denominator in Equation (5) are rewritten in polar form. Let the denominator be:

$$[eqn]$$

where

$$[eqn]$$ $$[eqn]$$

Let the numerator be:

$$[eqn]$$

where,

$$[eqn]$$ $$[eqn]$$

Substituting the polar forms into Equation (5), the complex modulus can be simplified as:

$$[eqn]$$

Substituting Equation (19) into Equation (12) and using the square-root operation rule for complex numbers:

$$[eqn]$$

Therefore, the attenuation coefficient $[eqn]$ is the imaginary part of the complex wavenumber $[eqn]$ :

$$[eqn]$$

Equation (21) contains the phase difference. To obtain a more intuitive physical interpretation, it is necessary to simplify it. According to the trigonometric identity $[eqn]$ , an explicit result can be obtained via the sine difference formula. Let $[eqn]$ , $[eqn]$ , and $[eqn]$ , then $[eqn]$ and $[eqn]$ . Using the trigonometric formulas for the sum and half of angles, Equation (21) can be shown to be equivalent to the following form:

$$[eqn]$$

where $[eqn]$ is the acoustic velocity. This equation quantitatively characterizes the coupling relationship among acoustic attenuation, frequency $[eqn]$ , and the fractional order $[eqn]$ .

A smaller $[eqn]$ leads to pronounced velocity dispersion, while as $[eqn]$ , the system approaches the elastic limit and the acoustic velocity becomes nearly independent of frequency.

In the fractional Zener model, the viscoelastic properties of the medium influence not only ultrasonic attenuation but also velocity dispersion. Figure 4a presents the variation in the normalized ultrasonic attenuation coefficient $[eqn]$ with normalized frequency $[eqn]$ for different fractional orders $[eqn]$ . When $[eqn]$ is small, the attenuation coefficient increases significantly with frequency, indicating strong energy scattering within the medium. As $[eqn]$ increases, the growth rate of attenuation is markedly reduced, and the medium gradually transitions toward a low-loss elastic state. This trend agrees with the theoretical prediction of Equation (22), demonstrating that the fractional parameter $[eqn]$ plays a key role in governing band-limited attenuation behavior.

Figure 4b shows the variation in the corresponding normalized phase velocity $[eqn]$ with frequency. The results indicate that the phase velocity is relatively low in the low-frequency region and gradually increases with frequency, eventually reaching a plateau at high frequencies. When $[eqn]$ is small, the medium exhibits pronounced dispersion; as $[eqn]$ increases, the dispersion progressively weakens, and in the limit $[eqn]$ the medium approaches an ideal elastic solid, for which the wave velocity is nearly independent of frequency. This phenomenon confirms that the fractional Zener model can accurately describe the velocity dispersion induced by band-limited attenuation, providing a theoretical basis for interpreting subsequent experimental results and for parameter inversion.

2.4. Model Degeneration Verification

Based on the fractional constitutive relation, this subsection derives the ultrasonic wave attenuation in a viscoelastic cementitious medium. The key idea is that material viscoelasticity introduces frequency-dependent energy dissipation, which manifests as a complex modulus and, consequently, a complex wavenumber. The imaginary part of the wavenumber directly governs the measurable ultrasonic attenuation coefficient.

Under the assumptions stated in Section 2.3, the ultrasonic wave propagation and attenuation are derived as follows.

To verify the correctness and general applicability of the proposed model, the model is reduced to two classical limiting cases: a viscous fluid ( $[eqn]$ ) and an elastic solid ( $[eqn]$ ). This demonstrates that the model is consistent with classical theory under limiting conditions.

In the viscous-fluid limit ( $[eqn]$ ), the equilibrium modulus $[eqn]$ of the material is much smaller than the instantaneous modulus $[eqn]$ , and the mechanical response is dominated by viscosity. To further degenerate the model into a purely viscous fluid (Newtonian fluid), one additional step is needed: letting $[eqn]$ . Through small-parameter asymptotic analysis, one can show that the attenuation coefficient $[eqn]$ is proportional to $[eqn]$ , which corresponds to the classical acoustic attenuation law of Newtonian fluids.

Considering Equation (22), the asymptotic behavior as $[eqn]$ is analyzed as follows: $[eqn]$ (retaining the first-order term of the Taylor expansion), $[eqn]$ (small-angle approximation), $[eqn]$ (small-angle approximation), Substituting these into the numerator gives: $[eqn]$ retaining terms up to first order in $[eqn]$ :

$$[eqn]$$

Substituting into the square-root term in the denominator:

$$[eqn]$$

Combining the constant term and the $[eqn]$ -dependent term yields:

$$[eqn]$$

Combining the numerator and denominator yields:

$$[eqn]$$

When $[eqn]$ , the $[eqn]$ -dependent terms in the denominator can be neglected:

$$[eqn]$$

In the viscous-fluid limit, the material response is predominantly viscous and elasticity can be neglected, which requires $[eqn]$ (the long-term modulus approaches zero). Meanwhile, the fractional parameter $[eqn]$ is linked to the viscosity η through the relaxation time τ $[eqn]$ . As $[eqn]$ , $[eqn]$ must remain finite, implying that $[eqn]$ tends to a constant. A more rigorous derivation shows that, under the double limit $[eqn]$ and $[eqn]$ , we have:

$$[eqn]$$

Here, $[eqn]$ denotes the acoustic velocity; in a purely viscous fluid, $[eqn]$ , where $[eqn]$ is the effective modulus.

When $[eqn]$ , the material behaves as an elastic solid. We directly substitute $[eqn]$ and consider the high-frequency approximation to demonstrate that the attenuation coefficient $[eqn]$ is proportional to $[eqn]$ , which is consistent with the acoustic attenuation characteristics of elastic solids.

Substituting the exact value $[eqn]$ into Equation (22) and simplifying yields:

$$[eqn]$$

In the elastic-solid limit ( $[eqn]$ ), the material response is dominated by elasticity and viscous effects become negligible. This corresponds to the high-frequency limit ( $[eqn]$ ) or ( $[eqn]$ ). When $[eqn]$ , we have $[eqn]$ . Substituting this approximation yields:

$$[eqn]$$

Since $[eqn]$ , substituting this relation yields the further simplified form:

$$[eqn]$$

To remain consistent with the classical form, we note that in an elastic solid the acoustic velocity is $[eqn]$ , and $[eqn]$ therefore:

$$[eqn]$$

Therefore, the above expression can be written as:

$$[eqn]$$

Here, $[eqn]$ adopts the long-term acoustic velocity $[eqn]$ .

The proposed fractional model correctly degenerates to classical theories under both limiting conditions. When $[eqn]$ , small-parameter expansion and double-limit analysis demonstrate that $[eqn]$ is proportional to $[eqn]$ , consistent with the acoustic attenuation law of Newtonian fluids. When $[eqn]$ , high-frequency approximation shows that $[eqn]$ is proportional to $[eqn]$ , matching the attenuation behavior of Hookean elastic solids.

In summary, the limiting cases $[eqn]$ and $[eqn]$ recover the classical viscous-dominated and elastic-dominated responses, respectively, confirming the physical consistency of the proposed formulation. These results demonstrate that the fractional model provides a continuous and unified description bridging the two extremes, while the detailed derivations are retained above for completeness.

3. Materials and Methods

This study employs fractional rheological theory to quantitatively characterize the microstructural evolution of cement-based materials during the setting process using ultrasonic monitoring. The experimental section systematically details the material preparation, ultrasonic testing system, experimental procedures, and data processing methods, ensuring reproducibility and precise control throughout the testing process.

3.1. Material Preparation and Specimen Design

The experiments employed P·O 42.5 ordinary Portland cement (Jiuchang Building Materials Business Department, Shandong, China), whose chemical composition is shown in Table 1 and complies with the GB/T 176-2017 standard [29]. A polycarboxylate-based superplasticizer (Dongrun Baisheng New Materials Co., Ltd., Sichuan, China) (PCE) was added at 0.2% of the cement mass. In addition to the chemical composition summarized in Table 1, the cement paste density and ultrasonic wave speed were measured at the initial stage to provide baseline material properties for subsequent analysis. The initial bulk density of the fresh cement paste was approximately 1900 kg/m^3^, and the corresponding longitudinal ultrasonic wave velocity was on the order of (1.2–1.8) × 10^3^ m/s, which is consistent with typical values reported for early-age cement pastes. These properties primarily affect the absolute magnitude of wave propagation parameters but do not alter the trend-based identification of the fractional parameter $[eqn]$ during hydration.

To investigate the influence of the water–cement ratio on rheological behavior, four levels were selected: 0.40, 0.45, 0.50, and 0.55. For each mixture, 500 mL of cement paste was prepared. The specific mix proportions were as follows: cement was fixed at 1000 g, and the corresponding water masses were 400 g, 450 g, 500 g, and 550 g, respectively. The PCE dosage was uniformly maintained at 2.0 g.

Mixing was carried out in a controlled-temperature environment (25 ± 1 °C). A planetary mixer was used to dry-mix cement and superplasticizer at 140 r/min for 1 min, after which water was added and wet-mixing was performed at 285 r/min for 3 min to ensure uniform dispersion without agglomeration. The fresh paste was poured into a custom acrylic container (150 mm × 50 mm × 50 mm), allowed to rest for 5 min to eliminate entrapped air bubbles, and then subjected to the monitoring stage.

3.2. Design of the Ultrasonic Monitoring System

The ultrasonic monitoring system was constructed based on a through-transmission configuration. The core transmitting and receiving components were a pair of 25 kHz piezoelectric ceramic transducers (model: DYW-H25-200CTY) (Dayu Electronics Technology Co., Ltd., Fujian, China), with a center frequency of 25 kHz ± 1 kHz. The transducers exhibit a broadband frequency response covering the excitation frequency range used in this study, ensuring stable signal transmission and reception during early-age monitoring. The choice of low-frequency PZT elements minimizes scattering effects and enhances penetration in highly attenuative fresh cement paste. The transducers were mounted on opposite sides of an acrylic container, with a fixed sensing distance of 150 mm.

The excitation signal was generated by an STM32F407 microcontroller (STMicroelectronics, Geneva, Switzerland) and shaped using a Hanning-windowed 5-cycle sinusoidal burst. After amplification by a high-voltage driver, a 50-Vpp excitation pulse was delivered to the transmitting transducer. The received signal first passed through a 40 dB preamplifier and was then sampled by a 16-bit ADC module at a sampling rate of 1 MHz.

To compensate for temperature-induced variations in acoustic velocity, the system incorporated an NTC thermistor connected via an SPI interface to monitor the ambient temperature (accuracy ±0.5 °C). The measured temperature was used to correct sound velocity according to the empirical relation $[eqn]$ where $[eqn]$ . Temperature effects were controlled by conducting all experiments in a temperature-stabilized laboratory environment. In addition, temperature was continuously monitored during the test period to ensure that fluctuations remained within a narrow range. The proposed framework relies on baseline-referenced attenuation trends rather than absolute signal amplitudes, which inherently reduces sensitivity to minor temperature-induced variations. The effectiveness of temperature compensation was verified by comparing attenuation trends under similar hydration conditions, confirming consistent evolution of the identified fractional parameter $[eqn]$ .

3.3. Experimental Protocol and Procedure

Each experiment lasted 120 min and began immediately after paste preparation and the subsequent resting period. Ultrasonic signals were collected every 5 min. At the initial stage, the system was calibrated using deionized water as the reference medium. During the experiments, the ambient temperature was maintained at 25 ± 1 °C and the relative humidity at 50 ± 5% RH. The baseline amplitude $[eqn]$ and temperature $[eqn]$ were recorded throughout the process. Figure 5 experimental setup for embedded ultrasonic sensing and the early-age evolution of cement paste observed from 0 to 120 min.

During formal measurements, three pulse signals were acquired in each cycle and averaged to suppress random noise. At the end of each test, all time-domain signals and temperature data were stored for subsequent analysis. To ensure data consistency, each water–cement ratio level was tested in triplicate, and the final results were obtained by averaging after removing outliers.

3.4. Data Processing and Parameter Inversion

Ultrasonic signal processing was performed using a time-domain attenuation approach. The raw signals were first corrected for DC offset to eliminate hardware-induced baseline drift. Subsequently, a Hilbert transform was applied to extract the amplitude envelope of the signal, representing the acoustic energy decay. To suppress noise interference, the envelope data were smoothed using a moving-average filter.

Based on the physical model of amplitude decay in a viscoelastic medium ( $[eqn]$ ), the natural logarithm of the envelope amplitude within the steady-attenuation region was taken and linearly fitted against propagation distance. To extract the ultrasonic attenuation coefficient, the received signal was first processed using an envelope detection method. The attenuation region was selected from the main transmitted wave packet after excluding the initial excitation transient and late-time reflections from specimen boundaries. Specifically, a time window corresponding to the stable propagation region with high signal-to-noise ratio was identified, and the exponential decay of the envelope within this region was fitted to determine the attenuation coefficient. This procedure ensures consistency across repeated measurements and minimizes the influence of noise and boundary effects. The slope of the fitted line corresponds to the decay rate $[eqn]$ (unit: s^−1^), and the attenuation coefficient $[eqn]$ (unit: Np/m) was obtained using $[eqn]$ where $[eqn]$ is the sound velocity.

The density $[eqn]$ of each sample was measured immediately after specimen preparation and used as a constant input parameter. The instantaneous elastic modulus $[eqn]$ was calculated from the propagation theory of elastic waves in solids, using the measured sound velocity $[eqn]$ and density $[eqn]$ .

The long-term elastic modulus $[eqn]$ was estimated using the effective modulus method, whose theoretical basis originates from effective medium theory in cementitious materials. This method simplifies long-term creep and viscoelastic effects into a reduced elastic modulus, generally expressed as $[eqn]$ where $[eqn]$ is the ultimate creep coefficient.

Considering the early-age characteristics of cementitious materials in this study, $[eqn]$ was parameterized using an empirical relation associated with the water–cement ratio: $[eqn]$ which enables rapid estimation and reflects the influence of mix proportion on flow and deformation behavior.

Finally, the coefficient of determination ( $[eqn]$ ) was calculated to evaluate the goodness of linear fits, ensuring the reliability of parameter inversion.

Each experiment was repeated three times under identical conditions to assess repeatability. Outliers were screened by comparing signal quality and attenuation trends; no abnormal measurements inconsistent with physical expectations were observed. For each time point, the mean attenuation coefficient is reported, while the experimental variability is represented by the minimum and maximum values among the three repeated measurements. This approach emphasizes repeatability and trend consistency rather than statistical significance testing, which is appropriate for early-age monitoring studies.

4. Results and Discussion

The core objective of this section is to integrate experimental data with the theoretical model to systematically analyze time-domain signal characteristics, ultrasonic attenuation behavior, variations in the fractional order $[eqn]$ , and multi-model comparisons. These analyses collectively reveal the evolution of the microstructure as it transitions from a viscous fluid to an elastic solid. Cement pastes with different water–cement ratios (Section 3.1) were tested using the ultrasonic monitoring system, and the results were compared with the theoretical predictions of the fractional constitutive model (Equation (5)), providing a quantitative basis for monitoring the rheological properties of cement-based materials.

4.1. Evolution Characteristics of Typical Time-Domain Signals of Cement Paste

As shown in Figure 6, taking the ultrasonic signal of cement paste with a water–cement ratio of 0.40 as an example, the time-domain responses at four key moments—10 min, 60 min, 90 min, and 120 min—were comparatively analyzed. At 10 min (a), the signal exhibits characteristics typical of a viscous-fluid-dominated state. The overall amplitude remains relatively high, but the envelope decays extremely rapidly, with oscillations vanishing within a very short duration and the waveform appearing irregular. This indicates that, during the early stage of setting, the dispersed cement particles are separated by free water, and the pronounced viscous dissipation causes ultrasonic energy to be quickly attenuated during propagation. At this moment, the material possesses very low structural strength, and its mechanical behavior closely resembles that of a Newtonian fluid. In terms of the fractional Zener model, this corresponds to a fractional order $[eqn]$ approaching zero, meaning that the constitutive response is dominated by viscous behavior.

As hydration progresses to 60 min (b) and 90 min (c), the time-domain signals exhibit clear transitional characteristics. Compared with the 10 min response, the most notable change is the increased persistence of oscillations. Beyond 500 μs, multiple cycles of regular, continuous high-frequency oscillations become visible. Although the peak amplitude decreases, the envelope attenuation rate is substantially reduced.

This transformation signifies fundamental microstructural evolution within the paste: hydration products (C–S–H gels) form progressively and interconnect to establish an initial three-dimensional flocculated network. This network enhances the continuity of the medium and strengthens its elastic response. Consequently, the ultrasonic energy-loss mechanism shifts from purely viscous dissipation toward viscoelastic damping, with elastic wave-propagation characteristics becoming increasingly pronounced.

In terms of fractional-order theory, this evolution is reflected by the gradual increase in the fractional order $[eqn]$ toward 1, indicating that the material’s “memory” and elastic contribution are continuously strengthened, and the constitutive behavior transitions from a fractional derivative toward an integer-order form.

By 120 min (d), the time-domain signal has evolved to exhibit propagation characteristics dominated by an elastic solid. The signal amplitude stabilizes at a relatively low level, and the waveform displays highly regular, continuous high-frequency oscillations, with only minimal attenuation throughout the entire observation window. This indicates that a dense and rigid spatial network structure has formed within the paste, where particles are interconnected through strong bonding forces.

In such a medium, ultrasonic energy loss arises primarily from elastic scattering caused by microstructural inhomogeneities (e.g., unhydrated particles, microvoids) rather than viscous friction. This stage corresponds to a further increase in the fractional order $[eqn]$ in the fractional model, reflecting the strengthening of elastic behavior and the near-solid-like mechanical response of the material.

4.2. Verification of the Fractional Model and Analysis of Attenuation Characteristics

To validate the effectiveness of the proposed fractional Zener model in describing ultrasonic attenuation during the setting of cement paste, this section systematically compares experimental measurements with model predictions for water–cement ratios (W/C) of 0.40, 0.45, 0.50, and 0.55. Figure 7 illustrate the evolution of the attenuation coefficient $[eqn]$ over time (0–120 min) for each water–cement ratio. Each figure presents the experimental mean value, the fluctuation range of the measurements, the theoretical prediction curve, and the corresponding evolution of the fractional order $[eqn]$ , which serves as a key parameter in the model.

As shown in the figures, for all water–cement ratios, the experimentally measured attenuation coefficients $[eqn]$ (blue circles) exhibit a variation trend that is highly consistent with the theoretical predictions based on the fractional model (Equation (22), red squares). Both curves increase monotonically with setting time, which corroborates the earlier time-domain signal analysis: as hydration products accumulate and the flocculated network develops, the dominant attenuation mechanism transitions from early-stage viscous dissipation to later-stage elastic scattering, resulting in a continuous increase in the macroscopic attenuation coefficient $[eqn]$ .

To quantitatively assess the prediction accuracy of the model, the relative error (RE) between the theoretical predictions and the experimental mean values was calculated over the entire setting process. The results show that the model performs well across all water–cement ratios, although the prediction accuracy varies:

1. W/C = 0.40 (a): The average relative error is 1.85%, with the maximum error (~7.96%) occurring at the early stage. This indicates that for denser pastes with lower water–cement ratios, the fractional model can capture the attenuation characteristics with very high precision.
2. W/C = 0.45 (b): The average relative error is 3.33%. The strong agreement between the predicted curve and the experimental mean further confirms the robustness of the model.
3. W/C = 0.50 (c): The average relative error is 3.35%. Although the fluctuation range of the experimental data (grey region) is larger—suggesting slightly poorer microstructural uniformity at this water–cement ratio—the model still effectively tracks the overall trend.
4. W/C = 0.55 (d): The average relative error is 7.11%, the highest among all groups. This suggests that at higher water–cement ratios, the greater amount of free water leads to more complex rheological behavior, introducing certain limitations to the model’s predictive capability.

As the key parameter of the proposed model, the fractional order $[eqn]$ provides a unique perspective on the microstructural evolution of cement paste. The temporal variation curves (green dash–dot lines) in all four figures exhibit a consistent trend: $[eqn]$ starts from a relatively low value (0.1–0.3) in the early stage of setting and increases approximately linearly as hydration progresses, reaching the range of 0.33–0.55 at 120 min.

At the beginning of setting, $[eqn]$ represents an ideal viscous-fluid response. In the experiments, the initial $[eqn]$ values are greater than zero but far below one, indicating that although the fresh paste is dominated by viscous behavior, it already exhibits a certain degree of “memory” behavior—distinct from Newtonian fluids—which reflects the presence of an initial flocculated microstructure.

The monotonic increase of $[eqn]$ directly captures the continuous formation and strengthening of the internal three-dimensional network. As hydration products such as C–S–H gel progressively interweave and interconnect, the solid-like (elastic) nature of the material increases, while the fluid-like (viscous) contribution diminishes. The inherent long-memory effect of fractional derivatives naturally describes this history-dependent microstructural reconstruction process.

Ultimately, the trend $[eqn]$ corresponds to the ideal Hookean elastic solid. Since $[eqn]$ does not reach 1 by the end of the 120 min observation window, this indicates that hydration remains incomplete and the microstructure is still evolving, with the material exhibiting characteristics of a typical viscoelastic solid.

As shown in Figure 8a, for the same setting time, a higher water–cement ratio corresponds to a lower value of $[eqn]$ . For example, at 120 min, $[eqn]$ for W/C = 0.40, whereas $[eqn]$ for W/C = 0.55. This quantitative relationship indicates that increasing the water–cement ratio dilutes the concentration of cement particles and slows down the formation rate of the flocculated network. Consequently, the paste exhibits stronger macroscopic fluidity (i.e., lower elasticity), which is fully consistent with traditional rheological understanding.

Figure 8b presents the temporal evolution of the residuals ( $[eqn]$ ) between the model predictions and the experimental measurements. For all water–cement ratios, the mean residuals remain very close to zero (ranging from −0.26 to 0.24), indicating that the model predictions are essentially unbiased. In addition, the residual sequences fluctuate randomly around zero throughout the entire setting period, with no observable systematic trends, and the fluctuation magnitude is limited (root-mean-square error RMS between 0.63 and 1.13). These characteristics collectively demonstrate that the fractional model achieves high predictive accuracy for the attenuation coefficients of cement pastes across different water–cement ratios.

To further examine whether the model exhibits systematic bias, Figure 8c analyzes the relationship between the model residuals and the intrinsic parameter $[eqn]$ . The scatter plots show that the residuals are randomly distributed around the zero line, and the trend lines for all water–cement ratios are nearly horizontal. Correlation analysis indicates that the absolute values of the correlation coefficients (r) are all below 0.3, and the corresponding (p)-values are much greater than the 0.05 significance threshold. This confirms that no statistically significant correlation exists between the residuals and the value of $[eqn]$ . In other words, the prediction deviations of the model are random across different structural states (represented by different $[eqn]$ values), rather than arising from any structural deficiency of the model itself.

Finally, autocorrelation analysis of the residual sequences (Figure 8d) shows that, for all cases, the autocorrelation coefficients decay and stabilize within the 95% confidence interval after lag 2. This indicates that the residuals exhibit no significant autocorrelation beyond the second lag and therefore satisfy the key characteristics of white noise. This result confirms that the fractional model has successfully captured the essential temporal dependencies in the evolution of the attenuation coefficient, with no important systematic time information left unmodeled.

Figure 6, Figure 7 and Figure 8 indicate that the evolution of the fractional parameter $[eqn]$ is generally monotonic but not strictly linear over time. Instead, $[eqn]$ exhibits stage-dependent trends that correspond well to the classical hydration process. During the dormant period, $[eqn]$ increases slowly, reflecting a predominantly viscous response of the fresh cement paste. In the acceleration stage, a more pronounced increase in $[eqn]$ is observed, associated with rapid formation of hydration products and the development of a load-bearing microstructure. As hydration progresses toward rigidity and early hardening, the growth rate of $[eqn]$ $[eqn]$ gradually decreases, indicating a transition toward a more elastic-dominated viscoelastic state.

For mixtures with higher water–cement ratios, the prediction error of ultrasonic attenuation becomes more pronounced. This behavior can be attributed to increased variability associated with a higher free-water content, which enhances scattering, local heterogeneity, and transient coupling effects during early hydration. In such systems, the ultrasonic attenuation is more sensitive to small fluctuations in microstructure and pore water distribution, leading to increased dispersion in the measured signals. Consequently, deviations between predicted and measured attenuation are more noticeable at higher water–cement ratios, particularly at very early ages.

The sensitivity of $[eqn]$ to measurement noise and algorithmic adjustments was examined by analyzing its temporal stability and repeatability across repeated measurements. Although small fluctuations in $[eqn]$ may arise from signal noise or envelope fitting procedures, the identified $[eqn]$ evolution is primarily governed by long-term trends rather than pointwise variations. Since $[eqn]$ is extracted through baseline-referenced attenuation trends, its estimation is inherently less sensitive to random noise compared with instantaneous ultrasonic indicators. The observed smooth evolution of $[eqn]$ across hydration stages suggests that the parameter predominantly reflects underlying physical dynamics rather than numerical artifacts.

In some cases, $[eqn]$ approaches a quasi-stationary zone or plateau at more advanced hydration ages. This behavior does not indicate a loss of sensitivity but rather reflects the stabilization of the viscoelastic microstructure once the primary hydration reactions have largely completed. At this stage, the cement paste exhibits a relatively stable elastic-dominated response, and further changes in intrinsic energy dissipation mechanisms become less pronounced. Therefore, the plateau of $[eqn]$ can be interpreted as a signature of early-age mechanical stabilization rather than a limitation of the sensing framework.

The practical use of $[eqn]$ and $[eqn]$ as hydration indicators is subject to several limitations. First, both parameters rely on baseline-referenced measurements and are therefore most effective during early-age hydration, when material evolution is monotonic and well controlled. Second, $[eqn]$ reflects cumulative attenuation effects and may be influenced by experimental conditions such as coupling and boundary effects, whereas $[eqn]$ is a model-derived parameter whose interpretation depends on the validity of the assumed viscoelastic framework. Consequently, $[eqn]$ and $[eqn]$ should be interpreted as complementary indicators of hydration progression rather than absolute material properties.

The proposed method exhibits limited sensitivity to random measurement noise, as both $[eqn]$ and $[eqn]$ are extracted from envelope-based attenuation trends and evaluated through temporal evolution rather than instantaneous values. Temperature fluctuations primarily affect absolute signal amplitudes and wave speeds; however, their influence on $[eqn]$ is mitigated by baseline referencing and trend-based interpretation. Transducer misalignment and coupling conditions may affect absolute attenuation levels, but their impact on the temporal evolution of $[eqn]$ is expected to be secondary, particularly when consistent sensor geometry and fixed excitation–reception paths are maintained.

While the present study focuses on early-age hydration, the proposed framework is, in principle, capable of detecting non-monotonic phenomena such as micro-cracking or segregation. Such processes are expected to introduce abrupt or irreversible changes in attenuation behavior and in the evolution of $[eqn]$ , in contrast to the smooth and monotonic trends observed during hydration. However, distinguishing these mechanisms from environmental or operational effects would require auxiliary sensing information and extended validation, which is beyond the scope of the current work.

The current formulation and experimental validation are limited to cement paste and therefore represent a matrix-level framework. In fiber-reinforced or highly heterogeneous concrete, additional effects such as multi-scattering, anisotropy, and guided-wave propagation are expected to influence ultrasonic attenuation. While the fractional viscoelastic model remains applicable at the effective-medium level, its extension to such systems would require incorporating appropriate homogenization or multi-scattering models. This represents an important direction for future work.

The fractional-order framework is inherently frequency-dependent and can be naturally extended to other ultrasonic frequencies or multi-frequency excitation schemes. Since the fractional constitutive model describes viscoelastic behavior across a broad frequency range, $[eqn]$ can be identified using different excitation frequencies, provided that the corresponding attenuation data are available. Multi-frequency measurements may further enhance parameter identifiability and robustness by constraining the inversion process and enabling frequency-domain consistency checks.

Taken together, the validation results for all four water–cement ratios demonstrate that the proposed fractional model not only provides high-accuracy predictions of ultrasonic attenuation during the setting of cement paste but also yields residuals that are genuinely random. This verifies the robustness and completeness of the model. Overall, the evolution of $[eqn]$ provides a physically consistent and stage-sensitive indicator of hydration progression, offering a robust descriptor of early-age material state that complements conventional ultrasonic attenuation analysis.

4.3. Comparative Analysis Between the Fractional Model and Other Viscoelastic Theories

After validating the fractional model developed in this study, we compare it with three representative viscoelastic models: the DFL (Decoupled Fractional Laplacian) model, the Fractal Modification model, and the classical Kelvin–Voigt model [30,31]. All of these models aim to describe wave propagation in viscoelastic media, yet they differ substantially in theoretical foundations and application focus.

The proposed fractional model introduces the fractional derivative order $[eqn]$ , which effectively captures the memory effects and gradual microstructural evolution occurring during the setting of cement paste. Its frequency-domain complex-modulus formulation provides a direct and intuitive representation of viscoelastic response.

The DFL model employs a decoupled fractional Laplacian operator and is particularly efficient for handling anisotropic attenuation in highly heterogeneous media. Its wavefield simulations allow clear separation of amplitude decay and phase dispersion.

The Fractal Modification model is derived from fractal theory, using parameters such as $[eqn]$ or $[eqn]$ as scaling variables. It is especially suited for mineral slurries where particle irregularity and clustering are prominent.

In contrast, the Kelvin–Voigt model—an integer-order viscoelastic representation composed of a parallel spring–dashpot system—is simple and intuitive but cannot adequately capture frequency-independent attenuation or the dynamic microstructural evolution characteristic of cement paste during hydration.

We performed a quantitative comparison between the proposed fractional model and three representative models—the DFL method, the Fractal Modification model, and the Kelvin–Voigt model—under different water–cement ratios (0.40, 0.45, 0.50, 0.55). The results are shown in Figure 9, and the statistical summary of average relative errors for each model is provided in the accompanying Table 2.

All models are able to capture the increasing trend of the attenuation coefficient over time, consistent with the physical evolution of the microstructure from a viscous fluid to a viscoelastic solid during cement hydration. However, the prediction accuracy varies markedly among the models. The fractional model proposed in this study exhibits the best performance across all water–cement ratios, with average relative errors of 1.85% (W/C = 0.40), 3.33% (W/C = 0.45), 3.54% (W/C = 0.50), and 7.11% (W/C = 0.55). It accurately reflects the temporal evolution of the attenuation coefficient.

In contrast, the prediction errors of the DFL method range from 10.61% to 12.32%, and those of the Fractal Modification model range from 8.07% to 13.16%, both significantly higher than those of the fractional model. As illustrated in (a) and (b), these models exhibit systematic overestimation in the early stage of setting (t < 60 min), particularly for W/C = 0.55, where the Fractal Modification model reaches an error of 13.16%. This behavior likely stems from the sensitivity of these models to the strong scattering effects induced by the rigid flocculated network formed at later hydration stages.

The classical Kelvin–Voigt model exhibits the highest instability among all models, with relative errors ranging from 9.27% to 18.42%. Notably, this model shows its largest error at a low water–cement ratio (W/C = 0.40, 18.42%), whereas the error becomes relatively smaller (9.27%) at a medium ratio (W/C = 0.50). This inconsistency reveals the difficulty of integer-order viscoelastic models in uniformly describing the complex viscoelastic behavior of cement pastes with different mix proportions.

Regarding the trend with water–cement ratio, prediction errors increase for all models as W/C increases from 0.40 to 0.55. For instance, the fractional model’s error increases from 1.85% to 7.11%, the DFL method from 10.61% to 10.71%, and the Fractal Modification model from 8.07% to 13.16%. This reflects that higher water–cement ratios introduce more free water and greater microstructural heterogeneity, thereby increasing the uncertainty in ultrasonic attenuation behavior. Despite this, the fractional model retains the highest prediction stability across all conditions: its maximum error (7.11%) is still lower than the minimum error of all other models, highlighting its advantage in modeling time-dependent material behavior through the fractional order $[eqn]$ .

Overall, these results demonstrate that the fractional model driven by ultrasonic sensing can accurately capture the microstructural evolution during cement paste setting. The fractional parameter $[eqn]$ serves as a real-time, sensitive indicator of the transition from viscous fluid to elastic solid. Compared with conventional rheometers—which require sample extraction and may disturb the hydration process—the proposed method enables continuous, non-invasive monitoring. Moreover, the model’s accurate description of low-frequency attenuation characteristics provides a solid theoretical basis for developing embedded real-time inversion algorithms.

From an SHM viewpoint, the monotonic and physically interpretable evolution of $[eqn]$ observed in the experiments highlights its potential as a state variable rather than a purely empirical fitting parameter. During early hydration, $[eqn]$ captures the transition from viscous-dominated to elastic-dominated behavior, while in later stages or under service conditions, changes in $[eqn]$ may reflect alterations in the internal stress state, microcrack density, or damage-induced viscoelasticity. Compared with traditional ultrasonic velocity or amplitude indicators, $[eqn]$ integrates attenuation behavior over a broad frequency band and is therefore less sensitive to measurement noise and boundary effects, which is advantageous for field-scale SHM applications.

The sensitivity of the proposed fractional-order ultrasonic model to aggregates, fibers, or reinforcement should be carefully discussed. In practical concrete, coarse aggregates introduce strong multiple scattering and impedance contrast, fibers induce anisotropic attenuation and dispersion, and steel reinforcement causes mode conversion and waveguiding effects. These factors primarily affect the propagation path and scattering characteristics of ultrasonic waves rather than the intrinsic viscoelastic behavior of the cementitious matrix.

From a modeling perspective, the fractional parameter $[eqn]$ represents an effective viscoelastic state of the load-bearing matrix at the scale of ultrasonic wavelengths. In heterogeneous composites, $[eqn]$ can therefore be interpreted as an equivalent or homogenized parameter that incorporates matrix behavior and mesoscale interactions. However, additional correction terms—such as aggregate-induced scattering losses or reinforcement-related boundary effects—would need to be explicitly incorporated into the attenuation model for quantitative interpretation in full concrete or reinforced concrete systems.

It should be acknowledged that the fractional parameter $[eqn]$ is influenced by multiple physical factors and should not be interpreted as a hydration-exclusive variable in a strict sense. In the present study, $[eqn]$ primarily reflects hydration-driven microstructural evolution because the experiments were conducted under well-controlled temperature and moisture conditions, with no external mechanical loading. Under such conditions, the dominant contribution to changes in $[eqn]$ originates from the progressive formation of the cementitious network and the associated viscoelastic transition.

In practical SHM applications, however, additional environmental and operational factors may affect the inferred value of $[eqn]$ . Temperature fluctuations can modify relaxation times and shift the frequency-dependent viscoelastic response, leading to apparent changes in $[eqn]$ . Moisture gradients may alter pore-fluid distribution and local damping behavior, while mechanical loading or microcracking can introduce irreversible changes in energy dissipation mechanisms. These effects are expected to influence ultrasonic attenuation through different physical pathways compared with hydration-induced structural build-up.

From an identifiability perspective, $[eqn]$ should therefore be regarded as an effective state parameter that integrates multiple contributions to viscoelastic attenuation. Decoupling hydration effects from environmental or damage-related influences requires auxiliary information and compensation strategies. For example, temperature effects can be mitigated through explicit temperature correction or multi-frequency normalization, while moisture-related variations may be constrained using humidity sensors or coupled poroelastic models. Load- or damage-induced changes in $[eqn]$ are typically non-monotonic and irreversible, in contrast to the monotonic and smooth evolution observed during early-age hydration, providing a potential criterion for distinguishing different mechanisms.

Consequently, the practical interpretation of $[eqn]$ in SHM should rely on baseline referencing and trend analysis rather than absolute values alone. When combined with environmental sensing, multi-parameter inversion, or digital twin frameworks, $[eqn]$ can be effectively decoupled from external influences and used as a robust indicator of material state evolution and structural integrity.

Building upon the identifiability framework discussed above, the practical deployment of the proposed sensing approach depends on the robustness and stability of the ultrasonic measurement system under realistic conditions. While the laboratory system demonstrates reliable performance under controlled environments, its field applicability requires careful consideration of long-term sensor durability, signal stability, and installation-related variability. In cementitious environments, the durability of embedded piezoelectric (PZT) transducers is mainly governed by mechanical protection, chemical compatibility, and moisture ingress; previous SHM studies have shown that appropriate encapsulation or inert housing can ensure stable electromechanical performance even under alkaline and moisture-rich conditions.

From a signal interpretation perspective, the proposed framework relies on baseline-referenced attenuation trends rather than absolute signal amplitudes, which inherently enhances robustness against slow signal drift caused by temperature fluctuations, curing-induced stiffness evolution, or gradual coupling changes over days or weeks. Provided that a sufficient signal-to-noise ratio is maintained, the fractional parameter $[eqn]$ can be consistently identified through trend analysis, supporting long-term monitoring applications. Although sensor coupling conditions and placement may affect absolute attenuation levels, their influence on the temporal evolution of $[eqn]$ is expected to be secondary when consistent excitation–reception paths, standardized embedding procedures, and network-based redundancy are employed. Overall, these considerations indicate that the sensing system and interpretation strategy are compatible with in situ SHM applications beyond laboratory settings.

5. Conclusions

This study presented a physics-informed, fractional-order ultrasonic sensing framework for monitoring the early-age evolution of cementitious materials. By integrating a fractional Zener viscoelastic model with ultrasonic attenuation theory, a quantitative and physically interpretable relationship was established between hydration-driven microstructural changes and measurable acoustic attenuation. The fractional parameter $[eqn]$ was shown to effectively characterize the continuous transition from viscous-dominated to elastic-dominated behavior during early hydration.

A custom ultrasonic measurement system, consisting of a low-frequency ultrasonic transducer, power amplifier, high-resolution data acquisition unit, and temperature compensation module, was developed to validate the proposed sensing approach. Experiments conducted on cement pastes with various water–cement ratios confirmed that the fractional-order model achieves significantly higher accuracy and robustness in capturing attenuation evolution compared with traditional integer-order and empirical models. The temporal evolution of $[eqn]$ provides a clear, physically grounded indicator of structural build-up and agrees well with the known rheological and microstructural transitions of cement hydration.

It should be emphasized that, while this work concentrates on early-age cementitious materials, the proposed fractional-order ultrasonic framework is not limited to short-term hydration monitoring. On the contrary, early-age sensing establishes a physically consistent reference state that is essential for long-term SHM applications. By providing an interpretable viscoelastic baseline, the fractional parameter $[eqn]$ enables the differentiation between intrinsic material evolution, environmental effects, and damage-induced anomalies during the service life of concrete structures. This capability is particularly valuable for intelligent construction, embedded sensing systems, and digital twin-based SHM, where material-scale information must be continuously linked to structural-scale performance assessment.

The proposed framework provides a promising foundation for next-generation, model-based ultrasonic sensing technologies aimed at quality assessment and structural state evaluation in high-end engineering applications, including large-scale concrete components, tunnel linings, marine foundations, and intelligent construction systems. By introducing a physically interpretable fractional parameter $[eqn]$ , the framework establishes a direct link between microstructural evolution and measurable ultrasonic attenuation, offering a compact and robust sensing strategy that is well suited for integration into sensor-empowered intelligent construction and digital twin frameworks for cement-based infrastructure.

Nevertheless, the present study is subject to certain limitations. All experimental investigations were conducted on cement paste, and the direct quantitative applicability of the results to concrete containing aggregates, fibers, or reinforcement has not yet been experimentally validated. In heterogeneous cementitious systems, coarse aggregates, fibers, and reinforcement are expected to introduce additional scattering, anisotropy, and multi-path propagation effects, which can significantly influence ultrasonic wave attenuation and dispersion. Consequently, the current framework should be regarded as a matrix-level sensing and modeling approach that captures the intrinsic viscoelastic evolution of cementitious materials, rather than a complete description of full concrete or reinforced concrete systems.

Future work will focus on extending the proposed framework by incorporating effective medium theories, multi-scattering correction models, and guided-wave formulations to explicitly account for aggregate- and reinforcement-induced effects. In parallel, efforts will be directed toward integrating the sensing model into embedded ultrasonic sensor networks and automated construction platforms, as well as developing real-time parameter inversion algorithms for field deployment. These extensions will enable the transition from controlled material-scale sensing to component- and structural-scale SHM, further supporting lifecycle monitoring and digital twin applications for concrete infrastructure.