Next Article in Journal
Using the Smith-Watson-Topper Parameter and Its Modifications to Calculate the Fatigue Life of Metals: The State-of-the-Art
Previous Article in Journal
Spatial Dispersion in Hypercrystal Distributed Feedback Lasing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Advances in Permanent Deformation Modeling of Asphalt Concrete—A Review

by
Mequanent Mulugeta Alamnie
1,*,
Ephrem Taddesse
1 and
Inge Hoff
2
1
Department of Engineering Science, University of Agder, 4879 Grimstad, Norway
2
Department of Civil and Environmental Engineering, Norwegian Science and Technology University (NTNU), 7491 Trondheim, Norway
*
Author to whom correspondence should be addressed.
Materials 2022, 15(10), 3480; https://doi.org/10.3390/ma15103480
Submission received: 12 April 2022 / Revised: 5 May 2022 / Accepted: 9 May 2022 / Published: 12 May 2022

Abstract

:
Permanent deformation is one of the dominant asphalt concrete damages. Significant progress has been made to realistically predict the damage. In the last decade, the mechanistic approach has been the focus of research, and the fundamental theories of viscoelasticity, viscoplasticity, continuum mechanics, and micromechanics are applied to develop the material laws (constitutive equations). This paper reviews the advancement of permanent deformation models including analogical, microstructural, and continuum-based methods. Pavement analysis using the nonlinear damage approach (PANDA) is the most comprehensive and theoretically sound approach that is available in the literature. The model coupled different damages and other phenomena (such as cracking, moisture, and phenomena such as healing, aging, etc.). The anisotropic microstructure approach can be incorporated into the PANDA approach for a more realistic prediction. Moreover, the interaction of fatigue and permanent deformation is the gap that is lacking in the literature. The mechanistic approaches have the capacity to couple these damages for unified asphalt concrete damage prediction.

1. Introduction

A flexible pavement, the longest continuous structure, comprises an asphalt concrete layer supported by unbound compacted layers (aggregate bases and subgrade). Asphaltic materials have been used for roadway construction since the end of the nineteenth century [1]. Asphalt concrete (also called bituminous mixture or hot mix asphalt) is a complex heterogeneous and three-phase material (aggregate matrix, mastic, and the air void). Such a material’s performance depends on the mixture composition, proportion, mechanical properties, and environmental conditions. It is characterized as a viscoelastic, viscoplastic, and time- and temperature-dependent material. Due to external loads and environmental factors, different distresses or damages occur in the asphalt layer. Permanent deformation or rutting is one of the primary distresses, making pavements rough and unsafe for driving, causing hydroplaning, etc. The rutting distress was noticed as a primary asphalt performance criterion [2]. The asphalt concrete’s susceptibility to permanent deformation is linked to material attributes and climatic and loading factors [3,4,5]. Material-related factors include excessive asphalt content, fine aggregate, high natural sand percentage, rounded aggregate particles, the moisture content in the mix, or granular materials and soils. From the asphalt concrete constituent properties, the chemistry of asphalt binder is the component (if not only) that makes bituminous mixtures a complex rate-dependent, nonlinear material. This nonlinear (viscous) behavior of the binder makes the permanent deformation evolution of asphalt mixtures a nonlinear mechanism [6]. For this reason, the rheological characteristics of binder (shear modulus, G*, and viscosity, η) are used to classify deformability properties of mixtures. The Strategic Highway Research Program (SHRP) rheological parameter, G*/sin δ (δthe phase angle), is the widely used criterion for rutting characterization [7]. The increased value of the criteria G*/sin δ leads to a reduced tendency to permanent deformation. However, the G*/sin δ showed limitations as it cannot predict modified binders. Thus, researchers have proposed modified rheological parameters for rutting [8]. Laukkanen et al. [8] conducted multiple stress creep-recovery (MSCR) tests on unmodified and modified binders. They concluded that the non-recoverable creep compliance parameter and accumulated strain at the end of the MSCR test showed a strong correlation and predicted mixture rutting performance compared to other rheological indicators. Meena et al. [9] investigated the rutting performance of asphalt mixture through the prediction of a resilient modulus (MR) based on the G*/sin δ rheological model. The second major constituent of asphalt concrete mixtures is the aggregate. As the primary load-carrying component, the aggregate gradation, property, angularity, texture, etc., have direct influence on the permanent deformation resistance of asphalt concrete. Research has shown that fine aggregates are better for rutting resistance [6]. Decreasing the maximum aggregate size is also good for rut resistance [10], and a fine aggregate texture is highly correlated to rutting [11]. Two theories have been assumed related to fine aggregate’s role for permanent deformation; first, fillers serve to fill the voids between aggregate particles, thereby increasing the density and strength of the compacted mixture, and secondly, the fine particles of the filler become suspended in the asphaltic binder, forming a mastic. The suspended filler particles absorb binder components, hence increasing the viscosity of the binder and, consequently, the toughness of the mixes. On the contrary, Kandhal et al. [12] reported that coarse- or fine-graded Superpave mixtures do not significantly differ in rutting resistance. In addition, the aggregate type and the chemical composition also play important roles for the creep-recovery behavior of asphalt mixtures [13]. For example, siliceous aggregate mixes show a better recovery property than do calcareous aggregates. Thus, permanent deformation is a complex phenomenon where aggregate, asphalt, and asphalt–aggregate interaction (adhesion) properties control the overall performance. These properties may change over time as a result of associated damages such as aging or moisture to the asphalt–aggregate interface and fatigue cracking. Moreover, temperature-susceptible asphalt concrete and cold weather paving, which leads to low density, are factors for permanent deformation. Other climatic factors that affect rutting are temperature, precipitation, duration, type of loading, and loading extent.
The rutting of asphalt concrete is generally related to three mechanisms. The first mechanism is related to wear rutting in the wheel path, mainly due to studded tire abrasion [14]. The second mechanism of rutting is due to the viscoplastic strain accumulation (permanent deformation) in the asphalt layer. This mechanism is caused by the densification (volume change) and shear flow at a high temperature and stress level. The third form of rutting is due to a substructural failure (subsidence) of the granular subbase, subgrade layers [15,16]. Furthermore, the development of permanent deformation is a gradual and simultaneous mechanism of densification (closing of voids), shearing (slippage due to loss of adhesion between aggregates and binder), and dilation [17] as well as the initiation and growth of micro-crack damage [18]. Again, the accumulation of permanent deformation (viscoplastic strain) within the microstructure of asphalt concrete involves three phenomena: (1) viscoplastic deformation associated with the asphalt binder, (2) rotation and slippage of aggregates (evolution of the microstructure), and (3) crack initiation and propagation (microcracks and macrocracks) [19,20,21,22]. The deformation resistance of asphalt concrete is derived from the aggregate matrix and the viscous asphalt mastic. The microstructure changes due to loading (such as air void reduction) and chemical transformation (such as aging) causes the continuous modification of aggregate matrix and asphalt mastic with time. In addition, the relaxation ability of the pavement upon load removal changes as the microstructure is continuously modified [23,24,25]. Furthermore, the growth of permanent deformation is highly dependent on the stress path and strain rate [26]. Excessive deviatoric stress and the environment cause a nonlinear, plastic, and viscous flow phenomenon where the stress–strain relationship shows strong nonlinearity [27] (especially at high strains where the linear thermo-rheological properties are not valid).
In the literature, both permanent deformation and rutting terminologies are used interchangeably. The term “Rutting” is used to describe the pavement surface roughness due to the vertical depression along the wheel path caused by the permanent deformation or wear in the asphalt layer. Permanent deformation is the accumulation of irrecoverable strain due to densification, shear deformation, and crack growth in asphalt concrete. In this paper, the term permanent deformation is used to refer to the plastic and viscoplastic strains (irrecoverable deformation). Therefore, rutting (RD, mm) is expressed as follows.
RD = i = 1 i h i ε v p , i
where h i is the ith layer’s thickness, and ε v p , i is the viscoplastic/permanent strain in the ith layer.

2. Objective and Scope

The aim of this review is to provide state-of-the art information on the developments of permanent deformation modeling for asphalt concrete. In the review, the main permanent deformation modeling theories, methods, models, and calibration tests are discussed in detail. The focus of the paper is studying the advancement of constitutive modeling approaches (analogical, microstructural, and continuum-based) and synthesize the capacity of the approaches, merits, and limitations. The robustness of permanent deformation models to account for/couple simultaneous damages such as moisture, fatigue cracking, etc., were also explored. The literature studied in this paper was collected using keywords (strings) such as viscoelastic, viscoplastic, viscodamage, permanent deformation, rutting, continuum damage, microstructure or micromechanics, mechanistic methods, creep-recovery, etc. The organization of the paper is depicted in the flow chart in Figure 1.

3. Permanent Deformation Prediction

3.1. Analytical Models

The evolution of the irrecoverable deformation of asphaltic material due to cyclic loading is described by three distinct stages, as shown in Figure 2a. The primary zone is described by the rapid accumulation of permanent deformation at a decreasing strain rate. In the secondary zone is the constant rate of permanent deformation with strain hardening, and the tertiary stage is characterized by the increasing rate of deformation accumulate and crack formation. The flow time (FT) or the flow number (FN) is defined as the time or number of loading cycles when shear deformation under constant volume commences. Several researchers have verified that the asphalt concrete’s deformation evolution showed all the three phases [28,29,30,31,32]. This three-stage deformation property is also regarded as asphalt concrete material property. The most common and simulative test used to characterize permanent deformation is the triaxial creep-recovery test. An example in Figure 2b shows that the creep–recovery deformation is dependent on confining stress. Confinement increases the friction between aggregates and increases the resistance to deformation. Although there is no conclusive research, the in situ confining pressure of asphalt concrete is approximated to be between 100 kPa and 225 kPa.
Several analytical models have been proposed over several decades to predict the permanent deformation of asphalt concrete. Sousa et al. [33] gave a summary of selected analytical models. Some of the commonly used permanent deformation models and the corresponding calibration tests are summarized in Table 1. These analytical models are used to predict rutting from typical laboratory experimental data (axial stress–strain test data). The models in Table 1 can be classified as empirical and mechanistic-empirical, which are calibrated using simulative laboratory or field data. The models give the macroscopic responses of the measured data and hardly relate the fundamental material properties. Other models are regression equations to fit rutting data and lack the explicit physical or material property for the modeling parameters.

3.2. Calibration Tests

