Next Article in Journal
Pentavalent Vanadium Species as Potential Corrosion Inhibitors of Al2Cu Intermetallic Phase in the Sulfuric(VI) Acid Solutions
Next Article in Special Issue
Double-Recycled Reclaimed Asphalt Pavement: A Laboratory Investigation at Low Temperatures Based on Different Mathematical Approaches
Previous Article in Journal
Palladium/Carbon Nanofibers by Combining Atomic Layer Deposition and Electrospinning for Organic Pollutant Degradation
Previous Article in Special Issue
Characterization of Bonding between Asphalt Concrete Layer under Water and Salt Erosion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

From Complex Modulus E* to Creep Compliance D(t): Experimental and Modeling Study

1
École de Technologie Supérieure (ÉTS), Construction Engineering, 1100 Notre-Dame Ouest, Montreal, QC H3C 1K3, Canada
2
Institut de Recherche en Constructibilité, Université Paris-Est, 94234 École Spéciale des Travaux Publics, 28 avenue du Président Wilson, 94234 Cachan CEDEX, France
*
Author to whom correspondence should be addressed.
Materials 2020, 13(8), 1945; https://doi.org/10.3390/ma13081945
Submission received: 11 January 2020 / Revised: 2 April 2020 / Accepted: 9 April 2020 / Published: 21 April 2020

Abstract

:
Creep compliance (D(t)) is a very important input for the thermal cracking resistance in the Mechanistic-Empirical Pavement Design Guide (MEPDG). The aim of the work presented here is to predict the results of creep compliance D(t) from the result of complex modulus E*(ω). The work plan is divided in two main parts: an experimental part consisting of creep tests, and a modeling part. Three configurations were compared together, namely direct tensile, direct compression and indirect tensile tests. The modelling part consists of using a 2S2P1D model coupled to Kopelman approximation to switch from the frequency domain to the time domain. Additionally, 2S2P1D was used to calibrate the generalized Kelvin–Voigt model and get the creep compliance directly from E* results. The experimental results show that D(t) from direct tensile and direct compression are the same in the viscoelastic domain and are greater than D(t) from the indirect tensile test. The indirect tensile test (IDT) seems to be very difficult to achieve compared to the other two variants. The converted results using the 2S2P1D model coupled to Kopelman approximation and the results from the GKV model describe the experimental points very well.

Highlights

  • The use of rheological models is a very simple and good method to predict creep compliance D(t) from complex modulus E(ω).
  • Creep compliance from direct tensile and direct compression tests are the same and are more representative to the field.
  • The IDT creep test configuration allow to measure parallel creep compliance which is not representative with thermal cracks observed in the real roads.
  • Direct creep compliance tests are easier to perform comparing to the IDT configuration proposed by AASHTO T322-07.

1. Introduction

