Next Article in Journal
Numerical Simulation Study of Gas-Solid Heat Transfer and Decomposition Processes of Limestone Calcined with Blast Furnace Gas in a Parallel Flow Regenerative Lime Kiln
Next Article in Special Issue
Self-Healing Performance Assessment of Bacterial-Based Concrete Using Machine Learning Approaches
Previous Article in Journal
Study on Simulation and Experiment of Cu, C-Doped Ag/Ni Contact Materials
Previous Article in Special Issue
Experiment on the Properties of Soda Residue-Activated Ground Granulated Blast Furnace Slag Mortars with Different Activators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi Expression Programming Model for Strength Prediction of Fly-Ash-Treated Alkali-Contaminated Soils

1
Department of Civil and Environmental Engineering, College of Engineering, King Faisal University (KFU), Al-Ahsa 31982, Saudi Arabia
2
Department of Civil Engineering, National Institute of Technology Warangal, Warangal 506004, India
3
Osaimi Geotechnic Company, Tabuk 31952, Saudi Arabia
4
State Key Laboratory of Ocean Engineering, Shanghai Key Laboratory for Digital Maintenance of Buildings and Infrastructure, School of Naval Architecture, Ocean & Civil Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
5
Department of Civil Engineering, University of Engineering & Technology (UET), Peshawar 39161, Pakistan
6
Department of Structural Engineering, Military College of Engineering (MCE), National University of Science and Technology (NUST), Islamabad 44000, Pakistan
7
Department of Civil Engineering, CECOS University of IT and Emerging Sciences, Peshawar 25000, Pakistan
8
Department of Mechanical Engineering, College of Engineering, King Faisal University (KFU), Al-Ahsa 31982, Saudi Arabia
*
Author to whom correspondence should be addressed.
Submission received: 21 April 2022 / Revised: 24 May 2022 / Accepted: 31 May 2022 / Published: 6 June 2022

Abstract

:
Rapid industrialization is leading to the pollution of underground natural soil by alkali concentration which may cause problems for the existing expansive soil in the form of producing expanding lattices. This research investigates the effect of stabilizing alkali-contaminated soil by using fly ash. The influence of alkali concentration (2 N and 4 N) and curing period (up to 28 days) on the unconfined compressive strength (UCS) of fly ash (FA)-treated (10%, 15%, and 20%) alkali-contaminated kaolin and black cotton (BC) soils was investigated. The effect of incorporating different dosages of FA (10%, 15%, and 20%) on the UCSkaolin and UCSBC soils was also studied. Sufficient laboratory test data comprising 384 data points were collected, and multi expression programming (MEP) was used to create tree-based models for yielding simple prediction equations to compute the UCSkaolin and UCSBC soils. The experimental results reflected that alkali contamination resulted in reduced UCS (36% and 46%, respectively) for the kaolin and BC soil, whereas the addition of FA resulted in a linear rise in the UCS. The optimal dosage was found to be 20%, and the increase in UCS may be attributed to the alkali-induced pozzolanic reaction and subsequent gain of the UCS due to the formation of calcium-based hydration compounds (with FA addition). Furthermore, the developed models showed reliable performance in the training and validation stages in terms of regression slopes, R, MAE, RMSE, and RSE indices. Models were also validated using parametric and sensitivity analysis which yielded comparable variation while the contribution of each input was consistent with the available literature.

1. Introduction

The natural soil–water system can be altered by the interaction of pollutants released from mining, agricultural, and industrial activities. Rapid industrialization has manifested in the form of an exponential rise in the usage of chemical agents which include hydroxides and carbonates and bicarbonates [1,2]. Among the various pollutants, the adverse effects of alkali contamination on the soil–water system have been well established by earlier researchers [3,4]. The lower concentrations of alkali can alter the structure of the soil [5], and at higher concentrations, alkali can augment the formation of new compounds in the form of zeolites and ettringite [4,6]. The alkali interactions with the expansive clays and subsequent swelling are a natural phenomenon primarily driven by the presence of smectite group clay minerals (montmorillonite, illite, and vermiculite) with expanding lattices [5,7]. However, Rao et al. [1] identified the heaving of foundations resting on non-swelling soils contributed by acid and alkali contamination. The incidence of geotechnical failures induced by alkali interaction reported by Rao et al. [1] confirmed the adverse effects of alkali contamination. Traditionally, the swelling of soil was addressed by the utilization of conventional binders such as lime and cement. However, the efficiency of traditional binders was contradicted by Hunter (1990), who reported 3.2 m heave of lime-treated soils. Further, with the growing impetus to promote sustainable binders with relatively lower carbon emissions, the usage of industrial by-products has gained momentum. Among the various industrial by-products, the efficiency of fly ash (FA) as a sustainable stabilizer/binder has been broached by recent researchers [8,9,10,11,12]. Fly ashes inherently do not exhibit unconfined compressive strength (UCS) due to lack/loss of cohesion in either dry or fully saturated conditions [13,14]. However, the UCS of FA is derived either from the difference in capillary action between the coarse and fine fractions [15,16,17] or due to internal friction [18]. Though the studies related to FA treatment on alkali-induced swelling are well documented, the attempts to quantify the variation in strength characteristics are sparse [4,19]. Hence an attempt has been made to evaluate the influence of alkali contamination on the UCS of two types of clay soils, i.e., expansive black cotton (BC) soil and non-expansive kaolin soil, for varying curing periods (1, 7, 14, and 28 days). Further, the effectiveness of FA (10%, 15%, and 20% dry weight) as a stabilizer to enhance the UCSkaolin and UCSBC soils was also evaluated for the aforementioned varying curing periods.
To overcome the limitation of the laborious and time-consuming nature of laboratory studies for the evaluation of engineering characteristics, there is a growing trend of developing numerical models for solving engineering problems [20,21,22] and prediction models for swift estimation of engineering characteristics of soils [23,24]. Initially, most of the prediction models were developed on the basis of regression analysis with relatively limited databases. To overcome this limitation, the advent of artificial intelligence (AI)-based models is being promoted primarily due to their ability to estimate/predict results even for larger databases. Sinha et al. [25] are among the early researchers to introduce artificial neural networks (ANNs), a machine learning language of MATLAB, for the estimation of compaction characteristics for a database of 55 soil samples. The advent of the ANNs has resulted in their consistent usage in dealing with complex physical and mathematical problems [26,27]. Along similar lines, there has been a recent advancement in genetic programming (GP)-based models in the form of gene algorithms (GAs) with the objective of identifying optimized solutions. Cramer (1985) was the first to introduce GP, which was subsequently improved by Koza [28] with varying sizes and shapes. The most advanced methods among the existing linear-based GP techniques are genetic expression programming (GEP) and multi expression programming (MEP); both are genotype computation programming methods that generate tree-like models/programs. The limitations of GAs and GP are addressed by both the methods, which manifests in their efficiency and swift execution (2 to 4 times faster) [29,30,31]. Their unique feature is their ability to learn and adapt by varying their shape, size, and composition, which is akin to living organisms [32,33]. However, unlike GEP, MEP adopts a demonstrative approach and utilizes linear chromosomes for the encoding of programs/solutions. However, the final solution (chromosome) is chosen based on the fitness value of an individual chromosome. In general, the governing parameters in MEP include subpopulation size and number, code length, function set, and crossover probability. Recently, this approach has gained greater applications in the geotechnical engineering field for addressing varied problems which include the prediction of compaction characteristics [34,35], compressive strengths [36,37], permeability and compressibility characteristics [38], deformation modulus [39], soil water characteristic curves, and peak ground acceleration [40]. However, attempts to generate models for the prediction of geotechnical characteristics of the contaminated soils are scarce. Moreover, considering the greater uncertainty associated with clayey soils, there is a need to develop a more realistic MEP model to estimate their strength characteristics. Thus, in the present study, attempts were made to develop an MEP model for the prediction of the UCS of alkali-contaminated clayey soils. For the development of the model, a comprehensive database of 384 soil samples was considered (part of the data was sourced from Ashfaq et al. [19]), and a brief description of the results and the contributing mechanism is presented in the following sections.

2. Materials and Methods

2.1. Laboratory Studies

The physical properties of the kaolin and BC soils are presented in Table 1, and it can be noted that inherently both the soils are classified as CH (highly plastic clays). The kaolin soil was commercially procured from Heilen Biopharm Pvt. Ltd. The FA considered in the study is sourced from Kakatiya Thermal Power Station (18°23′00.8″ N; 79°49′33.6″ E) Bhupalpally, Telangana, India, and its mineral composition is presented in Table 2 and is classified as class F (due to <30% CaO content). On the other hand, the BC soil was procured from the National Institute of Technology, Warangal (campus), by open excavation at a depth of 1 m (soil profile studies confirmed the presence of BC soil up to a depth of 1.8 m from the ground level), as shown in Table 1. The samples were oven-dried, and a 0.425 mm passing fraction was considered for the UCS tests. Firstly, the soil samples were inundated at varying concentrations of alkali for predetermined curing periods followed by UCS testing. Prior to the UCS testing, the samples were compacted to 95% of MDD of the respective soils, and the strain rate of 0.35 mm/min was applied in accordance with ASTM D2166. To maintain consistency and repeatability, the alkali dosage was maintained at OMC (presented in Table 3) for all the soil–alkali combinations. For the stabilized case, varying dosages of FA (10%, 15%, and 20% by dry weight) were added to the alkali–soil combinations.

2.2. MEP Model Development

As stated earlier, the MEP approach is among the most significant linear configurations of the genetic programming (GP) series since it has the ability to deliver simplistic mathematical formulae to forecast a particular prediction model [35,41]. Therefore, the formulization of the UCSkaolin and UCSBC soil was performed in the Multi-Expression Programming X (MEPX version 2021.08.28.0-beta) by incorporating experimental records, as shown in Table 3. Sufficient laboratory test data of 384 different soils for the UCS prediction of FA-treated alkali-contaminated soils were collected by performing an experimental study [42]. The MEP genes are the substrings of varying lengths that keep the chromosomal length constant and equivalent to the total genes on each chromosome. Each gene provides instructions for making a function or a terminal sign, whereas a gene encoding a function includes addresses to the function parameters. The function arguments always have lower parameter estimates than the location of the function on that chromosome [41]. A detailed methodology for generating equations is provided here, and details of the advanced GP approach (MEP simulation) can be found elsewhere [34,35,39,42,43,44,45].
Two-thirds of the entire data was considered for the MEP model development, whereas one-third was utilized to validate the formulated model. Table 4 shows the maximum and minimum values of the input (FA dosage, alkali concentration, and curing days) and output parameters (UCSkaolin and UCSBC) used to perform strength prediction of FA-treated alkali-contaminated soils in the case of both training and testing data. The maximum as well as the minimum values of all the input and output characteristics have been tabulated in Table 4. Figure 1 shows the frequency histograms (i.e., the scatter of the data) of the input attributes considered in the current study. The curves are smooth and uniformly distributed, which shows a good type of data. In addition, standard deviation (SD), kurtosis, and skewness for all the parameters are given. A smaller SD shows that the parameters are near the respective average value. The kurtosis value represents the sharpness of the peak of a frequency distribution curve. It clarifies the shape of probability distribution [34]. It is pertinent to mention that the kurtosis value is only useful when used in conjunction with the SD value. It is possible that an attribute might have a high kurtosis (bad), but the overall standard deviation is low (good). A kurtosis value of ±1 is considered very good for most psychometric uses, but ±2 is also usually acceptable. The kurtosis values of FA dosage, alkali concentration, curing days, and UCSkaolin approach zero and therefore represent a mesokurtic distribution which can be seen in the histogram plot, i.e., Figure 1. However, the kurtosis value in the case of UCSkaolin is comparatively higher, which represents a leptokurtic distribution (Figure 1). Lastly, the skewness depicts the extent to which a distribution of values deviates from symmetry around the mean. Bryne (2010) argued that data are considered to be normal if skewness is between −2 and +2.
Furthermore, Table 5 shows the Pearson correlation coefficient values (represented by ‘r’) for the input parameters and the two output parameters, i.e., UCSkaolin and UCSBC. It can be seen that the impact of all the three input parameters is linearly increasing (because r-values are positive). In the case of both UCSkaolin and UCSBC, the order of increasing impact of parameters follows the order: FA dosage > alkali concentration > curing age.
The details of 36 different trials (18 for each type of soil) undertaken to develop a model for the UCSkaolin and UCSBC with an optimal combination of hyperparameters are provided in Table 6. The set of hyperparameters (i.e., number of subpopulations, size of subpopulations, code length, tournament size, and number of generations) was varied in this study to achieve the optimal performance of models. A single parameter was modified whereas the rest were kept unchanged (as shown in Table 7) in an attempt to investigate the effect of different code settings on the correlation coefficient (R) and the mean squared error (MSE), as shown in Figure 2 and Figure 3, respectively.
In order to evaluate the UCSkaolin, the best performance was noted in the case of Trial 18 (R = 0.9465, Averaged MSE = 1245), wherein the number of subpopulations, size of subpopulation, code length, number of generations, and tournament size were kept as 20, 1000, 100, 150, and 6, respectively. On the other hand, in determining the UCSBC, the best performance was noted in the case of Trial 14 (R = 0.9672, Averaged MSE = 2220), wherein the number of subpopulations, size of subpopulation, code length, number of generations, and tournament size equal 100, 2000, 80, 150, and 2, respectively.
Using the above-mentioned adjusted hyperparameter setting, the simplified mathematical expressions given for the UCSkaolin (Equation (1)) and UCSBC (Equation (2)) were obtained, via C++ code, in order to predict the targeted UCS.
UCS Kaolin = 7     X 0 10     X 1 + X 2 ( 5     ( X 0 X 1 ) ) X 1 + 2     X 2 + ( X 1     X 2 + ( 2     X 0     X 1 ) X 2 ) X 0 X 2 +   X 0     X 1 + ( X 0     X 1     X 2 ) 8 + 258
UCS BC = 9     X 0     X 1   +   2     X 1 ( X 1 + X 2 + 1 )     ( 4     ( X 2   +   1 ) ) X 0     1 + X 1 ( X 1 + ( 4     ( X 2   +   1 ) ) X 0     1 + X 1 ( X 0     1 ) ) ( X 0     X 1 ) X 2   +   162
where X 0 , X 2 , and X 3 represent fly ash dosage (%), alkali concentration (N), and curing period (days), respectively.

3. Results and Discussion

3.1. Strength Characteristics

3.1.1. Effect of Alkali Contamination

The variation in the UCSkaolin and UCSBC with curing periods and alkali concentrations is presented in Figure 4. The UCSkaolin and UCSBC linearly decreased with the rise in alkali concentration, and BC soil exhibited a relatively greater fall in the UCS compared to kaolin soil. With the increase in curing periods, both the soils increased linearly for the controlled case. However, under the contaminated case, the kaolin soil exhibited a slight increase in the UCS, whereas the UCSBC remained constant at lower curing periods and was drastically reduced at higher curing periods. The variation in UCSkaolin and UCSBC in the untreated case can be attributed to their inherent mineralogical difference which enables the formation of primary hydration compounds [4]. Under a contaminated scenario, the linear decrease in UCSkaolin and UCSBC soils may be attributed to the increase in charge of clay particles with the pH of the soil. The rise in pH contributes to the subsequent dissolution of silica which varies with the size and crystallinity of quartz, a commonly found mineral in both the soils [4]. Furthermore, an increase in the UCSkaolin with the curing period may be attributed to the precipitation of hydration compounds such as nontronite and sodium silicate hydrate. Similar observations for clayey soils were made by Sivapullaiah and Reddy [4].

3.1.2. Effect of FA Dosage and Curing Period

The variation in the UCSkaolin and UCSBC with the alkali concentrations and FA dosage is presented in Figure 5. Considering brevity, the results pertaining to a 28-day curing period have been presented here. It is evident from the results that the FA addition has contributed sufficiently to the linear increase in the UCSkaolin and UCSBC for both the controlled and alkali-contaminated cases. In contrast to the contaminated case, the increase in the UCSBC is substantially higher with an increment of more than 900% compared to the 350% increase noted for kaolin soil. The increase in UCSkaolin and UCSBC is more pronounced at higher concentrations. The linear increase in UCS of both soils is attributed to the decrease in clay content with the FA addition [2]. The greater increment at higher concentration is attributed to the greater affinity of dissolved silica (due to higher concentration of alkali) to react with calcium from FA and subsequent formation of pozzolanic compounds. The pozzolanic compounds formed not only resist alkali attack on mineral phases of soil but also offer greater resistance to compressive loading which is manifested in the form of increased UCS [4,19].

3.2. Comparison between Experimental and Predicted Results

This segment focuses on efficacy examination and relative study of the MEP models generated to compute the UCSkaolin and UCSBC, using a variety of performance indices. To evaluate the prediction efficiency and accuracy of the proposed MEP models using MEPX software, eight analytical standard indicators, namely regression line slope, correlation coefficient (R), root mean squared error (RMSE), mean absolute error (MAE), root squared error (RSE), relative root mean square error (RRMSE), Nash–Sutcliffe efficiency (NSE) and performance index (ρ) were used in this study [46,47]. These performance measures are defined by the following Equation (3) to Equation (9):
                                            R = i = 1 n ( E i E ¯ i ) ( P i P ¯ i ) i = 1 n ( E i E i ) 2 i = 1 n ( P i P ¯ i ) 2
R M S E = i = 1 n ( E i P i ) 2 n
M A E = i = 1 n | E i P i | n
R S E = i = 1 n ( P i E i ) 2 i = 1 n ( E ¯ E i ) 2
R R M S E = 1 | E ¯ | i = 1 n ( E i P i ) 2 n
N S E = 1 i = 1 n ( E i P i ) 2 i = 1 n ( P i P ¯ i ) 2  
ρ = ( 1 | E ¯ | i = 1 n ( E i P i ) 2 n ) 1 + R
where Ei and Pi are the ith actual and predicted output values, respectively; E ¯ i and P ¯ i are the average values of the actual and predicted output values, respectively; and n is the number of samples. In addition, the objective function (OBF), as given in Equation (10), shall have a minimum value for better formulation of the model. A smaller OBF helps in overcoming the overfitting problem. A value approaching zero exhibits an excellent predictive capability.
O B F = ( n t r a i n n v a l i d n ) ρ t r a i n + 2 ( n v a l i d n ) ρ v a l i d
To construct accurate and robust AI-based predictive models, the ratio of experimental values and inputs in the experimental database (as shown in Table 3) must be greater than 3 and ideally greater than 5 [48]. In the current investigation, the prescribed ratio is 269/3 = 89.66 (training set) and 115/3 = 38.33 (validation set), which is significantly within safe limits and therefore depicts the robustness and superiority of the developed MEP models for the kaolin and BC treated soils. The observed (actual) and forecasted UCSkaolin and UCSBC in the training and validation phases, as well as the efficacy metrics (i.e., slope, R, RMSE, MAE, RSE, RRMSE, and ρ), are shown in Figure 6a,b, respectively. The 45° regression line with a horizontal axis depicts the ideally fit (1:1) line having an inclination corresponding to 1 [49,50]. For good, reliable, and highly correlated models, the dispersion pattern of the data points should be closer to the diagonal line crossing the origin, with a trend line of slope approximately equaling unity, R-value greater than 0.8, and reduced error measurements (i.e., R, RMSE, MAE, RSE, RRMSE, and ρ), as shown in Figure 6 and Table 8. For both the kaolin and BC soil, the slopes of the trend lines are closer to 1 (0.90: training, 1.01: validation; 0.97: training, 0.96: validation, respectively). In addition, the R is above 0.8 (closer to 1) for both types of soils, which reflects a reasonably strong correlation between the model predicted outputs (i.e., UCSkaolin and UCSBC) and experimental observations. Furthermore, the OBF value of kaolin soil was 0.025694481, whereas that of BC soil was 0.025050897 in the current study.
Only a higher R-value is not the sole indication of the reliability and accuracy of the machine learning models [34]. Therefore, a number of error measurements were considered to validate the robustness of the developed models. These error metrics include R, RMSE, MAE, RSE, RRMSE, and ρ. The optimizer of the MEP algorithm was set to minimize the MSE while increasing the R statistic. In each model (kaolin or BC soil), the MSE and MAE were relatively low as compared to the maximum expected output, while the RSE reached zero. The optimized UCSkaolin model has MSE and MAE equaling 1245 and 19.6 and 2220 MPa and 30 for the training and validation phases, respectively. Likewise, the discussed attributes are also lower for the optimized UCSBC model. Furthermore, for both optimized models, the RSE tends to approach zero in each phase (i.e., training and validation), confirming their superior functionality. The consistent and accurate performance of the developed models is due to the structural flow of the MEP algorithm. The MEP follows the reproduction procedure to move the relevant information to the subsequent generation and uses the mutation function for optimization inside the chosen chromosomes.
Thus, the predefined configuration of the function is not taken into consideration [51,52]. In addition, the MEP technique produces randomized functions and selects the one that best fits the experimental results [33,40,53].
It is essential to further validate the accuracy of the developed MEP models using the values of the residual error, i.e., the difference between the model-estimated and experimental UCS [54,55]. The positive/negative minimum and maximum error obtained for the UCSkaolin model (Figure 7a) are −160 kPa and 100 kPa, respectively, and are ±130 kPa for the UCSBC model (Figure 7b). The majority of the error readings run along the x-axis, indicating a significant frequency of low error values. In conjunction with significantly higher correlations and reduced error measurements, the proposed models could be advantageously employed for the prediction of UCSkaolin and UCSBC, assisting practitioners and designers to save time and skip costly laboratory tests.
The plot of actual experimental values and the ultimate response of the MEP model for estimation of UCSkaolin and UCSBC can be seen in Figure 8a,b, respectively. In each case, the modeled values of the training and validation phases almost go along the observed (experimental) output, which shows the efficiency and accuracy of the formulated MEP models.
For the kaolin soil, the MAE and RMSE of the validation dataset are 6.31% and 7.94% lesser than the training dataset, respectively, and 9.28% and 26.66% lesser, respectively, in the case of BC soil. The improved performance in the testing stage depicts that the proposed MEP models have effectively learned the non-linear relationships among the inputs and response parameters with considerably lower error statistics and higher generalization capability [56,57]. Thus, the proposed model can be used for the prediction of UCSkaolin and UCSBC soil, which will aid in avoiding the heavy testing process.

