Research on the Modulus Decay Model under a Three-Dimensional Stress State of Asphalt Mixture during Fatigue Damage

: Fatigue damage can develop within asphalt pavement due to the continuous impact of driving loads and natural elements. Understanding the process of asphalt mixtures’ fatigue damage is crucial for guiding the design, maintenance, and repair of asphalt pavement. This research aims to establish a model that characterizes the mixtures’ modulus decay under a three-dimensional (3-D) stress state. Firstly, asphalt mixes were subjected to direct tensile (DT), indirect tensile (IDT), unconﬁned compressive (UC) strength and fatigue tests, and the resulting data were analyzed. Then, modulus decay models under DT, IDT, and UC conditions were established, and the modulus decay patterns under the three loading modes were compared and analyzed. Finally, using the fatigue stress strength ratio ∆ (a fatigue resistance index for asphalt mixtures that takes into account the impacts of stress state and loading rate), a uniﬁed characterization model for asphalt mixes’ modulus decay under a 3-D stress state was created. According to the study’s ﬁndings, asphalt mixes’ modulus decay during fatigue damage exhibits obvious nonlinear characteristics. While the asphalt mixes’ modulus decay law with various loading modalities is similar under the same conditions, the decay rate may differ. Essentially, the speed of the modulus decay of a certain asphalt mixture primarily depends on the value of ∆ during service. A larger ∆ indicates a faster modulus decay. This study offers a theoretical foundation for the conversion from material fatigue damage to structural fatigue damage, which is vital for enhancing the asphalt pavements’ construction quality and longevity.


Introduction 1.Background
In the discipline of solid mechanics, known as fatigue damage mechanics [1][2][3], the evolution of damage factors [4], like internal stress, strain, and modulus of materials subjected to repeated loading over time, is primarily studied.Due to factors like vehicular loads and natural environmental conditions, asphalt pavement may suffer from fatigue damage within the structure [5,6].With the accumulation of damage, the materials and structural performance of the pavement will progressively deteriorate until they eventually fail [7].Optimizing the design of the pavement structure and carrying out the necessary maintenance procedures throughout its service life is crucial for improving the construction quality and durability of asphalt pavements.Further, investigating the development process of asphalt mixes' fatigue damage is critical for optimizing pavement structure design and making informed maintenance decisions [8,9].Therefore, numerous indoor experimental studies have been conducted by road researchers on asphalt mixtures [10][11][12][13][14][15], simulating various service conditions of asphalt pavement through diverse fatigue test methods.The asphalt mixes' fatigue performance was examined in these investigations, and fatigue damage models were suggested for several scenarios.Specifically, Zheng et al. [10] established a nonlinear fatigue damage equation for asphalt mixes based on the S-N fatigue equation and damage mechanics theory.Zhou et al. [11] simulated the asphalt mixes' fatigue damage process using three-point bending fatigue tests and a nonlinear finite element method, respectively.They predicted asphalt concrete specimens' fatigue life and crack length during unstable fracture.To evaluate the stiffness variations of two types of asphalt mixes at rest brought on by temperature, fatigue damage growth, and healing, Daniel et al. [12] used the impact resonance method.Shan et al. [13] utilized CT scanning technology to capture internal cross-sectional images of asphalt mixtures at various stages of fatigue damage during IDT fatigue testing.Subsequently, they examined the progression of fatigue damage over time.Based on four-point bending tests, Izaks et al. [14] found that developing High Modulus Asphalt Concrete (HMAC) mixtures with high fatigue resistance performance using RA content up to 70% by using polymer-modified bitumen binder is possible.Lv et al. [15] conducted small-scale accelerated pavement tests on asphalt mixtures using the Model Moving Load Simulator (MMLS3) to evaluate their fatigue damage performance.
It should be noted that one crucial mechanical parameter that plays a significant role in determining the asphalt pavement's service life is the modulus of pavement material [16][17][18].An essential factor to take into account in pavement design and maintenance is the decay pattern of the modulus, which can serve as a general measure of asphalt pavement's damage law [18][19][20][21][22].Thus, according to the strain equivalence assumption of viscoelastic damage mechanics theory, some scholars defined the reduction in modulus (stiffness) as a damage variable to study the asphalt mixes' fatigue damage features.For example, Zhou et al. [20] established a damage evolution equation considering the joint impact of fatigue and creep damage effects.Based on the small-scale accelerated pavement test results, Lv et al. [15] established a stiffness decrease equation for asphalt mixes and proposed the corresponding reduction law during fatigue.Lv et al. [21] also analyzed the modulus decay characteristics of asphalt mixtures with various loading circumstances.By utilizing DT fatigue testing, Liu et al. [22] revealed the attenuation law of remaining strength and dynamic modulus.Additionally, they established a coupled fatigue damage evolution model of remaining strength and dynamic modulus.
However, two aspects need to be noted when studying the modulus decay law of asphalt mixtures.Firstly, the modulus of asphalt pavement materials has a specific stress dependence.As asphalt pavement materials are viscoelastic and composed of multiple phases, scales, and components, their modulus is not an inherent property of asphalt mixtures and is influenced by stress states [23,24].Secondly, the modulus of asphalt pavement materials is influenced by the pavement's structure.Modulus levels and stress states change as a result of the different layers of asphalt mixture that make up the pavement structure.
During the actual service period, asphalt pavement structures experience a 3-D stress condition [25].Consequently, there may be significant discrepancies between the results of indoor material tests conducted using one particular stress condition and the actual pavement service state, which hinders the precise assessment of the asphalt pavement structures' fatigue performance and impedes scientific maintenance decision-making.Although full-scale accelerated loading tests [26][27][28] can provide a more realistic simulation of the actual service state of asphalt pavement, their high testing costs and lengthy testing cycles have limited their widespread use.Therefore, an urgent need is to develop a model that can accurately reflect the modulus decay of asphalt mixtures under actual 3-D service conditions through indoor experiments.
This study focuses on understanding the process of fatigue damage in asphalt mixtures under the three-dimensional service condition, which is crucial for guiding the design and maintenance of asphalt pavement structures.On the basis of the fatigue stress strength ratio ∆ as a measure of durability against fatigue [29], a unified characterization model for the asphalt mixes' modulus decay under a 3-D stress condition was developed.The findings of this study provide a theoretical foundation for the reasonable conversion of structural fatigue damage from material fatigue damage, allowing for a more realistic characterization of the asphalt pavement structures' fatigue damage patterns through indoor material tests, which is beneficial for improving asphalt pavements' construction quality and service life.

Objective
This paper aims to establish a model that characterizes the asphalt mixtures' modulus decay under a 3-D stress state during the fatigue damage.Strength and fatigue tests were conducted under the DT, IDT, and UC modes.The primary objectives of this research are as follows: (1) The relationship between the asphalt mixtures' strength and loading rate was established.(2) The conventional S-N fatigue equations for asphalt mixtures under the DT, IDT, and UC modes were obtained and analyzed.Additionally, the normalized fatigue equation for asphalt mixtures under the various modes was established based on the fatigue stress strength ratio ∆.
(3) Through stress-controlled fatigue tests, the modulus decay equations of asphalt mixtures under the DT, IDT, and UC modes were developed and compared.(4) A unified characterization model for the modulus decay of asphalt mixture under the 3-D stress state was established based on the fatigue stress strength ratio ∆.

Asphalt
Table 1 displays the conventional performance test outcomes of the high-content (30%) rubber-modified asphalt employed in this investigation.The test methods followed the specification JTG E20-2011 [30].  2 and 3, respectively.Table 4 shows the basic technical indicators of the filler.Each index complies with the JTG F40-2004 specification [31].The SAC-13 graded rubber asphalt mixture [32][33][34] was utilized in this study.The design principle of the SAC mix is to form a skeleton with coarse aggregates and fill the voids in the structure with fine aggregates and asphalt.SAC asphalt mixture has excellent high-temperature deformation resistance, water damage resistance, and slip resistance [35].Additionally, it has advantages such as small segregation, easy compaction, and asphalt conservation [36].
The grading design for SAC-13 [35,36] was conducted, which included a maximum particle size of 13.2 mm, a control percent passing of 30% for the 4.75 mm sieve, and 8% for the 0.075 mm sieve.The grading curve is plotted as shown in Figure 1, and the final asphalt mixture asphalt-aggregate ratio is 5.5%.Table 5 displays the Marshall performance indicators of the SAC-13 mixes at the optimal asphalt-aggregate ratio.

Experiment Approach
This study tested asphalt mixtures for strength and fatigue using the DT, IDT, and UC modes.The strength test adopts two methods: stress and displacement control.The fatigue test determines a sample's fatigue life by using a stress-control mode and counting the number of loading cycles till a specimen fails, either marked by obvious fracture cracking or deformation.Strength and fatigue tests were conducted on an MTS testing machine.LVDT sensors are utilized during the fatigue test to detect real-time specimen displacement, allowing for the analysis of modulus changes throughout the process.Testing is consistently conducted at 15 ± 1 • C, and the specimen should be insulated at that temperature for 4 h or more before the test.

Experiment Approach
This study tested asphalt mixtures for strength and fatigue using the DT, IDT, and UC modes.The strength test adopts two methods: stress and displacement control.The fatigue test determines a sample's fatigue life by using a stress-control mode and counting the number of loading cycles till a specimen fails, either marked by obvious fracture cracking or deformation.Strength and fatigue tests were conducted on an MTS testing machine.LVDT sensors are utilized during the fatigue test to detect real-time specimen displacement, allowing for the analysis of modulus changes throughout the process.Testing is consistently conducted at 15 ± 1 °C, and the specimen should be insulated at that temperature for 4 h or more before the test.

Specimen Preparation
The DT samples were created following the guidelines of specification JTG E20-2011 [30].Initially, a rectangular specimen (length and width of 300 mm, thickness of 50 mm) was molded using an Asphalt Mixer Wheel Grinding Machine.Subsequently, it was divided into beam samples (length of 250 mm, width and height of 50 mm).The asphalt mixture beam specimen is fixed between two circular iron plates with joints using epoxy stone dry-hanging adhesive, which has a higher tensile strength than the asphalt mixture and causes minimal deformation, ensuring accurate test results.Two displacement sensors are symmetrically attached on either side during fatigue testing of the asphalt mixture beam specimen.Refer to Figure 2a for the DT specimen.
The SGC rotary compactor was utilized to create cylindrical IDT specimens with a height and diameter of 100 ± 2 mm.The compaction process involved a vertical loading pressure of 600 ± 18 kPa, an effective internal rotation angle of 1.16 ± 0.02°, and a

Specimen Preparation
The DT samples were created following the guidelines of specification JTG E20-2011 [30].Initially, a rectangular specimen (length and width of 300 mm, thickness of 50 mm) was molded using an Asphalt Mixer Wheel Grinding Machine.Subsequently, it was divided into beam samples (length of 250 mm, width and height of 50 mm).The asphalt mixture beam specimen is fixed between two circular iron plates with joints using epoxy stone dry-hanging adhesive, which has a higher tensile strength than the asphalt mixture and causes minimal deformation, ensuring accurate test results.Two displacement sensors are symmetrically attached on either side during fatigue testing of the asphalt mixture beam specimen.Refer to Figure 2a for the DT specimen.
The SGC rotary compactor was utilized to create cylindrical IDT specimens with a height and diameter of 100 ± 2 mm.The compaction process involved a vertical loading pressure of 600 ± 18 kPa, an effective internal rotation angle of 1.16 ± 0.02 • , and a compaction speed of 30 ± 0.5 r/min.Afterward, the specimens were cut into 60 ± 2 mm height and 100 ± 2 mm diameter specimens for indirect tensile testing, as displayed in Figure 2b.
The UC samples, with a height and diameter of 100 ± 2 mm, were produced using the SGC rotary compactor.This method is consistent with the indirect stretching method discussed earlier.The unconfined compressive specimen is shown in Figure 2c.compaction speed of 30 ± 0.5 r/min.Afterward, the specimens were cut into 60 ± 2 mm height and 100 ± 2 mm diameter specimens for indirect tensile testing, as displayed in Figure 2b.
The UC samples, with a height and diameter of 100 ± 2 mm, were produced using the SGC rotary compactor.This method is consistent with the indirect stretching method discussed earlier.The unconfined compressive specimen is shown in Figure 2c.

Strength Test Approach
The MTS testing machines were employed to conduct the strength tests.For the stress-control mode DT, IDT, and UC strength tests, eight loading rates ranging from 0.02 MPa/s to 6 MPa/s were utilized.Before the fatigue test, displacement-control mode strength tests were performed.During the displacement-control mode, the loading rates for the DT, IDT, and UC strength tests were maintained at 5 mm/min, 50 mm/min, and 2 mm/min, respectively, following JTG E20-2011 standards [30].Each group of tests consisted of three parallel specimens, and the maximum load (N) was recorded upon specimen failure.The DT, IDT, and UC strength tests are shown in Figure 3.
The formula for calculating DT strength: The formula for calculating IDT strength: 0.006287 where IDT

S
represents IDT strength (MPa); h denotes specimen height (mm); P is the maximum test load (N).
The formula for calculating UC strength: where UC S represents UC strength (MPa); P is the maximum load (N) at the point of specimen failure; D denotes specimen diameter (mm); A is the cross-sectional area (mm 2 ).

Strength Test Approach
The MTS testing machines were employed to conduct the strength tests.For the stresscontrol mode DT, IDT, and UC strength tests, eight loading rates ranging from 0.02 MPa/s to 6 MPa/s were utilized.Before the fatigue test, displacement-control mode strength tests were performed.During the displacement-control mode, the loading rates for the DT, IDT, and UC strength tests were maintained at 5 mm/min, 50 mm/min, and 2 mm/min, respectively, following JTG E20-2011 standards [30].Each group of tests consisted of three parallel specimens, and the maximum load (N) was recorded upon specimen failure.The DT, IDT, and UC strength tests are shown in Figure 3.

Fatigue Test Approach
In this study, the selection of stress levels for the DT, IDT, and UC fatigue tests was based on the strength (static load strength) under displacement-control mode, respectively.Moreover, the selected stress level should not exceed the static load strength.Four specimens were tested in parallel for fatigue at each stress level, with a loading frequency of 10 Hz and a continuous haversine wave loading waveform.Refer to Figure 4 for the fatigue test schematic diagram.The formula for calculating DT strength: where S DT represents DT strength (MPa); P is the maximum load (N) at the point of specimen failure; a, b represents the length and width of the cross section (a = b = 50 mm); A is the cross-sectional area (mm 2 ).The formula for calculating IDT strength: where S IDT represents IDT strength (MPa); h denotes specimen height (mm); P is the maximum test load (N).
The formula for calculating UC strength: where S UC represents UC strength (MPa); P is the maximum load (N) at the point of specimen failure; D denotes specimen diameter (mm); A is the cross-sectional area (mm 2 ).

Fatigue Test Approach
In this study, the selection of stress levels for the DT, IDT, and UC fatigue tests was based on the strength (static load strength) under displacement-control mode, respectively.Moreover, the selected stress level should not exceed the static load strength.Four specimens were tested in parallel for fatigue at each stress level, with a loading frequency of 10 Hz and a continuous haversine wave loading waveform.Refer to Figure 4

Fatigue Test Approach
In this study, the selection of stress levels for the DT, IDT, and UC fatigue tests was based on the strength (static load strength) under displacement-control mode, respectively.Moreover, the selected stress level should not exceed the static load strength.Four specimens were tested in parallel for fatigue at each stress level, with a loading frequency of 10 Hz and a continuous haversine wave loading waveform.Refer to Figure 4

Damage Variables Based on Modulus
When subjected to the factors such as forces, temperature, and loading speed, microdefects can form within a material.These defects can nucleate, expand, and converge, ultimately leading to material failure [7].Variations in macroscopic physical variables like modulus, stress, and strength frequently accompany this process of material deterioration.Specific physical quantities can be defined as damage variables to study the microscopic changes within materials during the damage process, and corresponding fatigue damage models can be established to describe the damage process [4].
In 1963, the mechanist Rabotnov used the following damage factor D to describe the damage and study the constitutive equations of metal materials [37].

Derivation of Modulus Decay Equation 2.3.1. Damage Variables Based on Modulus
When subjected to the factors such as forces, temperature, and loading speed, microdefects can form within a material.These defects can nucleate, expand, and converge, ultimately leading to material failure [7].Variations in macroscopic physical variables like modulus, stress, and strength frequently accompany this process of material deterioration.Specific physical quantities can be defined as damage variables to study the microscopic changes within materials during the damage process, and corresponding fatigue damage models can be established to describe the damage process [4].
In 1963, the mechanist Rabotnov used the following damage factor D to describe the damage and study the constitutive equations of metal materials [37].
where A is the effective bearing area; A is the nominal area, which is the area in the undamaged state.The fatigue test of the mixes allows for continuous measurement of the modulus, enabling the asphalt mixes' modulus decay during the test to be used as a damage variable.
In this study, a damage factor formula based on the modulus was defined to describe damage in fatigue tests, as shown below.
where D is the damage variable of the specimen; E is the modulus of the material after fatigue; E 0 is the initial modulus of the specimen.

Nonlinear Fatigue Damage Equation for Asphalt Mixes
Creating a fatigue damage model is crucial in studying material fatigue damage, as it links the number of load actions to the resulting material damage.On the one hand, materials' nonlinear fatigue damage law is related to their initial damage state, and the amount of damage that occurs at different damage stages, even under the same load state, is not the same.On the other hand, the nonlinear fatigue damage law of materials varies with different load levels during cyclic loading.The most intuitive manifestation of the nonlinearity of material fatigue damage is that the image describing the material fatigue damage model is nonlinear.Equation ( 6) is a commonly used nonlinear fatigue damage model considering the influence of stress amplitude.
By integrating Equation ( 6), the fatigue damage variable formula can be derived as follows: Correspondingly, Equation ( 8) is the expression of the material fatigue equation: where α, γ and M are the material parameters linked to the average stress, temperature, and stress amplitude; D represents the damage variable; N stands for the number of cycles of load in the fatigue test; N f represents the fatigue life of the material.

Modulus Decay Equation
Fatigue tests allow for continuous modulus measurement without affecting asphalt mixture performance.The modulus steadily decreases as internal damage accumulates within the material, making it an effective measure for describing damage in asphalt mixtures.The simultaneous Equations ( 5) and ( 7) are as follows: The modulus decay equation can be obtained by eliminating the damage variable D in Equation (9).
where m is the modulus decay equation parameter, which can be fitted through fatigue tests; E is the effective modulus; E 0 is the initial modulus; N represents the number of load cycles in the fatigue test; N f denotes the fatigue life.Displacement-controlled strength tests for asphalt mixture were performed implementing the procedure in Section 2.2.2, and the experimental results are shown in Table 6.A stress-control mode was adopted to perform the DT, IDT, and UC strength tests at various loading rates to investigate the impact of the loading rate on asphalt mix strength.The experiment outcomes, presented in Table 7, provide insights into the viscoelastic behavior of asphalt mixtures.Equation (12) shows that the strength and loading rate were fitted using a power function equation.The curves are displayed in Figure 5.

Results and Discussion
where S is the strength of asphalt mixture; A and B are the equation parameters; v is the loading rate.Figure 5 demonstrates that the loading rate positively affects the strength of asphalt mixtures for all stress states, but the pace of increase reduces with higher loading rates.Unconfined compressive mode yields much higher strength values than the direct and indirect tensile modes, while the direct and indirect tensile modes produce similar strength values at the same loading rate.
where S is the strength of asphalt mixture; A and B are the equation parameters; v is the loading rate.
Figure 5 demonstrates that the loading rate positively affects the strength of asphalt mixtures for all stress states, but the pace of increase reduces with higher loading rates.Unconfined compressive mode yields much higher strength values than the direct and indirect tensile modes, while the direct and indirect tensile modes produce similar strength values at the same loading rate.

Analysis of Fatigue Test Results
Fatigue experiments were performed on asphalt mixtures by applying the experimental procedure outlined in Section 2.2.3.Before the fatigue test, a displacement-control mode strength test was conducted, and the outcomes are presented in Table 6.On this basis, five stress levels ranging from 0.2 MPa to 1 MPa were chosen for the DT and IDT fatigue tests, while UC fatigue experiments were performed at stress levels of 1.5 MPa, 2 MPa, 2.5 MPa, 3 MPa, and 3.5 MPa.The fatigue test data are summarized in Table 8.

Analysis of Fatigue Test Results
Fatigue experiments were performed on asphalt mixtures by applying the experimental procedure outlined in Section 2.2.3.Before the fatigue test, a displacement-control mode strength test was conducted, and the outcomes are presented in Table 6.On this basis, five stress levels ranging from 0.2 MPa to 1 MPa were chosen for the DT and IDT fatigue tests, while UC fatigue experiments were performed at stress levels of 1.5 MPa, 2 MPa, 2.5 MPa, 3 MPa, and 3.5 MPa.The fatigue test data are summarized in Table 8.Equation ( 13) is the S-N fatigue model based on phenomenological methods [38,39].
Equation ( 14) is the formula for calculating t n .
where N f is the fatigue life; k and c are the regression parameters; t n is the nominal stress ratio; S is the static load strength; σ is the stress levels of fatigue test.Figure 6 illustrates the fitting of the S-N fatigue formulas for asphalt mixtures under various stress states (loading modes) using Equation (12), founded on the fatigue experiment data obtained from the DT, IDT, and UC modes.where f N is the fatigue life; k and c are the regression parameters; n t is the nomi- nal stress ratio; S is the static load strength;  is the stress levels of fatigue test.
Figure 6 illustrates the fitting of the S-N fatigue formulas for asphalt mixtures under various stress states (loading modes) using Equation ( 12), founded on the fatigue experiment data obtained from the DT, IDT, and UC modes.Upon examining Figure 6, it is apparent that the S-N equation of asphalt concrete displays noteworthy discrepancies when subjected to the different loading modes.It is evident in the considerable variation in the fatigue equation parameter k values obtained from the three methods, resulting in a substantial separation between the fatigue equation curves under the DT, IDT, and UC modes.The fundamental reason for this phenomenon is that the nominal stress ratio t n , which serves as a fatigue resistance indicator, fails to account for the impact of the stress state and loading rate on asphalt mixtures as viscoelastic materials.Hence, to accurately assess the fatigue and damage properties of the mixtures, it is imperative to employ a fatigue resistance index that considers the material's 3-D stress condition and loading rate.This approach can significantly mitigate the influence of the experimental loading modes, specimen geometry, and size.

Fitting Results of Modulus Decay Equation with Various Stress States
This study defines the modulus' attenuation as a fatigue damage variable, and the focus of this section is to investigate asphalt mixes' modulus decay equation under various stress conditions during fatigue testing.A modulus value can be calculated in each load cycle of the fatigue test.The calculation method for modulus E is as follows.
where σ is the stress level of fatigue test; ε ei is the strain; E is the modulus of asphalt mixture.The Equation (11) derived in Section 2.  Figure 7 shows that the parameter m in the modulus decay equation positively correlates with the stress level.It means that as the stress level increases, the dynamic modulus decays faster.

Comparative Analysis of Modulus Decay Law under Different Stress States
Because the stress levels during fatigue testing differ between the UC mode and DT/IDT modes, the modulus decay equations for the three modes cannot be directly compared.To address this, the nominal stress ratio was obtained by dividing the stress levels equally by the static load strength for dimensionless processing.As a result, the modulus decay equations under the DT, IDT, and UC modes could be compared by examining the parameter m under the same nominal stress ratio.The m values and nominal stress ratios t n for each mode were fitted using a power function, and the resulting fitting curves are displayed in Figure 8.  pared.To address this, the nominal stress ratio was obtained by dividing the stress levels equally by the static load strength for dimensionless processing.As a result, the modulus decay equations under the DT, IDT, and UC modes could be compared by examining the parameter m under the same nominal stress ratio.The m values and nominal stress ratios n t for each mode were fitted using a power function, and the resulting fitting curves are displayed in Figure 8.

Comparative Analysis of Modulus Decay Law under Different Stress States
IDT: In an effort to contrast and analyze the modulus decay patterns with various stress states, it is essential to ensure that the fatigue test load levels remain consistent, i.e., the nominal stress ratio t n should be the same.Taking t n = 0.2, the modulus decay equations for the DT, IDT, and UC modes are: IDT : UC : The modulus decay curves were plotted according to Equations ( 19)-( 21), as shown in Figure 9. IDT: The modulus decay curves were plotted according to Equations ( 19)-( 21), as shown in Figure 9. Figure 9 suggests that the modulus decay patterns are generally similar under the different loading modes despite differences in the stress states.The decay of the modulus first occurs in the stable decay stage, followed by the rapid decay stage.However, the m value of the mixes' modulus decay curve varies under the different loading modes, even when subjected to the same nominal stress ratio.The modulus decay rate from fast to slow is IDT > DT > UC.
This result is likely because the n t does not consider how the loading rate and stress condition affect asphalt mixtures.To address this, the next step in this study will be to establish a unified model for describing asphalt mixes' modulus decay pattern under a 3-D stress condition, using the fatigue stress strength ratio Δ as a fatigue resistance index.Figure 9 suggests that the modulus decay patterns are generally similar under the different loading modes despite differences in the stress states.The decay of the modulus first occurs in the stable decay stage, followed by the rapid decay stage.However, the m value of the mixes' modulus decay curve varies under the different loading modes, even when subjected to the same nominal stress ratio.The modulus decay rate from fast to slow is IDT > DT > UC.

Unified Characterization Model for Modulus Decay of Asphalt Mixture under the 3-D Stress State
This result is likely because the t n does not consider how the loading rate and stress condition affect asphalt mixtures.To address this, the next step in this study will be to establish a unified model for describing asphalt mixes' modulus decay pattern under a 3-D stress condition, using the fatigue stress strength ratio ∆ as a fatigue resistance index.

Fatigue Stress Strength Ratio and Normalized Fatigue Equation
The stress condition inside asphalt pavement is highly complex under actual driving loads.Asphalt mixtures' strength and fatigue behavior are affected by the stress state and loading rate.Traditional phenomenological fatigue analysis methods fail to consider these factors, resulting in an unscientific approach.To overcome this, Lv et al. [29,39] proposed a novel approach for analyzing the fatigue characteristic of asphalt mix, which accounts for the actual pavement service conditions and utilizes the idea of a 3-D stress state strength yield surface.This compensates for the limitations of traditional fatigue performance analysis methods that do not consider the impacts of stress conditions and loading rates.The fatigue stress strength ratio ∆ is served as a fatigue resistance index, with Equation ( 22) defining its formula.
where σ ei represents the initial equivalent stress, while σ es represents the equivalent stress at the point of fatigue failure.The formula for equivalent stress is: where σ 1 is the first principal stress; σ 2 is the second principal stress; σ 3 is the third principal stress.Equation (22) shows the ratio of the stress level to material resistance during asphalt mixture fatigue testing.Using the ratio ∆, which considers the stress state and loading rate, a normalized fatigue equation is established as follows.
where N f represents fatigue life; γ represents the regression parameter.
The initial equivalent stress σ ei of the fatigue test and the equivalent stress during fatigue failure σ es can be calculated according to Tables 12 and 13, respectively.

Loading Modes
Stress Levels σ i (MPa) Table 13.Calculation method of σ es .

Loading Modes Stress-Controlled
In this study, the fatigue tests apply a haversine waveform load.The loading rate v can be calculated by the test's loading frequency f (loading cycle T) and stress level σ.
Based on Equation (24) and the power function equations established in Section 3.1.1for the strength and loading rate of asphalt mixtures in the DT, IDT, and UC modes, the ultimate strength at various stress levels can be determined.Thus, σ es can be calculated.Table 14 displays the calculated ∆ values.
Based on the data of fatigue life N f (in Section 3.1.2)and the fatigue stress strength ratio, the normalized fatigue equation for various stress conditions can be fitted by Equation (24).The resulting curve for the asphalt mix is displayed in Figure 10.
Figure 10 illustrates that the use of the fatigue stress strength ratio ∆ has successfully achieved the uniform description of asphalt mixtures' fatigue behavior with the various stress states.As a fatigue resistance indicator, the ratio ∆ considers the impact of the loading mode and rate on asphalt mixtures and correlates with the stress state and vehicle speed within a pavement structure.Based on the ratio ∆, this study established a modulus decay model for asphalt mixes under a 3-D stress condition.

Establishment of the Modulus Decay Model under the 3-D Stress State
Section 3.2 developed the modulus decay equations of asphalt mixes under various stress conditions corresponding to the different loading modes.In Section 3.3.1, the fatigue stress strength ratios were calculated for the DT, IDT, and UC modes at different stress levels.Table 15 displays the values of the modulus decay equation parameter m and the fatigue stress strength ratio ∆. m and ∆ are fitted using a power function equation, as shown in Figure 11.(24).The resulting curve for the asphalt mix is displayed in Figure 10.
Figure 10 illustrates that the use of the fatigue stress strength ratio Δ has successfully achieved the uniform description of asphalt mixtures' fatigue behavior with the various stress states.As a fatigue resistance indicator, the ratio Δ considers the impact of the loading mode and rate on asphalt mixtures and correlates with the stress state and vehicle speed within a pavement structure.Based on the ratio Δ, this study established a modulus decay model for asphalt mixes under a 3-D stress condition.
Therefore, a modulus decay model of asphalt mixture under the 3-D stress condition, Equation (27), was established by uniting Equations (11) and (26).Section 3.2 developed the modulus decay equations of asphalt mixes under various stress conditions corresponding to the different loading modes.In Section 3.3.1, the fatigue stress strength ratios were calculated for the DT, IDT, and UC modes at different stress levels.Table 15 displays the values of the modulus decay equation parameter m and the fatigue stress strength ratio Δ. m and Δ are fitted using a power function equa- tion, as shown in Figure 11.
where E represents the modulus during fatigue, while E 0 is its initial value; N represents the number of fatigue load cycles, while N f represents the fatigue life; m is a parameter of the modulus decay equation that can be determined through fatigue tests; ∆ is the fatigue stress strength ratio, which accounts for the effects of the stress state and loading rate, and represents the ratio of the stress level to material resistance during fatigue testing of asphalt mixtures; n is the regression parameter.
where E represents the modulus during fatigue, while 0 E is its initial value; N rep- resents the number of fatigue load cycles, while f N represents the fatigue life; m is a parameter of the modulus decay equation that can be determined through fatigue tests; Δ is the fatigue stress strength ratio, which accounts for the effects of the stress state and loading rate, and represents the ratio of the stress level to material resistance during fatigue testing of asphalt mixtures; n is the regression parameter.
Figure 12 shows the modulus decay surface under the 3-D stress state, which was drawn using MATLAB R2018a based on Equation (27).The surface function is 0.328 (1 ) , (0,1), (0,1) x Z y x y     . Figure 12 illustrates that the modulus decay model ( 27) can describe the modulus decay behavior under any service state (represented by fatigue stress strength ratio Δ).In asphalt pavement structures, the fatigue stress strength ratio at a given point continuously changes during service.Analyzing the variation law of the fatigue stress strength ratio at a particular point in the asphalt pavement structure under different fatigue life ratios, and using Equation ( 27), a structural modulus decay model can be established for that point in the asphalt pavement structure.Theoretically, the modulus decay curve of an asphalt pavement structure should be a curve on the asphalt mixes' modulus decay curved Figure 12 illustrates that the modulus decay model (27) can describe the modulus decay behavior under any service state (represented by fatigue stress strength ratio ∆).In asphalt pavement structures, the fatigue stress strength ratio at a given point continuously changes during service.Analyzing the variation law of the fatigue stress strength ratio at a particular point in the asphalt pavement structure under different fatigue life ratios, and using Equation ( 27), a structural modulus decay model can be established for that point in the asphalt pavement structure.Theoretically, the modulus decay curve of an asphalt pavement structure should be a curve on the asphalt mixes' modulus decay curved surface under the 3-D stress condition.As a result, the relationship between material modulus decay and structural modulus decay has been established.
This study establishes a theoretical foundation for the reasonable conversion of structural fatigue damage from material fatigue damage, making it possible to characterize the fatigue damage pattern of asphalt pavement structural layers through indoor experiments.

Conclusions
This study investigated the fatigue damage properties of asphalt mixtures under the DT, IDT, and UC modes.Firstly, the strength and fatigue experiment results were analyzed, followed by establishing the modulus decay models under the different loading modes.The modulus decay patterns were then compared and analyzed.Using the fatigue stress strength ratio ∆ as a measure of fatigue resistance, a unified description model for the modulus decay of asphalt mixtures under the 3-D stress condition was established.In summary, the study provides the following conclusions.
(1) The power function relationships between asphalt mixtures' strength and loading rate v are as follows.
S DT = 2.275v 0.213 , R 2 = 0.965 S IDT = 1.994v 0.195 , R 2 = 0.978 S UC = 11.602v0.154 , R 2 = 0.985 The strength of the asphalt mixture rises with the loading rate during testing.However, under various loading approaches, the asphalt mixes' strength varies even at the same loading rate.
(2) The conventional S-N fatigue equations of asphalt mix under the DT, IDT, and UC methods are as follows.
DT : N f = 2594.06(1/tn ) 3.457 , R 2 = 0.966 IDT : N f = 14.636(1/t n ) 4.469 , R 2 = 0.960 UC : N f = 109.08(1/tn ) 4.438 , R 2 = 0.894 The S-N fatigue equation of the mix exhibits significant differences under the DT, IDT, and UC modes.This phenomenon can be attributed to the fact that the nominal stress ratio t n , as a fatigue resistance indicator, does not account for the impacts of the stress conditions and loading rates on asphalt mixtures as viscoelastic materials.
(3) The modulus decay patterns are generally similar across the different loading modes.Modulus decay occurs first in the stable decay stage and then in the rapid decay stage.However, the parameter m of the modulus decay curve varies across the different loading modes under the same nominal stress ratio t n , with the modulus decay being fastest for IDT, followed by the DT and UC modes in descending order.
(4) Using the fatigue stress strength ratio ∆, the fatigue features of asphalt mixtures with various loading conditions can be uniformly described.The normalized fatigue equation for asphalt mixtures under different stress conditions was established.
N f = (1/∆) 5.109 , R 2 = 0.908 The fatigue resistance index ∆ considers the impact of the loading mode and rate on asphalt mixtures.It can correspond to the stress condition within the pavement structure and vehicle driving speed.
(5) A unified characterization model for the modulus decay of asphalt mixture under the 3-D stress state was established based on the ratio ∆.
This model comprehensively characterizes the asphalt mixes' modulus decay law under various service conditions, providing a reference to understanding the fatigue damage behavior of asphalt pavement structures.The modulus decay curve of an asphalt pavement structure should theoretically be a curve in the modulus decay curved surface of asphalt mixture under the 3-D stress condition, thus establishing a relationship between material modulus decay and structural modulus decay.
represents DT strength (MPa); P is the maximum load (N) at the point of specimen failure; , a b represents the length and width of the cross section ( 50 a b   mm); A is the cross-sectional area (mm 2 ).
for the fatigue test schematic diagram.
for the fatigue test schematic diagram.

Figure 5 .
Figure 5. Fitting results of strength and loading rate.

Figure 5 .
Figure 5. Fitting results of strength and loading rate.

Figure 6 .
Figure 6.S-N fatigue equations under different stress states.Upon examining Figure 6, it is apparent that the S-N equation of asphalt concrete displays noteworthy discrepancies when subjected to the different loading modes.It is evident in the considerable variation in the fatigue equation parameter k values ob- tained from the three methods, resulting in a substantial separation between the fatigue equation curves under the DT, IDT, and UC modes.The fundamental reason for this

Figure 6 .
Figure 6.S-N fatigue equations under different stress states.

3 . 3
was used to fit the modulus ratio E/E 0 and fatigue life ratio N/N f .The fitted modulus decay model parameter m of DT, IDT, and UC are summarized in Tables 9-11.Examples of modulus decay curves (one of the four parallel UC experiments) are shown in Figure7.

Figure 8 .
Figure 8. Fitting curves of parameter m and n t under the different loading modes.(a) DT.(b)

Figure 8
Figure 8 displays that the parameter m in the mixes' modulus decay equation rises as the nominal stress ratio n t increases.The corresponding fitting equations are pre-

Figure 8 .
Figure 8. Fitting curves of parameter m and t n under the different loading modes.(a) DT.(b) IDT.(c) UC.

Figure 9 .
Figure 9. Modulus decay curves under the different loading modes when n t = 0.2.

Figure 9 .
Figure 9. Modulus decay curves under the different loading modes when t n = 0.2.

Figure 10 .
Figure 10.The normalized fatigue equation curve of asphalt mixture.Figure 10.The normalized fatigue equation curve of asphalt mixture.

Figure 10 .
Figure 10.The normalized fatigue equation curve of asphalt mixture.Figure 10.The normalized fatigue equation curve of asphalt mixture.By observing Figure11, it is evident that the parameter m in the modulus decay equation increases as the ratio ∆ increases.The power function relationship equation between m and ∆, obtained by fitting the data, is presented below.m = ∆ 0.328 , R 2 = 0.910(26)

Figure 11 .
Figure 11.The fitting results of m and Δ.By observing Figure11, it is evident that the parameter m in the modulus decay equation increases as the ratio Δ increases.The power function relationship equation between m and Δ, obtained by fitting the data, is presented below.

Figure 11 .
Figure 11.The fitting results of m and ∆.

Figure 12 .
Figure 12.Modulus decay curved surface under the three-dimensional stress state.

Figure 12 .
Figure 12.Modulus decay curved surface under the three-dimensional stress state.

the
Postgraduate Scientific Research Innovation Project of Changsha University of Science and Technology (No. CXCLY2022012 and No. CLSJCX22016).

Table 1 .
Performance test outcomes of the asphalt in this investigation.
2.1.2.Aggregate and FillerLimestone aggregate and limestone mineral powder filler were utilized in this study.The aggregate density and mechanical index test results are presented in Tables

Table 2 .
Density and water absorption of limestone aggregate.

Table 3 .
Technical indexes of limestone aggregate.

Table 4 .
Technical indexes of filler.

Table 7 .
Strength test results by stress-controlled mode.

Table 8 .
Fatigue test results of asphalt mixture.
In the table: σ is the stress level (MPa); t n is the nominal stress ratio; N f is the mean fatigue life (cycles); CV is the coefficient of variation (%).

Table 9 .
Modulus decay equation parameters m fitting results of DT.

Table 10 .
Modulus decay equation parameters m fitting results of IDT.

Table 11 .
Modulus decay equation parameters m fitting results of UC.

Table 14 .
The calculated ∆ values.(in Section 3.1.2)and the fatigue stress strength ratio, the normalized fatigue equation for various stress conditions can be fitted by Equation f N

Table 15 .
The values of m and ∆.

Table 15 .
The values of m and Δ.