# 2S2P1D Model Calibration Error from User Panel for One Bitumen and One Bituminous Mixture

Authors: S. Mangiafico; C. Sauzéat; H. Di Benedetto
Journal: Advances in Materials Science and Engineering
 (2019)

---

## Abstract

The objective of this study is to analyse the differences between experimental LVE properties of both a straight-run bitumen and a bituminous mixture and simulations with analogical 2S2P1D (2 Springs, 2 Parabolic elements, and 1 Dashpot) model fitted by 14 different users. Data for the bitumen consisted of isotherms ofG∗ and φ obtained from DSR complex modulus tests at 12 different temperatures ranging from −29.9°C to 60.0°C and frequencies ranging from 6.3 to 40 Hz, for a total of 60 data points. Data for the bituminous mixture consisted of isotherms of E∗ and φ obtained from strain-controlled traction/compression complex modulus tests at 8 different temperatures ranging from −29.7°C to 38.8°C and frequencies ranging from 0.01 to 10 Hz, for a total of 55 data points. All users worked independently and for the same time duration of one hour to fit the 2S2P1D model on both sets of data. Successful simulations of experimental data of both bitumen and mixture were generally obtained by all the users over the whole range of frequencies and temperatures, regardless of their familiarity and experience with the model. The accuracy of the model to fit experimental data is all the more evident if the great spans of complex modulus (G∗ of the bitumen between 10−2 and 103 MPa, E∗ of the mixture between 10 and 40000 MPa) are considered. The obtained results highlight the convenience of 2S2P1D model to perform multiscale modelling of LVE behaviour of bituminous materials, from bitumens to mixtures.


---

## Body

## 1. Introduction

Modelling is an essential step in science. Several model-fitting methods, as well as sensitivity analysis methods, have been developed and are used in various scientific and engineering domains ([1–6], among others). The cited references offer a limited and nonexhaustive list of examples of existing methods, such as genetic algorithms, finite difference methods, direct differentiation methods, adjoint variable methods, etc.In the domain of civil engineering, a wide range of models are available to simulate the various types of mechanical behaviour of construction materials. The 2S2P1D (2 Springs, 2 Parabolic creep elements, and 1 Dashpot) is a linear viscoelastic (LVE) analogical model developed at the University of Lyon/ENTPE to simulate the behaviour of all types of bituminous materials (bitumens, mastics, and bituminous mixtures) [7]. More details on the model are given in Section 2.This model has gained attention in the scientific and technical community of the domain. The development of an accurate algorithm to automatically fit the model on a given set of experimental data by optimizing model constants is under way. Some studies are found in the literature where the model is calibrated by applying least square method to minimize the gap between experimental and simulated values of a specific LVE property (as an example, norm of complex modulus of binders was considered in [8]). However, its calibration is still generally performed manually by individual users. The model is fitted on experimental data by adapting the values of its constants through a trial and error procedure, until satisfactory simulations are obtained and visually judged [9].To the extent of the knowledge of the authors, no study has been published yet on the variability of 2S2P1D simulations of the same materials obtained by different users.The objective of this study was precisely to perform an analysis of the differences between experimental LVE properties of both a straight-run bitumen and a bituminous mixture (respectively,G∗ and E∗) and 2S2P1D simulations obtained by various users with different levels of expertise of the model.


## 2. 2S2P1D Analogical Model

As indicated by the acronym, the analogical 2S2P1D model (Figure1) is a partial derivative model consisting of an assembly of two linear elastic springs, two parabolic creep elements, and one Newtonian dashpot. For a given material with a LVE behaviour, the model can be used to simulate its complex modulus E∗ according to equation (1), where the variables are as follows: ω is the pulsation, related to frequency f as ω=2πf; E00 is the asymptotic static modulus, for ω⟶0; E0 is the asymptotic glassy modulus, for ω⟶∞; δ, k, and ℎ are dimensionless constants related to the two parabolic creep element, with k<h; β is a dimensionless constant related to the Newtonian dashpot as shown in equation (2); and τ is a characteristic time and the only constant of the model depending on temperature T.Figure 1
Analogical scheme of 2S2P1D model.If the time-temperature superposition principle (TTSP) is respected, temperature shift factorsaT at any temperature T can be obtained as in equation (3), where τ0 is the value of τ at the reference temperature Tref. Shift factors aT can be fitted with the Williams–Landel–Ferry equation (4), as a function of constants C1 and C2 and reference temperature Tref.(1)E∗=E00+E0−E001+δiωτ−k+iωτ−h+iωβτ−1,(2)η=E0−E00βτ,(3)aT=ττ0,(4)logaT=−C1T−TrefC2+T−Tref.The model can be used to simulate complex shear modulusG∗ of the same material, by replacing constants E00 and E0 with the corresponding constants G00 and G0 (Figure 1), respectively, equal to the static and glassy shear complex modulus, as in the following equation:(5)G∗=G00+G0−G001+δiωτ−k+iωτ−h+iωβτ−1.Complex modulusE∗ can be decomposed into its norm (E∗) and phase angle (φ), as well as its real (E1) and imaginary (E2) parts, as in equations (6) and (7). The same can be done for complex shear modulus G∗, as in equations (8) and (9).(6)E∗=E1+iE2=E∗cosφ+iE∗sinφ=E∗eiφ,(7)E∗=E12+E22,φ=arctanE2E1,(8)G∗=G1+iG2=G∗cosφ+iGsinφ=G∗eiφ,(9)G∗=G12+G22,φ=arctanG2G1.In this paper, complex shear modulusG∗ was used to simulate the LVE behaviour of a bitumen, while complex modulus E∗ was used for a bituminous mixture. An extension of the model has been proposed to take into account the tridimensional linear viscoelastic behaviour of materials through their complex Poisson’s ratio [10, 11]. In this paper, only the unidimensional formulation of the model has been studied.


## 3. Description of Input Experimental Data and Calibration Panel

The LVE behaviour of a straight-run unaged bitumen and of a bituminous mixture was simulated using the 2S2P1D model. The model was fitted on experimental data available for these two materials.Data for the bitumen consisted of isotherms ofG∗ and φ obtained from DSR complex modulus tests at 12 different temperatures ranging from −29.90°C to 60.00°C and frequencies ranging from 6.3 to 40 Hz, for a total of 60 data points (Figure 2). Data for the bituminous mixture consisted of isotherms of E∗ and φ obtained from complex modulus tests (strain-controlled traction/compression sinusoidal loading on cylindrical specimens with a 75 mm diameter and 150 mm length) at 8 different temperatures ranging from −29.65°C to 38.83°C and frequencies ranging from 0.01 to 10 Hz, for a total of 55 data points (Figure 3).Figure 2
Experimental data of the bitumen given to calibration panel: isotherms of norm (a) and phase angle (b) of complex shear modulusG∗.
(a)
(b)Figure 3
Experimental data of the bituminous mixture given to calibration panel: isotherms of norm (a) and phase angle (b) of complex modulusE∗.
(a)
(b)Data for the two materials were given to a calibration panel, composed of 14 users (including master and PhD students, postdoctoral fellows, senior academic researchers, and professionals in private companies) with different levels of knowledge, familiarity, and experience of the 2S2P1D model, from beginners (having learnt the theory of the model and practicing with it for less than one week) to experienced users (with more than 10 years of knowledge and extensive practice of the model). All users were given one hour to perform the calibration of the model for both materials (bitumen and mixture). All users worked independently, without any communication with the other components of the panel.Before fitting the 2S2P1D model on each material, every user had to shift the isotherms and obtain master curves of norm and phase angle of complex modulus (G∗ and E∗ for the bitumen and the mixture, respectively). The reference temperatures of 10.00°C for the bitumen and 9.53°C for the mixture were imposed. For each i-th user, values of temperature shift factors aTi at all temperatures were first obtained by manually shifting isotherms and visually evaluating the goodness of the overlap. The WLF equation was then fitted on all values of aTi obtained for each material, therefore obtaining the values of C1 and C2 at the imposed reference temperature. For this reason, for the same material, every user obtained a different set of WLF constants and, therefore, slightly different master curves. The fitting of the 2S2P1D was performed on master curves obtained using shift factors calculated with the WLF equation.At the end of the imposed time (one hour), each user was asked to give his/her values of the seven 2S2P1D constants (G00 or E00, G0 or E0, δ, k, h, β, and τ) and the two WLF constants (C1 and C2; Tref was imposed) for each material.The same Microsoft Excel® workbook, specifically designed to perform the fit, was given to all users. Several spreadsheets allow users to input data and compare experimental data points and simulated curves (WLF equation, Cole-Cole, Black and master curves).


