# On the Viscoelastic Parameters of Gussasphalt Mixture Based on Modified Burgers Model: Deviation and Experimental Verification

Authors: Faxiang Xie; Dengjing Zhang; Ao Zhou; Bohai Ji; Lin Chen
Journal: Advances in Materials Science and Engineering
 (2017)

---

## Abstract

Viscoelasticity is an important characteristic of gussasphalt mixtures. The aim of this study is to find the correct viscoelastic material parameters of the novel gussasphalt applied in the 4th Yangtze River Bridge based on the modified Burgers model. This study firstly derives the explicit Prony series form of the shear relaxation modulus of viscoelastic material from Laplace transformation, to fulfill the parameter inputting requirements of commonly used finite element software suites. Secondly, a kind of uniaxial penetration creep experiment on the gussasphalt mixtures is conducted. By fitting the creep compliance, the viscoelastic parameters characterized by the modified Burgers model are obtained. And thirdly, based on the viscoelastic test data of asphalt mixtures, the Prony series formula derived in this study is verified through the finite element simulation. The comparison results of the relative errors between the finite element simulation and the theoretical calculation confirm the reliability of the Prony series formulas deduced in this research. And finally, a stress-correcting method is proposed, which can significantly improve the accuracy of model parameters identification and reduce the relative error between the finite element simulation and the experimental data.


---

## Body

## 1. Introduction

Due to many eminent advantages, asphalt pavements have become the main structural form of Chinese highway. As a kind of viscoelastic material, asphalt mixtures produce viscous flows while bearing heavy loads at high temperature, thus producing nonrecoverable residual deformations [1]. With increasingly growing axle loads, rutting has become one of the most serious types of damage of asphalt pavements. Therefore, it is of great importance for the pavement engineers to describe the viscoelastic properties of asphalt mixtures while predicting rut formation and depth [2, 3]. Meanwhile, the new formulation of gussasphalt mixture and novel structural layer has been adopted in the construction of the 4th Nanjing Yangtze River Bridge which is 6th largest suspension bridge in the world [4]. However, to the best of the authors’ knowledge, the viscoelastic properties of the novel asphalt pavement which could be used to predict the ruts have never been measured, which has motivated research.Much research has investigated viscoelastic models of asphalt mixtures, such as the generalized model, rheological model, and master curve model. Among these models, the rheological Burgers model is extensively mentioned for their convenience in practical applications. Xu and Zhu [5] proposed the modified Burgers model and carried out the rutting experiments to calculate the rut depth which verified the effectiveness of the method. Chen et al. [6] transformed the Burgers model to Prony series and used finite element software to calculate the creep compliance parameters and compared the results with Maxwell model. Behzadfar and Hatzikiriakos [7] studied the viscoelastic properties of bitumen via rheological testing, where they proposed their rheological constitutive models as well.For general finite element analyze software such as ANSYS or ABAQUS, in which, the commonest way of setting the parameters of viscoelastic materials is entering the parameters as a Prony series [8, 9], describing the shear relaxation modulus thereof. However, they are inconsistent with parametric forms of the viscoelastic model. Saboo and Kumar [10] investigated the strain creep response of varying asphalt binders via four-element Burgers model under different temperature and stress amplitudes. The data of uniaxial creep tests with nonlinear least squares method was applied to fit the creep compliance at different temperatures. Compared with the Burgers model, the modified Burgers model can better explain the long-term deformation characteristics of asphalt mixtures more properly. Scholars have extensively studied the transformations and methods behind the Prony series version of the shear modulus based on the Burgers model [11, 12], while surprisingly the transformation method of the modified Burgers model has not yet been reported. Therefore, exploring the conversion method of the Prony series of the modified Burgers model is of greater theoretical and practical significance as it can better analyze the viscoelasticity of asphalt mixtures by using developed finite element analysis software, which has also motivated this research.Different experimental methods had been applied to evaluate the creep performance of the asphalt concrete. Airey and Rahimzadeh [13] designed a special tension-compression apparatus to test the viscoelastic properties of the asphalt mixtures. Diaz [11] applied modified cyclic creep test to evaluate the creep performance of the cold mix asphalt. Yao et al. [14] utilized flexural beam test to present the dynamic modulus and flexural stiffness of different surfacing materials on steel bridge decks including gussasphalt, epoxy asphalt concrete, which exhibited the difference between the dynamic modulus and the flexural stiffness. Yan et al. [15] made use of bending beam rheometer to determine the time-varying temperature rheological properties of waste tire rubber and reclaimed low density polyethylene asphalt mixtures. In the monograph [16], Sun presented extensive finite element analysis and experimental data showing that a kind of direct penetration method can reflect the actual stress state of the pavement properly. The static strength parameters including the shear strength, friction angle, and cohesive strength are also deducted. Thus, the uniaxial penetration method should also be suitable to measure the creep performance of the gussasphalt mixture. Meanwhile, the new formulation of gussasphalt mixture and novel structural layer has been utilized in the construction of the 4th Nanjing Yangtze River Bridge which is 6th largest suspension bridge in the world [4]. However, the viscoelastic properties of the novel asphalt pavement have never been measured, which motivated the penetration experiments.In this study, we devote to obtain the correct viscoelastic parameters of gussasphalt mixture via direct penetration method. First, the Prony series transformation method of the modified Burgers model is deduced and a simple numerical simulation is adopted to demonstrate the correctness of the proposed transformation equations. Then, to explore the viscoelastic properties of the gussasphalt, uniaxial penetration creep experiments under different temperatures are carried out. In the end, the stress-correction procedure to obtain the viscoelastic parameters is proposed and verified. The experimental data and simulation data are compared which shows the correctness of the method in this paper.


