Prediction of the hot asphalt mix properties using deep neural networks

Marshall design process is the most common method used for estimating the Optimum asphalt content (OAC) and this process is called the asphalt mix design. However, this method is time-consuming, labor-intensive, and its results are subjected to variations. This paper employs artificial neural network (ANN) for the estimation of Marshall test parameters (OAC, Stability, Flow, Air voids, Voids in mineral aggregate) using the aggregate gradation as the input of the prediction process. Multiple ANNs are tested in order to optimize the NN hyperparameters and produce accurate predictions. Different activation functions, number of hidden layers, and number of neurons per hidden layer are tested and heatmaps are generated to compare the performance of every ANN. Results show that the optimum ANN hyperparameters change depending on the predicted parameter. Finally, the deep NN can predict the OAC, stability, flow, density, air voids, and voids in mineral aggregate with R values of 0.91, 0.8, 0.53, 0.65, 0.77, and 0.66. The linear activation function is the most efficient activation function and generates more accurate results than the logistic and the hyperbolic tangent functions. Additionally, it is shown that the deep neural network approach represents a major innovative tool for the prediction of the asphalt mix properties as results of this approach outperforms results of the shallow ANN that consists of a single hidden layer which is the only approach used in the literature. Thus, the use of the deep ANN can be useful during the phase of the design of the asphalt mix process because of its ability to predict variables with high accuracy. For example, the ANN with 3 hidden layers and 16 neurons per layer with the linear activation function can predict the OAC with high accuracy (R = 0.91), which can be helpful in the design process as the ANN can be employed for the prediction of the OAC of the asphalt mix.


Background
Asphalt concrete was originally invented with the purpose of carrying the high pressure and heavy loads generated from aircrafts [1]. In general, pavement consists of multiple layers to transfer the heavy loads to the soil without causing soil failure [2]. The process of asphalt pavement design consists of two main processes. The first process is the thickness design which focuses on estimating the required thickness of each pavement layer to transfer traffic load safely to the soil. The second process, which is the main focus of this paper, is called the asphalt mix design and it focuses on estimating the optimum aggregate and bitumen characteristics in the mix [3][4][5][6][7][8]. Usually, the asphalt mix design process is carried out through laboratory tests with the goal of estimating the optimum asphalt content (OAC) in the asphalt mix [9] as the OAC has a significant influence on the final performance of the mix [10]. Low asphalt content leads to harsh mix which causes durability issues, while high

Open Access
Beni-Suef University Journal of Basic and Applied Sciences asphalt content leads to rutting, flushing, and insufficient air voids [11].
Marshall mix design method was firstly proposed in 1939 by Bruce Marshall in the Mississippi State Highway Department [1]. Marshall test was widely adopted in 1948 with slight modifications from country to country [9]. Currently, Marshall test remains the most common way used for estimating the optimum asphalt content [12,13]. Marshall test procedure requires the preparation of at least fifteen samples for five different asphalt contents (three samples for each) then draw the design curves [14,15] and estimate the OAC that satisfies a predefined criterion which depends on the Maximum stability, Maximum density or unit weight, predefined air voids percent, and a min value for voids in mineral aggregate [9]. Finally, the OAC is estimated as the average value that corresponds to the maximum stability, maximum unit weight, and prespecified air voids percentage [14,15]. As a result, the estimation of the OAC is subjected to significant deviations, as it relies on the average value of a group of different ACs. Additionally, Marshall design method requires significant time for sample preparation and testing [14,16]. Consequently, a large number of studies focus on developing alternative methods for Marshall test to save time.
AC is a composite material that consists of aggregate and bitumen. Aggregate makes up a high proportion of volume and mass of mixtures (around 95% of the mix weight), hence it is considered as the most important constituent of asphalt concrete [17] so characteristics of the asphalt mix mainly depend on the aggregate used and its gradation [18]. Thus, it is predicted that the aggregate properties have an enormous impact on the mixture properties. Aggregate gradation can be defined as the distribution of particle sizes expressed as a percentage of the total weight [18,19] and considered as the centerpiece property of aggregate which needs careful consideration due to its effect on mix properties and performance of HMA mixtures, including air void, stability, stiffness, durability, permeability, workability, fatigue resistance, frictional resistance, resistance to moisture damage [18,20] and also rutting resistance of asphalt concrete under traffic and environmental loads [21]. As a result, this parameter is considered as a very important parameter in the process of mixture design. For example, changing the aggregate gradation changes the required asphalt content because the asphalt content is responsible for coating the aggregate surface and filling the voids between the aggregate particles [16]. Moreover, the main sources of the mix strength or stability are the friction between the aggregate particles, interlocking resistance between the aggregate particles and the consistency of the bitumen used [14,17]. Therefore, the main objective of this paper is to test the effectiveness of using ANN for accurate and fast prediction of the asphalt mix properties (mechanical and volumetric properties) to enhance the mix design process in the laboratory. The ANNs developed in this study were trained and tested on the basis of laboratory data of HMAs for two types of aggregate. The tests were carried out in the Highway and Airport Engineering Research Laboratory, Cairo University, Egypt for the most common aggregate (3D, 4C) and bitumen (60/70) used in Egypt.
Over the past few years, ANNs have been used in the prediction problem and indeed in the pavement engineering field. The study by Kaseko and Ritchie (1993) is one of the oldest studies that employed ANN for the detection of pavement cracks [22]. This was followed by Gagarin, Flood, and Albrecht (1994) study that employed ANN for the estimation of different truck attributes from the strain response readings from bridges [23]. In 1995, Cal employed an ANN that takes plasticity index, water content, and liquid limit as an input for the prediction of the soil classification [24]. In 1998, Roberts and Attoh-Okine employed an ANN for the prediction of the pavement condition [25]. Similarly, the study by Attoh-Okine (2001) proposed an ANN that uses the pavement characteristics and the surrounding conditions such as weather, age, and traffic condition in the prediction of the pavement condition [26]. In 2010, Tapkın et al. proposed an ANN for estimating Marshal test results for dense asphalt mixes that are modified with different types and percentages of an additive which is the polypropylene [10]. In 2014, Khuntia et al. used ANN for the prediction of Marshall test results for the polyethylene modified mixes using two techniques which are ANNs and support vector machine [27]. In 2016, Ozturk et al. proposed an ANN for the prediction of the volumetric properties of the asphalt mix [28]. In 2017, Ivica and Ivan proposed an ANN for the prediction of the air voids and the asphalt content of asphalt mixes in Croatia [29]. In 2018, Baldo et al. used the ANN technique for the prediction of the mechanical properties of asphalt mixes (stability and flow) [9]. In 2019, Nguyen et al. used artificial intelligence techniques for the prediction of the characteristics of the stone matrix asphalt [30].
Since 2006, when Hinton proposed several techniques for deep learning structures, configuring deep networks with more than three layers have shown widespread success in training neural networks. Artificial Neural Networks (ANN) is a learning algorithm that imitates the human neural system. An ANN consists of multiple nodes, called neurons, that communicate through synapses. Typically, there are three sets of layers: input layer, hidden layer, and output layer and each type of layer plays unique roles. The Input layer receives input information, the output layer yields output signals, and the intermediate layer (hidden layer) receives signals from the input layer and manipulates the information to give results to the output layer. An ANN model can have multiple intermediate layers [31]. In general, ANNs can be defined as parallel distributed processors that have the potential to store the learned knowledge and use it in the future. The fundamental unit of any ANN is the artificial neuron. The function of the neuron is to process the input signals and modulate its own response through an activation function or sometimes called the transfer function. The activation function determines the interruption or transmission of the outgoing impulse. Each neuron computes a weighted sum of elements of the input vector (Xs) through weights associated with the connections (W) [32]. Then, the neuron output value is calculated by applying the assigned transfer function to the weighted sum as follows: There are many types of neural networks such as feedforward and feedback. Also, there are many types of training techniques deepening on the data such as supervised and unsupervised learning. In this study, the supervised learning is employed for the proposed ANN. Additionally, the backpropagation learning algorithm is the most common training algorithm used for the training of ANN, so the backpropagation technique is used for the training of the proposed ANN in this paper. Additionally multiple ANN will be tested to optimize the ANN architecture.

Material selection
In general, two aggregate gradations are commonly used in Egypt (4C and 3D) [10][11][12][13][14][15][16][17][18]. Between 2015 and 2016, the highway and research laboratory, Cairo University, Egypt collected 2073 asphalt mix samples from different locations all over Egypt. Out of these 2073 samples, 1186 are wearing course samples, and 887 are binder course samples. Results of laboratory tests show that gradation 4C is the most common gradation used in the wearing course, and gradation 3D is the most common gradation used in the binder course samples followed by grade gradation 4C as shown in Fig. 1. On the other hand, the figure shows the poor quality control in sites as most of the collected samples do not satisfy any aggregate gradations according to the Egyptian code for urban and rural roads.
For the bitumen or binder used in Egypt, the most common binder grade used in Egypt is binder 60-70. These results comply with the results of two recent studies by Mousa et al. (2018) and Othman and Abdelwahab (2021) [10,18] in Egypt which showed that binder grade 60-70 and aggregate 4C and 3D are the most common binder and aggregate gradations used in Egypt [10].

Characteristics of the materials used
In this study, 102 asphalt mix samples are prepared and tested using bitumen 60-70 and dolomite aggregate gradations 4C and 3D. Properties of the aggregate used (3D and 4C) are summarized in Table 1 and tier limits are shown in a sieve size distribution graph in Fig. 2. Additionally, properties of the bitumen used (60/70) are summarized in Table 2.

Marshall test
For every asphalt mix, five binder contents were chosen, and three samples were prepared and tested for every single binder content. Thus, for each mix, at least 15 samples were prepared and tested in the standard Marshall test to estimate the flow, stability values, and the volumetric characteristics. Then, for every binder content, the average value of these three samples was used as an indicator for the stability, flow, and the remaining volumetric properties. At this point, the five design curves were drawn to estimate the OAC that corresponds to the maximum stability, maximum density, and specific percentage of air voids in the mix as follows: In this paper, the standard Marshall test is used to estimate the optimum asphalt content for 102 asphalt mixes samples prepared using bitumen 60-70 and aggregates (3D and 4C). Out of these 102 samples, 63 were prepared using aggregate grade 3D and 39 were prepared using aggregate grade 3C.

ANN structure
The accuracy of any ANN depends on setting the hyperparameters and selecting the appropriate ANN architecture. However, till the moment, there is no scientific way for the selection of the hyperparameters or the architecture. Thus, the optimum ANN architecture is usually selected following the trial-and-error technique. Thus, in this study multiple ANN are tested using different activation functions, different number of hidden layers, and different number of neurons per hidden layer.
The activation function or the transfer function may take multiple analytical expressions. In this study, three different activation functions are tested. These activation functions are shown in Fig. 3 and can be expressed as follows:  Figure 4 shows the general architecture of the ANNs tested in this study. Additionally, while initial weights of the ANN significantly influence the final performance of the ANN and the speed of the training process, there is no specific approach for estimating the best initial weights. Thus, initial weights are assigned randomly at the beginning of the training process.
For the training of the ANN, the dataset was randomly divided into two sets: training set (61%-62 samples) and testing set (49%-40 samples) while the backpropagation technique was used to train the data for 500 iterations. For the testing set, 50% of the testing samples were used to monitor the performance of the training process to avoid overfitting and this segment of the dataset is called the validation set. Additionally, the early stopping technique was employed to avoid overfitting. Overfitting occurs when the training set error declines while the validation or the testing set error increases. Thus, the early stopping technique stops the training once the validation set error starts to increase [34,35]. However, the training of the neural network is a stochastic process, and the error of the validation set might go up or down at any point. Thus, the first overfitting point might not be the perfect point to stop the training. As a result, the performance of the ANN is monitored for a number of iterations and if the validation set error keeps increasing the training process should stop. On the other hand, if the validation set error declines, the training process should continue [36]. In this study, 100 iterations are used to monitor the validation set error. After these 100 iterations, if the validation set error declined again the training process proceeds, other than that the training process is halted.

Results
The selection of the optimum ANN hyperparameters should be based on the evaluation of the ANNs' ability to predict the outputs [37,38], which through statistical indicators. The main goal of the evaluation process or testing process is to estimate the ANN generalization error. The term generalization can be defined as the ability of the model to return high accuracy results on unknown data that is not used in the training process. The most common statistical indicator is the coefficient of correlation (R) that measures the linear relation between the ANN predictions (X prediction ) and the target values of the actual values ( X actual ) and it is calculated as follows: The coefficient of correlation of the different ANNs for the testing set are shown in Tables 3, 4, 5, 6, 7 and 8, and for visualization purposes, the cells are highlighted from red which indicates the lowest value to green which indicates the highest value according to their R-value. Table 3 reports the R values in the prediction of the OAC for the different ANNs. Results show that the logistic activation function produces the highest R-value, and the highest R-value (= 0.91) is achieved using three hidden layers and 16 neurons per layer. Additionally, for the three activation functions, no specific pattern can be observed with respect to the number of neurons per hidden layer or the number of hidden layers in the ANN. Table 4 reports the R values in the prediction of the stability value for the different ANNs. Results show that, in general, the R values of the linear activation function are better than the R-values produced using the other activation functions. On the other hand, the highest R-value (= 0.8) is achieved using an ANN with three hidden layers, 9 neurons per layer, and the hyperbolic tanh activation function. Additionally, for the logistic and tanh activation functions, no specific pattern can be observed with respect to the number of neurons per hidden layer or the number of hidden layers in the ANN. However, for the linear activation function, the performance of the ANN gets worth with the increase in the number of the hidden layers. Table 5 reports the R-values in the prediction of the flow value for the different ANNs. Results show that, in general, the R-values of the linear activation function are better than the R values produced using the other activation functions. On the other hand, the highest R-value (= 0.53) is achieved using an ANN with two hidden layers, 15 neurons per layer, and the logistic activation function. Additionally, for the three activation functions, the performance of the ANN, in general, gets worth with the increase in the number of the hidden layers. Table 6 reports the R-values in the prediction of the density value for the different ANNs. Results show that, in general, the R-values of the linear activation function are better than the R values produced using the other activation functions, and the highest R-value (= 0.65) is achieved using an ANN with two hidden layers, 11 neurons per layer. Additionally, for the three activation functions, the performance of the ANN, in general, gets worth with the increase in the number of the hidden layers.

Discussion
It can be seen that the linear activation function produces the best R-values; however, the linear activation function does not provide the highest R-value for all the outputs. Additionally, the optimum ANN architecture changes according to the desired output. Thus, in this section, a search for the optimum ANN that produces the best results for all outputs will be investigated.

Optimum ANN architecture for all output variables
Instead of choosing the optimum ANN that optimizes the performance for a single desired output. This section focuses on choosing the optimum ANN architecture that optimizes that performance for all predictions. In other words, this section focuses on choosing the optimum ANN architecture that produces the best predictions for all the outputs. Thus, a weight is assigned to each output R-value and then the weighted sum is calculated for every signal ANN. In this case, the new R-value is calculated as follows: The updated R-values are shown in Table 9. It can be seen that the linear activation function, in general, provides the best results for all the predictors, followed by the logistic activation function, then the hyperbolic tanh activation function. From the updated table, the optimum ANN that provides the best predictions for all the six predictions has two hidden layers, 11 neurons per layer, and using the linear activation function as shown in Fig. 5. A comparison between the best possible ANN architecture and the optimum ANN for all predictions is shown in Table 10. It is clear that the difference between the two ANN is not significant for the OAC, stability, flow, density, and AV; however, there is a considerable error in the predictions of the VMA for the optimum ANN for all predictions when compared with the optimum ANN for the VMA.

Optimum ANN for the prediction of the asphalt mix mechanical characteristics (stability and flow)
Similar to the previous section, this section focuses on choosing the optimum ANN architecture that optimizes that ANN predictions of the asphalt mix mechanical properties. Thus, a weight is assigned to each output R-value and then the weighted sum is The updated R-values are shown in Table 11. It can be seen that the linear activation function, in general, provides the best results for all the predictors, followed by the logistic activation function, then the hyperbolic tanh activation function. From the updated table, results show that the optimum ANN that provides the best predictions of the asphalt mix mechanical properties is the same as the optimum ANN for all predictions. This ANN has two hidden layers, 11 neurons per layer, and using the linear activation function.

Optimum ANN for the prediction of the asphalt mix volumetric characteristics (AV, VMA, and density)
Similar to the previous two section, this section focuses on choosing the optimum ANN architecture that optimizes that ANN predictions of the asphalt mix volumetric properties. Thus, a weight is assigned to each output R-value and then the weighted sum is calculated for every signal ANN. In this case, the new R-value is calculated as follows: The updated R-values are shown in Table 12. It can be seen that the linear activation function, in general, provides the best results for all the predictors, followed by the logistic activation function, then the hyperbolic tanh activation function. From the updated table, results show that the optimum ANN that provides the best predictions of the asphalt mix volumetric properties is the same as the optimum ANN for all predictions and the same as the optimum ANN for the prediction of the asphalt mix mechanical properties. This ANN has two hidden layers, 11 neurons per layer, and using the linear activation function.

Conclusions
This study focuses on the prediction of the OAC, mechanical, and volumetric properties of the hot asphalt mix from the aggregate gradation using ANN. 204 different ANN are tested using different activation functions, number of hidden layers, and number of neurons per layer in order to choose the ANN that optimizes the prediction accuracy based on results of 102 asphalt mix samples that were prepared and tested in the standard Marshall test. The main outcomes of this study can be summarized as follows: -The optimum ANN architecture changes depending on the desired output as shown in Table 13.

Number of Hidden Layers
Linear acƟvaƟon funcƟon L ogisƟc acƟvaƟon funcƟon T anh acƟvaƟon funcƟon  1  2  3  4  1  2  3  4  1  2  -The use of the hyperbolic tangent function provides the worst level of accuracy, followed by the logistic activation function. On the other hand, the linear activation function is the most efficient activation function. -The optimum ANN for the prediction of all the six output variables consists of two hidden layers, every hidden layer has elven neurons, and with a linear activation function. This architecture is also the best architecture for the prediction of either the mechanical properties of the mix, or the volumetric properties.
-Results of this study show that the development of deep ANN represents a major innovative tool for the prediction of the asphalt mix properties as it is verified that the deep ANN is useful for the predictions of any of the desired outputs as summarized in Table 1. -The use of the ANN can be useful during the phase of the design of the asphalt mix process because of its ability to predict variables with high accuracy. For example, the ANN with 3 hidden layers and 16 neurons per layer with the linear activation function can predict the OAC with high accuracy (R = 0.91), which  can be helpful in the design process as the ANN can be employed for the design of the hot asphalt mix. In this case, the ANN can be used for the prediction of the OAC instead of the standard Marshall test procedures that require the preparation of fifteen samples. Then, only three samples are prepared and tested to estimate the design parameters and make sure they match the design criteria. This approach saves time, resources, and the required effort to estimate the OAC.