The thermomechanical behavior of bituminous asphalt mixes is very complex. According to Benedetto [1] four behaviors are observed depending on the strain amplitude (|ε|) and the number of cycles (N) of the applied loading. For small strain amplitudes (|ε| < 10−4) and a few hundred loading cycles, the behavior can be considered as linear viscoelastic (LVE). For tens of thousands of cycles, the phenomenon of fatigue manifests. Rutting is created by the accumulation of irreversible permanent deformations caused by near-breaking amplitudes. The variation of the temperature can cause deformations of a few percent (|ε| >> 10−4), the observed behavior is strongly nonlinear.
The viscoelastic linear behavior of asphalt mixes is mainly characterized by two types of tests: quasi-static loading tests in the time domain (relaxation and creep) and cyclic loading tests in the frequency domain (complex modulus). Relaxation is the evolution of the stress in time for the application of a constant deformation, and creep is the evolution of the strain for a constant stress applied. As for the complex modulus, it is the expression of rigidity of viscoelastic materials under cyclic loading; it represents the ratio between the cyclic stress by the cyclic strain.
Laplace–Carson transform is a mathematical tool allowing substituting variables, in particular to go from the time domain to the frequency domain. Because integral equations describing material LVE behavior in a time domain are difficult to apply, the Laplace–Carson transform is used to turn integral equations into algebraic, therefore simplifying calculations (Equation (1)) [1,2].
g * ( p ) = p 0 + e p t g ( t ) d t
where g(t) is a time dependent function, g*(p) is function dependent on p, p is a complex variable, corresponding to time in the transform domain. If p = iω (ω is the angular frequency) we have:
E * ( ω ) × D *   ( ω ) = 1
where E* is the complex modulus (relaxation in frequency domain) and D* is the complex compliance (creep compliance in frequency domain).
Creep compliance (ratio between deformations in time by the maximum stress applied) is a very important input for low temperature or thermal cracking module in the Mechanistic-Empirical Pavement Design Guide (MEPDG) software [2,3]. Too low creep may lead to cracking at low temperature and too high creep may lead to high potential of rutting at service temperatures [4]. This type of test (creep) is also a very important element for characterization of bitumen at low temperature; the PG High-Low temperature system is based on a creep test (bending beam rheometer, BBR) for the determination of the low temperature of bitumen [5,6,7,8,9].
The American Association of State Highway and Transportation Officials (AASHTO T322-07) describes a test protocol for the creep compliance [10]. The standard proposes to work in indirect tensile test (IDT) configuration on a cylindrical test specimen of 150 mm in diameter and 38 mm to 50 mm of thickness. The specimen compaction is performed by the gyratory compactor. The test is generally carried out at three different temperatures ranging from 10 °C to −30 °C depending on the grade of the bitumen. A static load is applied along a diametral axis in order to have an indirect tensile stress. The creep compliance is a non-destructive test; the applied loading must always be in the linear viscoelastic domain. However, the load must be high enough to cause sufficient horizontal deformation so the noise in the data acquisition process is insignificant. During the test, vertical and horizontal deformations are measured and the creep compliance represents the time-dependent horizontal deformation by the applied stress, which can be calculated by the given Equation (3):
D ( t ) = Δ X t m , t × D a v g × b a v g P a v g × G L × C c m p l
C c m p l = 0.6354 × ( X Y ) 1 0.332
where:
D(t):creep compliance at time t (1/kPa).
GL:gauge length in meters (mm).
Davg:average diameter of all specimens (mm).
bavg:average thickness of all specimens (mm).
Pavg:average creep load (kN).
ΔXtm,t:trimmed mean of the normalized horizontal deformations of all specimen faces at time t (mm/mm).
Ccmpl:correction factor.
X Y :absolute value of the ratio of the normalized, trimmed mean of the horizontal deformations (i.e., ΔXtm,t) to the normalized, trimmed mean of the vertical deformations (i.e, ΔYtm,t) at a time corresponding to ½ the total creep test time for all specimen faces.
The National Cooperative Highway Research Program (NCHRP) in its report 530 presents a laboratory test program to evaluate and compare measurement of creep compliance using IDT configuration, uniaxial compression, and uniaxial tension tests [11,12]. The objective was to review the test procedure described in AASHTO T322-07. The results show that the behavior is very anisotropic. The compliance measured with the IDT test was lower than the compliance determined in uniaxial compression, which in turn was lower than the uniaxial tension compliance. The NCHRP chose to validate the test procedure described in AASHTO T322-07 because the compliance measured in perpendicular plane in relation to the compaction is generally higher than that determined parallel to compaction. However, we note that thermal cracks at low temperatures are generally perpendicular to the direction of compaction of the pavement. Lowering road temperature creates contraction stress; the contraction is blocked in the longitudinal direction of the road which leads to transverse cracking [4,13]. The IDT configuration allows determining parallel creep compliance contrary to the direct tensile and direct compression configurations. Additionally, we note that compaction procedure with gyratory compactor is not very representative of field compaction.
The viscoelastic theory offers the possibility to get creep compliance results directly from complex modulus results. Olard et al. conducted a study at binder scale for the prediction of BBR results (creep compliance test) from shear complex modulus results (G*) [14]. Complex modulus tests were performed by Métravib (tensile/compression configuration) then the results were converted to shear complex modulus using Poisson’s ratio (|E*| = 2 (1 + ν) |G*|). incompressibility was supposed, so that ν = 0.5 and the approximation |E*| = 3 | G*|was used [4,14]. To change from the frequency domain to the time domain, Olard et al. used four approximate equations namely Kopelman equation (1958), Christensen equation (1982), Ninomiya and Ferry equation (1959) and Schwarzl and Struik equation (1968) [14,15,16]. Creep compliance results obtained from BBR data were compared with the converted results from the approximate equations (prediction). The predicted results describe the experimental data well.
Any linear viscoelastic body can be modelled by the generalized Kelvin–Voigt model (GKV) (Figure 1a) [17]. This model is an assembly of simple rheological elements like springs and dashpots, the complex modulus and the creep compliance are given by Equations (5) and (6).
E * ( ω ) = ( 1 E 00 + 1 i ω η 00 + i = 1 n 1 E i + i ω η i ) 1
D ( t ) = 1 E 00 + i = 1 n 1 1 E i   ( 1 e t τ i ) + 1 η 00
where:
E*(ω):complex modulus (MPa).
E00 and η00value of the complex modulus and the viscosity when ωτ → 0.
ω:is the angular frequency, so that ω/2π is equal to the frequency.
ηi and Ei:viscosity and elastic modulus of each Kelvin Voigt branch.
D(t):creep compliance (1/MPa).
τ:is the relaxation time for single branch.
The unidimensional rheological model with continuous spectrum called 2S2P1D (Figure 1b) [18] is widely used for the prediction of the LVE behavior of bituminous materials; it consists of two elastic (spring) elements, two parabolic creep elements and one viscous (dashpot) element. The complex modulus is given by the Equation (7). However, no analytical expression of the creep compliance function of this model is available in the time domain.
E * ( ω ) = E 00 + E 0 E 00 1 + δ ( i ω t ) k + ( i ω t ) h + ( i ω β t ) 1
where:
  • E0 is the value of the complex modulus when ωτ →∞.
  • δ, k and h are dimensionless constants and β is a dimensionless parameter, related to Newtonian viscosity η.