## 2. The Explicit Prony Series Form of the Modified Burgers Model

### 2.1. Prony Series Derivation of the Modified Burgers Model

Asphalt mixture is a kind of viscoelastic material, which exhibits creep properties under various loads, where the recoverable and nonrecoverable creeps of hardening mixture materials can be described by using Kelvin model and Maxwell model. The Burgers model is connected by the Kelvin model and the Maxwell model in series with four-unit-four-parameters to predict the creep behaviors of mixtures. Based on the model, Xu and Zhu [5], proposed a five-unit-four-parameters viscoelastic model, as shown in Figure 1, which is called the modified Burgers model.Figure 1
Modified Burgers model.Creep compliance is defined as the strain function changing with time under the effect of unit stress. According to the constitutive model in Figure1, the creep compliance of the modified Burgers model is expressed as(1)Jt=1E1+1ab1-e-bt+1E21-e-E2/η2t.In general finite element analysis software, the Prony series of the shear relaxation modulus of relevant materials is generally input by defining the properties of nonlinear viscoelastic materials. The form of the Prony series used for the shear relaxation modulus is(2)Gt=G0+G1∑inaie-t/ti,where G1, ai, and ti represent the initial shear modulus, the relative shear modulus, and the relaxation time, respectively.In the expression of creep compliance equation (1), the five viscoelastic parameters cannot be directly accessed as the Prony series, while relevant viscoelastic parameters should be converted into the parameters relating to the Prony parameters of the shear relaxation modulus.The three-dimensional (3D) form is utilized to represent (1):(3)Jt=121G1+1AB1-e-Bt+1G21-e-G2/n2t,where G1 and G2 stand for the shear moduli corresponding to E1 and E2 and n2 indicates the 3D viscosity coefficient corresponding to η2. In addition, A and B corresponding to a and b represent the 3D parameters of the modified Burgers model.When the viscoelastic volume deformation is not considered, the shear modulus shows the following relationships with the elastic modulus and the five coefficients are(4)G1=E121+μ,G2=E23,n2=η23,A=a3,B=b.Concerning the convolution relationship between the creep complianceJ(t) and the relaxation modulus Y(t) of viscoelastic materials, it can be found that [12](5)∫0tJτYt-τdτ=∫0tJτYt-τdτ=t.The Laplace transformation is applied to both sides of (5), leading to(6)Js¯∗Ys¯=1s2,where Js¯ and Ys¯ are Laplace transformation of J(s) and Y(s).It can be seen, from the Laplace transformation of (3), that(7)Js¯=121G1s+1Ass+B+1sn2s+G2.Then, the Laplace transformation of the relaxation modulus is(8a)Ys¯=1Js¯∗s2=2s1G1+1As+B+1n2s+G2.To prepare the system for inverse Laplace transformation, formula (8a) is rewritten as(8b)Ys¯=2∗q1s+q2s+α+q3s+β,where the parameters are (9)α,β=-12∗-ABn2+G2+G1n2+G1AAn2±ABn2+G2+G1n2+G1AAn22-4∗ABG2+G1G2+G1ABAn2,q1=ABG1G2ABG2+G1G2+G1AB,q2=N-αMβ-α,q3=βM-Nβ-αand the inner parametersM and N are(10)M=q2+q3=G12AB+G2ABG1+G2+G1G2,N=q2β+q3α=AB2n2+G22G12n2ABG1+G2+G1G2.An inverse Laplace transformation is performed to the both sides of (8b), which states that(11)Yt=2∗q1+q2e-αt+q3e-βt.Then the shear modulusG(t) can be obtained:(12)Gt=12Yt=q1+q2e-αt+q3e-βt.In this way, the parameters of Prony series of materials require to be input to the chosen finite element software as(13)a1=q2G1,t1=1α,a2=q3G1,t2=1β.From (13), the explicit parameters for Prony series of the modified Burgers model are obtained, which can be utilized directly to the finite element soft suites such as ANSYS.


### 2.2. Verification of Prony Series Transformation Equation

To verify the correctness of the derived results of the Prony series of the modified Burgers model, this study establishes the ANSYS finite element model, as shown in Figure2. The SOLID185 element with the size being 1 mm ∗ 1 mm ∗ 1 mm is used, and there is totally 1 element with 8 nodes. The material properties are taken from the properties of the Modified Burgers model in Table 1. The boundary conditions are shown in Figure 2, and the uniform load of 0.7 MPa is applied to the top surface of the element.Table 1
Material properties of the calculation example.
Parameters of modified Burgers model
Parameters of Prony series (shear modulus)
E
1 (MPa)
E
2 (MPa)
η
2 (Mpa∗s)
a (Mpa∗s)
b (/s)
a
1
t
1
a
2
t
2
87.04
20.87
5.70 × 104
1.44 × 104
5.04 × 10−3
0.2065
980.9978
0.6592
74.1388Figure 2
Schematic diagram of the ANSYS model.The nodal strain data of the nodes obtained by the ANSYS simulation are transformed into the creep compliance and compared with the theoretical calculation results of the modified Burgers model, as demonstrated in Table2. It can be found that the finite element simulation results are consistent with the theoretical calculation results, and the relative error is less than 0.70%. So the simulation results show that the Prony series results are reliable. By employing this method, the material parameters of the viscoelastic material can be accurately entered to the finite element software. However, the numerical simulation unit force is uniform, which has difference with the actual creep penetration experiments.Table 2
The comparison between equation calculation results and the ANSYS simulation.
Time (s)
Creep compliance of equation calculation (MPa−1)
Strain of ANSYS
Creep compliance of ANSYS (MPa−1)
Relative error (%)
1
1.16
E
-
02
8.07
E
-
03
1.15
E
-
02
-
0.37
101
1.87
E
-
02
1.32
E
-
02
1.88
E
-
02
0.50
201
2.37
E
-
02
1.67
E
-
02
2.38
E
-
02
0.67
301
2.72
E
-
02
1.92
E
-
02
2.74
E
-
02
0.69
601
3.41
E
-
02
2.40
E
-
02
3.43
E
-
02
0.60
901
3.86
E
-
02
2.72
E
-
02
3.88
E
-
02
0.54
1801
4.84
E
-
02
3.40
E
-
02
4.86
E
-
02
0.43
2701
5.54
E
-
02
3.89
E
-
02
5.56
E
-
02
0.36
3601
6.04
E
-
02
4.24
E
-
02
6.05
E
-
02
0.30


## 2.1. Prony Series Derivation of the Modified Burgers Model