The most common laboratory test protocols used cylindrical specimens of dimensions (diameter by height) 100 mm by 150 mm for creep-recovery, or 150 mm by 50–70 mm for creep. The shear strain from simple shear tests is also used to model permanent deformation [16,34]. The coupling of shear and axial strain components in permanent deformation modeling has not been performed yet in previous studies [35]. The indirect tensile test is also used for permanent deformation with different specimen dimensions of 150 or 100 mm diameter by 50 to 70 mm thickness [36]. Moreover, the wheel tracking test is used to simulate permanent deformation [37,38,39].
Table 1. Some selected permanent deformation analytical models.
Table 1. Some selected permanent deformation analytical models.
Model (Equation)VariablesDescription ReferenceCalibration Test
ε p = aN b a, bAccurate for secondary stages, small stress/strain deformation Creep, creep-recovery
ε p = AN B + C ( e DN 1 )
Where, A = 115 ( σ 1 σ 3 ) E *
B = 0.182 + 0.294 ( σ VM σ VL )
A, B, C, DMost widely used analytical model for permanent deformation in all three creep stages; the first part is a power function (for low stresses). and the second part is for high stresses (tertiary stage). ( σ VM —maximum stress, σ VL – plastic failure threshold; σ 1 , σ 3 —axial and lateral stress, E * —stiffness)Francken [40]Repeated triaxial compression
ε p N = AN m , A = J ( M r σ a ) S J , S, mAnalytical Power Model based on dissipated energy rate; A is a function of resilient modulus and applied stress Khedr Safwan [41]Multiple step dynamic test
ε p = ε o e ( ρ N ) β ε o , β, ρ Widely used analytical model to fit all creep stagesTseng and Lytton [42]Triaxial creep-recovery test
ε p = δ 1 ( 1 e δ 2 N ) + δ 3 ( e δ 4 N 1 ) δ 1 ,   δ 2 ,
δ 3 , δ 4
Analytical model for three creep stages, mainly developed for unbound materials ( δ 1 , δ 3 —scale primary and tertiary strain; δ2, δ 4 —rate parameters)Wilshire and Evans [43]Creep tests
ε p ε r = 10 k 1 T k 2 N k 3 k 1 ,   k 2 , k 3 Mechanistic-empirical (MEPDG) Model
ε r —resilient strain, T—temperature
Triaxial creep-recovery test
ε p I = aN I b ;
ε p II = ε p I + c ( N II N I ) ;
ε p III = ε p II + d ( e k ( N N II ) 1 )
a, b, c, d, kThree-stage rutting model (Modified Francken model) for accurate flow number identification Zhou et al. [28]Triaxial creep-recovery test
ε p = A + BN Ce DN A, B, C, DTwo phase, Linear exponent model (mainly for unbound granular materials)Cerni et al. [44]Triaxial creep-recovery test
ε p = A + BN ( C + N ) α A, B, C, αIncremental model: mechanistic based as a function of viscoplastic hardening (H and α), loading time, deviatoric stress, and rest period (A, C contain the parameters related to initial behavior of permanent deformation)Choi et al. [45]Triaxial creep-recovery
ε p = A + BN r e d ( C + N r e d ) α
N r e d = N × 10 α t o t
α t o t = α ξ p + α σ d
α ξ p = a 1 ξ p a 2 + a 3
α σ d = b 1 ( σ d P a ) b 2 + b 3
a 1 , a 2 ,   a 3 , b 1 , b 2 , b 3 , A, B, C, α A mechanistic shift model based on the load time and stress–shift function (master curve). P a is atmospheric pressure, ξ p is reduced loading time, and σ d is deviatoric stressChoi et al. [46]Triaxial creep-recovery
ε p = a ( N b + e c N ) a, b, cA three-stage model modified from Francken modelFang et al. [4]Wheel tracking, Uniaxial cyclic compression
As shown in Figure 3, different models have different accuracy for the same permanent deformation data. It is evident that the fitting accuracy is variable especially in the primary and tertiary stages of deformation. The incremental or Choi, Tseng–Lytton, and Francken models showed close predictions of the measured data. The Francken model is the most widely used for permanent deformation modeling [30].
Each model has a different accuracy of fitting all the creep stages of permanent deformation. Some of the models presented above have clear limitations such as the implicit empiricism, being unable to model load history and hardening–relaxation behavior, a lack of capacity to couple other simultaneous damages, etc. As discussed in the next sections, mechanistic permanent deformation prediction methods are aimed to resolve the limitations of empirical/mechanistic-empirical models. The latest mechanistic models have relied on rigorous material models, which are based on fundamental theories of mechanics, stress–strain relationships, and environmental factors.

4. Permanent Deformation Damage Modeling

4.1. Stress–Strain Response

The stress-, time-, and temperature-dependent viscoelastic, viscoplastic, and viscodamage properties of asphalt concrete material offered considerable challenges to accurately model the response under variable loading conditions. The response of asphalt concrete is stress path-dependent [26]. It is also a rate- and history-dependent material [47,48]. The linearity limits of asphalt concrete are 150 and 100 micro-strain in compression and tension, respectively [49], but others suggest 122 micro-strain as a limit for linear response [50]. Permanent deformation damage is induced on the asphalt beyond this strain limit at high temperatures. Moreover, the stress–strain evolution is highly dependent on the stress and strain levels, number of loading cycles (loading time), and temperature range. In a typical creep-recovery test, the stress–strain hysteresis loops evolve nonlinearly as shown in Figure 4a. The loop has the recovery and non-recovery (permanent deformation) parts. In each creep–recovery cycle, the deviatoric stress causes a non-recoverable strain and creates an open stress–strain hysteresis loop. Figure 4b shows the stress–strain responses of a constant rate compressive strength test (without recovery time). Asphalt concrete undergoes creep deformation during the load phase and a delayed recovery upon the load removal during the rest period. Traditionally, the additive decomposition of strain is applied to separate the permanent strain and recoverable strain from the total strain.
The schematic Figure 5 shows the strain components of a single pulse creep-recovery loading. The strain components can be separated into four strain components [51,52] (elastic ε e , viscoelastic ε ve , plastic ε p , and viscoplastic ε vp ). The elastic (time-independent) and viscoelastic (time-dependent) are recoverable, while plastic (time-independent) and viscoplastic (time-dependent) strains are non-recoverable parts of total strain. The total strain is expressed as follows.
ε tot = ε e + ε ve + ε p + ε vp
The permanent strain ( ε vp = ε p + ε vp ) is obtained by subtracting the viscoelastic strain ( ε ve = ε e + ε ve ) components from the total strain.
However, the strain decomposition approach is questioned when the hardening–relaxation mechanism of asphalt concrete is considered. In a loading–unloading cycle, the recoverable strain is dependent on the rest period. For tests with short rest periods, the computed viscoplastic strain can be overestimated [53]. The limitation of strain decomposition is from the inherent interaction between the viscoelastic and viscoplastic strain (the viscoelastic response is also a function of viscoplastic deformation history).

4.2. Constitutive Models

The computational modeling of asphalt concrete poses difficulties mainly due to the material nonlinearity, complexity to characterize under repeated and moving loads, and variable environmental conditions (temperature, moisture, etc.) [35,54]. The constitutive equations for the linear viscoelastic strain ( ε ve ) in undamaged conditions is defined by the Boltzmann superposition principle.
ε ve ( t ) = 0 t D ( t τ ) d σ ( τ ) d τ d τ
σ ( t ) = 0 t E ( t τ ) d ε ve d τ d τ  
where t is time; E ( t ) and D ( t ) are relaxation and creep compliance moduli, respectively; and τ is the integration variable. Prony series forms of the creep compliance and relaxation moduli are expressed as follows.
D ( t ) = D o + i = 1 N D i [ 1 exp ( t τ i ) ]
E ( t ) = E + i = 1 M E i [ exp ( t ρ i ) ]
where N and M are the total numbers of Prony terms; D o , D i , and τ i are creep compliance model coefficients; and E , E i , and ρ i are relaxation model coefficients. Often, D ( t ) is obtained from E ( t ) data via the interconversion technique [55,56] ( 0 t E ( t τ ) dD ( τ ) d τ d τ = 1 ). Once the creep compliance function is defined, the viscoelastic strain can be determined, and the viscoplastic strain is calculated from the total strain by the additive decomposition technique. A sufficient rest period is necessary to completely remove the delayed recovery strain from viscoplastic strain evolution. Cao and Kim [53] showed that 99% of viscoelastic strain is recovered for the 0.4 s pulse period and 100 s rest duration in the first cycle and 98% in cycle number 10 within 60 s. Therefore, to obtain a true viscoplastic strain, about 100 s rest period is required [45]. Experimental observations showed that the total deformation in each cycle decreases due to hardening as the number of load cycles increases. Since the viscoelastic strain is obtained from a separate dynamic modulus test (i.e., a constant viscoelastic strain), subtracting a constant cyclic viscoelastic strain from a decreasing total strain can result in a negative viscoplastic strain. That means a decreasing viscoelastic deformation model should be proposed. It is the microstructural change due to viscoplastic deformation that causes a change in viscoelastic deformation. This interaction between viscoelastic and viscoplastic deformation is referred to as viscoelastic–viscoplastic coupling, according to [53]. There is no available literature that couples the two deformations.

Analogical Models

The family of different analogical models has been used to model the viscoelastic–viscoplastic response of time-dependent materials [57,58,59,60]. The common classic mechanical models are spring, dashpot, slip device, pot, parabolic elements, etc., and the combination of these elements to analogs. The mechanical elements are advantageous to visualizing the stress and strain responses using the analogs. The Maxwell model (for the viscoelastic model), Kelvin model (for creep response), and the Burgers model are used to model the viscoelastic and viscoplastic strain (Figure 6a,b). The governing differential equations (viscoelastic constitutive equations) were developed from a number of springs and dashpots arranged in series and parallel. In the generalized Maxwell model, the same strain is shared across all elements, and the stress is additive, while in the generalized Burgers model the strains are additive, and the stress is the same for each element. It can be noted here that the generalized Burgers model shares the same framework as classical viscoplasticity models and allows nonlinearities based on stress to be accommodated more easily [57]. The viscoelastic and viscoplastic components can be calculated using the hereditary integral formulation as follows.
ε ve ( t ) = D ve ( 0 ) σ ( t ) + 0 t σ ( τ ) dD ve ( t τ ) d ( t τ ) d τ
ε vp ( t ) = D vp ( 0 ) σ ( t ) + 0 t σ ( τ ) dD vp ( t τ ) d ( t τ ) d τ
D ve and D vp are viscoelastic and viscoplastic creep compliance. The hereditary integrals in Equations (7) and (8) are different from the one in Equation (3). A formulation based on stress and the rate of compliance rather than a formulation based on the rate of stress and compliance is advantageous to avoid problems due to the sudden application of a stress in which the rate of stress can be extremely high (e.g., in a creep test). The first derivatives of the viscoelastic and viscoplastic creep compliance for a generalized Burgers model are dD ve ( t τ ) d ( t τ ) = i = 1 N 1 λ i e ( t τ ) / τ i and dD vp ( t τ ) d ( t τ ) = 1 λ , D ve ( 0 ) = D vp ( 0 ) = 0 , τ i = λ i / E i , where λ i is the viscosity of the ith Voigt element; E i is the modulus of elasticity; and λ is the viscosity of viscoplastic element. A power function for creep compliance is also used for small stress cases. Moreover, the model in Figure 7c is a modified Burgers model with a plasticity element for asphalt mixture [60,61]. The additional elastoplastic network composed of the spring and slider in parallel is used. The limit stress in the slider modeling plasticity is denoted by σ o . The authors also extended the fractional rheological model for nonlinear elastic, nonlinear viscous, and plastic properties and formulated a differential equation to characterize the viscoelastic–plastic response of asphalt concrete. Other similar analogical models such as the 2S2P1D (two springs, two parabolic. and one dashpot elements) in Figure 7a and DBN (Di Benedetto–Neifar) in Figure 7b are also frequently used to predict linear viscoelastic and creep responses (for a small number of cycles) for binders and bituminous mixtures [58,59,62]. The DBN model is a special case of the Kelvin–Voigt model where the DBN model has an elastoplastic (EP) element instead of an elastic element only.
From these models, the governing differential equation for creep compliance or stress and strain functions is derived to predict the permanent deformation response of asphalt concrete.
One can note that the parabolic elements (k, h) in 2S2P1D model, the slider element in (Figure 6c), and the spring-pot element in the fractional model (Figure 6d) have similar arrangements. The parabolic creep elements in 2S2P1D (Figure 7a), the elastic-plastic (EP) elements in DNB (Figure 7b) and the spring-pot elements in fractional model (Figure 6d) have also similar functions to model the elastoplastic response of asphalt concrete. The slip device shown in Figure 7c is placed in parallel with linear viscoelastic (LVE) element which has similar property as the fractional model in Figure 6c. The slip device functions as an irreversible deformation, which means that whenever the LVE device relaxes during unloading, the slip device locks, thereby disallowing strain recovery. Based on this phenomenology, the viscoelastic integrals are proposed for hardening and permanent deformation, considering that viscoelastic deformation, viscoplastic deformation, and hardening function are history dependent. Subramanian et al. [48] proposed a viscoelastic-like viscoplastic constitutive model for the permanent deformation of asphalt concrete. The proposed model takes the following form in Macaulay brackets.
ε vp ( t ) = 0 t D 1 ( t τ ) d σ vp ( τ ) d τ d τ
H ( t ) = H o + 0 t D 2 ( t τ ) d σ H ( τ ) d τ d τ
where H ( t ) is the material hardening variable, H o denotes the initial hardening state, σ H ( t ) and σ vp ( t ) are functions of deviatoric stress for hardening and viscoplastic deformation calculation, and D 1 ( t ) and D 2 ( t ) are compliance functions. The two stress terms are approximated by power functions, as follows.
σ H ( t ) = H 1 ( σ d ( t ) ) q 1
σ vp ( t ) = ( G 1 σ d ( t ) ) p 1 + ( G 2 σ d ( t ) ) p 2 ( H ( t ) ) α
where σ d ( t ) is the deviatoric stress history; and H 1 , q 1 , G 1 , G 2 , p 1 , p 2 , and α are parameters. This model considers only hardening during loading pulses and ignored the softening mechanism. Based on the model in [48], Cao and Kim [53] proposed a viscoplastic model using “internal stress” as the hardening variable. They hypothesized that as soon as the deformation of the LVE device is constrained (in Figure 7c), the internal stress inside the device starts to develop due to stress relaxation. As illustrated in Figure 8, the internal stress decreases with time due to LVE device relaxation before the next load cycle but has a direction opposite to the external load. Once the applied stress rises above the level of the concurrent resisting internal stress in the LVE device, the slip device is unlocked and becomes frictionless, allowing the overall deformation of the mechanical analog to increase in a viscoelastic fashion. The proposed viscoelastic-type viscoplastic model takes the following form [53].
ε vp ( t ) = 0 t D ( t τ ) d σ d ( τ ) d τ d τ
This model introduces the coupling of viscoelastic and viscoplastic responses using the internal stress as a hardening variable and viscoelastic-like hardening–relaxation spectrum.

4.3. Continuum Based Models

4.3.1. Damage Density

Continuum mechanics is a standalone and widely applied theory for damage formulation. The concept of continuum damage mechanics (CDM) was pioneered by Kachanov [63], who introduced a scalar measure called damage variable or damage density ϕ , which is defined as follows.
ϕ = 1 A A d = A d A A d
where A—real (intact) area, A d —damaged area, A d A is the area of micro-damage, ϕ = 0 means the initial state, and ϕ = 1 mean complete rupture. Based on the damage density function and effective area, the effective stress concept in CDM is defined as follows.
σ ¯ ij = σ ij 1 ϕ
where σ ¯ ij is the effective stress tensor in an undamaged configuration; and σ ij is the nominal Cauchy tensor in damage configuration. For more accurate modeling, the damage evolution is modified as follows [64].
σ ¯ ij = σ ij ( 1 ϕ ) 2
The classic Kachanov–Robotnov damage models [63,65] were extensively used for creep damage ( ϕ c ) modeling for different materials.
Kachanov :   ϕ c ˙ = G ( σ A ) r ( 1 ϕ c ) k
Robotnov :   ϕ c ˙ = C 1 σ γ ( 1 ϕ c ) η ,   ϕ c ˙ = C 2 exp ( k ε ) ( 1 ϕ c ) η
where A, G, r, k, and C 1 , γ , η , C 2 , k are material constants. Recently, Darabi et al. [64] developed a continuum viscodamage model using the effective total strain ( ε ef T ), viscoplastic hardening, and temperature coupling functions.
ϕ c ˙ = Γ o φ [ Y ¯ ( 1 ϕ c ) 2 Y 0 ] q exp ( k ε ef T ) G ( T )
where Γ o φ is reference damage viscosity, Y 0 is reference damage force, Y ¯ is damage driving force in effective configuration, G ( T ) is temperature coupling term, and k and q are constants.

4.3.2. Viscoplasticity

The classic Perzyna viscoplastic hardening rule [66] assumes a constant hardening variable for a cyclic creep-recovery load, defined as
ε ˙ vp = Γ vp f F ( σ ) σ
where Γ vp is the viscoplastic fluidity parameter such that 1 / Γ vp is a measure of viscoplastic viscosity; and F ( σ ) σ is a measure of the direction of viscoplastic strain. The classical hardening assumes that the viscoplastic strain rate decreases progressively with an increase in loading time. However, the hardening function is not constant due to the hardening–relaxation behavior [23]. Different researchers pointed out that the Perzyna-type rate models have limitations such as that the model cannot capture the load history effect, the relaxation or softening behavior during the rest period is ignored, and it assumes a constant hardening parameter [23,47,48]. As illustrated in Figure 9, the viscoplastic strain rate is no longer a decreasing function. The quantity qvp is the hardening–relaxation internal state variable that memorizes the maximum experienced viscoplastic strain for which the hardening recovery has occurred.

4.4. Mechanistic Methods

4.4.1. Pavement Analysis Using the Nonlinear Damage Approach (PANDA)

Pavement analysis using the nonlinear damage approach (PANDA) is the latest generation of mechanistic pavement design approach [67]. The PANDA is a mechanistic-based pavement analysis method that is founded on three classic theories: (1) Schapery’s [68] nonlinear viscoelasticity, (2) Perzyna’s [66] viscoplasticity, and (3) Darabi’s [64] viscodamage constitutive relationship. Based on the three constitutive equations, the PANDA approach has an unlimited capacity to couple different damage mechanisms of pavement structures. The approach coupled the temperature, rate, and time-dependent viscoelastic and viscoplastic models to predict the permanent deformation of asphalt concrete. For example, healing, aging, hardening–relaxation, moisture-induced damage, and other behaviors are conveniently incorporated into the PANDA approach [69,70,71,72,73]. In Figure 10, the development of the PANDA model and the constitutive equations are summarized.
(1)
First, the linear and nonlinear viscoelastic variables are obtained from dynamic modulus (for linear viscoelastic) and creep-recovery tests (nonlinear viscoelastic). The nonlinear viscoelastic strain is formulated using the well-known Schapery’s viscoelastic constitutive equation [68].
(2)
Secondly, the viscoelastic strain is deducted from the total strain to extract the viscoplastic strain from the same creep-recovery test data using the strain decomposition principle. Then, the classic Perzyna’s viscoplasticity [66] theory is adopted to predict the viscoplastic strain evolution. The Drucker–Prager yield surface function is often used [23,74].
ε ˙ vp = Γ vp f σ y o N F ( σ ) σ
(3)
The third foundation of the PANDA constitutive model is the viscodamage model [64] using the continuum damage mechanics (CDM) theory. The effective strain is used in effective configuration. The viscodamage model mainly predicts the permanent deformation in the tertiary creep.
Therefore, the PANDA model encompasses nonlinear viscoelastic, viscoplastic (hardening), and viscodamage responses. The thermo-piezo-rheological viscoelastic properties [75] coupled with the viscoplastic yield criteria (the Drucker–Prager yield surface) is an integral part of the PANDA method [76]. Several material parameters need to be optimized to calibrate the viscoelastic–viscoplastic–viscodamage, the hardening–relaxation, moisture damage, and healing responses of asphalt concrete. The parameters and their physical meaning are presented in Table 2 Despite the unlimited capacity of the PANDA approach, calibrating the large number of model parameters is a laborious task. Hence, a systematic procedure is followed to obtain material parameters with a smaller number of tests. Once robust mechanistic constitutive models are developed and calibrated, the numerical implementation (finite element modeling) is performed using the user-defined material (UMAT) tool to define material laws in commercial software such as ABAQUS. Finite element modeling (FEM) is conducted at the desired modeling space (2D or 3D), and realistic tire–pavement contact [77], traffic, and full pavement structure (asphalt, base, subbase, and subgrade), etc., can be constructed [67,76,78,79,80]. Although the PANDA approach is still at a research stage, it is evident that advantages as well as limitations can be listed. In Table 3, some of the merits and limitations of the PANDA are described. The calibration of PANDA models used uniaxial creep, uniaxial constant stress creep-recovery, the crosshead strain rate test, and multiple stress creep tests, etc. The influence of confining pressure on the linear viscoelastic as well as nonlinear viscoelastic responses is significant [75,81]. Most studies that tried to calibrate the PANDA models used uniaxial test data, or some used a single confining pressure.

4.4.2. Microstructural Based Models

The micromechanics approach is probably the best way to account the effects of individual mixture constituents and their interactions and the anisotropy of heterogenous asphalt mixture. The microstructure change in asphalt concrete is mainly due to the friction between the aggregates and interlocking bond breakage. This mechanism is responsible for the accumulation of permanent deformation rather than the recoverable part of the deformation. The continuous increase in the resistance of the material due to the permanent microstructure rearrangement is physically related to the strain hardening. The hardening parameter reflects the combined effect of the cohesion of asphalt binders, the adhesion properties between binder and aggregate, and the frictional properties of the aggregate structure. The fabric of granular media refers to the size, shape, and arrangement of the solid particles and the associated voids. The scalar quantity, like void ratio, is not capable of characterizing the directional nature of fabric and describing the state of packing of granular materials [85]. The microstructure approach is necessary to consider the directional nature of granular fabric. The approach is capable of modeling nonlinearities such as heterogeneity, aggregate distribution, anisotropy, crack, and air void. The micromechanics coupled with the continuum damage approach is a powerful way to model the permanent deformation of granular asphalt concrete with an appropriate representative volume element (RVE) [19,20,21,22]. The fundamental element of the granular microstructure is a directional unit vector m and the vector magnitude Δ (Figure 11).
Δ = 1 M [ k = 1 M ( cos 2 θ k ) 2 + k = 1 M ( sin 2 θ k ) 2 ] 1 2  
where M is the total number of objects analyzed in an image; and θ k is the orientation of the unit vector n and ranges between −90 and +90°. Theoretically, the value of Δ is between zero and unity; Δ = 0 indicates that objects are completely randomly distributed, which is analogous to isotropic materials; and Δ = 1 indicates that objects are entirely oriented in one direction, which is analogous to perfectly transverse anisotropic materials. Oda and Nakayama [86] introduced a symmetrical microstructure tensor F i j that gives a measure of the two-dimensional anisotropy produced by the preferred orientation of non-spherical particles.
F i j = ( ( 1 Δ ) / ( 3 + Δ ) 0 0 0 ( 1 + Δ ) / ( 3 + Δ ) 0 0 0 ( 1 Δ ) / ( 3 + Δ ) )
For simplicity the anisotropy value can be taken as constant (the initial value). For isotropic elements, F 11 = ( 4 Δ ) / 3 ( 3 + Δ ) and F 22 = F 33 = ( 2 Δ ) / 3 ( 3 + Δ ) .
The microstructure tensor is incorporated in the Drucker–Prager yield function by modifying the first and second invariants. The classical Perzyna’s viscoplastic and continuum damage models are then used to formulate the microstructure-based viscoplastic strain model to characterize permanent deformation of asphalt concrete. The detail derivation of the model can be found in the study [19,20,86]. The final form of the viscoplastic model in a triaxial compression test (axial strain) based on micro-structural anisotropy takes the following form.
ε ˙ vp = [ X β ( 1 Y ) ] 3 2 X + 3 β Y X + 3 β 2 + 3 2 β 2 Y 2 ε ˙ ij vp ε ˙ ij vp = ε ˙ 11 vp  
where X = 1 3 4 9 24 μ ( Δ 3 + Δ ) 2 , Y = 4 3 24 λ ( Δ 3 + Δ ) 2 , and λ and μ are anisotropy coefficients that reflect the effect of the aggregate anisotropic distribution on the confining and shear stresses, respectively. The viscoplastic model considers phenomena including the effect of the binder fluidity ( Γ ); confinement and aggregate friction ( α ); aggregate interlocking and dilation ( β ); binder cohesion and its adhesion to the aggregates ( κ ); anisotropy of aggregate distribution ( Δ ); and microstructure damage ( ξ ). As shown in Figure 12, the anisotropy parameter significantly contributes to compressive viscoplastic behavior and has little effect or no effect on the tensile and linear viscoelastic property [22]. In a conventional continuum method (isotropic mode), phenomenologically motivated, microstructural-based viscoplastic models have been developed, for example in [23,25] and others.
Microstructure modeling is directly integrated with the utilization of digital technologies to capture the particles arrangement, deformation, etc., in the granular materials packing. Digital image correlations, X-ray chromatography, and other tools were used for asphalt concrete. Coleri et al. [88] used the X-ray computed tomography (CT) method to model asphalt concrete rutting in a full-scale heavy vehicle simulation (HVS) test site. The object segmentation and processing procedure using a digital camera is shown in Figure 13, and an example of the discretization phases is shown in Figure 14. Using the digital technology, parts of a heterogenous mixture can be easily dissected, modeled, and analyzed into the FEM and other post processes.
In general, the microstructural based micromechanics modeling is ideal to consider the effect of individual constituent particles in viscoplastic evolution [89]. However, this approach needs image analyzing tools, discretization, and finite element modeling. Several material behaviors such as hardening–relaxation, healing, stress history, viscodamage, etc., need to be integrated for accurate formulation. Limited literature can be found in this area. With the advancement of morphological image analyzing tools, the microstructural approach will be the next active research area to characterize permanent deformation.