Several studies present a method to get viscosities (ηi) and modulus (Ei) of the generalized Kelvin–Voigt (GKV) model from 2S2P1D model [19,20,21]. After the calibration of the 2S2P1D model in the frequency domain (E*), the method consists of minimizing the difference between the GKV model results and the results given by a reference viscoelastic model (2S2P1D).

2. Methods and Materials

2.1. Methods

The objective of this study is to predict the results of creep compliance in time domain D(t) directly from the results of complex modulus.
The methodology is divided into two main parts: an experimental part and a modelling part. The experimental part consists of realization of complex modulus tests and creep compliance tests (Figure 2).
Complex modulus tests were performed in traction-compression configuration according to Quebec’s standard LC26-700. Seven test temperatures ranging from −20 °C to 40 °C with a step of 10 °C and a frequency sweep from 0.1 Hz to 10 Hz were used. The imposed strain was ±50 µstrain. Complex modulus results were fitted using 2S2P1D model, then two different methods were used to switch from the frequency domain to the time domain for creep compliance. First, 2S2P1D results were coupled with the approximation proposed by Kopelman which is mathematically simpler than the other three equations proposed in Olard’s et al. [14] study. Second, the GKV model was calibrated according to the minimization with 2S2P1D model proposed by Tiouajni et al. [19]. and the creep compliance in time domain is given by Equation (6).
Creep tests were performed using three configurations namely: direct tensile configuration, direct compression configuration and indirect tensile configuration. Direct tensile and direct compression tests were achieved on cylindrical specimens with a diameter of 75 ± 1 mm and a height of 135 ± 2 mm. A static traction or compression stress load was applied on the samples for 300 s; it took 0.1 s to reach the imposed load, without impact, from a contact force of 50 N. The measurement of deformation was insured by three extensometers of 25 mm in length. Four test temperatures were used ranging from −20 °C to 10 °C with a step of 10 °C. The amplitude of the imposed stress was fixed by taking in to account complex modulus results. For each test temperature, before to start the creep test and at the end of it, a partial complex modulus test was done in order to check if the material was still in the viscoelastic behavior domain. Indirect tensile creep tests were performed on cylindrical samples with a diameter of 150 mm and a thickness of 38 mm. AASHTO T322-07 has not set a load application speed or a maximum amplitude, so, in order to have the same test parameters, the IDT loading was applied in 0.1 s and the amplitude was fixed manually in order to have a minimum 50 µstrain at 0.1 s.

2.2. Materials

A 0–10 mm surface mix, namely ESG-10, was used. The gradation is given in Figure 3. The bitumen is an unmodified PG 58-34, the penetration and the ring and ball temperature of this bitumen are respectively 153 ± 1, 1/10 mm and 43 ± 0.1 °C. The mixing temperature used was 150 °C, and the compaction temperature 135 °C. The aggregates were heated 15 °C more than the bitumen as indicated in the LC (26-003) test method, i.e., 165 °C.
The cylindrical samples for complex modulus tests, direct tensile and direct compression creep tests were cored from slabs (100 mm × 180 mm × 500 mm) compacted in laboratory using MLPC (Materiel Laboratoire Ponts et Chaussées) wheel compactor according to Quebec’s LC test method (Figure 4). For indirect tensile creep tests, the specimens were cut from cylindrical samples compacted by gyratory compactor (Figure 5). Before doing tests, the air voids content was measured using hydrostatic weighing and the maximum density of asphalt mix (Dmm).