Asphalt mixture is a kind of viscoelastic material, which exhibits creep properties under various loads, where the recoverable and nonrecoverable creeps of hardening mixture materials can be described by using Kelvin model and Maxwell model. The Burgers model is connected by the Kelvin model and the Maxwell model in series with four-unit-four-parameters to predict the creep behaviors of mixtures. Based on the model, Xu and Zhu [5], proposed a five-unit-four-parameters viscoelastic model, as shown in Figure 1, which is called the modified Burgers model.Figure 1
Modified Burgers model.Creep compliance is defined as the strain function changing with time under the effect of unit stress. According to the constitutive model in Figure1, the creep compliance of the modified Burgers model is expressed as(1)Jt=1E1+1ab1-e-bt+1E21-e-E2/η2t.In general finite element analysis software, the Prony series of the shear relaxation modulus of relevant materials is generally input by defining the properties of nonlinear viscoelastic materials. The form of the Prony series used for the shear relaxation modulus is(2)Gt=G0+G1∑inaie-t/ti,where G1, ai, and ti represent the initial shear modulus, the relative shear modulus, and the relaxation time, respectively.In the expression of creep compliance equation (1), the five viscoelastic parameters cannot be directly accessed as the Prony series, while relevant viscoelastic parameters should be converted into the parameters relating to the Prony parameters of the shear relaxation modulus.The three-dimensional (3D) form is utilized to represent (1):(3)Jt=121G1+1AB1-e-Bt+1G21-e-G2/n2t,where G1 and G2 stand for the shear moduli corresponding to E1 and E2 and n2 indicates the 3D viscosity coefficient corresponding to η2. In addition, A and B corresponding to a and b represent the 3D parameters of the modified Burgers model.When the viscoelastic volume deformation is not considered, the shear modulus shows the following relationships with the elastic modulus and the five coefficients are(4)G1=E121+μ,G2=E23,n2=η23,A=a3,B=b.Concerning the convolution relationship between the creep complianceJ(t) and the relaxation modulus Y(t) of viscoelastic materials, it can be found that [12](5)∫0tJτYt-τdτ=∫0tJτYt-τdτ=t.The Laplace transformation is applied to both sides of (5), leading to(6)Js¯∗Ys¯=1s2,where Js¯ and Ys¯ are Laplace transformation of J(s) and Y(s).It can be seen, from the Laplace transformation of (3), that(7)Js¯=121G1s+1Ass+B+1sn2s+G2.Then, the Laplace transformation of the relaxation modulus is(8a)Ys¯=1Js¯∗s2=2s1G1+1As+B+1n2s+G2.To prepare the system for inverse Laplace transformation, formula (8a) is rewritten as(8b)Ys¯=2∗q1s+q2s+α+q3s+β,where the parameters are (9)α,β=-12∗-ABn2+G2+G1n2+G1AAn2±ABn2+G2+G1n2+G1AAn22-4∗ABG2+G1G2+G1ABAn2,q1=ABG1G2ABG2+G1G2+G1AB,q2=N-αMβ-α,q3=βM-Nβ-αand the inner parametersM and N are(10)M=q2+q3=G12AB+G2ABG1+G2+G1G2,N=q2β+q3α=AB2n2+G22G12n2ABG1+G2+G1G2.An inverse Laplace transformation is performed to the both sides of (8b), which states that(11)Yt=2∗q1+q2e-αt+q3e-βt.Then the shear modulusG(t) can be obtained:(12)Gt=12Yt=q1+q2e-αt+q3e-βt.In this way, the parameters of Prony series of materials require to be input to the chosen finite element software as(13)a1=q2G1,t1=1α,a2=q3G1,t2=1β.From (13), the explicit parameters for Prony series of the modified Burgers model are obtained, which can be utilized directly to the finite element soft suites such as ANSYS.


## 2.2. Verification of Prony Series Transformation Equation

To verify the correctness of the derived results of the Prony series of the modified Burgers model, this study establishes the ANSYS finite element model, as shown in Figure2. The SOLID185 element with the size being 1 mm ∗ 1 mm ∗ 1 mm is used, and there is totally 1 element with 8 nodes. The material properties are taken from the properties of the Modified Burgers model in Table 1. The boundary conditions are shown in Figure 2, and the uniform load of 0.7 MPa is applied to the top surface of the element.Table 1
Material properties of the calculation example.
Parameters of modified Burgers model
Parameters of Prony series (shear modulus)
E
1 (MPa)
E
2 (MPa)
η
2 (Mpa∗s)
a (Mpa∗s)
b (/s)
a
1
t
1
a
2
t
2
87.04
20.87
5.70 × 104
1.44 × 104
5.04 × 10−3
0.2065
980.9978
0.6592
74.1388Figure 2
Schematic diagram of the ANSYS model.The nodal strain data of the nodes obtained by the ANSYS simulation are transformed into the creep compliance and compared with the theoretical calculation results of the modified Burgers model, as demonstrated in Table2. It can be found that the finite element simulation results are consistent with the theoretical calculation results, and the relative error is less than 0.70%. So the simulation results show that the Prony series results are reliable. By employing this method, the material parameters of the viscoelastic material can be accurately entered to the finite element software. However, the numerical simulation unit force is uniform, which has difference with the actual creep penetration experiments.Table 2
The comparison between equation calculation results and the ANSYS simulation.
Time (s)
Creep compliance of equation calculation (MPa−1)
Strain of ANSYS
Creep compliance of ANSYS (MPa−1)
Relative error (%)
1
1.16
E
-
02
8.07
E
-
03
1.15
E
-
02
-
0.37
101
1.87
E
-
02
1.32
E
-
02
1.88
E
-
02
0.50
201
2.37
E
-
02
1.67
E
-
02
2.38
E
-
02
0.67
301
2.72
E
-
02
1.92
E
-
02
2.74
E
-
02
0.69
601
3.41
E
-
02
2.40
E
-
02
3.43
E
-
02
0.60
901
3.86
E
-
02
2.72
E
-
02
3.88
E
-
02
0.54
1801
4.84
E
-
02
3.40
E
-
02
4.86
E
-
02
0.43
2701
5.54
E
-
02
3.89
E
-
02
5.56
E
-
02
0.36
3601
6.04
E
-
02
4.24
E
-
02
6.05
E
-
02
0.30


## 3. Penetration Experiments and Verification

### 3.1. Penetration Test Schemes of the Asphalt Mixture

To obtain the viscoelastic parameters of the practical gussasphalt mixture, the uniaxial penetration experiments, which had been practically proved effective to evaluate the shear resistance of the asphalt mixture by Sun [16], was conducted. The uniaxial penetration tests in this study are performed on rutting specimens (the gradation composition is shown in Table 3) with the dimension of 30 cm × 28 cm × 5 cm. To avoid size effection of the coarse aggregates, the dimension of the pressure head is chosen as 100 mm which is about 3 times larger than the largest coarse aggregate in the gussasphalt.Table 3
Gradation composition of the asphalt mixture.
Component
Grading/asphaltcontent
Mineral aggregate (basalt)
2# (4.75~9.5 mm)
27
3# (2.36~4.75 mm)
23
Natural sand
24
Mineral dust
26
High elastic modified bitumen
8.3%In these experiments, an instantaneous load is applied on the center of the specimen and kept unchanged, during which the time-varying displacement of the pressure head was measured, which can be utilized to obtain the timely strain of the specimen. The only measurement in the experiment is the displacement of the pressure head. The constant pressure is controlled by the material test system and the exampling frequency is 30 Hz. The whole process of the uniaxial penetration creep experiment is simple and very easy to apply, which is an obvious advantage of the uniaxial penetration method.WANCE material test system is employed in the tests, equipped with the pressure head to apply a constant pressure of 0.7 MPa for 3600 seconds at constant temperatures of 40°C, 50°C, and 60°C, respectively. To keep the specimens working in constant temperatures, a water tank controlled by a constant temperature bath is applied, shown in Figure3(a). Figure 3(a) also shows the whole configuration of the penetration experiments and Figure 3(b) presents a specimen after the experiment.Figure 3
Uniaxial penetration tests on the asphalt rutting board.
(a)
Test equipment
(b)
The asphalt rutting specimens


### 3.2. Procedure to Obtain Viscoelastic Parameters via Experimental Data and Finite Element Method

To fit the compliance curve of the gussasphalt, the stress as well as the strain should be known. In the finite element simulation mentioned above, the stress of the specimen is assumed to be uniform and known. However, in the actual experiments, the stress status of the gussasphalt is unknown. Thus, after obtaining the measured data from the experiments, to find the proper viscoelastic parameters of the gussasphalt mixtures, three steps are essential in this paper.Firstly, we assume that the vertical stress in the asphalt mixture is constant and equals the pressure applied by the head, which is 0.7 MPa. Based on the assumption, the creep compliance curve is fitted by utilizing the modified Burgers model, and each parameter in the model can be calculated by using Matlab via (12) and (13). The initial fitting results and the Prony series are listed in Table 4 by assuming Poisson’s ratio of the mixture 0.25, which are regarded as the initial parameters.Table 4
Fitting results of experimental data.
Temperature
Parameters of modified Burgers model
Parameters of Prony series (shear modulus)
E
1 (MPa)
E
2 (MPa)
η
2 (Mpa∗s)
a (Mpa∗s)
b (/s)
a
1
t
1
a
2
t
2
40
88.34
34.85
8.03
E
+
04
1.81
E
+
04
3.62
E
-
03
0.1432
1167.1852
0.6801
96.3290
50
72.02
25.95
5.94
E
+
04
1.28
E
+
04
5.47
E
-
03
0.2107
984.6694
0.6095
76.4215
60
58.92
21.22
4.53
E
+
04
9.66
E
+
03
5.47
E
-
03
0.1933
947.9533
0.6303
72.6607Secondly, using the initial material parameters in Table4, the finite element model of the uniaxial penetration creep test is established where the ANSYS element SOLID 185 is adopted, as shown in Figure 4. The circular loading area is located in the center of the model with a diameter of 10 cm and a transitional zone between the loading areas. The model boundary is 18 cm. The boundary conditions include the vertically constrained bottom surface and the bottom with four corners fixed in three directions. A constant load of 0.7 MPa is applied to the central circular loading area. The calculation results under temperature 40°C are shown in Figure 5 as an example.Figure 4
The model of uniaxial penetration creep test finite element.Figure 5
Comparison between the finite element results and the experimental results before correction (40°C).Figure5 clearly shows that the simulation results and experimental results are far from matching. This phenomenon is predictable because the inner vertical stress of the rutting plate is not the same value applied by the pressure head. By analyzing the vertical stress of the test example, Figures 6 and 7 show the contour of the vertical stress inside the plate. We can find that the vertical stress inside the asphalt material changes along the height of the specimen. If we take the average vertical stress (Z direction) of the asphalt cylinder under the pressure head as the representative value (the vertical stress points are chosen as illustrated in Figure 7 and the stress points under the pressure head line along the dimension due to the symmetry), it is significantly smaller than the value applied by the pressure head in the tests. Meanwhile, although the mean value of the vertical principal stress varies with time, as shown in Figure 8, the changing amplitude is tiny and less than 0.1% of the vertical stress during the test time. Thus, the average vertical stress of the gussasphalt mixture under the pressure head in the end time of the experiment can be approximately taken as the mean vertical stress of the specimen along the height, which is 0.5929 Mpa in this example. It is much smaller and more precise than the initial pressure 0.7 MPa to calculate the creep compliance curve. Here a new parameter α is introduced to update the model parameters of the modified Burgers model.(14)α=p2p1,where p1 is the pressure stress before updating and p2 is the updated stress. If the parameter α equals 1 or α-1<ε (where ε is a small value) then we can conclude that the results have converged. Using the definition, the correction coefficient α in current example can be calculated (where α is 0.7/0.5929 = 1.18). Then by turning the creep compliance J(t) into αJ(t) and refit the curve, the updated material parameters can be achieved, as shown in Table 5.Table 5
Experimental data fitting results after stress correction.
Temperature
Parameters of modified Burgers model
Parameters of Prony series(shear modulus)
Stress correctionα
Iteration number
E
1 (MPa)
E
2 (MPa)
η
2 (Mpa∗s)
a (Mpa∗s)
b (/s)
a
1
t
1
a
2
t
2
40
74.83
29.52
6.80
E
+
04
1.53
E
+
04
3.62
E
-
03
0.1432
1167.1852
0.6801
96.3290
1.18
1
50
61.00
21.98
5.03
E
+
04
1.08
E
+
04
5.47
E
-
03
0.2107
984.6694
0.6095
76.4215
1.18
1
60
49.91
17.97
3.84
E
+
04
8.18
E
+
03
5.47
E
-
03
0.1933
947.9533
0.6303
72.6607
1.18
1Figure 6
Vertical stress cloud of central section.Figure 7
Distribution of vertical stress under the loading area.Figure 8
Time history of mean values of the vertical stress.Thirdly, theupdated material parameters are utilized to resimulate the experiment. The results before and after correction as well as the relative errors are shown in Figures 9 and 10. The same procedure can be applied to the experimental specimens under temperature 50°C and 60°C. Figure 9 demonstrates the displacements measured in the penetration experiments as well as predicted results before and after stress correction by finite element analysis under different temperatures. Figure 10 shows the comparison of the relative errors between the simulated and experimental results before and after stress correction. From the comparison of these figures, it can be clearly seen that the accuracy of finite element simulation of the penetration tests at different temperatures is greatly improved after the introduction of the correction coefficient. The maximum error of the finite element simulation is about 5% comparing with the experimental results. More iterations by following the procedure in Figure 11 can still be implemented to obtain more precise results if necessary. Then the rheological parameters of the gussasphalt mixtures can be obtained successfully.Figure 9
Comparison between the finite element results and the experimental results before and after correction at 40°C (a), 50°C (b), or 60°C (c).
(a)
(b)
(c)Figure 10
Comparison of relative error before and after correction at 40°C (a), 50°C (b), or 60°C (c).
(a)
(b)
(c)Figure 11
Flow chart of stress-correction method.


## 3.1. Penetration Test Schemes of the Asphalt Mixture

To obtain the viscoelastic parameters of the practical gussasphalt mixture, the uniaxial penetration experiments, which had been practically proved effective to evaluate the shear resistance of the asphalt mixture by Sun [16], was conducted. The uniaxial penetration tests in this study are performed on rutting specimens (the gradation composition is shown in Table 3) with the dimension of 30 cm × 28 cm × 5 cm. To avoid size effection of the coarse aggregates, the dimension of the pressure head is chosen as 100 mm which is about 3 times larger than the largest coarse aggregate in the gussasphalt.Table 3
Gradation composition of the asphalt mixture.
Component
Grading/asphaltcontent
Mineral aggregate (basalt)
2# (4.75~9.5 mm)
27
3# (2.36~4.75 mm)
23
Natural sand
24
Mineral dust
26
High elastic modified bitumen
8.3%In these experiments, an instantaneous load is applied on the center of the specimen and kept unchanged, during which the time-varying displacement of the pressure head was measured, which can be utilized to obtain the timely strain of the specimen. The only measurement in the experiment is the displacement of the pressure head. The constant pressure is controlled by the material test system and the exampling frequency is 30 Hz. The whole process of the uniaxial penetration creep experiment is simple and very easy to apply, which is an obvious advantage of the uniaxial penetration method.WANCE material test system is employed in the tests, equipped with the pressure head to apply a constant pressure of 0.7 MPa for 3600 seconds at constant temperatures of 40°C, 50°C, and 60°C, respectively. To keep the specimens working in constant temperatures, a water tank controlled by a constant temperature bath is applied, shown in Figure3(a). Figure 3(a) also shows the whole configuration of the penetration experiments and Figure 3(b) presents a specimen after the experiment.Figure 3
Uniaxial penetration tests on the asphalt rutting board.
(a)
Test equipment
(b)
The asphalt rutting specimens


## 3.2. Procedure to Obtain Viscoelastic Parameters via Experimental Data and Finite Element Method

To fit the compliance curve of the gussasphalt, the stress as well as the strain should be known. In the finite element simulation mentioned above, the stress of the specimen is assumed to be uniform and known. However, in the actual experiments, the stress status of the gussasphalt is unknown. Thus, after obtaining the measured data from the experiments, to find the proper viscoelastic parameters of the gussasphalt mixtures, three steps are essential in this paper.Firstly, we assume that the vertical stress in the asphalt mixture is constant and equals the pressure applied by the head, which is 0.7 MPa. Based on the assumption, the creep compliance curve is fitted by utilizing the modified Burgers model, and each parameter in the model can be calculated by using Matlab via (12) and (13). The initial fitting results and the Prony series are listed in Table 4 by assuming Poisson’s ratio of the mixture 0.25, which are regarded as the initial parameters.Table 4
Fitting results of experimental data.
Temperature
Parameters of modified Burgers model
Parameters of Prony series (shear modulus)
E
1 (MPa)
E
2 (MPa)
η
2 (Mpa∗s)
a (Mpa∗s)
b (/s)
a
1
t
1
a
2
t
2
40
88.34
34.85
8.03
E
+
04
1.81
E
+
04
3.62
E
-
03
0.1432
1167.1852
0.6801
96.3290
50
72.02
25.95
5.94
E
+
04
1.28
E
+
04
5.47
E
-
03
0.2107
984.6694
0.6095
76.4215
60
58.92
21.22
4.53
E
+
04
9.66
E
+
03
5.47
E
-
03
0.1933
947.9533
0.6303
72.6607Secondly, using the initial material parameters in Table4, the finite element model of the uniaxial penetration creep test is established where the ANSYS element SOLID 185 is adopted, as shown in Figure 4. The circular loading area is located in the center of the model with a diameter of 10 cm and a transitional zone between the loading areas. The model boundary is 18 cm. The boundary conditions include the vertically constrained bottom surface and the bottom with four corners fixed in three directions. A constant load of 0.7 MPa is applied to the central circular loading area. The calculation results under temperature 40°C are shown in Figure 5 as an example.Figure 4
The model of uniaxial penetration creep test finite element.Figure 5
Comparison between the finite element results and the experimental results before correction (40°C).Figure5 clearly shows that the simulation results and experimental results are far from matching. This phenomenon is predictable because the inner vertical stress of the rutting plate is not the same value applied by the pressure head. By analyzing the vertical stress of the test example, Figures 6 and 7 show the contour of the vertical stress inside the plate. We can find that the vertical stress inside the asphalt material changes along the height of the specimen. If we take the average vertical stress (Z direction) of the asphalt cylinder under the pressure head as the representative value (the vertical stress points are chosen as illustrated in Figure 7 and the stress points under the pressure head line along the dimension due to the symmetry), it is significantly smaller than the value applied by the pressure head in the tests. Meanwhile, although the mean value of the vertical principal stress varies with time, as shown in Figure 8, the changing amplitude is tiny and less than 0.1% of the vertical stress during the test time. Thus, the average vertical stress of the gussasphalt mixture under the pressure head in the end time of the experiment can be approximately taken as the mean vertical stress of the specimen along the height, which is 0.5929 Mpa in this example. It is much smaller and more precise than the initial pressure 0.7 MPa to calculate the creep compliance curve. Here a new parameter α is introduced to update the model parameters of the modified Burgers model.(14)α=p2p1,where p1 is the pressure stress before updating and p2 is the updated stress. If the parameter α equals 1 or α-1<ε (where ε is a small value) then we can conclude that the results have converged. Using the definition, the correction coefficient α in current example can be calculated (where α is 0.7/0.5929 = 1.18). Then by turning the creep compliance J(t) into αJ(t) and refit the curve, the updated material parameters can be achieved, as shown in Table 5.Table 5
Experimental data fitting results after stress correction.
Temperature
Parameters of modified Burgers model
Parameters of Prony series(shear modulus)
Stress correctionα
Iteration number
E
1 (MPa)
E
2 (MPa)
η
2 (Mpa∗s)
a (Mpa∗s)
b (/s)
a
1
t
1
a
2
t
2
40
74.83
29.52
6.80
E
+
04
1.53
E
+
04
3.62
E
-
03
0.1432
1167.1852
0.6801
96.3290
1.18
1
50
61.00
21.98
5.03
E
+
04
1.08
E
+
04
5.47
E
-
03
0.2107
984.6694
0.6095
76.4215
1.18
1
60
49.91
17.97
3.84
E
+
04
8.18
E
+
03
5.47
E
-
03
0.1933
947.9533
0.6303
72.6607
1.18
1Figure 6
Vertical stress cloud of central section.Figure 7
Distribution of vertical stress under the loading area.Figure 8
Time history of mean values of the vertical stress.Thirdly, theupdated material parameters are utilized to resimulate the experiment. The results before and after correction as well as the relative errors are shown in Figures 9 and 10. The same procedure can be applied to the experimental specimens under temperature 50°C and 60°C. Figure 9 demonstrates the displacements measured in the penetration experiments as well as predicted results before and after stress correction by finite element analysis under different temperatures. Figure 10 shows the comparison of the relative errors between the simulated and experimental results before and after stress correction. From the comparison of these figures, it can be clearly seen that the accuracy of finite element simulation of the penetration tests at different temperatures is greatly improved after the introduction of the correction coefficient. The maximum error of the finite element simulation is about 5% comparing with the experimental results. More iterations by following the procedure in Figure 11 can still be implemented to obtain more precise results if necessary. Then the rheological parameters of the gussasphalt mixtures can be obtained successfully.Figure 9
Comparison between the finite element results and the experimental results before and after correction at 40°C (a), 50°C (b), or 60°C (c).
(a)
(b)
(c)Figure 10
Comparison of relative error before and after correction at 40°C (a), 50°C (b), or 60°C (c).
(a)
(b)
(c)Figure 11
Flow chart of stress-correction method.


## 4. Conclusion

In this paper, to obtain the viscoelastic parameters of gussasphalt mixtures, the explicit parameters of modified Burgers model based on Prony series are derived and the uniaxial penetration creep experiments of gussasphalt mixture are carried out. The parameter identification procedure by stress correcting is proposed to find accurate model parameters of gussasphalt mixture and the effectiveness of the method is validated by FEM simulation and experimental data.Based on the work, the following conclusions can be drawn:(1)
The explicit form Prony series of the modified Burgers model is obtained. The viscoelastic properties of gussasphalt mixtures can be easily input into the commercial finite element software to explore the nonlinear properties. The correctness of the transformation equations of the modified Burgers model is verified by the finite element model.(2)
The viscoelastic parameter identification procedure combined with the stress-correcting method of the gussasphalt mixture based on uniaxial penetration creep experiment under different temperatures is carried out. After stress correcting, the experimental data can be well simulated by FEM software, following the presented procedure.(3)
The penetration experiment can imitate the stress status of the pavement and is relatively easy to realize. The method combined with the Prony series via the modified Burgers model as well as stress-correcting iteration proposed in this paper can be extended to other materials with viscoelastic properties.



---

*Source: 4324765-2017-03-07.xml*