3.3. Model Validity

The validity of the model is an important aspect of the AI modeling process. The model may perform better during the training stage for one set of data, whereas it may yield decreased performance for a new dataset. Therefore, the AI model shall be validated using an unused dataset to investigate the accuracy of the developed model for future applications [46,58,59]. As described in Section 3.2, the developed MEP model was validated using 30% of the experimental data; however, for further validation, the simulated dataset was created to evaluate the effect of contributing parameters on UCSkaolin and UCSBC shown as sensitivity and parametric analysis.

Sensitivity Analysis and Parametric Study of MEP Model

The testing of ML-based simulations is critical to ensuring that the recommended models are trustworthy and continue to perform well over a variety of datasets. The goal of sensitivity and parametric research is to confirm the efficacy of the proposed MEP models in terms of their interdependency on physical events [60,61,62]. The sensitivity analysis (SA) of the models on the complete dataset demonstrates how sensitive a generated model is to a change in the input variable in question [57,61,63]. The SA is being used to evaluate the impact of the input factors employed in this study on the anticipated UCS of contaminated soils.
For a specific independent variable ( Y i ) , the SA is carried out with Equations (11) and (12) for the overall experimental database considered in the current research. This means that one of the independent variables was changed between its extreme values while keeping the remaining input variables at their average values, and the output was recorded in the form of f(Yi). Next, the second independent variable was changed and the output was monitored.
R k = f m a x ( Y k ) f m i n ( Y k )
R e l a t i v e   I m p o r t a n c e   ( % ) = S A   ( % ) = R k n j = 1 R j 100
f m i n ( Y k ) and f m a x ( Y k ) are the minimum and maximum values of the anticipated results based on kth domain of the input variable in the preceding equations, with the remaining inputs kept at their mean. The results of SA can be observed in Figure 9, which shows that the curing period and alkali concentration have almost equal contributions to yielding UCS. The FA dosage contributes 13.37% among the three attributes. In the case of BC soil, the curing period significantly outperforms the other two variables; however, alkali concentration and fly ash dosage also contribute 22.96 and 7.73%, which is an important aspect to consider when investigating FA-incorporated alkali-contaminated soil.
To begin, Figure 10 visually represents the parametric analysis of the inputs used in this work (FA dosage, alkali content, and curing duration) for the prediction of UCS of kaolin and BC soil. The UCS of kaolin soil varies linearly with the amount of fly ash and alkali content, and a second-order polynomial trend is detected for the curing duration. In the case of BC soil, a straightforward linearly rising trend is found for each input (FA dose, alkali content, and curing duration). The increase in UCS of both types of soils with the curing duration is obviously according to the physical process involved, correctly captured by the MEP models, thus validating the models in this respect. Zha et al. [64] found a significant increase in UCS up to 28 days of curing while investigating the stabilization of metal-contaminated soil by alkaline residue. A similar increasing trend in UCS with an increase in curing duration was observed by Fasihnikoutalab et al. [65]. An increase in UCS was also observed with a rise in alkali content and fly ash while investigating alkali-activated geopolymer-incorporated kaolin soil [65]. Therefore, the developed models are deemed validated on the basis of the new experimental dataset and simulated data, which shows the model behavior similar to that of the physical process involved in alkali-activated kaolin and black cotton soil.

4. Conclusions

In the present experimental-cum-modeling study, the effect of alkali contamination on the strength characteristics of two clayey soils (i.e., kaolin and BC soil) has been evaluated. The efficiency of FA in remediating the alkali-induced effects was also assessed. Finally, the results were utilized by formulating MEP-based computational prediction models for computing the UCS of both types of the soils (UCSkaolin and UCSBC) to overcome the demerits of laborious laboratory testing, cost, and time. The following conclusions can be drawn from the study:
  • The inundation of kaolin and BC soils in alkali solution caused the UCS property to decrease. The higher concentrations posed a significant impact in lowering the UCSkaolin and UCSBC. On the contrary, the FA treatment of alkali-contaminated soils resulted in a linear increase in the UCSkaolin and UCSBC, and an increase of 7-fold was witnessed for the BC soil. Hence, it is concluded that the alkali contamination acted as an activator for a subsequent pozzolanic reaction when FA was incorporated.
  • In order to obtain the optimal MEP model for predicting the UCSkaolin and UCSBC, a total of 18 trials (each) were undertaken while considering the variation in (a) number of subpopulations, (b) subpopulation size, (c) code length, (d) tournament size, and (e) number of generations. The corresponding performance of all the trials was evaluated using a variety of performance indices, i.e., correlation coefficient and averaged MSE value. The best MEP model (kaolin and BC soil) was achieved in the case of 20 and 70 subpopulations, 1000 and 50 subpopulation size, 100 each code length, 6 each tournament size, 150 and 100 number of generations, 0.9465 and 0.9538 R-value, and 1245 kPa and 4400 kPa averaged MSE value, respectively.
  • Simple regression equations developed in this study (Equations (1) and (2)) for kaolin and BC contaminated soils can readily be used to forecast the UCS property. The equations have been generated from relatively high accuracy models evaluated using R, MAE, RMSE, and RSE (0.937, 19.6, 18.271, 0.128 and 0.956, 30, 17.151, 0.108) for the training data of kaolin and BC soils, respectively.
  • The generated models were evaluated using parametric and sensitivity analysis as second-level validation. The results obtained from the parametric study manifested a variation in UCS conforming to the literature for kaolin and BC soil with the change in the given input parameters. The sensitivity analysis of kaolin soil showed that curing period and alkali concentration had comparable contributions, followed by the FA dosage, whereas for BC, soil the following increasing trend was observed: curing period > alkali concentration > FA dosage.

Author Contributions

Conceptualization, K.K., M.I. and F.I.S.; Data curation, M.I., M.A.K. and F.E.J.; Formal analysis, M.A.K., M.I.F. and F.E.J.; Funding acquisition, K.K., M.N.A., F.I.S. and F.E.J.; Investigation, M.A. and M.I.; Methodology, M.A. and M.A.K.; Project administration, K.K. and M.N.A.; Resources, K.K., M.N.A. and F.I.S.; Software, M.A. and F.E.J.; Supervision, K.K.; Validation, F.I.S. and M.I.F.; Visualization, M.N.A. and M.I.F.; Writing—original draft, K.K., M.A., M.I. and F.E.J.; Writing—review & editing, K.K. and M.I. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (GRANT556).

Acknowledgments

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

Conflicts of Interest

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; or in the decision to publish the results.

References

  1. Rao, S.; Rao, K.S. Ground Heave from Caustic Soda Solution Spillage—A Case Study. Soils Found. 1994, 34, 13–18. [Google Scholar] [CrossRef] [Green Version]
  2. Sivapullaiah, P.V.; Manju, M. Effect of zeolitization on the volume-change behaviour of kaolinitic soils. Can. Geotech. J. 2006, 43, 969–978. [Google Scholar] [CrossRef]
  3. Sivapullaiah, P.V. Manju Ferric chloride treatment to control alakli induced heave in weathered red earth. Geotech. Geol. Eng. 2006, 24, 1115–1130. [Google Scholar] [CrossRef]
  4. Sivapullaiah, P.; Reddy, P.H.P. Fly ash to control alkali-induced volume changes in soils. Proc. Inst. Civ. Eng.-Ground Improv. 2009, 162, 167–173. [Google Scholar] [CrossRef]
  5. Mitchell, J.K.; Soga, K. Fundamentals of Soil Behavior; John Wiley & Sons: New York, NY, USA, 2005; Volume 3. [Google Scholar]
  6. Sivapullaiah, P.V.; Sankara, G.; Allam, M.M. Mineralogical changes and geotechnical properties of an expansive soil interacted with caustic solution. Environ. Earth Sci. 2010, 60, 1189–1199. [Google Scholar] [CrossRef]
  7. Jalal, F.E.; Xu, Y.; Jamhiri, B.; Memon, S.A. On the Recent Trends in Expansive Soil Stabilization Using Calcium-Based Stabilizer Materials (CSMs): A Comprehensive Review. Adv. Mater. Sci. Eng. 2020, 2020, 1510969. [Google Scholar] [CrossRef] [Green Version]
  8. Moghal, A.A.B. State-of-the-Art Review on the Role of Fly Ashes in Geotechnical and Geoenvironmental Applications. J. Mater. Civ. Eng. 2017, 29, 04017072. [Google Scholar] [CrossRef]
  9. Ashfaq, M.; Moghal, A.A.B.; Basha, B.M. Reliability-Based Design Optimization of Chemically Stabilized Coal Gangue. J. Test. Eval. 2021, 51. [Google Scholar] [CrossRef]
  10. Indiramma, P.; Sudharani, C.; Needhidasan, S. Utilization of fly ash and lime to stabilize the expansive soil and to sustain pollution free environment—An experimental study. Mater. Today Proc. 2020, 22, 694–700. [Google Scholar] [CrossRef]
  11. Rakhimova, N.; Rakhimov, R.Z. Toward clean cement technologies: A review on alkali-activated fly-ash cements incorporated with supplementary materials. J. Non-Cryst. Solids 2019, 509, 31–41. [Google Scholar] [CrossRef]
  12. Firat, S.; Khatib, J.; Yilmaz, G.; Comert, A.T. Effect of curing time on selected properties of soil stabilized with fly ash, marble dust and waste sand for road sub-base materials. Waste Manag. Res. J. A Sustain. Circ. Econ. 2017, 35, 747–756. [Google Scholar] [CrossRef]
  13. DiGioia, A.M.; Nuzzo, W.L. Fly ash as structural fill. J. Power Div. 1972, 98, 77–92. [Google Scholar] [CrossRef]
  14. Gray, D.H.; Lin, Y.-K. Engineering properties of compacted fly ash. J. Soil Mech. Found. Div. 1972, 98, 361–380. [Google Scholar] [CrossRef]
  15. Leonards, G.A.; Bailey, B. Pulverized coal ash as structural fill. J. Geotech. Eng. Div. 1982, 108, 517–531. [Google Scholar] [CrossRef]
  16. Singh, M.; Garg, M. Cementitious binder from fly ash and other industrial wastes. Cem. Concr. Res. 1999, 29, 309–314. [Google Scholar] [CrossRef]
  17. Antiohos, S.; Tsimas, S. Activation of fly ash cementitious systems in the presence of quicklime: Part I. Compressive strength and pozzolanic reaction rate. Cem. Concr. Res. 2004, 34, 769–779. [Google Scholar] [CrossRef]
  18. Singh, S.; Panda, A. Utilization of fly ash in geotechnical construction. Proc. Indian Geotech. Conf. 1996, 1, 547–550. [Google Scholar]
  19. Ashfaq, M.; Heeralal, M.; Reddy, P.H.P. A Study on Strength Behavior of Alkali-Contaminated Soils Treated with Fly Ash. In Recycled Waste Materials; Springer: Berlin/Heidelberg, Germany, 2019; pp. 137–143. [Google Scholar]
  20. Bhatti, U.N.; Bashmal, S.; Khan, S.; Ali, S. Design and Optimization of 6-DOF Platform Top Plate under Realistic Joint Conditions. Int. J. Simul. Model. 2020, 19, 53–64. [Google Scholar] [CrossRef]
  21. Bhatti, U.N.; Khan, S.; Ali, S.; Horoub, M.M. A Numerical Modeling Study of the Effects of Various Joint Boundary Conditions on Stiffness Behavior of 6DOF Platform’s Top Plate. In Proceedings of the 2018 9th International Conference on Mechanical and Aerospace Engineering (ICMAE), Budapest, Hungary, 10–13 July 2018; pp. 243–248. [Google Scholar]
  22. Ali, S.; Kamran, M.A.; Khan, S. Effect of baffle size and orientation on lateral sloshing of partially filled containers: A numerical study. Eur. J. Comput. Mech. 2017, 26, 584–608. [Google Scholar] [CrossRef]
  23. Onyelowe, K.C.; Jalal, F.E.; Onyia, M.E.; Onuoha, I.C.; Alaneme, G.U. Application of Gene Expression Programming to Evaluate Strength Characteristics of Hydrated-Lime-Activated Rice Husk Ash-Treated Expansive Soil. Appl. Comput. Intell. Soft Comput. 2021, 2021, 6686347. [Google Scholar] [CrossRef]
  24. Al-Busultan, S.; Aswed, G.K.; Almuhanna, R.R.A.; Rasheed, S.E. Application of Artificial Neural Networks in Predicting Subbase CBR Values Using Soil Indices Data. In Proceedings of the IOP Conference Series: Materials Science and Engineering, Chennai, India, 16–17 September 2020; IOP Publishing: Bristol, UK, 2020; Volume 671, p. 012106. [Google Scholar]
  25. Sinha, S.K.; Wang, M.C. Artificial Neural Network Prediction Models for Soil Compaction and Permeability. Geotech. Geol. Eng. 2007, 26, 47–64. [Google Scholar] [CrossRef]
  26. Naqvi, S.R.; Tariq, R.; Hameed, Z.; Ali, I.; Taqvi, S.A.; Naqvi, M.; Niazi, M.B.; Noor, T.; Farooq, W. Pyrolysis of high-ash sewage sludge: Thermo-kinetic study using TGA and artificial neural networks. Fuel 2018, 233, 529–538. [Google Scholar] [CrossRef]
  27. Bong, J.T.; Loy, A.C.M.; Chin, B.L.F.; Lam, M.K.; Tang, D.K.H.; Lim, H.Y.; Chai, Y.H.; Yusup, S. Artificial neural network approach for co-pyrolysis of Chlorella vulgaris and peanut shell binary mixtures using microalgae ash catalyst. Energy 2020, 207, 118289. [Google Scholar] [CrossRef]
  28. Koza, J. On the programming of computers by means of natural selection. Genet. Program. 1992. [Google Scholar]
  29. Ferreira, C. Gene Expression Programming: Mathematical Modeling by an Artificial Intelligence; Springer: Berlin/Heidelberg, Germany, 2006; Volume 21. [Google Scholar]
  30. Iqbal, M.; Zhang, D.; Jalal, F.E.; Javed, M.F. Computational AI prediction models for residual tensile strength of GFRP bars aged in the alkaline concrete environment. Ocean Eng. 2021, 232, 109134. [Google Scholar] [CrossRef]
  31. Khan, S.; Horoub, M.M.; Shafiq, S.; Ali, S.; Bhatti, U.N. Optimization of vehicle suspension system using genetic algorithm. In Proceedings of the 2019 IEEE 10th International Conference on Mechanical and Aerospace Engineering (ICMAE), Brussels, Belgium, 22–25 July 2019; pp. 203–207. [Google Scholar]
  32. Iqbal, M.; Zhao, Q.; Zhang, D.; Jalal, F.E.; Jamal, A. Evaluation of tensile strength degradation of GFRP rebars in harsh alkaline conditions using non-linear genetic-based models. Mater. Struct. 2021, 54, 1–20. [Google Scholar] [CrossRef]
  33. Jalal, F.E.; Xu, Y.; Iqbal, M.; Javed, M.F.; Jamhiri, B. Predictive modeling of swell-strength of expansive soils using artificial intelligence approaches: ANN, ANFIS and GEP. J. Environ. Manag. 2021, 289, 112420. [Google Scholar] [CrossRef] [PubMed]
  34. Jalal, F.E.; Xu, Y.; Iqbal, M.; Jamhiri, B.; Javed, M.F. Predicting the compaction characteristics of expansive soils using two genetic programming-based algorithms. Transp. Geotech. 2021, 30, 100608. [Google Scholar] [CrossRef]
  35. Wang, H.-L.; Yin, Z.-Y. High performance prediction of soil compaction parameters using multi expression programming. Eng. Geol. 2020, 276, 105758. [Google Scholar] [CrossRef]
  36. Çanakcı, H.; Baykasoğlu, A.; Güllü, H. Prediction of compressive and tensile strength of Gaziantep basalts via neural networks and gene expression programming. Neural Comput. Appl. 2009, 18, 1031. [Google Scholar] [CrossRef]
  37. Mousavi, S.M.; Aminian, P.; Gandomi, A.H.; Alavi, A.H.; Bolandi, H. A new predictive model for compressive strength of HPC using gene expression programming. Adv. Eng. Softw. 2012, 45, 105–114. [Google Scholar] [CrossRef]
  38. Johari, A.; Habibagahi, G.; Ghahramani, A. Prediction of SWCC using artificial intelligent systems: A comparative study. Sci. Iran. 2011, 18, 1002–1008. [Google Scholar] [CrossRef] [Green Version]
  39. Alavi, A.H.; Gandomi, A.H.; Sahab, M.G.; Gandomi, M. Multi expression programming: A new approach to formulation of soil classification. Eng. Comput. 2009, 26, 111–118. [Google Scholar] [CrossRef]
  40. Cabalar, A.F.; Cevik, A. Genetic programming-based attenuation relationship: An application of recent earthquakes in turkey. Comput. Geosci. 2009, 35, 1884–1896. [Google Scholar] [CrossRef]
  41. Oltean, M.; Grosan, C. A comparison of several linear genetic programming techniques. Complex Syst. 2003, 14, 285–314. [Google Scholar]
  42. Oltean, M.; Dumitrescu, D. Multi expression programming. J. Genet. Program. Evolvable Mach. Kluwer Second. Tour Rev. 2002. [Google Scholar] [CrossRef]
  43. Gandomi, A.; Faramarzifar, A.; Rezaee, P.G.; Asghari, A.; Talatahari, S. New design equations for elastic modulus of concrete using multi expression programming. J. Civ. Eng. Manag. 2015, 21, 761–774. [Google Scholar] [CrossRef] [Green Version]
  44. Fallahpour, A.; Olugu, E.U.; Musa, S.N. A hybrid model for supplier selection: Integration of AHP and multi expression programming (MEP). Neural Comput. Appl. 2015, 28, 499–504. [Google Scholar] [CrossRef]
  45. Fallahpour, A.; Wong, K.Y.; Rajoo, S.; Tian, G. An evolutionary-based predictive soft computing model for the prediction of electricity consumption using multi expression programming. J. Clean. Prod. 2020, 283, 125287. [Google Scholar] [CrossRef]
  46. Biswas, R.; Bardhan, A.; Samui, P.; Rai, B.; Nayak, S.; Armaghani, D.J. Efficient soft computing techniques for the prediction of compressive strength of geopolymer concrete. Comput. Concr. 2021, 28, 221–232. [Google Scholar]
  47. Biswas, R.; Rai, B. Effect of cementing efficiency factor on the mechanical properties of concrete incorporating silica fume. J. Struct. Integr. Maint. 2020, 5, 190–203. [Google Scholar] [CrossRef]
  48. Iqbal, M.F.; Javed, M.F.; Rauf, M.; Azim, I.; Ashraf, M.; Yang, J.; Liu, Q.-F. Sustainable utilization of foundry waste: Forecasting mechanical properties of foundry sand based concrete using multi-expression programming. Sci. Total Environ. 2021, 780, 146524. [Google Scholar] [CrossRef] [PubMed]
  49. Khan, M.A.; Zafar, A.; Akbar, A.; Javed, M.; Mosavi, A. Application of Gene Expression Programming (GEP) for the Prediction of Compressive Strength of Geopolymer Concrete. Materials 2021, 14, 1106. [Google Scholar] [CrossRef]
  50. Azim, I.; Yang, J.; Iqbal, M.F.; Javed, M.F.; Nazar, S.; Wang, F.; Liu, Q.-F. Semi-analytical model for compressive arch action capacity of RC frame structures. Structures 2020, 27, 1231–1245. [Google Scholar] [CrossRef]
  51. Faradonbeh, R.S.; Hasanipanah, M.; Amnieh, H.B.; Armaghani, D.J.; Monjezi, M. Development of GP and GEP models to estimate an environmental issue induced by blasting operation. Environ. Monit. Assess. 2018, 190, 351. [Google Scholar] [CrossRef] [PubMed]
  52. Emamgolizadeh, S.; Bateni, S.; Shahsavani, D.; Ashrafi, T.; Ghorbani, H. Estimation of soil cation exchange capacity using Genetic Expression Programming (GEP) and Multivariate Adaptive Regression Splines (MARS). J. Hydrol. 2015, 529, 1590–1600. [Google Scholar] [CrossRef]
  53. Nazari, A.; Torgal, F.P. Modeling the compressive strength of geopolymeric binders by gene expression programming-GEP. Expert Syst. Appl. 2013, 40, 5427–5438. [Google Scholar] [CrossRef] [Green Version]
  54. Babanajad, S.K.; Gandomi, A.H.; Alavi, A.H. New prediction models for concrete ultimate strength under true-triaxial stress states: An evolutionary approach. Adv. Eng. Softw. 2017, 110, 55–68. [Google Scholar] [CrossRef]
  55. Alavi, A.H.; Gandomi, A.H.; Nejad, H.C.; Mollahasani, A.; Rashed, A. Design equations for prediction of pressuremeter soil deformation moduli utilizing expression programming systems. Neural Comput. Appl. 2013, 23, 1771–1786. [Google Scholar] [CrossRef]
  56. Gandomi, A.; Alavi, A.H.; MirzaHosseini, M.R.; Nejad, F.M. Nonlinear Genetic-Based Models for Prediction of Flow Number of Asphalt Mixtures. J. Mater. Civ. Eng. 2011, 23, 248–263. [Google Scholar] [CrossRef]
  57. Iqbal, M.F.; Liu, Q.-F.; Azim, I.; Zhu, X.; Yang, J.; Javed, M.F.; Rauf, M. Prediction of mechanical properties of green concrete incorporating waste foundry sand based on gene expression programming. J. Hazard. Mater. 2020, 384, 121322. [Google Scholar] [CrossRef]
  58. Bardhan, A.; Gokceoglu, C.; Burman, A.; Samui, P.; Asteris, P.G. Efficient computational techniques for predicting the California bearing ratio of soil in soaked conditions. Eng. Geol. 2021, 291, 106239. [Google Scholar] [CrossRef]
  59. Kaloop, M.R.; Bardhan, A.; Kardani, N.; Samui, P.; Hu, J.W.; Ramzy, A. Novel application of adaptive swarm intelligence techniques coupled with adaptive network-based fuzzy inference system in predicting photovoltaic power. Renew. Sustain. Energy Rev. 2021, 148, 111315. [Google Scholar] [CrossRef]
  60. Azim, I.; Yang, J.; Javed, M.F.; Iqbal, M.F.; Mahmood, Z.; Wang, F.; Liu, Q.-F. Prediction model for compressive arch action capacity of RC frame structures under column removal scenario using gene expression programming. Structures 2020, 25, 212–228. [Google Scholar] [CrossRef]
  61. Azim, I.; Yang, J.; Iqbal, M.F.; Mahmood, Z.; Javed, M.F.; Wang, F.; Liu, Q.-F. Prediction of Catenary Action Capacity of RC Beam-Column Substructures under a Missing Column Scenario Using Evolutionary Algorithm. KSCE J. Civ. Eng. 2021, 25, 891–905. [Google Scholar] [CrossRef]
  62. Liu, Q.-F.; Iqbal, M.F.; Yang, J.; Lu, X.-Y.; Zhang, P.; Rauf, M. Prediction of chloride diffusivity in concrete using artificial neural network: Modelling and performance evaluation. Constr. Build. Mater. 2020, 268, 121082. [Google Scholar] [CrossRef]
  63. Hanandeh, S.; Ardah, A.; Abu-Farsakh, M. Using artificial neural network and genetics algorithm to estimate the resilient modulus for stabilized subgrade and propose new empirical formula. Transp. Geotech. 2020, 24, 100358. [Google Scholar] [CrossRef]
  64. Zha, F.; Pan, D.; Xu, L.; Kang, B.; Yang, C.; Chu, C. Investigations on Engineering Properties of Solidified/Stabilized Pb-Contaminated Soil Based on Alkaline Residue. Adv. Civ. Eng. 2018, 2018, 8595419. [Google Scholar] [CrossRef]
  65. Fasihnikoutalab, M.H.; Asadi, A.; Unluer, C.; Huat, B.K.; Ball, R.J.; Pourakbar, S. Utilization of Alkali-Activated Olivine in Soil Stabilization and the Effect of Carbonation on Unconfined Compressive Strength and Microstructure. J. Mater. Civ. Eng. 2017, 29, 06017002. [Google Scholar] [CrossRef]
Figure 1. Frequency histograms of the input and output parameters: (a) fly ash dosage (%), (b) alkali concentration (N), (c) curing age (days), (d) UCSkaolin, and (e) UCSBC.
Figure 1. Frequency histograms of the input and output parameters: (a) fly ash dosage (%), (b) alkali concentration (N), (c) curing age (days), (d) UCSkaolin, and (e) UCSBC.
Materials 15 04025 g001
Figure 2. Comparison of normalized averaged MSE and correlation for the developed MEP model for kaolin soil.
Figure 2. Comparison of normalized averaged MSE and correlation for the developed MEP model for kaolin soil.
Materials 15 04025 g002
Figure 3. Comparison of normalized averaged MSE and correlation for the developed MEP model for BC soil.
Figure 3. Comparison of normalized averaged MSE and correlation for the developed MEP model for BC soil.
Materials 15 04025 g003
Figure 4. The variation in UCSkaolin and UCSBC with curing period and concentration of alkali after a 28-day curing period.
Figure 4. The variation in UCSkaolin and UCSBC with curing period and concentration of alkali after a 28-day curing period.
Materials 15 04025 g004
Figure 5. The UCS variation of both the soils with the concentration of alkali and fly ash dosage.
Figure 5. The UCS variation of both the soils with the concentration of alkali and fly ash dosage.
Materials 15 04025 g005
Figure 6. Comparison of experimental and predicted results to evaluate the UCS in the case of (a) kaolin soil and (b) BC soil.
Figure 6. Comparison of experimental and predicted results to evaluate the UCS in the case of (a) kaolin soil and (b) BC soil.
Materials 15 04025 g006
Figure 7. Error analysis of the developed models to evaluate the UCS in the case of (a) kaolin soil and (b) BC soil.
Figure 7. Error analysis of the developed models to evaluate the UCS in the case of (a) kaolin soil and (b) BC soil.
Materials 15 04025 g007
Figure 8. Tracing of experimental results by predicted values to evaluate the UCS in the case of (a) kaolin soil and (b) BC soil.
Figure 8. Tracing of experimental results by predicted values to evaluate the UCS in the case of (a) kaolin soil and (b) BC soil.
Materials 15 04025 g008
Figure 9. Sensitivity analysis of the developed MEP models: (a) kaolin soil; (b) BC soil.
Figure 9. Sensitivity analysis of the developed MEP models: (a) kaolin soil; (b) BC soil.
Materials 15 04025 g009
Figure 10. Parametric study of input variables for kaolin and BC soil MEP models: (a,d) fly ash dosage, (b,e) alkali concentration, and (c,f) curing period.
Figure 10. Parametric study of input variables for kaolin and BC soil MEP models: (a,d) fly ash dosage, (b,e) alkali concentration, and (c,f) curing period.
Materials 15 04025 g010
Table 1. Physical properties of kaolin and BC soils used in the current study.
Table 1. Physical properties of kaolin and BC soils used in the current study.
PropertyKaolin SoilBC Soil
Specific gravity2.562.65
pH7.37.1
USCS classificationCHCH
Liquid limit (%)4162
Plasticity index (%)1928
Optimum moisture content (%)1723
Maximum dry density (g/cc)1.811.67
Table 2. Chemical composition of fly ash in the current study.
Table 2. Chemical composition of fly ash in the current study.
Chemical ConstituentsValue (%)
Silica (SiO2)62.9
Alumina (Al2O3)21.7
Ferric oxide (Fe2O3)4.5
Calcium oxide (CaO)6.8
Magnesia (MgO)1.08
Titanium (TiO2)0.06
Potash (K2O)0.04
Sulfur (SO3)0.7
Loss on ignition2.21
Table 3. Experimental database of the input parameters and output parameters in the current study.
Table 3. Experimental database of the input parameters and output parameters in the current study.
S. No.Fly Ash Dosage (%)Alkali Concentration (N)Curing Age (Days)UCSBC (kPa)UCSkaolin (kPa)
1001280255
2001271261
3001269259
4001286275
5001278268
6001288277
7007272262
8007274264
9007281270
10007286275
11007289278
12007280269
130014300265
140014280262
150014298279
160014285266
170014301281
180014296277
190028310270
200028286267
210028296277
220028308288
230028301281
240028296276
25011267236
26011260231
37820414631851
37920428729998
38020428732987
38120428750996
382204287451040
383204287321012
384204287401004
Table 4. Statistical description of input and output parameters used for MEP modeling.
Table 4. Statistical description of input and output parameters used for MEP modeling.
Fly Ash Dosage (%)Alkali Concentration (N)Curing Age (Days)UCSBC (kPa)UCSkaolin (kPa)
Minimum00144119
Maximum204287501040
Mean11.251.7512.5379.51369.06
Median12.51.510.5365325.5
SD7.401.4810.06109.14183.30
Kurtosis−1.1537−1.1537−1.14271.13742.2691
Skewness−0.43640.43640.50250.86671.4534
Table 5. Pearson correlation coefficient values for the input parameters and the UCS of alkali-contaminated soils.
Table 5. Pearson correlation coefficient values for the input parameters and the UCS of alkali-contaminated soils.
Fly Ash Dosage (%)Alkali Concentration (N)Curing (Days)UCSkaolin,BC (kPa)
Fly ash dosage (%)1
Alkali concentration (N)01
Curing age (days)001
UCSkaolin (kPa)0.5899060.5083030.1851891
UCSBC (kPa)0.7248090.2704960.3219861
Table 6. Parameter setting for MEP algorithm settings for strength prediction of fly-ash-treated alkali-contaminated soils.
Table 6. Parameter setting for MEP algorithm settings for strength prediction of fly-ash-treated alkali-contaminated soils.
ParametersKaolin SoilBC Soil
Number of subpopulations20100
Subpopulation size10002000
Code length10080
Crossover probability0.90.9
Crossover typeUniform
Mutation probability0.001
Tournament size2
Operators0.5
Variables0.5
Constants0
Number of generations150
Function set+, −, ×, /
Terminal setProblem input
Replication number10
Error measureMean squared error
Problem typeRegression
SimplifiedYes
Random seed0
Number of runs10
Number of threads1
Table 7. Details of trials undertaken in selecting the best MEP models.
Table 7. Details of trials undertaken in selecting the best MEP models.
MEP TrialNo. of SubpopulationSubpopulation SizeCode LengthNo. of GenerationsTournament SizeR2RAvg. MSETime (min)
Kaolin Soil
11010020100268.5482.7981481
220 69.2883.2374891
370 66.2781.4171772
4100 66.2781.4144363
5200 64.5780.3652096
6100500 77.2087.86293725
7 1000 79.0988.93252148
8 1500 78.8988.82256272
9 2000 80.3489.63248585
10 30 82.6090.882109130
11 50 83.6691.471951220
12 80 87.1993.381527300
13 100 87.6593.621474429
14 150 88.9894.331455667
15 200 88.0093.811315925
16201000 150 87.1993.37155140
17 489.3394.511895106
18 689.5894.651245102
BC Soil
11010020100272.4585.1214,6971
220 278.0288.3310,9801
370 277.8188.2111,1872
4100 277.5688.0795783
5200 276.3987.4098048
6100500 279.1988.99873323
7 1000 280.2689.59848652
8 1500 281.1390.078105100
9 2000 280.8889.938026145
10 30 279.2689.037993190
11 50 278.8088.777256330
12 80 280.5589.756592357
13 100 280.0089.447633393
14 80150293.5496.722220552
15 200292.1996.022638549
1670500100100270.1183.73845090
17 487.6693.635976110
18 690.9795.384400112
Table 8. Performance index values of the final MEP models for alkali-activated soils.
Table 8. Performance index values of the final MEP models for alkali-activated soils.
DatasetPerformance IndexKaolin SoilBC Soil
TrainingR0.937130.95661
RMSE18.27117.151
MAE19.630.0
RSE0.12800.1078
RRMSE0.05430.0564
NSE0.87200.8922
ρ0.02800.02882
TestingR0.900140.96243
RMSE21.98722.995
MAE30.554.7
RSE0.19720.0841
RRMSE0.04580.0441
NSE0.80280.9159
ρ0.02410.0225
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khan, K.; Ashfaq, M.; Iqbal, M.; Khan, M.A.; Amin, M.N.; Shalabi, F.I.; Faraz, M.I.; Jalal, F.E. Multi Expression Programming Model for Strength Prediction of Fly-Ash-Treated Alkali-Contaminated Soils. Materials 2022, 15, 4025. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15114025

AMA Style

Khan K, Ashfaq M, Iqbal M, Khan MA, Amin MN, Shalabi FI, Faraz MI, Jalal FE. Multi Expression Programming Model for Strength Prediction of Fly-Ash-Treated Alkali-Contaminated Soils. Materials. 2022; 15(11):4025. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15114025

Chicago/Turabian Style

Khan, Kaffayatullah, Mohammed Ashfaq, Mudassir Iqbal, Mohsin Ali Khan, Muhammad Nasir Amin, Faisal I. Shalabi, Muhammad Iftikhar Faraz, and Fazal E. Jalal. 2022. "Multi Expression Programming Model for Strength Prediction of Fly-Ash-Treated Alkali-Contaminated Soils" Materials 15, no. 11: 4025. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15114025

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop