Next Article in Journal
Micropipes in SiC Single Crystal Observed by Molten KOH Etching
Next Article in Special Issue
Fatigue Models for Airfield Concrete Pavement: Literature Review and Discussion
Previous Article in Journal
A Brief Photocatalytic Study of ZnO Containing Cerium towards Ibuprofen Degradation
Previous Article in Special Issue
Evaluation of the Aging of Styrene-Butadiene-Styrene Modified Asphalt Binder with Different Polymer Additives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Viscoelastic Damage Characteristics of Asphalt Mixtures Using Fractional Rheology

1
School of Transportation, Southeast University, Nanjing 211189, China
2
National Demonstration Center for Experimental Road and Traffic Engineering Education (Southeast University), Nanjing 211189, China
*
Author to whom correspondence should be addressed.
Submission received: 11 September 2021 / Revised: 2 October 2021 / Accepted: 4 October 2021 / Published: 8 October 2021
(This article belongs to the Special Issue Experimental Testing and Constitutive Modelling of Pavement Materials)

Abstract

:
The mechanical behavior of asphalt mixtures at high stress levels are characterized by non-linear viscoelasticity and damage evolution. A nonlinear damage constitutive model considering the existence of creep hardening and creep damage mechanisms in the entire creep process is proposed in this study by adopting the fractional rheology theory to characterize the three-stage creep process of mixtures. A series of uniaxial compressive creep tests under various stresses were conducted at different temperatures to verify the model. The results indicated that the model predictions were in good agreement with the creep tests. The relationship between the model parameters and applied stresses was established, and the stress range in which the mixture exhibited only creep consolidation was obtained. The damage to the asphalt mixture was initiated in the steady stage; however, it developed in the tertiary stage. A two-parameter Weibull distribution function was used to describe the evolution between the damage values and damage strains at different stress levels and temperatures. The correlation coefficients were greater than 0.99 at different temperatures, indicating that a unified damage evolution model could be established. Thus, the parameters of the unified model were related to material properties and temperature, independent of the stress levels applied to the mixtures.

1. Introduction

Asphalt mixtures are considered as heterogeneous pavement materials typically composed of asphalt, aggregate, admixtures, and air voids. Furthermore, asphalt mixtures are the primary materials used worldwide for constructing pavement layers because of their remarkable advantages, such as short construction periods, long service lives, recyclability, ease of maintenance, and comfortable driving. Although aggregates, which are a type of rigid material, account for approximately 90% of all components within the asphalt mixtures, such mixtures exhibit complex rheological properties. Asphalt mixtures exhibit both viscous (fluid) and elastic (solid) behavior [1,2,3], and their mechanical properties are highly dependent on temperature, loading stress level, loading time, and loading speed [4,5]. The theoretical and experimental studies on the mechanical behaviors of asphalt mixtures provide a crucial basis for designing pavement layers [6,7]. Creep behavior is amongst the mechanical characteristics of asphalt mixtures [8], and it is strongly related to the rutting formation of the pavement structure [3,9]. Therefore, it is important to establish a constitutive model that can accurately predict the creep behavior of asphalt mixtures.
The improvement of constitutive models is gaining increasing interest in the field of rheology in order to better describe the mechanical behaviors of complex rheological materials. Furthermore, several rheological models of asphalt mixtures were proposed in the past [10,11,12,13,14,15,16,17,18,19]. These models can be divided into integral constitutive and differential constitutive models [15,20]. The integral viscoelastic constitutive model, which considers the loading history and memory of materials, is developed based on methods such as the Boltzmann superposition principle or the Leaderman modified superposition principle. Compared with multiple integral models, single integral models have fewer material parameters and simpler forms, which makes them easier to use in practical engineering. Schapery’s model, which is widely applied to characterize the viscoelastic properties of asphalt mixtures, is a single integral constitutive model [21]. Im et al. [22] described the viscoelastic and viscoplastic behavior of mixtures by combining Schapery’s viscoelastic model and Perzyna viscoplasticity with a generalized Drucker–Prager yield surface.
Classical integer differential constitutive models, such as the Maxwell model, Kelvin model, Burgers model, and generalized models, comprise springs and dashpots in series or in parallel. The mathematical expressions of the differential constitutive models can correspond to the mechanical models that have advantages of intuitionism and simplicity. The classical integer derivative constitutive models are widely used to analyze the creep behavior of asphalt mixtures. The Burgers model provides a better description of mechanical behavior than the Maxwell and Kelvin–Voigt models [23]. However, Cheng et al. [24] reported that, compared to the Burgers model, the generalized Kelvin (GK) model and the generalized Maxwell (GM) model have a better performance in describing the creep performance of asphalt mixtures. Although the classical constitutive models can describe the viscoelastic properties of asphalt mixtures, they yield a large error in describing the initial stage of the creep performance [13]. With the development of viscoelasticity theory, the fractional derivative operator that can describe the historical memory dependence and spatially wide relevance of a material meets the research requirements for examining the complex viscoelastic, non-Newtonian fluid and the porous media mechanics [25]. Furthermore, the fractional calculus theory achieves great success in describing the rheological behavior of asphalt materials [13,14,15,26,27]. Lagos-Varas et al. [13] established a viscoelastic model using fractional-order derivatives, and unlike the Burgers model, the fractional order model can predict the elastic jump observed at the beginning of the creep modulus. Celauro et al. [28] demonstrated that the fractional-order Burgers model can accurately characterize the creep/recovery response of asphalt mixtures. In addition, the fractional derivative model has the advantage of exhibiting a simpler expression, utilizing fewer parameters, and providing more accurate results. Furthermore analysis revealed that a fractional derivative model comprising two elements of fractional derivatives in series with a Maxwell element can be used to model the rheological behavior of asphalt binders. This includes the dynamic viscoelastic behavior, static creep, and relaxation characteristics [27].
The entire creep process of asphalt mixtures is characterized into three stages: primary (decelerated stage), steady (stationary stage), and tertiary (accelerated stage). Furthermore, different deformation mechanisms occur throughout the creep process of asphalt mixtures at various stresses. At a lower stress level, the asphalt mixture is gradually compacted and the aggregates are gradually stacked to form a stable skeleton structure, which results in a “consolidation effect” [19]. This deformation mechanism is known as the creep hardening mechanism [18], and the asphalt mortar resists deformation together with the skeleton structure during the procedure. However, at higher stress levels, the stress in the asphalt mixture exceeds the frictional resistance between the aggregates and the bonding between the aggregate and asphalt mortar, which results in the fragmentation of the asphalt mixture and the destruction of the skeleton. This consequently accelerates the viscoelastic flow rate and creep failure of the asphalt mixture in the tertiary stage. Here, the deformation mechanism is called the damage softening mechanism [18]. Linear models (e.g., the ordinary Burgers model, generalized Kelvin model, and fractional-order Burgers model) cannot describe the tertiary stage of creep.
One potential approach is to add a plastic element into the classic differential models to construct viscoelastoplastic models to accurately describe the three-stage creep behavior of asphalt mixtures at high stresses. The Nishihara model, which adds a plastic element to the Burgers model, was developed to describe the tertiary stage of creep and was used in the mechanical response analysis of bituminous materials [19]. A nonlinear viscoelastic-plastic creep model for asphalt mixtures was developed using variable-order fractional calculus, where time-varying viscoplastic elements were used to describe the tertiary stage of creep performance [15]. In addition, damage occurs in the asphalt mixture during the creep process. Therefore, another feasible approach is coupling a continuum damage evolution law with a linear viscoelastic model. Sun et al. [29] developed a damage constitutive model on the basis of the generalized Burgers model and the Perzyna viscoplastic flow theory, combined with a modified Rabotnov damage theory, which can better reflect the three stages of deformation in creep testing and the hardening and softening effects in the constant strain rate compression testing. Zhang et al. [18] suggested that the characteristics of the three-stage permanent deformation were attributed to a competition between the damage softening and strain hardening effect, and they proposed a viscoelastoplastic damage mechanics model wherein the damage and hardening variables were introduced to modify the Burger’s model for describing the deformation of asphalt mixtures.
Most existing creep damage models are formed by the integral order difference constitutive model coupled directly with the damage factor. These models have the disadvantage of many material parameters and overly complex expressions, which are not conducive to their generalization in practical engineering. Furthermore, these models directly couple the strain caused by damage (called the damage strain) to the viscoelastic strain of the asphalt mixture, which is not conducive to studying the strain development mechanism and damage characteristics during the creep process. Therefore, this study aims to propose a creep damage model based on the fractional derivative operator that considers both the creep hardening mechanism and damage softening mechanism to characterize and analyze the creep process of asphalt mixtures at various stress levels and temperatures. Simultaneously, a unified damage evolution model is constructed for different stress levels, and the development of the strain and damage characteristics of the asphalt mixture throughout the creep process is analyzed. Finally, a method for the statistical quantification of the damage evolution of the asphalt mixture is proposed.
The paper is structured in the following way: a nonlinear viscoelastic creep damage model adopting the fractional rheology theory is developed in Part 2. The compression creep tests are conducted on AC-13 asphalt mixtures in Part 3. In Part 4, the experimental results are described and the proposed model is experimentally validated. Meanwhile, the model parameters and the creep damage evolution of the asphalt mixtures are analyzed in Part 4. Finally, several conclusions are drawn in Part 5.