3. Results and Discussion

3.1. Complex Modulus Results and 2S2P1D Calibration

Complex modulus results are presented in Figure 6 and Figure 7. The calibration of the 2S2S1D model was performed in Cole–Cole curve and black space in order to fix the shape parameters (k, h, δ and β). Asymptotic moduli E00 and E0 were fixed by a linear extrapolation respectively at very low frequencies and very high frequencies. Characteristic time τ depends on temperature as follows:
a T = τ τ 0
where τ0 is the characteristic time at the reference temperature (Tref) and aT is the shift factor. aT values may be approximated by the Williams–Landel–Ferry (WLF) equation (Equation (9)) when Time-Temperature superposition is valid [22,23,24].
log ( a T ) = C 1   ( T T r e f ) C 2 + ( T T r e f )
where C1 (23.90) and C2 (148.90) are empirical constants depending on the material. The calibration of the GKV model was performed according to the minimization proposed by [19,20], 19 branches of Kelvin–Voigt model were used.
Master curve is presented in Figure 7 at Tref = 0 °C. As it is noted, the 2S2P1D model and GKV model describe very well the experimental points. The error between the experimental points and 2S2P1D model are presented in Figure 7. The majority of the points are located within an error interval of ± 5% for the norm of complex modulus (|E*|) and ± 2° for phase angle (φ). The maximum error noted is 10.5% for |E*| and 3.3° for φ. The maximum error between the GKV model and the 2S2P1D model is 14% for |E*| and 3° for φ.
The calculation of creep compliance is done using two different methods. First, using the GKV model, the creep compliance of this model is given by Equation (6). The second method consists of coupling the 2S2P1D results with the approximation proposed by Kopelman. The 2S2P1D model doesn’t have a creep function in frequency domain while using Kopelman Equation (15) can allow to switch from the frequency domain to the time domain. The explanations of the calculations are given below:
E * = E 1 + i E 2
D * = D 1 i D 2
E * ( ω ) × D * ( ω ) = 1
D 1 = E 1 E 1 2 + E 2 2
D 2 = E 2 E 1 2 + E 2 2
D ( t ) = | D ( ω ) | ω = 1 / t
where:
  • E1: the real part of the complex modulus (storage modulus) (MPa).
  • E2: the imaginary part of the complex modulus (loss modulus) (MPa).
  • D1: the real part of the complex compliance (1/MPa).
  • D2: the imaginary part of complex compliance (1/MPa).
  • i: the imaginary unit.
Creep compliance from the GKV model and from the 2S2P1D model coupled with Kopelman approximation are presented in Figure 8. The two curves are very close. The error is calculated from 10−5 s to 10+5 s using Equation (16).
%   E r r o r = D ( t ) 2 S 2 P 1 D + K o p e l m a n D ( t ) G K V D ( t ) 2 S 2 P 1 D + K o p e l m a n
The results are presented in Figure 9. The error is between −2% to −15%.

3.2. Direct Tensile and Direct Compression Creep Tests Results

Creep tests were performed using three configurations namely direct tensile test, direct compression test and indirect tensile test. As shown in Figure 10, for direct tensile tests and direct compression tests, a static stress was applied to the cylindrical samples. This configuration allows to have a homogeneous stress distribution, so the creep compliance can be calculated directly using the Equation (17).
The load was applied in 0.1 s. The amplitude of force and stress applied are presented in Table 1. For the first 10 to 15 s, the applied load was not stable. Creep compliance was calculated (Equation (17)) when the load is stabilized at ± 2% of the target value at t equal to 0.1, 1, 3, 10, 30, 100 and 300 s.
D ( t ) = ε ( t ) σ m a x
The results of creep compliance are presented in Figure 11 at different test temperatures (isotherms at −20 °C, −10 °C, 0 °C and 10 °C). According to the results, we notice that the creep compliances from direct tensile and from direct compression are the same. The majority of points are superposed on the equality line with an error interval between ±3% and a maximum error of 18% at high creep compliance value.
This results are contradictory with the results presented by NCHRP [10]. The explanation could be that the results presented by NCHRP were out of the VEL domain. As mentioned before, to verify that our results are inside the linear viscoelastic domain, partial complex modulus tests were performed before and after direct traction and compression tests. Exceeding the viscoelastic limit is characterized by a loss of |E*|. Table 2 presents the percentage of loss of |E*| at test temperatures and frequencies. These percentages were calculated using Equation (18). The non-significant losses observed confirm that viscoelastic limit was not exceeded.
%   o f   l o s s = | E b e f o r e   d i r e c t   c r e e p   t e s t s * | | E a f t e r   d i r e c t   c r e e p   t e s t s * | | E b e f o r e   d i r e c t   c r e e p   t e s t * |
The application of the time temperature superposition principle (TTSP) allows having a unique curve. The WLF equation is widely used to determine the shift factor (at) for bituminous materials. However, this shift factor is used in the frequency domain. In order to apply the TTSP to creep compliance results, two methods are used to determine the shift factor in time domain. First a graphical method is used (experimental: EXP) and second the inverse of at WLF calculated from complex modulus test were used. The results are presented in Figure 12. It can be noted that the graphical exp shift factor can be approximated by the inverse of WLF equation.
The master curve built with the experimental points is presented at 0 °C (Figure 13) comparatively with the creep compliance predicted with GKV model and 2S2P1D model coupled to Kopelman approximation. Both of models seem to describe very well the experimental points. These results confirm that it is possible to predict creep compliance (static load test) results directly from complex modulus (cyclic load test) results.

3.3. Indirect Tensile Creep Test Results

IDT creep tests were carried out by applying a static load on a diametral axis on cylindrical specimens. The deformations were measured on both sides of the specimens. The applied loads as well as the deformations on both faces are presented in Figure 14.
The difference in strain measured on both sides increased by increasing test temperature. At −20 °C the variation of horizontal strain between face 1 and face 2 was 2% ± 5% and for vertical strain the difference was 30% ± 10%. At −10 °C the gap for horizontal strain was 4% ± 3% and vertical strain it was 5% ± 1%. Finally, for 0 °C the difference was very significant for horizontal strain 37% ± 3% and 16% ± 1% for vertical strain. Additionally, it can be noted that the maximal vertical strain at 0 °C was greater than the limit of viscoelastic domain found in the literature. However, because of the shape of the specimen, it was not possible to test the complex modulus before and after the static load. The gap between strains measured on both faces can be an indicator that the viscoelastic limit was exceeded; for that reason, the test was stopped after 100 s and not carried out at 10 °C.
The ratio between the mean horizontal strain (X) and the mean vertical strain (Y) was measured, the results are presented in Figure 15. The points follow a polynomial law of order two in the range of temperature tested.
IDT creep results were calculated using the equation proposed by AASHTO. The results (isotherms) are presented in Figure 16. The values appear lower than direct tensile (DT) and direct compression (DC) tests except for the test points at 0 °C, which may confirm that the linear viscoelastic limit was exceeded. At −20 °C, creep compliance from IDT test was lower than direct tensile and direct compression creeps with factor ranging between 15.0 to 1.05. At −10 °C there is a factor of 3.0 to 1.05 with two points at 0.85. At 0 °C, the factor ranges from 0.75 to 0.35.
The master curves of IDT creep compliance test are presented at 0 °C in Figure 17 where the creep compliance from GKV model and the creep compliance from the 2S2P1D model coupled to Kopelman approximation and from direct tensile and compression creep test are compared. To build the master curve for IDT, test shift factors were determined graphically and seem to be ten times bigger the inverse of WLF law using the same C1 and C2 from complex modulus test. It is difficult to know the reason of that, C1 and C2 are empirical constants depending on the material. However, this may be related to the correction factor proposed by the AASHTO.
The results show that creep compliance determined from IDT (corrected) test appears lower than the creep predicted from GKV model and 2S2P1D model coupled to Kopelman approximation except for the points determined at 0 °C. This could be explained by the fact that the viscoelastic domain was exceeded and therefore the creep was much more pronounced in this case. Additionally, Figure 17 shows the non-corrected IDT creep compliance (without applying Ccmpl). In this case, to build the master curve, the shift factors were determined using the inverse of WLF law.

4. Conclusions

