Machine Learning Aided Design and Prediction of Environmentally Friendly Rubberised Concrete

: Not only can waste rubber enhance the properties of concrete (e.g., its dynamic damping and abrasion resistance capacity), its rational utilisation can also dramatically reduce environmental pollution and carbon footprint globally. This study is the world’s ﬁrst to develop a novel machine learning-aided design and prediction of environmentally friendly concrete using waste rubber, which can drive sustainable development of infrastructure systems towards net-zero emission, which saves time and cost. In this study, artiﬁcial neuron networks (ANN) have been established to determine the design relationship between various concrete mix composites and their multiple mechanical properties simultaneously. Interestingly, it is found that almost all previous studies on the ANNs could only predict one kind of mechanical property. To enable multiple mechanical property predictions, ANN models with various architectural algorithms, hidden neurons and layers are built and tailored for benchmarking in this study. Comprehensively, all three hundred and ﬁfty-three experimental data sets of rubberised concrete available in the open literature have been collected. In this study, the mechanical properties in focus consist of the compressive strength at day 7 (CS7), the compressive strength at day 28 (CS28), the ﬂexural strength (FS), the tensile strength (TS) and the elastic modulus (EM). The optimal ANN architecture has been identiﬁed by customising and benchmarking the algorithms (Levenberg–Marquardt (LM), Bayesian Regularisation (BR) and Scaled Conjugate Gradient (SCG)), hidden layers (1–2) and hidden neurons (1–30). The performance of the optimal ANN architecture has been assessed by employing the mean squared error (MSE) and the coefﬁcient of determination ( R 2 ) . In addition, the prediction accuracy of the optimal ANN model has ben compared with that of the multiple linear regression (MLR).


Introduction
Rubber or elastomer is a common material and is widely used as an essential material in the manufacture of tires. Because the demand for rubber has continued to increase over time, the global consumption of rubber in 2017 was 13,225 thousand metric tons of natural rubber and 15,189 thousand metric tons of synthetic rubber [1]. The generation of waste rubber in the EU is estimated to be more than 1.43 billion tons per year and has been growing at a rate comparable to the EU's economic growth. Nearly 5 billion tires, including stacked tires, will have been discarded by 2030 [2]. Thus, the utilisation of waste rubber resources is seen as an effective method for reducing their adverse effects on the environment, maintaining natural resources and reducing the demand for storage space [3]. At present, the main methods for the disposal of waste rubber are incineration and burial. There is a detrimental effect on the environment when waste rubber is burned because of the emissions of carbon dioxide and cyanide. According to the American Rubber Manufacturers Association Report, only approximately 5.5% of waste rubber is used for civil engineering. If more waste rubber is reused, more resources can be saved and negative effects on the environment can be reduced. As the primary source of waste rubber, scrapped tires, being an important waste material, have been studied and examined in relation to the field of construction [4]. The application of discarded automobile rubber tires in the civil engineering industry can be traced back to the last century. Waste rubber tires mainly ended up in landfill or were used as cushioning materials until the 1960s when rubber tires began to be used on a large scale in the construction industry due to the increasing amount of waste and enhanced environmental protection plans [5]. Therefore, waste rubber is now used as a substitute for aggregates (as either fine or coarse aggregates). When using discarded rubber tires to wholly or partially replace fine aggregates, the resulting concrete is lighter in weight. By using rubber tires instead of coarse aggregates, the elasticity and energy absorption capacity of the concrete are increased accordingly [6,7]. Moreover, the shrinkage of rubberised concrete increases with the increase in rubber sand content [8].
Research articles have shown a growth in the noise reduction coefficient with an increase in rubber replacing sand. Therefore, replacing aggregates with waste rubber in concrete at an opportune ratio not only supports the improvement of the performance of the concrete but also avoids the environmental pollution and waste of resources caused by conventional treatments [9,10]. The use of waste rubber can help engineers and stakeholders achieve a system's zero emission status in a faster and safer manner. This practice will fundamentally underpin the United Nation's sustainable development goals (SDGs) as well as the "race to zero" campaign.
Although there is a huge potential in the substitution of aggregates with waste rubber, the weakness in the mechanical properties and the durability of concrete due to the poor performance of rubber in bonding cement particles has been mentioned [11][12][13][14]. Thus, various methods for improving the durability and the mechanical properties of concrete have been proposed. For instance, the following treatments of rubber have been employed to improve the compressive strength of rubberised concrete: water soaking and washing, utilising rubber particles with large sizes, NaOH treatment, treatment with acetone or ethanol [15,16]. Furthermore, an innovative method, compression concrete casting, has been proposed to improve the compressive strength and elastic modulus of rubberised concrete. When 20% of the coarse aggregates in compressed concrete samples was replaced with rubber particles, the compressive strength and elastic modulus of the concrete were enhanced by 35% and 29%, respectively [17,18]. Regarding the flexural strength, the splitting tensile strength and the elastic modulus, the methods of water washing, water soaking and coating with cement paste were seen to enhance these properties. In order to obtain rubberised concrete with promising workability, rubber particles were added into a mixer with other concrete components, with the exception of water for dry mixing, and then water was added as the other components were mixed homogeneously [16].
Machine learning is often referred to as being part of the artificial intelligence used to analyse data to make smart decisions [19]. It is a method of realising artificial intelligence which has the ability to learn and predict data [20]. By adopting machine learning, it is possible to predict the performance of rubberised concretes that have different compositions. Machine learning can be developed by a variety of algorithms which are commonly classified into four different learning types according to their learning style: supervised, unsupervised, semi-supervised and reinforcement learning [21]. Supervised learning is suitable for data that have features and labels. In other words, data are provided to predict the labels. Unsupervised learning is only used in features with no labels, which means that data are provided to look for hidden structures. The difference between the above two styles is that supervised learning only uses labelled sample sets for learning, while unsupervised learning only uses unlabeled sample sets. For semi-supervised learning, some of the data are unlabelled, but most of it is labelled. Compared to supervised learning, the cost of semi-supervised learning is lower, but it can achieve higher accuracy. Reinforcement learning also uses unlabeled data, but it is possible to see whether it is getting closer or further away from the correct answer. In engineering design, computeraided methods such as machine learning and data statistics have been effectively used and can provide powerful benefits [22], especially when dealing with materials with complex variables and high uncertainty, such as composite materials [23]. A number of studies have shown that machine learning models have been widely applied and used as valuable tools for the prediction of the mechanical properties of concrete [24][25][26].
An artificial neuron network (ANN), a type of machine learning, is a simplified mathematical model that can simulate the function of natural biological neural networks to learn from past experience for solving new problems [27,28]. Since a large amount of data, such as compositions and properties of concrete, needs to be processed, the ordinary statistical methods cannot be sufficiently applied to the prediction of concrete properties. Furthermore, the prediction accuracy of the ordinary statistical methods may not be satisfied without proper algorithmic support. Based on previous experimental results, ANN models were used for predicting the mix proportion of polymer concrete to demonstrate the potential in saving time and costs [27,29]. Moreover, an ANN model with one hidden layer and 11 hidden neurons was utilised for predicting the compressive strength of concrete containing silica fume; 0.9724 of the R 2 value indicating the high prediction accuracy of the ANN model was obtained. An acceptable MSE value of the ANN model employed for estimating the compressive strength and elastic modulus of lightweight concrete was also acquired [30]. Therefore, it is clear that ANN can be selected as the suitable training method for this study. According to the previous studies related to the ANN for concrete property prediction, it can be found that the ANN was employed for predicting only one output in almost all previous studies. However, in reality, multiple mechanical properties are often required as part of engineering design and technical specifications. If an ANN model can predict multiple outputs (i.e., mechanical properties) with a satisfactory prediction accuracy, the optimal design and prediction of the concrete properties can be established, resulting in a reduction in material wastes and unnecessary costs. Thus, it is essential to automate the concrete design and prediction, which can improve waste management strategies towards net-zero built environments.
A typical structure of the ANN is demonstrated in Figure 1. It usually consists of three parts, the input layer, the output layer, and one or more hidden layers. Each layer has different numbers of neurons linked together by connections. In order to improve the accuracy of an ANN model, it is generally recommended to set one and two hidden layers containing multiple neurons. Moreover, the weight consists of the sum of regression coefficients and bias. The corresponding weight of layers can be added to each connection. An ANN model can be optimised by adjusting the weight during the training process until the error is reduced to an acceptable level [31]. Furthermore, the sigmoid function can be applied as the active function to analyse the effect of input elements and the weight on this element being processed [32]. For an ANN model, the number of hidden layers, connections and neurons are confirmed by the complexity of the raw data. The more complex the raw data are, the more hidden layers and neurons there are [33]. In order to obtain a high-accuracy model, the number of hidden layers and neurons of the model can be modified and compared in this study.
The aim of this research is to establish a novel machine learning apporach capable of designing and predicting multiple mechanical properties of rubberised concrete with various compositions. The ANN models are found to be capable of managing the complicated relationships between the inputs and outputs and of designing rubberised concrete that enhances resource conservation and environment protection by decreasing the experimental cost. Firstly, 335 experimental data sets of rubberised concrete properties with different compositions have been collected from published articles in open literature. Subsequently, the ANN models with different architectures (1-2 hidden layers and 1-30 hidden neurons) have been designed utilising MATLAB. Then, the optimal ANN architecture can be determined, followed by the performance evaluation. Finally, the prediction accuracy of multiple linear regression (MLR) conducted in the comparative analysis section is compared with that of the optimal ANN architecture capable of designing and predicting multiple mechanical properties of rubberised concrete.

Data Collection
The data used in this research have been collected from published articles in the open literature. In order to import data into the machine learning model in MATLAB, two items of the data sets for the pre-treatment method of rubber have been replaced by digitisation numbers described in Table 1. The data set size of the cement type and the ratio of rubber replacement are listed in Tables 2 and 3, respectively. Table 1. The representation of numbers for the pre-treatment rubber method.

1
No special treatment 2 Pre-treated with NaOH 3 Pre-coated with limestone  The collected rubberised concrete data in this research have been mainly categorised into three aspects: mandatory elements, characteristic elements and output elements as described below.

•
Mandatory Elements (ME) In this research, ME includes the percentage of rubber replacement (RR), the particle size of rubber (PSR), the proportion of fine aggregates (FA), the moisture content of fine aggregates (MCFA), the particle size of fine aggregates (PSFA), the proportion of rubber (R), the pre-treatment method of rubber (PR), the proportion of cement (C), the cement type (CT), the proportion of water (W), the proportion of water-reducing admixture (WRM), the proportion of coarse aggregates (CA), the particle size of coarse aggregates (CAPS), and the water-cement ratio (WCR).

•
Characteristic Elements (CE) CE indicates the parameters which are not included in all data sets, such as the proportion of slag (SG), the proportion of fly ash (FA) and the proportion of silica fume (SF).

• Output Elements (OE)
In this research, compressive strength at day 7 (CS7) and compressive strength at day 28 (CS28) of rubberised concrete, flexural strength (FS), splitting tensile strength (STS) and elastic modulus (EM) are considered as OE.
According to the aforementioned data sets classification, ME, CE and OE are the inputs and outputs of the ANN models accordingly. Table 4 shows the range of these parameters. In this study, data normalisation is proposed to reduce the negative influence of singular sample data in the intermited data clusters. Moreover, implementing data normalisation can avoid the overfitting problem. The reason is that different variables contain different dimensions, which may generate impacts on the data analysis [61]. Applying the data normalisation is able to limit data values within the range between zero and one that can enhance the comparability of data. The inputs and outputs in this research have been processed with the data normalisation method by utilising Equation (1) [62].
where, X i denotes the normalised data, X indicates the experimental data, X max and X min denote the maximum and minimum experimental data. By using the function, mapminmax, installed in MATLAB, the data normalisation was conducted.

Data Importation
Three hundred and fifty-three data sets are introduced to the ANN models for predicting mechanical properties of rubberised concrete. All imported data sets are randomly divided into three parts for training, validation and testing, respectively. The fixed allocation ratios of data sets between training, validation and testing aspects are 70%, 15% and 15%. The data sets for training are utilised for training models by modifying weights. The validation sets are used to adapt the model selection, that is, to do the final optimisation and determination of models, such as choosing the number of hidden neurons and hidden layers; while the testing set is purely to prove the generalisation of the trained models.
The raw data have been divided into the inputs and outputs data and then were imported into a workspace. Two data files are created in the workspace in this research, with "Input data" being a 17 × 353 matrix representing static data which has been composed of 353 samples of 17 elements. For those data sets without characteristic elements, the missed data are replaced by zero. Meanwhile, "Output data" is a 5 × 353 matrix, representing static data constituted of 353 samples of 5 elements.

Toolbox Selection
The neural net fitting toolbox installed in MATLAB has been selected as the training application in this study. It is noted that the output elements are continuous variables, which implies that a regression model is very suitable.

Hidden Layers and Neurons Determination
The numbers of hidden layers and neurons are essential to improve the accuracy of ANN models. Redundant and insufficient hidden neurons may cause overfitting and underfitting issues, respectively, due to inappropriate estimation of the relationship between the inputs and the outputs [30]. Moreover, according to previous research, there are three rules to assist in determining the appropriate number of hidden neurons [63,64].

•
The number of hidden neurons should be between the size of the input layer and the size of the output layer.

•
The number of hidden neurons should be 2/3 of the size of the input layer plus 2/3 of the size of the output layer.

•
The number of hidden neurons should be less than double the size of the input layer. • Besides, another method for determining neurons of ANN models is proposed in Equations (2)-(4) [65][66][67].
where N h denotes the neurons in the hidden layers and N i indicates the neurons in the numbers of inputs. Based on the rules above, the number of hidden neurons in this research has been set to 1, 5, 10, 15, 20, 25, 30 for each hidden layer and they are, respectively, substituted into the ANN models. For the number of hidden layers, there is not any guidance on how to specify it. Seven pertinent published articles related to the ANN architecture are listed in Table 5. It can be observed that 1-3 hidden layers and 5-40 hidden neurons were employed in the ANN models for predictions of concrete properties, respectively. It is clear that the customisation of ANN models is necessary. Thus, one and two hidden layers have been utilised in this research accordingly.

Algorithms Selection
Three algorithms are employed in this study to compare and indicate an optimal performance of ANN models, such as the accuracy and the time consumption of trained models. Levenberg-Marquardt (LM), Bayesian Regularisation (BR) and Scaled Conjugate Gradient (SCG) algorithms are utilised in this study and the detailed information of the algorithms are listed as follows.
• LM LM is an algorithm that provides a solution of the numerical nonlinear minimisation. The significance of LM algorithm is that it can simultaneously achieve the advantages of the Gauss-Newton method and the gradient descent algorithm by changing parameters. Furthermore, LM algorithm can improve the shortcomings of both algorithms. The LM algorithm is a type of upgraded Newton method shown in Equation (5) [74][75][76].
where, I indicates the identity matrix, e represents the vector, J is a Jacobian matrix, x k denotes the weight at epoch K, and u is a damping factor. In order to get more accurate models, u can be increased or dropped according to the success or failure of steps, and then the performance function can be enhanced.
• BR BR algorithm is capable of reaching the generalisation by applying an excellent combination of weights and square errors on the basis of LM optimisation. Equation (6) can be written to explain the objective function by employing the weights of networks [77].
where, E D denotes the value of errors, E w indicates the value of weights, and α and β represent the function parameters. Moreover, in order to determine the optimal α and β parameters, Equation (7) is proposed [78].
where, D stands for the distributed weight, M denotes the optimum architecture of networks, P (D|α, β, M), P (D|M) and P (α, β| M) indicate the likelihood function, the normalisation parameter and the initial regularisation factor, respectively. The operating processes of BR are as follows. Firstly, the optimum values of α and β are determined by Equation (7). Thereafter, the parameters, α and β, are confirmed by employing Bayes' theorem. Moreover, the optimum α and β are determined when the P (D|α, β, M) reaches the maximum value. Subsequently, the optimum weight is confirmed according to the minimum value of the Hessian matrix in LM operation. Finally, the values of α and β keep changing simultaneously until the convergent of models is reached [74,78].
• SCG SCG algorithm is a combination of LM and CG approaches [79]. With respect to the gradient descent algorithms, a costly linear searching direction must be determined. The response analysis of all data sets has to be repetitively conducted for multiple calculations during each direction searching. The application of SCG does not require the implementation of a linear search in each iteration, which can dramatically save time [78,80].
By following the selection approach of hidden layers, neutrons and algorithms selection, a variety of tailored architectures of the ANN models in this study can be listed in Table 6.

Performance Evaluation of ANN Architectures
In this study, customised combinations of different hidden neurons (1,5,10,15,20,25,30), hidden layers (1-2) and LM, BR and SCG algorithms have been employed. Furthermore, each pattern of ANN architectures has been trained for five times to improve the accuracy of models. In order to investigate the performance of the trained models in this study, two statistical analyses named MSE and R 2 are employed. MSE is the average value of the cost function for minimising the sum of squared errors (SSE) during the linear regression model fitting process. This represents the mean square error between the predicted and the actual value. MSE value is calculated according to Equation (8) [68]. The lower MSE value indicates a trained model with higher accuracy.
where n is the number of samples, and (y i − y i ) is the result of the experimental value minus the predicted value on the testing sets being processed. Moreover, R 2 value is employed as an assistance method to determine the performance of trained models defined in Equation (9). R 2 value exhibits the percentage of real value changes, which can be influenced by the variation of the predicted value. The range of R 2 value is from zero to one. Considering Equation (9), the numerator part represents the sum of the squared difference between the real value and the predicted value, similar to MSE. The denominator part represents the sum of the squared difference between the real value and the mean [70]. If the result is 0, it means that the model fits poorly. If the result is 1, it means that the model is error-free. Generally, the larger the R 2 value is, the better the model fitting effect is.
According to the aforementioned sections, the main steps to attain the optimum ANN architecture selection are demonstrated in Figure 2.

Optimum ANN Architecture Determination
The optimal ANN architecture is confirmed by embarking on the ANN architecture selection steps mentioned in the previous section. The MSE and R 2 values of the testing sets are defined as the performance evaluation index of all ANN architectures with 1-2 hidden layers and 1-30 hidden neurons, which are listed in Table A1. For instance, the ANN architecture LM-17-1-1-5, can be explained as the architecture employing the LM algorithm consisting of 17 input elements, five output elements, and one hidden neuron in two hidden layers. Figures 3 and 4 Figure 3b, the MSE values of all ANN architectures with two hidden layers go through slight drops within the first 15 hidden neurons in general. Thereafter, the MSE value of the ANN models with LM, BR and SCG algorithms decreases from 21.5373 to 7.2420, 28.2334 to 9.1582 and 25.4727 to 11.6450 when the hidden neurons are changed from one to ten, one to five and one to fifteen, respectively. Subsequently, the MSE values of the ANN with LM and BR algorithms soar between 7.2420 and 30.3840, and between 9.1582 and 54.9310 from ten to fifteen hidden layers and five to twenty hidden layers, respectively. Meanwhile, the MSE value of the ANN model with SCG slightly increases from 11.6450 to 22.2361 between 15 and 30 hidden neurons. Afterwards, the MSE values of the ANN models with LM and BR algorithms moderately drop from 30.3840 to 19.0260, and from 54.9310 to 39.5577 when the hidden neurons are changed from 15 to 30, and from 20 to 30, respectively. It can be concluded that the minimum MSE values of ANN architectures with one and two hidden layers are 8.7459 and 7.2420 derived from LM-17-15-5 and LM-17-10-10-5, respectively. The increasing MSE values of ANN models with various architectures can be attributed to two problems: (i) the underfitting problem, and (ii) the overfitting problem. Considering the underfitting problem, the learning capacity of ANN models is relatively low because the numbers of the hidden layers or neurons of ANN models are rather insufficient to stimulate the complicated relationships between the inputs and the outputs, then resulting in a weak generalisation capacity of ANN models. On the contrary, the overfitting problem occurs when the learning capacity of ANN models is overly high. In other words, every single datum can be captured by ANN models. Thus, the MSE value of ANN models soars because the generalisation capacity of ANN models decreases owing to the overfitting problem.
With regard to Figure 4a,b, the trend of changes in R 2 value is observed to be the same as the MSE value shown in Figure 3a,b. The maximum R 2 value of ANN architectures with one hidden layer and two hidden layers are 0.9626 and 0.9710 at LM-17-15-5 and LM-17-10-10-5, respectively. It can be summarised that the ANN architecture with LM algorithm, which consists of two hidden layers with ten neurons in each layer, shows the lowest MSE value but the highest R 2 value of 7.2420 and 0.9710, respectively. Thus, the best ANN architecture for predicting multiple mechanical properties of rubberised concrete is LM-17-10-10-5.
The optimal ANN architecture for predicting multiple mechanical properties of rubberised concrete has been determined by utilising different algorithms (LM, BR and SCG), hidden layers (one and two layers) and hidden neurons (1,5,10,15,20,25,30) in this study. Each of the ANN architectures contains 17 inputs and five outputs. The comparison of each experimental and predicted mechanical property by employing the optimal ANN architecture, LM-17-10-10-5, is exhibited in Figure 5. The R 2 value indicates the capacity of ANN models for predicting each mechanical property of rubberised concrete. The line charts demonstrate the difference between the experimental and the predicted value. Regarding Figure 5a-d,i,j, the R 2 values of CS7, CS28 and EM are 0.9552, 0.9641 and 0.9576, respectively. The R 2 value of these mechanical properties is recognised as a high prediction accuracy. The R 2 value of the ANN model, which predicts FS, is 0.8493 that is relatively lower than that of the ANN models for predicting CS7, CS28 and EM. Moreover, the R 2 of the ANN model for predicting STS, 0.6545, is the lowest among all ANN models in this study. Based on Figure 5h, the fitting line cannot coincide completely with the Y = T line. Namely, there is a certain degree of deviation between the two lines. This phenomenon also can be found in Figure 5g. For instance, the predicted STS is 5.70 MPa and 4.54 MPa, which are significantly lower than the experimental STS of 8.00 MPa and 7.00 MPa at sample 256 and sample 225, respectively. This phenomenon can be attributed to the underfitting of the model. The ANN model is too simple to explain the relationship between the inputs and STS.

Comparative Analysis
In the previous section, the optimal ANN architecture has been selected as LM-17-10-10-5 by comparisons using the MSE and R 2 values. Linear regression has been commonly utilised to simulate the complicated relationship between the dependent and independent factors. To inspect the prediction accuracy of the optimal ANN architecture, MLR has been adopted in this study. Therefore, the prediction accuracy of both methods is defined by employing R 2 . The higher the R 2 value, the better the prediction accuracy. The complicated relationship between 17 inputs and five outputs is calculated by employing MLR, explained in Equation (10) [30].
where, y indicates the predicted mechanical properties, x m denotes the independent variables, a 0 denotes the y-intercept, and a m indicates the regression coefficients. The regression coefficients are similar to the traditional regression models by applying the least square method shown in Equation (10). Therefore, the best function can be identified by minimising the sum of squared errors between experimental and predicted values. MLR has been conducted by employing "SPSS Statistics R27" in this study. The data sets are split at a ratio 8:2 for training and testing. The R 2 value of each predicted mechanical property derived from MLR and the optimal ANN model are demonstrated in Table 7. According to Table 7, the R 2 values of MLR for predicting the CS7, CS28, FS, STS and EM of rubberised concrete are 0.660, 0.673, 0.601, 0.460 and 0.773, respectively. It is evident that the R 2 value of MLR is relatively lower than that of the optimal ANN model. It can be interpreted that the prediction accuracy of MLR for predicting mechanical properties of rubberised concrete is lower than that of the ANN models. Wherein, the attention should be paid to STS prediction of MLR. Note that the R 2 value of MLR for STS prediction (around 0.460) is much lower than those of the other prediction models for mechanical properties. The reason can be attributed to the fact that the relationship between the inputs and the tensile strength is nonlinear. Of interest, the same phenomenon occurs when STS is predicted by employing ANN models.

Conclusions
This study is the world's first to establish a novel machine learning-aided design and prediction of eco-friendly rubberised concrete, enhancing engineering applications for sustainable infrastructures towards net-zero emission. The advanced machine learning approach is capable of designing and predcting multiple attributes (i.e., mechanical properties) simultaneously, which is a key novelty of this study. This approach is more rational and practical since, in reality, engineers need to satisfy all limit states criteria (for both serviceability and ultimate conditions). To enable the study, a comprehensive collection of 353 data sets consisting of 17 input elements of pertinent rubberised concrete including its five mechanical properties as outputs have been collected from all reputable sources in the open literature, and have been processed for training and testing a variety of ANN models. ANN models with 1-2 hidden layers, 1-30 hidden neurons and three types of algorithms (LM, BR and SCG) have been designed, validated and evaluated by benchmarking the R 2 values and MSE values. Subsequently, the optimal ANN architecture, which best predicts the outcome, has been customised and obtained as LM-17-10-10-5. Then, the R 2 value acting as the prediction accuracy index of MLR was compared with that of the optimal ANN model. The conclusions can be drawn as follows: • The ANN architecture with LM algorithm, two hidden layers and ten hidden neurons in each hidden layer is the optimal option for simultaneously predicting multiple mechanical properties of eco-friendly rubberised concrete. • Based on the MSE (7.2420) and R 2 (0.9710) values of the optimal ANN architecture, excellent prediction accuracy of the machine learning can be attained.

•
The R 2 value of MLR is relatively lower than that of the optimal ANN model. This traditionally implies that the prediction accuracy of the ANN model is relatively higher than that of MLR.