2. Theoretical Background

2.1. Modelling of Creep Hardening Mechanism

In the classical integer-order differential viscoelastic model, the elastic mechanical behavior of the material is modeled with a Hook spring, and the viscous behavior is modeled with a Newton dashpot. The viscoelastic behavior of the asphalt mixtures is achieved by various combinations of Hook springs and Newton dashpots in series or in parallel. Gu et al. [26] adopted the Abel spring-pot element to describe the viscoelastic behavior based on the fractional order calculus theory for better describing the complex viscoelastic mechanical behavior of the materials. The fractional order viscoelastic element (Abel spring-pot element) is shown in Figure 1 and the stress–strain relation of the Abel spring-pot element is expressed as in [26]:
σ t = ξ D r ε t 0 r 1 ,
where σ, t, ε, ξ, and Dr denote the nominal stress, function variable, strain, material coefficient, and differential operator of the Riemann–Liouville fractional calculus [30], respectively, which is defined as:
D r f t = d r f t d t r = 1 Γ 1 r d d t 0 t f τ t τ r d τ 0 r 1 ,
where r denotes the fractional order of differentiation, Г(r) denotes the gamma function, and Γ r = 0 t r 1 e t d t . According to Equations (1) and (2), the spring-pot reduces to the Hook spring when r = 0, whereas it reduces to the Newton dashpot when r = 1. Thus, the spring-pot can be regarded as a type of fractional differential viscoelastic element.
The strain can be obtained from Equation (1) as:
ε t = 1 ξ D - r σ t ,
where D-r denotes the integral operator of the Riemann–Liouville fractional calculus [30], which is defined as:
D - r f t = d - r f t d t - r = 0 t t τ r - 1 Γ r f τ d τ 0 r 1 ,
In the case of the creep, under a constant stress level σ t = σ 0 H ( t ) , the creep strain can be obtained from Equations (3) and (4) as:
ε t = σ 0 t r ξ Γ r + 1 = σ 0 K t r   0 r 1 ,
where K is related to the coefficient ξ and the fractional order r of the material, and K = ξ Γ r + 1 . H(t) denotes the Heaviside step function, that is, H t = 1 t > 0 0 t 0 . According to Equation (5), K is related to the deformation of the material and referred to here as the deformation factor.
Find the derivative of Equation (5) with respect to t:
ε ˙ = d ε d t = σ 0 r K t r 1 0 r 1 .
Equation (6) indicates that the strain rate of the material, which decreases as time increases, and is consistent with the decrease in the strain rate at the decelerating stage of creep. When time t is sufficiently large, the strain rate converges to zero, which is consistent with the consolidation effect at the lower stress levels [19]. Figure 2 shows the Abel model used to describe the consolidation effect and the decelerated stage of creep. (The test data in the figure were randomly selected from the creep tests in Section 3 to illustrate the applicability of the model.) Figure 2 shows that the Abel model can characterize both the hardening effect and the primary stage of creep performance very effectively. Thus, the Abel spring-pot model can be applied to characterize the creep hardening mechanism of asphalt mixtures.
The load is kept constant ( σ t = σ 0 H ( t ) ) for an assigned time t ¯ , where the viscoelastic material enters the creep phase. When t > t ¯ , the load is removed and the material enters the recovery phase, and the stress is zero. The stress in the creep recovery phase can be expressed as:
σ t = σ 0 H t σ 0 H t t ¯ .
where t ¯ denotes the loading time.
The creep-recovery strain can be obtained from Equations (5) and (7) as:
ε t = σ 0 Κ t r H t t ¯ t t ¯ r   0 r 1 .
Equation (8) shows that for aviscoelastic materials, there are always irrecoverable strains in the material, no matter how much stress is applied. In other words, viscoelastic materials have permanent strains at any time after being stressed.

2.2. Modeling of Damage Softening Mechanism

