# Finite Element Analysis on the Creep Constitutive Equation of High Modulus Asphalt Concrete

Authors: Xiu-shan Wang; Tao-tao Fan; Lum Kitmeng
Journal: Advances in Materials Science and Engineering
 (2015)

---

## Abstract

In order to obtain the viscoelastic constitutive relation of High Modulus Asphalt Concrete (HMAC), the constitutive relationship is used to carry on the various numerical calculation, utilize integral transforming in its Bailey-Norton creep rule, and establish new practical model. Combining with the collection data from creep tests of uniaxial compression, regression calculation is carried out with 1STOPT fitting software; thus the creep parameters of High Modulus Asphalt Concrete in different temperatures were obtained, which can be used in HMAC creep model.


---

## Body

## 1. Introduction

High Modulus Asphalt Concrete (HMAC) is originally a kind of asphalt concrete mixture agitated by hard asphalt, a certain grading of stone, and additives. Under 15°C, 10 Hz test conditions, its dynamic modulus can reach over 14000 MPa [1–6], which fully exhibited such advantages as high modulus, good rutting resistance, low temperature cracking, and low thermal fatigue cracking sensitivity. As the viscoelastic material, HMAC has unique constitutive relation different from other elastic, elastic-plastic materials; thus it is more difficult to grasp its mechanical characteristics. However, the development of viscoelastic theory provides mechanical analysis of HMAC with an effective tool finite element method. Since there is no complex analytic derivation in the finite element method, it is convenient to simulate viscoelastic of road materials and analyze the actual stress state of road surface [7–11].Due to the simplicity and high controllability of creep test method, many researchers have conducted a lot of laboratory tests to reveal the actual stress condition of road surface through the asphalt concrete material stress-strain relationship. But it is still very difficult to obtain viscoelastic parameters from these asphalt stress-strain data and then use computer software to simulate its mechanical characteristics. There are two reasons: one is the low precision degree of data fitting; the other is the inconsistency of formula in data fitting and computer simulation software. Therefore, though many researchers have viscoelastic data, they adopt elastic or elastic-plastic method instead to calculate the stress state of the approaching road surface [12, 13]. In order to better utilize these data, this paper analyzed the creep model in ABAQUS suitable for high modulus asphalt and carried on integration processing on Bailey-Norton creep law to establish new practical model. Based on the numerical simulation of uniaxial compression creep test, the creep parameters of High Modulus Asphalt Concrete in different temperatures were obtained by regression calculation with software 1STOPT of high efficiency and user-friendly control.


## 2. The Creep Constitutive Equations of Asphalt Concrete

The finite element analysis on material nonlinear problem has two aspects, namely, the establishment of the constitutive formula and the solution to nonlinear equations. Material constitutive relations are usually divided into two categories: the full amount of the constitutive equations and the incremental constitutive equations. Since the plastic deformation is unrecoverable inelastic deformation, stress state cannot be determined by the current state of deformation but rather by the route and deformation history. To accommodate this situation, the incremental constitutive equations are usually adopted in finite element analysis.Generally, creep effect is time-dependent; that is, under constant load conditions, the deformation of material increases with time. With constant load and displacement, creep effect has two stages: redistribution of stress within the structure as the first stage and the steady stress state of the structure as the second. In the process from transient creep state to steady creep state, the first stage is named as transient creep stage and the second stage is named as steady state creep stage. The steady creep stage can be analyzed by the full amount of the finite element analysis. But, for the transient creep, the load of a given change over time, and the structural creep effect under displacement, an incremental finite element method analysis is used usually along with thermal elastic-plastic incremental analysis, known as thermal elastic-plastic creep analysis.Asphalt concrete, which is the material dependent on time, temperature, and stress, presents deformation response in elasticity, plasticity, viscoelasticity, and viscoplasticity under repeated loads. The strain parameters can be expressed as(1)
ε
=
ε
e
+
ε
p
+
ε
ve
+
ε
vp
.In (1),  ε is the total strain varying with time, ε
e is elastic strain (recoverable and independent of time), ε
p is plastic strain (unrecoverable and independent of time), ε
ve is viscoelastic strain (recoverable and time-dependent), and ε
v
p is viscoplastic strain (unrecoverable and time-dependent).Obviously, such plastic nature of the asphalt mixture as plastic strain and viscoplastic strain can generate permanent deformation, and those plastic strains are cumulative under repeated loads. So only viscoplastic strain remains to be calculated. However, it is very hard to differentiate viscoelastic part and viscoplastic part especially when the two are changeable under repeated loads. On the other hand, creep tests can easily identify elastic strain and inelastic strain, that is, creep strain, which should include viscoplastic strain and a part of viscoelastic strain at a certain point of time. Since the current test method can only measure the combined effects of the two, creep and plasticity cannot be treated respectively.The creep deformation of the material can be expressed as a function of temperatureT, stress σ, and time t; that is, ε
c
=
f
(
T
,
σ
,
t
) which can be used to analyze the creep deformation [14–17]. Therefore, Bailey-Norton creep model in ABAQUS can be used to simulate the nonlinearity of the asphalt concrete layer. Equation (2) is expressed in the form of creep strain rate, as follows: (2)
ε
c
r
=
C
1
σ
C
2
t
C
3
e
-
C
4
/
T
.In (2), ε
c
r is creep strain rate; e is 2.718; t is time; T is temperature; C
1, C
2, C
3, and C
4 are material parameters.In order to simplify the model and enhance its usability, the Bailey-Norton creep law was retransformed by carrying on integral calculation on both sides simultaneously. The results are as follows:(3)
∫
ε
c
ε
c
r
d
ε
c
r
=
∫
0
t
C
1
σ
C
2
t
C
3
e
-
C
4
/
T
d
t
.When the temperature is constant, C
4
=
0, and (3) can be converted into(4)
ε
c
r
=
C
1
σ
C
2
t
C
3
+
1
C
3
+
1
.Equation (4) is the time hardening creep model (expressed by creep rate) in the converted ABAQUS; C
1, C
2, C
3 are model parameters dependent on temperature, which can be determined by material testing.


## 3. Creep Tests on High Modulus Asphalt Concrete

Since the High Modulus Asphalt Concrete is mainly used in sections of large traffic volume and harsh stress environment, grade A asphalt in specification [18] was chosen as asphalt material in the High Modulus Asphalt Concrete. In this study, two kinds of asphalt 50# and 70#, respectively, from Zhonghai Company were chosen based on the market survey, and five kinds of asphalt 70#-1 with additive from French PR company were chosen according to its mixed recommended dosage. They are 50#-1, 50#-2, 70#-1, 70#-2, and 70#-3. The performances of them can all meet the norms. The optimum asphalt content is shown in Table 1. Aggregate grading was selected from specification AC-20 median grading and the grading composition is shown in Table 2.Table 1
The optimum asphalt content.
Mixture type
50#-1
50#-2
70#-1
70#-2
70#-3
Best oilstone ratio
4.1
4.4
4.0
4.1
4.3Table 2
Grading composition in AC-20 median asphalt concrete.
Sieve sizegrading composition
Pass rate of hole sieve (%)
31.5
26.5
19
16
13.2
9.5
4.75
2.36
1.18
0.6
0.3
0.15
0.075
AC-20
100
100
95
85
71
61
41
30
22.5
16
11
8.5
5Use MTS810 material testing machine imported from America to conduct the test with uniaxial compression creep method. The specimen size was 100 mm in diameter, its height was 100 mm, and test temperature was 20°C, 40°C, and 60°C, respectively. A series of pressure (0.1 MPa, 0.2 MPa, 0.3 MPa, 0.4 MPa, and 0.5 MPa) was loaded on the specimen for 60 minutes and then unloaded for 10 minutes before starting processing data. Through these tests, the changes of vertical accumulative strain over time for specimens of different type produced at different temperatures during creep process were obtained. Due to space limitations, only two graphs are presented here with Figure1 showing creep curves for different types of asphalt mixture at temperature 40°C and Figure 2 showing the creep curves of 70#-1 asphalt mixture at different temperatures.Figure 1
Creep curves for different types of asphalt mixture.Figure 2
Creep curves of 70#-1 asphalt at different temperatures.From Figure1, it can be seen that, at loading stage, the instantaneous deformation of the specimen first occurs under the initial load, and then, with continuing role of dead load, the specimen deformation keeps increasing until the deformation increments are steady. After unloading, the elastic deformation resumes immediately, viscoelastic deformation gradually recovers over time, and plastic deformation remains deformed permanently because of its irreversibility. Asphalt mixtures present the following rule according to their cumulative creep strains: 70#-3 < 50#-2 < 70#-2 < 70#-1 < 50#-1. The cumulative strain of asphalt concrete 50#-1 is the maximum, indicating its weakest capacity to resist deformation under sustained loads; concrete 70#-3, mixed with 0.7% additive, has the minimum cumulative strain, indicating its strongest capacity to resist high temperature deformation.From Figure2, it can be seen that, under the same stress conditions, as the temperature increases, both the accumulated deformation and the residual deformation after relaxation are increased. At the high temperature of 40°C and 60°C, instantaneous deformation of High Modulus Asphalt Concrete with 0.7% additive is consistent with the ordinary asphalt concrete, while the cumulative deformation of the former is far less than the latter with the growth of loading time, indicating a significant improvement in high temperature deformation assistance of the asphalt mixture with admixture.


## 4. Creep Model Parameters Results

According to the creep performance test results of High Modulus Asphalt Concrete under different temperature conditions as well as treatment and regression conducted by professional fitting software 1STOPT, creep parametersC
1, C
2, and C
3 of 5 different types of asphalt concrete are obtained, respectively, as shown in Table 3.Table 3
Creep parameters for asphalt concrete.
Mixture type
Creep parameters
C
1
C
2
C
3
20°C
40°C
60°C
20°C
40°C
60°C
20°C
40°C
60°C
50#-1
2.98e − 15
1.06e − 16
1.85e − 21
2.815
3.101
3.860
−0.809
−0.820
−0.914
50#-2
6.23e − 16
2.22e − 17
1.27e − 22
1.808
2.369
2.898
−0.646
−0.696
−0.631
70#-1
8.67e − 16
3.09e − 17
5.39e − 22
2.026
2.655
3.248
−0.690
−0.712
−0.670
70#-2
6.52e − 16
2.33e − 17
4.06e − 22
1.921
2.305
2.469
−0.690
−0.663
−0.670
70#-3
1.26e − 6
1.66e − 7
1.87e − 8
0.378
0.517
0.565
−0.820
−0.875
−0.720From Table3, the following can be seen.At the same temperature, creep model coefficientC
1 has the biggest amplitude of variation and exhibits the change in the order of magnitudes with additive, while partial stress index C
2 is reduced in five asphalt concrete creep parameters of different types. Time index C
3 is always negative, its absolute value being less than 1, and changes little with additive variation.As for creep parameters at different temperatures, whatever the material is, the amplitude of variation in the creep model coefficientC
1 is always the biggest with the series of change in the gap, and the higher the temperature is, the smaller the value of C
1 is. As for partial stress index C
2, its value increases along with the temperature increase. The time index C
3, always being −1 and 0, shows no significant change regardless of temperature change.At different temperatures, compared with the other four mixtures, the creep coefficient modelC
1 of High Modulus Asphalt Concrete (70#-3) shows increases in series, while the partial stress index C
2 decreases by 4 to 6 times.


## 5. Conclusions

(1)
By integration process of the time hardening creep model in ABAQUS, a new simplified creep model is established.
(2)
From laboratory High Modulus Asphalt Concrete creep test, it can be learned that the size of creep accumulative strain of asphalt mixtures shows the following rule: 70#-3 < 50#-2 < 70#-2 < 70#-1 < 50#-1. Among all the five types, high modulus asphalt mixture with added additive shows the strongest high temperature deformation resistance.
(3)
Creep parameters of different types of High Modulus Asphalt Concrete at different temperatures are obtained with the professional fitting software 1STOPT on the basis of numerical simulation of uniaxial compression creep test.



---

*Source: 860454-2015-10-11.xml*