Effect of Weather Aging on Viscoelasticity and Fatigue Performance of Asphalt Mastic

The long-term effect of climate factors, such as sunlight, oxygen, and water, leads to the performance degradation of the asphalt mastic, which is the binding part in the asphalt mixture. It is not conducive to satisfy the long-term performance requirements of long-life asphalt pavement. In this study, five kinds of base asphalt mastic and styrene-butadiene-styrene (SBS) modified asphalt mastic were prepared with the filler-asphalt ratio of 0.6, 0.8, 1.0, 1.2, and 1.4. The indoor simulated weather aging tests were carried out considering multi-factors including sunlight, oxygen, and water. The master curves of the complex shear modulus and phase angle of the asphalt mastic with different aging degrees were obtained by the frequency sweep test. The curves of fatigue damage characteristics and fatigue life were fitted based on the viscoelastic continuum damage (VECD) model. The influence of weather aging on the viscoelasticity and fatigue performance of asphalt mastic were analyzed. Results indicated that the effect of weather aging increases the elastic component and decreases the viscous component. The fatigue performance of SBS modified asphalt mastic was better than that of base asphalt mastic. As the aging degree deepens, the brittle failure characteristics of asphalt mastic with a higher filler–asphalt ratio were more obvious. The base asphalt mastic becomes more sensitive to the strain level due to weather aging, and its fatigue life increased under the low strain loading and decreased under the high strain loading. The fatigue performance of SBS modified asphalt mastic was less sensitive to the strain level. The fatigue life reduced after aging under low and high strain load. Taking the impact of weather aging on the fatigue performance into consideration, the optimal filler–asphalt ratios of the base asphalt mastic SBS modified asphalt mastic are 1.0 and 1.2, respectively.


Introduction
For long-life asphalt pavement with a service life of 40 to 60 years, the long-term effect of the weather aging will gradually lead to the performance degradation of the asphalt mastic, which is the binding part in the asphalt pavement. It will cause irreversible performance degradation and reduce the service life. Therefore, it is necessary to study the influence of aging on the long-term performance of asphalt mastic [1]. During the service of asphalt pavement, many factors could lead to the aging of asphalt mastic, such as sunlight, oxygen, and moisture. These factors react with the asphalt mastic and then resulting in aging, which include thermal oxidative aging, photo-oxidative aging, and water-damage aging.
For thermo-oxygen aging, thermal oxidation aging occurs during asphalt pavement construction including asphalt mixture mixing, laying and compaction. Thermal oxidation aging is mainly due to thermal oxidation and light component loss of asphalt binders at high temperatures [2,3]. Thermal oxidation aging is usually simulated using the rolling thin film oven test (RTFOT) based on ASTM D2872-04 and pressure aging vessel test (PAV) based on ASTM D6521. For photo-oxygen aging, as early as 1961, Traxler found that light could accelerate the aging of asphalt pavement materials [4].
Further study showed that the ultraviolet radiation will break the energy-rich bond, such as C-C and C-H, of the polymer chain in the asphalt. A series of photo-oxidative reactions will occur, thus changing the composition and structure of the asphalt and deteriorating its mechanical properties [5,6]. It is general to adopt the ultraviolet lamp or the high-pressure mercury lamp as the light source. In addition, the xenon lamp is also used as the light source because its spectral lines are similar with those of sunlight [7,8]. The radiation intensity and light duration are generally set according to the actual solar radiation in the target region. For water-damage aging, indoor simulation of water damage aging methods include boiling method, immersion method, and freezing and thawing cycle method. By these methods, researchers found that water damage of asphalt binder mainly exists in two aspects. One is that the water diffuses and permeates in the asphalt binder, and reacts with water-sensitive components and polar groups. This leads to the gradual deterioration of the asphalt [9][10][11]. The other is that the deterioration of the asphalt performance attenuates the adhesion between the binder and the aggregates [12][13][14][15].
However, the long-term effects of light, oxygen, humidity and heat jointly lead to the aging of asphalt pavement materials and degradation of pavement performance. In order to simulate the actual long-term aging process of asphalt pavement materials, it is necessary to consider the multi-factor coupling aging test of asphalt and asphalt mixture, as the climate aging mode based on the actual environment.
Additionally, asphalt mastic is a kind of mixture of asphalt and mineral powder with a particle size less than 0.075 mm, which binds coarse and fine aggregates into a whole, so it plays a crucial role in the macroscopic properties of asphalt concrete (such as viscoelasticity and fatigue) [16]. In the early years, based on the dynamic shear rheometer, the fatigue factor of asphalt binder at 1% strain level and 10 Hz loading frequency was used as the evaluation index of asphalt fatigue performance [17]. Later, time scanning tests for controlling stress or strain were widely used [18]. In recent years, Johnson et al. proposed linear amplitude scanning test (LAS) for asphalt cement based on viscoelastic continuous damage theory (VECD) to solve the defect of time-consuming time scanning test [19,20]. Due to the differences between asphalt and asphalt mastic, specific test scheme setting and evaluation index selection of asphalt mortar still need to be modified.
Therefore, based on the specific climate of northeast China, the long-term effects of multi-factors, such as sun-light, oxygen, and water were considered comprehensively. The base asphalt mastic and SBS modified asphalt mastic were selected as the object, which were made of limestone powder mixed with 70# base asphalt and SBS modified asphalt, respectively. An indoor accelerated weather aging test method was designed, which could simulate the environmental conditions of the seasonal frozen region in northeast China. The viscoelasticity and fatigue performance of asphalt mastic with different aging degrees were studied. It is attempted to provide theoretical and experimental basis for optimizing the selection of asphalt pavement materials and the determination of filler-asphalt ratio in the seasonal frozen region of northeast China.

Raw Materials
The raw material used in this paper was Pen 70 matrix asphalt produced by Inner Mongolia Luant Asphalt Hi-Tech Co.,Wuxi, China SBS modified asphalt was the commercial product with PG70-28 that supplied by Hubei Guochuang Hi-Tech Materials Co, Wuhan, China. In accordance with the testing methods in the "Standard Test Methods of Bitumen and Bituminous Mixtures for Highway Engineering" (JTG E20-2011), several technical indicators of base asphalt, modified asphalt and limestone powder were tested. The testing results are shown in Tables 1 and 2. All the indicators meet the technical requirements in the "Technical Specification for Construction of Highway Asphalt Pavements" (JTG F40-2004). Limestone powder was obtained by grinding Guilin limestone finely with a Materials 2021, 14, 6163 3 of 22 density of 2.750 g/cm 3 and a specific surface area of 0.582. X-ray fluorescence spectroscopy (Quant X type, Thermo Fisher Scientific Inc, Massachusetts, USA.) was used to analyze the raw material raw samples of limestone powder for compositional analysis, and the results are shown in the Table 3.

Preparation of Asphalt Mastic
Five kinds of base asphalt mastic and SBS modified asphalt mastic were prepared with the filler-asphalt ratio of 0.6, 0.8, 1.0, 1.2, and 1.4. The preparation process was as follows: (1) The asphalt was heated to a fluid state in the oven. The heating temperature of the base asphalt and SBS modified asphalt is 140 • C and 170 • C, respectively. A certain quality of the asphalt required for further tests was poured into the stainless steel cups; (2) The limestone powder was heated in the oven at 140 • C for 2 h to make itself dry and the internal temperature uniform. The required quality of the limestone powder depended on the asphalt quality and the filler-asphalt ratio; (3) The stainless steel cups containing the asphalt were heated on the electric stove to keep the asphalt in a fluid state during the mixing process. The mixing temperature was maintained at 140 • C ± 5 • C (base asphalt), 165 • C ± 5 • C (SBS modified asphalt); (4) The limestone powder was added to the asphalt in batches, and mixed with the high-speed shearing machine. The speed was controlled at 1000 r/min and each mixing time was 3 min. After all the limestone powder was added, the asphalt mastic was mixed for another 5 min to ensure the even distribution of the limestone powder in the mixture, as shown in Figure 1; (5) The evenly mixed asphalt mastic was poured into the silastic mold, making two kinds of samples with a thickness of 2 mm, a diameter of 8 mm and a thickness of 1 mm, a diameter of 25 mm. The asphalt mastic was also poured into the weather aging tray with 50 g per tray for short-term aging and weather aging tests.  (5) The evenly mixed asphalt mastic was poured into the silastic mold, making two kinds of samples with a thickness of 2 mm, a diameter of 8 mm and a thickness of 1 mm, a diameter of 25 mm. The asphalt mastic was also poured into the weather aging tray with 50 g per tray for short-term aging and weather aging tests.

Weather Aging Methods
In order to realize the aging effect of climate factors, such as sunlight, oxygen, temperature, and precipitation on asphalt mixture, this study adopted the weather aging box (Figure 2), referred to the artificial weathering aging method of vulcanized rubber in the standard GB/T 15255-2015 [21] and the accelerated aging test method for asphalt in ASTM D4798 [22], and comprehensively considered the natural environmental conditions in the seasonal frozen region of Northeast China. The conditions of the indoor accelerated weather aging test were as follows: (1) The irradiation intensity of the xenon lamp was 600 W; (2) The temperature of the black sheet was (65 ± 2.5) °C; (3) The relative humidity was controlled at (65 ± 5) %; (4) The rainfall intensity was 25 mm/h. An aging cycle was 1 h, including 54 min of light and 6 min of rainfall. Moreover, the accelerated rate of asphalt aging was determined to be 3.84 by using a UV intensity meter to measure the UV radiation intensity of 135.6 W/m 2 outdoors and 520 W/m 2 inside the aging chamber, respectively. Meanwhile, taking into account the effects of water, temperature and airflow rate in the experiment, a correction factor 1.25 times was chosen to correct the accelerated rate of asphalt aging to 4.8. The average annual sunshine in the northeastern region is 485.5 h. Combining this with the accelerated rate, it was found that 100 h of simulated aging in the laboratory is equivalent to 1 year of filed aging of the asphalt, and in this experiment the asphalt was aged for 100, 200, and 300 h to characterizes 1, 2, and 3 years of aging, and the specimens were turned over every 10 h.
The asphalt mastic with different weather aging time was shown in Figure 3. It could be seen that as the aging time increased, a small amount of mineral powder precipitated in the asphalt mastic with a high filler-asphalt ratio. Generally, the precipitation phenomenon in the base asphalt mastic was more obvious. It means that in the environment of light, oxygen, and water, the stability of the asphalt mastic is affected by the content of the mineral powder. With the increase in the filler-asphalt ratio, the stability becomes worse.

Weather Aging Methods
In order to realize the aging effect of climate factors, such as sunlight, oxygen, temperature, and precipitation on asphalt mixture, this study adopted the weather aging box (Figure 2), referred to the artificial weathering aging method of vulcanized rubber in the standard GB/T 15255-2015 [21] and the accelerated aging test method for asphalt in ASTM D4798 [22], and comprehensively considered the natural environmental conditions in the seasonal frozen region of Northeast China. The conditions of the indoor accelerated weather aging test were as follows: (1) The irradiation intensity of the xenon lamp was 600 W; (2) The temperature of the black sheet was (65 ± 2.5) • C; (3) The relative humidity was controlled at (65 ± 5) %; (4) The rainfall intensity was 25 mm/h. An aging cycle was 1 h, including 54 min of light and 6 min of rainfall. Moreover, the accelerated rate of asphalt aging was determined to be 3.84 by using a UV intensity meter to measure the UV radiation intensity of 135.6 W/m 2 outdoors and 520 W/m 2 inside the aging chamber, respectively. Meanwhile, taking into account the effects of water, temperature and airflow rate in the experiment, a correction factor 1.25 times was chosen to correct the accelerated rate of asphalt aging to 4.8. The average annual sunshine in the north-eastern region is 485.5 h. Combining this with the accelerated rate, it was found that 100 h of simulated aging in the laboratory is equivalent to 1 year of filed aging of the asphalt, and in this experiment the asphalt was aged for 100, 200, and 300 h to characterizes 1, 2, and 3 years of aging, and the specimens were turned over every 10 h.
The asphalt mastic with different weather aging time was shown in Figure 3. It could be seen that as the aging time increased, a small amount of mineral powder precipitated in the asphalt mastic with a high filler-asphalt ratio. Generally, the precipitation phenomenon in the base asphalt mastic was more obvious. It means that in the environment of light, oxygen, and water, the stability of the asphalt mastic is affected by the content of the mineral powder. With the increase in the filler-asphalt ratio, the stability becomes worse.

Dynamic Frequency Sweep Test
When the asphalt mastic is in the range of linear viscoelasticity, the stress-strain relationship is independent of the test temperature and loading frequency. The master curves of the viscoelastic parameters can be determined according to the time-temperature equivalence principle. Although in the range of non-linear viscoelasticity, the principle is no longer applicable. So, the dynamic rheological tests of asphalt mastic should be carried out within the range of linear viscoelasticity. Firstly, strain sweep is performed on the asphalt mastic to determine the linear viscoelasticity range. With the increase in the

Dynamic Frequency Sweep Test
When the asphalt mastic is in the range of linear viscoelasticity, the stress-strain relationship is independent of the test temperature and loading frequency. The master curves of the viscoelastic parameters can be determined according to the time-temperature equivalence principle. Although in the range of non-linear viscoelasticity, the principle is no longer applicable. So, the dynamic rheological tests of asphalt mastic should be carried out within the range of linear viscoelasticity. Firstly, strain sweep is performed on the asphalt mastic to determine the linear viscoelasticity range. With the increase in the

Dynamic Frequency Sweep Test
When the asphalt mastic is in the range of linear viscoelasticity, the stress-strain relationship is independent of the test temperature and loading frequency. The master curves of the viscoelastic parameters can be determined according to the time-temperature equivalence principle. Although in the range of non-linear viscoelasticity, the principle is no longer applicable. So, the dynamic rheological tests of asphalt mastic should be carried out within the range of linear viscoelasticity. Firstly, strain sweep is performed on the asphalt mastic to determine the linear viscoelasticity range. With the increase in the mineral powder content in the asphalt mastic, the structural asphalt increases, and the thickness of the asphalt film on the surface of the mineral powder decreases. The volume filling and hardening effect of the mineral powder on the asphalt mastic is enhanced, which macroscopically shows that the linear viscoelasticity range of the asphalt mastic decreases [23]. Meanwhile, the effect of weather aging will harden the asphalt mastic and reduce the linear viscoelasticity range. Therefore, base asphalt mastic and SBS modified asphalt mastic with the filler-asphalt ratio of 1.4 and aged for 400 h were selected for the strain sweep test. Figure 4 showed the sweep results of base asphalt mastic and SBS modified asphalt mastic at 20 • C and 10 Hz. According to the SHRP definition of the ultimate linear viscoelastic strain of asphalt materials, the linear viscoelastic limit is the strain at which the maximum complex shear modulus reduces to 90%. As shown in Figure 4, the complex shear modulus of base and modified asphalt mastic reduce to 90%Gmax were at 0.12% and 0.11%, respectively. Therefore, the linear viscoelastic strain range of the base asphalt mastic is 0-0.12% under the test conditions of 20 • C and 10 Hz, and that of the SBS modified asphalt mastic is 0-0.11%.
Materials 2021, 14, x FOR PEER REVIEW 6 of 23 mineral powder content in the asphalt mastic, the structural asphalt increases, and the thickness of the asphalt film on the surface of the mineral powder decreases. The volume filling and hardening effect of the mineral powder on the asphalt mastic is enhanced, which macroscopically shows that the linear viscoelasticity range of the asphalt mastic decreases [23]. Meanwhile, the effect of weather aging will harden the asphalt mastic and reduce the linear viscoelasticity range. Therefore, base asphalt mastic and SBS modified asphalt mastic with the filler-asphalt ratio of 1.4 and aged for 400 h were selected for the strain sweep test. Figure 4 showed the sweep results of base asphalt mastic and SBS modified asphalt mastic at 20 °C and 10 Hz. According to the SHRP definition of the ultimate linear viscoelastic strain of asphalt materials, the linear viscoelastic limit is the strain at which the maximum complex shear modulus reduces to 90%. As shown in Figure 4, the complex shear modulus of base and modified asphalt mastic reduce to 90%Gmax were at 0.12% and 0.11%, respectively. Therefore, the linear viscoelastic strain range of the base asphalt mastic is 0-0.12% under the test conditions of 20 °C and 10 Hz, and that of the SBS modified asphalt mastic is 0-0.11%.
(a) (b) Based on the linear viscoelasticity range determined by the strain sweep test, the strain level of the frequency sweep test is set to 0.01%, and the loading frequency is 0.2 Hz, 0.4 Hz, 0.6 Hz, 0.8 Hz, 1 Hz, 2 Hz, 4 Hz, 6 Hz, 8 Hz, 10 Hz, 20 Hz, and 30 Hz. The test temperature is 5 °C, 15 °C, 25 °C, and 35 °C. Within the range of 10 °C to 30 °C, samples with a diameter of 8 mm and a height of 2 mm are used. When the test temperature is above 30 °C, samples with a diameter of 25 mm and a height of 1 mm are used. In order to analyze the viscoelastic characteristics of the asphalt mastic in a wider range of temperature and frequency, the WLF equation (Equation (1)) is adopted on the basis of the timetemperature equivalence principle. Taking 20 °C as the reference temperature, the shift factor is calculated. The complex shear modulus and phase angle at other temperatures and frequencies are shifted to the reference curve. Thus, the master curves of the complex shear modulus and phase angle of the asphalt mastic are obtained.
where, T α is the shift factor at the reference temperature; T is the test temperature; T0 is the reference temperature; C1 and C2 are the fitting parameters.
Referring to the existing studies, the Sigmoidal, CA, and CAM models have relatively high fitting accuracy to the master curves of the complex shear modulus of the asphalt  Based on the linear viscoelasticity range determined by the strain sweep test, the strain level of the frequency sweep test is set to 0.01%, and the loading frequency is 0.2 Hz, 0.4 Hz, 0.6 Hz, 0.8 Hz, 1 Hz, 2 Hz, 4 Hz, 6 Hz, 8 Hz, 10 Hz, 20 Hz, and 30 Hz. The test temperature is 5 • C, 15 • C, 25 • C, and 35 • C. Within the range of 10 • C to 30 • C, samples with a diameter of 8 mm and a height of 2 mm are used. When the test temperature is above 30 • C, samples with a diameter of 25 mm and a height of 1 mm are used. In order to analyze the viscoelastic characteristics of the asphalt mastic in a wider range of temperature and frequency, the WLF equation (Equation (1)) is adopted on the basis of the time-temperature equivalence principle. Taking 20 • C as the reference temperature, the shift factor is calculated. The complex shear modulus and phase angle at other temperatures and frequencies are shifted to the reference curve. Thus, the master curves of the complex shear modulus and phase angle of the asphalt mastic are obtained.
where, α T is the shift factor at the reference temperature; T is the test temperature; T 0 is the reference temperature; C 1 and C 2 are the fitting parameters.
Referring to the existing studies, the Sigmoidal, CA, and CAM models have relatively high fitting accuracy to the master curves of the complex shear modulus of the asphalt mastic [21]. Due to the clear physical meaning, the CAM model is adopted to fit the master curve of complex shear modulus, as shown in Equation (2).
where, G g is the glassy shear modulus; ω c is the crossover frequency; ω r is the reduced frequency; v = log 2 R , R is the ratio of G g and |G * | at the crossover frequency; v is the parameter determining the shape of master curve.
Since the deviation in the low frequency range is large when the CAM model is used to fit the phase angle master curve of the asphalt mastic, the double-logistic model (Equation (3)) is adopted to fit the master curve of the phase angle [24]. The shift factor is the same as that of the complex shear modulus master curve.
where, δ P is the plateau phase angle; f red is the reduced frequency; f P is the frequency at the plateau; H( f red − f P ) and H( f P − f red ) are two heaviside step functions; S R stands for the rise (or fall) on the right side and δ L for the left side of plateau δ P .

Fatigue Test
A new test based on the theory of viscoelastic continuum damage is adopted. According to the "Estimating damage tolerance of asphalt binders using the linear amplitude sweep" (AASHTO TP101-14), the test consists of two parts. Firstly, frequency sweep is carried out on the asphalt mastic at the designed temperature. The strain is 0.01% and the sweep frequency is 0.2 Hz, 0.4 Hz, 0.6 Hz, 0.8 Hz, 1 Hz, 2 Hz, 4 Hz, 6 Hz, 8 Hz, 10 Hz, 20 Hz, and 30 Hz. Based on complex shear modulus and phase angle, the undamaged material parameter α in the further fatigue damage analysis is calculated. Secondly, at a loading frequency of 10 Hz, the asphalt mastic is swept with a linear increase in strain amplitude from 0.1% to 30% within 310 s. The linear amplitude sweep test at 15 • C was performed on the asphalt mastic. It was found that the adhesive failure occurred quickly between the asphalt mastic and the upper and lower rotors. Therefore, the temperature of the linear amplitude sweep test in this study is set to 20 • C.

S-VECD Model
The establishment of the simplified viscoelastic continuum damage (S-VECD) model is based on three principles.
(1) Elastic-viscoelastic correspondence principle based on the pseudo strain. The viscoelastic problem can be simplified into the elastic problem with this principle; (2) Work potential theory based on the continuous damage mechanics, and it can quantify the influence of internal microdamage on the macroscopic mechanical behavior; (3) Time-temperature equivalence principle. It includes the influence of loading time and temperature on the physical and mechanical behavior of materials [25].
(1) Elastic-viscoelastic correspondence principle based on the pseudo strain; The linear viscoelastic materials are dependent on time and temperature. The machinal response is not only related to the current load, but also the entire loading history. The constitutive relation of viscoelastic materials can be represented by Equations (4) and (5).
where, E(t) and D(t) are, respectively, the relaxation modulus and creep compliance of materials; τ is the integration variable. The Schapery elastic-viscoelastic correspondence principle replaces the viscoelastic strain of materials with the pseudo strain. The viscoelastic problem is transformed into the elastic problem, as shown in Equation (6).
where, ε R is the pseudo strain; ε is the physical strain of materials; E R is the specific reference modulus. The constitutive relation of viscoelastic materials based on the pseudo strain can be obtained by bringing Equation (4) into Equation (6), as shown in Formula (7). For viscoelastic materials, the modulus attenuation under load is due to the combined effect of viscoelastic behavior and internal damage. The constitutive relation based on the pseudo strain can separate the viscoelastic behavior and the damage, thus effectively analyzing the damage development of materials under load.
(2) Work potential theory based on the continuous damage mechanics; Since it is difficult to define and quantify the microdamage inside materials, the continuous damage mechanics describes the damage evolution process of materials by reducing the effective stiffness resulted from the microstructure change. The VECD model adopts Schapery work potential theory to quantify the damage. It defines the damage as the internal state variable of materials, and quantitatively analyzes the fatigue damage process.
The work potential principle of viscoelastic materials is shown in Equation (8).
where, W R is the pseudo strain energy density; S is the internal state variable of materials characterizing the damage. Equation (9) expresses the damage rate of materials.
where, α is the material parameter under undamaged state, α = 1 m ; m is the fitted slope of composite shear modulus master curve of materials in linear viscoelasticity range.
The pseudo strain energy density of materials can be calculated by the following Equation (10).
where, C(S) is the pseudo modulus, which is the function of material damage parameter S. It can be presented by Equations (11)- (13).
where, τ P is the peak shear stress in any loading cycle; ε R P is peak pseudo shear strain in any loading cycle; dynamic modulus ratio (DMR) is the specimen variability compensation parameter, which generally needs to be limited with 0.9~1.1; |G * | 0 is the dynamic modulus at the simulated temperature and frequency; |G * | I is the initial dynamic modulus.
The expression of damage variable S can be derived from Equations (9)- (12).
where, i is the time interval selected for damage calculation.
The characteristic curve formula of material damage can be obtained by using power function to fit the pseudo modulus C and damage variable S of materials, as shown in Equation (15).
where, C 1 and C 2 are the fitting parameters.
where, f is the loading frequency. Substituting the cumulative damage variable S f at fatigue failure into Equation (16), the number of cycles at failure N f can be estimated by the Equation (17):

Fatigue Failure Criterion of Asphalt Mastic
The fatigue failure of materials should be determined according to the internal damage distribution. Existing methods make it difficult to obtain the physical condition change of material interior during the fatigue process. In general, the fatigue failure point of materials is defined based on the phenomenological method or energy method. In the time sweep test of asphalt, the main measures determining the fatigue life are as follows.
(1) Cycle number when the dynamic modulus is reduced to 50% of the initial modulus (N f50 ); (2) Cycle number at the peak of |G * | × i, which i is the loading times (N NM ); (3) Cycle number at the inflection point of the change rate of |G * |·(N G * ); (4) Cycle number at the peak of energy dissipation rate (N DER ); (5) The fatigue life is determined on the basis of the cumulative dissipation energy ratio (CDER). Specifically, the relationship between CDER and the loading cycle is linear at the beginning. When the relation curve between CDER and the loading cycle deviates from the straight line by 20%, the material is considered to have fatigue failure. The cycle number at this time is defined as the fatigue life.
In the researches on linear amplitude sweep (LAS) test of asphalt, a variety of measures can determine the material fatigue failure based on the phenomenological method or the energy method. AASHTO TP101 (2012) took the shear stress peak point as the fatigue failure point of asphalt. Although in the new version (2018), the fatigue failure was defined at the point when |G * | sin δ decreases to 65% of the initial value. Wang believed that although phenomenological indicators, such as peak shear stress and peak phase angle, can be used to define the fatigue failure of asphalt binders, it was necessary to establish a fatigue failure index based on the energy method, which could be used for VECD model prediction. The index related to pseudo strain energy was proposed [26]. Under the effect of load, the pseudo strain energy of viscoelastic materials W R t includes the stored pseudo strain energy W R S and released pseudo strain energy W R r , which are presented in Equations (17)- (20). Figure 5 shows the W R S and W R r curves of asphalt mastic in the LAS test. It can be seen that W R r increases with the increase in the pseudo strain, and W R S drops rapidly after reaching the peak. With the initiation and development of damage, the slope of W R r curve increases gradually. It indicates that the energy storage capacity of materials decreases as the damage becomes worse. Therefore, the peak value of W R S represents the maximum energy that the material can store, which can be used as the definition of asphalt fatigue failure. Recently, the fatigue failure criterion of asphalt binder based on W R S has been applied to some extent, which proves its applicability in the asphalt fatigue behavior [27].  Figure 6 was the curves of complex shear modulus, shear stress, and storage pseudo strain energy for base asphalt mastic and SBS modified asphalt mastic. As could be seen, when the complex shear modulus dropped to 50% of the initial value, the shear stress, and stored pseudo strain energy of the asphalt mastic had not reached the peak value. It is conservative to use W . It means that at the beginning of the shear stress attenuation, the energy storage capacity of the asphalt mastic is still in the enhancement stage. The material still has a strong carrying capacity. Therefore, compared with indexes related to shear stress and complex shear modulus, it is more suitable to use the maximum value of the stored pseudo strain energy to determine the fatigue failure of the asphalt mastic.  Figure 6 was the curves of complex shear modulus, shear stress, and storage pseudo strain energy for base asphalt mastic and SBS modified asphalt mastic. As could be seen, when the complex shear modulus dropped to 50% of the initial value, the shear stress, and stored pseudo strain energy of the asphalt mastic had not reached the peak value. It is conservative to use 0.5|G * | 0 as the fatigue failure index. If the inflection point of the complex shear modulus curve is taken as the fatigue failure point, the time at the inflection point of base asphalt mastic is basically the same as the time when W R S reaches the maximum value (Figure 6a). Although, for SBS modified asphalt mastic, the complex shear modulus curve is still in a stable attenuation stage when W R S reaches the maximum value (Figure 6b). The inflection point fails to occur by the end of the test. At this time, the attenuation amplitude of the shear stress curve is large and the obvious fatigue failure occurs. It is unreasonable to use the inflection point of the complex shear modulus curve as the fatigue failure index of SBS modified asphalt mastic. For two kinds of asphalt mastic, the peak of the shear stress occurs slightly earlier than that of W R S . It means that at the beginning of the shear stress attenuation, the energy storage capacity of the asphalt mastic is still in the enhancement stage. The material still has a strong carrying capacity. Therefore, compared with indexes related to shear stress and complex shear modulus, it is more suitable to use the maximum value of the stored pseudo strain energy to determine the fatigue failure of the asphalt mastic. tic, the peak of the shear stress occurs slightly earlier than that of R s W . It means that at the beginning of the shear stress attenuation, the energy storage capacity of the asphalt mastic is still in the enhancement stage. The material still has a strong carrying capacity. Therefore, compared with indexes related to shear stress and complex shear modulus, it is more suitable to use the maximum value of the stored pseudo strain energy to determine the fatigue failure of the asphalt mastic.

Influence of Weather Aging on Viscoelasticity Performance
Taking the base asphalt mastic and SBS modified asphalt mastic with the filler-asphalt ratio of 0.6 and 1.4 as examples, the influence of weather aging on the viscoelasticity performance of the asphalt mastic was studied. Figure 7 showed the complex shear modulus master curves of asphalt mastic with different aging degrees. It was observed that with the increment of aging time, the complex shear modulus of the two kinds of asphalt mastic tended to increase. Under the condition of low-frequency loading, the complex shear modulus was sensitive to the effect of weather aging. The increment was significant. Although under high-frequency loading, the complex shear modulus tended to the same asymptotic line, less sensitive to weather aging. It is consist with the results of Liang et al., which studied the master curve of modified asphalt with different aging time and found [28]. The combined aging effect of light, oxygen and water increased the stiffness of the asphalt mastic. Due to the time limit of the aging test, the aging cycle interval was small. Thus, the variation of the complex shear modulus master curve with the aging time was relatively small. Figure 8 was the phase angle master curve of base asphalt mastic and SBS modified asphalt mastic with different aging degrees. As the aging degree deepened, the phase angle of base asphalt mastic dropped. The elastic component increased while the viscous component decreased. For SBS modified asphalt mastic, the peak value of phase angle increased after short-term aging. The loading frequency corresponding to the peak phase angle decreased and the viscous component increased. In the weather aging stage, the peak value of phase angle dropped with the increment of aging time. According to the description of the aging mechanism of SBS modified bitumen, there are two parts of SBS modified bitumen aging, the degradation of the modifier and the change of the chemical composition of the asphalt [29]. Therefore, after conducted the short-term aging, SBS modifier was partially broken and decomposed into small molecules due to thermal oxidative aging and photo-oxidative aging. In this stage, the decomposition effect of SBS modifier is greater than the hardening effect caused by the oxidative dimerization of base asphalt. It increases the viscosity of SBS modified asphalt and the peak value of phase angle. During the stage of weather aging, with the continuity of thermal oxidative aging and photo-oxidative aging, the decomposition products and the small molecules in base asphalt polymerize to macromolecules. This gradually hardens SBS modified asphalt and decreases the phase angle. The elasticity of asphalt mastic increases and the viscosity drops.
Although under high-frequency loading, the complex shear modulus tended to the same asymptotic line, less sensitive to weather aging. It is consist with the results of Liang et al., which studied the master curve of modified asphalt with different aging time and found [28]. The combined aging effect of light, oxygen and water increased the stiffness of the asphalt mastic. Due to the time limit of the aging test, the aging cycle interval was small. Thus, the variation of the complex shear modulus master curve with the aging time was relatively small.  aging and photo-oxidative aging. In this stage, the decomposition effect of SBS modifier is greater than the hardening effect caused by the oxidative dimerization of base asphalt. It increases the viscosity of SBS modified asphalt and the peak value of phase angle. During the stage of weather aging, with the continuity of thermal oxidative aging and photooxidative aging, the decomposition products and the small molecules in base asphalt polymerize to macromolecules. This gradually hardens SBS modified asphalt and decreases the phase angle. The elasticity of asphalt mastic increases and the viscosity drops.

Influence of Weather Aging on Fatigue Performance
The influence of weather aging on the fatigue performance of base asphalt mastic and SBS modified asphalt mastic was analyzed in this section. It should be noted that after aging for 400 h, the adhesive failure occurred quickly between the rotors and the base asphalt mastic with the filler-asphalt ratio of 1.2 and 1.4. This demonstrates that the adhesion performance of the base asphalt mastic with deep aging degree is reduced. The interface cracking is more likely to occur.
The stress-strain curves of base asphalt mastic and SBS modified asphalt mastic with different aging degrees were presented in Figure 9. It showed that as the aging degree deepened, the peak shear stress of asphalt mastic increased. Comparing Figure 9a,b and Figure 9c,d, it could be seen that for the asphalt mastic with different mineral powder content, the shear strain corresponding to the peak shear stress varied differently with the aging degree. When the filler-asphalt ratio of asphalt mastic was 0.6, the shear strain corresponding to the peak shear stress was basically the same under different aging degrees. The shear strain when the fatigue failure occurred was also consistent. It indicated that the weather aging had less influence on the deformability of the asphalt mastic with a low filler-asphalt ratio. When the filler-asphalt ratio was 1.4, the shear strain corresponding to the peak shear stress decreased gradually with the deepening of aging degree. As the aging degree deepened, the brittle failure characteristics became more obvious and the failure strain was smaller. So, the deformability of the asphalt mastic with a high filler-asphalt ratio was greatly affected by weather aging. Figure 9c,d, it could be seen that for the asphalt mastic with different mineral powder content, the shear strain corresponding to the peak shear stress varied differently with the aging degree. When the filler-asphalt ratio of asphalt mastic was 0.6, the shear strain corresponding to the peak shear stress was basically the same under different aging degrees. The shear strain when the fatigue failure occurred was also consistent. It indicated that the weather aging had less influence on the deformability of the asphalt mastic with a low filler-asphalt ratio. When the filler-asphalt ratio was 1.4, the shear strain corresponding to the peak shear stress decreased gradually with the deepening of aging degree. As the aging degree deepened, the brittle failure characteristics became more obvious and the failure strain was smaller. So, the deformability of the asphalt mastic with a high fillerasphalt ratio was greatly affected by weather aging.  Figure 10 showed the C S − curves of base asphalt mastic and SBS modified asphalt mastic with different aging degrees. As can be seen in Figure 10a,b, within 300 h before weather aging, the C S − curves of base asphalt mastic with the filler-asphalt ratio of 0.6 and 1.4 rose gradually as the aging time increased. Under the same damage variable S, the material integrity of asphalt mastic was higher with the longer aging time. After 400 h of weather aging, the material integrity attenuation rate C of JZ0.6 asphalt mastic increased significantly in the initial damage state when S was less than 200. Once S was greater than 200, the attenuation rate C of JZ0.6 asphalt mastic was lower than that of asphalt mastic with low aging degree. In Figure 10c, within 300 h before weather aging, the C S − curves of SBS0.6 asphalt mastic showed a gradual rise trend. After 400 h of aging, the curves dropped in the late damage state. In Figure 10d, within 300 h before aging, the C S − curves of SBS1.4 asphalt mastic rose gradually with the increase in aging time. After being aged for 400 h, the curves dropped in the initial damage state when S was less than 250. The attenuation rate of material integrity increased. When S was greater than 250, the curves rose relatively and the material integrity attenuation rate decreased. It displayed significant damage resistance ability in the later stage.  Figure 10 showed the C − S curves of base asphalt mastic and SBS modified asphalt mastic with different aging degrees. As can be seen in Figure 10a,b, within 300 h before weather aging, the C − S curves of base asphalt mastic with the filler-asphalt ratio of 0.6 and 1.4 rose gradually as the aging time increased. Under the same damage variable S, the material integrity of asphalt mastic was higher with the longer aging time. After 400 h of weather aging, the material integrity attenuation rate C of JZ0.6 asphalt mastic increased significantly in the initial damage state when S was less than 200. Once S was greater than 200, the attenuation rate C of JZ0.6 asphalt mastic was lower than that of asphalt mastic with low aging degree. In Figure 10c, within 300 h before weather aging, the C − S curves of SBS0.6 asphalt mastic showed a gradual rise trend. After 400 h of aging, the curves dropped in the late damage state. In Figure 10d, within 300 h before aging, the C − S curves of SBS1.4 asphalt mastic rose gradually with the increase in aging time. After being aged for 400 h, the curves dropped in the initial damage state when S was less than 250. The attenuation rate of material integrity increased. When S was greater than 250, the curves rose relatively and the material integrity attenuation rate decreased. It displayed significant damage resistance ability in the later stage.
greater than 200, the attenuation rate C of JZ0.6 asphalt mastic was lower than that of asphalt mastic with low aging degree. In Figure 10c, within 300 h before weather aging, the C S − curves of SBS0.6 asphalt mastic showed a gradual rise trend. After 400 h of aging, the curves dropped in the late damage state. In Figure 10d, within 300 h before aging, the C S − curves of SBS1.4 asphalt mastic rose gradually with the increase in aging time. After being aged for 400 h, the curves dropped in the initial damage state when S was less than 250. The attenuation rate of material integrity increased. When S was greater than 250, the curves rose relatively and the material integrity attenuation rate decreased. It displayed significant damage resistance ability in the later stage.   The fatigue life Equation (17) can be represented in another way: The fatigue life equation at the undamaged interface is: The fatigue life Equation (17) can be represented in another way: where, A = f ·2 α ·S 1−αC 2 +α f (1−αC 2 +α)(C 1 C 2 ) α (|G * | 0 ) 2α and B = −2α. The fatigue life equation at the undamaged interface is: The fatigue life equation at the damaged interface is: Let N f h = N f 0 , the calculation equation of the cross fatigue strain is obtained: When the strain load is γ c , the fatigue life of the aged and unaged asphalt mastic is the same. When the strain load level is lower than γ c , the fatigue life of the aged asphalt mastic is higher than that of the unaged asphalt mastic. It means that under low strain load, environment aging has a positive impact on the fatigue performance of asphalt mastic. When the strain load level is higher than γ c , the fatigue life of the aged asphalt mastic is lower than that of the unaged asphalt mastic. It indicates that under high strain load, the influence of environmental aging on the fatigue performance of asphalt mastic is negative. According to the viscoelastic characteristic analysis of the asphalt mastic, the complex shear modulus of the aged asphalt mastic is relatively high. It shows high stiffness in terms of mechanical behavior. Under low strain load, the internal damage of the asphalt mastic is slight and the fatigue life is long due to the high strength. The deformation and bearing capacity of the aged asphalt mastic is relatively poor. The high strain load leads to the great internal damage of the asphalt mastic and the attenuation of fatigue life. Table 4 summarized the γ c values of fatigue life curves for base asphalt mastic with different aging degrees. With the increment of the filler-asphalt ratio, γ c values of the aged and unaged asphalt mastic first increased and then decreased. The cross fatigue strain of the asphalt mastic with the filler-asphalt ratio of 0.6, 1.2, and 1.4 was relatively small. Under high strain load, weather aging had an adverse impact on the fatigue performance of these three kinds of asphalt mastic. The cross fatigue strain of the asphalt mastic with the filler-asphalt ratio of 0.8 and 1.0 was greater than 10%. Within the strain load range of 0.01% to 10%, the fatigue life of the two kinds of asphalt mastic increased after aging. The influence of weather aging on the fatigue life of the base asphalt mastic with middle filler-asphalt ratio was positive. Referring to Figure 12, the fatigue life of the unaged base asphalt mastic increased with the increase in filler-asphalt ratio. When the aging degree was the same, the aged base asphalt mastic with the filler-asphalt ratio of 1.0 had the maximum fatigue life under the strain load of 0.1%, 1%, and 5%. Therefore, based on the fatigue performance of the aged asphalt mastic, the best filler-asphalt ratio of the base asphalt mastic was 1.0.     Figure 13 depicted the fatigue life of SBS modified asphalt mastic with different aging degrees. It pointed out that as the aging degree deepened, the fatigue life of SBS modified asphalt mastic was reduced within the strain load range of 0.01% to 10%. The slope of the fatigue life curves increased, which meant that the fatigue life became more sensitive to the strain load. Compared with base asphalt mastic, the influence of weather aging on the strain load sensitivity of SBS modified asphalt mastic was less. Within the strain load range, the cross fatigue strain failed to occur. Figure 14 presented the fatigue life of SBS modified asphalt mastic under the strain load of 0.1%, 0.5%, and 5%. When the filler-asphalt is less than 1.2, the fatigue life of asphalt mastic increased with the increase in the value of filler-asphalt ratio, and then reduced when the filler-asphalt ratio is 1.4. It indicated SBS modified asphalt mastic when the filler-asphalt ratio is 1.2 possesses the maximum fatigue life regardless of the degree of aging. Therefore, the fatigue performance of SBS modified asphalt mastic with different aging degrees was optimal when the filler-asphalt ratio was 1.2.

Conclusions
(1) Comparing the curves of complex shear modulus, shear stress, and storage pseudo strain energy for matrix asphalt mastic and SBS modified asphalt mastic, it is more suitable to use the maximum value of the stored pseudo strain energy to determine the fatigue failure of the asphalt mastic than shear stress and complex shear modulus; (2) As the aging time increases, the complex shear modulus of asphalt mastic increases significantly under low-frequency loading. Although under high-frequency loading, the complex shear modulus tends to the same asymptotic line. With the deepening of aging degree, the phase angle of base asphalt mastic and SBS modified asphalt mastic drops. The elastic component increases while the viscous component decreases; (3) According to the relationship between of fatigue life and the filler-asphalt ratio, the best filler-asphalt ratio of base asphalt and SBS modified asphalt are 1.0 and 1.2, respectively; (4) The base asphalt mastic becomes more sensitive to the strain level after aging than SBS modified asphalt mastic. The fatigue life of base asphalt mastic increases under the low strain load and decreases under the high strain load, but the fatigue life of SBS modified asphalt mastic gradually decreases within the strain load range of 0.1% to 10% as the aging degree deepens.