Next Article in Journal
Model and Simulation of Benefit Distribution of Collaborative Cooperation in the Supply Chain of General Contracting Projects
Previous Article in Journal
Factors Hindering Solar Photovoltaic System Implementation in Buildings and Infrastructure Projects: Analysis through a Multiple Linear Regression Model and Rule-Based Decision Support System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Ensemble Prediction Method for Reclaimed Asphalt Pavement (RAP) Mixtures Containing Different Constituents

1
Faculty of Civil Engineering, Sahand University of Technology, Tabriz 51335-1996, Iran
2
Department of Civil and Environmental Engineering, Tarbiat Modares University, Tehran 14115-111, Iran
3
Department of Civil, Geological, and Mining Engineering, Polytechnique Montréal, Montreal, QC H3T 1J4, Canada
4
Department of Civil Engineering, University of Science and Culture, Tehran 14619-68151, Iran
*
Author to whom correspondence should be addressed.
Buildings 2023, 13(7), 1787; https://doi.org/10.3390/buildings13071787
Submission received: 8 June 2023 / Revised: 30 June 2023 / Accepted: 8 July 2023 / Published: 13 July 2023
(This article belongs to the Section Building Materials, and Repair & Renovation)

Abstract

:
Fatigue and rutting are two common damage types in asphalt pavements. Reclaimed asphalt pavement (RAP), as a sustainable approach in the pavement industry, deals with the foregoing damage. Fatigue and rutting characteristics of asphalt pavement are generally assessed using laboratory tests, taking a long time and consuming significant amounts of raw material. This study aims to propose a novel approach for predicting fatigue and rutting performance of RAP mixtures. A new ensemble prediction method, named COA-KNN, is introduced by combining the coyote optimization algorithm and K-nearest neighbor to increase the accuracy of fatigue and rutting prediction. In order to evaluate the accuracy, the proposed method was compared against robust prediction methods, including random forest (RF), gradient boosting (GB), decision tree regression (DT), and multiple linear regression (MLR). Afterward, the influence of each variable on the mentioned damages is examined, and the variables are ranked based on their relative influence on the mentioned damages. The results suggest that COA-KNN outperformed other prediction techniques when comparing different performance indicators. Total binder content in asphalt mixes and the PG span of the virgin binder added to the recycled asphalt mixture had the highest relative influence on fatigue and rutting performance, respectively.

1. Introduction

Transportation is a critical system for every communication as most human and goods transmissions are conducted using infrastructure [1,2]. As an enormous consumer of natural resources, the transportation industry contributes 25% and 22% of global energy usage and fossil fuels burned across the world, respectively. Accordingly, 22.7% of global greenhouse gas (GHG) production is generated in this industry. A total of 7% of the mentioned GHG is emitted by pavement which is responsible for transferring 60% of freight and 80% of passenger traffic [3]. Furthermore, the pavement industry consumes a significant amount of energy. For instance, aggregate and binder production and hot mix asphalt (HMA) operation to produce asphalt consume 54, 5810, and 275 MJ/ton energy, respectively [4]. Thus, the pavement sector can be considered a harmful industry to the environment, which is a significant concern.
The concerning environmental impacts and the high cost of crude oil have caused the pavement industry to develop a method to recycle waste pavement materials [5]. Consequently, many methods have been suggested as sustainable approaches, including the application of recycled materials such as reclaimed asphalt pavement (RAP) [6,7] or warm and cold mix asphalt technologies [8]. Warm and cold mix technologies still face problems, such as moisture damage potential and durability due to water in the binder emulsion, long curing time, and uncertainty of long-term performance [9]. Consequently, among these sustainable options, asphalt recycling is the most common method currently used in this industry, as it creates a cycle of reusing natural resources and reduces the need for virgin binder and material [10]. However, RAP application provides the mentioned environmental benefits, increasing its content in asphalt mixtures, and the effect of the degree of blending between virgin and RAP binders may deteriorate some characteristics of asphalt pavements [11], such as cracking resistance, moisture damage, etc. Thus, there has been a dire need to come up with new solutions to increase RAP content without the mentioned problems [11].
To this end, laboratory tests are generally executed to evaluate the mentioned mechanical behavior of asphalt mixes as a conventional technique. Two major damage types of concern in asphalt pavements are fatigue cracking and rutting. Thus, previous studies investigated these failures in mixtures containing RAP using laboratory tests. In order to evaluate fatigue cracking, the resilient modulus, fracture energy, dynamic modulus, etc., have been investigated using various tests such as four-point bending beam, dynamic modulus, semi-circular bending, and dissipated creep strain energy tests [12,13]. Similarly, regarding the rutting resistance of asphalt mixtures, rut depth is a parameter that is used for evaluating this performance. Accordingly, rutting has been investigated in various studies with tests including wheel tracking and asphalt pavement analyzer tests [14,15,16,17,18,19]. These laboratory tests generally showed that increasing the RAP content in mixtures increased the stiffness of mixtures and improved their rutting resistance but clearly did not specify the effects of the mixture components on the behavior of RAP mixes.
All the previous laboratory studies experimentally evaluated the performance of asphalt mixtures containing RAP, though the optimal mixtures are selected between the produced samples. That is, a limited number of specimens have been generally examined in experimental laboratory experiments. Moreover, laboratory testing takes a long time to prepare the specimens, equilibrate the test condition, and perform the tests. Additionally, these complex procedures take expert technicians, expensive equipment, and even more energy and material. Furthermore, laboratory specimens harm the environment because they should be stored in landfills after performing the tests. To diminution these issues, researchers began to apply soft computing (i.e., machine learning and optimization) techniques to predict various asphalt damages [20]. As such, Garbowski and Pożarycki [21] used an optimization framework to determine the thickness and stiffness of the pavement layers through a multi-level inverse approach. Sarkhani Benemaran et al. [22] applied extreme gradient boosting to predict the resilient modulus of flexible pavement foundations. Esmaeili-Falak and Sarkhani Benemaran [23] developed a new prediction model using extreme gradient boosting to predict the resilient modulus of modified base materials.
Ghasemi et al. [24] applied artificial neural networks (ANN) and linear regression to predict the dynamic modulus of asphalt, considering volumetric and particle size gradation as input variables. Only nine mixture proportions were taken as the dataset. The results suggested that ANN outperforms linear regression when comparing the prediction accuracy. Behnood and Daneshvar [20] estimated the dynamic modulus of asphalt mixes using the M5P algorithm and a dataset including 4022 asphalt mixture samples. The dataset included different input variables, including aggregate gradation, volumetric properties, binder characteristics, etc., and one output variable (i.e., dynamic modulus). The mentioned model obtained an R2 of 0.919. Behnood and Golafshani [25], in another study with the same dataset and variables, used the biogeography-based optimization (BBO) algorithm, which enhanced the accuracy with an R2 of 0.9601.
Likewise, prediction techniques have been applied to predict the rutting of asphalt pavements. Initially, regression models were used to predict the rutting of asphalt pavements [26]. However, recently, researchers have started to apply powerful prediction techniques to predict rutting. As such, Ullah and Zainab [27] studied the rutting performance with a dataset including 0 to 60% RAP. The input variables were loading cycles, RAP percentage, RAP binder content, etc. Using an ANN-based model, they predicted the rutting behavior with a testing data R2 of 0.997. Majidifard et al. [28] predicted the rut depth of asphalt mixes using genetic expression programming. The dataset consisted of 96 test samples. The used parameters were mix and binder characteristics, aggregate gradation, and the content of additives in the mixtures, which led to an R2 of 0.84.
Previous studies mainly focused on dynamic modulus and rut depth as indices for fatigue and rutting resistance, respectively, which are the mechanical properties of asphalt mixtures. Since the aim of this study is the evaluation of fatigue and rutting characteristics of these mixes in the field, two general indexes (fatigue life and rutting resistance index) are chosen for the prediction procedure. In addition, the used algorithm is a key element of the prediction procedure, and hence, an accurate and robust algorithm must be developed to achieve reliable results. Moreover, a major part of the dataset used for the prediction procedure was concentrated on low percentages of RAP. The lack of data for mixes with high percentages of RAP affects the prediction procedure resulting in inaccurate models.
In this study, unlike the discussed studies, two general indices were used for fatigue and rutting resistance of asphalt pavement instead of dynamic modulus and rut depth to better indicate the field behavior of asphalt mixtures. Meanwhile, in this study, the database included mixes containing high percentages of RAP. An appropriate prediction algorithm should be applied to the mentioned prediction problem to reach the highest possible accuracy. In this regard, a new robust ensemble method was developed to predict the fatigue and rutting performance of asphalt mixes accurately. The mentioned method was developed by merging a machine learning technique with a metaheuristic optimization algorithm. The introduced prediction method was then compared against conventional prediction methods in order to evaluate its performance and prioritize the performance of all the prediction algorithms applied in this research on fatigue and rutting prediction problems. Further, the impact of each input variable on the outputs of the models was investigated. That is, the features (variables) were ranked based on their relative influence on fatigue and rutting resistance. Therefore, researchers could select the optimal parameters to improve the fatigue and rutting performance of asphalt mixtures containing high RAP content.

2. Methodology

This study attempts to predict the fatigue and rutting performance of asphalt mixtures containing RAP. To this end, two general indices for fatigue and rutting (i.e., Nf and RRI, which will be further explained in the following sections) are used since these indices better indicate the fatigue and rutting life of asphalt. The Table of abbreviations with their explanations are shown in Appendix A. In order to achieve higher accuracy in the predicted values, a new ensemble method is developed. The performance of the introduced method is further evaluated by comparing its accuracy with random forest, gradient boosting, decision tree, and multiple linear regression. Ultimately, the effect of each input variable (e.g., RAP percentage and total binder percentage) on the fatigue and rutting prediction models were assessed and prioritized. The methodology flowchart is indicated in Figure 1.
The following sub-sections are presented in this section:
  • The process of data preparation
  • The applied algorithms for the prediction procedure
  • The description of machine learning performance indicators used to evaluate the accuracy of prediction algorithms
  • The methods applied to analyze the effects of the model’s inputs on output variables

2.1. The Data Preparation Process

An accurate model needs to contemplate all affecting features on the intended parameters, which are fatigue and rutting indices, in this study. As this study aims at the fatigue and rutting prediction of asphalt mixes containing RAP, two datasets are extracted and generated from authentic international publications. A dataset including 227 asphalt mixture samples is prepared for the fatigue life modeling [29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66]. The model’s input variables (features) include the RAP percentage, the span of the PG of the virgin binder, the intermediate-temperature performance grade (PG) of the virgin binder, the rejuvenator content, the virgin asphalt content, the asphalt content in RAP, the total asphalt content, the percentage of aggregate smaller than 4.75 mm (fine aggregate), the percentage of aggregate larger than 4.75 mm divided by the percentage of aggregate smaller than 4.75 mm (coarse aggregate/fine aggregate ratio), and the dynamic modulus test temperature. It should be noted that since fatigue life (Nf) is very skewed, the natural logarithmic scale of it is preferred to be considered for the model’s output variable since the logarithmic transformation increases the fatigue life prediction accuracy [67]. The fatigue life of the mixes is evaluated to assess the fatigue performance of asphalt pavements. Equation (1) is applied to calculate the fatigue life of the asphalt mixture [68].
N f = 0.0796 ε t 3.291 | E | 0.854
where Nf is the allowable cycles of fatigue life, εt is the tensile strain at the bottom of asphalt layers under the wheel (strain) which has been considered 0.0002 strains [69], and |E*| is the dynamic modulus of asphalt mixture (MPa) extracted from dynamic modulus test (AASHTO TP62).
The description of the fatigue prediction input and output variables and their descriptive statistics are given in Table 1. The dataset contains 12 input variables and one output variable.
Likewise, a dataset including 226 data samples of asphalt mixes containing RAP is employed for rutting resistance performance prediction [14,15,16,19,68,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88]. Different input variables, including the RAP content, the span of PG of the virgin binder, the high-temperature PG of the virgin binder, the rejuvenator content, the virgin binder content, the binder content in RAP, the total binder content, the nominal maximum aggregate size, the percentage of fine aggregate size, and the coarse aggregate to fine aggregate ratio, are taken into consideration as the model’s input features. Meanwhile, the rutting resistance index (RRI) is considered the model’s output variable since RRI implies the rutting resistance of asphalt pavement based on the rut depth and the loading cycles reported from the laboratory tests. The RRI can be calculated from Equation (2) [89].
R R I = N × ( 1 R D )
where RRI is the rutting resistance index, N signifies the number of cycles at the end of the test, and RD is the rut depth (inch) (AASHTO T324). In this equation, the maximum allowable rut depth is 1 in [89].
The description of the input and output features and their descriptive statistics are presented in Table 2 for rutting. The dataset contains ten input variables and one output variable. Further, the distributions of fatigue and rutting indices are shown in Figure 2.
It is noteworthy that 75% of data samples are used as training data. Training data is used to tune hyperparameters and train the model. That is, K-fold cross-validation (considering K = 5) and grid search methods are simultaneously used to tune the hyperparameters of prediction techniques. Then, the other 25% of (testing data) is used to evaluate the performance of prediction techniques on unseen data samples.

2.2. Prediction Algorithms

As described before, in order to increase the accuracy of the fatigue and rutting prediction, a new hybrid algorithm is proposed, and to assess the performance of this newly developed algorithm, four conventional methods are employed. Before running the machine learning algorithms, the input variables are scaled using the standard scaler in the python SKlearn library. Standard scaler applies Equation (3) to scale the range of all variables.
S i = x i u s t d
where, Si implies the scale value of the data sample i, xi is the raw data value of the data sample I, and u and std signify the mean and standard deviation of data samples.
In this section, the utilized algorithms are explained.

2.2.1. Decision Tree

The decision tree (DT) is a prediction technique capable of modeling regression and classification problems. DT simulates the structure of trees, particularly their leaves, and nodes, to divide the dataset into categories with similar characteristics. DT is a straightforward prediction method, and its outcomes are easy to interpret. DT has been widely used for prediction purposes since it can be modeled as a white-box prediction technique. It was demonstrated that DT could face large-scale datasets with a significant number of samples. DT generally works based on some decision rules, determining feature cut-off thresholds. Initially, it is assumed that all data samples constitute a tree. Then the algorithm moves forward and splits samples into different sub-samples using a feature threshold in each leaf node. Subsequently, each sub-sample is divided into smaller sub-samples using another feature threshold. This data sample division is kept on until the model obtains pure sub-samples or termination criteria are met. Afterward, DT employs a back-forward pruning technique to cut nonessential splits (branches) to enhance the model’s efficiency and decrease the computational complexity [90].

2.2.2. Random Forest

Random forest (RF) is an ensemble prediction technique applied for regression, classification, and outlier detection. RF mainly obtains high prediction accuracy, performs well in large-scale datasets, and reduces the likelihood of over-fitting. It was demonstrated that RF is a computationally efficient prediction method and can perform well even if the number of data samples is inadequate. RF generally employs different randomly generated decision trees. Then, the generated models are aggregated using the bagging technique. That is, a powerful ensemble method is created using bootstrap aggregation decision trees [91]. Consequently, decision tree algorithms are run and predict the model’s dependent variable. RF considers the prediction value of all decision trees in the prediction process. Hence, the average value of decision trees’ outputs is considered the RF prediction value [3].

2.2.3. Gradient Boosting

Gradient boosting regression is an ensemble prediction method that can be used for regression and classification problems. Gradient boosting regression (GB) combines various decision trees to create a powerful prediction model. GB generates basic decision trees iteratively and combines these weak regressions through boosting process. Additional decision trees are added in each iteration of the boosting process. Hence, a new regression model is generated by GB at each iteration, and the combination method is optimized using a loss function (e.g., absolute error or squared error) and training data. Accordingly, the model seeks to add basic trees that can reduce the loss function to the prediction model. It has been demonstrated that GB is a powerful method for databases containing unbalanced data and outliers, and it can guarantee prediction efficiency for these databases [92].

2.2.4. Multiple Linear Regression

Multiple linear regression (MLR) is a conventional technique applied for prediction, data representation, and an indication of relevant feature relationships in cases where data follow a linear pattern. MLR is a quick and straightforward method that is easy to interpret since it is a white-box prediction technique. However, the regression structure in MLR should be determined before running the model, and the predefined structure may reduce the model’s flexibility and prediction accuracy [93].

2.2.5. The Proposed Prediction Method (COA-KNN)

As mentioned, a new ensemble regression method is introduced in this study to predict asphalt’s fatigue index and rutting index accurately. The proposed model is generated using the coyote optimization algorithm (COA) as a robust metaheuristic algorithm and K-nearest neighbors regression (KNN) as a powerful prediction technique. COA and KNN are described in this section, and afterward, the proposed ensemble regression method is presented.
COA is a novel swarm intelligence optimization algorithm proposed by Pierezan and Coelho. This algorithm mimics the social behaviors and interactions of coyote Canis latrans. COA assumes that each solution vector is a coyote in the feasible region (coyote society). The coyote’s behavior is modeled to investigate the feasible region and finds the optimal solution. Each coyote is a vector representing the values for each independent variable of the optimization problem. First, some coyotes are generated by allocating random values to all independent variables. Afterward, the fitness values of coyotes are calculated using the problem’s objective function. Consequently, the coyotes are categorized into different herds. That is, solution vectors are divided into various groups to investigate different parts of the feasible region at the same time and avoid sticking to locally optimal solutions. Subsequently, the coyotes are compared based on their fitness value, and the most valuable coyotes (highest fitness value) in each group are called alpha [94].
Then, the transfer culture operation is performed, and all coyotes are attracted by their groupmates and their groups’ alpha. Hence, they are transferred to points, which are located near their alpha and their groupmates. The coyotes with higher fitness values have higher attraction than others. Therefore, more coyotes are accumulated in areas with better previous experiences. Meanwhile, some solution vectors are transferred between groups in each iteration. This process reduces the likelihood of sticking to local optimal solutions by scattering some solutions in the problem’s feasible region. Ultimately, the death and birth operator is performed to enhance the algorithm’s performance. In this operation, the weakest solution vectors are removed from society (death), and new random solution vectors are created (birth) to investigate new areas in the feasible region. These processes are performed in all iterations, and after meeting the termination criteria, the algorithm stops working, and the solution vector with the highest fitness value (best coyote) is reported as the optimization problem’s optimal solution [95].
K-nearest neighborhood regression (KNN) is a robust regression model applied for statistical analysis and prediction since the 1970s. KNN is a non-parametric regression model that predicts the output values considering the resembling data points in the data set. That is to say; all training data points are plotted in an F-dimension space according to the prediction problem’s input values, where F implies the number of input values in the prediction problem. Subsequently, the distance of each data point to all data samples is computed. Then, the output value of K nearest data samples is employed to predict the output value of the mentioned data point. Among K nearest data samples, the data points with lower distances affect the output value more than the data samples with longer distances. However, the model is highly sensitive to the value of K, and the model may not obtain accurate predictions if appropriate values are not considered for K [96]. To this end, the proposed model combines different KNN models with different K values to prevail in the mentioned deficiency.
As previously mentioned, COA and KNN are combined to propose a new ensemble prediction method. In this regard, an optimization problem is modeled, which is indicated in Equations (4) and (5).
min | y i p r e y i exp | i { 1 , 2 , , I }
s . t : y i p r e = c 0 + k = 1 K c k × y i K N N k i { 1 , 2 , , I } , k { 1 , 2 , , K }
where yipre and yiexp signify the prediction value and real value of validation data sample i. yiKNNk denotes the prediction value of validation data sample i predicted by the KNN model k. c0 and ck are the constant value and the coefficient of the KNN model k. That is, c0 and ck are the optimization decision variables that should be optimized using COA. I is the number of validation data, and K is the number of KNN models applied in the proposed ensemble method.
Equation (4) is the objective function of the optimization problem. As can be seen, the introduced model aims to minimize the mean absolute error of validation data by finding the optimal coefficient of each KNN model. Training data are used to run the model. Then, the trained model predicts the response variable (i.e., rutting or fatigue) for the validation data (yiKNNk). Afterward, the real validation data and their predicted values are used to evaluate the mean absolute error. The optimization model aims to minimize the mean absolute error by finding optimal values of decision variables. Equation (5) is the constraint of the optimization problem, in which the prediction value is calculated using the values predicted by different KNN models, coefficients, and the optimal constant value. In this investigation, ten KNN models with different K values from 1 to 10 are modeled. However, it is recommended that the number of KNN is increased in the cases where the number of data points is higher than that of the current study. Afterward, these ten KNN models are run. Consequently, the validation data are employed in the models, and the validation data output value for all KNN models is estimated. Then, the optimal constant value and the optimal coefficients for each KNN model are obtained by solving the optimization problem (using COA). The KNN models with the coefficient of zero are removed from the ensemble model, and the prediction values are computed using the KNN models with their non-zero coefficients and optimal constant value. This process removes the models with inappropriate K values, and only robust KNN models remain in the proposed prediction method. Ultimately, testing data are applied to assess the prediction accuracy of the proposed model and other conventional prediction techniques. Since COA is a metaheuristic optimization algorithm, it was run five times to find the optimal structure of KNN models [97].

2.3. The Evaluation of Algorithms Performance Using the Performance Indicators

The accuracy of prediction algorithms is evaluated using machine learning performance indicators. Consequently, these performance indicators are utilized to compare the performance of prediction algorithms and detect the most accurate model. In this study, coefficient of determination (R2), mean absolute error (MAE), mean absolute percentage error (MAPE), mean squared error (MSE), variance account for (VAF), and A10-index are considered machine learning performance indicators [98]. These indicators are presented in this part.
R2 is the proportion of the variance in the dependent variable that is predictable from the independent variables, calculated using Equation (6) [99].
R 2 = 1 i = 1 m ( X i Y i ) 2 i = 1 m ( Y ¯ Y i ) 2
where m equals the number of values (data samples), Xi equals the predicted ith value, Yi equals the real ith value, and Y ¯ equals the average of the real values.
Mean absolute error (MAE) is an indicator used for measuring the arithmetic average of deviations between predicted and real values. MAE is estimated using Equation (7) [100].
M A E = 1 m i = 1 m | X i Y i |
Mean absolute percentage error (MAPE) implies the mean percentage deviation between predicted and real values. MAPE can be calculated based on Equation (8).
M A P E = 1 m i = 1 m | Y i X i Y i |
Mean squared error (MSE) is the average squares of the differences between predicted and real values. The formula of MSE is shown in Equation (9) [101].
M S E = 1 m i = 1 m ( X i Y i ) 2
Variance accounts for (VAF) and A10-index are also applied to compare the performance of prediction techniques:
V A F = [ 1 var ( Y i X i ) var ( Y i ) ] × 100
A 10 i n d e x = n u m b e r   o f   d a t a   w i t h   v a l u e s   o f   r a t e   a c t u a l / p r e d i c t e d   v a l u e s   ( r a n g i n g   f r o m   0.9   t o   1.1 )   n × 100

2.4. Prioritizing Features on Asphalt Mixture Characteristics

A good practice in modeling a prediction technique is to analyze the sensitivity of potential variables to the prediction model [102]. In ensemble learning techniques, the important weight represents the sensitivity of potential variables to the prediction model. In this study, the effect of each input variable (e.g., RAP content, total binder content, the span of the PG of the virgin binder, etc.) on the fatigue and rutting prediction models is investigated. In this regard, prediction models are run for each dataset. Afterward, the models with the highest accuracy are detected. Subsequently, the average importance weight of the most accurate models is used to prioritize each input feature for each dataset. The ranking of input variables on fatigue and rutting performance can determine the essential parameters in optimizing fatigue and rutting performance of asphalt mixes containing RAP, which is an optimal approach to modify the mixture proportion or mixture conditioning considering the most important input variables.

3. Results and Discussion

As previously mentioned, this study aims to propose a novel approach to accurately predict two major damages (i.e., fatigue and rutting) in asphalt mixtures containing different RAP contents using two general indices. To this end, various robust prediction techniques, including random forest, gradient boosting, decision tree, and multiple linear regression, are employed. Moreover, a new hybrid machine learning method is introduced using a combination of COA and KNN. In this section, the following steps are described, respectively:
  • The performance of prediction algorithms on fatigue and rutting prediction is analyzed.
  • The effects of each variable on fatigue and rutting characteristics are presented, and then the mentioned variables are prioritized based on their effects on the prediction performance.
  • The error histogram of the fatigue and rutting models is presented.

3.1. Performance of Algorithms

As mentioned, four conventional machine learning algorithms are used for the prediction process. Moreover, a new ensemble method, called COA-KNN, is proposed to increase the prediction accuracy. The developed model and the conventional prediction algorithms are then run on each dataset (i.e., fatigue and rutting) individually in order to predict the fatigue and rutting behavior of asphalt mixes containing RAP. In this section, the performance of prediction techniques on fatigue and rutting prediction problems are discussed. The performance of prediction techniques is compared through the machine learning performance indicators.

3.1.1. Performance of Machine Learning Algorithms—Fatigue

The obtained prediction model for fatigue, using multiple linear regression, is shown in Equation (12). It should be noted that the scaled values of the coefficients are given in this equation, so S refers to the parameters’ scaled value.
ln ( F a t i g u e ) = 0.121 × S R A P + 0.072 × S P G S p a n 0.094 × S P G i n t e r + 0.040 × S R e j u v e n a t o r 0.056 × S A C V i r g i n 0.023 × S A C R A P + 0.137 × S A C T o t a l + 0.072 × S N M A S + 0.043 × S A g g . b e l o w # 4 + 0.010 × S A g g . a b o v e # 4 / A g g . b e l o w # 4 + 0.034 × S T e m p e r a t u r e + 17.844
where ln (Fatigue) represents the natural logarithm predicted value of the fatigue index. SRAP, SPGspan, SPGinter, SRejuvenator, SACvirgin, SACRAP, SACTotal, SNMAS, SAgg.below#4, SAgg.above#4/Agg.below#4, and STemperature signify the scaled values of the RAP content, the span of PG of virgin added binder, the intermediate-temperature PG of virgin binder, the rejuvenator content, the virgin binder content, the binder content in RAP, the total binder content, the nominal maximum aggregate size in gradation, the percentage of aggregate smaller than 4.75 mm (fine aggregate), the percentage of aggregate larger than 4.75 mm divided by the percentage of aggregate smaller than 4.75 mm (coarse aggregate/fine aggregate), and the dynamic modulus test temperature, respectively.
In multiple linear regression, the coefficient of each feature indicates the importance and the effect of the feature on the damage. That is, a greater coefficient implies that the feature has more impact on the model’s output (i.e., fatigue). According to Equation (12), since the total binder content coefficient is 0.137, it is the most important feature of fatigue. In addition, the coarse aggregate/fine aggregate, with a coefficient value of 0.010, has the lowest impact on fatigue performance. Based on the results, the most effective feature on fatigue performance of asphalt mixtures containing RAP is the virgin binder content, followed by RAP content, the intermediate-temperature PG of virgin binder, the span of PG of virgin added binder, the nominal maximum aggregate size in gradation, the virgin binder content, the fine aggregate, the rejuvenator content, the dynamic modulus test temperature, the binder content in RAP, and the coarse aggregate/fine aggregate, respectively.
In the COA-KNN algorithm, to achieve the final predicted values, different KNN models with different values for K (the number of neighbors in KNN) are applied simultaneously. The result of each KNN model is multiplied by an optimal coefficient to present the ultimate prediction value. The optimal coefficients of the mentioned ensemble method are determined by an optimization algorithm (i.e., COA).
The predicted values of KNN models for fatigue performance are given in Equation (13). As can be seen, only four KNN models are available in the final equation (Equation (13)), and this is because of allocating a coefficient of zero to other KNN models in the optimization process.
ln ( F a t i g u e ) =   0.079 × K N N ( 5 ) + 0.004 × K N N ( 8 ) + 0.918 × K N N ( 10 )
where KNN (i) signifies the KNN model considering the number of neighbors (K) is equal to i. To predict the fatigue performance, the predicted value by KNN(5) is multiplied by 0.079 and is summed up with the predicted value by KNN(8) multiplied by 0.004, then again is summed up with the predicted value by KNN(10) multiplied by 0.918. As mentioned, the coefficient of other KNN models is zero, i.e., they are not chosen by the model. That is, three KNN models predict fatigue performance, and COA-KNN applies these predicted values to predict the ultimate prediction value.
The machine learning performance indicators of prediction algorithms on the fatigue index prediction problem for testing and training data are indicated in Figure 3. The machine learning performance indicators were used to evaluate the performance of the developed model in comparison to the conventional methods. Regarding Figure 3, the testing data R2 of COA-KNN are 0.07 more than the gradient boosting, 0.10 more than the random forest, 0.26 more than the decision tree, and 0.35 more than multiple linear regression. Thus, these testing data R2 obtained using COA-KNN are more than the other conventional techniques, indicating the proposed method’s higher accuracy. Based on these testing data MAE performance indicators, the lowest MAE of test data belongs to COA-KNN, followed by gradient boosting, random forest, decision tree, and multiple linear regression. The MAEs of COA-KNN test data are 0.01, 0.02, 0.04, and 0.07 ln (cycles), lower than the MAE test data in gradient boosting, random forest, decision tree, and multiple linear regression, respectively. Hence, it can be theorized that COA-KNN outperforms other prediction techniques in terms of MAE. Similarly, the MAPE COA-KNN test data are 0.04% lower than gradient boosting, 0.11% lower than random forest, 0.21% lower than decision tree, and 0.38% lower than multiple linear regression. Thus, it can be postulated that the introduced prediction model (COA-KNN) is more accurate than the other methods regarding MAPE. The same trend can be observed from the testing data MSE. That is, the MSE test data obtained using COA-KNN are 0.01, 0.01, 0.03, and 0.03 (ln (cycles))2 lower than gradient boosting, random forest, decision tree, and multiple linear regression, respectively. The VAF test data of COA-KNN are 5.77%, 5.97%, 14.85%, and 47.47% higher than gradient boosting, random forest, decision tree, and multiple linear regression, in the order given. However, the highest testing data A10-index is reached by decision tree, which is 5.7% higher than COA-KNN.
Considering all performance indicators, it can be deduced that the introduced ensemble regression method (i.e., COA-KNN) outperforms all conventional prediction techniques. In addition, among the conventional models, random forest and gradient boosting have acceptable performances. Nonetheless, linear regression and decision tree regression may not be qualified to predict fatigue damages.

3.1.2. Performance of Machine Learning Algorithms—Rutting

The multiple linear regression prediction equation for rutting performance is presented in Equation (14).
R u t t i n g = 593.19 × S R A P 1431.64 × S P G S p a n + 1007.17 × S P G H i g h 607.97 × S R e j u v e n a t o r + 2031.25 × S A C V i r g i n + 1375.87 × S A C R A P 1704.87 × S A C T o t a l + 197.98 × S N M A S 554.33 × S A g g . b e l o w # 4 2681.74 × S A g g . a b o v e # 4 / A g g . b e l o w # 4 + 12765
where Rutting represents the predicted value of the rutting index, SRAP, SPGspan, SPGhigh, SRejuvenator, SACvirgin, SACRAP, SACTotal, SNMAS, SAgg.below#4, and SAgg.above#4/Agg.below#4 imply the scaled value of the RAP content, the span of PG of virgin added binder, the high-temperature PG of virgin binder, the rejuvenator content, the virgin binder content, the binder content in RAP, the total binder content, the nominal maximum aggregate size in gradation, the percentage of fine aggregate, and the coarse aggregate/fine aggregate value.
Based on Equation (14), it can be observed that the rutting performance is highly affected by the coarse aggregate/fine aggregate value with a coefficient of −2684.74, indicating the corresponding impact is negative. Meanwhile, the NMAS coefficient in aggregate gradation shows that aggregate gradation has the lowest effect on the rutting index, with a value of 197.983. Other effective input features on the rutting performance of asphalt mixtures containing RAP are the virgin binder content, the total binder content, the span of PG of virgin added binder, the binder content in RAP, the high-temperature PG of virgin binder, the rejuvenator content, the RAP content, and the percentage of fine aggregate, respectively.
The KNN prediction model for rutting performance is shown in Equation (15).
R u t t i n g =   0.749 × K N N ( 1 ) + 0.076 × K N N ( 6 ) + 0.174 × K N N ( 7 ) + 51.562
where KNN (i) is the KNN model in which the number of neighbors (K) is equal to i. According to Equation (15), the predicted value by KNN (1) multiplied by 0.749, summed up with 0.076 and 0.174 times the predicted values by KNN (6) and KNN (7), respectively. The obtained value is summed up with the optimal constant (i.e., 51.562), and the ultimate value is the COA-KNN rutting performance prediction.
The machine learning performance indicators of prediction algorithms on the rutting index prediction problem for testing and training data are indicated in Figure 4. As can be seen, the testing data R2 reaches the maximum value using COA-KNN. The testing data R2 of COA-KNN is 0.050 more than random forest, 0.091 more than gradient boosting, 0.222 more than decision tree, and 0.806 more than multiple linear regression. Hence, COA-KNN outweighs other prediction algorithms based on testing data R2. Considering the MAE test data, the MAE test data of COA-KNN are 361.44 inches lower than random forest, 890.44 inches lower than gradient boosting, 897.21 inches lower than decision tree, and 3123.89 inches lower than multiple linear regression. Accordingly, COA-KNN is the most accurate algorithm in terms of error reduction. Similarly, the MAPE test data of COA-KNN are 1.44%, 5.06%, 5.70%, and 22.60% lower than that of random forest, decision tree, gradient boosting, and multiple linear regression, in the order mentioned. Regarding the MSE, the lowest MSE test data are related to COA-KNN. The MSE test data of COA-KNN are significantly lower than conventional methods, which are 1,297,190 inches2 less than random forest, 2,506,610 inches2 less than gradient boosting, 7,206,180 inches2 less than decision tree, and 19,897,280 inches2 less than multiple linear regression. The VAF test data of COA-KNN are 6.6%, 11.72%, 28.97%, and 60.87% higher than random forest, gradient boosting, decision tree, and multiple linear regression. Likewise, the highest A10 index is obtained using COA-KNN, with a value of 98.73%.
Thus, it can be postulated that COA-KNN performs better than conventional prediction techniques. To this end, the introduced ensemble method predicts the rutting damage with the highest accuracy, while the multiple linear regression has the worst performance. Among the other conventional algorithm, the prediction models of random forest and gradient boosting are more accurate.

3.2. Relative Influence of Variables

It is essential to assess the effects of each input variable (e.g., RAP percentage) on the prediction model’s output (i.e., fatigue and rutting) to understand better how asphalt performance can be improved. Since the results show that multiple linear regression is the worst algorithm based on performance indicators, the outcomes of Equations (12) and (14) are not trustable, and a practical approach should be considered to analyze model parameters meticulously. To this end, the importance weight values presented by Random Forest and Gradient Boosting (the most accurate conventional techniques) are used to prioritize models’ features. Hence, the average value of importance weights presented by the mentioned algorithms is employed in order to consider both algorithms. That is, the average importance weights values of the two models are calculated, and subsequently, the features are prioritized based on it. The feature with the largest average important weight value is considered a variable that affects the model’s output the most.

3.2.1. Fatigue Input Feature Performance

Table 3 shows the important weight values of the random forest and gradient boosting for fatigue performance of asphalt mixtures containing RAP. As can be seen from Table 3, the total binder content in asphalt mixes affects the fatigue performance the most since the binder content is responsible for the stiffness of asphalt mixes which is in line with the results presented by Sreedhar and Coleri [103]. The second important feature in fatigue damage is the virgin binder added to the recycled asphalt mixtures due to the fact that the virgin binder lessens the stiffness of the mixture containing stiff binder in RAP; therefore, it increases the fatigue life of the mixture [17]. The intermediate-temperature PG of virgin binder is the next feature that can affect fatigue cracking since fatigue cracking occurs in the intermediate temperature of the asphalt pavement environment [104]. It is worth noting that RAP content is the fourth important feature of fatigue performance of mixtures containing RAP as it makes the mixture stiffer [17]. Other input features affecting fatigue cracking in asphalt pavements are the dynamic modulus test temperature, the asphalt content in RAP, the span of PG of virgin added binder, the percentage of aggregate larger than 4.75 mm, and the coarse aggregate/fine aggregate value, the nominal maximum aggregate size in gradation, the percentage of aggregate smaller than 4.75 mm and the rejuvenator content, in the order given. Hence, it can be concluded that the total binder content, the virgin binder added to the recycled asphalt mixtures, and the intermediate-temperature PG of the virgin binder are the most important parameters affecting the fatigue index. In this regard, considering these parameters to enhance recycled asphalt’s fatigue performance can be an optimal approach.

3.2.2. Rutting Input Feature Performance

The importance weight values of random forest and gradient boosting and their average value for input variables in the rutting performance of asphalt mixtures containing RAP are in Table 4. As can be perceived, the most influential parameter on rutting performance is the PG span of the virgin binder added to the recycled asphalt mixture. This rheological property of binder (i.e., PG span) plays the main role in asphalt pavement rutting performance, and this result is in harmony with the outcomes of Taher et al. [105]. The following important feature in rutting is the total binder content in the mixture followed by the coarse aggregate/fine aggregate value, as the binder in mixes affects the friction between aggregates and causes these deformations [106]. Furthermore, investigations have shown that the coarser gradations in asphalt mixes perform better in rutting than fine gradations [107], which is in accordance with the results of the present study. RAP content in recycled asphalt mixes is the fourth feature affecting the rutting damage, which indicates the improvement in rutting performance due to the stiff binder of RAP [31]. The binder content in RAP, the percentage of fine aggregate, the high-temperature PG of added virgin binder, the virgin binder content, the nominal maximum aggregate size in gradation, and the rejuvenator content, are the following essential parameters on rutting, respectively. Therefore, optimizing the PG span of the virgin binder, the total binder content in the mixture, and the coarse aggregate/fine aggregate value should be considered the priority in the cases where rutting performance enhancement is investigated.

3.3. Error Histogram

The performance of prediction algorithms can be evaluated by comparing the predicted values with the measured values at laboratories. The lower distance of data points from the quality lines indicates higher accuracy and less error in the method [25]. It is noteworthy that the data points used for the charts belong to the training and testing datasets.

3.3.1. Accuracy of the Machine Learning Methods—Fatigue Prediction

Figure 5 depicts the predicted values versus the experimental values of the natural logarithm of the fatigue index. As can be seen from Figure 5, in COA-KNN, approximately all data points are on or close to the quality line, which means the appropriate correlation between the predicted and experimental values. Random forest and gradient boosting have acceptable performances compared to other conventional algorithms. The data points in decision tree show more distance from the quality line followed by multiple linear regression, indicating a higher error in contrast to other prediction algorithms. Overall, it can be concluded that the introduced ensemble model outperforms the previously developed methods with significantly higher accuracy and lower error.

3.3.2. Accuracy of the Machine Learning Methods—Rutting Prediction

The graphs of the predicted data points versus measured data points for the rutting performance of asphalt mixtures containing RAP are represented in Figure 6. As can be seen, the best prediction performance is for the newly developed algorithm (COA-KNN) due to the accumulation of all data points on or around the quality line except for just two data samples. Additionally, among the conventional algorithms, random forest and gradient boosting perform admissibly as most data points are near the quality line. The data points in the decision tree graph show further distance from the quality line, indicating a higher error in comparison to the foregoing methods. Likewise, except for one point, none of the data points in multiple linear regression are on the quality line, implying the predicted values of the model are not reliable. To this end, it can be concluded that the proposed ensemble regression method exceeds in performance compared to the mentioned techniques.

4. Conclusions

This study aims to accurately predict the fatigue and rutting performances of asphalt mixtures containing different contents of RAP. Contrary to previous investigations, fatigue and rutting indices are used instead of dynamic modulus and rut depth, as these indices directly indicate the fatigue and rutting life of asphalt pavements. Random forest, gradient boosting, decision tree, and multiple linear regression were utilized as conventional prediction tools. Meanwhile, a new hybrid ensemble method, called COA-KNN, was developed to improve the prediction accuracy, and the results were compared against conventional prediction algorithms. COA-KNN was further applied to two databases of fatigue and rutting performances. The following conclusions can be drawn from the results of this investigation:
  • The R2 values of COA-KNN test data are 0.07, 0.10, 0.26, and 0.35 more than that of the gradient boosting, random forest, decision tree, and multiple linear regression, respectively, in the fatigue performance prediction.
  • Applying COA-KNN to the fatigue database reduces the MAE test data by 0.01, 0.02, 0.04, and 0.07 ln (cycles) compared to gradient boosting, random forest, decision tree, and multiple linear regression, respectively.
  • COA-KNN MAPE test data are 0.04%, 0.11%, 0.21%, and 0.38% lower than gradient boosting, random forest, decision tree, and multiple linear regression, respectively, in fatigue performance prediction.
  • The MSE test data obtained using COA-KNN for fatigue prediction are 0.01, 0.01, 0.03, and 0.03 (ln (cycles))2 lower than gradient boosting, random forest, decision tree, and multiple linear regression, respectively.
  • The R2 values of test data attained using COA-KNN for rutting performance are 0.050, 0.091, 0.222, and 0.806 more than random forest, gradient boosting, decision tree, and multiple linear regression, respectively.
  • The MAE COA-KNN test data are 361.44, 890.44, 897.21, and 3123.89 inches lower than random forest, gradient boosting, decision tree, and multiple linear regression, respectively, in rutting performance prediction.
  • Replacing COA-KNN on the fatigue database with random forest, decision tree, gradient boosting, and multiple linear regression can reduce the MAPE test data by 1.44%, 5.06%, 5.70%, and 22.60%, respectively.
  • The MSE test data of COA-KNN are 1297190, 2506610, 7206180, and 19897280 inches2 less than the random forest, gradient boosting, decision tree, and multiple linear regression, respectively, for the rutting database.
  • Considering the performance indicators, COA-KNN outperforms the other conventional algorithms in fatigue and rutting predictions. Moreover, random forest and gradient boosting methods had an appropriate accuracy, compared to other methods, in both datasets.
  • Based on the importance weights of random forest and gradient boosting, the most effective features on fatigue performance of asphalt mixes containing RAP are the total binder content in the mix, the virgin binder content, and the intermediate-temperature PG of the virgin binder, respectively.
  • According to the ranking of the variables on the rutting characteristic, the PG span of the virgin binder, total binder content, and the coarse-to-fine aggregate ratio are the most effective parameters on the rutting damage of recycled asphalt pavements.
  • One of the limitations of this study is to apply a limited number of machine learning algorithms. Hence, it is recommended that other machine learning techniques, such as artificial neural networks or Gaussian Processes [108], will be applied to predict rutting and fatigue in future studies, and their performance will be compared with the proposed method in the current investigation.

Author Contributions

Methodology, H.N., H.J. and M.M.K.; Software, S.G. and H.N.; Validation, S.G., H.N., H.J. and M.M.K.; Formal analysis, S.G.; Resources, Z.A.; Data curation, Z.A.; Writing—original draft, Z.A.; Writing—review and editing, H.N., H.J. and M.M.K.; Supervision, H.J. and M.M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The table of abbreviations.
Table A1. The table of abbreviations.
RAPReclaimed asphalt pavement
COA-KNNCoyote Optimization Algorithm and K-Nearest Neighbor
RFRandom Forest regression
GBGradient Boosting regression
DTDecision Tree regression
MLRMultiple Linear Regression
GHGGlobal Greenhouse Gas
HMAHot Mix Asphalt
ANNArtificial Neural Networks
BBOBiogeography-Based Optimization algorithm
RRIRutting Resistance Index
PGPerformance Grade
RDRut Depth
MAEMean Absolute Error
MAPEMean Absolute Percentage Error
MSEMean Squared Error
NMASNominal Maximum Aggregate Size
RASReclaimed Asphalt Shingles
CRCrumbed Rubber

References

  1. Jahanbakhsh, H.; Karimi, M.M.; Jahangiri, B.; Nejad, F.M. Induction Heating and Healing of Carbon Black Modified Asphalt Concrete under Microwave Radiation. Constr. Build. Mater. 2018, 174, 656–666. [Google Scholar] [CrossRef]
  2. Karimi, M.M.; Jahanbakhsh, H.; Jahangiri, B.; Moghadas Nejad, F. Induced Heating-Healing Characterization of Activated Carbon Modified Asphalt Concrete under Microwave Radiation. Constr. Build. Mater. 2018, 178, 254–271. [Google Scholar] [CrossRef]
  3. Naseri, H.; Waygood, E.O.D.; Wang, B.; Patterson, Z.; Daziano, R.A. A Novel Feature Selection Technique to Better Predict Climate Change Stage of Change. Sustainability 2021, 14, 40. [Google Scholar] [CrossRef]
  4. Jahanbakhsh, H.; Karimi, M.M.; Naseri, H.; Nejad, F.M. Sustainable Asphalt Concrete Containing High Reclaimed Asphalt Pavements and Recycling Agents: Performance Assessment, Cost Analysis, and Environmental Impact. J. Clean. Prod. 2020, 244, 118837. [Google Scholar] [CrossRef]
  5. Karimi, M.M.; Amani, S.; Jahanbakhsh, H.; Jahangiri, B.; Alavi, A.H. Induced Heating-Healing of Conductive Asphalt Concrete as a Sustainable Repairing Technique: A Review. Clean. Eng. Technol. 2021, 4, 100188. [Google Scholar] [CrossRef]
  6. Kamali, Z.; Karimi, M.M.; Dehaghi, E.A.; Jahanbakhsh, H. Using Electromagnetic Radiation for Producing Reclaimed Asphalt Pavement (RAP) Mixtures: Mechanical, Induced Heating, and Sustainability Assessments. Constr. Build. Mater. 2022, 321, 126315. [Google Scholar] [CrossRef]
  7. Xing, C.; Li, M.; Liu, L.; Lu, R.; Liu, N.; Wu, W.; Yuan, D. A Comprehensive Review on the Blending Condition between Virgin and RAP Asphalt Binders in Hot Recycled Asphalt Mixtures: Mechanisms, Evaluation Methods, and Influencing Factors. J. Clean. Prod. 2023, 398, 136515. [Google Scholar] [CrossRef]
  8. Praticò, F.G.; Giunta, M.; Mistretta, M.; Gulotta, T.M. Energy and Environmental Life Cycle Assessment of Sustainable Pavement Materials and Technologies for Urban Roads. Sustainability 2020, 12, 704. [Google Scholar] [CrossRef] [Green Version]
  9. Diab, A.; Sangiorgi, C.; Ghabchi, R.; Zaman, M.; Wahaballa, A.M. Warm Mix Asphalt (WMA) Technologies: Benefits and Drawbacks—A Literature Review. In Functional Pavement Design; CRC Press: Boca Raton, FL, USA, 2017; pp. 1145–1154. [Google Scholar] [CrossRef]
  10. Abdollahi, S.F.; Karimi, M.M.; Jahanbakhsh, H.; Tabatabaee, N. Cracking Performance of Rubberized RAP Mixtures with Sasobit. Constr. Build. Mater. 2022, 319, 126090. [Google Scholar] [CrossRef]
  11. Ziari, H.; Aliha, M.R.M.; Moniri, A.; Saghafi, Y. Crack Resistance of Hot Mix Asphalt Containing Different Percentages of Reclaimed Asphalt Pavement and Glass Fiber. Constr. Build. Mater. 2020, 230, 117015. [Google Scholar] [CrossRef]
  12. Yousefi, A.; Behnood, A.; Nowruzi, A.; Haghshenas, H. Performance Evaluation of Asphalt Mixtures Containing Warm Mix Asphalt (WMA) Additives and Reclaimed Asphalt Pavement (RAP). Constr. Build. Mater. 2021, 268, 121200. [Google Scholar] [CrossRef]
  13. Jahanbakhsh, H.; Hosseini, P.; Nejad, F.M.; Habibi, M. Intermediate Temperature Fracture Resistance Evaluation of Cement Emulsified Asphalt Mortar. Constr. Build. Mater. 2019, 197, 1–11. [Google Scholar] [CrossRef]
  14. Behbahani, H.; Ayazi, M.J.; Moniri, A. Laboratory Investigation of Rutting Performance of Warm Mix Asphalt Containing High Content of Reclaimed Asphalt Pavement. Pet. Sci. Technol. 2017, 35, 1556–1561. [Google Scholar] [CrossRef]
  15. Colbert, B.; You, Z. The Determination of Mechanical Performance of Laboratory Produced Hot Mix Asphalt Mixtures Using Controlled RAP and Virgin Aggregate Size Fractions. Constr. Build. Mater. 2012, 26, 655–662. [Google Scholar] [CrossRef]
  16. Tran, N.H.; Taylor, A.; Willis, R. Effect of Rejuvenator on Performance Properties of HMA Mixtures with High RAP and RAS Contents. NCAT Rep. 2012, 1, 5–12. [Google Scholar]
  17. Izaks, R.; Haritonovs, V.; Klasa, I.; Zaumanis, M. Hot Mix Asphalt with High RAP Content. Procedia Eng. 2015, 114, 676–684. [Google Scholar] [CrossRef] [Green Version]
  18. Guduru, G.; Goli, A.K.; Matolia, S.; Kuna, K.K. Chemical and Performance Characteristics of Rejuvenated Bituminous Materials with High Reclaimed Asphalt Content. J. Mater. Civ. Eng. 2021, 33, 04020434. [Google Scholar] [CrossRef]
  19. Miró, R.; Valdés, G.; Martínez, A.; Segura, P.; Rodríguez, C. Evaluation of High Modulus Mixture Behaviour with High Reclaimed Asphalt Pavement (RAP) Percentages for Sustainable Road Construction. Constr. Build. Mater. 2011, 25, 3854–3862. [Google Scholar] [CrossRef]
  20. Behnood, A.; Daneshvar, D. A Machine Learning Study of the Dynamic Modulus of Asphalt Concretes: An Application of M5P Model Tree Algorithm. Constr. Build. Mater. 2020, 262, 120544. [Google Scholar] [CrossRef]
  21. Garbowski, T.; Pożarycki, A. Multi-Level Backcalculation Algorithm for Robust Determination of Pavement Layers Parameters. Inverse Probl. Sci. Eng. 2017, 25, 674–693. [Google Scholar] [CrossRef]
  22. Sarkhani Benemaran, R.; Esmaeili-Falak, M.; Javadi, A. Predicting Resilient Modulus of Flexible Pavement Foundation Using Extreme Gradient Boosting Based Optimised Models. Int. J. Pavement Eng. 2022, 1–20. [Google Scholar] [CrossRef]
  23. Esmaeili-Falak, M.; Benemaran, R.S. Ensemble Deep Learning-Based Models to Predict the Resilient Modulus of Modified Base Materials Subjected to Wet-Dry Cycles. Geomech. Eng. 2023, 32, 583–600. [Google Scholar]
  24. Ghasemi, P.; Aslani, M.; Rollins, D.K.; Williams, R.C. Principal Component Neural Networks for Modeling, Prediction, and Optimization of Hot Mix Asphalt Dynamics Modulus. Infrastructures 2019, 4, 53. [Google Scholar] [CrossRef] [Green Version]
  25. Behnood, A.; Mohammadi Golafshani, E. Predicting the Dynamic Modulus of Asphalt Mixture Using Machine Learning Techniques: An Application of Multi Biogeography-Based Programming. Constr. Build. Mater. 2021, 266, 120983. [Google Scholar] [CrossRef]
  26. Kaloush, K.E.; Witczak, M.W. Development of a Permanent to Elastic Strain Ratio Model for Asphalt Mixtures; Development of the 2002 Guide for the Design of New and Rehabilitated Pavement Structures. NCHRP 1-37 A, Inter Team Technical Report; Scientific Research Publishing: Wuhan, China, 2000. [Google Scholar]
  27. Ullah, S.; Tanyu, B.F.; Zainab, B. Development of an Artificial Neural Network (ANN)-Based Model to Predict Permanent Deformation of Base Course Containing Reclaimed Asphalt Pavement (RAP). Road Mater. Pavement Des. 2021, 22, 2552–2570. [Google Scholar] [CrossRef]
  28. Majidifard, H.; Jahangiri, B.; Rath, P.; Urra Contreras, L.; Buttlar, W.G.; Alavi, A.H. Developing a Prediction Model for Rutting Depth of Asphalt Mixtures Using Gene Expression Programming. Constr. Build. Mater. 2021, 267, 120543. [Google Scholar] [CrossRef]
  29. Mogawer, W.S.; Austerman, A.J.; Bonaquist, R. Determining the Influence of Plant Type and Production Parameters on Performance of Plant-Produced Reclaimed Asphalt Pavement Mixtures. Transp. Res. Rec. 2012, 2268, 71–81. [Google Scholar] [CrossRef]
  30. Elkashef, M.; Williams, R.C. Improving Fatigue and Low Temperature Performance of 100% RAP Mixtures Using a Soybean-Derived Rejuvenator. Constr. Build. Mater. 2017, 151, 345–352. [Google Scholar] [CrossRef]
  31. Mogawer, W.; Austerman, A.; Mohammad, L.; Kutay, M.E. Evaluation of High RAP-WMA Asphalt Rubber Mixtures. Road Mater. Pavement Des. 2013, 14, 129–147. [Google Scholar] [CrossRef]
  32. Valdés, G.; Pérez-Jiménez, F.; Miró, R.; Martínez, A.; Botella, R. Experimental Study of Recycled Asphalt Mixtures with High Percentages of Reclaimed Asphalt Pavement (RAP). Constr. Build. Mater. 2011, 25, 1289–1297. [Google Scholar] [CrossRef] [Green Version]
  33. Nabizadeh, H.; Haghshenas, H.F.; Kim, Y.R.; Aragão, F.T.S. Effects of Rejuvenators on High-RAP Mixtures Based on Laboratory Tests of Asphalt Concrete (AC) Mixtures and Fine Aggregate Matrix (FAM) Mixtures. Constr. Build. Mater. 2017, 152, 65–73. [Google Scholar] [CrossRef]
  34. Bonicelli, A.; Calvi, P.; Martinez-Arguelles, G.; Fuentes, L.; Giustozzi, F. Experimental Study on the Use of Rejuvenators and Plastomeric Polymers for Improving Durability of High RAP Content Asphalt Mixtures. Constr. Build. Mater. 2017, 155, 37–44. [Google Scholar] [CrossRef]
  35. Mogawer, W.S.; Booshehrian, A.; Vahidi, S.; Austerman, A. Evaluating the Effect of Rejuvenators on the Degree of Blending and Performance of High RAP, RAS, RAP/RAS Mixtures. Assoc. Asph. Paving Technol. AAPT 2013, 82, 403–436. [Google Scholar] [CrossRef]
  36. Tran, N.; Xie, Z.; Julian, G.; Taylor, A.; Willis, R.; Robbins, M.; Buchanan, S. Effect of a Recycling Agent on the Performance of High-RAP and High-RAS Mixtures: Field and Lab Experiments. J. Mater. Civ. Eng. 2017, 29, 04016178. [Google Scholar] [CrossRef]
  37. Abreu, L.P.F.; Oliveira, J.R.M.; Silva, H.M.R.D.; Fonseca, P.V. Recycled Asphalt Mixtures Produced with High Percentage of Different Waste Materials. Constr. Build. Mater. 2015, 84, 230–238. [Google Scholar] [CrossRef] [Green Version]
  38. Cooper, S.B.; Rouge, B. Characterization of Hma Mixtures Containing High Recycled Asphalt Pavement Content with Crumb Rubber Additives; Louisiana State University and Agricultural & Mechanical College: Baton Rouge, LA, USA, 2008. [Google Scholar]
  39. Bueche, N.; Dumont, A.G.; Vanelstraete, A.; De Visscher, J.; Vansteenkiste, S.; Vervaecke, F.; Gaspar, L.; Thogersen, F. Laboratory and Alt-Evaluation of High Stiffness Underlayers With High Percentage of Re-Use as Developed in the NR2C-Project. In Proceedings of the 4th Eurasphalt and Eurobitume Congress, Copenhagen, Denmark, 21–23 May 2008. [Google Scholar]
  40. Büchler, S.; Falchetto, A.C.; Walther, A.; Riccardi, C.; Wang, D.; Wistuba, M.P. Wearing Course Mixtures Prepared with High Reclaimed Asphalt Pavement Content Modified by Rejuvenators. Transp. Res. Rec. 2018, 2672, 96–106. [Google Scholar] [CrossRef]
  41. Sabouri, M. Evaluation of Performance-Based Mix Design for Asphalt Mixtures Containing Reclaimed Asphalt Pavement (RAP). Constr. Build. Mater. 2020, 235, 117545. [Google Scholar] [CrossRef]
  42. Arshadi, A.; Steger, R.; Ghabchi, R.; Zaman, M.; Hobson, K.; Commuri, S. Performance Evaluation of Plant-Produced Warm Mix Asphalts Containing RAP and RAS. Asph. Paving Technol. Assoc. Asph. Paving Technol.-Proc. Tech. Sess. 2017, 86, 403–425. [Google Scholar] [CrossRef]
  43. Daniel, J.S.; Mogawer, W.S. Determining the Effective PG Grade of Binder in RAP Mixes; The New England Transportation Consortium: Concord, NH, USA, 2010. [Google Scholar]
  44. Tavakol, M.; Hossain, M.; Heptig, B. Minimum Virgin Binder Content Needed in Recycled Superpave Mixtures to Resist Fatigue Cracking and Moisture Damage. J. Mater. Civ. Eng. 2018, 30, 04018126. [Google Scholar] [CrossRef]
  45. Vahidi, S.; Mogawer, W.S.; Booshehrian, A. Effects of GTR and Treated GTR on Asphalt Binder and High-RAP Mixtures. J. Mater. Civ. Eng. 2014, 26, 721–727. [Google Scholar] [CrossRef]
  46. Pasetto, M.; Giacomello, G.; Pasquini, E. Effectiveness of Rejuvenators for Asphalt Mixtures with High Reclaimed Asphalt Pavement Content in Cold Climates. In Proceedings of the International Symposium on Asphalt Pavement & Environment, Padua, Italy, 11–13 September 2019; pp. 3–13. [Google Scholar]
  47. Ozer, H.; Al-Qadi, I.; Kanaan, A.; Lippert, D. Performance Characterization of Asphalt Mixtures at High Asphalt Binder Replacement with Recycled Asphalt Shingles. Transp. Res. Rec. 2013, 2371, 105–112. [Google Scholar] [CrossRef]
  48. Haghshenas, H.; Nabizadeh, H.; Kim, Y.-R.; Santosh, K. Research on High-RAP Asphalt Mixtures with Rejuvenators and WMA Additives; University of Nebraska: Lincoln, NE, USA, 2019. [Google Scholar]
  49. Sabahfar, N.; Ahmed, A.; Aziz, S.R.; Hossain, M. Cracking Resistance Evaluation of Mixtures with High Percentages of Reclaimed Asphalt Pavement. J. Mater. Civ. Eng. 2017, 29, 06016022. [Google Scholar] [CrossRef]
  50. Ahmed, A. Evaluation of Cracking Potential of Superpave Mixtures with High Reclaimed Asphalt Pavement Content. Master’s Thesis, Kansas State University, Manhattan, KS, USA, 2015. [Google Scholar]
  51. Cong, P.; Zhang, Y.; Liu, N. Investigation of the Properties of Asphalt Mixtures Incorporating Reclaimed SBS Modified Asphalt Pavement. Constr. Build. Mater. 2016, 113, 334–340. [Google Scholar] [CrossRef]
  52. Mogawer, W.S.; Fini, E.H.; Austerman, A.J.; Booshehrian, A.; Zada, B. Performance Characteristics of High Reclaimed Asphalt Pavement Containing Bio-Modifier. Road Mater. Pavement Des. 2016, 17, 753–767. [Google Scholar] [CrossRef]
  53. Zhang, J.; Guo, C.; Chen, T.; Zhang, W.; Yao, K.; Fan, C.; Liang, M.; Guo, C.; Yao, Z. Evaluation on the Mechanical Performance of Recycled Asphalt Mixtures Incorporated with High Percentage of RAP and Self-Developed Rejuvenators. Constr. Build. Mater. 2021, 269, 121337. [Google Scholar] [CrossRef]
  54. Ma, T.; Ding, X.; Zhang, D.; Huang, X.; Chen, J. Experimental Study of Recycled Asphalt Concrete Modified by High-Modulus Agent. Constr. Build. Mater. 2016, 128, 128–135. [Google Scholar] [CrossRef]
  55. Kocak, S.; Kutay, M.E. Use of Crumb Rubber in Lieu of Binder Grade Bumping for Mixtures with High Percentage of Reclaimed Asphalt Pavement. Road Mater. Pavement Des. 2017, 18, 116–129. [Google Scholar] [CrossRef]
  56. Baaj, H.; Ech, M.; Tapsoba, N.; Sauzeat, C.; Di Benedetto, H. Thermomechanical Characterization of Asphalt Mixtures Modified with High Contents of Asphalt Shingle Modifier (ASM®) and Reclaimed Asphalt Pavement (RAP). Mater. Struct./Mater. Constr. 2013, 46, 1747–1763. [Google Scholar] [CrossRef]
  57. Hossain, M.; Musty, H.Y.; Sabahfer, N. Use of High-Volume Reclaimed Asphalt Pavement (RAP) for Asphalt Pavement Rehabilitation Due to Increased Highway Truck Traffic from Freight Transportation; University of Nebraska: Lincoln, NE, USA, 2012. [Google Scholar]
  58. Sanchez, X.; Tighe, S.L. Steps towards the Detection of Reclaimed Asphalt Pavement in Superpave Mixtures. Road Mater. Pavement Des. 2019, 20, 1201–1214. [Google Scholar] [CrossRef]
  59. Ozer, H.; Al-Qadi, I.L.; Carpenter, S.H.; Aurangzeb, Q.; Roberts, G.L.; Trepanier, J. Evaluation of RAP Impact on Hot-Mix Asphalt Design and Performance. Asph. Paving Technol. Assoc. Asph. Paving Technol. Proc. Tech. Sess. 2009, 78, 317–348. [Google Scholar]
  60. Safi, F.R.; Al-Qadi, I.L.; Hossain, K.; Ozer, H. Total Recycled Asphalt Mixes: Characteristics and Field Performance. Transp. Res. Rec. 2019, 2673, 149–162. [Google Scholar] [CrossRef]
  61. Zhu, J.; Ma, T.; Fan, J.; Fang, Z.; Chen, T.; Zhou, Y. Experimental Study of High Modulus Asphalt Mixture Containing Reclaimed Asphalt Pavement. J. Clean. Prod. 2020, 263, 121447. [Google Scholar] [CrossRef]
  62. Boriack, P.; Katicha, S.W.; Flintsch, G.W. Laboratory Study on Effects of High Reclaimed Asphalt Pavement and Binder Content. Transp. Res. Rec. J. Transp. Res. Board 2014, 2445, 64–74. [Google Scholar] [CrossRef]
  63. Pérez-martínez, M.; Marsac, P.; Gabet, T.; Hammoum, F.; Lopes, M.; Pouget, S. Durability Analysis of Different Warm Mix Asphalt Containing Reclaimed Asphalt Pavement. In Proceedings of the XXVth WRC Seoul Seoul, Seoul, Republic of Korea, 2–6 November 2015; pp. 1–15. [Google Scholar]
  64. Xie, Z.; Tran, N.; Julian, G.; Taylor, A.; Blackburn, L.D. Performance of Asphalt Mixtures with High Recycled Contents Using Rejuvenators and Warm-Mix Additive: Field and Lab Experiments. J. Mater. Civ. Eng. 2017, 29, 04017190. [Google Scholar] [CrossRef]
  65. Tran, N.; Taylor, A.; Turner, P.; Holmes, C.; Porot, L. Effect of Rejuvenator on Performance Characteristics of High RAP Mixture. Road Mater. Pavement Des. 2017, 18, 183–208. [Google Scholar] [CrossRef]
  66. Al-Qadi, I.L.; Qazi, A.; Carpenter, S.H. Impact of High RAP Content on Structural and Performance Properties of Asphalt Mixtures; Research Report FHWA-ICT-12-002; University of Illinois, Urbana-Champaign: Champaign, IL, USA, 2012; pp. 1–107. [Google Scholar]
  67. Daneshvar, D.; Behnood, A. Estimation of the Dynamic Modulus of Asphalt Concretes Using Random Forests Algorithm. Int. J. Pavement Eng. 2020, 23, 250–260. [Google Scholar] [CrossRef]
  68. Huang, Y.H. Pavement Analysis and Design (Second Edition); Pearson: London, UK, 2004; ISBN 0131424734. [Google Scholar]
  69. European Standard EN 12697-24; Bituminous Mixtures—Test Methods for Hot Mix Asphalt—Part 24: Resistance to Fatigue. British Standards Institution: London, UK, 2005.
  70. Silva, H.M.R.D.; Oliveira, J.R.M.; Jesus, C.M.G. Are Totally Recycled Hot Mix Asphalts a Sustainable Alternative for Road Paving? Resour. Conserv. Recycl. 2012, 60, 38–48. [Google Scholar] [CrossRef]
  71. Zaumanis, M.; Arraigada, M.; Wyss, S.A.; Zeyer, K.; Cavalli, M.C.; Poulikakos, L.D. Performance-Based Design of 100% Recycled Hot-Mix Asphalt and Validation Using Traffic Load Simulator. J. Clean. Prod. 2019, 237, 117679. [Google Scholar] [CrossRef]
  72. Pradyumna, T.A.; Mittal, A.; Jain, P.K. Characterization of Reclaimed Asphalt Pavement (RAP) for Use in Bituminous Road Construction. Procedia Soc. Behav. Sci. 2013, 104, 1149–1157. [Google Scholar] [CrossRef] [Green Version]
  73. Cooper, S.B.; Mohammad, L.N.; Elseifi, M.A. Laboratory Performance of Asphalt Mixtures Containing Recycled Asphalt Shingles, Reclaimed Asphalt Pavement, and Recycling Agents. J. Mater. Civ. Eng. 2017, 29, D4016001. [Google Scholar] [CrossRef]
  74. Pradhan, S.K.; Sahoo, U.C. Influence of Softer Binder and Rejuvenator on Bituminous Mixtures Containing Reclaimed Asphalt Pavement (RAP) Material. Int. J. Transp. Sci. Technol. 2021, 11, 46–59. [Google Scholar] [CrossRef]
  75. Pradhan, S.K.; Sahoo, U.C. Evaluation of Recycled Asphalt Mixtures Rejuvenated with Madhuca Longifolia (Mahua) Oil. Int. J. Pavement Res. Technol. 2021, 14, 43–53. [Google Scholar] [CrossRef]
  76. Celauro, C.; Bernardo, C.; Gabriele, B. Production of Innovative, Recycled and High-Performance Asphalt for Road Pavements. Resour. Conserv. Recycl. 2010, 54, 337–347. [Google Scholar] [CrossRef]
  77. Maupin, G.W.; Diefenderfer, S.D.; Gillespie, J.S. Evaluation of Using Higher Percentages of Recycled Asphalt Pavement in Asphalt Mixes in Virginia; Virginia Transportation Research Council: Charlottesville, VA, USA, 2008. [Google Scholar]
  78. Widyatmoko, I. Mechanistic-Empirical Mixture Design for Hot Mix Asphalt Pavement Recycling. Constr. Build. Mater. 2008, 22, 77–87. [Google Scholar] [CrossRef]
  79. Hill, B.; Behnia, B.; Buttlar, W.G.; Reis, H. Evaluation of Warm Mix Asphalt Mixtures Containing Reclaimed Asphalt Pavement through Mechanical Performance Tests and an Acoustic Emission Approach. J. Mater. Civ. Eng. 2013, 25, 1887–1897. [Google Scholar] [CrossRef]
  80. Vargas, A. Evaluation of the Use of Reclaimed Asphalt Pavement in Stone Matrix Asphalt Mixtures. Master’s Thesis, Auburn University, Auburn, AL, USA, 2007. [Google Scholar]
  81. Pradhan, S.K.; Sahoo, U.C. Effectiveness of Pongamia Pinnata Oil as Rejuvenator for Higher Utilization of Reclaimed Asphalt (RAP) Material. Innov. Infrastruct. Solut. 2020, 5, 92. [Google Scholar] [CrossRef]
  82. Gedafa, D.S.; Berg, A.; Karki, B.; Saha, R.; Melaku, R.S. Cracking and Rutting Performance of Field and Laboratory HMA Mixes. In Proceedings of the Airfield and Highway Pavements 2019: Testing and Characterization of Pavement Materials—Selected Papers from the International Airfield and Highway Pavements Conference 2019, Chicago, IL, USA, 18 July 2019; pp. 12–19. [Google Scholar] [CrossRef]
  83. Aurangzeb, Q.; Al-Qadi, I.L.; Pine, W.J.; Trepanier, J.S.; Abuawad, I.M. Thermal Cracking Potential in Asphalt Mixtures with High RAP Contents. In Proceedings of the 7th RILEM International Conference on Cracking in Pavements, Delft, The Netherlands, 20–22 June 2012; Springer: Berlin/Heidelberg, Germany, 2012; Volume 4, pp. 1271–1280. [Google Scholar] [CrossRef]
  84. Al-Saffar, Z.H.; Yaacob, H.; Mohd Satar, M.K.I.; Mohd Usak, S.N.; Jaya, R.P.; Hassan, N.A.; Radeef, H.R.; Warid, M.N.M. Evaluating the Performance of Reclaimed Asphalt Pavement Incorporating PelletRAP as a Rejuvenator. In Proceedings of the 4th National Conference on Wind & Earthquake Engineering, Putrajaya, Malaysia, 16–17 October 2020; IOP Publishing: Bristol, UK, 2021; Volume 682. [Google Scholar] [CrossRef]
  85. Kim, S.; Shen, J.; Myung Jeong, M. Effects of Aggregate Size on the Rutting and Stripping Resistance of Recycled Asphalt Mixtures. J. Mater. Civ. Eng. 2018, 30, 04017280. [Google Scholar] [CrossRef]
  86. Pradhan, S.K.; Sahoo, U.C. Effectiveness of Polanga Oil as Rejuvenator for Asphalt with High RAP Content. In Proceedings of the Airfield and Highway Pavements 2019: Testing and Characterization of Pavement Materials—Selected Papers from the International Airfield and Highway Pavements Conference 2019, Chicago, IL, USA, 18 July 2019; pp. 114–126. [Google Scholar] [CrossRef]
  87. Jamal; Hafeez, I.; Yaseen, G.; Aziz, A. Influence of Cereclor on the Performance of Aged Asphalt Binder. Int. J. Pavement Eng. 2020, 21, 1309–1320. [Google Scholar] [CrossRef]
  88. Al-Saffar, Z.H.; Yaacob, H.; Mohd Satar, M.K.I.; Putra Jaya, R. The Tailored Traits of Reclaimed Asphalt Pavement Incorporating Maltene: Performance Analyses. Int. J. Pavement Eng. 2020, 23, 1800–1813. [Google Scholar] [CrossRef]
  89. Zhang, W.; Shen, S.; Wu, S.; Mohammad, L.N. Prediction Model for Field Rut Depth of Asphalt Pavement Based on Hamburg Wheel Tracking Test Properties. J. Mater. Civ. Eng. 2017, 29, 04017098. [Google Scholar] [CrossRef]
  90. Rau, C.S.; Wu, S.C.; Chien, P.C.; Kuo, P.J.; Chen, Y.C.; Hsieh, H.Y.; Hsieh, C.H. Prediction of Mortality in Patients with Isolated Traumatic Subarachnoid Hemorrhage Using a Decision Tree Classifier: A Retrospective Analysis Based on a Trauma Registry System. Int. J. Environ. Res. Public Health 2017, 14, 1420. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Naseri, H.; Shokoohi, M.; Jahanbakhsh, H.; Karimi, M.M. A Novel Soft Computing Approach to Better Predict Flexible Pavements Roughness. SSRN 4191842. Available online: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=4191842 (accessed on 2 May 2023).
  92. Wei, Z.; Meng, Y.; Zhang, W.; Peng, J.; Meng, L. Downscaling SMAP Soil Moisture Estimation with Gradient Boosting Decision Tree Regression over the Tibetan Plateau. Remote Sens. Environ. 2019, 225, 30–44. [Google Scholar] [CrossRef]
  93. Naseri, H.; Jahanbakhsh, H.; Moghadas Nejad, F. Developing a Novel Machine Learning Method to Predict the Compressive Strength of Fly Ash Concrete in Different Ages. AUT J. Civ. Eng. 2019, 4, 423–436. [Google Scholar] [CrossRef]
  94. Pierezan, J.; Maidl, G.; Massashi Yamao, E.; dos Santos Coelho, L.; Cocco Mariani, V. Cultural Coyote Optimization Algorithm Applied to a Heavy Duty Gas Turbine Operation. Energy Convers. Manag. 2019, 199, 111932. [Google Scholar] [CrossRef]
  95. Naseri, H.; Ehsani, M.; Golroo, A.; Moghadas Nejad, F. Sustainable Pavement Maintenance and Rehabilitation Planning Using Differential Evolutionary Programming and Coyote Optimisation Algorithm. Int. J. Pavement Eng. 2021, 23, 2870–2887. [Google Scholar] [CrossRef]
  96. Thanh Noi, P.; Kappas, M. Comparison of Random Forest, k-Nearest Neighbor, and Support Vector Machine Classifiers for Land Cover Classification Using Sentinel-2 Imagery. Sensors 2017, 18, 18. [Google Scholar] [CrossRef] [Green Version]
  97. Naseri, H.; Shokoohi, M.; Jahanbakhsh, H.; Golroo, A.; Gandomi, A.H. Evolutionary and Swarm Intelligence Algorithms on Pavement Maintenance and Rehabilitation Planning. Int. J. Pavement Eng. 2022, 23, 4649–4663. [Google Scholar] [CrossRef]
  98. Naseri, H.; Jahanbakhsh, H.; Hosseini, P.; Moghadas Nejad, F. Designing Sustainable Concrete Mixture by Developing a New Machine Learning Technique. J. Clean. Prod. 2020, 258, 120578. [Google Scholar] [CrossRef]
  99. Naseri, H.; Hosseini, P.; Jahanbakhsh, H.; Hosseini, P.; Gandomi, A.H. A Novel Evolutionary Learning to Prepare Sustainable Concrete Mixtures with Supplementary Cementitious Materials. Environ. Dev. Sustain. 2022, 25, 5831–5865. [Google Scholar] [CrossRef]
  100. Naseri, H.; Jahanbakhsh, H.; Khezri, K.; Shirzadi Javid, A.A. Toward Sustainability in Optimizing the Fly Ash Concrete Mixture Ingredients by Introducing a New Prediction Algorithm. Environ. Dev. Sustain. 2022, 24, 2767–2803. [Google Scholar] [CrossRef]
  101. Naseri, H.; Jahanbakhsh, H.; Foomajd, A.; Galustanian, N.; Karimi, M.M.; Waygood, E.O.D. A Newly Developed Hybrid Method on Pavement Maintenance and Rehabilitation Optimization Applying Whale Optimization Algorithm and Random Forest Regression. Int. J. Pavement Eng. 2022, 1–13. [Google Scholar] [CrossRef]
  102. Mrówczyński, D.; Gajewski, T.; Garbowski, T. Sensitivity Analysis of Open-Top Cartons in Terms of Compressive Strength Capacity. Materials 2023, 16, 412. [Google Scholar] [CrossRef]
  103. Sreedhar, S.; Coleri, E. Effects of Binder Content, Density, Gradation, and Polymer Modification on Cracking and Rutting Resistance of Asphalt Mixtures Used in Oregon. J. Mater. Civ. Eng. 2018, 30, 04018298. [Google Scholar] [CrossRef]
  104. FHA. Administration Superpave Fundamentals; FHA: Washington, DC, USA, 2008. [Google Scholar]
  105. Taher, B.M.; Mohamed, R.K.; Mahrez, A. A Review on Fatigue and Rutting Performance of Asphalt Mixes. Sci. Res. Essays 2011, 6, 670–682. [Google Scholar]
  106. Golalipour, A.; Jamshidi, E.; Niazi, Y.; Afsharikia, Z.; Khadem, M. Effect of Aggregate Gradation on Rutting of Asphalt Pavements. Procedia Soc. Behav. Sci. 2012, 53, 440–449. [Google Scholar] [CrossRef] [Green Version]
  107. Ahmed, M.A.; Attia, M.I.E. Impact of Aggregate Gradation and Type on Hot Mix Asphalt Rutting In Egypt. Int. J. Eng. Res. Appl. (IJERA) 2013, 3, 2249–2258. [Google Scholar]
  108. Garbowski, T. Stochastic Model Reduction Applied to Inverse Analysis. In Proceedings of the VI International Conference on Adaptive Modeling and Simulation ADMOS 2013, Lisbon, Portugal, 3–5 June 2013. [Google Scholar]
Figure 1. The methodology flowchart of this investigation.
Figure 1. The methodology flowchart of this investigation.
Buildings 13 01787 g001
Figure 2. Distribution of (a): fatting index (b) rutting index.
Figure 2. Distribution of (a): fatting index (b) rutting index.
Buildings 13 01787 g002
Figure 3. Performance of different prediction algorithms on the fatigue index prediction problem: (a) R2, (b) MAE, (c) MAPE, (d) MSE, (e) VAF, and (f) A10-index.
Figure 3. Performance of different prediction algorithms on the fatigue index prediction problem: (a) R2, (b) MAE, (c) MAPE, (d) MSE, (e) VAF, and (f) A10-index.
Buildings 13 01787 g003aBuildings 13 01787 g003bBuildings 13 01787 g003c
Figure 4. Performance of different prediction algorithms on the rutting index prediction problem: (a) R2, (b) MAE, (c) MAPE, (d) MSE, (e) VAF, and (f) A10-index.
Figure 4. Performance of different prediction algorithms on the rutting index prediction problem: (a) R2, (b) MAE, (c) MAPE, (d) MSE, (e) VAF, and (f) A10-index.
Buildings 13 01787 g004aBuildings 13 01787 g004bBuildings 13 01787 g004c
Figure 5. Predicted vs. experimental fatigue index values of (a) Random Forest, (b) Gradient Boosting, (c) Decision Tree, (d) Multiple Linear Regression, and (e) COA-KNN.
Figure 5. Predicted vs. experimental fatigue index values of (a) Random Forest, (b) Gradient Boosting, (c) Decision Tree, (d) Multiple Linear Regression, and (e) COA-KNN.
Buildings 13 01787 g005aBuildings 13 01787 g005bBuildings 13 01787 g005c
Figure 6. Predicted vs. experimental rutting index values of (a) Random Forest, (b) Gradient Boosting, (c) Decision Tree, (d) Multiple Linear Regression, and (e) COA-KNN.
Figure 6. Predicted vs. experimental rutting index values of (a) Random Forest, (b) Gradient Boosting, (c) Decision Tree, (d) Multiple Linear Regression, and (e) COA-KNN.
Buildings 13 01787 g006aBuildings 13 01787 g006bBuildings 13 01787 g006c
Table 1. Description of the variables and their descriptive statistics for fatigue index.
Table 1. Description of the variables and their descriptive statistics for fatigue index.
VariableTypeDescriptionDataPearson Correlation with Ln (Nf)MaximumMinimumAverageVarianceStandard ErrorKurtosis
NfFatigue
Index
Calculated fatigue indexTraining- 1.25 × 10 8 2.43 × 10 7 5.88 × 10 7 1.91 × 10 7 1.52 × 10 6 0.91
Testing- 1.09 × 10 8 3.09 × 10 7 6.14 × 10 7 1.87 × 10 7 3.15 × 10 6 0.14
Ln (Nf)Fatigue
Index
Ln of Calculated
fatigue index
Training118.6417.0117.840.310.02−0.01
Testing118.5017.2517.890.290.05−0.10
TemperatureTest
Condition
Dynamic modulus test temperature (°C)Training0.00525.004.0019.413.150.256.72
Testing0.00525.0015.0019.892.300.391.70
RAPRAPRAP content (%)Training−0.292100.000.0031.7621.951.751.20
Testing−0.38470.000.0030.7418.683.16−0.31
PGSpanBinderSpan of PG of binderTraining0.23298.0062.0086.479.280.740.39
Testing0.12798.0062.0087.718.761.480.60
PGInterBinderIntermediate-temperature PG of binderTraining−0.24640.0010.0022.924.670.371.90
Testing−0.27837.0010.0022.435.460.920.90
RejuvenatorRejuvenatorRejuvenator
content (%)
Training0.05213.800.001.303.250.264.45
Testing0.01012.000.001.453.160.533.87
ACVirginVolumetricVirgin asphalt content (%)Training0.3266.740.003.691.310.100.48
Testing0.3978.001.003.621.290.222.98
ACRAPVolumetricRAP asphalt
content (%)
Training−0.1177.900.004.271.930.150.98
Testing−0.2458.900.004.391.990.341.61
ACTotalVolumetricTotal asphalt
content (%)
Training0.3668.003.705.320.760.060.54
Testing0.4288.004.005.280.750.133.54
NMASGradationNominal maximum aggregate size (mm)Training−0.21720.004.7513.623.550.28−0.13
Testing−0.14520.004.7513.873.560.60−0.01
Fine agg.GradationAggregate smaller than 4.75 mm (%)Training0.21593.006.1057.5413.171.050.68
Testing−0.14876.8033.9056.0711.902.01−0.97
Course agg./Fine agg.GradationAggregate larger than 4.75 mm/Aggregate smaller than 4.75 mmTraining0.14913.290.061.691.390.1131.89
Testing−0.2143.310.511.470.740.120.24
Table 2. Description of the variables and their descriptive statistics for the rutting index.
Table 2. Description of the variables and their descriptive statistics for the rutting index.
VariableTypeDescriptionDataPearson Correlation with RRIMaximumMinimumAverageVarianceStandard ErrorKurtosis
RRIRutting
Index
Calculated
rutting index
Training137,007.904094.4912,767.915666.04450.771.88
Testing127,401.606305.5112,577.054975.55853.300.23
RAPRAPRAP content (%)Training0.008100.000.0036.0424.971.990.56
Testing−0.058100.000.0035.9618.273.133.70
PGSpanBinderSpan of PG of binderTraining−0.14698.0062.0082.669.380.75−0.26
Testing−0.09698.0068.0084.508.041.380.13
PGHighBinderHigh-temperature PG of binderTraining−0.08476.0046.0060.247.870.63−0.29
Testing−0.03076.0046.0060.748.181.40−0.01
RejuvenatorRejuvenatorRejuvenator content (%)Training−0.09115.000.001.883.400.272.08
Testing−0.0819.280.000.842.410.417.55
ACVirginVolumetricVirgin asphalt content (%)Training−0.0388.000.003.891.380.110.71
Testing−0.1245.700.513.751.010.172.33
ACRAPVolumetricRAP asphalt content (%)Training0.0637.900.003.941.850.150.59
Testing−0.2656.200.004.551.310.226.97
ACTotalVolumetricTotal asphalt content (%)Training−0.1918.003.705.540.760.060.64
Testing−0.4026.604.005.400.610.100.04
NMASGradationNominal
maximum aggregate size (mm)
Training0.09125.004.7514.695.270.42−0.47
Testing0.34625.004.7513.413.790.651.85
Fine agg.GradationAggregate smaller than 4.75 mm (%)Training0.08788.9010.0053.7416.631.32−0.51
Testing−0.13093.0025.2052.1615.402.640.32
Course agg./Fine agg.GradationAggregate larger than 4.75 mm
/Aggregate smaller than 4.75 mm
Training−0.2129.000.121.100.950.0829.89
Testing−0.0852.970.081.110.690.121.21
Table 3. Ranking of input features in fatigue.
Table 3. Ranking of input features in fatigue.
Random ForestGradient BoostingAverage
Input
Features
Importance WeightRankingInput
Features
Importance WeightRankingInput
Features
Importance WeightRanking
ACTotal0.2141ACTotal0.3131ACTotal0.2641
ACVirgin0.1672RAP0.1232ACVirgin0.1322
PGInter0.1163PGInter0.1053PGInter0.1103
Temperature0.0914ACVirgin0.0974RAP0.1044
PGSpan0.0865Temperature0.0865Temperature0.0885
RAP0.0856ACRAP0.0786ACRAP0.0766
ACRAP0.0747NMAS0.0497PGSpan0.0677
Course agg./Fine agg.0.0518PGSpan0.0478Course agg./Fine agg.0.0488
Fine agg.0.0489Course agg./Fine agg.0.0449NMAS0.0439
NMAS0.03810Fine agg.0.03210Fine agg.0.04010
Rejuvenator0.03111Rejuvenator0.02411Rejuvenator0.02711
Table 4. Ranking of input features in rutting.
Table 4. Ranking of input features in rutting.
Random ForestGradient BoostingAverage
Input
Features
Importance WeightRankingInput
Features
Importance WeightRankingInput
Features
Importance WeightRanking
PGSpan0.2531PGSpan0.2991PGSpan0.2761
Course agg./Fine agg.0.1542ACTotal0.1842ACTotal0.1642
ACTotal0.1453RAP0.1363Course agg./Fine agg.0.1323
ACRAP0.1404Course agg./Fine agg.0.1104RAP0.1144
RAP0.0925Fine agg.0.0965ACRAP0.0925
Fine agg.0.0606PGHigh0.0716Fine agg.0.0786
PGHigh0.0517ACRAP0.0457PGHigh0.0617
ACVirgin0.0508ACVirgin0.0428ACVirgin0.0468
NMAS0.0499Rejuvenator0.0119NMAS0.0289
Rejuvenator0.00610NMAS0.00610Rejuvenator0.00810
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ghavami, S.; Alipour, Z.; Naseri, H.; Jahanbakhsh, H.; Karimi, M.M. A New Ensemble Prediction Method for Reclaimed Asphalt Pavement (RAP) Mixtures Containing Different Constituents. Buildings 2023, 13, 1787. https://doi.org/10.3390/buildings13071787

AMA Style

Ghavami S, Alipour Z, Naseri H, Jahanbakhsh H, Karimi MM. A New Ensemble Prediction Method for Reclaimed Asphalt Pavement (RAP) Mixtures Containing Different Constituents. Buildings. 2023; 13(7):1787. https://doi.org/10.3390/buildings13071787

Chicago/Turabian Style

Ghavami, Sadegh, Zeynab Alipour, Hamed Naseri, Hamid Jahanbakhsh, and Mohammad M. Karimi. 2023. "A New Ensemble Prediction Method for Reclaimed Asphalt Pavement (RAP) Mixtures Containing Different Constituents" Buildings 13, no. 7: 1787. https://doi.org/10.3390/buildings13071787

APA Style

Ghavami, S., Alipour, Z., Naseri, H., Jahanbakhsh, H., & Karimi, M. M. (2023). A New Ensemble Prediction Method for Reclaimed Asphalt Pavement (RAP) Mixtures Containing Different Constituents. Buildings, 13(7), 1787. https://doi.org/10.3390/buildings13071787

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