There is a damage effect throughout the creep of the asphalt mixture, and the damage evolution results in a nonlinear mechanical behavior. Herein, the damage variable D is used to describe the deterioration of the Abel spring-pot element, and damage evolution is related to the stress levels and its own distributed microdefects in the asphalt mixture. Based on the Kachanov damage theory [31], the damage evolution can be defined as:
d D d t = C σ q 1 D v ,
where C, u, and v represent the temperature-dependent material parameters, and 0 D 1 .
Taking the integral of Equation (9) with the initial condition D = 0 at t = 0+ and the critical condition D = 1 at t = tf yields the damage evolution as [19]:
D = 1 1 1 + v 0 t C σ u d t 1 / 1 + v   0 t t f ,
where tf denotes the damage-induced failure time.
For σ t = σ 0 H ( t ) , the damage evolution of creep can be written as:
D = 1 1 t / t f 1 / ( v + 1 ) ,
where tf can be expressed as:
t f = 1 C ( 1 + v ) σ 0 u ,
According to the principle of the strain equivalence of the continuum damage mechanics [32], a constitutive relationship of the damaged material is obtained by replacing the nominal stress σ with the effective stress σ ¯ . The relationship between σ and σ ¯ is expressed as:
σ ¯ = σ 1 D ,
Combining Equations (5), (11), and (13), the creep strain, considering the damage evolution based on the fractional derivative operator, can be expressed as:
ε t = σ 0 K 1 D t r = σ 0 t r K 1 t / t f 1 / ( v + 1 ) = σ 0 K ' t r 0 r 1 ,
where K ' = K 1 D . As can be seen from Equation (14), the damage causes the deformation coefficient to no longer be constant, and as the damage accumulates, the deformation coefficient deteriorates and the strain increases until the specimen is damaged or fails.
In this study, the creep strain, when considering the damage, is decomposed into a strain caused by the damage and the strain from the creep hardening mechanism in the undamaged state. The strain induced by the damage (referred to as the damage strain in this study) can be expressed as:
ε d t = σ 0 t r K 1 D σ 0 t r K = 1 1 D 1   σ 0 t r K 0 r 1 ,
where εd denotes the damage strain, and this can be applied to characterize the damage softening mechanism of the asphalt mixtures.

2.3. Modeling of Fractional Derivative Creep Damage Model

A simple creep damage model is proposed by combining a fractional derivative model that describes the creep hardening mechanism with a fractional derivative damage strain that describes the damage deterioration mechanism, as shown in Figure 3. Luo et al. [33] described the asymmetry of the dynamic viscoelastic properties of asphalt mixtures using a modified fractional Zener model, where the two spring-pot elements have different fractional orders r and the same material coefficient ξ. In this study, it is suggested that not only should the two spring-pot elements describing the two mechanisms have different r, the coefficients ξ should also not be the same to better study the creep performance of asphalt mixtures. The fractional derivative creep damage model can then be expressed as:
ε t = σ 0 t r K 1 + σ 0 t α K 2 1 t / t f 1 / ( v + 1 ) 1 ,
where r and α denote the fractional orders of the two spring-pot elements, and ξ1 and ξ2 denote the coefficients of the two spring pots.

3. Uniaxial Compression Static Creep Test

3.1. Material

An AH-70 virgin binder was employed to form the AC-13 asphalt mixture, which contained basalt aggregates. The gradation of the asphalt mixture was designed according to the Standard Test Methods of Bitumen and Bituminous Mixtures for Highway Engineering (JTG E20-2011), as shown in Figure 4. The performance tests of the asphalt binder and basalt aggregate were carried out according to JTG E20-2011 and JTG E42-2005 (Test Methods of Aggregate for Highway Engineering), respectively. The performances of the asphalt binder and basalt aggregates are presented in Table 1 and Table 2. The optimum asphalt–aggregate mass ratio was determined to be 4.9% according to JTG E20-2011; the air voids in this study were set to 4.0 ± 0.2%. The Saint Venant principle clearly indicated that the larger the ratio of the height to the diameter of the specimen subjected to uniaxial stress, the greater the area beyond the end of the specimen that approached the uniaxial state. Therefore, this study adopted a gyratory compaction method to prepare a specimen with a height to diameter ratio of 1.5 for attenuating the effect of end friction on the stress state of the specimen. First, a certain amount of the mixture was compacted into a cylindrical specimen (diameter = 150 mm, height = 170 mm) using a gyratory compactor. The amount of asphalt mixture required to form the cylindrical specimen was calculated based on the maximum theoretical density, void ratio and specimen volume, and the amount should be such that the height of the specimen formed reached plus or minus 3 mm of the desired height. In this study, the amount of the mixture required to form a specimen was 7694 g. The parameters of the gyratory compactor were set according to the Superpave design method, where the internal rotation angle is 1.16° ± 0.02°, the vertical pressure was 600 ± 18 kPa, and the rotation rate was 30 r/min ± 0.5 r/min. Meanwhile, this study chose to set the required number of rotational compaction (100) as the end condition of rotational compaction. After compaction, the mold was demolded to obtain the specimen. Then, the formed specimen was core drilled, and the ends of the specimen were removed to obtain a cylindrical specimen (diameter = 100 mm, height = 150 mm).

3.2. Methodology

Creep tests were conducted on AC-13 asphalt mixtures at various stress levels and temperatures to comprehensively analyze the viscoelastic damage behavior of asphalt mixes. Various researchers selected different test temperatures and test stress levels. Ghorbani et al. [7,9] assessed the temperature effects of mixtures at 5, 20, 35, and 50 °C. Cheng et al. [24] analyzed the viscoelastic properties of different mixtures from 10 to 50 °C at 10 °C increments. Luo et al. [15] studied the viscoelastic properties of asphalt mixtures at a room temperature of 25 °C only. This study was based on the work of Cheng et al. [24] who analyzed the creep properties of the AC-13 asphalt mixture from 10 to 50 °C at 20 °C increments. Considering the temperature, the loading capacity of the instrument, and the loading time, the stress levels selected at different temperatures were different for exhibiting both the consolidation effect and the tertiary stage of creep performance at different temperatures. For example, at 50 °C and 0.8 MPa, the asphalt mixture exhibited the three stages of creep, while, at 10 °C, the asphalt mixture only exhibited a consolidation effect at 0.8 MPa. To allow the asphalt mixture to exhibit the three stages of creep, it required a larger stress level. Therefore, in this study, in order to study the change of viscoelasticity under different stresses, and to consider the comparison of viscoelasticity under the same stress at different temperatures, the creep test conditions shown in Table 3 were selected.
The specimen was kept under the testing temperature for more than 4 h before the experiment to ensure a consistent temperature throughout the specimen and to eliminate the effects of temperature inhomogeneity. Before the creep test, the specimen was preloaded with 0.05 MPa and held for 60 s to eliminate the effect of mechanical errors. Double-greased membranes was used to reduce the effect of friction between the clamps and the specimen. All tests were carried out with the UTM-25 tester (IPC Australia), which could apply a maximum stress of 2.8 MPa to the specimens. All preset stress levels were maintained for 3600 s at 30 and 50 °C. However, at 10 °C, the stress duration was increased to 10,000 s to ensure that the specimen exhibited the tertiary stage of creep at 2.8 and 2.5 MPa.

4. Test Results and Analysis

4.1. Test Results

Three tests were performed at each stress level, and the results were averaged. The averaged creep strain vs. time curves for the applied stress levels at different temperatures are shown in Figure 5.
Figure 5 shows that the creep deformation of the asphalt mixtures is stress- and temperature-dependent. Here, the creep curve at 30 °C and 0.8 MPa is used as an example illustrate the three stages of creep and the time at which each stage begins. The curve of three stages of creep and the curve of creep rate are shown in Figure 6a,b, respectively. The creep rate curve was obtained by deriving the creep curve versus time. In Figure 6, t1 and t2 represent the start times of the steady stage and tertiary stage, respectively. Additionally, during this time, the creep rate remains almost constant and the value reaches a minimum. Therefore, the change in creep rate can be used to determine which creep stage the asphalt mixture is in. Combining Figure 5 and Figure 6, at all temperatures, the creep process under larger stress levels exhibited all three stages of creep, whereas under smaller stress levels it only exhibited decelerated and stationary stages (i.e., consolidation effects). During the primary stage, the strain rate of the asphalt mixture gradually decreased with increasing in loading time; this is referred to as creep hardening. Creep hardening is different from strain hardening, which refers to the phenomenon wherein the material enters plastic deformation and the strength increases with an increase in deformation. The prerequisite for strain hardening is that material stresses should reach the yield strength and produce irreversible plastic deformation, i.e., the hardening is attributed to the strengthening of the yield strength. However, even under lower stress, the asphalt mixture exhibited varying degrees of decay creep; thus, creep hardening in the decelerating stage was not attributed to strain hardening, but to the compaction effect of the asphalt mixture under compressive stress. In this stage, the asphalt mixture was gradually compacted under the compressive stress, the voids were reduced and the strength gradually increased. In the accelerated stage, a rapid increase in deformation in addition to an increase in the strain rate with time was observed because of the presence and evolution of the damage, and this eventually led to a creep failure in the asphalt mixture. The strain rate at the accelerated creep stage was related to the stress and temperature: the higher the stress, the greater the strain rate; the higher the temperature, the greater the strain rate.
Thus, based on the testing results of the creep performance, this study considers that there are two mechanisms in the entire creep process based on the creep testing results: creep hardening and creep damage deterioration. When the creep hardening mechanism is dominant, the strain rate of the asphalt mixture gradually decreases with an increase in loading time, and the creep decays. When the damage deterioration mechanism is dominant, the strain rate increases with an increase in loading time, and the creep exhibits an accelerated state. When the two mechanisms occur in proximity, the strain rate remains constant and asphalt mixture is in a stable state. The point at which the steady stage enters the accelerated stage is called the flow time. The flow time depends on the level of stress applied; the higher the stress level, the earlier the damage deterioration mechanism occurs, and thus the flow time is shorter.

4.2. Model Verification and Parameter Determination

We verify whether the proposed creep damage model can accurately describe and reproduce all creep behavior presented in the tests. Equation (16) is determined by a curve fitting to the creep test data using the Levenberg–Marquardt optimization algorithm in 1stop. The parameters of the proposed model are listed in Table 4; the testing data of creep tests under various conditions and the simulation results of the proposed model are illustrated in Figure 7. The parameters of the correlation coefficients R2 are greater than 0.998, which indicates that the fractional derivative creep damage model can accurately characterize the three-stage creep process of asphalt mixtures under various conditions.
Table 4 indicates that, at lower temperatures (10 °C), when the asphalt mixtures are subjected to smaller stress levels (0.8–2.0 MPa), the deformation factor K1, which describes the creep hardening mechanism, remains approximately constant and its value is much greater than at the other temperatures. This indicates that the deformation factor K1 is independent of the stress levels at lower stress levels (0.8–2.0 MPa) and the asphalt mixture has greater strength at lower temperatures (10 °C). This is most likely due to the fact that asphalt mortars have a greater strength and exhibit elastic behavior at lower temperatures. Additionally, in the smaller stress range, the asphalt mixture has a smaller strain, and the compaction within the asphalt mixture is not sufficient to change the skeletal structural strength of the mixture, so the deformation factor K1 remains constant. As only consolidation effects are exhibited within the asphalt mix at these small stresses, these smaller stress ranges (0.8–2.0 MPa at 10 °C and <0.2 MPa at 30 °C) are referred to as creep consolidation stress (CCT) ranges. The creep consolidation stress range is temperature-dependent. At higher temperatures (30 °C) the asphalt mixtures exhibit a creep consolidation only at a smaller stress level (0.2 MPa); while, at lower temperatures (10 °C), the asphalt mixtures exhibit creep consolidation behavior over a wide stress range (0.8–2.0 MPa).
In the creep consolidation state, the fractional order r increases as the stress increases and the parameters (α, K2), describing the damage mechanism, are close to zero. It follows from Equation (16) that when the deformation factor K2 tends to be zero, the damage must be sufficiently small to enable the proposed constitutive model to describe the creep properties. The strain induced by the damage is too small relative to the strain of the hardening mechanism. Therefore, the Abel model can be applied to characterize the entire creep process of an asphalt mixture under smaller stress levels. This once again shows that the Abel model can describe the consolidation effect of the asphalt mixes. Under other creep conditions, both the fractional order r and deformation factor K1 increase as the stress increases for the hardening mechanism. The fractional order increases because the flow rate of the asphalt mortar in the asphalt mixture increases as the stress increases, which causes the asphalt mixture to gradually exhibit a viscous behavior. The increase in the deformation factor indicates that the resistance to deformation increases with the increasing stress. As the stress increases in the hardening mechanism, the faster the asphalt mixture is compacted, the faster the increase in the skeletal strength between the aggregates, and the smaller the deformation. For the softening mechanism, the fractional order α increases with the increasing stress, whereas the deformation factor K2 decreases with the increasing stress. The fractional order increases because the flow rate of the asphalt mortar increases as the stress increases. The decrease in the deformation factor is attributed to the presence of damage in the softening mechanism; the greater the damage to the inter-aggregate asphalt mortar as the stress increases, and the greater the deformation. For the same stress level, asphalt mixtures exhibit different viscoelastic properties at different temperatures. For example, r increases and K1 decreases as the temperature increases for the stress of 0.8 MPa. This is because, at higher temperatures, the asphalt mortar exhibits a stronger viscosity, while at lower temperatures the asphalt exhibits a stronger elasticity. When the model is applied to describe the three stages of creep, i.e., when the hardening and deterioration mechanisms are present, the fractional order of the damage softening mechanism is larger compared to the hardening mechanisms at higher temperatures (30 °C and 50 °C). This is due to the presence of damage, which makes the asphalt mixture exhibit stronger viscosity. The deformation factor of the damage softening mechanism is also larger compared to the hardening mechanism. This is because the asphalt mixture is compacted in the primary stage, which makes its resistance to deformation increase when it enters the softening mechanism. However, at 10 °C, the parameter α describing the damage mechanism is small, indicating that the damage of the asphalt mixture at this temperature tends to occur more closely to the elastic damage.

4.3. Analysis of the Relationship between Model Parameters and Stress Levels

As can be seen from Table 4, several parameters in the fractional derivative creep damage model are related to both loading stress levels and temperatures. Figure 8, Figure 9 and Figure 10 illustrate the distribution of the parameters of various stress levels at different temperatures. Table 4 indicates that the material parameter v is essentially equal for different stresses at the same temperature; therefore, parameter v at a given temperature is considered the average of the parameter v at different stresses, as indicated in Figure 8f, Figure 9f and Figure 10f. When the asphalt mixtures were subjected to linear viscoelastic stress levels, the parameters K1 remained approximately constant, and parameters α and K2 were close to zero. Thus, the parameters K1, α, and K2 were considered constants in the range of the linear viscoelastic stresses. Under other stress levels, fractional orders r and α showed a good first-order exponential decay distribution regarding the stress levels, and the deformation factors, K1 and K2, showed a good linear distribution with the stress levels, as indicated in Figure 8, Figure 9 and Figure 10. The linear viscoelastic stress range could also be obtained using K1, as shown in Figure 9b and Figure 10b. The intersection point of the two straight lines, which described the relationship between K1 and the applied stress, was the point of demarcation between the linear viscoelastic range and the nonlinear viscoelastic range. When the stress level is less than the intersection point, the asphalt mixture is in the linear viscoelastic range and vice versa in the nonlinear viscoelastic range. The creep failure times at different stresses were fitted using Equation (12) to obtain the material parameters that describe the damage at different temperatures; the results are shown in Figure 8e, Figure 9e and Figure 10e.
Based on the relationship between the model parameters and stresses given above, the creep curves could be predicted for any stress level at the three temperatures. Generally, the standard load was taken as 0.7 MPa when designing the pavement structure. Therefore, the creep behavior of the asphalt mixture was predicted and compared with the experimental data at 30 °C and 0.7 MPa, as shown in Figure 11. It was revealed that the predicted values were in good agreement with the experimental results.

4.4. Analysis of Damage Evolution

It is believed that the fractional derivative creep damage model provides a good description of the three-stage creep process in asphalt mixtures. Once the parameters of the proposed constitutive model are identified, the proposed damage evolution model can be applied to calculate the damage value D. The damage variation curves with a loading time at various stresses and temperatures are shown in Figure 12.
As indicated in Figure 12, at a given temperature, under higher stresses, the damage value D increases in an approximately linear fashion during the steady stage and rapidly during the tertiary stage. In contrast, at lower stresses, the damage value D remains at virtually zero throughout the creep process. Regardless of the applied stresses, the asphalt mixture exhibited negligible damage during the primary stage [16,34]. This is because, in the initial stages of loading, the asphalt mixture is compacted, which makes the mixture structurally hardened, at which point no damage occurs to the asphalt mixture. On a macroscopic scale, asphalt mixtures show a decrease in the void fraction and an increase in the resistance to deformation; this indicates creep hardening. Over time, the asphalt mixture structure continues to stabilize, and the flow of the asphalt mortar tends to slow down. However, the asphalt mortar and coarse aggregate interface are gradually damaged and develop into microcracks. Macroscopically, the asphalt mixture exhibits a stable void ratio and constant strain rate. As the asphalt mortar flow deformation continues to accumulate, microcracks coalesce to form macro cracks, the mineral skeleton of the asphalt mixture begins to destabilize, and the asphalt mixture enters an accelerated creep phase. Macroscopically, it is characterized by cracks in the asphalt mixture and a rapid reduction in load-bearing capacity until the mixture breaks down in the creep phase; in general, damage to the asphalt mixture is initiated in the steady stage and developed in the tertiary stage, which is consistent with the works reported by Zeng et al. [16] and Al-rub et al. [34]. The experimental data and model damage analysis show the feasability of the damage evolution proposed in this study.
Equation (11) considers only the time dependence of the damage evolution, assuming that the damage only develops over time; it does not consider the relationship between the damage evolution and the relevant mechanical quantities such as stress or strain. Equations (11) and (16) show that each moment t corresponds to a damage value and a creep strain value. Therefore, through creep tests, the relationship curve between the damage value and strain can be established for different stress levels. At different stress levels, the damage values corresponding to the strain values at each moment of the decelerated stage are small and essentially close to zero. The strain caused by the damage (damage strain, εd) is obtained by subtracting the strain in the decelerated stage from the total creep strain in the tests. According to Equation (16), εd can be expressed as:
ε d t = ε t σ 0 t r K 1 .
Based on the time dependence of the damage evolution and damage values obtained from the creep experiments, the relationship curves between the damage value D and εd for different stress levels at different temperatures is established in this study, as shown in Figure 13.
Figure 13 shows that the relationship curves between the creep D and the εd of the asphalt mixture at different stress levels overlap. Therefore, a unified damage evolution model between D and εd can be developed, where the damage parameters in the damage evolution relationship are only material- and temperature-dependent, independent of the magnitude of the stress levels. Katsuki et al. [4] used the Weibull distribution function to describe the variation in the damage with the strain, and verified the feasibility of the model using the experimental computed tomography images from Tashman at al. [35]. Therefore, in this study, the Weibull distribution function is used to describe the evolution of the damage with εd, and the model expression is shown in Equation (18). By adopting the proposed damage model (Equation (18)) to fit the experimental data under different levels of stress for different temperatures in Figure 13, the damage model curve shown in Figure 13 is obtained. The parameters are listed in Table 5. The parameters of the correlation coefficient R2 are greater than 0.992, which shows that the Weibull distribution function can be used to describe the damage evolution and demonstrate the feasibility of establishing a unified damage evolution model for different stress levels. Statistically speaking, the microelement strength of randomly distributed microdefects within the material obeys the Weibull distribution under εd:
D = 1 exp ε d p q ,
where p and q are temperature-dependent material parameters, respectively.
The nonlinear damage constitutive model proposed in this study has the advantages of fewer parameters and a clear physical meaning of parameters relative to the creep model proposed by Zeng et al. [16]. Relative to the model proposed by Zhang et al. [19], there is no need to assume that there is critical stress. Additionally, relative to other damage models [16,17,19], this study decouples the damage strain from the creep strain and establishes a unified damage evolution model for different stress levels, whose parameters are only related to material properties and temperature, and are independent of the magnitude of the applied stress.

5. Conclusions

A fractional viscoelastic creep damage model was proposed to characterize the creep process of asphalt mixtures at different stress levels and temperatures based on the fact that the asphalt mixtures exhibited different deformation mechanisms throughout the entire creep process at different stress levels. The following conclusions were drawn:
(1)
The compressive creep tests showed that the creep processes of asphalt mixtures at three different temperatures exhibited a consolidation effect at low stress levels and a tertiary stage at high stress levels. Based on the compressive creep tests, this study considered that the creep hardening and creep damage deterioration mechanisms were at work throughout the creep process. A nonlinear fractional creep damage model was proposed by combining the two models. The results of the compression creep tests were in good agreement with the proposed creep damage model, which confirmed the effectiveness of the model in describing the three-stage creep process of the asphalt mixture at different stress levels and temperatures.
(2)
The parameters of the proposed model were explicit in a physical sense, and relationships between the parameters and the applied stresses were established. The creep consolidation stress range could be obtained through the relationship between K1 and the applied stresses.
(3)
The damage value was negligible in the primary stage; it increased in an approximately linear fashion in the steady stage and rapidly during the tertiary stage. The damage was initiated in the steady stage and developed in the tertiary stage. The damage value remained at virtually zero throughout the creep process at low stresses; however, it increased rapidly with the increasing stress.
(4)
Based on the statistical quantification of the asphalt mixture damage evolution, a unified creep damage evolution model, which represented the relationship between the damage evolution and damage strain, was established. The damage parameters in this damage evolution relationship were related to the material properties and temperature and were independent of the magnitude of the stress levels. The model described well the asphalt mixture damage evolution process, and the damage obeyed the Weibull distribution.
It should be noted that only the relationship between the model parameters and applied stresses was established in the work. In order to generalize the proposed nonlinear damage constitutive model, it is also necessary to provide a way to construct a relationship between the model parameters, temperature, and the applied stress. The establishment of the relationship between the three will be considered in future studies.

Author Contributions

The authors confirm the contributions to the paper as follows: Conceptualization: X.G., Q.D.; data curation: Z.Y.; analysis and interpretation of results: Q.Z.; funding acquisition: X.G.; draft manuscript preparation: Q.Z.; review and editing: J.L. All authors reviewed the results and approved the final version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 51878162), the National Key Research and Development Program of China (Grant No. 2017YFF0205600).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ahmed, I.; Thom, N.; Bilal Ahmed Zaidi, S.; Carvajal-Munoz, J.S.; Rahman, T.; Dawson, A. Application of a novel linear-viscous approach to predict permanent deformation in simulative inverted pavements. Constr. Build. Mater. 2021, 267, 120681. [Google Scholar] [CrossRef]
  2. Doll, B.; Ozer, H.; Rivera-Perez, J.J.; Al-Qadi, I.L.; Lambros, J. Investigation of viscoelastic fracture fields in asphalt mixtures using digital image correlation. Int J. Fract. 2017, 205, 37–56. [Google Scholar] [CrossRef]
  3. Dong, N.; Wang, D.; Zhang, S.; Chen, Z.; Liang, H.; Ni, F.; Yu, J.; Yu, H. Exploring creep and recovery behavior of hot mix asphalt field cores with multi-sequenced repeated load test. Constr. Build. Mater. 2021, 279, 122435. [Google Scholar] [CrossRef]
  4. Katsuki, D.; Gutierrez, M. Viscoelastic damage model for asphalt concrete. Acta Geotech. 2011, 6, 231–241. [Google Scholar] [CrossRef]
  5. Nguyen, Q.T.; Di Benedetto, H.; Nguyen, Q.P.; Hoang, T.T.N.; Bui, V.P. Effect of time–temperature, strain level and cyclic loading on the complex Poisson’s ratio of asphalt mixtures. Constr. Build. Mater. 2021, 294, 123564. [Google Scholar] [CrossRef]
  6. Zbiciak, A.; Brzeziński, K.; Michalczyk, R. Constitutive models of pavement asphaltic layers based on mixture compositions. J. Civ. Eng. Manag. 2017, 23, 378–383. [Google Scholar] [CrossRef]
  7. Ghorbani, B.; Arulrajah, A.; Narsilio, G.; Horpibulsuk, S. Experimental and ANN analysis of temperature effects on the permanent deformation properties of demolition wastes. Transp. Geotech. 2020, 24, 100365. [Google Scholar] [CrossRef]
  8. Alrashydah, E.A.I.; Abo-Qudais, S.A. Modeling of creep compliance behavior in asphalt mixes using multiple regression and artificial neural networks. Constr. Build. Mater. 2018, 159, 635–641. [Google Scholar] [CrossRef]
  9. Ghorbani, B.; Arulrajah, A.; Narsilio, G.; Horpibulsuk, S.; Win Bo, M. Thermal and mechanical properties of demolition wastes in geothermal pavements by experimental and machine learning techniques. Constr. Build. Mater. 2021, 280, 122499. [Google Scholar] [CrossRef]
  10. Bai, F.; Yang, X.; Zeng, G. Stochastic Viscoelastic–Viscoplastic Response of Asphalt Mixture under Uniaxial Compression. J. Eng. Mech. 2017, 143, 04017049. [Google Scholar] [CrossRef]
  11. Gao, D.; Wang, P.; Li, M.; Luo, W. Modelling of nonlinear viscoelastic creep behaviour of hot-mix asphalt. Constr. Build. Mater. 2015, 95, 329–336. [Google Scholar] [CrossRef]
  12. Huang, C.; Wang, F.; Gao, T.; Gao, D.; Kachanov, L.M. A New Viscoelastic Mechanics Model for the Creep Behaviour of Fibre Reinforced Asphalt Concrete. Frat. Integrita Strut. 2018, 12, 108–120. [Google Scholar] [CrossRef]
  13. Lagos-Varas, M.; Movilla-Quesada, D.; Arenas, J.P.; Raposeiras, A.C.; Castro-Fresno, D.; Calzada-Pérez, M.A.; Vega-Zamanillo, A.; Maturana, J. Study of the mechanical behavior of asphalt mixtures using fractional rheology to model their viscoelasticity. Constr. Build. Mater. 2019, 200, 124–134. [Google Scholar] [CrossRef]
  14. Lagos-Varas, M.; Raposeiras, A.C.; Movilla-Quesada, D.; Arenas, J.P.; Castro-Fresno, D.; Muñoz-Cáceres, O.; Andres-Valeri, V.C. Study of the permanent deformation of binders and asphalt mixtures using rheological models of fractional viscoelasticity. Constr. Build. Mater. 2020, 260, 120438. [Google Scholar] [CrossRef]
  15. Luo, W.; Li, B.; Zhang, Y.; Yin, B.; Dai, J. A Creep Model of Asphalt Mixture Based on Variable Order Fractional Derivative. Appl. Sci. 2020, 10, 3862. [Google Scholar] [CrossRef]
  16. Zeng, G.; Yang, X.; Bai, F.; Gao, H. Visco-elastoplastic damage constitutive model for compressed asphalt mastic. J. Cent. South. Univ. 2014, 21, 4007–4013. [Google Scholar] [CrossRef]
  17. Zhang, J.; Li, Z.; Chu, H.; Lu, J. A viscoelastic damage constitutive model for asphalt mixture under the cyclic loading. Constr. Build. Mater. 2019, 227, 116631. [Google Scholar] [CrossRef]
  18. Zhang, J.; Zhu, C.; Li, X.; Pei, J.; Chen, J. Characterizing the three-stage rutting behavior of asphalt pavement with semi-rigid base by using UMAT in ABAQUS. Constr. Build. Mater. 2017, 140, 496–507. [Google Scholar] [CrossRef]
  19. Zhang, Y.; Liu, X.; Yin, B.; Luo, W. A Nonlinear Fractional Viscoelastic-Plastic Creep Model of Asphalt Mixture. Polymers 2021, 13, 1278. [Google Scholar] [CrossRef] [PubMed]
  20. Luo, W.; Jazouli, S.; Vu-Khanh, T. Modeling of Nonlinear Viscoelastic Creep of Polycarbonate. e-Polymers 2007, 7, 1. [Google Scholar] [CrossRef] [Green Version]
  21. Darabi, M.K.; Huang, C.-W.; Bazzaz, M.; Masad, E.A.; Little, D.N. Characterization and validation of the nonlinear viscoelastic-viscoplastic with hardening-relaxation constitutive relationship for asphalt mixtures. Constr. Build. Mater. 2019, 216, 648–660. [Google Scholar] [CrossRef]
  22. Im, S.; You, T.; Ban, H.; Kim, Y.-R. Multiscale testing-analysis of asphaltic materials considering viscoelastic and viscoplastic deformation. Int. J. Pavement. Eng. 2017, 18, 783–797. [Google Scholar] [CrossRef]
  23. Chang Kuo-Neng, G.; Meegoda Jay, N. Micromechanical Simulation of Hot Mix Asphalt. J. Eng. Mech. 1997, 123, 495–503. [Google Scholar] [CrossRef]
  24. Cheng, Y.; Li, H.; Li, L.; Zhang, Y.; Wang, H.; Bai, Y. Viscoelastic Properties of Asphalt Mixtures with Different Modifiers at Different Temperatures Based on Static Creep Tests. Appl. Sci. 2019, 9, 4246. [Google Scholar] [CrossRef] [Green Version]
  25. Xu, X.; Cui, Z. Investigation of a fractional derivative creep model of clay and its numerical implementation. Comput. Geotech. 2020, 119, 103387. [Google Scholar] [CrossRef]
  26. Gu, L.; Zhang, W.; Ma, T.; Qiu, X.; Xu, J. Numerical Simulation of Viscoelastic Behavior of Asphalt Mixture Using Fractional Constitutive Model. J. Eng. Mech. 2021, 147, 04021027. [Google Scholar] [CrossRef]
  27. Xu, Y.; Shan, L.; Tian, S. Fractional Derivative Viscoelastic Response Model for Asphalt Binders. J. Mater. Civil. Eng. 2019, 31, 04019089. [Google Scholar] [CrossRef]
  28. Celauro, C.; Fecarotti, C.; Pirrotta, A.; Collop, A.C. Experimental validation of a fractional model for creep/recovery testing of asphalt mixtures. Constr. Build. Mater. 2012, 36, 458–466. [Google Scholar] [CrossRef]
  29. Sun, L.; Zhu, H.; Zhu, Y. Two-Stage Viscoelastic-Viscoplastic Damage Constitutive Model of Asphalt Mixtures. J. Mater. Civil. Eng. 2013, 25, 958–971. [Google Scholar] [CrossRef]
  30. Koeller, R.C. Applications of Fractional Calculus to the Theory of Viscoelasticity. J. Appl. Mech. 1984, 51, 299–307. [Google Scholar] [CrossRef]
  31. Kachanov, L.M. Rupture Time Under Creep Conditions. Int J. Fract. 1999, 97, 11–18. [Google Scholar] [CrossRef]
  32. Lemaitre, J. How to use damage mechanics. Nucl. Eng. Des. 1984, 80, 233–245. [Google Scholar] [CrossRef]
  33. Luo, W.; Liang, S.; Zhang, Y. Fractional Differential Constitutive Model for Dynamic Viscoelasticity of Asphalt Mixture. China J. Highw. Transp. 2020, 33, 34–43. [Google Scholar]
  34. Al-rub, R.K.A.; You, T.; Masad, E.A.; Little, N. Mesomechanical modeling of the thermo-viscoelastic, thermo-viscoplastic, and thermo-viscodamage response of asphalt concrete. Int J. Eng. Sci. 2011, 3, 14–33. [Google Scholar] [CrossRef]
  35. Tashman, L.; Masad, E.; Little, D.; Zbib, H. A microstructure-based viscoplastic model for asphalt concrete. Int J. Plast. 2005, 21, 1659–1685. [Google Scholar] [CrossRef]
Figure 1. Schematic of fractional viscoelastic element (Abel spring-pot element).
Figure 1. Schematic of fractional viscoelastic element (Abel spring-pot element).
Materials 14 05892 g001
Figure 2. Consolidation effect and decelerating stage of the asphalt mixture: Abel model vs. test.
Figure 2. Consolidation effect and decelerating stage of the asphalt mixture: Abel model vs. test.
Materials 14 05892 g002
Figure 3. Schematic of fractional creep damage model.
Figure 3. Schematic of fractional creep damage model.
Materials 14 05892 g003
Figure 4. Gradations of the AC-13 asphalt mixture.
Figure 4. Gradations of the AC-13 asphalt mixture.
Materials 14 05892 g004
Figure 5. Creep strain curves of asphalt mixtures under different stress levels at (a) 10 °C, (b) 30 °C, and (c) 50 °C.
Figure 5. Creep strain curves of asphalt mixtures under different stress levels at (a) 10 °C, (b) 30 °C, and (c) 50 °C.
Materials 14 05892 g005
Figure 6. Schematic of the three stages (a) curve of three stages of creep, (b) curve of creep rate of asphalt mixture.
Figure 6. Schematic of the three stages (a) curve of three stages of creep, (b) curve of creep rate of asphalt mixture.
Materials 14 05892 g006
Figure 7. Comparison between the tests data of the creep tests and the fractional derivative creep damage model predicted creep curves at (a) 10 °C, (b) 30 °C, and (c) 50 °C.
Figure 7. Comparison between the tests data of the creep tests and the fractional derivative creep damage model predicted creep curves at (a) 10 °C, (b) 30 °C, and (c) 50 °C.
Materials 14 05892 g007
Figure 8. Variation of the proposed model parameters with stress levels at 50 °C (a) r, (b) K1, (c) α, (d) K2, (e) tf, and (f) v.
Figure 8. Variation of the proposed model parameters with stress levels at 50 °C (a) r, (b) K1, (c) α, (d) K2, (e) tf, and (f) v.
Materials 14 05892 g008
Figure 9. Variation of the proposed model parameters with stress levels at 30 °C (a) r, (b) K1, (c) α, (d) K2, (e) tf, and (f) v.
Figure 9. Variation of the proposed model parameters with stress levels at 30 °C (a) r, (b) K1, (c) α, (d) K2, (e) tf, and (f) v.
Materials 14 05892 g009
Figure 10. Variation of the proposed model parameters with stress levels at 10 °C (a) r, (b) K1, (c) α, (d) K2, (e) tf, and (f) v.
Figure 10. Variation of the proposed model parameters with stress levels at 10 °C (a) r, (b) K1, (c) α, (d) K2, (e) tf, and (f) v.
Materials 14 05892 g010
Figure 11. Comparison between the test data of the creep and predicted values at 30 °C and 0.7 MPa.
Figure 11. Comparison between the test data of the creep and predicted values at 30 °C and 0.7 MPa.
Materials 14 05892 g011
Figure 12. Variation curves of damage values with time under different compressive stress levels at (a) 10 °C, (b) 30 °C, and (c) 50 °C.
Figure 12. Variation curves of damage values with time under different compressive stress levels at (a) 10 °C, (b) 30 °C, and (c) 50 °C.
Materials 14 05892 g012
Figure 13. The relation curve between creep damage value and damage strain of asphalt mixture under different stress levels at (a) 10 °C, (b) 30 °C, and (c) 50 °C.
Figure 13. The relation curve between creep damage value and damage strain of asphalt mixture under different stress levels at (a) 10 °C, (b) 30 °C, and (c) 50 °C.
Materials 14 05892 g013
Table 1. Performances of asphalt binder.
Table 1. Performances of asphalt binder.
TypePenetration
(25 °C, 0.1 mm)
Ductility (15 °C, cm)Softening Point (°C)Brinell Viscosity (60 °C) Elastic Recovery (%)
Virgin binder70.1604728510
Requirements60–80≥40≥45≥160-
SpecificationT0604 2011T0605 2011T0606 2011T0625 2011T0662 2000
Table 2. Performances of basalt aggregates.
Table 2. Performances of basalt aggregates.
TypeCrushing Value (%)Los Angeles Abrasion Value (%)Needle and Flaky Particle Content (%)Water Absorption Rate (%)Apparent Density (g/cm3)
Basalt aggregates13.514.29.80.42.93
Requirements≤26≤28≤15≤2.0>2.6
SpecificationT0316 2005T0317 2005T0312 2005T0307 2005T0304 2005
Table 3. Creep experiment conditions of asphalt mixtures.
Table 3. Creep experiment conditions of asphalt mixtures.
Testing Temperature/°CStress Levels /MPa
100.8, 1.4, 2.0, 2.5, 2.8
300.2, 0.5, 0.65, 0.8, 0.7 (validation)
500.2, 0.35, 0.5, 0.8
Table 4. The parameters of the proposed constitutive models.
Table 4. The parameters of the proposed constitutive models.
Testing Tmperature
/°C
Stress
/MPa
rK1/MPa·srαK2/MPa·sαtf /svR2
102.800.2711453.2110.02918.5897350.63410.0720.999
2.500.2441316.9230.02319.867610,172.40510.1030.998
2.000.2191113.4683.849 × 10−64.027 × 10−48.577 × 1089.9870.999
1.400.2161084.97001.407 × 10−73.792 × 10139.7790.999
0.800.2121116.15607.878 × 10−127.661 × 101810.0850.999
300.800.225141.9680.396310.243746.4925.1740.999
0.650.165125.0240.350352.4322407.4425.2670.999
0.500.131104.1580.326386.82610,525.5285.3610.998
0.200.08073.06402.1961.8014 × 1065.2480.999
500.800.254144.6830.517455.56755.5522.5580.999
0.500.147106.8030.441542.917278.3712.5990.999
0.350.11596.7040.419587.8241643.6182.5830.999
0.200.08079.4320.411624.14513,912.0832.5600.998
Table 5. The parameters of damage evolution model.
Table 5. The parameters of damage evolution model.
Testing Temperature/°CpqR2
100.1730.9640.992
300.0550.7220.995
500.0210.6870.995
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, Q.; Gu, X.; Yu, Z.; Liang, J.; Dong, Q. Viscoelastic Damage Characteristics of Asphalt Mixtures Using Fractional Rheology. Materials 2021, 14, 5892. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14195892

AMA Style

Zhang Q, Gu X, Yu Z, Liang J, Dong Q. Viscoelastic Damage Characteristics of Asphalt Mixtures Using Fractional Rheology. Materials. 2021; 14(19):5892. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14195892

Chicago/Turabian Style

Zhang, Qipeng, Xingyu Gu, Zilu Yu, Jia Liang, and Qiao Dong. 2021. "Viscoelastic Damage Characteristics of Asphalt Mixtures Using Fractional Rheology" Materials 14, no. 19: 5892. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14195892

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop