Abstract

This paper proposes using the master-curve band (MCB) to incorporate the unavoidable measurement errors and modeling uncertainty into the bitumen master-curve construction. In general, the rheological property of bitumen within the linear viscoelastic region is characterized by the master curve of modulus and/or phase angle, provided that the bitumen complies with the time-temperature superposition principle (TTSP). However, the master-curve construction is essentially a mathematical fitting process regardless of whether or not the original data is perfect enough to fit. For this reason, the MCB was introduced to consider the uncertainty information instead of a single master curve. Rheological data of four kinds of bitumen including unaged and aged bitumen were used to construct the MCBs. The results indicated that the generalized sigmoidal model showed the widest master-curve band, followed by Christensen-Anderson-Marasteanu (CAM) and CAM () models. The width of MCB was a useful tool to identify the sensitivity of bitumen to rheological models. The sensitivity of bitumen to rheological models is associated with the number of active parameters in rheological models and model parameters’ confidence intervals. The construction of an MCB was beneficial to select the rheological models. Accordingly, the CAM () model is proved to be the best to analyze the aging effects.

1. Introduction

Most of the distresses that occur in asphaltic pavements are associated with the binder phase in asphalt mixtures, such as adhesion failure (moisture damage) [1, 2], cohesion failure within bitumen (cracking) [3], fatigue failure [4, 5], and permanent deformation (rutting) [6]. In the characterization of bitumen’s physical properties, rheological behavior was stressed due to its close connection to the above distresses [7]. The terminology rheology was conceptualized as the study of deformation and flow properties of materials [8]. Bituminous material is manifested as being viscoelastic, showing typical time- and temperature-dependent properties [9]. By far, the rheological properties of bitumen are conventionally represented by two primary parameters, which are complex shear modulus (G) and phase angle (δ) measured through the dynamic shear rheometer (DSR).

It is worth noting that the complex shear modulus and phase angle of bitumen only make sense with a given temperature and loading frequency. However, the property of bitumen under a single condition of temperature and loading frequency is not sufficient to reveal the entire performance of bitumen in a broad range of temperature and frequency. For this reason, Airey [8] proposed the black diagram to identify the rheological performance without considering the temperature and frequency. To date, the black diagram has been successfully used in many studies of bitumen rheology [7, 10, 11]. However, the black diagram is not able to reveal the temperature or frequency sensitivity of bitumen. Alternatively, the master curve of modulus or phase angle was proposed to handle the rheological data [1214]. The master-curve construction introduced a time-temperature shift factor to bring rheological data measured at different temperatures and frequencies to one overall continuous curve. The preparatory work required for constructing the master curve is to determine the shift factor as a temperature function. Many approaches were proposed to determine the temperature dependence of shift factors [15], such as a manual procedure [16], the Arrhenius equation [17], the Williams-Landel-Ferry (WLF) formula [18], logarithmic-linear relationship with temperature [19], the Fox method [20], and viscosity-temperature susceptibility (VTS) method [21].

In addition to the determination of shift factors, rheological models are usually used to achieve a perfectly smooth master curve. Yusoff et al. [22] classified rheological models into three categories: nonlinear multivariable models, empirical algebraic equations, and mechanical element models. The nonlinear multivariable models are not usually used as they are empirical models presented in the form of a figure. The empirical algebraic equations, such as the sigmoidal model and Christensen-Anderson-Marasteanu (CAM) model, and mechanical element models, such as the 2S2P1D model, were extensively adopted in many kinds of research. The measured rheological data is mathematically fitted through an optimization process, converting into the master curve at the reference temperature. Herein, it is noticeable that the achieved mater curve, to some extent, mathematically approximates the experimental data. However, uncertainty in the above process could cause significant variation between experimental and predicted data. Figure 1 illustrates the potential uncertainty factors that could give rise to variation between measurement and prediction.

Among these potential uncertainty factors shown in Figure 1, sample variation can be eliminated by measuring multiple parallel samples. As for the testing protocols, some efforts have been made toward accurate measurement for decades. Carswell et al. [23] discussed testing protocols that might cause measurement variation: temperature control, sample geometry, sample preparation, equipment calibration, linear range, and so forth. In their discussion, temperature control and sample geometry are two critical sources of errors. With the development of the device and correlation method, temperature control is no longer an issue if the standard procedure is followed. The geometry size effect can be eliminated by using appropriate testing plates [8]. The bitumen sample showed high modulus in some conditions, and the compliance error could induce inaccurate measurement. Correspondingly, Wu et al. [24] proposed an efficient method to correct rheological data’s compliance error using the 2S2P1D model.

Nonetheless, the measurement uncertainty caused by laboratory or technician difference cannot be avoided (referring to measurement error in Figure 1). A study launched by RILEM investigated the repeatability of binder fatigue tests [25]. This study invited varying laboratories to measure the same type of bituminous materials. In this report, the coefficient of variation measured for complex modulus was up to 15.6%. It is illustrated that the measurement error of rheological properties is practically unavoidable.

On the other hand, in the practical realization of the master curve, the determination of shift factors is flexible toward a goodness of fit of rheological models. Therefore, the variation of shift factors usually depends on the rheological model selection. In this case, unavoidable uncertainty mainly rests on the measurement error, model variation, and fitting variation.

Based on the above analysis, finding an approach that incorporates uncertainty information into the master curve could benefit a more accurate estimation of bitumen rheological properties. For this reason, this study proposed a method to consider the uncertainty that exists in master-curve construction by extending a master curve to a master-curve band (MCB). The MCB has the following potential advantages compared with the master curve: (1) the uncertainty information is included in the master-curve band; (2) the sensitivity of bitumen to the rheological model can be quantitatively featured; and (3) the MCB can help select appropriate rheological models to study the aging or modification effect.

2. The Development of the Master-Curve Band

2.1. Data Variation Measured by DSR

In the first phase of the RILEM Round Robin Binder Fatigue Test mentioned above, it was found that the coefficient of variation (COV) measured for complex modulus was 15.6%, and, for phase angle, the maximum coefficient of variation observed was 3.7% [25]. After improving the test protocols, in the second phase of their study, the variation coefficient changed to 6.7% and 4.8% for complex modulus and phase angle, respectively [25]. Although this study was carried out to identify reasonable fatigue test methods, the results indicated that DSR measurement showed unavoidable measurement errors for complex modulus and phase angle. Besides, the measurement error of the complex modulus is much more significant than that of the phase angle. Unfortunately, a single master-curve cannot include the data variation in its current presentation form.

2.2. Model Variation Caused by Model Types and Fitting Process

Many linear viscoelastic rheological models have been developed to construct the master curve of bitumen. Empirical equations such as CAM and mechanical models such as the 2S2P1D model have been successfully applied to characterize the bitumen viscoelastic properties. Nur et al. compared several mathematical models applied to the unaged, aged, and modified bitumen binders [26]. The results showed that the generalized sigmoidal model performed the best in terms of the goodness of fit, while the CAM model was reported to be senior to other models. Nonetheless, modeling sensitivity has not been thoroughly investigated. One rheological model might be powerful enough to satisfy measurement data. However, it cannot exclude the possibility of overfitting. Therefore, a good quality model should show the goodness of fit and efficiently identify the variation caused by data variation and modeling uncertainty.

2.3. Quantitative Description of Data Variation and Modeling Uncertainty

This section introduced how to consider the data variation and modeling uncertainty quantitatively. As previously noted, the unavoidable data variation could reach up to 15.6%. Therefore, the measured data was extended by multiple corresponding coefficients of variation, as shown in Figure 2(a). Seven data groups, including the measured one, were calculated.

Seven data groups were then subjected to the construction of a master curve using a specific rheological model. In this case, seven groups of model parameters can be obtained. It is time-consuming to conclude the distribution of model parameters experimentally. Alternatively, in previous literature, it was reported that measurements of asphalt mixtures reasonably approximate a normal distribution curve [27]. Therefore, it was assumed that the distribution of model parameters caused by modeling uncertainty also complies with the normal distribution, as shown in Figure 2(b).

Moreover, this study considered three confidence levels (95%, 90%, and 85%) to study the model sensitivity of bitumen. It should be stressed that this study only focuses on the development of a master-curve band for complex modulus. However, the master-curve band’s construction for phase angle can also be achieved following the identical procedure proposed in this study.

2.4. The Mathematical Expression of a Master-Curve Band

For example, the rheological model, namely, the generalized sigmoidal model and the Williams-Landel-Ferry (WLF) formula, was used to construct the master curve. The WLF equation is as follows [15]:where C1 and C2 are two constants relative to material properties.

The generalized sigmoidal model was proposed to describe the stiffness of asphalt mixtures. It has been proven applicable to the bituminous binder and other materials, as shown in the following equation [28]:where G is the complex modulus and f is the reduced frequency. δ, α, β, λ, and γ are model parameters.

According to the time-temperature superposition principle, the complex modulus measured at temperature T and frequency f can be obtained at reference temperature Tref. The reduced frequency is determined through the shift factor, as expressed in the following equation [29]:

The construction of the master curve was performed with the aid of the Solver function in Excel. Combining equation 1 with equation 3, the master curve can be derived by optimizing the parameters with nonlinear least-square regression methods. In particular, this method minimized the sum of square errors (SSE) between measured data and predicted data, which can be written as follows [30]:

The model parameters were achieved by solving equation 4, consisting of a parameter vector x = (x1, x2, ⋯, xn). It is hypothesized that individual parameters are independent of each other. As introduced previously, each parameter should not be a specific value but a range, as shown in equation 5. The normal distribution fitting can determine the lower and upper values for a single parameter.

Based on equation 5, the parameter vector x can be extended to data space U, as expressed in the following equation:

Accordingly, the single master curve is extended to a master-curve band. Using a generalized sigmoidal model, the master-curve band can be induced as the two following equations:

2.5. Procedure to Construct the Master-Curve Band

This section explained how to draw a master-curve band step by step, as shown in Figure 3. After the rheological test, reference temperature was firstly determined. In what follows, the shift factors are determined based on the WLF or Arrhenius equations. The rheological model’s selection is subjective, which could be any empirical or mechanical models. Once the rheological model is determined, model parameters can be obtained by solving the target minimization function. The measured isochronal plots are multiplied by the coefficients from 85% to 115% for the master-curve band construction. In total, six other groups of isochronal plots are calculated.

The normal distribution fit is then performed separately at 85%, 90%, and 95% confidence levels. The data space U in equation 6 for each parameter is obtained at a specific confidence level. According to equation 7, the master-curve band is generated by performing a calculating experiment using MATLAB code. All possible values in data space U are considered. Figure 4 illustrates an example of extending from master curve to master-curve band using a generalized sigmoidal model.

2.6. Data Collection and Selected Rheological Models

Four kinds of bitumen were subjected to DSR measurements to obtain rheological data. Two of them were unaged and unmodified bitumen, marked as U15 and U70. The numbers following “U” indicated the penetration grade. One bitumen was the aged U15, noted as A15. The other one was the SBS modified bitumen based on U70, coded as SBS. All the samples were subjected to the DSR frequency sweep following the standard procedure.

Two rheological models were used to construct the master-curve band: generalized sigmoidal model (introduced above, GS) and CAM model. According to a previous study, bitumen’s theoretical glassy modulus is supposed to be 1 GPa (shear test). Therefore, the CAM model was utilized in two ways. One refers to the CAM model with the glassy modulus fixed at 1 GPa. The other one with the glassy modulus unfixed was noted as CAM () model.

The CAM model is mathematically represented as follows [31]:where represents the glassy modulus, ωc is the crossover frequency, and υ and are model parameters.

3. Basic Properties of the Master-Curve Band

The master-curve band is conceptualized by incorporating the measurement and/or modeling uncertainty into the master curve. However, very few previous studies introduced the master-curve band as well as the basic properties of the master-curve band. Besides, the rheological model’s effect on the master-curve band is still unknown. The basic properties varied with different confidence levels and the bitumen types are vague. Therefore, this section discusses the basic properties of the master-curve band regarding the effect of rheological models, bitumen types, confidence intervals, and bitumen aging, respectively.

3.1. Effect of Rheological Models

Figure 5 presents the master-curve bands of bitumen U15 constructed using three different models at a 95% confidence level. It can be seen that, despite using the same data, different models can result in different master-curve bands. The generalized sigmoidal model showed the widest band, followed by the CAM and CAM () models. It indicates that the CAM () model performed the best in modeling U15 as the least deviation caused by modeling. The only difference between CAM and CAM () models rests on whether the glassy modulus was fixed. However, the CAM model, which has one parameter less than the CAM () model, showed a wider band. It indicates that the fix of the glassy modulus would change the values of other parameters. In this sense, the CAM model and CAM () model are essentially two different models. Table 1 presents the parameter intervals calculated at a confidence level of 95%. It can be seen that the GS model owns five model parameters, which contributed to its wider master-curve band. Although the CAM model has only three model parameters, the interval of was from 17.61 to 28.32, while in the CAM () model the value was only from 1.23 to 1.42. Therefore, it is not true that the band width solely depends on the number of model parameters. Parameter interval is another critical factor that influences the width of the master-curve band.

On the other hand, it can be concluded that and are two model parameters that will not significantly change due to data variation. In this study, parameters with such kind of character are called inert parameters. In contrast, parameters with a considerable change due to data variation are called active parameters. The inert and active parameters can be distinguished by calculating the variation coefficient with a cut-off value of 20%. In this case, the numbers of active parameters in the GS, CAM (), and CAM models are four, two, and one, respectively.

3.2. Effect of Bitumen Types

Figure 6 shows the master-curve bands of different bitumen binders using the GS model at a 95% confidence level. Various bitumen binders responded differently to the identical construction procedure of master-curve bands. As shown in Figure 6, the MCB of U70 showed the widest band, followed by U15 and SBS. The band width reflects the model sensitivity of bitumen. Herein, modified bitumen investigated in this study showed the smallest sensitivity to the data variation. Provided that an identical measurement variation occurred, the master curve of modified bitumen could present the smallest deviation to the correct master curve. Similarly, hard and unmodified bitumen showed less deviation compared to softer and unmodified bitumen in this study.

To further explain the sensitivity of bitumen to rheological models, the confidence intervals for different bitumen are presented in Table 2. For U15, as the confidence level increased, the interval was enlarged except for the parameter γ. It indicated that the parameter γ is relatively stable when using the GS model to construct U15 even in the presence of data variation. However, for other parameters, data variation could result in a change of parameter value.

In terms of U70, the model parameters can be divided into two groups. The first group (active group), including δ, α, and λ, complied with the tendency found in U15. The other group (including β and γ) is the inert group, in which the interval remained almost unchanged as the confidence level increased. In this case, two model parameters would not be significantly altered due to the data variation.

Nevertheless, the variation of U70 between the lower and upper values is relatively more considerable than that of U15, which is ascribed to a wider master-curve band. In the case of SBS, the parameter λ was also classified into the inert group. However, the parameter interval in the active group was the smallest compared to that of U70 or U15.

It can be concluded that the model sensitivity of bitumen is correlated with not only the type of rheological model but also the type of bitumen.

3.3. Effect of Confidence Levels

Figure 7 presents the MCBs of U15 constructed at different confidence levels. When the rheological model was the GS model, the increase of confidence level increased MCB width. As the confidence level increased, the confidence interval of model parameters correspondingly grew. However, in the CAM and CAM () models, MCB width change with the increase of confidence level was not significant. It illustrated that although the confidence level was altered, the constructed master-curve bands were not sensitive to the change of confidence level. In this condition, given that a modeling error existed, the utilization of the CAM and CAM () models would not significantly change the final constructed master curve. This modeling stability, to some extent, can alleviate the modeling uncertainty of the master curve.

Nevertheless, this finding brings about some other issues when comparing master curves between similar bitumen types. Assuming that bitumen was aged insignificantly, the master curve of aged bitumen was similar to its origin bitumen. In this condition, the master-curve comparison cannot reveal the essential variation caused by aging, considering the modeling uncertainty.

3.4. Effect of Bitumen Aging

Figure 8 compares the master-curve bands of unaged and aged bitumen (U15 and A15) using different rheological models. When the rheological model was selected as the generalized sigmoidal model, the master-curve band of A15 overlapped with U15, especially at high frequencies. The master-curve band refers to all possible constructed master curves if modeling uncertainty existed. Therefore, the overlapped area claimed the possibility that the master curves of A15 and U15 were close to each other. In this case, the effect of aging on rheological properties of bitumen might be enervated. Therefore, it is not appropriate for researchers to identify the aging effect using a generalized sigmoidal model to model uncertainty.

As can be seen from Figures 8(b) and 8(c), the overlapping area of master-curve bands for A15 and U15 can be considerably reduced using the CAM or CAM () models. It can be found that the overlapping area was relevant to the master-curve band’s width. As aforementioned, bitumen showed the smallest sensitivity to the CAM () model. Therefore, the width of the master-curve band tended to present in a small value. A small value of band width indicated a small deviation caused by the modeling uncertainty. Hence, it can be concluded that the modeling uncertainty can be eliminated as much as possible if an appropriate rheological model is used. In previous studies of rheological modeling of bitumen binders, only the goodness of fit was concerned. The introduction of master-curve bands revealed that the rheological model’s selection should consider satisfying fitting results and the resistance ability to the modeling uncertainty.

4. Conclusion and Outlook

4.1. Conclusion

In the engineering field, measurement and modeling uncertainty represent a critical topic. However, there is a lack of studies considering the uncertainty in the construction of the master curve. For this reason, this study proposed a method to conclude the uncertainty property. The fundamental technology includes two parts, which are separately relative to the measurement error and modeling uncertainty. For the measurement error, the coefficient of variation was used to describe the error quantitatively. The normal distribution fitting was adopted to obtain the parameter intervals at a specific confidence level for the modeling uncertainty. Consequently, the master-curve band can be generated through a calculating experiment via MATLAB. Through investigating the fundamental properties of the master-curve band, some critical findings were concluded as follows:(1)Ignorance of modeling uncertainty in the estimation of aging or modification would possibly hinder an appropriate conclusion to be achieved, especially when the aging or modification degree is relatively small.(2)Master-curve bands constructed for the same bitumen were different on account of various rheological models. The width of the master band is dominated by the sensitivity of bitumen to the rheological model.(3)The sensitivity of bitumen to rheological models was determined by two factors. The first is the number of active parameters in the selected rheological model. The second is related to the confidence interval of active parameters.(4)In this study, the generalized sigmoidal model showed the widest master-curve band, while the CAM () model performed the best. Therefore, it is recommended to use the CAM () model rather than the generalized sigmoidal model to analyze the aging effect.

4.2. Outlook

It is undeniable that this proposed method could not be an efficient approach to characterize the master-curve construction’s uncertainty information. This study’s primary purpose is to introduce uncertainty into the master curve to strengthen the reliability of the master curve. It is believed that much work should be done to modify the master-curve band toward an efficient and thoughtful approach to include the uncertainty information in the master curve. From the evaluation of the authors, the master-curve band can be modified from the following perspectives:(1)Notably, a detailed standard procedure to construct the master-curve band was not recommended. The reason is that limited types of bitumen and rheological models were investigated in this study. Therefore, an investigation with more samples and models is expected to validate some arguments in constructing the master-curve band. In addition, more potential influence factors of the master-curve band will need to be discussed.(2)The coefficient of variation used in this study was based on the RILEM research. However, it will be more rational when the coefficient of variation is determined by repeatedly conducting the DSR measurements.(3)Normal distribution fitting was adopted to feature the modeling uncertainty, and all the parameters were assumed to be independent of each other. However, the data space of model parameters can be established by other advanced uncertainty analysis methods.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to acknowledge the National Natural Science Foundation of China (no. 52078190) and the German Research Foundation (O.E. 514/10-1). The authors also gratefully acknowledge financial support from the China Scholarship Council (CSC no. 201706710009).