4.4.3. Finite Element Simulation

Even if pavement is considered a homogeneous, isotropic, elastic, multi-layered system, the calculated stress–strain distribution in the structure under the simple action of a passing wheel is very complex. For a realistic modeling of pavements, each layer and material constituent should be modeled with the appropriate material constitutive law, and the interaction of each layer should be analyzed as a pavement system. The finite element method is the most versatile approach in the mechanistic pavement design approaches. The FEM is an integral part of the PANDA approach; the modeling space (2D or 3D), wheel loading configuration, constitutive law, etc., are of great importance [67,69,76,79]. Collop et al. [57] implemented Burgers analogical constitutive model into FEM to predict viscoplastic deformation. A comprehensive FEM study by [80] investigated the effect of loading scenario in 2D and 3D modeling spaces and implemented the PANDA constitutive law. They investigated 11 different loading modes (given in Table 4) and compared three different material constitutive equations: elasto-viscoplastic, viscoelastic–viscoplastic, and viscoelastic–viscoplastic–viscodamage or PANDA. Based on wheel tracking test data, they concluded that the 2D model significantly overestimates permanent deformation (rutting) compared to the 3D moving load case (which is the most realistic case). Moreover, they found that the viscoelastic–viscoplastic–viscodamage model gives higher rutting damage predictions compared to elasto-viscoplastic and viscoelastic–viscoplastic models. The viscoelastic–viscoplastic and elasto-viscoplastic models give close predictions, where the form gives a still higher rutting prediction. On the contrary, another study using the viscoelastic–viscoplastic–viscodamage constitutive model and 3D FEM modeling of tire–pavement interaction showed that the pulsatile and equivalent loading assumptions overestimated rutting compared to the realistic moving load [90]. The FEM simulation has been applied for different viscoelastic, viscoplastic, and crack modeling. The linear viscoelastic simulation has been modeled using the analogical models (e.g., Maxwell model), the Prony series, and time–temperature superposition principles in commercial software such as ABAQUS.

5. Permanent Deformation and Fatigue Interaction Damage

Permanent deformation and fatigue cracking are the two dominant pavement damage mechanisms. Extensive research has been conducted to characterize the two damages independently. Advanced mechanistic models were proposed [64,83,84,91]. Fatigue damage occurred due to the formation and propagation of cracks due to repetitive loads. The bottom-up cracking of the asphalt concrete layer was the traditional fatigue mechanism. However, experimental and field observations showed that top-down cracking due to tire compression is also the cause of fatigue damage [92,93]. The crack initiation in compression occurs when the viscoplastic strain hardening reaches saturation at the flow number. The extra energy at maximum saturation (hardening) is consumed to initiate microcracks and increase phase angle. When asphalt with pre-existing cracks is subjected to a compressive load, wing cracks develop and propagate parallel to the load direction [94,95], contrary to tensile loading (cracks grow perpendicular to the stress direction). When pavement temperature is considered, fatigue cracking is critical at low temperatures, and permanent deformation is a high-temperature damage. At elevated temperatures, the critical energy threshold (the mixture relaxes faster) becomes more resistant to micro-cracks and needs more energy to initiate cracks [96]. On the other hand, deformation via plastic flow (aggregate re-orientation) is dominant at high temperatures. When asphalt concrete is subjected to repetitive loading, energy is dissipated by viscous flow and/or plastic flow, leading to fatigue cracking and/or permanent deformation [97], and some part of the energy is transferred into heat [50]. The energy dissipation caused material ductility exhaustion, hardening, and viscoplastic flow. The classic energy balance principle states that the decreasing rate of potential energy (stored and recoverable) during crack initiation is equal to the dissipated energy rate due to plastic/viscoplastic deformation and crack opening, and different failure criteria were proposed [98,99]. The dissipated energy (W) in a cyclic load is the area under the stress–strain hysteresis and expressed as follows.
W = 0 τ σ ( t ) ε ( τ ) τ d τ
where τ is the integration variable, σ is the stress function, and ε is the strain function. One hysteresis cycle consists of the total strain energy, which is the sum of damage-causing dissipated energy (the hysteresis area) and the recovery strain energy (stored energy).
The classic approach considered rutting and fatigue separately (cf. Figure 15a) [100]. However, both fatigue and permanent deformation damages are caused by the same load in pavements. The difference is which damage mode dominates at different temperatures. Therefore, a realistic damage modeling should consider the simultaneous damage evolution of fatigue and permanent deformation. An interaction domain of rutting and fatigue can be shown in Figure 15b [101]. For an interaction damage condition, the dissipated energy (in Equation (26)) can be decomposed further into permanent deformation and fatigue damage-causing energies. As discussed above, the existing models for permanent deformation or fatigue damages are based on the idealized hysteresis loops of pure creep and pure fatigue. Only very few attempts can be found; for example, Luo et al. [99] investigated the permanent deformation behavior of pre-cracked asphalt concrete.
In a cyclic creep-recovery test, material hardening grows, and the hysteresis loops shift horizontally due to the accumulation of irrecoverable viscoplastic strain (open loops). On the other hand, in the idealized fatigue test (either stress- or strain-controlled), cyclic load is applied that causes stiffness reduction and a phase angle increment with a negligible viscoplastic strain (hysteresis loops do not shift horizontally or loops are closed). In simultaneous fatigue-permanent deformation damage, the material hardening contributes to initiate fatigue cracking at intermediate temperatures, and the fatigue cracks accelerate the microstructure change and contribute to viscoplastic deformation. The viscous flow creates a plastic zone (the crack initiation point), and fatigue cracking can evolve without plastic flow (aggregate movement and re-orientation). That means the simultaneous occurrence of the two damages cannot be ignored as they are interdependent. Thus, the hysteresis loops can be modified by superposing the two loops (Figure 4a and Figure 16a) to account for the interaction of the two damages, as shown in Figure 16b. Similar fatigue–creep interaction loops can be found elsewhere for steel [102,103]. A detailed discussion of the interaction of the two damages is presented in next sections.
One of the limitations for interaction damage modeling is the lack of integrated testing protocol for permanent deformation and fatigue. Some attempts were made, such as Indirect Tension (IDT) testing for deformation and fracture and applying a haversine loading waveform in fatigue tests [104,105]. The haversine load is used with the assumption that the pulse can be decomposed into pure creep and pure sinusoidal components. However, the creep-recovery behavior cannot be captured with such approaches. In another study, Zhang et al. [106] conducted a destructive dynamic modulus test to simultaneously characterize permanent deformation and fracture properties using 16 different asphalt mixtures. In their work, they successfully modeled viscoelastic, viscoplastic, and viscofracture properties using a new destructive dynamic modulus test. Another reason for the independent treatment of fatigue and rutting damages on asphalt concrete damage is the perception that permanent deformation (creep) is expected in the early life of pavements, while fatigue is high-cycle damage after asphalt concrete accumulates sufficient hardening and age. However, this precedence is dependent on several factors such as temperature, loading range (loading time, rest period, deviatoric stress), material stiffness, number of cycles, and other climatic factors. For example, fatigue cracking can evolve before creep damage in stiff and thick pavements. The top-down cracking due to tire compression can develop without permanent deformation. Some field and laboratory observations have also shown that fatigue cracking accompanied rutting [92]. Therefore, both the creep–fatigue and fatigue–creep damage sequences can occur in asphalt concrete pavements [99].
The main objective in the design and service life estimation of asphalt concrete is to estimate the number of cycles ( N f ) to initiate fatigue cracking and/or the critical strain ( ε c —rupture strain) that causes the flow or rutting of asphalt concrete. The simplest way to evaluate the interaction damage of creep and fatigue is separately calculating the creep and fatigue damages and adding them together.
n N f + ε ε c = 1
By the life fraction rule, the sum of fatigue damage ( ϕ f ) and creep damage ( ϕ c ) equals a certain damage density value, ϕ . The nonlinear summative decomposition of damage can be expressed as follows for a more general representation.
d ϕ = d ϕ c + d ϕ f = f c ( ε , σ ,   T , ϕ c + ϕ f ) dN c + f f ( ε , σ ,   T , ϕ c + ϕ f ) dN f
The classic continuum damage models were used to define the damage rates for creep [63,65] and fatigue damage [107,108]. Lemaitre et al. [109] proposed sequential damage interaction, and the total damage accumulation during one cycle of creep followed by fatigue and fatigue followed by a creep sequence can be expressed as follows (respectively).
ϕ c _ f = 0 n c δ ϕ c δ N δ N   + n c n c + n f δ ϕ f δ N δ N
ϕ f _ c = 0 n f δ ϕ c δ N δ N   + n f n f + n c δ ϕ f δ N δ N
For time-dependent materials such as asphalt concrete, damage densities ϕ c _ f and ϕ f _ c will not be the same for the same number of cycles in both sequences. This is because of the different modes of damage formation in creep and fatigue and the possible interactive damage one on another. A parameter can be defined for “interactive” damage. Skelton et al. [101] presented analytical expressions to allow creep to be modified by fatigue and fatigue to be modified by creep. The combined equation for the “creep–fatigue” and “fatigue–creep” interactive damages can take the following form.
ϕ f 1 I c f ϕ c + ϕ c 1 I f c ϕ f = 1
The interaction coefficients I cf (creep on fatigue) and I fc (fatigue on creep) can take any value between zero and unity. The damage density for fatigue and creep can be determined by rearranging terms and solving a quadratic solution. In the creep–fatigue sequence, it is assumed that pure creep is followed by fatigue damage. The fatigue damage part is modified to account for the pre-existing creep damage. Therefore, the creep–fatigue and fatigue–creep interactive damages can be expressed as follows, respectively.
ϕ c + ϕ f 1 I c f ϕ c = 1 , ϕ f + ϕ c 1 I f c ϕ f = 1
Re-arranging and solving for ϕ c and ϕ f gives the following expressions.
ϕ c = ( 1 + I cf ) ( I cf 1 ) 2 + 4 I cf ϕ f 2 I cf ,   ϕ f = ( 1 + I fc ) ( I fc 1 ) 2 + 4 I fc ϕ c 2 I fc
The creep–fatigue and fatigue–creep interaction coefficients I cf and I fc are non-zero and take any value between −1 and 1. Figure 17 shows the relationship between creep and fatigue damage densities at random values of creep–fatigue interaction coefficients.

6. Summary, Conclusions, Further Research

6.1. Summary

The practice is always trailing the theory in pavement design, and the case of the permanent deformation prediction is no different. Many highway agencies and road administrations are practicing an empirical method of permanent deformation characterizations (for example, the Marshall test). Some progresses are being made to transition from empirical to mechanistic-empirical design methods. The primary motivation of this literature study was to explore the advancement of asphalt concrete permanent deformation characterization, constitutive modeling, and application. The current state of research is focused on the formulation of a mechanistic method that applies the fundamental theories of mechanics and materials to predict permanent deformation damage. In the last decade, promising advancements have been made in the development of comprehensive and coupled permanent deformation modeling. Pavement analysis using the nonlinear damage approach (also known as PANDA) constitutive model is one of the notable progresses that is gaining wide acceptance (as the next generation mechanistic pavement design approach). Another promising area of progress is the microstructural approach aided with 3D digital image analysis equipment and finite element modeling. This digital technology-based permanent deformation modeling is also the future prospect for accurate permanent deformation modeling and prediction of different asphalt mixtures. The mechanistic methods offer unlimited potential to expand the modeling parameters, and different damages and phenomena such as aging, healing, moisture damage, pre-crack, etc., can be coupled to unify damage prediction. The drawbacks of mechanistic models are sophistication, requiring extensive testing, and the calibration of several modeling variables. For example, the viscoelastic–viscoplastic–viscodamage model requires more than 21 model variables to be optimized using at least two different experiments and several test repetitions at different temperatures stresses and strains. The improvement of the latest models from the classic viscoplastic strain hardening model is the consideration of cyclic hardening and relaxation mechanisms and the viscodamage of asphalt concrete. Although the mechanistic method is theoretically appealing, the calibration cost and rigorous equations can be considered the limitations. The micromechanics approach considers the evolution of permanent deformation related to changes in the microstructure of asphalt concrete constituents (aggregates and mastic). Thus, the micromechanics method is regarded as the most realistic way of modeling heterogeneous materials such as asphalt concrete. Based on the extensive literature study, the permanent deformation prediction and modeling approaches can be categorized into four aspects, as shown in Table 5.

6.2. Conclusions

A review study was conducted to explore the state of the art on permanent deformation prediction from the 1960s’ pure empirical to latest mechanistic (from 2011) methods. The review study revealed that the latest constitutive models integrate and couple different theories, i.e., continuum mechanics (1958), nonlinear viscoelastic (1969), viscoplasticity (1971), and viscodamage (2011), along with the crucial time–temperature superposition principle. Such coupling techniques offered advantages of integrating different asphalt concrete damages and opened the possibility of a unified asphalt damage model in the future. The PANDA model is one of the most comprehensive permanent deformation modeling approaches available in the literature (the next generation to the mechanistic-empirical method). The calibration and/or validation tests are reliant upon the conventional creep or creep-recovery tests in either confined or unconfined modes. The computation and experiment cost of mechanistic methods are the limitations. The practical application of the mechanistic models is very limited at the moment. Moreover, the fatigue–permanent deformation (rutting) interaction is often ignored in the existing (studied) literature. It is inferred that both damages can evolve simultaneously as the same load caused both damages. The mechanistic approach has the potential to couple the two predominant pavement damages. From the extensive study, it can be synthesized that a unified permanent deformation damage model can be developed by integrating continuum damage and microstructure approaches and coupling fatigue, moisture, healing, aging, and other physical and chemical phenomena in asphalt concrete.

6.3. Further Research

The mechanistic method is “universally” applicable to predict damage regardless of climatic conditions, stress state, or material type. This characteristic presents wide, open research questions, for example, (1) developing a unified pavement damage performance prediction model, (2) a coupled model for fatigue and rutting damages using continuum mechanics and viscoelastic and viscoplastic theories, (3) considering (coupling) different strains such as shear and axial strains in permanent deformation prediction models, and (4) developing simplified (unified) asphalt concrete test methods to characterized different damages simultaneously (fatigue, rutting, moisture, etc.). The authors of this paper are conducting research on the simultaneous creep–fatigue damage evolution in a sequential manner.

Author Contributions

Conceptualization, M.M.A.; methodology, M.M.A.; investigation, M.M.A.; writing—original draft preparation, M.M.A.; writing—review and editing, M.M.A., E.T. and I.H.; supervision, E.T. and I.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Roberts, F.L.; Mohammad, L.N.; Wang, L.B. History of Hot Mix Asphalt Mixture Design in the United States. J. Mater. Civ. Eng. 2002, 14, 279–293. [Google Scholar] [CrossRef]
  2. Vesic, A.S.; Domaschuk, L. Theoretical Analysis of Structural Behavior of Road Test Flexible Pavements; National Academy of Science: Washington, DC, USA, 1964. [Google Scholar]
  3. Archilla, A.R.; Madanat, S. Development of a Pavement Rutting Model from Experimental Data. J. Transp. Eng. 2000, 126, 291–299. [Google Scholar] [CrossRef]
  4. Fang, H.; Liu, Q.; Mo, L.; Javilla, B.; Shu, B.; Wu, S. Characterization of three-stage rutting development of asphalt mixtures. Constr. Build. Mater. 2017, 154, 340–348. [Google Scholar] [CrossRef]
  5. Celauro, C. Influence of the Hourly Variation of Temperature on the Estimation of Fatigue Damage and Rutting in Flexible Pavement Design. Int. J. Pavement Eng. 2004, 5, 221–231. [Google Scholar] [CrossRef]
  6. Delgadillo, R.; Bahia, H.U. The Relationship between Nonlinearity of Asphalt Binders and Asphalt Mixture Permanent Deformation. Road Mater. Pavement Des. 2010, 11, 653–680. [Google Scholar] [CrossRef]
  7. Kennedy, T.W.; Huber, G.A.; Harrigan, E.T.; Cominsky, R.J.; Hughes, C.S.; Von Quintus, H.; Moulthrop, J.S. Superior Performing Asphalt Pavements (Superpave): The Product of the SHRP Asphalt Research Program; Strategic Highway Research Program, National Research Council: Washington, DC, USA, 1994. [Google Scholar]
  8. Laukkanen, O.-V.; Soenen, H.; Pellinen, T.; Heyrman, S.; Lemoine, G. Creep-recovery behavior of bituminous binders and its relation to asphalt mixture rutting. Mater. Struct. 2014, 48, 4039–4053. [Google Scholar] [CrossRef]
  9. Meena, S.; Biligiri, K.P. Binder rheological predictive models to evaluate rutting performance of asphalt mixtures. Constr. Build. Mater. 2016, 111, 556–564. [Google Scholar] [CrossRef]
  10. İskender, E. Rutting evaluation of stone mastic asphalt for basalt and basalt–limestone aggregate combinations. Compos. Part B Eng. 2013, 54, 255–264. [Google Scholar] [CrossRef]
  11. Masad, E.; Olcott, D.; White, T.; Tashman, L. Correlation of Fine Aggregate Imaging Shape Indices with Asphalt Mixture Performance. Transp. Res. Rec. J. Transp. Res. Board 2001, 1757, 148–156. [Google Scholar] [CrossRef]
  12. Kandhal, P.S.; Cooley, L. Coarse- Versus Fine-Graded Superpave Mixtures: Comparative Evaluation of Resistance to Rutting. Transp. Res. Rec. J. Transp. Res. Board 2002, 1789, 216–224. [Google Scholar] [CrossRef] [Green Version]
  13. Singh, M.; Kumar, P.; Anupam, A.K. Effect of type of aggregate on permanent deformation of bituminous concrete mixes. Road Mater. Pavement Des. 2016, 17, 417–433. [Google Scholar] [CrossRef]
  14. Lundberg, J.; Janhäll, S.; Gustafsson, M.; Erlingsson, S. Calibration of the Swedish studded tyre abrasion wear prediction model with implication for the NORTRIP road dust emission model. Int. J. Pavement Eng. 2019, 22, 432–446. [Google Scholar] [CrossRef] [Green Version]
  15. Ahmed, A.; Erlingsson, S. Evaluation of permanent deformation models for unbound granular materials using accelerated pavement tests. Road Mater. Pavement Des. 2013, 14, 178–195. [Google Scholar] [CrossRef]
  16. Deacon, J.A.; Harvey, J.T.; Guada, I.; Popescu, L.; Monismith, C.L. Analytically Based Approach to Rutting Prediction. Transp. Res. Rec. J. Transp. Res. Board 2002, 1806, 9–18. [Google Scholar] [CrossRef]
  17. Song, J.; Pellinen, T. Dilation Behavior of Hot Mix Asphalt Under Triaxial Loading. Road Mater. Pavement Des. 2007, 8, 103–125. [Google Scholar] [CrossRef]
  18. Park, S.W.; Richard Kim, Y.; Schapery, R.A. A viscoelastic continuum damage model and its application to uniaxial behavior of asphalt concrete. Mech. Mater. 1996, 24, 241–255. [Google Scholar] [CrossRef]
  19. Tashman, L.; Masad, E.; Zbib, H.; Little, D.; Kaloush, K. Microstructural Viscoplastic Continuum Model for Permanent Deformation in Asphalt Pavements. J. Eng. Mech. 2005, 131, 48–57. [Google Scholar] [CrossRef]
  20. Yu, H.; Muhunthan, B.; Shen, S. Anisotropic Nonlinear Elastoviscoplastic Model for Rutting of Asphalt Mixtures. J. Eng. Mech. 2014, 140, 242–249. [Google Scholar] [CrossRef]
  21. Kim, Y.-R.; Lutif, J.E.S.; Allen, D.H. Determining Representative Volume Elements of Asphalt Concrete Mixtures without Damage. Transp. Res. Rec. J. Transp. Res. Board 2009, 2127, 52–59. [Google Scholar] [CrossRef]
  22. Underwood, S.; Heidari, A.H.; Guddati, M.; Kim, Y.R. Experimental Investigation of Anisotropy in Asphalt Concrete. Transp. Res. Rec. J. Transp. Res. Board 2005, 1929, 238–247. [Google Scholar] [CrossRef]
  23. Darabi, M.K.; Abu Al-Rub, R.K.; Masad, E.A.; Huang, C.-W.; Little, D.N. A modified viscoplastic model to predict the permanent deformation of asphaltic materials under cyclic-compression loading at high temperatures. Int. J. Plast. 2012, 35, 100–134. [Google Scholar] [CrossRef]
  24. Darabi, M.K.; Abu Al-Rub, R.K.; Masad, E.A.; Little, D.N. Cyclic Hardening-Relaxation Viscoplasticity Model for Asphalt Concrete Materials. J. Eng. Mech. 2013, 139, 832–847. [Google Scholar] [CrossRef]
  25. Karimi, M.M.; Tabatabaee, N.; Jahangiri, B.; Darabi, M.K. Constitutive modeling of hardening-relaxation response of asphalt concrete in cyclic compressive loading. Constr. Build. Mater. 2017, 137, 169–184. [Google Scholar] [CrossRef]
  26. Blanc, J.; Gabet, T.; Hornych, P.; Piau, J.-M.; Di Benedetto, H. Cyclic triaxial tests on bituminous mixtures. Road Mater. Pavement Des. 2014, 16, 46–69. [Google Scholar] [CrossRef] [Green Version]
  27. Airey, G.D.; Rahimzadeh, B.; Collop, A.C. Viscoelastic linearity limits for bituminous materials. Mater. Struct. 2003, 36, 643–647. [Google Scholar] [CrossRef]
  28. Zhou, F.; Scullion, T.; Sun, L. Verification and Modeling of Three-Stage Permanent Deformation Behavior of Asphalt Mixes. J. Transp. Eng. 2004, 130, 486–494. [Google Scholar] [CrossRef]
  29. Zhou, F.; Scullion, T. Discussion: Three Stages of Permanent Deformation Curve and Rutting Model. Int. J. Pavement Eng. 2002, 3, 251–260. [Google Scholar] [CrossRef]
  30. Biligiri, K.P.; Kaloush, K.E.; Mamlouk, M.S.; Witczak, M.W. Rational Modeling of Tertiary Flow for Asphalt Mixtures. Transp. Res. Rec. J. Transp. Res. Board 2007, 2001, 63–72. [Google Scholar] [CrossRef]
  31. Ameri, M.; Sheikhmotevali, A.H.; Fasihpour, A. Evaluation and comparison of flow number calculation methods. Road Mater. Pavement Des. 2014, 15, 182–206. [Google Scholar] [CrossRef]
  32. Javilla, B.; Mo, L.; Hao, F.; An, S.; Wu, S. Systematic comparison of two-stage analytical rutting models of asphalt mixtures. Constr. Build. Mater. 2017, 153, 716–727. [Google Scholar] [CrossRef]
  33. Sousa, J.B.; Craus, J.; Monismith, C.L. Summary Report on Permanent Deformation in Asphalt Concrete; University of California: Berkeley, CA, USA, 1991. [Google Scholar]
  34. Uzan, J.; Motola, Y. Damage Evaluation in Simple Shear Tests with and without Stress Reversal of Asphalt Concrete. Road Mater. Pavement Des. 2006, 7, 71–86. [Google Scholar] [CrossRef]
  35. Li, Q.; Lee, H.J.; Lee, S.Y. Permanent Deformation Model Based on Shear Properties of Asphalt Mixtures: Development and Calibration. Transp. Res. Rec. 2011, 2210, 81–89. [Google Scholar] [CrossRef]
  36. Zieliński, P. Indirect tensile test as a simple method for rut resistance evaluation of asphalt mixtures—Polish experience. Road Mater. Pavement Des. 2020, 23, 112–128. [Google Scholar] [CrossRef]
  37. Perraton, D.; Di Benedetto, H.; Sauzéat, C.; De La Roche, C.; Bankowski, W.; Partl, M.; Grenfell, J. Rutting of bituminous mixtures: Wheel tracking tests campaign analysis. Mater. Struct. 2011, 44, 969–986. [Google Scholar] [CrossRef]
  38. Mohammad, L.N.; Elseifi, M.; Cao, W.; Raghavendra, A.; Ye, M. Evaluation of various Hamburg wheel-tracking devices and AASHTO T 324 specification for rutting testing of asphalt mixtures. Road Mater. Pavement Des. 2017, 18, 128–143. [Google Scholar] [CrossRef]
  39. Ahmed, A.W.; Erlingsson, S. Evaluation of a permanent deformation model for asphalt concrete mixtures using extra-large wheel-tracking and heavy vehicle simulator tests. Road Mater. Pavement Des. 2014, 16, 154–171. [Google Scholar] [CrossRef]
  40. Francken, L. Permanent deformation law of bituminous road mixes in repeated triaxial compression. In Proceedings of the 4th International Conference on Structural Design of Asphalt Pavements, Ann Arbor, MI, USA, 22–26 August 1977; pp. 483–496. [Google Scholar]
  41. Khedr Safwan, A. Deformation Mechanism in Asphaltic Concrete. J. Transp. Eng. 1986, 112, 29–45. [Google Scholar] [CrossRef]
  42. Tseng, K.-H.; Lytton, R.L. Prediction of Permanent Deformation in Flexible Pavement Materials. In Implication of Aggregates in the Design, Construction and Performance of Flexible Pavements; ASTM International: West Conshohocken, PA, USA, 1989; Volume 1016, pp. 154–172. [Google Scholar]
  43. Wilshire, B.; Evans, R.W. Acquisition and analysis of creep data. J. Strain Anal. Eng. Des. 1994, 29, 159–165. [Google Scholar] [CrossRef]
  44. Cerni, G.; Cardone, F.; Virgili, A.; Camilli, S. Characterisation of permanent deformation behaviour of unbound granular materials under repeated triaxial loading. Constr. Build. Mater. 2012, 28, 79–87. [Google Scholar] [CrossRef]
  45. Choi, Y.T.; Subramanian, V.; Guddati, M.N.; Kim, Y.R. Incremental Model for Prediction of Permanent Deformation of Asphalt Concrete in Compression. Transp. Res. Rec. 2012, 2296, 24–35. [Google Scholar] [CrossRef]
  46. Choi, Y.-T.; Kim, Y.R. Development of characterisation models for incremental permanent deformation model for asphalt concrete in confined compression. Road Mater. Pavement Des. 2013, 14, 266–288. [Google Scholar] [CrossRef]
  47. Yun, T.; Kim, Y.R. A viscoplastic constitutive model for hot mix asphalt in compression at high confining pressure. Constr. Build. Mater. 2011, 25, 2733–2740. [Google Scholar] [CrossRef]
  48. Subramanian, V.; Guddati, M.N.; Richard Kim, Y. A viscoplastic model for rate-dependent hardening for asphalt concrete in compression. Mech. Mater. 2013, 59, 142–159. [Google Scholar] [CrossRef] [Green Version]
  49. Levenberg, E.; Uzan, J. Triaxial Small-Strain Viscoelastic-Viscoplastic Modeling of Asphalt Aggregate Mixes. Mech. Time-Depend. Mater. 2004, 8, 365–384. [Google Scholar] [CrossRef]
  50. Di Benedetto, H.; Nguyen, Q.T.; Sauzéat, C. Nonlinearity, Heating, Fatigue and Thixotropy during Cyclic Loading of Asphalt Mixtures. Road Mater. Pavement Des. 2011, 12, 129–158. [Google Scholar] [CrossRef]
  51. Perl, M.; Uzan, J.; Sides, A. Visco-elasto-plastic constitutive law for a bituminous mixture under repeated loading. Transp. Res. Rec. 1983, 911, 20–26. [Google Scholar]
  52. Sides, A.; Uzan, J.; Perl, M. A comprehensive viscoelasto-plastic characterization of sand-asphalt compressive and tensile cyclic loading. J. Test. Eval. 1985, 13, 49–59. [Google Scholar]
  53. Cao, W.; Kim, Y.R. A viscoplastic model for the confined permanent deformation of asphalt concrete in compression. Mech. Mater. 2016, 92, 235–247. [Google Scholar] [CrossRef]
  54. Uzan, J. Permanent Deformation in Flexible Pavements. J. Transp. Eng. 2004, 130, 6–13. [Google Scholar] [CrossRef]
  55. Park, S.; Kim, Y. Interconversion between relaxation modulus and creep compliance for viscoelastic solids. J. Mater. Civ. Eng. 1999, 11, 76–82. [Google Scholar] [CrossRef]
  56. Mun, S.; Chehab, G.R.; Kim, Y.R. Determination of Time-domain Viscoelastic Functions using Optimized Interconversion Techniques. Road Mater. Pavement Des. 2007, 8, 351–365. [Google Scholar] [CrossRef]
  57. Collop, A.C.; Scarpas, A.; Kasbergen, C.; De Bondt, A. Development and Finite Element Implementation of Stress-Dependent Elastoviscoplastic Constitutive Model with Damage for Asphalt. Transp. Res. Rec. J. Transp. Res. Board 2003, 1832, 96–104. [Google Scholar] [CrossRef]
  58. Di Benedetto, H.; Mondher, N.; Sauzéat, C.; Olard, F. Three-dimensional Thermo-viscoplastic Behaviour of Bituminous Materials: The DBN Model. Road Mater. Pavement Des. 2007, 8, 285–315. [Google Scholar] [CrossRef]
  59. Olard, F.; Di Benedetto, H. General “2S2P1D” Model and Relation Between the Linear Viscoelastic Behaviours of Bituminous Binders and Mixes. Road Mater. Pavement Des. 2003, 4, 185–224. [Google Scholar] [CrossRef]
  60. Grzesikiewicz, W.; Wakulicz, A.; Zbiciak, A. Non-linear problems of fractional calculus in modeling of mechanical systems. Int. J. Mech. Sci. 2013, 70, 90–98. [Google Scholar] [CrossRef]
  61. Zbiciak, A. Mathematical description of rheological properties of asphalt-aggregate mixes. Bull. Pol. Acad. Sci. Tech. Sci. 2013, 61, 65–72. [Google Scholar] [CrossRef] [Green Version]
  62. Di Benedetto, H.; Olard, F.; Sauzéat, C.; Delaporte, B. Linear viscoelastic behaviour of bituminous materials: From binders to mixes. Road Mater. Pavement Des. 2004, 5, 163–202. [Google Scholar] [CrossRef]
  63. Kachanov, L.M. On time to rupture in creep conditions. Izv. Akad. Nauk SSSR 1958, 8, 26–31. (In Russian) [Google Scholar]
  64. Darabi, M.K.; Abu Al-Rub, R.K.; Masad, E.A.; Huang, C.-W.; Little, D.N. A thermo-viscoelastic–viscoplastic–viscodamage constitutive model for asphaltic materials. Int. J. Solids Struct. 2011, 48, 191–207. [Google Scholar] [CrossRef] [Green Version]
  65. Rabotnov, I.N.; Rabotnov, I.U.N.; Rabotnov, I.U.r.N.; Rabotnov, J.N.; Rabotnov, Y.N.; Leckie, F.A. Creep Problems in Structural Members; North-Holland Publishing Company: Amsterdam, The Netherlands, 1969; Volume 7. [Google Scholar]
  66. Perzyna, P. Thermodynamic Theory of Viscoplasticity. In Advances in Applied Mechanics; Elsevier: Amsterdam, The Netherlands, 1971; pp. 313–354. [Google Scholar]
  67. Rushing, J.F.; Darabi, M.K.; Rahmani, E.; Little, D.N. Comparing rutting of airfield pavements to simulations using Pavement Analysis Using Nonlinear Damage Approach (PANDA). Int. J. Pavement Eng. 2017, 18, 138–159. [Google Scholar] [CrossRef]
  68. Schapery, R.A. On the characterization of nonlinear viscoelastic materials. Polym. Eng. Sci. 1969, 9, 295–310. [Google Scholar] [CrossRef]
  69. Huang, C.-W.; Abu Al-Rub, R.K.; Masad, E.A.; Little, D.N.; Airey, G.D. Numerical implementation and validation of a nonlinear viscoelastic and viscoplastic model for asphalt mixes. Int. J. Pavement Eng. 2011, 12, 433–447. [Google Scholar] [CrossRef]
  70. Abu Al-Rub, R.K.; Darabi, M.K.; Kim, S.-M.; Little, D.N.; Glover, C.J. Mechanistic-based constitutive modeling of oxidative aging in aging-susceptible materials and its effect on the damage potential of asphalt concrete. Constr. Build. Mater. 2013, 41, 439–454. [Google Scholar] [CrossRef]
  71. Abu Al-Rub, R.K.; Darabi, M.K.; Little, D.N.; Masad, E.A. A micro-damage healing model that improves prediction of fatigue life in asphalt mixes. Int. J. Eng. Sci. 2010, 48, 966–990. [Google Scholar] [CrossRef]
  72. Shakiba, M.; Darabi, M.K.; Abu Al-Rub, R.K.; You, T.; Little, D.N.; Masad, E.A. Three-dimensional microstructural modelling of coupled moisture–mechanical response of asphalt concrete. Int. J. Pavement Eng. 2015, 16, 445–466. [Google Scholar] [CrossRef]
  73. Huang, C.-W.; Abu Al-Rub, R.K.; Masad, E.A.; Little, D.N. Three-Dimensional Simulations of Asphalt Pavement Permanent Deformation Using a Nonlinear Viscoelastic and Viscoplastic Model. J. Mater. Civ. Eng. 2011, 23, 56–68. [Google Scholar] [CrossRef] [Green Version]
  74. Darabi, M.K.; Kola, R.; Little, D.N.; Garg, N. Time-dependent Drucker-Prager-Cap model coupled with PANDA (Pavement Analysis Using Nonlinear Damage Approach) to predict rutting performance of flexible pavements. Constr. Build. Mater. 2020, 244, 118326. [Google Scholar] [CrossRef]
  75. Rahmani, E.; Darabi, M.K.; Abu Al-Rub, R.K.; Kassem, E.; Masad, E.A.; Little, D.N. Effect of confinement pressure on the nonlinear-viscoelastic response of asphalt concrete at high temperatures. Constr. Build. Mater. 2013, 47, 779–788. [Google Scholar] [CrossRef]
  76. Darabi, M.K.; Kola, R.; Little, D.N.; Rahmani, E.; Garg, N. Predicting Rutting Performance of Flexible Airfield Pavements Using a Coupled Viscoelastic-Viscoplastic-Cap Constitutive Relationship. J. Eng. Mech. 2019, 145, 04018129. [Google Scholar] [CrossRef]
  77. Behnke, R.; Wollny, I.; Hartung, F.; Kaliske, M. Thermo-mechanical finite element prediction of the structural long-term response of asphalt pavements subjected to periodic traffic load: Tire-pavement interaction and rutting. Comput. Struct. 2019, 218, 9–31. [Google Scholar] [CrossRef]
  78. Kim, S.-M.; Darabi, M.K.; Little, D.N.; Abu Al-Rub, R.K. Effect of the Realistic Tire Contact Pressure on the Rutting Performance of Asphaltic Concrete Pavements. KSCE J. Civ. Eng. 2018, 22, 2138–2146. [Google Scholar] [CrossRef]
  79. You, T.; Im, S.; Kim, Y.-R.; Little, D.N. Mechanistic Modeling to Evaluate Structural Performance of Bituminous Pavements with Inelastic Deformation and Fatigue Damage of Mixtures. J. Eng. Mech. 2017, 143, 04016126. [Google Scholar] [CrossRef]
  80. Abu Al-Rub, R.K.; Darabi, M.K.; Huang, C.-W.; Masad, E.A.; Little, D.N. Comparing finite element and constitutive modelling techniques for predicting rutting of asphalt pavements. Int. J. Pavement Eng. 2012, 13, 322–338. [Google Scholar] [CrossRef]
  81. Mulugeta Alamnie, M.; Taddesse, E.; Hoff, I. Thermo-piezo-rheological characterization of asphalt concrete. Constr. Build. Mater. 2022, 329, 127106. [Google Scholar] [CrossRef]
  82. Lai, J.; Bakker, A. 3-D Schapery representation for non-linear viscoelasticity and finite element implementation. Comput. Mech. 1996, 18, 182–191. [Google Scholar] [CrossRef]
  83. Darabi, M.K.; Abu Al-Rub, R.K.; Masad, E.A.; Little, D.N. Constitutive modeling of fatigue damage response of asphalt concrete materials with consideration of micro-damage healing. Int. J. Solids Struct. 2013, 50, 2901–2913. [Google Scholar] [CrossRef] [Green Version]
  84. Shakiba, M.; Darabi, M.K.; Abu Al-Rub, R.K.; Masad, E.A.; Little, D.N. Microstructural modeling of asphalt concrete using a coupled moisture–mechanical constitutive relationship. Int. J. Solids Struct. 2014, 51, 4260–4279. [Google Scholar] [CrossRef] [Green Version]
  85. Yideti, T.F.; Birgisson, B.; Jelagin, D.; Guarin, A. Packing theory-based framework to evaluate permanent deformation of unbound granular materials. Int. J. Pavement Eng. 2013, 14, 309–320. [Google Scholar] [CrossRef]
  86. Oda, M.; Nakayama, H. Yield Function for Soil with Anisotropic Fabric. J. Eng. Mech. 1989, 115, 89–104. [Google Scholar] [CrossRef]
  87. 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]
  88. Coleri, E.; Harvey, J.T.; Yang, K.; Boone, J.M. A micromechanical approach to investigate asphalt concrete rutting mechanisms. Constr. Build. Mater. 2012, 30, 36–49. [Google Scholar] [CrossRef]
  89. Chehab, G.R.; Seo, Y.; Kim, Y.R. Viscoelastoplastic damage characterization of asphalt-aggregate mixtures using digital image correlation. Int. J. Geomech. 2007, 7, 111–118. [Google Scholar] [CrossRef]
  90. Shojaeifard, M.; Baghani, M.; Shahsavari, H. Rutting investigation of asphalt pavement subjected to moving cyclic loads: An implicit viscoelastic–viscoplastic–viscodamage FE framework. Int. J. Pavement Eng. 2018, 21, 1393–1407. [Google Scholar] [CrossRef]
  91. Underwood, B.S. A continuum damage model for asphalt cement and asphalt mastic fatigue. Int. J. Fatigue 2016, 82, 387–401. [Google Scholar] [CrossRef]
  92. Pellinen, T.K.; Christensen, D.W.; Rowe, G.M.; Sharrock, M. Fatigue-Transfer Functions: How Do They Compare? Transp. Res. Rec. J. Transp. Res. Board 2004, 1896, 77–87. [Google Scholar] [CrossRef]
  93. Ling, M.; Luo, X.; Chen, Y.; Gu, F.; Lytton, R.L. Mechanistic-empirical models for top-down cracking initiation of asphalt pavements. Int. J. Pavement Eng. 2018, 21, 464–473. [Google Scholar] [CrossRef]
  94. Zhang, Y.; Luo, X.; Luo, R.; Lytton, R.L. Crack initiation in asphalt mixtures under external compressive loads. Constr. Build. Mater. 2014, 72, 94–103. [Google Scholar] [CrossRef] [Green Version]
  95. Dyskin, A.V.; Sahouryeh, E.; Jewell, R.J.; Joer, H.; Ustinov, K.B. Influence of shape and locations of initial 3-D cracks on their growth in uniaxial compression. Eng. Fract. Mech. 2003, 70, 2115–2136. [Google Scholar] [CrossRef]
  96. Onifade, I.; Birgisson, B.; Balieu, R. Energy-based damage and fracture framework for viscoelastic asphalt concrete. Eng. Fract. Mech. 2015, 145, 67–85. [Google Scholar] [CrossRef]
  97. Widyatmoko, I.; Ellis, C.; Read, J.M. Energy dissipation and the deformation resistance of bituminous mixtures. Mater. Struct. 1999, 32, 218–223. [Google Scholar] [CrossRef]
  98. Ghuzlan, K.A.; Carpenter, S.H. Energy-Derived, Damage-Based Failure Criterion for Fatigue Testing. Transp. Res. Rec. J. Transp. Res. Board 2000, 1723, 141–149. [Google Scholar] [CrossRef]
  99. Luo, X.; Luo, R.; Lytton, R.L. Energy-based mechanistic approach for damage characterization of pre-flawed visco-elasto-plastic materials. Mech. Mater. 2014, 70, 18–32. [Google Scholar] [CrossRef]
  100. Di Benedetto, H.; Partl, M.N.; Francken, L.; De La Roche Saint André, C. Stiffness testing for bituminous mixtures. Mater. Struct. 2001, 34, 66–70. [Google Scholar] [CrossRef]
  101. Skelton, R.P.; Gandy, D. Creep-fatigue damage accumulation and interaction diagram based on metallographic interpretation of mechanisms. Mater. High Temp. 2008, 25, 27–54. [Google Scholar] [CrossRef]
  102. Chaboche, J.L. Constitutive equations for cyclic plasticity and cyclic viscoplasticity. Int. J. Plast. 1989, 5, 247–302. [Google Scholar] [CrossRef]
  103. Barbera, D.; Chen, H.; Liu, Y. On Creep Fatigue Interaction of Components at Elevated Temperature. J. Press. Vessel Technol. 2016, 138, 041403. [Google Scholar] [CrossRef] [Green Version]
  104. Liu, H.; Yang, X.; Jiang, L.; Lv, S.; Huang, T.; Yang, Y. Fatigue-creep damage interaction model of asphalt mixture under the semi-sine cycle loading. Constr. Build. Mater. 2020, 251, 119070. [Google Scholar] [CrossRef]
  105. Gupta, R.; Atul Narayan, S.P. Tensile creep of asphalt concrete in repeated loading tests and its effect on energy dissipation. Int. J. Pavement Eng. 2019, 22, 927–939. [Google Scholar] [CrossRef]
  106. Zhang, Y.; Luo, R.; Lytton, R.L. Characterizing Permanent Deformation and Fracture of Asphalt Mixtures by Using Compressive Dynamic Modulus Tests. J. Mater. Civ. Eng. 2012, 24, 898–906. [Google Scholar] [CrossRef]
  107. Chaboche, J.L.; Lesne, P.M. A Non-Linear Continuous Fatigue Damage Model. Fatigue Fract. Eng. Mater. Struct. 1988, 11, 1–17. [Google Scholar] [CrossRef]
  108. Cheng, G. A fatigue damage accumulation model based on continuum damage mechanics and ductility exhaustion. Int. J. Fatigue 1998, 20, 495–501. [Google Scholar] [CrossRef]
  109. Lemaitre, J.; Plumtree, A. Application of Damage Concepts to Predict Creep-Fatigue Failures. J. Eng. Mater. Technol. 1979, 101, 284–292. [Google Scholar] [CrossRef]
Figure 1. Flow chart of the study structure.
Figure 1. Flow chart of the study structure.
Materials 15 03480 g001
Figure 2. (a) Schematic three-stage creep deformation and (b) example of a creep-recovery test result at different confining pressures (100 kPa and 200 kPa).
Figure 2. (a) Schematic three-stage creep deformation and (b) example of a creep-recovery test result at different confining pressures (100 kPa and 200 kPa).
Materials 15 03480 g002
Figure 3. Example—permanent deformation prediction of some models using test data: (a) logarithmic scale and (b) normal scale.
Figure 3. Example—permanent deformation prediction of some models using test data: (a) logarithmic scale and (b) normal scale.
Materials 15 03480 g003
Figure 4. Stress–strain state in compression: (a) hysteresis loop under creep-recovery and (b) crosshead strain rate test.
Figure 4. Stress–strain state in compression: (a) hysteresis loop under creep-recovery and (b) crosshead strain rate test.
Materials 15 03480 g004
Figure 5. Schematic—typical creep-recovery strain components.
Figure 5. Schematic—typical creep-recovery strain components.
Materials 15 03480 g005
Figure 6. (a) Maxwell’s viscoelastic model, (b) Burgers viscoelastic model, (c) viscoelastic–plastic model with slider, and (d) elastic-visco-plastic nonlinear fractional (spring-pot) model.
Figure 6. (a) Maxwell’s viscoelastic model, (b) Burgers viscoelastic model, (c) viscoelastic–plastic model with slider, and (d) elastic-visco-plastic nonlinear fractional (spring-pot) model.
Materials 15 03480 g006
Figure 7. (a) The 2S2P1D rheological model (h and k are two parabolic creep elements), (b) DBN model for bituminous mixtures, and (c) slip device for the viscoplastic model with viscoelastic component.
Figure 7. (a) The 2S2P1D rheological model (h and k are two parabolic creep elements), (b) DBN model for bituminous mixtures, and (c) slip device for the viscoplastic model with viscoelastic component.
Materials 15 03480 g007
Figure 8. Schematic of representation of internal stress evolution during loading and unloading. Adapted with permission from [53]. 2016, Mechanics of Materials, Elsevier.
Figure 8. Schematic of representation of internal stress evolution during loading and unloading. Adapted with permission from [53]. 2016, Mechanics of Materials, Elsevier.
Materials 15 03480 g008
Figure 9. Schematic—the evolution of the effective viscoplastic strain and the hardening–relaxation during loading–unloading cycle. ( ε e vp q vp the driving force for viscoplastic softening or recovery in the viscoplastic hardening). Adapted with permission from [23], 2012, International Journal of Plasticity, Elsevier.
Figure 9. Schematic—the evolution of the effective viscoplastic strain and the hardening–relaxation during loading–unloading cycle. ( ε e vp q vp the driving force for viscoplastic softening or recovery in the viscoplastic hardening). Adapted with permission from [23], 2012, International Journal of Plasticity, Elsevier.
Materials 15 03480 g009
Figure 10. Development of the PANDA mechanistic permanent deformation model.
Figure 10. Development of the PANDA mechanistic permanent deformation model.
Materials 15 03480 g010
Figure 11. Two-dimensional particles orientation for micromechanics modeling (orientation,   θ and vector, m).
Figure 11. Two-dimensional particles orientation for micromechanics modeling (orientation,   θ and vector, m).
Materials 15 03480 g011
Figure 12. Effect of anisotropy on the stress–strain relationship predicted by the model. Reprinted with permission from [87]. 2005, International Journal of Plasticity, Elsevier.
Figure 12. Effect of anisotropy on the stress–strain relationship predicted by the model. Reprinted with permission from [87]. 2005, International Journal of Plasticity, Elsevier.
Materials 15 03480 g012
Figure 13. Digital image analysis process for microstructure modeling.
Figure 13. Digital image analysis process for microstructure modeling.
Materials 15 03480 g013
Figure 14. The procedure followed for aggregate domain creation: (a) unprocessed 2D X-ray CT image used for the development of image stacks (left: before HVS trafficking and right: after HVS trafficking), (b) 3D images developed from X-ray CT image stacks (left: before HVS trafficking and right: after HVS trafficking), and (c) 3D aggregate volumes (left: before HVS trafficking and right: after HVS trafficking). Note: direction of HVS traffic is out of the page. Reprinted with permission from [88]. 2012, Construction and Building Materials, Elsevier.
Figure 14. The procedure followed for aggregate domain creation: (a) unprocessed 2D X-ray CT image used for the development of image stacks (left: before HVS trafficking and right: after HVS trafficking), (b) 3D images developed from X-ray CT image stacks (left: before HVS trafficking and right: after HVS trafficking), and (c) 3D aggregate volumes (left: before HVS trafficking and right: after HVS trafficking). Note: direction of HVS traffic is out of the page. Reprinted with permission from [88]. 2012, Construction and Building Materials, Elsevier.
Materials 15 03480 g014
Figure 15. (a) Damage domains of bituminous mixtures (ε—strain amplitude and N—number of cycles) and (b) fatigue and permanent deformation interactive damage for asphalt concrete (NC, NIC, Nf, and NIf failure lines for creep and fatigue domains).
Figure 15. (a) Damage domains of bituminous mixtures (ε—strain amplitude and N—number of cycles) and (b) fatigue and permanent deformation interactive damage for asphalt concrete (NC, NIC, Nf, and NIf failure lines for creep and fatigue domains).
Materials 15 03480 g015
Figure 16. (a) Cyclic fatigue hysteresis and (b) creep–fatigue interaction stress–strain loop.
Figure 16. (a) Cyclic fatigue hysteresis and (b) creep–fatigue interaction stress–strain loop.
Materials 15 03480 g016
Figure 17. Parametric study of creep–fatigue damage interaction coefficient.
Figure 17. Parametric study of creep–fatigue damage interaction coefficient.
Materials 15 03480 g017
Table 2. Summary of material parameters of the constitutive equations and physical meaning in Figure 10.
Table 2. Summary of material parameters of the constitutive equations and physical meaning in Figure 10.
ParameterPhysical MeaningTheory/DomainCalibration Test
Viscoelastic Model Parameters [68,82]
D o Instantaneous creep compliance. Characterizes the instantaneous elastic part of the viscoelastic strain Linear viscoelasticCreep compliance or
Dynamic Modulus
λ n nth retardation time. A measure of the required time for the viscoelastic material to relax the induced stress
D n nth coefficient of the Prony series associated with the nth retardation time kn. These parameters characterize the transient compliance of the material
g o , g 1 , g 2 Nonlinear viscoelastic parameters, where g o measures the reduction or the increase in the instantaneous compliance; g 1 defines the nonlinearity effects in the transient compliance; and g 2 is the nonlinear parameter accounting for the loading rate effects on the creep responseSchapery’s Nonlinear viscoelasticCreep-recovery
Viscoplastic Model Parameters [64,66,69]
Γ v p Viscoplastic fluidity parameter, such that 1 / Γ v p   is a measure of viscoplastic viscosityPerzyna’s Viscoplastic (with Drucker–Prager yield criteria)Creep, creep-recovery (creep part), constrain strain rate test
κ o , κ 1 , κ 2 Isotropic hardening parameters, where κ o defines the initial yield strength; κ o + κ 1 defines the saturated limit of the hardening function; and κ 2 defines the hardening rate and controls the shape of the hardening function versus the effective viscoplastic strain ( ε ˙ e v p )
N Viscoplastic rate sensitivity exponent and describes the nonlinear rate dependency of viscoplastic response
α   , β Govern the pressure sensitivity of the yield surface and plastic potential functions. Related to the angle of friction in the asphalt concrete
d vp Model parameter distinguishing viscoplastic responses in extension and contraction modes of loading
Visco-damage Model Parameters [64]
Γ v d Visco-damage fluidity parameter ( 1 / Γ v d is a measure of damage viscosity)CDM (Tertiary creep)Creep tests
q Stress dependency parameter. Defines the sensitivity of the damage evolution due to stress level
Y o The reference damage force obtained at a reference stress of a creep test
k Strain exponent parameter. Defines the sensitivity of the damage evolution due to strain level
d vd Model parameter distinguishing visco-damage responses in extension and contraction modes of loading
Micro-damage healing model Parameters [83]
Γ h Micro-damage healing fluidity parameter, such that 1 / Γ h is a measure of healing viscosityCDM for healing
k 1 , k 2 Healing model parameters that describe the effect of the damage and healing histories on the healing evolution
Temperature coupling terms parameters [64]
θ 1 , θ 2 , θ 3 Temperature sensitivity model parameters for viscoplastic, viscodamage and microdamage healing, respectivelyTime–temperature superposition Dynamic Modulus, Creep compliance
T o Reference temperature
Hardening–relaxation Model Parameters [24,25]
Γ h r The hardening–relaxation fluidity parameter, such that 1/ Γ h r is the hardening–relaxation retardation time controlling the rate of the hardening–relaxation Creep-recovery (at different rest periods)
S 1 , S 2 , S 3 Hardening–relaxation rate-sensitivity parameters that describe the relaxation behavior of viscoplastic hardening due to recoverable potential during rest period
Moisture damage Model Parameters [84]
p Adhesion or cohesion moisture damage parameterCDM for moisture damagePull-off test
δ Parameter that describes the moisture damage history
Table 3. Advantages and limitations of the PANDA Model.
Table 3. Advantages and limitations of the PANDA Model.
AdvantagesDisadvantages
The approach is based on fundamental theories of mechanics (mechanistic)
Unlimited modeling capacity
The material properties and temperature coupling are integrated
The approach has the capacity to couple different damage types (such as moisture, fatigue cracking, etc.) with permanent deformation
It enables realistic rutting prediction with full pavement structure modeling using finite element method (3D modeling)
Moving loads can be modeled, which was not possible in traditional methods
etc.
Complex constitutive equations
Large number of modeling parameters
It requires many different calibration tests
It is computationally and experimentally costly
It is at development stage and some theories have limitations (e.g., classic viscoplasticity theory)
The numerical implementation is based on user material (UMAT) model (not standalone)
etc.
Table 4. The commonly used FEM simulation loading assumptions.
Table 4. The commonly used FEM simulation loading assumptions.
Loading Approach
ModeLoad ConfigurationPulse
Loading
Equivalent
Loading
Moving
Loading
2DSingle wheel (plane strain)
Single wheel (axisymmetric)
Single wheel (moving loading)
Multiple wheel (moving loading)
3DSingle wheel (rectangular)
Whole wheel path
Single wheel (circular)
Single wheel (moving)
Multiple wheel
FEM simulation applied.
Table 5. Summary of permanent deformation modeling approaches.
Table 5. Summary of permanent deformation modeling approaches.
Approach
Properties, Theories, MethodsPANDAMicrostructure AnalogicalEmpirical
Continuum damage mechanics--
Viscoelastic-
Viscoplasticity
Micromechanics--
Finite element simulation-
Time–temperature superposition-
Coupling (healing, moisture)--
Coupling fatigue–rutting damage----
Full pavement deformation model--
Hardening–relaxation mechanism--
Uniaxial and triaxial repeated load test
√ Incorporated, † Not yet incorporated.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alamnie, M.M.; Taddesse, E.; Hoff, I. Advances in Permanent Deformation Modeling of Asphalt Concrete—A Review. Materials 2022, 15, 3480. https://doi.org/10.3390/ma15103480

AMA Style

Alamnie MM, Taddesse E, Hoff I. Advances in Permanent Deformation Modeling of Asphalt Concrete—A Review. Materials. 2022; 15(10):3480. https://doi.org/10.3390/ma15103480

Chicago/Turabian Style

Alamnie, Mequanent Mulugeta, Ephrem Taddesse, and Inge Hoff. 2022. "Advances in Permanent Deformation Modeling of Asphalt Concrete—A Review" Materials 15, no. 10: 3480. https://doi.org/10.3390/ma15103480

APA Style

Alamnie, M. M., Taddesse, E., & Hoff, I. (2022). Advances in Permanent Deformation Modeling of Asphalt Concrete—A Review. Materials, 15(10), 3480. https://doi.org/10.3390/ma15103480

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