## 4. Analysis of Results from 2S2P1D Calibration Panel

In this section, experimental data points and 2S2P1D simulations of all users are compared, for both materials. A quantitative estimation of the relative (forG∗ and E∗) and absolute (for φ) differences between data points and simulations is performed.


### 4.1. 2S2P1D Fitting for a Bitumen

Table1 reports values of WLF and 2S2P1D constants obtained by all users of the calibration panel for the bitumen. Figure 4 shows a comparison of the simulations of the LVE behaviour of the bitumen obtained by all the users and the experimental data points. In particular, Figures 4(a) and 4(b) show, respectively, Black and Cole-Cole spaces, Figures 4(c) and 4(d) show master curves of, respectively, G∗ and φ, and Figure 4(e) shows the various WLF curves obtained by the users, together with the average WLF curve, according to the constants shown in Table 1. Experimental points shown in master curves of Figures 4(c) and 4(d) are shifted using shift factors calculated with the mentioned average WLF curve, with the only purpose of having a graphical representation.Table 1
WLF and 2S2P1D constants obtained for the bitumen by all the users of the calibration panel (Tref = 10.00°C).
C
1 (−)
C
2 (°C)
G
00 (MPa)
G
0 (MPa)
δ (−)
k (−)
h (−)
β (−)
τ (s)
log
τ (−)
User 1
22.98
165.49
0
877
2.38
0.223
0.560
115
4.10 × 10−4
−3.39
User 2
16.89
120.95
0
900
1.90
0.195
0.500
200
1.50 × 10−4
−3.82
User 3
22.56
169.33
0
900
1.79
0.195
0.510
250
1.90 × 10−4
−3.72
User 4
24.16
178.59
0
888
1.95
0.205
0.515
230
2.00 × 10−4
−3.70
User 5
19.05
141.39
0
900
2.10
0.204
0.550
200
2.70 × 10−4
−3.57
User 6
18.30
137.00
0
870
2.60
0.230
0.600
100
5.50 × 10−4
−3.26
User 7
28.83
200.18
0
880
2.30
0.217
0.560
140
3.00 × 10−4
−3.52
User 8
30.86
196.09
0
900
2.20
0.210
0.550
200
3.40 × 10−4
−3.47
User 9
32.31
235.19
0
900
1.88
0.200
0.500
300
1.80 × 10−4
−3.74
User 10
32.17
223.30
0
893
2.00
0.204
0.550
220
2.50 × 10−4
−3.60
User 11
18.03
124.75
0
862
2.70
0.235
0.608
80
6.00 × 10−4
−3.22
User 12
25.37
179.43
0
870
2.53
0.226
0.580
115
4.80 × 10−4
−3.32
User 13
23.15
170.33
0
870
2.60
0.230
0.580
120
5.00 × 10−4
−3.30
User 14
26.05
203.21
0
900
1.85
0.197
0.505
400
2.00 × 10−4
−3.70
Mean
24.34
174.66
0
886
2.20
0.212
0.548
191
3.30 × 10−4
−3.52
St. dev.
5.25
35.11
0
14
0.32
0.014
0.037
88
1.52 × 10−4
0.20
COV
21.56%
20.10%
—
1.62%
14.48%
6.72%
6.73%
46.32%
45.91%
−5.68%
Note. Coefficient of variation (COV): standard deviation/mean.Figure 4
Results of 2S2P1D calibration panel for the bitumen: (a) Black space; (b) Cole-Cole plan; (c) master curves ofG∗ at 10.00°C; (d) master curves of φ at 10.00°C; (e) WLF equation used by users to obtain master curves in (c and d) and average WLF curve.Apart fromG00, imposed equal to zero by all user, the constant with the lowest coefficient of variation (COV, defined as the ratio between standard deviation and mean) is G0 (1.62%), followed by the two constants k and h associated with the parabolic elements (lower than 7%). Although a high COV is observed for the characteristic time τ, the variability of this constant is more rigorously evaluated by considering its logarithm log τ (whose COV is lower than 6% in absolute value), as it can be easily understood by taking into account its relationship with temperature and viscosity (equations (2) and (3)).Generally satisfactory approximations of experimental points are obtained with 2S2P1D simulations of all users. At low frequency/high temperature, slight underestimation ofG∗ and overestimation of φ are observed for all users. However, these differences between the model and the experimental points are negligible in absolute value (less than 100 kPa for G∗).In order to have a quantitative estimation of the differences between the 2S2P1D simulation obtained by eachi-th user and experimental points, relative errors δG1i,δG2i and δG∗i (respectively for G1, G2 and G∗) and absolute error δφi (for φ) were calculated for each equivalent frequency aTiω (60 data points) according to the following equations:(10)δG1iaTiω=G1,expaTiω−G1,2S2P1DiaTiωG1,expaTiω,(11)δG2iaTiω=G2,expaTiω−G2,2S2P1DiaTiωG2,expaTiω,(12)δG∗iaTiω=Gexp∗aTiω−G2S2P1Di∗aTiωGexp∗aTiω,(13)δφiaTiω=φexpaTiω−φ2S2P1DiaTiω,where G1,expaTiω, G2,expaTiω, Gexp∗aTiω, and φexpaTiω are experimental values of, respectively, G1, G2, G∗, and φ at the equivalent frequency aTiω. G1,2S2P1DiaTiω, G2,2S2P1DiaTiω, G2S2P1Di∗aTiω, and φ2S2P1DiaTiω are values of, respectively, G1, G2, G∗, and φ at the equivalent frequency aTiω, calculated with the 2S2P1D according to the constants of the i-th user reported in Table 1.Master curves ofδG1i,δG2i, δG∗i, and δφi are shown in Figure 5. Bold lines represent 10% and 5° limits for, respectively, relative and absolute errors. In order to avoid confusion, it is very important to highlight that errors δG1i,δG2i, δG∗i, and δφi were calculated for each i-th user at equivalent frequencies obtained with shift factors aTi calculated according to his/her WLF constants shown in Table 1. For this reason, the total range of equivalent frequencies is slightly different for each user. The results confirm the qualitative judgement made from Figure 4. The relative errors observed for equivalent frequencies lower than approximately 101 Hz are negligible because of the low values of G1, G2, and G∗ at these frequencies (lower than 100 kPa). The error for phase angle φ is generally lower than 5°, apart from some exceptions for some users.Figure 5
Errors between 2S2P1D simulations of all users and experimental data for the bitumen (Tref = 10.00°C): relative errors for G1 (a), G2 (b), and G∗ (c) and absolute error for φ (d). Bold lines represent 10% and 5° limits for, respectively, relative and absolute errors.
(a)
(b)
(c)
(d)The accuracy of the model to fit experimental data is all the more evident if the great spans of equivalent frequencies (between 10−6 and 1010 Hz at the reference temperature of 10°C) and G∗ (between 10−2 and 103 MPa) are taken into consideration. The three most experienced users (2, 5, and 12) of the panel obtained particularly accurate approximations of experimental data, as shown by their corresponding master curves (at 10°C) of relative errors for G∗ and absolute errors for φ in Figure 6. This observation confirms the ability of the model to accurately simulate LVE behaviour of bitumens.Figure 6
Errors between 2S2P1D simulations of the three most experienced users of the panel and experimental data for the bitumen (Tref = 10.00°C): relative error G∗ (a) and absolute errors for φ (b). Bold lines represent 10% and 5° limits for, respectively, relative and absolute errors.
(a)
(b)For each user, global error parameters over the whole range of frequencies were calculated for each of the considered relative and absolute errors. Equations (14)–(17) were used to calculate average values (δG1i¯, δG2i¯, δG∗i¯, and δφi¯) and standard deviations (μG1i, μG2i, μG∗i, and μφi) of relative and absolute errors of each user at the N considered equivalent frequencies aTi (N = 60), plotted in histograms of Figure 7. Error bars are plotted according to standard deviation values.(14)δG1i¯=∑aTiωδG1iaTiωN,μG1i=∑aTiωδG1iaTiω−δG1i¯2N,(15)δG2i¯=∑aTiωδG2iaTiωN,μG2i=∑aTiωδG2iaTiω−δG2i¯2N,(16)δG∗i¯=∑aTiωδG∗iaTiωN,μG∗i=∑aTiωδG∗iaTiω−δG∗i¯2N,(17)δφi¯=∑aTiωδφiaTiωN,μφi=∑aTiωδφiaTiω−δφi¯2N.Figure 7
Histograms of global errors of 2S2P1D model fitting of all users for the bitumen (Figure4), calculated according to equations (14)–(17): (a) relative errors for G1 and G2 (δG1i¯ and δG2i¯, respectively); (b) relative error for G∗ (δG∗i¯) and absolute error for φ (δφi¯). Error bars are plotted according to standard deviations μG1i, μG2i, μG∗i, and μφi.
(a)
(b)A further error estimation was performed by considering absolute values of errorsδG1i,δG2i, δG∗i, and δφi. Therefore, equations (18)–(21) were used to calculate average values (ΔG1i¯, ΔG2i¯, ΔG∗i¯, and Δφi¯) and standard deviations (MG1i, MμG2i, MμG∗i, and Mφi) of absolute values of relative and absolute errors of each user at the N considered equivalent frequencies aTi (N = 60), plotted in histograms of Figure 8. As in Figure 7, error bars are plotted according to standard deviation values.(18)ΔG1i¯=∑aTiωδGiaTiωN,MG1i=∑aTiωδG1iaTiω−ΔG1i¯2N,(19)ΔG2i¯=∑aTiωδG2iaTiωN,MG2i=∑aTiωδG2iaTiω−ΔG2i¯2N,(20)ΔG∗i¯=∑aTiωδG∗iaTiωN,MG∗i=∑aTiωδG∗iaTiω−ΔG∗i¯2N,(21)Δφi¯=∑aTiωδφiaTiωN,Mφi=∑aTiωδφiaTiω−Δφi¯2N.Figure 8
Histograms of global errors of 2S2P1D model fitting of all users for the bitumen (Figure4), calculated according to equations (18)–(21): (a) relative errors for G1 and G2 (ΔG1i¯ and ΔG2i¯, respectively); (b) relative error for G∗ (ΔG∗i¯) and absolute error for φ (Δφi¯). Error bars are plotted according to standard deviations MG1i, MμG2i, MμG∗i, and Mφi.
(a)
(b)Apart from some exceptions regarding few users, global error parameters for phase angle are generally lower than 5° even considering the corresponding standard deviations. Higher error parameters are found forG1, G2, and G∗, especially when absolute values of errors are considered. However, these parameters were calculated by arithmetically averaging all the 60 data points over the whole frequency range available. For this reason, the global error parameters are affected by the high relative errors found at low frequency/high temperature. As already discussed, in this part of the frequency spectrum, the relative errors are not particularly important, given the low values of G1, G2, and G∗ (lower than 100 kPa).Finally, overall error parameters were calculated for the whole calibration panel by averaging errors obtained for the 14 users, according to equations (22)–(29). The obtained parameters (δG1^, δG2^, δG∗^ and δφ^; ΔG1^, ΔG2^, ΔG∗^, and Δφ^ for absolute values of errors) and corresponding standard deviations (μG1^, μG2^, μG∗^, and μφ^; MG1^, MG2^, MG∗^, and Mφ^ for absolute values of errors) are reported in Table 2.(22)δG1^=∑i∑aTiωδG1iaTiω14N,μG1^=∑i∑aTiωδGiaTiω−δG1i^214N,(23)δG2^=∑i∑aTiωδG2iaTiω14N,μG2^=∑i∑aTiωδG2iaTiω−δG2i^214N,(24)δG∗^=∑i∑aTiωδG∗iaTiω14N,μG∗^=∑i∑aTiωδG∗iaTiω−δG∗i^214N,(25)δφ^=∑i∑aTiωδφiaTiω14N,μφ^=∑i∑aTiωδφiaTiω−δφi^214N,(26)ΔG1^=∑i∑aTiωδG1iaTiω14N,MG1^=∑i∑aTiωδG1iaTiω−ΔG1i^214N,(27)ΔG2^=∑i∑aTiωδG2iaTiω14N,MG2^=∑i∑aTiωδG2iaTiω−ΔG2i^214N,(28)ΔG∗^=∑i∑aTiωδG∗iaTiω14N,MG∗^=∑i∑aTiωδG∗iaTiω−ΔG∗i^214N,(29)Δφ^=∑i∑aTiωδφiaTiω14N,Mφ^=∑i∑aTiωδφiaTiω−Δφi^214N.Table 2
Overall global errors of 2S2P1D model fitting for a bitumen, calculated for the whole calibration panel according to equations (22)–(29).
δ
G
1
^
μ
G
1
^
δ
G
2
^
μ
G
2
^
δ
G
∗
^
μ
G
∗
^
δ
φ
^
μ
φ
^
Δ
G
1
^
M
G
1
^
Δ
G
2
^
M
G
2
^
Δ
G
∗
^
M
G
∗
^
Δ
φ
^
M
φ
^
−12.46%
56.13%
−6.27%
37.86%
−9.05%
40.49%
0.60°
3.40°
35.80%
44.99%
25.36%
28.80%
26.39%
32.02%
2.49°
2.39°


### 4.2. 2S2P1D Fitting for a Bituminous Mixture

Table3 reports values of WLF and 2S2P1D constants obtained by all users of the calibration panel for the bituminous mixture. Figure 9 shows a comparison of the simulations of the LVE behaviour of the mixture obtained by all the users and the experimental data points. In particular, Figures 9(a) and 9(b) show, respectively, Black and Cole-Cole spaces, Figures 9(c) and 9(d) show master curves of, respectively, E∗ and φ, and Figure 9(e) shows the various WLF curves obtained by the users, together with the average WLF curve, according to the constants shown in Table 3. As for bitumen in Section 4.1, experimental points shown in master curves of Figures 9(c) and 9(d) are shifted using shift factors calculated with the mentioned average WLF curve, in order to have a graphical representation.Table 3
WLF and 2S2P1D constants obtained for the bituminous mixture by all the users of the calibration panel (Tref = 9.53°C).
C
1 (−)
C
2 (°C)
E
00 (MPa)
E
0 (MPa)
δ (−)
k (−)
h (−)
β (−)
τ (s)
log
τ (−)
User 1
23.42
154.38
26
40700
2.26
0.189
0.593
180
2.10 × 10−1
−0.68
User 2
26.20
170.02
20
40900
2.00
0.185
0.545
400
1.00 × 10−1
−1.00
User 3
30.63
203.51
20
41300
1.95
0.178
0.563
200
9.00 × 10−2
−1.05
User 4
23.93
158.02
20
40850
2.13
0.185
0.565
230
1.00 × 10−1
−1.00
User 5
22.75
156.29
13
40500
2.30
0.195
0.600
800
1.50 × 10−1
−0.82
User 6
19.80
135.00
25
40500
2.10
0.190
0.550
180
9.50 × 10−2
−1.02
User 7
31.07
199.78
18
40800
2.15
0.188
0.560
230
9.00 × 10−2
−1.05
User 8
30.86
196.09
14
41000
2.10
0.179
0.610
200
1.40 × 10−1
−0.85
User 9
35.17
238.76
23
41000
2.00
0.180
0.550
180
9.00 × 10−2
−1.05
User 10
37.95
238.01
23
41000
2.15
0.185
0.580
230
2.50 × 10−1
−0.60
User 11
25.41
170.77
18
40850
2.15
0.189
0.595
250
1.40 × 10−1
−0.85
User 12
25.67
172.20
19
40500
2.30
0.193
0.593
240
1.30 × 10−1
−0.89
User 13
20.28
134.42
15
40700
2.30
0.190
0.620
1000000
1.70 × 10−1
−0.77
User 14
27.89
189.50
16
41000
2.28
0.188
0.595
320
1.10 × 10−1
−0.96
Mean
27.22
179.77
19
40829
2.16
0.187
0.580
71689a
1.33 × 10−1
−0.90
St. dev.
5.36
33.00
4
233
0.12
0.005
0.024
267186b
4.89 × 10−2
0.14
COV
19.70%
18.36%
20.68%
0.57%
5.56%
2.71%
4.18%
372.7%c
36.70%
−15.96%
Note. Coefficient of variation (COV): standard deviation/mean. a280, if the unusually high value of β of user 13 is ignored. b168, if the unusually high value of β of user 13 is ignored. c60.0%, if the unusually high value of β of user 13 is ignored.Figure 9
Results of 2S2P1D calibration panel for a mixture: (a) Cole-Cole plan; (b) Black space; (c) master curves ofG∗ at 9.53°C; (d) master curves of φ at 9.53°C; (e) WLF equation used by users to obtain master curves in (c and d) and average WLF curve.A very low value (0.57%) is found for the COV ofE00. Also, similarly to what was found for the bitumen, COV of constants k and h is quite low (lower than 5%). If the unusually high value of β chosen by user 13 is neglected, the COV for this constant is equal to 60.0%, of the same order of magnitude of the COV for the same constant found for the bitumen. Concerning this point, it should be highlighted that the value of β is directly related to the viscosity of the linear dashpot of the analogical model. For this reason, its variations should be considered in logarithmic scale; that is, a variation of this constant from 100 to 200 has a more important effect on the LVE simulation than a variation from 900 to 1000. In the case of the 2S2P1D fitting carried out by user 13, it can be easily verified that a relatively small variation occurs to the corresponding simulation if the value of β is fixed to 1000 instead of 1000000.Overall, one important conclusion that can be drawn from the results of the calibration panel is that the variability (expressed by the COV) of constantsG00 (for the bitumen), E00 and E0 (for the mixture), δ, k, h, and log τ (for both) is approximately equal or lower than 20%.2S2P1D simulations of all users satisfactorily approximate experimental points. At very low frequency/high temperature (less than 10−3 Hz), a negligible underestimation of E∗ and φ is observed for all users. However, these differences are negligible in absolute value (less than 20 MPa for E∗) and concern only very few experimental points of the isotherms obtained at 38.83°C. As in the case of the bitumen, satisfactory fitting of experimental data is observed over the whole spans of equivalent frequencies (between 10−6 and 1010 Hz at the reference temperature of 9.5°C) and E∗ (between 10 and 40000 MPa).Similarly to what was done for the bitumen, relative errorsδE1i,δE2i, and δE∗i (respectively, for E1, E2, and E∗) and absolute error δφi (for φ) were calculated for each equivalent frequency aTiω (55 data points) according to the following equations:(30)δE1iaTiω=E1,expaTiω−E1,2S2P1DiaTiωE1,expaTiω,(31)δE2iaTiω=E2,expaTiω−E2,2S2P1DiaTiωE2,expaTiω,(32)δE∗iaTiω=Eexp∗aTiω−E2S2P1Di∗aTiωEexp∗aTiω,(33)δφiaTiω=φexpaTiω−φ2S2P1DiaTiω,where the variables are listed as follows: E1,expaTiω, E2,expaTiω, Eexp∗aTiω, and φexpaTiω are experimental values of, respectively, E1, E2, E∗, and φ at the equivalent frequency aTiω; E1,2S2P1DiaTiω, E2,2S2P1DiaTiω, E2S2P1Di∗aTiω, and φ2S2P1DiaTiω are values of, respectively, E1, E2, E∗, and φ at the equivalent frequency aTiω, calculated with the 2S2P1D according to the constants of the i-th user reported in Table 1.Master curves ofδE1i,δE2i, δE∗i, and δφi are shown in Figure 10 (10% and 5° limits for, respectively, relative and absolute errors are represented by bold lines). As in the case of the bitumen, errors δE1i, δE2i, δE∗i, and δφi were calculated for each i-th user at equivalent frequencies obtained with shift factors aTi calculated according to his/her WLF constants shown in Table 3. Hence, master curves corresponding to different users have slightly different total ranges of equivalent frequencies. Differences greater than 5° between experimental and simulated φ (parameter δφi in Figure 10(d)) are observed only for the last two points of the isotherm at 38.83°C. Concerning errors in the simulation of E1, E2, and E∗, for the majority of users, values greater than 10% of parameters δE1i, δE2i, and δE∗i are observed only for equivalent frequencies lower than approximately 10−2 Hz. Only for two users, values of 2S2P1D model simulations are significantly higher than experimental data for frequencies between 10−3 and 102 Hz. The important relative errors observed for E2 at frequencies higher than 105 Hz can be neglected because of the very low values of E2 in this frequency range (lower than 2000 MPa) with respect to E∗ (higher than 30000 MPa). Similarly to what was done for the bitumen (Section 4.1), the master curves (at 9.5°C) of δE∗i and δφi obtained for the mixture by the three most experienced users (2, 5, and 12) are plotted in Figure 11. Excellent correspondences between simulations and experimental data are observed, confirming the ability of the model to simulate LVE behaviour of mixtures.Figure 10
Errors between 2S2P1D simulations of all users and experimental data for the bituminous mixture (Tref = 9.53°C): relative errors for E1 (a), E2 (b), and E∗ (c) and absolute error for φ (d). Bold lines represent 10% and 5° limits for, respectively, relative and absolute errors.
(a)
(b)
(c)
(d)Figure 11
Errors between 2S2P1D simulations of the three most experienced users of the panel and experimental data for the bituminous mixture (Tref = 9.53°C): relative error for E∗ (a) and absolute error for φ (b). Bold lines represent 10% and 5° limits for, respectively, relative and absolute errors.
(a)
(b)Average values (δE1i¯, δE2i¯, δE∗i¯, and δφi¯) and standard deviations (μE1i, μE2i, μE∗i, and μφi) of relative and absolute errors δE1i,δE2i, δE∗i, and δφi of each user at the N considered equivalent frequencies aTi (N = 55) were calculated according to equations (34)–(37) and plotted in histograms of Figure 12. Error bars are plotted according to standard deviation values.(34)δE1i¯=∑aTiωδE1iaTiωN,μE1i=∑aTiωδE1iaTiω−δE1i¯2N,(35)δE2i¯=∑aTiωδE2iaTiωN,μE2i=∑aTiωδE2iaTiω−δE2i¯2N,(36)δE∗i¯=∑aTiωδE∗iaTiωN,μE∗i=∑aTiωδE∗iaTiω−δE∗i¯2N,(37)δφi¯=∑aTiωδφiaTiωN,μφi=∑aTiωδφiaTiω−δφi¯2N.Figure 12
Histograms of global errors of 2S2P1D model fitting of all users for the bituminous mixture (Figure9), calculated according to equations (34)–(37): (a) relative errors for E1 and E2 (δE1i¯ and δE2i¯, respectively); (b) relative error for E∗ (δE∗i¯) and absolute error for φ (δφi¯). Error bars are plotted according to standard deviations μE1i, μE2i, μE∗i, and μφi.
(a)
(b)Equations (38)–(41) were used to calculate average values (ΔE1i¯, ΔE2i¯, ΔE∗i¯, and Δφi¯) and standard deviations (ME1i, MμE2i, MμE∗i, and Mφi) of absolute values of relative and absolute errors δE1i,δE2i, δE∗i, and δφi of each user at the N considered equivalent frequencies aTi (N = 55), plotted in histograms of Figure 13. As shown in Figure 12, error bars are plotted according to standard deviation values.(38)ΔE1i¯=∑aTiωδE1iaTiωN,ME1i=∑aTiωδE1iaTiω−ΔE1i¯2N,(39)ΔE2i¯=∑aTiωδE2iaTiωN,ME2i=∑aTiωδE2iaTiω−ΔE2i¯2N,(40)ΔE∗i¯=∑aTiωδE∗iaTiωN,ME∗i=∑aTiωδE∗iaTiω−ΔE∗i¯2N,(41)Δφi¯=∑aTiωδφiaTiωN,Mφi=∑aTiωδφiaTiω−Δφi¯2N.Figure 13
Histograms of global errors of 2S2P1D model fitting of all users for the bituminous mixture (Figure9), calculated according to equations (38)–(41): (a) relative errors for E1 and E2 (ΔE1i¯ and ΔE2i¯, respectively); (b) relative error for E∗ (ΔE∗i¯) and absolute error for φ (Δφi¯). Error bars are plotted according to standard deviations ME1i, MμE2i, MμE∗i, and Mφi.
(a)
(b)Global error parametersδE∗i¯ and ΔE∗i¯ (for E∗) are generally lower than 10% and parameters δφi¯ and Δφi¯ (for φ) are generally lower than 5° (even considering the corresponding standard deviations), except for the same few users. These estimates show that very satisfactory simulations were obtained by the vast majority of the users. This observation is all the more relevant because the error parameters were calculated by arithmetically averaging all the 55 data points over the whole frequency range available.Overall error parameters were calculated for the whole calibration panel by averaging errors obtained for the 14 users, according to equations (42)–(49). The obtained parameters (δE1^, δE2^, δE∗^, and δφ^; ΔE1^, ΔE2^, ΔE∗^, and Δφ^ for absolute values of errors) and corresponding standard deviations (μE1^, μE2^, μE∗^, and μφ^; ME1^, ME2^, ME∗^, and Mφ^ for absolute values of errors) are reported in Table 4.(42)δE1^=∑i∑aTiωδE1iaTiω14N,μE1^=∑i∑aTiωδE1iaTiω−δE1i^214N,(43)δE2^=∑i∑aTiωδE2iaTiω14N,μE2^=∑i∑aTiωδE2iaTiω−δE2i^214N,(44)δE∗^=∑i∑aTiωδE∗iaTiω14N,μE∗^=∑i∑aTiωδE∗iaTiω−δE∗i^214N,(45)δφ^=∑i∑aTiωδφiaTiω14N,μφ^=∑i∑aTiωδφiaTiω−δφi^214N,(46)ΔE1^=∑i∑aTiωδE1iaTiω14N,ME1^=∑i∑aTiωδE1iaTiω−ΔE1i^214N,(47)ΔE2^=∑i∑aTiωδE2iaTiω14N,ME2^=∑i∑aTiωδE2iaTiω−ΔE2i^214N,(48)ΔE∗^=∑i∑aTiωδE∗iaTiω14N,ME∗^=∑i∑aTiωδE∗iaTiω−ΔE∗i^214N,(49)Δφ^=∑i∑aTiωδφiaTiω14N,Mφ^=∑i∑aTiωδφiaTiω−Δφi^214N.Table 4
Global errors of 2S2P1D model fitting for a bituminous mixture, calculated for the whole calibration panel according to equations (42)–(49).
δ
E
1
^
μ
E
1
^
δ
E
2
^
μ
E
2
^
δ
E
∗
^
μ
E
∗
^
δ
φ
^
μ
φ
^
Δ
E
1
^
M
E
1
^
Δ
E
2
^
M
E
2
^
Δ
E
∗
^
M
E
∗
^
Δ
φ
^
M
φ
^
1.70%
12.65%
7.06%
15.71%
2.77%
12.94%
1.02°
3.14°
8.10%
9.86%
10.86%
13.38%
8.15%
10.42%
1.48°
2.95°


## 4.1. 2S2P1D Fitting for a Bitumen

Table1 reports values of WLF and 2S2P1D constants obtained by all users of the calibration panel for the bitumen. Figure 4 shows a comparison of the simulations of the LVE behaviour of the bitumen obtained by all the users and the experimental data points. In particular, Figures 4(a) and 4(b) show, respectively, Black and Cole-Cole spaces, Figures 4(c) and 4(d) show master curves of, respectively, G∗ and φ, and Figure 4(e) shows the various WLF curves obtained by the users, together with the average WLF curve, according to the constants shown in Table 1. Experimental points shown in master curves of Figures 4(c) and 4(d) are shifted using shift factors calculated with the mentioned average WLF curve, with the only purpose of having a graphical representation.Table 1
WLF and 2S2P1D constants obtained for the bitumen by all the users of the calibration panel (Tref = 10.00°C).
C
1 (−)
C
2 (°C)
G
00 (MPa)
G
0 (MPa)
δ (−)
k (−)
h (−)
β (−)
τ (s)
log
τ (−)
User 1
22.98
165.49
0
877
2.38
0.223
0.560
115
4.10 × 10−4
−3.39
User 2
16.89
120.95
0
900
1.90
0.195
0.500
200
1.50 × 10−4
−3.82
User 3
22.56
169.33
0
900
1.79
0.195
0.510
250
1.90 × 10−4
−3.72
User 4
24.16
178.59
0
888
1.95
0.205
0.515
230
2.00 × 10−4
−3.70
User 5
19.05
141.39
0
900
2.10
0.204
0.550
200
2.70 × 10−4
−3.57
User 6
18.30
137.00
0
870
2.60
0.230
0.600
100
5.50 × 10−4
−3.26
User 7
28.83
200.18
0
880
2.30
0.217
0.560
140
3.00 × 10−4
−3.52
User 8
30.86
196.09
0
900
2.20
0.210
0.550
200
3.40 × 10−4
−3.47
User 9
32.31
235.19
0
900
1.88
0.200
0.500
300
1.80 × 10−4
−3.74
User 10
32.17
223.30
0
893
2.00
0.204
0.550
220
2.50 × 10−4
−3.60
User 11
18.03
124.75
0
862
2.70
0.235
0.608
80
6.00 × 10−4
−3.22
User 12
25.37
179.43
0
870
2.53
0.226
0.580
115
4.80 × 10−4
−3.32
User 13
23.15
170.33
0
870
2.60
0.230
0.580
120
5.00 × 10−4
−3.30
User 14
26.05
203.21
0
900
1.85
0.197
0.505
400
2.00 × 10−4
−3.70
Mean
24.34
174.66
0
886
2.20
0.212
0.548
191
3.30 × 10−4
−3.52
St. dev.
5.25
35.11
0
14
0.32
0.014
0.037
88
1.52 × 10−4
0.20
COV
21.56%
20.10%
—
1.62%
14.48%
6.72%
6.73%
46.32%
45.91%
−5.68%
Note. Coefficient of variation (COV): standard deviation/mean.Figure 4
Results of 2S2P1D calibration panel for the bitumen: (a) Black space; (b) Cole-Cole plan; (c) master curves ofG∗ at 10.00°C; (d) master curves of φ at 10.00°C; (e) WLF equation used by users to obtain master curves in (c and d) and average WLF curve.Apart fromG00, imposed equal to zero by all user, the constant with the lowest coefficient of variation (COV, defined as the ratio between standard deviation and mean) is G0 (1.62%), followed by the two constants k and h associated with the parabolic elements (lower than 7%). Although a high COV is observed for the characteristic time τ, the variability of this constant is more rigorously evaluated by considering its logarithm log τ (whose COV is lower than 6% in absolute value), as it can be easily understood by taking into account its relationship with temperature and viscosity (equations (2) and (3)).Generally satisfactory approximations of experimental points are obtained with 2S2P1D simulations of all users. At low frequency/high temperature, slight underestimation ofG∗ and overestimation of φ are observed for all users. However, these differences between the model and the experimental points are negligible in absolute value (less than 100 kPa for G∗).In order to have a quantitative estimation of the differences between the 2S2P1D simulation obtained by eachi-th user and experimental points, relative errors δG1i,δG2i and δG∗i (respectively for G1, G2 and G∗) and absolute error δφi (for φ) were calculated for each equivalent frequency aTiω (60 data points) according to the following equations:(10)δG1iaTiω=G1,expaTiω−G1,2S2P1DiaTiωG1,expaTiω,(11)δG2iaTiω=G2,expaTiω−G2,2S2P1DiaTiωG2,expaTiω,(12)δG∗iaTiω=Gexp∗aTiω−G2S2P1Di∗aTiωGexp∗aTiω,(13)δφiaTiω=φexpaTiω−φ2S2P1DiaTiω,where G1,expaTiω, G2,expaTiω, Gexp∗aTiω, and φexpaTiω are experimental values of, respectively, G1, G2, G∗, and φ at the equivalent frequency aTiω. G1,2S2P1DiaTiω, G2,2S2P1DiaTiω, G2S2P1Di∗aTiω, and φ2S2P1DiaTiω are values of, respectively, G1, G2, G∗, and φ at the equivalent frequency aTiω, calculated with the 2S2P1D according to the constants of the i-th user reported in Table 1.Master curves ofδG1i,δG2i, δG∗i, and δφi are shown in Figure 5. Bold lines represent 10% and 5° limits for, respectively, relative and absolute errors. In order to avoid confusion, it is very important to highlight that errors δG1i,δG2i, δG∗i, and δφi were calculated for each i-th user at equivalent frequencies obtained with shift factors aTi calculated according to his/her WLF constants shown in Table 1. For this reason, the total range of equivalent frequencies is slightly different for each user. The results confirm the qualitative judgement made from Figure 4. The relative errors observed for equivalent frequencies lower than approximately 101 Hz are negligible because of the low values of G1, G2, and G∗ at these frequencies (lower than 100 kPa). The error for phase angle φ is generally lower than 5°, apart from some exceptions for some users.Figure 5
Errors between 2S2P1D simulations of all users and experimental data for the bitumen (Tref = 10.00°C): relative errors for G1 (a), G2 (b), and G∗ (c) and absolute error for φ (d). Bold lines represent 10% and 5° limits for, respectively, relative and absolute errors.
(a)
(b)
(c)
(d)The accuracy of the model to fit experimental data is all the more evident if the great spans of equivalent frequencies (between 10−6 and 1010 Hz at the reference temperature of 10°C) and G∗ (between 10−2 and 103 MPa) are taken into consideration. The three most experienced users (2, 5, and 12) of the panel obtained particularly accurate approximations of experimental data, as shown by their corresponding master curves (at 10°C) of relative errors for G∗ and absolute errors for φ in Figure 6. This observation confirms the ability of the model to accurately simulate LVE behaviour of bitumens.Figure 6
Errors between 2S2P1D simulations of the three most experienced users of the panel and experimental data for the bitumen (Tref = 10.00°C): relative error G∗ (a) and absolute errors for φ (b). Bold lines represent 10% and 5° limits for, respectively, relative and absolute errors.
(a)
(b)For each user, global error parameters over the whole range of frequencies were calculated for each of the considered relative and absolute errors. Equations (14)–(17) were used to calculate average values (δG1i¯, δG2i¯, δG∗i¯, and δφi¯) and standard deviations (μG1i, μG2i, μG∗i, and μφi) of relative and absolute errors of each user at the N considered equivalent frequencies aTi (N = 60), plotted in histograms of Figure 7. Error bars are plotted according to standard deviation values.(14)δG1i¯=∑aTiωδG1iaTiωN,μG1i=∑aTiωδG1iaTiω−δG1i¯2N,(15)δG2i¯=∑aTiωδG2iaTiωN,μG2i=∑aTiωδG2iaTiω−δG2i¯2N,(16)δG∗i¯=∑aTiωδG∗iaTiωN,μG∗i=∑aTiωδG∗iaTiω−δG∗i¯2N,(17)δφi¯=∑aTiωδφiaTiωN,μφi=∑aTiωδφiaTiω−δφi¯2N.Figure 7
Histograms of global errors of 2S2P1D model fitting of all users for the bitumen (Figure4), calculated according to equations (14)–(17): (a) relative errors for G1 and G2 (δG1i¯ and δG2i¯, respectively); (b) relative error for G∗ (δG∗i¯) and absolute error for φ (δφi¯). Error bars are plotted according to standard deviations μG1i, μG2i, μG∗i, and μφi.
(a)
(b)A further error estimation was performed by considering absolute values of errorsδG1i,δG2i, δG∗i, and δφi. Therefore, equations (18)–(21) were used to calculate average values (ΔG1i¯, ΔG2i¯, ΔG∗i¯, and Δφi¯) and standard deviations (MG1i, MμG2i, MμG∗i, and Mφi) of absolute values of relative and absolute errors of each user at the N considered equivalent frequencies aTi (N = 60), plotted in histograms of Figure 8. As in Figure 7, error bars are plotted according to standard deviation values.(18)ΔG1i¯=∑aTiωδGiaTiωN,MG1i=∑aTiωδG1iaTiω−ΔG1i¯2N,(19)ΔG2i¯=∑aTiωδG2iaTiωN,MG2i=∑aTiωδG2iaTiω−ΔG2i¯2N,(20)ΔG∗i¯=∑aTiωδG∗iaTiωN,MG∗i=∑aTiωδG∗iaTiω−ΔG∗i¯2N,(21)Δφi¯=∑aTiωδφiaTiωN,Mφi=∑aTiωδφiaTiω−Δφi¯2N.Figure 8
Histograms of global errors of 2S2P1D model fitting of all users for the bitumen (Figure4), calculated according to equations (18)–(21): (a) relative errors for G1 and G2 (ΔG1i¯ and ΔG2i¯, respectively); (b) relative error for G∗ (ΔG∗i¯) and absolute error for φ (Δφi¯). Error bars are plotted according to standard deviations MG1i, MμG2i, MμG∗i, and Mφi.
(a)
(b)Apart from some exceptions regarding few users, global error parameters for phase angle are generally lower than 5° even considering the corresponding standard deviations. Higher error parameters are found forG1, G2, and G∗, especially when absolute values of errors are considered. However, these parameters were calculated by arithmetically averaging all the 60 data points over the whole frequency range available. For this reason, the global error parameters are affected by the high relative errors found at low frequency/high temperature. As already discussed, in this part of the frequency spectrum, the relative errors are not particularly important, given the low values of G1, G2, and G∗ (lower than 100 kPa).Finally, overall error parameters were calculated for the whole calibration panel by averaging errors obtained for the 14 users, according to equations (22)–(29). The obtained parameters (δG1^, δG2^, δG∗^ and δφ^; ΔG1^, ΔG2^, ΔG∗^, and Δφ^ for absolute values of errors) and corresponding standard deviations (μG1^, μG2^, μG∗^, and μφ^; MG1^, MG2^, MG∗^, and Mφ^ for absolute values of errors) are reported in Table 2.(22)δG1^=∑i∑aTiωδG1iaTiω14N,μG1^=∑i∑aTiωδGiaTiω−δG1i^214N,(23)δG2^=∑i∑aTiωδG2iaTiω14N,μG2^=∑i∑aTiωδG2iaTiω−δG2i^214N,(24)δG∗^=∑i∑aTiωδG∗iaTiω14N,μG∗^=∑i∑aTiωδG∗iaTiω−δG∗i^214N,(25)δφ^=∑i∑aTiωδφiaTiω14N,μφ^=∑i∑aTiωδφiaTiω−δφi^214N,(26)ΔG1^=∑i∑aTiωδG1iaTiω14N,MG1^=∑i∑aTiωδG1iaTiω−ΔG1i^214N,(27)ΔG2^=∑i∑aTiωδG2iaTiω14N,MG2^=∑i∑aTiωδG2iaTiω−ΔG2i^214N,(28)ΔG∗^=∑i∑aTiωδG∗iaTiω14N,MG∗^=∑i∑aTiωδG∗iaTiω−ΔG∗i^214N,(29)Δφ^=∑i∑aTiωδφiaTiω14N,Mφ^=∑i∑aTiωδφiaTiω−Δφi^214N.Table 2
Overall global errors of 2S2P1D model fitting for a bitumen, calculated for the whole calibration panel according to equations (22)–(29).
δ
G
1
^
μ
G
1
^
δ
G
2
^
μ
G
2
^
δ
G
∗
^
μ
G
∗
^
δ
φ
^
μ
φ
^
Δ
G
1
^
M
G
1
^
Δ
G
2
^
M
G
2
^
Δ
G
∗
^
M
G
∗
^
Δ
φ
^
M
φ
^
−12.46%
56.13%
−6.27%
37.86%
−9.05%
40.49%
0.60°
3.40°
35.80%
44.99%
25.36%
28.80%
26.39%
32.02%
2.49°
2.39°


## 4.2. 2S2P1D Fitting for a Bituminous Mixture

Table3 reports values of WLF and 2S2P1D constants obtained by all users of the calibration panel for the bituminous mixture. Figure 9 shows a comparison of the simulations of the LVE behaviour of the mixture obtained by all the users and the experimental data points. In particular, Figures 9(a) and 9(b) show, respectively, Black and Cole-Cole spaces, Figures 9(c) and 9(d) show master curves of, respectively, E∗ and φ, and Figure 9(e) shows the various WLF curves obtained by the users, together with the average WLF curve, according to the constants shown in Table 3. As for bitumen in Section 4.1, experimental points shown in master curves of Figures 9(c) and 9(d) are shifted using shift factors calculated with the mentioned average WLF curve, in order to have a graphical representation.Table 3
WLF and 2S2P1D constants obtained for the bituminous mixture by all the users of the calibration panel (Tref = 9.53°C).
C
1 (−)
C
2 (°C)
E
00 (MPa)
E
0 (MPa)
δ (−)
k (−)
h (−)
β (−)
τ (s)
log
τ (−)
User 1
23.42
154.38
26
40700
2.26
0.189
0.593
180
2.10 × 10−1
−0.68
User 2
26.20
170.02
20
40900
2.00
0.185
0.545
400
1.00 × 10−1
−1.00
User 3
30.63
203.51
20
41300
1.95
0.178
0.563
200
9.00 × 10−2
−1.05
User 4
23.93
158.02
20
40850
2.13
0.185
0.565
230
1.00 × 10−1
−1.00
User 5
22.75
156.29
13
40500
2.30
0.195
0.600
800
1.50 × 10−1
−0.82
User 6
19.80
135.00
25
40500
2.10
0.190
0.550
180
9.50 × 10−2
−1.02
User 7
31.07
199.78
18
40800
2.15
0.188
0.560
230
9.00 × 10−2
−1.05
User 8
30.86
196.09
14
41000
2.10
0.179
0.610
200
1.40 × 10−1
−0.85
User 9
35.17
238.76
23
41000
2.00
0.180
0.550
180
9.00 × 10−2
−1.05
User 10
37.95
238.01
23
41000
2.15
0.185
0.580
230
2.50 × 10−1
−0.60
User 11
25.41
170.77
18
40850
2.15
0.189
0.595
250
1.40 × 10−1
−0.85
User 12
25.67
172.20
19
40500
2.30
0.193
0.593
240
1.30 × 10−1
−0.89
User 13
20.28
134.42
15
40700
2.30
0.190
0.620
1000000
1.70 × 10−1
−0.77
User 14
27.89
189.50
16
41000
2.28
0.188
0.595
320
1.10 × 10−1
−0.96
Mean
27.22
179.77
19
40829
2.16
0.187
0.580
71689a
1.33 × 10−1
−0.90
St. dev.
5.36
33.00
4
233
0.12
0.005
0.024
267186b
4.89 × 10−2
0.14
COV
19.70%
18.36%
20.68%
0.57%
5.56%
2.71%
4.18%
372.7%c
36.70%
−15.96%
Note. Coefficient of variation (COV): standard deviation/mean. a280, if the unusually high value of β of user 13 is ignored. b168, if the unusually high value of β of user 13 is ignored. c60.0%, if the unusually high value of β of user 13 is ignored.Figure 9
Results of 2S2P1D calibration panel for a mixture: (a) Cole-Cole plan; (b) Black space; (c) master curves ofG∗ at 9.53°C; (d) master curves of φ at 9.53°C; (e) WLF equation used by users to obtain master curves in (c and d) and average WLF curve.A very low value (0.57%) is found for the COV ofE00. Also, similarly to what was found for the bitumen, COV of constants k and h is quite low (lower than 5%). If the unusually high value of β chosen by user 13 is neglected, the COV for this constant is equal to 60.0%, of the same order of magnitude of the COV for the same constant found for the bitumen. Concerning this point, it should be highlighted that the value of β is directly related to the viscosity of the linear dashpot of the analogical model. For this reason, its variations should be considered in logarithmic scale; that is, a variation of this constant from 100 to 200 has a more important effect on the LVE simulation than a variation from 900 to 1000. In the case of the 2S2P1D fitting carried out by user 13, it can be easily verified that a relatively small variation occurs to the corresponding simulation if the value of β is fixed to 1000 instead of 1000000.Overall, one important conclusion that can be drawn from the results of the calibration panel is that the variability (expressed by the COV) of constantsG00 (for the bitumen), E00 and E0 (for the mixture), δ, k, h, and log τ (for both) is approximately equal or lower than 20%.2S2P1D simulations of all users satisfactorily approximate experimental points. At very low frequency/high temperature (less than 10−3 Hz), a negligible underestimation of E∗ and φ is observed for all users. However, these differences are negligible in absolute value (less than 20 MPa for E∗) and concern only very few experimental points of the isotherms obtained at 38.83°C. As in the case of the bitumen, satisfactory fitting of experimental data is observed over the whole spans of equivalent frequencies (between 10−6 and 1010 Hz at the reference temperature of 9.5°C) and E∗ (between 10 and 40000 MPa).Similarly to what was done for the bitumen, relative errorsδE1i,δE2i, and δE∗i (respectively, for E1, E2, and E∗) and absolute error δφi (for φ) were calculated for each equivalent frequency aTiω (55 data points) according to the following equations:(30)δE1iaTiω=E1,expaTiω−E1,2S2P1DiaTiωE1,expaTiω,(31)δE2iaTiω=E2,expaTiω−E2,2S2P1DiaTiωE2,expaTiω,(32)δE∗iaTiω=Eexp∗aTiω−E2S2P1Di∗aTiωEexp∗aTiω,(33)δφiaTiω=φexpaTiω−φ2S2P1DiaTiω,where the variables are listed as follows: E1,expaTiω, E2,expaTiω, Eexp∗aTiω, and φexpaTiω are experimental values of, respectively, E1, E2, E∗, and φ at the equivalent frequency aTiω; E1,2S2P1DiaTiω, E2,2S2P1DiaTiω, E2S2P1Di∗aTiω, and φ2S2P1DiaTiω are values of, respectively, E1, E2, E∗, and φ at the equivalent frequency aTiω, calculated with the 2S2P1D according to the constants of the i-th user reported in Table 1.Master curves ofδE1i,δE2i, δE∗i, and δφi are shown in Figure 10 (10% and 5° limits for, respectively, relative and absolute errors are represented by bold lines). As in the case of the bitumen, errors δE1i, δE2i, δE∗i, and δφi were calculated for each i-th user at equivalent frequencies obtained with shift factors aTi calculated according to his/her WLF constants shown in Table 3. Hence, master curves corresponding to different users have slightly different total ranges of equivalent frequencies. Differences greater than 5° between experimental and simulated φ (parameter δφi in Figure 10(d)) are observed only for the last two points of the isotherm at 38.83°C. Concerning errors in the simulation of E1, E2, and E∗, for the majority of users, values greater than 10% of parameters δE1i, δE2i, and δE∗i are observed only for equivalent frequencies lower than approximately 10−2 Hz. Only for two users, values of 2S2P1D model simulations are significantly higher than experimental data for frequencies between 10−3 and 102 Hz. The important relative errors observed for E2 at frequencies higher than 105 Hz can be neglected because of the very low values of E2 in this frequency range (lower than 2000 MPa) with respect to E∗ (higher than 30000 MPa). Similarly to what was done for the bitumen (Section 4.1), the master curves (at 9.5°C) of δE∗i and δφi obtained for the mixture by the three most experienced users (2, 5, and 12) are plotted in Figure 11. Excellent correspondences between simulations and experimental data are observed, confirming the ability of the model to simulate LVE behaviour of mixtures.Figure 10
Errors between 2S2P1D simulations of all users and experimental data for the bituminous mixture (Tref = 9.53°C): relative errors for E1 (a), E2 (b), and E∗ (c) and absolute error for φ (d). Bold lines represent 10% and 5° limits for, respectively, relative and absolute errors.
(a)
(b)
(c)
(d)Figure 11
Errors between 2S2P1D simulations of the three most experienced users of the panel and experimental data for the bituminous mixture (Tref = 9.53°C): relative error for E∗ (a) and absolute error for φ (b). Bold lines represent 10% and 5° limits for, respectively, relative and absolute errors.
(a)
(b)Average values (δE1i¯, δE2i¯, δE∗i¯, and δφi¯) and standard deviations (μE1i, μE2i, μE∗i, and μφi) of relative and absolute errors δE1i,δE2i, δE∗i, and δφi of each user at the N considered equivalent frequencies aTi (N = 55) were calculated according to equations (34)–(37) and plotted in histograms of Figure 12. Error bars are plotted according to standard deviation values.(34)δE1i¯=∑aTiωδE1iaTiωN,μE1i=∑aTiωδE1iaTiω−δE1i¯2N,(35)δE2i¯=∑aTiωδE2iaTiωN,μE2i=∑aTiωδE2iaTiω−δE2i¯2N,(36)δE∗i¯=∑aTiωδE∗iaTiωN,μE∗i=∑aTiωδE∗iaTiω−δE∗i¯2N,(37)δφi¯=∑aTiωδφiaTiωN,μφi=∑aTiωδφiaTiω−δφi¯2N.Figure 12
Histograms of global errors of 2S2P1D model fitting of all users for the bituminous mixture (Figure9), calculated according to equations (34)–(37): (a) relative errors for E1 and E2 (δE1i¯ and δE2i¯, respectively); (b) relative error for E∗ (δE∗i¯) and absolute error for φ (δφi¯). Error bars are plotted according to standard deviations μE1i, μE2i, μE∗i, and μφi.
(a)
(b)Equations (38)–(41) were used to calculate average values (ΔE1i¯, ΔE2i¯, ΔE∗i¯, and Δφi¯) and standard deviations (ME1i, MμE2i, MμE∗i, and Mφi) of absolute values of relative and absolute errors δE1i,δE2i, δE∗i, and δφi of each user at the N considered equivalent frequencies aTi (N = 55), plotted in histograms of Figure 13. As shown in Figure 12, error bars are plotted according to standard deviation values.(38)ΔE1i¯=∑aTiωδE1iaTiωN,ME1i=∑aTiωδE1iaTiω−ΔE1i¯2N,(39)ΔE2i¯=∑aTiωδE2iaTiωN,ME2i=∑aTiωδE2iaTiω−ΔE2i¯2N,(40)ΔE∗i¯=∑aTiωδE∗iaTiωN,ME∗i=∑aTiωδE∗iaTiω−ΔE∗i¯2N,(41)Δφi¯=∑aTiωδφiaTiωN,Mφi=∑aTiωδφiaTiω−Δφi¯2N.Figure 13
Histograms of global errors of 2S2P1D model fitting of all users for the bituminous mixture (Figure9), calculated according to equations (38)–(41): (a) relative errors for E1 and E2 (ΔE1i¯ and ΔE2i¯, respectively); (b) relative error for E∗ (ΔE∗i¯) and absolute error for φ (Δφi¯). Error bars are plotted according to standard deviations ME1i, MμE2i, MμE∗i, and Mφi.
(a)
(b)Global error parametersδE∗i¯ and ΔE∗i¯ (for E∗) are generally lower than 10% and parameters δφi¯ and Δφi¯ (for φ) are generally lower than 5° (even considering the corresponding standard deviations), except for the same few users. These estimates show that very satisfactory simulations were obtained by the vast majority of the users. This observation is all the more relevant because the error parameters were calculated by arithmetically averaging all the 55 data points over the whole frequency range available.Overall error parameters were calculated for the whole calibration panel by averaging errors obtained for the 14 users, according to equations (42)–(49). The obtained parameters (δE1^, δE2^, δE∗^, and δφ^; ΔE1^, ΔE2^, ΔE∗^, and Δφ^ for absolute values of errors) and corresponding standard deviations (μE1^, μE2^, μE∗^, and μφ^; ME1^, ME2^, ME∗^, and Mφ^ for absolute values of errors) are reported in Table 4.(42)δE1^=∑i∑aTiωδE1iaTiω14N,μE1^=∑i∑aTiωδE1iaTiω−δE1i^214N,(43)δE2^=∑i∑aTiωδE2iaTiω14N,μE2^=∑i∑aTiωδE2iaTiω−δE2i^214N,(44)δE∗^=∑i∑aTiωδE∗iaTiω14N,μE∗^=∑i∑aTiωδE∗iaTiω−δE∗i^214N,(45)δφ^=∑i∑aTiωδφiaTiω14N,μφ^=∑i∑aTiωδφiaTiω−δφi^214N,(46)ΔE1^=∑i∑aTiωδE1iaTiω14N,ME1^=∑i∑aTiωδE1iaTiω−ΔE1i^214N,(47)ΔE2^=∑i∑aTiωδE2iaTiω14N,ME2^=∑i∑aTiωδE2iaTiω−ΔE2i^214N,(48)ΔE∗^=∑i∑aTiωδE∗iaTiω14N,ME∗^=∑i∑aTiωδE∗iaTiω−ΔE∗i^214N,(49)Δφ^=∑i∑aTiωδφiaTiω14N,Mφ^=∑i∑aTiωδφiaTiω−Δφi^214N.Table 4
Global errors of 2S2P1D model fitting for a bituminous mixture, calculated for the whole calibration panel according to equations (42)–(49).
δ
E
1
^
μ
E
1
^
δ
E
2
^
μ
E
2
^
δ
E
∗
^
μ
E
∗
^
δ
φ
^
μ
φ
^
Δ
E
1
^
M
E
1
^
Δ
E
2
^
M
E
2
^
Δ
E
∗
^
M
E
∗
^
Δ
φ
^
M
φ
^
1.70%
12.65%
7.06%
15.71%
2.77%
12.94%
1.02°
3.14°
8.10%
9.86%
10.86%
13.38%
8.15%
10.42%
1.48°
2.95°


## 5. Conclusions

Complex modulus test results obtained for a straight-run bitumen and a bituminous mixture (respectively,G∗ and E∗) were simulated using 2S2P1D model by 14 users with different levels of expertise of the model. All users worked independently and for the same time duration of one hour to fit the model on both sets of data, providing their own values for all the model constants (including WLF constants at the imposed reference temperatures, 10°C for the bitumen and 9.5°C for the mixture). The following conclusions can be drawn from the analysis of relative and absolute differences between experimental data of the LVE behaviour of both materials and 2S2P1D simulations obtained by all the users:(i)
2S2P1D simulations of LVE behaviour of both the bitumen and the bituminous mixture obtained by the users of the calibration panel satisfactorily approximate experimental data of complex modulus (G∗ for the bitumen and E∗ for the mixture) over the whole range of frequencies and temperatures. The accuracy of the model to fit experimental data is all the more evident if the great spans of equivalent frequencies (between 10−6 and 1010 Hz at the reference temperatures of 10°C for the bitumen and 9.5°C for the mixture) and complex modulus (G∗ of the bitumen comprised between 10−2 and 103 MPa; E∗ of the mixture comprised between 10 and 40000 MPa) are taken into consideration.(ii)
Very low values of COV were found for glassy moduliG0 (1.62%) and E0 (0.57%), for the bitumen and the mixture, respectively. Generally, the variability (expressed by the COV) of constants E00 (for the mixture), δ, k, h, and log τ (for both bitumen and mixture) is approximately equal or lower than 20%.(iii)
Acceptable errors were generally observed over the whole range of frequencies and temperatures for the simulations obtained for both bitumen and mixture by all the users of the panel, regardless of their familiarity and experience with the model. The three most experienced users of the panel obtained particularly accurate approximations of experimental data, confirming the ability of the model to simulate LVE behaviour of both bitumens and mixtures.



---

*Source: 6547025-2019-05-26.xml*