Predictive Modeling of Compression Strength of Waste PET/SCM Blended Cementitious Grout Using Gene Expression Programming

<jats:p>The central aim of this study is to evaluate the effect of polyethylene terephthalate (PET) alongside two supplementary cementitious materials (SCMs)—i.e., fly ash (FA) and silica fume (SF)—on the 28-day compressive strength (CS28d) of cementitious grouts by using. For the gene expression programming (GEP) approach, a total of 156 samples were prepared in the laboratory using variable percentages of PET and SCM (0–10%, each). To achieve the best hyper parameter setting of the optimized GEP model, 10 trials were undertaken by varying the genetic parameters while observing the models’ performance in terms of statistical indices, i.e., correlation coefficient (R), root mean squared error (RMSE), mean absolute error (MAE), comparison of regression slopes, and predicted to experimental ratios (ρ). Sensitivity analysis and parametric study were performed on the best GEP model (obtained at; chromosomes = 50, head size = 9, and genes = 3) to evaluate the effect of contributing input parameters. The sensitivity analysis showed that: CS7d (30.47%) &gt; CS1d (28.89%) &gt; SCM (18.88%) &gt; Flow (18.53%) &gt; PET (3.23%). The finally selected GEP model exhibited optimal statistical indices (R = 0.977 and 0.975, RMSE = 2.423 and 2.531, MAE = 1.918 and 2.055) for training and validation datasets, respectively. The role of PET/SCM has no negative influence on the CS28d of cementitious grouts, which renders the PET a suitable alternative toward achieving sustainable and green concrete. Hence, the simple mathematical expression of GEP is efficacious, which leads to saving time and reducing labor costs of testing in civil engineering projects.</jats:p>

section_content: Plastic production and increased usage have led to significant environmental damage. Plastic waste generation is increasing rapidly with increases in the global population, and it is predicted that this worldwide generation may be doubled by 2050 [1]. An exponential increase in the production of plastics has been reported such that, during 1950 and 2015, it increased from 2. 3 million tons to 448 million tons [2]. Another study showed that about half of all plastic has been produced in the last 15 years. Of this plastic waste, approximately 8 million tons have been recorded as washing into the oceans every year [2]. It has been estimated that more than 8. 3 billion tons of plastic has been manufactured since 1950s and around 60% of that plastic has ended up in landfills and the natural environment [3]. Therefore, to cope with this global environmental issue, it is immensely pertinent to explore the recycling of waste plastics. In this scenario, replacement of virgin materials used in construction industries with waste materials is desideratum [4]. One such alternative is recycling and reusing waste plastic in the construction of roads/highways. Waste plastic has widely been used as a substitute for bitumen or as an aggregate replacement in asphalt concrete pavements. Significant improvement in the performance properties of binder and asphalt mixtures has been reported previously [5, 6]. A significant contribution to the economy and environment can be attained by recycling single-use waste plastic on a wide scale. 

Recycling of waste polyethylene terephthalate (PET) in cementitious materials and its effect on mechanical properties have also been investigated in a few studies. In the concrete industry, waste plastics are used as a fiber or replacement of fine and coarse aggregates [7]. They have been used to replace aggregates in order to produce lightweight concrete elements [8]. Studies have shown an increase in the tensile strength and crack resistance of concrete by inclusion of recycled PET as a fiber in the concrete mixture [9, 10]. However, the incorporation of PET as an aggregate or sand replacement with a volume of more than 10% causes a significant reduction in compressive strength compared with that of normal strength concrete [11, 12]. Moreover, partial replacement of cement with plastic waste has also been studied, where significant reduction in the compressive strength was noticed. The reduction in compressive strength may reach up to 23% and 72% at 5% and 20% replacement of the cement, respectively [13]. Schaefer et al. [14] also used waste PET as a partial replacement of cement in mortars and observed reduction in compressive strength. However, using gamma irradiated PET instead of regular PET and additional replacement of FA or SF in cement mortar resulted in improved strength properties [14]. Some other studies also reported that the addition of gamma irradiated waste plastic ameliorated the strength properties of cementitious grouts [15] [16] [17] [18]. 

The utilization of natural resources can be reduced while using waste plastic instead of aggregates or sand, which in return reduces environmental pollution and the cost of construction alongside contributing to better engineering properties of asphalt concrete [19]. Moreover, recycling waste plastic would reduce the burden of landfilling and save marine life [20]. Furthermore, industrial wastes and/or byproducts (such as FA and SF) are increasing due to increase in populations and industrialization all over the world. Due to shortage of land, many countries are facing problems of waste disposal which is a serious concern affecting the ecosystem. Therefore, reuse and recycling of these wastes and/or byproducts in the construction industries could be a viable alternative [21]. The perpetual and rapid utilization of natural resources is polluting the environment and deteriorating the surrounding environment [22]. FA and SF are categorized as supplementary cementing materials (SCM) owing to their pozzolanic properties, and therefore they contribute to the improvement in performance properties of cement concrete. The use of materials or mineral admixtures to replace cement in mortar and concrete is almost unavoidable in the effort to achieve sustainability, superior performance, and financial benefits [23, 24]. Therefore, the recycling of municipal wastes (such as glass, plastic, rubber, wood, etc. ) and industrial wastes/byproducts (such as FA, SF, ground granulated blast-furnace slag, etc. ) as replacements for cement, sand, or aggregates could be beneficial to the environment in terms of reducing the usage of non-renewable natural resources and to the construction industry in terms of cost and enhancing the properties of concrete [21, 25]. 

Artificial Intelligence (AI) approaches are becoming increasingly popular as a result of their improved predictive capability opposed to earlier techniques, and they are being utilized to simulate the complicated behavior of a range of structural engineering problems [26]. Data mining for processes in chemistry, materials, and especially civil engineering, has been frequently documented since 2000, which is largely attributed to the rapid growth of machine learning (ML) algorithms [27] [28] [29]. In the past few years, commonly deployed AI techniques employing conventional statistical methods in civil engineering include, artificial neural networks (ANNs) [30], genetic algorithm (GA) [31], genetic expression programming (GEP) [32], multi expression programming (MEP) [33], support vector machines (SVMs) [34], alternate decision trees (ADTs) [35], ensemble random forest (RF) regression [36], etc. Note that Giustolisi and Doglioni [37] advocated using white, black, and grey colors to categorize distinct mathematical models. The known variables and parameters of the first kind-the white-box model-are based on physical laws, which build correct physical connections and provide maximum transparency. Second, black-box models rely on regressive data-driven systems with unknown functional forms of connections among variables that must be approximated (e. g., ANN, ANFIS, etc. ). Finally, grey-box models represent logical systems where a mathematical framework better assesses the behavior of the system (e. g., GEP, MEP, etc. ). GEP can be classified as a 'grey-box model' because of its conceptualization of physical phenomena in a symbolic and easy way [38] [39] [40]. The GEP-based models are found to yield better results in contrast to other genetic-based and ANN-based approaches [29, 41, 42]. It is also associated with the fact that a typical GEP chromosome consists of head and tail that comprise particular symbols, thus providing a better way to encode syntactically correct computer programs as compared with the MEP approach [43]. In the past decade, GEP has been broadly and efficaciously used in addressing specific structural engineering issues, for instance, evaluating the compression strength of concrete incorporating (fly ash admixture, geopolymer, eco-efficient GGBSbased geopolymer, nanomaterials, etc. ), the splitting tensile strength from the compressive strength [44], and the post-fire compressive strength of recycled PET aggregate concrete reinforced with steel fibers [45], among others. The GEP approach is widely deployed in solving regressions, modeling functions, predicting, detecting, and data mining [46, 47]. In addition, the formulated genetic programming models are efficient due to absence of the assumption of fixed connection for developing a GEP model [48]. 

Furthermore, Ferreira and Jalali [49] predicted the early age compressive strength using a technique based on the nondestructive testing data. According to [50], models based on multi-layer feed-forward neural networks (MFNNs) may be built to forecast concrete's 28-day compression strength on the basis of its selective influencing parameters. Rafi and Nasir [51] suggested an analytical approach for predicting concrete's 28-day strength (CS 28d ) based on its 7-day strength (CS 7d ) by proposing a mathematical equation which was further expanded to forecast the strength of concrete constructed using Pakistani cements. To date, no attempts were made by previous researchers to formulate empirical models using the GEP approach for forecasting the 28-day compression strength (CS 28d ) of waste PET/SCM blended cementitious grout despite of their practical significance. With these uncertainties in consideration, the present research incorporates the GEP method to compute the compression strength of waste PET/SCM blended cementitious grout to compare the behavior of a variety of formulated GEP models.
section_name: Introduction

section_name: Materials and Methods

section_content: Materials used in the experimental program were collected from local vendors, which include ordinary Portland cement, waste PET, SF, FA, and superplasticizer. Waste PET (particle size less than 150 µm) was used as a replacement of cement. The waste PET was obtained from a plastics factory, where PET flakes are recycled into plastic-based products. This process of recycling PET results in a sufficient amount of powder being discarded for disposal. For use in cementitious grout, the PET powder obtained was sieved so that the particle size was less than 150 microns. Table 1 shows the basic characteristics of waste PET.
section_name: Materials

section_content: The mixing of cementitious grouts were performed in accordance with ASTM C305 [52]. The required quantity of cement, waste PET, and SCMs (i. e., FA and SF) were initially dry mixed followed by adding two-thirds water and further mixing. The remaining water and superplasticizer were added and mixed to ensure the homogeneity of cement grouts. The mixed proportion of grouts containing PET and SCM are presented in Table 2. The fresh cement grout was checked for the flow value using Malaysian flow cone apparatus [53]. According to specifications, 11-16 s are required for 1 L of fresh grout to flow-out of the cone [53]. The purpose of achieving such flow value (or flowability) is required to ensure proper filling of voids in an open-graded asphalt mixture. Moreover, mold (50 × 50 × 50 m 3 ) were also prepared from all combination of grouts and tested for evaluating the compressive strength. After demolding, the specimens were kept in water bath for curing until test date. On completion of curing period (i. e., 1, 7, and 28 d), specimens were tested using a universal testing machine (UTM) with a capacity 3000 kN and at a rate of 0. 90 kN/second for the determination of compressive strength (CS 1d, CS 7d, CS 28d, respectively) of hardened cement grouts following ASTM C109 [54]. The values of the input and output parameters are given in Table S1 (Supplementary Materials).
section_name: Cement Grout Preparation and Testing Methods

section_content: A GEP approach based on Darwin's evolution theory and Mendel's genetic theory is one of the most intellectually compelling computational intelligence approaches [55, 56]. With GEP, there are two languages: the language of genes and the language of expression trees (ETs), and understanding the sequence or structure of one is equivalent to understanding the other [57]. A typical GEP model involves the following steps, as shown in Figure 1a. a. 

In GEP modeling, chromosome numbers are randomly generated for designated numbers, and the Karva language (which represents symbols) is used for introducing the chromosomes. The chromosomes and genes typically consist of a head and a tail. A head consists of either a function or a set of terminal symbols, while a tail consists of only terminal symbols. A model's head size depends on the complexity of each parameter, whereas the number of genes determines the number of sub-ETs. 

b. 

The length of chromosomes is fixed, so the chromosomes can be transformed into an algebraic expression [58] as illustrated in Figure 1b. GEP genes contain a fixed list of terms, so each term represents a function; for example, arithmetic operations (+, -, ×, /), Boolean logic functions (AND, OR, NOT, etc. ), mathematical functions (cos, sin, ln), conditional functions (IF, THEN, ELSE), or other function categories [59]. c. 

The chromosomes are then represented by ETs, which come in diverse shapes and sizes. After that, primary genetic operators-such as crossover, mutation, transposition, and gene recombination (1-point, 2-point, and gene recombination)-are applied to each chromosome according to their ratios [60]. Figure 1b depicts a typical expression tree (ET), illustrating how crossovers and mutations work. In Equation ( 1 ), we can also see that the ET is expressed by means of a Karva notation or K-expression (Figure 1b ) [61]. 

Once the stopping condition (highest number of generations or best solution) is reached, the process is complete [62]. 

d. When the maximum number of iterations or favorite fitness value is not achieved, the roulette wheel method (as well as greedy over the selection, ranking selection, tournament selection, and elite selection [27]) is adopted, which selects the viable chromosomes of first-generation for continuation to the next. Herein, the process is rewound for a specified number of generations or until the selection of the optimal solution is made [63]. 

b. The length of chromosomes is fixed, so the chromosomes can be transformed into an algebraic expression [58] as illustrated in Figure 1b. GEP genes contain a fixed list of terms, so each term represents a function; for example, arithmetic operations (+, -, ×, /), Boolean logic functions (AND, OR, NOT, etc. ), mathematical functions (cos, sin, ln), conditional functions (IF, THEN, ELSE), or other function categories [59]. c. The chromosomes are then represented by ETs, which come in diverse shapes and sizes. After that, primary genetic operators-such as crossover, mutation, transposition, and gene recombination (1-point, 2-point, and gene recombination)-are applied to each chromosome according to their ratios [60]. Figure 1b depicts a typical expression tree (ET), illustrating how crossovers and mutations work. In Equation ( 1 ), we can also see that the ET is expressed by means of a Karva notation or K-expression (Figure 1b ) [61]. 

Once the stopping condition (highest number of generations or best solution) is reached, the process is complete [62]. d. When the maximum number of iterations or favorite fitness value is not achieved, the roulette wheel method (as well as greedy over the selection, ranking selection, tournament selection, and elite selection [27]) is adopted, which selects the viable chromosomes of first-generation for continuation to the next. Herein, the process is rewound for a specified number of generations or until the selection of the optimal solution is made [63].
section_name: Overview of GEP Method

section_name: Results and Discussion

section_content: While exercising modeling using the GEP model, genetic parameters-namely, num-
section_name: Effect of Genetic Parameters on Performance of Models

section_content: While exercising modeling using the GEP model, genetic parameters-namely, number of chromosomes, head size, and number of genes-were varied alongside the evaluation of models in the form of statistical equations for the correlation coefficient (R), mean absolute error (MAE), and root mean square error (RMSE) given as Equation (2) to Equation ( 4 ) in accordance with the previous research [65] [66] [67] [68] [69] [70]. Initially, number of chromosomes were changed from 30 to 150 to select the best number of chromosomes for the next model generation. The best model for variable setting chromosomes was based on the overall MAE and R 2 for the developed models. Subsequently, head size was changed between 8 and 10. Finally, number of genes was varied from 3 to 5. Further increase in the number of genes may affect the performance of the model; however, it also complexifies the mathematical equation. Figure 2 displays the effect of setting parameters on the performance of the models. Figure 2a -f depict the performance of the models for the training data and validation data with respect to number of chromosomes, head size, and number of genes, respectively. It can be observed that increasing the number of chromosomes from 30 to 50 increased the performance in the sense of increasing R and decreasing MAE and RMSE for both the training and validation data. In contrast, an additional increase in the number of chromosomes from 50 to 100 decreased the performance of the models. The best number of chromosomes was achieved at 50; therefore, more trials on developing the best models by changing head size were exercised with 50 chromosomes. 

Pearson's correlation coefficient, abbreviated as r, is amongst the most widely employed metrics of depicting relationship [71]. Interdependency as well as multicollinearity are better recognized to be examined since they present issues with the understanding of a given AI model [72]. Therefore, the Pearson's correlation matrix was computed for the experimental database of waste PET/SCM blended cementitious grout considered here which comprise five input parameters (PET, SCM, Flow, CS1d, and CS7d) for determination of the output parameter (i. e., CS28d), as shown in Table 3. The correlation matrix is defined with the help of a square symmetrical H × H matrix such that the (uv)th element equals the correlation coefficient Ruv among the (u)th and the (v)th variable. Note that the diagonal members (correlations of the considered input parameters with each other) are always 1 [42, 73]. The r varies from -1 and 1 (0 means zero correlation, while, ±1 depicts strong positive and negative correlation, respectively). The correlation degree follows the order: CS7d > CS1d > PET > SCM > Flow (r = 0. 93, 0. 85, -0. 67, 0. 64, and -0. 52, respectively) which indicates a major contribution of CS7d, CS1d, and PET on the CS28d of waste PET/SCM blended cementitious grout. Since the r values of all parameters are exceedingly higher, they were therefore considered while formulating the GEP models. To see the impact of changing head size on the performance of the developed models, it can be seen that increasing the head size from 8 to 9 increased the values of R from 0. 969 to 0. 977 for the training data; whereas, for the validation data, the value of R remained unchanged for the 8 and 9 head size. Similarly, the values of RMSE and MAE decreased from 3. 267 and 2. 218 MPa to 2. 432 and 1. 918 MPa for the training data. The counterpart values of the validation data somehow increased; however, based on overall value of R, the trials carried out with the 9 head size performed better. Trails 7, 8, and 9 were interpreted with comparatively lower performance. The variation in number of genes from 3 to 5 did not yield the fruitful performance of the model. In summary, the best setting parameters in this comprehensive exercise were obtained as 50, 9, and 3 number of chromosomes, head size, and number of genes, respectively. 

In the above equations, E k and P k represent the actual observations and estimated outcomes, respectively; E k and P k denote the average of the actual observations and estimated outcomes, respectively; and m is the total number of records.
section_name: Effect of Genetic Parameters on Performance of Models

section_name: Prediction Performance of Models

section_content: Five input parameters (PET, SCM, Flow, CS 1d, and CS 7d ) were selected in the current study for determination of the output parameter (i. e., CS 28d ). Although the parameters affecting the strength of cement-based materials are considered to be material mixing ratios-such as the water/binder ratio and the ratio of fine aggregates used-are also desideratum; however, these parameters were kept constant in this study, whereas the aforementioned parameters were varied. The prediction was made based on the varying parameters and this study can be extrapolated in future regarding CS prediction of PET/SCM blended cementitious grout for the design mix used in the study as presented in Section 2. 2. The formulation of CS 28d of waste PET/SCM blended cementitious grout is given by Equation ( 5 ) below. 

Pearson's correlation coefficient, abbreviated as r, is amongst the most widely employed metrics of depicting relationship [71]. Interdependency as well as multicollinearity are better recognized to be examined since they present issues with the understanding of a given AI model [72]. Therefore, the Pearson's correlation matrix was computed for the experimental database of waste PET/SCM blended cementitious grout considered here which comprise five input parameters (PET, SCM, Flow, CS 1d, and CS 7d ) for determination of the output parameter (i. e., CS 28d ), as shown in Table 3. The correlation matrix is defined with the help of a square symmetrical H × H matrix such that the (uv)th element equals the correlation coefficient Ruv among the (u)th and the (v)th variable. Note that the diagonal members (correlations of the considered input parameters with each other) are always 1 [42, 73]. The r varies from -1 and 1 (0 means zero correlation, while, ±1 depicts strong positive and negative correlation, respectively). The correlation degree follows the order: CS 7d > CS 1d > PET > SCM > Flow (r = 0. 93, 0. 85, -0. 67, 0. 64, and -0. 52, respectively) which indicates a major contribution of CS 7d, CS 1d, and PET on the CS 28d of waste PET/SCM blended cementitious grout. Since the r values of all parameters are exceedingly higher, they were therefore considered while formulating the GEP models. The summary of descriptive statistics of the input (PET, SCM, Flow, CS 1d, and CS 7d ) and output parameter (CS 28d ) are given in Table 4. The minimum and maximum limits, standard deviation (SD), kurtosis, as well as skewness values for all these parameters in the considered database have been tabulated. A smaller value of SD suggests that the values are mostly nearing the average (PET, SF, and Flow), while a larger SD represents comparatively higher spread out (CS 1d, CS 7d ) [74]. Skewness (value could be positive, zero, negative, or undefined) helps in determining the magnitude of asymmetry of the probability distribution in case of a real-valued arbitrary parameter from the standpoint of its average value [75]. Furthermore, kurtosis generally ranges between -10 (heavytailed) and + 10 (light-tailed), which helps in the elucidation of the shape of a probability distribution, as described by Brown and Greene [76]. The kurtosis values for PET, SF, CS 1d, and CS 7d are negative and range between -0. 3 and -1. 6 (follow mesokurtic distribution), whereas the only positive value is obtained in case of Flow, i. e., 0. 66 (follow leptokurtic distribution) [77, 78]. Ten trials were run in order to determine the parametric combination with the overall R (average of training and validation values) and MAE value as well as for CS 28d of waste PET/SCM blended cementitious grout (Table 5 ). The details of these undertaken trials for identification of optimal combination of the GEP parameters yielding better performance are presented (Table 4 ). The GEP algorithm is kept running indefinitely so that the correlation, as well as fitness functions (RMSE in this case), do not alter their respective values. Furthermore, while an arithmetic function is selected to connect them, the genes alongside their corresponding head sizes are continuously increased (addition function is taken into consideration here). The process is performed numerous times until the best model is established. Table 6 shows the best settings for the GEP model utilized in this study. The drop recorded in model error is determined by the R value, that is connected to the inclusion of predictor parameters. Note that, the overall R was highest-i. e., 0. 995while overall MAE was 2. 059 for Trial 2 at default settings (chromosomes = 30, head size = 8, genes number = 3). However, in case of Trial 5, the R value was 0. 953 while its MAE was 1. 9865, because of which the latter trial was considered to be the final optimized model. Further changes in the GEP settings resulted in lower R values, higher R 2, and lower MAE.
section_name: Statistical Evaluation

section_content: Figure 3a -j presents the comparison of the regression slopes between modeled and experimental results for all the 10 GEP trials (training and validation datasets are represented by blue and red, respectively). It can be seen in Figure 3 that the CS28 is accurately predicted for almost all the proposed models wherein the lowest-maximum regression slopes for training and testing datasets are 0. 9468 (Trial 1)-0. 9728 (Trial 6) and 0. 947 (Trial 8)-0. 9973 (Trial 4), respectively. In the previous studies, slopes ranging between 0. 80 and 0. 90 were considered indicators of the achievement of higher performance, when modeled with GEP approach [29, 79, 80]. Based on their findings, it can be stated that all of the formulated GEP trials exhibited comparatively better performance. Trial 5 was selected as the optimal trial on the basis of overall R and MAE values (in the previous section). In addition, the nearness of the points to the ideal fit (1:1 line) and the presence of the majority of the points inside the allowable confidence interval elucidates the validity of GEP models [42]. The points which are closer to the regression line depict better prediction performance of a formulated GEP model [81]. While considering the plots of regression slopes (in Figure 3 ), it can be observed that Trial 6 and Trial 4 possessed the maximum slopes for training and validation datasets, respectively. On the contrary, the values of slope in the case of Trial 5 were recorded as 0. 957 (training) and 0. 961 (validation). Upon comparing the overall MAE values of the three trials, it is found that Trial 5 has the lowest value (1. 9865) in contrast to Trial 4 (2. 096) and Trial 6 (2. 0325). Moreover, the comparatively higher R value (0. 953) in case of Trial 5 further makes it a final choice to be considered as the most optimal model. Note that, an R value exceeding 0. 8 illustrates better performance of the GEP model [79, 82]. Therefore, the prediction accuracy of Trial 5 is higher.
section_name: Comparison of the Regression Slopes between Experimental and Predicted Results

section_content: A variety of statistical methods are used to investigate the performance of developed models. The values obtained from the predicted models divided by the experimental results-expressed as predicted-to-experimental ratio (ρ)-is also employed to evaluate the performance of the developed model [83]. The ratios obtained from dividing predicted results by experimental values are plotted in the form of Table 7, taking the bin range as 0. 2. It can be observed that all models interpreted 100% ρ values in between 0. 8 and 1. 2. This suggests that the accuracy of all the developed models is within ±20% error. The results of ρ for Trial 5-which is concluded as the best model achieved from evaluation in terms of statistical indices and comparison of regression slopes-was plotted in the form of a histogram, decreasing the bin range from 0. 2 to 0. 1 to see the accuracy of the respective model in a closer view shown in Figure 4. ρ values of 98. 17% and 91. 49% o were obtained within 0. 9 to 1. 1 for the training and validation datasets, respectively. This manifests the model error within ±10%, reflecting a more robust performance for the selected optimum model. A variety of statistical methods are used to investigate the performance of developed models. The values obtained from the predicted models divided by the experimental results-expressed as predicted-to-experimental ratio (ρ)-is also employed to evaluate the performance of the developed model [83]. The ratios obtained from dividing predicted results by experimental values are plotted in the form of Table 7, taking the bin range as 0. 2. It can be observed that all models interpreted 100% ρ values in between 0. 8 and 1. 2. This suggests that the accuracy of all the developed models is within ±20% error. The results of ρ for Trial 5-which is concluded as the best model achieved from evaluation in terms of statistical indices and comparison of regression slopes-was plotted in the form of a histogram, decreasing the bin range from 0. 2 to 0. 1 to see the accuracy of the respective model in a closer view shown in Figure 4. ρ values of 98. 17% and 91. 49% o were obtained within 0. 9 to 1. 1 for the training and validation datasets, respectively. This manifests the model error within ±10%, reflecting a more robust performance for the selected optimum model.
section_name: Predicted to Experimental Ratio (ρ)

section_content: The expression trees shown in Figure 5 and the MATLAB code obtained from the GEP model were used to extract mathematical equation (Equation ( 6 )) to predict the compressive strength of PET containing SCM based grout in accordance with the previous literature [84]. The prediction equations were obtained for the best model, Trial 5, in this study. The complexity of the mathematical equation depends upon the number of sub-ETs which depends on the number of genes. It is evident that Trial 5 contains three genes in its setting parameters; hence, four sub-ETs can be observed in Figure 5. Furthermore, Figure 5 shows linking functions (+, -, /, ×, cubic root), constants (C1 to C8), and variables in the form of d0 to d4 denoting PET percentage, percentage of SCM, flow, 1-day compressive strength, and 7-day compressive strength of cementitious grout, respectively, used in the formulation of GEP mathematical model. 

where, X 0 = PET, X 1 = SCM, X 2 = Flow, X 3 = CS 1d, and X 4 = CS 7d. Secondly, the parametric study in Figure 6b -f shows that with an increase in the dosage levels of the PET (0-10%), proportioning of SCM (0-10%), CS1d, and CS7d led to a positive growth in the CS28d of cementitious grout mixes. Generally, the increase in PET leads to a decrease in the compression strength; however, in the current study, it was found that with increasing PET the CS28d increased very slightly (up to 3. 61%). This can be attributed to the inclusion of the SCMs (i. e., fly ash, silica fume, and metakaolin). It is pertinent to mention that when the percentage of recycled PET fiber is higher, the slump test and compressive strength values decreased [89]. Generally, the inclusion of plastic fragments to
section_name: GEP Formulation

section_content: It is critical to conduct a number of assessments over AI formulated models to guarantee that the suggested models possess robustness and that they perform well for different unseen data. Sensitivity and parametric tests demonstrate their robustness [72, 85]. The sensitivity analysis (SA) on the entire database determines how susceptible a constructed model is to changes in the variables under consideration [79, 86]. The relative contributions of the input factors (PET, SCM, Flow, CS 1d, and CS 7d ) are taken into consideration here to forecast the CS 28d of waste PET/SCM blended cementitious grout by performing the SA. 

For a particular input variable, Si, the SA can be computed using the Equations ( 7 ) and (8), respectively (such that t = 156, in the present study). 

where fmax (S i ) and fmin (S i ) denotes, respectively, the maximum and minimum of forecasted CS 28d on the basis of the ith input domain, whereas the rest of the input variables remain constant at their mean. 

The SA value ranges from 0 to 1, indicating the relative contribution of each input variable to the predicted output compression strength (SA = 1 indicates higher impact, SA = 0 indicates least effect) [87, 88]. 

Firstly, the results of sensitivity analysis in Figure 6a shows a rising trend of input parameter contributions to determine the CS 28d of waste PET/SCM blended cementitious as CS 7d (30. 47%) > CS 1d (28. 89%) > SCM (18. 88%) > Flow (18. 53%) > PET (3. 23%). Therefore, CS 7d is most substantial factor that influences the long term CS 28d of cementitious grout mixes, followed by CS 1d, flow, and proportion of the SCMs. content [91]. Inclusion of PET fibers alongside silica fume-cemented sand increases the unconfined compression strength, as well as the energy absorption capacity [92] because silica fume tends to decrease the detrimental consequences on strength characteristics while increasing the compression strength [93, 94]. On the other hand, the increase in flow values negatively impacted the CS28d of cementitious grout mixes, which agrees with the findings of the past researchers-i. e., when the slump flow of concrete decreases, the CS increases [95] [96] [97] [98].
section_name: Sensitivity and Parametric Analysis

section_content: This study aims to determine the impact of PET/SCM (polyethylene terephthalate (PET) and two supplementary cementitious materials (SCMs)-i. e., fly ash (FA) and silica fume (SF)-on the 28-day compressive strength (CS28d) of cementitious grouts using GEP. A new predictive model in the form of a simple mathematical expression has been formulated to compute the CS28d of cementitious grout mixes using easily determinable input parameters, such as, PET, SCM, flow, CS1d, and CS7d. Based on this investigation, the following conclusions can be drawn:

1. Ten GEP trials were conducted, and several performance indices were recorded Secondly, the parametric study in Figure 6b -f shows that with an increase in the dosage levels of the PET (0-10%), proportioning of SCM (0-10%), CS 1d, and CS 7d led to a positive growth in the CS 28d of cementitious grout mixes. Generally, the increase in PET leads to a decrease in the compression strength; however, in the current study, it was found that with increasing PET the CS 28d increased very slightly (up to 3. 61%). This can be attributed to the inclusion of the SCMs (i. e., fly ash, silica fume, and metakaolin). It is pertinent to mention that when the percentage of recycled PET fiber is higher, the slump test and compressive strength values decreased [89]. Generally, the inclusion of plastic fragments to concrete reduces the fresh and dry density, thus decreasing the weight of the generated concrete [90] ; however, since SCMs were used in conjunction, the CS 28d of cementitious grout mixes increased with higher dosages of PET. In addition, the cementitious materials (fly ash as well as silica fume) compensate the loss of CS 28d due to addition of the PET content [91]. Inclusion of PET fibers alongside silica fume-cemented sand increases the unconfined compression strength, as well as the energy absorption capacity [92] because silica fume tends to decrease the detrimental consequences on strength characteristics while increasing the compression strength [93, 94]. On the other hand, the increase in flow values negatively impacted the CS 28d of cementitious grout mixes, which agrees with the findings of the past researchers-i. e., when the slump flow of concrete decreases, the CS increases [95] [96] [97] [98].
section_name: Conclusions

section_content: This study aims to determine the impact of PET/SCM (polyethylene terephthalate (PET) and two supplementary cementitious materials (SCMs)-i. e., fly ash (FA) and silica fume (SF)-on the 28-day compressive strength (CS 28d ) of cementitious grouts using GEP. A new predictive model in the form of a simple mathematical expression has been formulated to compute the CS 28d of cementitious grout mixes using easily determinable input parameters, such as, PET, SCM, flow, CS 1d, and CS 7d. Based on this investigation, the following conclusions can be drawn:

1. Ten GEP trials were conducted, and several performance indices were recorded for each trial. While investigating the effect of varying genetic setting parameters (chromosomes, head size, genes) to evaluate the CS 28d of cementitious grout mixes, it was noticed that the best hyper parameter for a particular GEP model strongly depended on trial and access method. The optimal genetic parameter for one GEP model for computation of CS 28d may be completely different from the optimum parameter for another model. Depending on variety of indices, the best GEP model was obtained for number of chromosomes = 50, head size = 9, and number of genes = 3. 

2. The optimum statistical indices obtained in the case of the finally selected optimal model (Trial 5) were obtained as RMSE (2. 423 and 2. 531), MAE (1. 918 and 2. 055), and R (0. 977 and 0. 975) for training and validation datasets, respectively. In addition, the MAE values depict 3. 6% (training) and 3. 8% (validation) averaged error in the developed model. These values are significantly lower, which indicates accuracy and robust performance in terms of experimental and modeled values of the CS 28d of cementitious grout mixes in the formulated GEP model. 

3. The GEP model performance was further supplemented with the help of other statistical evaluating indices, such as: (i) slope of regression line between experimental and predicted results; and (ii) predicted to experimental ratios (ρ) for all the models. It it worth noting that the best model yielded regression slopes of 0. 966 (training) and 0. 985 (validation), which are more proximal to unity (i. e., ideal slope). The ρ values for all the 10 trials interpreted 100% of the values lying within ±20%. Furthermore, the optimal model resulted in 98. 17% and 91. 49% of the ρ values within ±10% error, which further confirms the final selection of the optimum model. 

4. The MATLAB code extracted from the final GEP model was used to form a mathematical equation comprising easily determinable input parameters to evaluate the CS 28d of cementitious grout mixes, thus avoiding the laborious, time-consuming, and costly testing of the samples and thereby improving the cost-effectiveness of civil engineering projects. 

5. The sensitivity analysis revealed that CS 7d is the most significant parameter which impacts the long term CS 28d of cementitious grout mixes followed by CS 1d, flow, and proportion of the SCM. Alongside incorporating the different SCMs, it is imperative to mention that the negative influence of the PET dosage on CS 28d was substantially neutralized (depicting only a slight increase). Moreover, the parametric study revealed that increasing the PET and proportioning of SCM (0 to 10% each), CS 1d, and CS 7d resulted in a positive increase in the CS 28d of cementitious grout mixes. Conversely, the increase in flow values negatively impacted the CS 28d. The current GEP model may be effectively deployed for future purposes to evaluate the 28-day compression strength of cementitious grout mixes.
section_name: Conclusions

section_content: Acknowledgments: This research was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (project no. AN000316 ). The authors wish to express their gratitude for the financial support that has made this study possible.

section_content: Data Availability Statement: All the data utilized in this research are available on a reasonable request from the corresponding author.

section_name: Conflicts of Interest:

section_content: Author Contributions: K. K., F. E. J. and M. I. K. contributed to the design of this research. M. I. K. and M. A. A. -F. contributed to performing the experiments. K. K., F. E. J. and M. I. K. contributed to writing the 'Introduction' and 'Materials and Methods' sections of the manuscript. K. K. and M. N. A. provided visualizations and contributed to funding acquisition. F. E. J., M. I. and M. N. A. contributed to the modeling and prepared all Figures and Tables.

section_content: The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; nor in the decision to publish the results.
section_name: Conflicts of Interest: