Performance Evaluation of Soft Computing for Modeling the Strength Properties of Waste Substitute Green Concrete

: The waste disposal crisis and development of various types of concrete simulated by the construction industry has encouraged further research to safely utilize the wastes and develop accurate predictive models for estimation of concrete properties. In the present study, sugarcane bagasse ash (SCBA), a by-product from the agricultural industry, was processed and used in the production of green concrete. An advanced variant of machine learning, i.e., multi expression programming (MEP), was then used to develop predictive models for modeling the mechanical properties of SCBA substitute concrete. The most signiﬁcant parameters, i.e., water-to-cement ratio, SCBA replacement percentage, amount of cement, and quantity of coarse and ﬁne aggregate, were used as modeling inputs. The MEP models were developed and trained by the data acquired from the literature; furthermore, the modeling outcome was validated through laboratory obtained results. The accuracy of the models was then assessed by statistical criteria. The results revealed a good approximation capacity of the trained MEP models with correlation coefﬁcient above 0.9 and root means squared error (RMSE) value below 3.5 MPa. The results of cross-validation conﬁrmed a generalized outcome and the resolved modeling overﬁtting. The parametric study has reﬂected the effect of inputs in the modeling process. Hence, the MEP-based modeling followed by validation with laboratory results, cross-validation, and parametric study could be an effective approach for accurate modeling of the concrete properties.


Introduction
The cement production process results in the depletion of limited natural resources and also emits a significant amount of carbon dioxide to the atmosphere [1][2][3]. Moreover, due to rapid industrialization, a large amount of waste and by-products are generated whose proper disposal otherwise remains a problem. In this domain, the concept of green concrete is becoming popular to condense the adverse effects of concrete and effectively utilize the waste products. Utilization of alternative materials as a cement replacement Sustainability 2021, 13, 2867 3 of 20 certain shortcomings due to its incapability to consider a few deviating datasets. For better modeling output, the deviating datasets need to be removed from the data [5].
Considering the drawbacks of several machine learning models, the current study employed an advanced algorithm, i.e., multiple expression programming (MEP), to model the mechanical properties of SCBA concrete. The MEP has the ability to encode multiple chromosomes to a single program. The best possible chromosome can be selected based on evaluating the fitness [35]. MEP is considered as an improved version of the GEP having the capacity to estimate accurate results, given the complexity of the target, as compared to other modeling techniques [36]. It is capable of highlighting and excluding errors from the final expression form and also pre-specifying of final expression is not compulsory. MEP has been rarely used in the civil engineering field despite its distinguished features and could be an effective tool for solving complex problems. In this study, the mechanical properties of SCBA concrete, i.e., compressive strength (CS), splitting tensile strength (ST), and flexural strength (FS), were modeled employing MEP with the most influential parameters as modeling inputs. A large dataset was collected from the literature. The robustness of the developed MEP was validated using the mechanical properties derived from laboratory testing. Furthermore, the behavior of the models, during both training and testing, was judged using various statistical indicators and the modeling output was cross-validated through a 10-fold method. Parametric analysis was performed to find out the effect and most influential input parameters.

Multi Expression Programming
John Holland (1975) was the first one who introduced the genetic algorithm (GA) motivated by the evolutionary theory. Similarly, the genetic programming (GP) was initially proposed by Nichael Cramer (1985) [37][38][39]. The GP is mainly an extended form of GA which can be distinguished by description of the solution. Various genetic operators such as crossover, reproduction, and mutation are brought together in GP with little variation [36,40]. An arbitrary population of the parse trees is created to increase the diversity of the model [41]. Based on fitness criteria, some portion of the trees are killed in the implementation stage. Furthermore, the best individuals are selected as a parent to create a new generation of individuals based on the error assessment criteria. This cyclic process of selecting the best individuals continues when termination criteria is reached or a specified number of iterations are attained [38]. However, the GP algorithm is accompanied by certain shortcomings. Firstly, the tree crossover genetic operators are utilized which creates a huge population of parse trees, resulting in an increased simulation time and space requirement [38]. Secondly, the non-linear structure of GP is based on phenotype and genotype, thus making it difficult to develop practical empirical equations for the output [42]. A single expression of pre-specified length is implemented in the GP process regardless of the problem complexity. Moreover, the process escalates with the number of genes in chromosomes to a certain value and it is practically impossible to enforce a complex chromosome to produce simple expression [35].
Recently, a linear variant of machine learning, i.e., multi expression programming (MEP), has been proposed with individual entities as a variable length [35,43]. The linear variants permit MEP to extricate the genotype and phenotype [44]. This technique is considered advantageous over other techniques due to linear chromosomes and permitting multiple solutions in an individual chromosome. This unique feature allows searching in a wider space to gather the best possible solution. In comparison to the GP technique, the MEP applies simple decoding procedures and is given particular importance in case of unknown complexity of targeted expression [45]. MEP can handle exceptions such as invalid expressions, division by zero. The gene is responsible for generating an exception, then it alters to an arbitrarily terminal symbol; thereby, no infertile individuals enter the next generation. This produces a margin in the structure of chromosomes during the evaluation process [35]. The methodology flowchart of MEP is presented in Figure 1. Which Sustainability 2021, 13, 2867 4 of 20 is adapted from [46]. The MEP is decoded similarly to the C and pascal compiler empirical relationship to machine coding. The outcome of the MEP process is a combination of mathematical operators or variables in a linear string of instructions [46]. The chromosome length is governed by several genes per chromosome, whereas the gene encoded the elements in function and terminal set. exceptions such as invalid expressions, division by zero. The gene is responsible for generating an exception, then it alters to an arbitrarily terminal symbol; thereby, no infertile individuals enter the next generation. This produces a margin in the structure of chromosomes during the evaluation process [35]. The methodology flowchart of MEP is presented in Figure 1. Which is adapted from [46]. The MEP is decoded similarly to the C and pascal compiler empirical relationship to machine coding. The outcome of the MEP process is a combination of mathematical operators or variables in a linear string of instructions [46]. The chromosome length is governed by several genes per chromosome, whereas the gene encoded the elements in function and terminal set. In order to develop a generalized relationship, the MEP model requires several fitting parameters. In this study, basic mathematical operators were considered to get simple expressions and a trial and error method was used to get the fitting parameters [47]. Population size was used to specify the number of programs required as the high population size can lead to complex and long convergence times, thus often causing overfitting problems. Moreover, the error can be reduced by running algorithms for an enormous number of generations. Several combinations of parameters were initiated, and the best possible grouping was selected based on model performance. The parameter setting for MEP modeling is presented in Table 1. In order to develop a generalized relationship, the MEP model requires several fitting parameters. In this study, basic mathematical operators were considered to get simple expressions and a trial and error method was used to get the fitting parameters [47]. Population size was used to specify the number of programs required as the high population size can lead to complex and long convergence times, thus often causing over-fitting problems. Moreover, the error can be reduced by running algorithms for an enormous number of generations. Several combinations of parameters were initiated, and the best possible grouping was selected based on model performance. The parameter setting for MEP modeling is presented in Table 1. A large number of studies have been conducted during the last decade that has focused on neural network and GEP models. The abovementioned advantages of MEP over other methods can lead to accurate and reliable models in many fields. The MEP has been applied in a few research studies to estimate elastic modulus of high and normal strength concrete [44], formulation of compressive strength of Portland cement [48], soil deformation modulus [49], consolidation depth of soil layer [46], and strength of polymer confined concrete columns [36]. In the present study, the SCBA concrete properties have been considered to evaluate the performance of the MEP.

Modeling Dataset
For this study, a comprehensive dataset of SCBA concrete properties was obtained from the published literature [4,[11][12][13]15,16, and was analyzed to get the most effective parameters affecting the properties of SCBA concrete. Literature study revealed that the mechanical properties of normal concrete are influenced by mix composition. For concrete with admixtures and supplementary cementitious materials, the type and dosage of the admixtures affect the properties of concrete [72]. Therefore, in this study, the water to cement ratio (w/c), the proportion of bagasse ash (SCBA%), fine aggregate content (FA), cement content (CC), and quantity of coarse aggregate (CA) are chosen as input for the prediction of SCBA concrete mechanical properties. Frequency histograms of the modeling input parameters are shown in Figure 2 to visualize the distribution [73]. The statistical parameters, i.e., mean, minimum, maximum, and standard deviation of input data used are also given in Table 2. A large number of studies have been conducted during the last decade that has focused on neural network and GEP models. The abovementioned advantages of MEP over other methods can lead to accurate and reliable models in many fields. The MEP has been applied in a few research studies to estimate elastic modulus of high and normal strength concrete [44], formulation of compressive strength of Portland cement [48], soil deformation modulus [49], consolidation depth of soil layer [46], and strength of polymer confined concrete columns [36]. In the present study, the SCBA concrete properties have been considered to evaluate the performance of the MEP.

Modeling Dataset
For this study, a comprehensive dataset of SCBA concrete properties was obtained from the published literature [4,[11][12][13]15,16, and was analyzed to get the most effective parameters affecting the properties of SCBA concrete. Literature study revealed that the mechanical properties of normal concrete are influenced by mix composition. For concrete with admixtures and supplementary cementitious materials, the type and dosage of the admixtures affect the properties of concrete [72]. Therefore, in this study, the water to cement ratio (w/c), the proportion of bagasse ash (SCBA %), fine aggregate content (FA), cement content (CC), and quantity of coarse aggregate (CA) are chosen as input for the prediction of SCBA concrete mechanical properties. Frequency histograms of the modeling input parameters are shown in Figure 2 to visualize the distribution [73]. The statistical parameters, i.e., mean, minimum, maximum, and standard deviation of input data used are also given in Table 2.

Cross-Validation with k-Fold Algorithm
Evaluating the performance of machine learning models can be a difficult task since the learning models usually fail to produce generalizable results based on data it has not been trained on. Usually, the data are divided into train and test set for models training and testing, respectively, and performance is then evaluated with statistical criteria. However, this method does not go well if a larger dataset is not available. Moreover, overfitting of the data is another problem in machine learning-based modeling. To ensure that the model can perform well on unseen data and to address the modeling overfitting, a resampling technique is used, called k-fold cross-validation. This method distributes the actual data to k subclasses [74]. The enhanced performance and efficiency of the 10-fold cross-validation method is reported in the literature [75]. In the current study, the 10-fold cross-validation is adopted by randomly dividing the dataset into ten subsets. Among all the ten subsets, each subset is used for validation to examine the grouping model and the same method is repeated for all the remaining subsets. The accuracy of the final model is then expressed in terms of mean accuracy obtained by the ten-fold method in ten individual rounds.

Performance Measures
The developed models are evaluated employing statistical checks, i.e., correlation coefficient (R), Nash Sutcliff efficiency (NSE), mean absolute error (MAE), relative root mean squared error (RRMSE), relative squared error (RSE), performance index (ρ), and root means squared error (RMSE) [44,76]. The high values of R and NSE and low values of RMSE and MAE indicate better performance. The indicator R quantifies the linear relationship between input and output [30] and an R-value of more than 0.8 signifies an excellent correlation between predicted and actual data [77]. Despotovic et al., (2016) [78] reported that a model can be considered excellent if the value of RRMSE is between 0 and 0.10; good if the values are between 0.11 and 0.20. The expressions for the statistical indicators are shown below (Equations (1)- (7)).
where n, P i , M i , P i , and M i , show a number of data points, predicted data, measured data, mean of predicted data, and mean of measured data, respectively. The T and TE represent the training and testing datasets, respectively.

Sugarcane Bagasse Ash (SCBA) Characterization
The waste SCBA was obtained from sugar industry located in Khyber Pakhtunkhwa, Pakistan. It has been reported in the literature that fibrous and unburnt content in waste bagasse ash reduces the pozzolanic activity. Removal of unburnt particles, proper characterization, and grinding up to cement fineness considerably increases its pozzolanic activity [4,79]. Therefore, the obtained SCBA was sieved from the #200 sieve and afterward grinded to reduce the particle size. Grinding was carried out in a ball mill machine and a grinding media (ceramic balls)-to-ash ratio of 5 by weight was kept constant. The grinding was done for different time periods, i.e., 15,30,45, and 60 min. The surface area of SCBA was determined according to ASTM C204 for each successive grinding period. It can be seen in Figure 3 that increasing the duration to 60 min, the specimen SCBA gives more surface area compared to cement.
where , , , , and , show a number of data points, predicted data, measured data, mean of predicted data, and mean of measured data, respectively. The T and TE represent the training and testing datasets, respectively.

Sugarcane Bagasse Ash (SCBA) Characterization
The waste SCBA was obtained from sugar industry located in Khyber Pakhtunkhwa, Pakistan. It has been reported in the literature that fibrous and unburnt content in waste bagasse ash reduces the pozzolanic activity. Removal of unburnt particles, proper characterization, and grinding up to cement fineness considerably increases its pozzolanic activity [4,79]. Therefore, the obtained SCBA was sieved from the #200 sieve and afterward grinded to reduce the particle size. Grinding was carried out in a ball mill machine and a grinding media (ceramic balls)-to-ash ratio of 5 by weight was kept constant. The grinding was done for different time periods, i.e., 15,30,45, and 60 min. The surface area of SCBA was determined according to ASTM C204 for each successive grinding period. It can be seen in Figure 3 that increasing the duration to 60 min, the specimen SCBA gives more surface area compared to cement. The chemical composition and amorphous nature of SCBA was evaluated through X-ray fluorescence (XRF) and X-ray diffraction (XRD). The chemical composition of SCBA is given in Table 3. According to the available literature [80,81], the pozzolanic activity mainly depends on the silica and alumina content in a pozzolanic material. It can be inferred from XRF results that the quantity of silica, alumina, and iron oxide is above 70%, which qualifies the ASTM C618-05 standard for pozzolan. The chemical composition and amorphous nature of SCBA was evaluated through X-ray fluorescence (XRF) and X-ray diffraction (XRD). The chemical composition of SCBA is given in Table 3. According to the available literature [80,81], the pozzolanic activity mainly depends on the silica and alumina content in a pozzolanic material. It can be inferred from XRF results that the quantity of silica, alumina, and iron oxide is above 70%, which qualifies the ASTM C618-05 standard for pozzolan. The XRD pattern of SCBA exhibits the amorphous nature of ash with few pronounced peaks. The SCBA was already processed and fibrous and unburnt contents were removed which enhances its properties. One peak at 2θ = 26 • was observed and is termed as cristobalite [4]. This is formed due to high and long temperature burning and is classified as a crystalline polymorph of silica. XRD pattern of SCBA is shown in Figure 4.  The XRD pattern of SCBA exhibits the amorphous nature of ash with few pronounced peaks. The SCBA was already processed and fibrous and unburnt contents were removed which enhances its properties. One peak at 2θ = 26° was observed and is termed as cristobalite [4]. This is formed due to high and long temperature burning and is classified as a crystalline polymorph of silica. XRD pattern of SCBA is shown in Figure 4. Scanning electron micrographs of SCBA are shown in Figure 5. X500 to X5000. In the micrographs, the heterogeneous nature of SCBA is prominent. The SCBA exhibited various shape particles such as elongated, needle, flat, oval, and irregular. The size of needle and oval shape particles are found up to 50 μm. Likewise, 10 μm size of irregular particles can also be observed. The irregular shape particles are mainly silica-rich [82]. Furthermore, voids over the surface are also visible in the form of black spots and size ranges from 5 to 10 μm. The voids in the SCBA particles may be due to excessive grinding. Overall, the size of different particles ranges from 5 to 50 μm which can be observed in Figure  5. Scanning electron micrographs of SCBA are shown in Figure 5. X500 to X5000. In the micrographs, the heterogeneous nature of SCBA is prominent. The SCBA exhibited various shape particles such as elongated, needle, flat, oval, and irregular. The size of needle and oval shape particles are found up to 50 µm. Likewise, 10 µm size of irregular particles can also be observed. The irregular shape particles are mainly silica-rich [82]. Furthermore, voids over the surface are also visible in the form of black spots and size ranges from 5 to 10 µm. The voids in the SCBA particles may be due to excessive grinding. Overall, the size of different particles ranges from 5 to 50 µm which can be observed in Figure 5.

Mix Proportions and Specimen Preparation
Experimental testing of SCBA concrete was conducted to verify the behavior of the model through the validation requirement. SCBA was incorporated in concrete as 0, 10, 20, 30, and 40% by replacing cement. Concrete specimens were casted, and comparison of hardened properties was made between SCBA modified matrix (BC) and control specimens (CM). The water-to-cement ratio for all the specimens was kept constant to compare the modified matrix with a non-modified matrix. The detailed mix design with the formulation is depicted in Table 4. Concrete specimens were prepared and tested for compressive strength (CS), splitting tensile strength (ST), and flexure strength (FS) at the curing age of 28 days in accordance with ASTM C39, ASTM C496, and ASTM C293, respectively.
The results obtained from laboratory tested concrete were used for validation of the MEP models.

Mechanical Properties of SCBA Concrete
Concrete specimens were prepared in the laboratory with 0, 10, 20, 30 and 40% cement replaced with SCBA to determine the compressive, splitting tensile and flexural

Mix Proportions and Specimen Preparation
Experimental testing of SCBA concrete was conducted to verify the behavior of the model through the validation requirement. SCBA was incorporated in concrete as 0, 10, 20, 30, and 40% by replacing cement. Concrete specimens were casted, and comparison of hardened properties was made between SCBA modified matrix (BC) and control specimens (CM). The water-to-cement ratio for all the specimens was kept constant to compare the modified matrix with a non-modified matrix. The detailed mix design with the formulation is depicted in Table 4. Concrete specimens were prepared and tested for compressive strength (CS), splitting tensile strength (ST), and flexure strength (FS) at the curing age of 28 days in accordance with ASTM C39, ASTM C496, and ASTM C293, respectively. The results obtained from laboratory tested concrete were used for validation of the MEP models.

Mechanical Properties of SCBA Concrete
Concrete specimens were prepared in the laboratory with 0, 10, 20, 30 and 40% cement replaced with SCBA to determine the compressive, splitting tensile and flexural strength. The results of mechanical properties are given in Table 5. It can be observed ( Table 5) that the strength of modified concrete (BC) increases up to 10% SCBA replacement and consistently decreases afterward. The maximum strength gained is at 10% SCBA replacement and may be due to the fine SCBA particles dispersed all over the concrete mix. Moreover, the silica content reacts with lime (resulted from cement hydration) and produces more calcium silicate hydrate (CSH) [53,83]. The decrease in strength as compared to CM is 6.5, 17.3, and 30.3% for 20BC, 30BC, and 40BC, respectively. This may be attributed to non-availability of adequate calcium hydroxide, which consequently reduced the formation of secondary CSH. Similar results were observed for ST and FS as shown in Table 5. In comparison with CM, the increase in ST is 25.3 and 15.8% for 10BC and 20BC, respectively. Moreover, the decrease in ST of SCBA concrete is 7.9 and 23.8% for 30BC and 40BC, respectively. For the case of FS, the maximum strength is achieved by 10BC concrete. When compared to CM, the reduction in FS for 20BC, 30BC and 40BC is 15.5, 28.8 and 42.5%, respectively. The decrease in ST and FS for 30% and 40% replacement level maybe due to less formation of CSH gel and calcium aluminate. The maximum strength gain by 10BC may be due to the microfibrous nature of SCBA which is related to secondary CSH formation. In addition, the formation of aluminates take place which interlocks between hydrated pastes and directly contributed to ST and FS [84,85]. In a nutshell, the improved properties of SCBA concrete may be due to the formation of more hydrated products, improved interfacial transition zone, and filling effect [84,86,87].

Formulation of Compressive Strength (CS)
The output of MEP model (based on five input variables) for compressive strength (CS) is given in Figure 6 for models training, testing and validation phase. The modeling results were decoded to get the empirical expressions for the prediction of the CS. The developed equation is shown Equation (8) and can be used to predict the 28 days CS. Figure 6 graphically demonstrates the comparison of model predicted and actual CS for training, testing, and validation. The expressions for regression line are also shown in the graph. It is known that the slope of the line should be nearly equal to 1 for an ideal situation. It can be deduced that the model took the effect of the inputs and retains a strong relationship as obvious from the slope, i.e., 0.8951, 0.9315, and 0.9014 for training, testing, and validation, respectively. The statistical indicators for CS model are shown in Table 6. The value of indicator, i.e., R, is 0.91 and 0.94 during training and testing, respectively. The RMSE and MAE are considerably low which highlighted the high accuracy of the model. From statistical analysis, it can be also inferred that the model retains a high generalized relationship.

Formulation of Splitting Tensile Strength (ST)
The splitting tensile strength (ST) of SCBA concrete at 28 days curing age can be predicted using Equation (9). A comparative assessment graph between MEP prediction and actual ST is presented in Figure 7. An excellent correlation exists among experimental and predicted data to estimate ST. The respective slope of the regression line is close to the ideal fit, i.e., 0.9351, 0.8903, and 0.9273 for training, testing, and validation, respectively. The model for ST performs exceptionally well on the training dataset. Table 6

Formulation of Splitting Tensile Strength (ST)
The splitting tensile strength (ST) of SCBA concrete at 28 days curing age can be predicted using Equation (9). A comparative assessment graph between MEP prediction and actual ST is presented in Figure 7. An excellent correlation exists among experimental and predicted data to estimate ST. The respective slope of the regression line is close to the ideal fit, i.e., 0.9351, 0.8903, and 0.9273 for training, testing, and validation, respectively. The model for ST performs exceptionally well on the training dataset. Table 6 shows the performance measures of the developed MEP model for ST. The NSE and R values for ST model are above 0.9 for model training as well as testing. The minimum and maximum value of RRMSE is 0.09 and 0.10, respectively. The performance index (ρ) is nearly equal to zero for ST model, hence revealing that the issue of over-fitting by the model has also been reduced to a large extent. It has been reported in the literature [88] that the performance of a model depends on the number of data points incorporated in the modeling process. In this study, 110 ST data points have been selected to model the tensile strength. Therefore, an accurate model with minimum error has been achieved.
Sustainability 2021, 13, 2867 12 of 21 model are above 0.9 for model training as well as testing. The minimum and maximum value of RRMSE is 0.09 and 0.10, respectively. The performance index (ρ) is nearly equal to zero for ST model, hence revealing that the issue of over-fitting by the model has also been reduced to a large extent. It has been reported in the literature [88] that the performance of a model depends on the number of data points incorporated in the modeling process. In this study, 110 ST data points have been selected to model the tensile strength.
Therefore, an accurate model with minimum error has been achieved. where

Formulation of Flexural Strength (FS)
The MEP modeling output for flexural strength (FS) is graphically presented in Figure 8. The slopes of the regression line are observed to be 0.9494, 0.9026, and 0.9332 for training, testing and validation data, respectively. FS can also be predicted using Equation (10) derived by MEP. Comparing with the models for compressive and tensile strength, an excellent estimate is observed between experimental and MEP model predicted data for FS as illustrated in Figure 8. The MEP model for FS also performs exceptionally well as illustrated in Table 6

Models Error Assessment
To infer the absolute error between the model predicted and actual data, the data points are plotted in Figure 9 showing the absolute error.

Models Error Assessment
To infer the absolute error between the model predicted and actual data, the data points are plotted in Figure 9 showing the absolute error.

Model Cross-Validation Results
The results of cross-validation are graphically shown in Figures 10-12

Model Cross-Validation Results
The results of cross-validation are graphically shown in Figures 10-12

Parametric Analysis
In machine learning-based modeling, it is essential to perform several checks and analyses to ensure that the models have considered the effect of all the physical phenomena and inputs. Parametric analysis has been recommended in literature studies to evaluate the robustness and accuracy of the final models. In the current study, the same method has been used by plotting the variation of the model output with a single input variable. The same procedure is applied to the individual input parameter. Figure 13 illustrates the parametric study results for the developed compressive strength (CS) model. The same patterns were observed for ST and FS as well. Therefore, parametric study results are discussed in detail for CS only in the below section.

Parametric Analysis
In machine learning-based modeling, it is essential to perform several checks and analyses to ensure that the models have considered the effect of all the physical phenomena and inputs. Parametric analysis has been recommended in literature studies to evaluate the robustness and accuracy of the final models. In the current study, the same method has been used by plotting the variation of the model output with a single input variable. The same procedure is applied to the individual input parameter. Figure 13 illustrates the parametric study results for the developed compressive strength (CS) model. The same patterns were observed for ST and FS as well. Therefore, parametric study results are discussed in detail for CS only in the below section. ena and inputs. Parametric analysis has been recommended in literature studies to evaluate the robustness and accuracy of the final models. In the current study, the same method has been used by plotting the variation of the model output with a single input variable. The same procedure is applied to the individual input parameter. Figure 13 illustrates the parametric study results for the developed compressive strength (CS) model. The same patterns were observed for ST and FS as well. Therefore, parametric study results are discussed in detail for CS only in the below section.  An increase in w/c decreases the strength parameters of concrete [88,89]. It can be observed in Figure 13a that the increased w/c reduced the concrete strength. The decrease in CS with increased w/c may be attributed to the dilution of the cement matrix. Moreover, the available literature revealed that the properties of SCBA concrete depend on its oxide composition and SCBA replacement percentage [4]. Figure 13b shows that the compressive strength of SCBA concrete increases with an increase in SCBA replacement up to 10%. Afterward, the strength decreases for higher replacement level. Maximum compressive strength has been achieved at 10% SCBA replacement. The same results were obtained during experimental analysis of SCBA concrete where maximum strength gain was attained by 10BC. Therefore, the output of the parametric study is parallel with laboratory testing and same findings were reported in the literature [4,12,13]. Figure 13c shows the variation of compressive strength with increasing cement content. The increasing trend can be observed as the higher cement content produces more calcium silicate hydrate which forms a dense structure, thereby directly increasing the compressive strength. Figure 13d,e illustrate the variation of compressive strength with fine and coarse aggregate content, respectively. Both the figures demonstrate that CS decreases with an increase in FA and CA content. As the aggregates are inert materials and used as a filler in concrete production, when keeping all the parameters constant, the properties of concrete decrease with the increase in aggregate quantity.

Conclusions
The present study conducted experimental testing on sugarcane bagasse ash (SCBA)incorporated sustainable concrete and adopted a machine learning technique to model the strength properties of SCBA concrete. The specific conclusions of this study are as follows.

1.
The SCBA showed good pozzolanic properties when processed, i.e., passed from sieve #200 and grinded up to cement fineness.

2.
Microfibrous structure and irregular shape particles were observed in the SEM images of processed SCBA.

3.
The concrete showed maximum strength when cement was replaced with 10% SCBA. Afterward, strength reduction was observed for higher replacement levels.

4.
The multi expression programming (MEP) was found to be very efficient in modeling the strength properties of SCBA concrete. The parametric study showed that the developed MEP models for SCBA concrete are accurate and revealed the effect of input parameters in the modeling output.
Moreover, the work presented in the current study has certain shortcomings. The key focus of this study was to examine the effect of concrete constituents on the mechanical properties of SCBA concrete. Indeed, other important factors also need to be investigated that are important to mechanical viewpoints such as curing conditions, type of cement, reactivity, type of ash, and testing conditions. It is recommended that further research should be conducted with an extensive dataset for model training and testing. Moreover, some deep learning techniques, i.e., convolutional neural network, neuro-fuzzy system, and ensemble modeling, should be considered for comparative analysis and assessment of concrete properties.