4.1. Construction of Deduction Model
Equation (1) may be transformed into Equation (7). Parameters
and
, representing mixture gradation, asphalt volume, and porosity, are delineated in Equations (8) and (9). In the prediction of the dynamic modulus of an operational asphalt pavement, these parameters are presumed to be constant, and the data are utilized to establish the fitting parameters.
The in-service asphalt pavement modulus
has the following relationship to the original pavement modulus
.
Equation (7) can be transformed into Equation (11).
When utilizing FWD data for calculating the dynamic modulus of in-service asphalt pavement, the relevant parameters are substituted into Equation (11). Equation (11) may be re-expressed as Equation (12).
In this context,
refers to the asphalt viscosity, which is related to the internal temperature of the pavement during the FWD test.
signifies the frequency at which the FWD test is administered. Based on the outcomes of Loulizi’s investigations, the FWD test frequency is consistently maintained at 33 Hz within the pavement structure [
25]. Consequently, this constancy permits the simplification of Equation (12) to Equation (13).
Equation (13) can be further converted to Equation (14).
Upon substitution of Equation (14) into Equation (11), Equation (15) is consequently obtained.
Dynamic modulus data at the same point, different temperatures, and different frequencies are formed in the dynamic modulus database constructed above, so
and
can use dynamic modulus data of the same point and different parameters.
can be calculated by composite viscosity directly, in which the temperature corresponding to composite viscosity is the internal temperature of the pavement structure corresponding to the data
.
can be calculated by composite viscosity directly, in which the temperature corresponding to the composite viscosity is the internal temperature of the pavement structure corresponding to the data
. Substitute the data into Equation (15), the input parameters are shown in
Table 6, and the fitting results are shown in
Table 7,
Figure 5, and Equation (16).
Dynamic modulus data previously established for various temperatures and frequencies at the specified test point enable the direct acquisition of parameters
and
. Since
and
are data associated with temperature and frequency, respectively, the number of fitting data can be effectively expanded by utilizing the correspondence of data at different temperatures and frequencies when constructing the fitting data. This can effectively improve the fitting effect. Parameters
and
, indicative of the composite viscosity, can be ascertained post-asphalt extraction from operational asphalt pavements, if it is available. Within this study, composite viscosity values derived from Equation (5) are utilized for parameters
and
, representing the composite viscosity at the internal temperature of the pavement structure aligned with the
data, while corresponding to the composite viscosity at the internal temperature in line with the
data. Upon the substitution of the data into Equation (15), the input parameters are presented in
Table 6, and the fitting results are depicted in
Table 7, accompanied by
Figure 5. The examination of
Table 7 and
Figure 5 indicates that the extrapolated model has an R-squared (R
2) value of 0.76. A small number of points exhibit deviation from the empirical data.
The deduction model is obtained after the deduction of the basic model. Although the basic model has been verified by a large number of data, it is also an empirical model obtained through a large number of data fitting, and it will have certain biases. In addition, in the fitting data, there is a certain deviation between the FWD modulus back-calculation results and the indoor dynamic modulus test results, and the viscosity data of asphalt is calculated through combination, so the overall data have a certain error. Although the deduction model has a certain reliability, there is a certain deviation from the basic model. Moreover, there are some errors in the fitting data, which leads to some deviation in the obtained dynamic modulus prediction modulus.
4.2. Construction of Gene Expression Programming Model
An exploratory analysis was conducted on the GEP model using the same dataset, with 80% allocated for training and the remaining 20% utilized for validation purposes. The parameter configurations are detailed in
Table 8. The quantity of genes correlates with the complexity of the predictive model. Within the constraints of limited data, an increase in gene count typically yields higher predictive accuracy, albeit at the cost of increased model complexity [
39]. To enhance the model’s applicability and streamline its structure, the gene number parameter was established at 2. Both the population size and gene head length serve as complexity coefficients for individual genes [
40,
41]. Augmented values for population size and gene head length parameters elevate computational complexity and predictive model accuracy while concurrently extending the duration of computation. Based on preliminary computational trials, the population size was configured at 30, and the gene head length was determined to be 12. The derived model is presented in Equation (17), and the computational outcomes for both the training and validation sets are depicted in
Figure 6 and
Table 9.
The inspection of
Figure 6 and
Table 9 reveals that the GEP model exhibits superior predictive capabilities for the dynamic modulus of in-service asphalt pavement. The R
2 for the GEP model stands at 0.88 in the training set and 0.83 in the validation set.
Figure 7 illustrates the deviations between the predicted and actual values of the GEP model within the validation dataset. A majority of the deviations are below 5%, with only 7.1% of the observations exhibiting a deviation exceeding 5%. The maximum deviation is 18.8% and the average deviation is 2.3%. These values suggest that the predictive model, which utilizes gene expression analysis, exhibits a high degree of accuracy. This is mainly due to the reliability of the GEP algorithm, which greatly improves the effectiveness of model construction. In addition, different from the deduction of the model itself, GEP fully utilizes the fitting data to construct the dynamic modulus prediction model, so that a better fitting effect can be obtained.
Figure 8 demonstrates the predictive performance of both the deduction model and the GEP model. The R
2 of the GEP model attained 0.88, in contrast to the deduction model which registered a lower R
2 of 0.76. Within the validation dataset, 92.9% of the data points for the GEP model maintained a deviation below 5%. This indicates that the GEP model exhibits a high level of fit. Utilizing the GEP model allows for the precise calculation of the dynamic modulus of in-service pavement across varying temperatures and frequencies.