The objective of the work presented here was to predict creep compliance D(t) results from complex modulus E*(ω) results. Several methods are proposed in the literature allowing to go from the frequency domain to the time domain and vice-versa such as the inverse of the transform of Carson-Laplace and by more or less complex numerical methods. Rheological models, if properly calibrated, seem to be a very good and simple method to go from E* (cyclic test) to D(t) (static test). To achieve this objective, the proposed working method consists of two main parts: an experimental part and a modelling part. For the experimental part, two types of test were performed namely: cyclic test–E*(ω) and creep static test–D(t). Complex modulus tests were done using direct tensile and compression tests. For creep test, a comparison was done between indirect tensile configuration (proposed by AASHTO) and direct tensile and direct compression tests.
The results show that D(t) from direct tensile and direct compression tests are the same in the viscoelastic behavior domain and are greater than D(t) from indirect tensile test. This can be explained by the direction of the measured creep modulus. D(t) perpendicular to the direction of compaction are higher compared to the parallel one. Compared to the field, the cracks observed by lowering of temperature are in an axis perpendicular to the direction of compaction of the roads. Moreover, the compaction method proposed by AASHTO (with the gyratory compactor) is not very representative compared to compaction on site. For direct tensile and direct compression tests, the samples were compacted using MLPC wheel compactor. Additionally, it is complicated to compare creep test in IDT configuration to direct compression or tensile tests. For IDT, the results show significant difference between the strain measured in both faces. Further, it’s very difficult to know if the load applied remains inside the viscoelastic behavior domain for IDT configuration. The proposed method for converting E*(ω) to D(t) consists in using existing rheological models. 2S2P1D is very commonly used for modelling of asphalt materials but this model does not have a creep function in time domain. To change from the frequency domain to the time domain, we propose to couple the results from 2S2P1D to the approximation proposed by Kopelman. The generalized Kelvin–Voigt model can predict the creep compliance in time domain, but it is more difficult to calibrate. A minimization method is proposed in literature to calibrate GKV model from 2S2P1D model and therefore having D(t) results directly from E* results. The two methods were compared, and the results show a very good concordance between them.

Author Contributions

Conceptualization, methodology and validation, A.D. (Abdeldjalil Daoudi), D.P., A.D. (Anne Dony) and A.C.; writing—original draft preparation, A.D. (Abdeldjalil Daoudi); writing—review and editing, D.P., A.D. (Anne Dony) and A.C. All authors have read and agree to the published version of the manuscript.

Funding

This research was partially funded by Transport Quebec.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Di Benedetto, H.; Corté, J.-F. Matériaux routiers bitumineux: Constitution et propriétés thermomécaniques des mélanges. Hermès Sci. Publ. 2004, 2, 75–215. [Google Scholar]
  2. Huang, Y.H. Pavement Analysis and Design, 2nd ed.; Pearson: Upper Saddle River, NJ, USA, 1993; Appendix A; pp. 655–670. [Google Scholar]
  3. Puppala, A.J. Estimating Stiffness of Subgrade and Unbound Materials for Pavement Design; Transportation Research Board: Washington, DC, USA, 2008; Volume 382. [Google Scholar]
  4. Olard, F. Comportement Thermomécanique des Enrobés Bitumineux à Basses Températures. Relat. Entre Les Propriétés Du Liant Et De L’enrobé. Ph.D. Thesis, The Institut National des Sciences Appliquées de Lyon, Lyon, France, 2003. [Google Scholar]
  5. Olidis, C.; Hein, D. Guide for the Mechanistic-Empirical Design of New and Rehabilitated Pavement Structures Materials Characterization: Is Your Agency Ready. In Proceedings of the 2004 Annual Conference of the Transportation Association of Canada, Québec City, QC, Canada, 19–22 September 2004. [Google Scholar]
  6. Al-Qadi, I.L.; Loulizi, A.; Aref, S.; Masson, J.F.; McGhee, K.M. Modification of bending beam rheometer specimen for low-temperature evaluation of bituminous crack sealants. Transp. Res. Rec. 2005, 1933, 96–106. [Google Scholar] [CrossRef]
  7. Anderson, D.A.; Kennedy, T.W. Development of SHRP binder specification. J. Assoc. Asph. Paving Technol. 1993, 62, 481–501. [Google Scholar]
  8. Bahia, H.U.; Anderson, D.A. The development of the bending beam rheometer; basics and critical evaluation of the rheometer. Phys. Prop. Asph. Cem. Binders. Astm Int. 1995, 28–50. [Google Scholar] [CrossRef]
  9. Bahia, H.U.; Anderson, D.A.; Christensen, D.W. The bending beam rheometer; a simple device for measuring low-temperature rheology of asphalt binders (with discussion). J. Assoc. Asph. Paving Technol. 1992, 61, 117–153. [Google Scholar]
  10. T322-07–Standard Method of Test for Determining the Creep Compliance and Strength of Hot-Mix Asphalt (HMA) Using the Indirect Tensile Test Device; American Association of State Highway and Transportation Officials: Washington, DC, USA, 2007.
  11. Christensen, D.W.; Bonaquist, R.F. Evaluation of Indirect Tensile Test (IDT) Procedures for Low-Temperature Performance of Hot Mix Asphalt; Transportation Research Board: Washington, DC, USA, 2004; Volume 530. [Google Scholar]
  12. Richardson, D.N.; Lusher, S.M. Determination of Creep Compliance and Tensile Strength of Hot-Mix Asphalt for Wearing Courses in Missouri. Missouri. Dept. Transp. 2008, RI05-052, 75. [Google Scholar]
  13. Olard, F.; Di Benedetto, H.; Eckmann, B.; Triquigneaux, J.P. Linear viscoelastic properties of bituminous binders and mixtures at low and intermediate temperatures. Road Mater. Pavement Des. 2003, 4, 77–107. [Google Scholar] [CrossRef]
  14. Olard, F.; Di Benedetto, H.; Eckmann, B. Rhéologie des bitumes: Prédiction des résultats des tests de fluage BBR à partir des résultats de module complexe. Bull. Des. Lab. Des. Ponts Et Chaussées 2004, 252–253, 3–15. [Google Scholar]
  15. Olard, F.; Di Benedetto, H.; Dony, A.; Vaniscote, J.C. Properties of bituminous mixtures at low temperatures and relations with binder characteristics. Mater. Struct. 2005, 38, 121–126. [Google Scholar] [CrossRef]
  16. Marasteanu, M.O. Inter-Conversions of the Linear Viscoelastic Functions Used for the Rheological Characterization of Asphalt Binders. Ph.D. Thesis, Pennsylvania State University, Philadelphia, PA, USA, 1999. [Google Scholar]
  17. Marques, S.P.C.; Creus, G.J. Computational Viscoelasticity; Springer Science & Business Media: Berlin, Germany, 2012; pp. 11–28. [Google Scholar]
  18. 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]
  19. Tiouajni, S.; Di Benedetto, H.; Sauzéat, C.; Pouget, S. Approximation of linear viscoelastic model in the 3 dimensional case with mechanical analogues of finite size: Application to bituminous materials. Road Mater. Pavement Des. 2011, 12, 897–930. [Google Scholar] [CrossRef]
  20. Carret, J.C.; Falchetto, A.C.; Marasteanu, M.O.; Di Benedetto, H.; Wistuba, M.P.; Sauzeat, C. Comparison of rheological parameters of asphalt binders obtained from bending beam rheometer and dynamic shear rheometer at low temperatures. Road Mater. Pavement Des. 2015, 16, 211–227. [Google Scholar] [CrossRef]
  21. 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]
  22. 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]
  23. Ferry, J.D. Viscoelastic Properties of Polymers; John Wiley & Sons: Somerset, NJ, USA, 1980. [Google Scholar]
  24. Perraton, D.; Di Benedetto, H.; Sauzéat, C.; Hofko, B.; Graziani, A.; Nguyen, Q.T.; Grenfell, J. 3Dim experimental investigation of linear viscoelastic properties of bituminous mixtures. Mater. Struct. 2016, 49, 4813–4829. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) The generalized Kelvin–Voigt model (GKV); (b) 2S2P1D model.
Figure 1. (a) The generalized Kelvin–Voigt model (GKV); (b) 2S2P1D model.
Materials 13 01945 g001
Figure 2. The work plan.
Figure 2. The work plan.
Materials 13 01945 g002
Figure 3. The gradation of asphalt mix (ESG-10).
Figure 3. The gradation of asphalt mix (ESG-10).
Materials 13 01945 g003
Figure 4. Sample preparation for complex modulus test, direct tensile test and direct compression test.
Figure 4. Sample preparation for complex modulus test, direct tensile test and direct compression test.
Materials 13 01945 g004
Figure 5. Sample preparation for indirect tensile test.
Figure 5. Sample preparation for indirect tensile test.
Materials 13 01945 g005
Figure 6. Complex modulus results in (a) Cole–Cole representation and (b) black space.
Figure 6. Complex modulus results in (a) Cole–Cole representation and (b) black space.
Materials 13 01945 g006
Figure 7. (a) Master curve of complex modulus. (b) the error between 2S2P1D model and the experimental points.
Figure 7. (a) Master curve of complex modulus. (b) the error between 2S2P1D model and the experimental points.
Materials 13 01945 g007
Figure 8. Creep compliance predicted with the GKV model and with the 2S2P1D coupled to Kopelman approximation.
Figure 8. Creep compliance predicted with the GKV model and with the 2S2P1D coupled to Kopelman approximation.
Materials 13 01945 g008
Figure 9. The error between D(t)GKV model and D(t) from 2S2P1D coupled to Kopelman approximation.
Figure 9. The error between D(t)GKV model and D(t) from 2S2P1D coupled to Kopelman approximation.
Materials 13 01945 g009
Figure 10. Creep loading for direct tensile and direct compression tests.
Figure 10. Creep loading for direct tensile and direct compression tests.
Materials 13 01945 g010
Figure 11. (a) Results of direct tensile and direct compression creep tests and (b) error analysis.
Figure 11. (a) Results of direct tensile and direct compression creep tests and (b) error analysis.
Materials 13 01945 g011aMaterials 13 01945 g011b
Figure 12. EXP shift factor vs. the inverse of Williams–Landel–Ferry (WLF) equation.
Figure 12. EXP shift factor vs. the inverse of Williams–Landel–Ferry (WLF) equation.
Materials 13 01945 g012
Figure 13. Master curve of creep compliance at 0 °C.
Figure 13. Master curve of creep compliance at 0 °C.
Materials 13 01945 g013
Figure 14. Creep loading for indirect tensile test at (a) −20 °C, (b) −10 °C, (c) 0 °C, (d) IDT configuration
Figure 14. Creep loading for indirect tensile test at (a) −20 °C, (b) −10 °C, (c) 0 °C, (d) IDT configuration
Materials 13 01945 g014
Figure 15. X/Y vs. test temperature.
Figure 15. X/Y vs. test temperature.
Materials 13 01945 g015
Figure 16. (a) Isotherms of IDT creep compliance. (b) IDT creep compliance vs. direct creep compliance.
Figure 16. (a) Isotherms of IDT creep compliance. (b) IDT creep compliance vs. direct creep compliance.
Materials 13 01945 g016aMaterials 13 01945 g016b
Figure 17. Master curve of creep compliance at 0 °C.
Figure 17. Master curve of creep compliance at 0 °C.
Materials 13 01945 g017
Table 1. Forces and stresses amplitude applied for direct tensile and direct compression tests.
Table 1. Forces and stresses amplitude applied for direct tensile and direct compression tests.
Temperature (°C)Force Max (N)Stress Max (MPA)
−2018000.41
−1015000.34
012000.27
105000.11
Table 2. calculation of loss in |E*|.
Table 2. calculation of loss in |E*|.
Modulus Loss
−20 °C−10 °C0 °C10 °C
10 Hz0.32%0.02%−0.13%1.50%
3 Hz0.19%0.51%0.51%2.41%
1 Hz0.39%0.53%0.17%2.11%
0.3 Hz0.05%0.00%−1.37%2.17%

Share and Cite

MDPI and ACS Style

Daoudi, A.; Perraton, D.; Dony, A.; Carter, A. From Complex Modulus E* to Creep Compliance D(t): Experimental and Modeling Study. Materials 2020, 13, 1945. https://doi.org/10.3390/ma13081945

AMA Style

Daoudi A, Perraton D, Dony A, Carter A. From Complex Modulus E* to Creep Compliance D(t): Experimental and Modeling Study. Materials. 2020; 13(8):1945. https://doi.org/10.3390/ma13081945

Chicago/Turabian Style

Daoudi, Abdeldjalil, Daniel Perraton, Anne Dony, and Alan Carter. 2020. "From Complex Modulus E* to Creep Compliance D(t): Experimental and Modeling Study" Materials 13, no. 8: 1945. https://doi.org/10.3390/ma13081945

APA Style

Daoudi, A., Perraton, D., Dony, A., & Carter, A. (2020). From Complex Modulus E* to Creep Compliance D(t): Experimental and Modeling Study. Materials, 13(8), 1945. https://doi.org/10.3390/ma13081945

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