Shear strength characteristics of binary mixture sand-carpet fibre using experimental study and machine learning

ABSTRACT Mixing carpet fibre in sand offers great potential for enhancing soil properties. In this study, direct shear tests were conducted on two different sands mixed with varying carpet fibre percentages to investigate the effects on soil strength, stiffness, and deformation. Artificial intelligence techniques were used to analyse the data and develop predictive models, including an empirical equation that predicts the shear strength. The results showed that the addition of carpet fibre improved soil properties, with increased strength, stiffness, and reduced deformation. The AI models, including the empirical equation, accurately predicted the mixture's shear strength. Furthermore, this study investigated the importance of each input parameter in predicting the mixture's shear strength. The input parameters are normal stress, void ratio, mean particle size, and the ratio of carpet fibre content to specific gravity. According to the results, normal stress is the most important parameter, and mean particle size is the least important.


Introduction
Fibres have become increasingly significant in geotechnical engineering as a result of their ability to improve the mechanical properties of soil and other geotechnical materials (Almeida et al. 2023).Fibres are substances that are longer than they are wide, and they can be classified into two categories: natural and synthetic.There are many natural fibres found in nature, such as coconut fibre, palm fibre, jute, flax, bamboo, cane, sisal, and barely straw (Hejazi et al. 2012).In contrast, synthetic fibres are manmade fibres, such as polypropylene, polyester, polyethylene, glass, and nylon (Hejazi et al. 2012).The use of fibres in geotechnical engineering is often used to strengthen soil.Tables 1 and 2 list the mechanical and physical properties of natural and synthetic fibres.
Carpet fibre has recently emerged as a potential reinforcement material for geotechnical engineering applications (Gowthaman, Nakashima, and Kawasaki 2018).The use of carpet fibres as a reinforcement material can provide a lowcost, environmentally friendly solution to improve the mechanical properties of soil and other geotechnical materials (Rahman, Siddiqua, and Cherian 2022).Carpet fibre falls under the category of synthetic or man-made fibres (Oksman 2000).Carpet fibres are non-degradable and are made to endure overuse and are resistant to mechanical wear which also makes them difficult to recycle.Carpet production requires large amounts of chemicals, the majority of which are derived from non-renewable petroleum.These chemicals may also be released during carpet disposal emitting harmful volatile organic components (VOCs) into the air, making them significant pollutants to the environment (Realff et al. 2005).
In 2004, the worldwide fibre production has reached 64 million tonnes per year (Wang 2006).On the contrary, a vast amount of waste carpet fibre is discarded into landfills annually (Wang 2007).The worldwide waste carpet disposal rate is 4-6 million tons per year (Fashandi, Pakravan, and Latifi 2019) which makes recycling them crucial to contribute to a more sustainable world.Within the engineering field, sustainability is defined as projects that are environmentally friendly causing minimal damage to the environment over their life cycle and also consider factors including safety, durability, efficiency, ecology, environmental protection, carbon reduction, landscape, energy saving, and waste reduction (Chou et al. 2020;Liu 2020;Liu, Chen, and Chou 2019;Nguyen et al. 2023).One of the ancient ideas for reusing these waste carpet fibre materials was in soil reinforcements.Soil reinforcement has gained increased attention throughout the years, and it can be defined as a technique to improve the stability and enhance the engineering characteristics of soil to develop parameters including, shear strength, compressibility, density, and hydraulic conductivity (Kazemian et al. 2010).Fibre reinforced soil has appealing properties including its tensile strength, increased flexibility, greater weather resistance, drainage capacity, its low cost, and its increased peak strength and stiffness properties (Harrison 2006;Kumar and Roy 2022).All these properties make fibres advantageous in many engineering applications.The application of soil reinforcement can be implemented in geo projects including retaining walls, railway embankments, protection of slopes, earthquakes, and soil-foundation engineering (Hejazi et al. 2012).Numerous studies investigated the behaviour of soil reinforced with fibre and found that by combining fibre with soil, the shear strength of the mixture increased due to the tensile resistance of the carpet fibre that intersects with the shear failure and holds the mixture together (Hejazi et al. 2012).
Research by Ghiassian et al. (2004), Sadek et al. (2010), and Nguyen (2022) has demonstrated the importance of inclusion fibres in improving different soil properties.Using carpet waste as a soil reinforcement material, Ghiassian et al. (2004) found that carpet strips were able to significantly increase the strength and ductility of fine sand, particularly at higher aspect ratios (the ratio of fibre length to diameter), and therefore recycling carpet waste suggests a promising and environmentally friendly alternative to landfilling.Sadek et al. (2010), in the same way, emphasized the impact of fibre length on the shear strength of fibre reinforced sand mixtures.Their research discovered that, even with the same aspect ratio, using longer fibres resulted in a significant improvement in the shear strength of the mixture.Nguyen (2022) investigated the mechanical characteristics of clayey sand soils reinforced with Polypropylene (PP) fibres in the same way.They performed compaction, direct shear, and triaxial compression experiments with different fibre lengths (10 mm, 15 mm, and 20 mm) and fibre contents (0%, 0.5%, 1%, and 1.5%).Their findings indicated that fibre reinforcement significantly enhanced the mechanical characteristics of the soil.The study found that fibre lengths of 20 mm and 15 mm, with a fibre content of 1%, produced the best results in terms of maximum dry density and better cohesiveness.The friction angle, on the other hand, did not improve at the same rate.Nguyen (2022) also proposed that fibre reinforcement efficiency could be further maximized by increasing the compaction energy.These studies demonstrate that fibre reinforcements, whether recycled carpet strips, longer fibres, or synthetic PP fibres, could significantly improve the mechanical and structural properties of soils.
The distribution of fibre for reinforcement with sand is a huge area of study that has attracted increasing attention in geotechnical applications.Many researchers dedicated their lab studies to investigate the shear strength of fibres in different distributions including random and discrete distribution.Gray and Ohashi (1983) used the concept of root reinforcement of soil (sand) where the fibres were tested in both perpendicular and inclined orientations to the shear plane and concluded that the greatest shear strength was achieved with fibre orientation of 60 degrees with respect to the shear surface.Diambra et al. (2010) undertook an experimental investigation using polypropylene fibres with Hostun RF sand.The mixture was tested in conventional triaxial compression and extension.It was noticed that the contribution of fibre to the strength was significant in compression but limited in extension confirming that it depends primarily on their orientation with respect to tensile strains.
Despite the importance of fibre carpets in soil, a comprehensive model for predicting their shear strength has not yet been developed.There are several reasons for this, including the multiplicity of effective factors and the non-linearity of the relationships between them.Over the last two decades, artificial intelligence (AI) methods have attracted the attention of engineering and geotechnical researchers.Various AI methods have been successfully applied to slope stability (Baghbani et al. 2022;Bardhan and Samui 2022), soil mechanics (Daghistani andAbuel-Naga 2023 Inazumi et al. 2020;Singh et al. 2021), soil dynamics (Lee, Lee, and Kim 2003;Baghbani et al. 2023a,;Gatto and Montrasio 2023), recycles materials (Baghbani et al. 2023;Daghistani et al. 2023;Saad et al. 2023;Baghbani et al., 2022) and soil cracking (Baghbani et al. 2022;Xu et al. 2022;Zhang et al. 2023).In spite of this, AI methods have not yet been used to predict the shear strength of soil and carpet fibre mixtures and thus form a gap in the history of studies that must be addressed.AI models were used in several studies to predict shear strength of soil.For example, Zhang et al., (2021a) addressed the accurate assessment of undrained shear strength (USS) in soft sensitive clays.It applied novel data-driven methods, including extreme gradient boosting (XGBoost) and random forest (RF), to capture the relationships between USS and various soil parameters.The models were developed using a general approach with five key feature variables and are compared with other machine learning methods.The XGBoost and RF models outperformed the alternatives, providing higher predictive accuracy and feature importance rankings.These models offered promising tools for USS prediction and enhance the interpretability of geotechnical parameters.In another study, Zhang et al., (2021b) explored the application of deep learning (DL) in geotechnical engineering.DL algorithms have powerful feature learning capabilities and are widely used in various fields.Zhang et al., (2021a) discussed the major DL algorithms and their applications in geotechnical engineering.It provided a summary of published literature, reference cases, and adopted DL algorithms in this field.The challenges and future prospects of DL in geotechnical engineering were also highlighted.
After conducting and analysing a series of direct shear tests on the mixture of sand and carpet fibre, two statistical methods and three AI methods were used to predict the shear strength of the mixture.The statistical methods were multiple linear regression (MLR) and partial least squares (PLS), and the AI methods were classification and regression random forest (CRRF), genetic programming (GP), and Genetic Algorithm-Emotional Neural Network (GA-ENN).Among the strengths of this study is the use of both a black box AI method, i.e.CRRF, and a white box AI method, i.e.GP.Methods of white box AI produce equations or trees that readers can use for a database.Based on four inputs of carpet fibre content/specific gravity, normal stress, void ratio, and mean particle size, this study presents a new equation for predicting the shear strength of carpet fibre and sand mixture using the GP method.After reviewing the results of the mathematical models, the best models are presented and the importance of the parameters in different models is discussed.

Material and method
Two types of sand were used, namely L-sand and G-sand, with mean particle sizes (D 50 ) 0.52 and 0.25 mm, respectively.These sands have coefficients of uniformity (C u ) of 2.04 and 1.62, respectively, and coefficients of curvature (C c ) of 0.96.It appears that both sands are poorly graded, and their sieve particle sizes represent in Figure 1.The sieve analysis was conducted according to the Australian standard (AS1774.192003).The specific gravity of the L-sand was 2.57, while that of the G-sand was 3.79.These specific gravity measurements were taken according to the Australian standard (AS1289.3.5.1 2006).In addition, scanning electron microscopes (SEM) were used to look for the sand particle shape.The SEM results based on 50 particles revealed that sands had sphericity of 0.91 of L-sand and 0.88 of G-sand, where one means a complete sphere.The sphericity is measured using Equation (1).
where S is the sphericity, and Re is circle with area equal to projected particle area.
The carpet fibres employed in this study were unused and shredded into lengths ranging from 15 mm to 25 mm, with an average length of 20 mm (Figure 2).The tensile strength and density, as shown in the literature, are 800 MPa and 0.9 g/cm3, respectively (Sriram and Sidhaarth 2022).The carpet fibre content (CFC) in the mixture was measured with Equation (2).
The mixture contents of carpet fibres in the sand were varied at levels of 0%, 0.25%, 0.5%, and 1% of sample dry weight.A direct shear apparatus was used in the experiments, with a mould of 60 × 60 mm.Fifty-two tests were conducted according to the Australian standard (AS1289.6.2.2 2020), where each sample was subjected to four normal stresses of 25, 50, 100, and 200 kPa and a shear rate of 1 mm/min.
In term of mathematical models, five different statistical techniques (MLR, PLS, CRRF, GP, GA-ENN) were considered in machine learning.Here is a brief identification about the using techniques.The MLR (multiple linear regression) is a statistical model that can predict dependent variables based on multiple independent variables.PLS (partial least squares) is another statistical method for handling highly correlated data and reducing dimensionality, CRRF (classification and regression random forest) is an AI black-box method that combine decision trees and random forest to classify and predict variables, GP (genetic programming) is an AI grey-box technique that utilizes natural selection and genetic principles to evolve computer programs and model, and GA-ENN (Genetic Algorithm-Emotional Neural Network) is a new hybrid method that Genetic Algorithm is employed to optimize the architecture and ENN is a modified version of ANN as it includes neural parameters which enhance the network learning process.

Shear strength, dilation, and contraction
As displayed in Figure 3, the shear strength of a sand-carpet fibre mixture can be influenced by the mixture density, CFC, and normal stress.Higher density mixtures tend to have higher shear strength due to increased interlocking of the sample particles which strengthens the mixture.The CFC can enhance shear strength by improving stability and rigidity, which can resist shear deformation.The normal stress also affects the shear strength of the mixture by increasing interlocking between sand particles and carpet fibres, making the mixture more stable and durable.
In a sand-carpet fibre mixture, dilation and contraction during shear can occur due to the changes in the size of the sample when it is subjected to shear forces.As in the dense states, the sand-carpet fibre mixture has a low void ratio, leading the sand particles to move apart and causing the sample to dilate.Conversely, in a loose state, the sample may contract due to the high void ratio within it.In this scenario, the carpet fibres may tend to align themselves in the direction of the applied shear force.According to Figure 4, when the mixture is in its loose state, it exhibits contraction during shear while subjected to normal stresses of 200, 100, and 50 kPa.This contraction phenomenon becomes more pronounced with an increase in the CFC, which indicates that the carpet fibres increase the void in the sample.In terms of shear stress response, the peak and residual shear strengths are very similar in the loose state, and the shear strength of the mixture increases as CFC increases.
In the dense state, under the normal stresses 200, 100, and 50 kPa and as showing in Figure 5.The mixture slightly contracts in the beginning where the sand particle adjusted in the existed void and then dilated during the shear, and it is showing that as the CFC increases, the mixture had more void and takes longer time to dilate comparing with less CFC mixture.In terms of shear strength characteristics, it is evident that the dense state mixture exhibits higher shear strength compared to the loose state.Additionally, as the CFC increases, the mixture tends to show a transition towards behaving more like a loose state, where there is not a significant change in the shear strength.Thus, as the carpet fibre content increases in the mixture, the voids increase, resulting in a looser mixture that experiences more contraction.Furthermore, the dilation and contraction of a sand-carpet fibre mixture during shear are influenced by CFC, mixture density, and normal stresses.According to Figure 6, as the CFC increases in the mixture, it takes longer to reach the maximum shear strength.Also, a loose state takes longer to reach its maximum shear strength than a dense one.Meanwhile, normal stress has a slight influence compared to other factors.The optimal values for these factors will depend on the intended use of the material and the shear loads that it is expected to experience.

Mixture density and compressibility
The sand-carpet fibre mixture was prepared with two densities: a loose state and a dense state, where the carpet fibre content was 0%, 0.25%, 0.5%, and 1% of the sample dry weight.The relative density values (D r ) for the loose state were 9.09%, 7.41%, 5.56%, and 3.23%, respectively.In contrast, for the dense state, the D r values were 90.91%, 92.59%, 94.44%, and 96.77%, respectively.The inclusion of carpet fibres in the sand mixture has a noticeable effect on the dry density and void ratio, as observed in Figure 7.The inclusion of carpet fibres decreases the mixture's density due to their lower specific gravity compared to sand, and they increase the void ratio in the sample by filling the space between sand particles.
The coefficient of volume compressibility (m v ) is an important parameter for evaluating the settlement and consolidation behaviour of the mixture under different loading conditions.The m v decreases with an increase in dry density and a decrease in CFC, as shown in Figure 8.This means that the mixture becomes more resistant to compression as it becomes denser and has less carpet fibre.

Data-driven modelling results
Statistical data from direct shear tests (Table 3), including 52 observations of five variables: shear strength, normal stress, void ratio, mean particle size (D 50 ), and CFC/G s .Mean and standard deviation values are provided for each variable, which are useful for developing an AI model to predict mixture shear strength based on normal stress, void ratio, mean particle size, and CFC/G s .
The experimental database was divided into training (80%) and testing (20%) sets to develop an AI model.Tables 4 and 5 provide statistical information about these sets, and their similarity is important to improve model accuracy.The similar values of means and standard deviations indicate that the training set is representative of the testing set, which suggests that the AI model will perform well on new data.MATLAB (R2021a: The MathWorks Inc., Natick, MA, U.S.A.) was used for all mathematical modelling in this study.

Multiple linear regression (MLR)
After performing a series of MLR modelling, the best model was identified.In Figure 9, the values predicted by the MLR model are compared to the actual values of shear strength in direct shear tests.According to the results, the MLR model has a moderate level of accuracy.
The accuracy of MLR model for predicting shear strength was evaluated using a training database of 41 observations and a testing database of 11 observations.The model's accuracy was measured (Table 6) using several statistical metrics, including MAE, RMSE, RMSLE, and R 2 .The MAE for the testing database was slightly higher at 21.724 compared to 17.369 for the training database.Similarly, the RMSE for the testing database was higher at 22.833 than the RMSE for the training database at 19.785.The R 2 value for the training database was 0.787, indicating that the model explains 78.7% of the data's variability.However, the R 2 value for the testing database was slightly lower at 0.653, suggesting that the model may not perform as well on unseen data as on training data.Overall, the MLR model demonstrated moderate accuracy in predicting shear strength, but its performance on new data may vary.

Partial least squares (PLS)
To find the best PLS model, several simulations were carried out under different conditions.Figure 10 illustrates the results of the best PLS model in terms of values predicted by the PLS model versus actual values obtained from the direct shear test.
Based on the results, PLS model has been found to be more accurate at predicting the shear strength of sand and carpet fibre mixtures.
The performance of a PLS model for predicting the shear strength of sand and carpet fibre mixtures is presented in Table 7, using 41 observations in the training dataset and 11 observations in the testing dataset.The model's accuracy is evaluated using statistical metrics such as MAE, RMSE, RMSLE, and R 2 .The PLS model exhibits a mean absolute error of 10.662 and 15.952 for the training and testing databases, respectively.Additionally, the RMSE values for the training and testing databases are 14.941 and 16.966, respectively.The model's RMSLE values for the training and testing databases are 0.228 and 0.438, respectively, showing smaller errors on a logarithmic scale than the RMSE values.The R 2 values for the training and testing datasets are 0.878 and 0.808, respectively, indicating that the model explains around 80% to 88% of the variance in shear strength values.Overall, the PLS model performs well in estimating shear strength, but its accuracy may vary.For higher accuracy, better models like AI models may be needed.

Classification and regression random forest (CRRF)
The specifications of the best CRRF model based on a thorough evaluation process are presented in Table 8.The CRRF model is a machine learning algorithm that combines the strengths of regression and random forest models, and specific parameters are optimized through a search process.The CRRF trees have a minimum node size of 2 and a minimum son size of 1, which determine the minimum number of observations required to split a node and create a new child node, respectively.The maximum depth of the tree is set at 7, indicating the maximum number of splits.The CRRF forest parameter has a M try value of 4 and a CP value of 0.00001, indicating the number of predictor variables randomly selected at each split and the minimum decrease in the mean squared error required to split a node, respectively.Each tree in the forest is constructed by random sampling with replacement, with a sample size of 41.The final model consists of 200 trees.
Based on the best CRRF model, Figure 11 shows the shear strength values predicted by the CRRF model versus the actual shear strength values obtained by the direct shear test.According to the results, the CRRF model predicted shear strength values more accurately than the two statistical models.
The performance of the CRRF model is presented in Table 9.The CRRF model was trained and tested on datasets with 41 and 11 observations, respectively.The MAE for the testing dataset was 4.857, indicating a difference of approximately 4.8 units between the predicted and actual values.The model's RMSE for the testing dataset was 6.695, which is lower than the RMSE of the training dataset at 9.324, suggesting better performance on unseen data.The RMSLE values for the testing and training datasets were 0.118 and 0.098, respectively, indicating high accuracy in predicting shear strength.
The R 2 values for the testing and training datasets were 0.970 and 0.953, respectively, demonstrating a good fit between the model and the data.Overall, the CRRF model accurately predicts shear strength, but slightly higher error metrics for the training dataset suggest overfitting.

Genetic programming
After having successful results of using a black box AI method (CRRF), a white box AI method (GP), was used.Figure 12 presents a comparison of predicted and actual shear strength values for the best GP model, which was determined after conducting numerous GP modelling analyses.According to the obtained results, the GP model is more accurate than the CRRF model.
The performance of the GP model in predicting the shear strength of sand and carpet fibre mixtures is presented in Table 10, including both training and testing databases, evaluated using MAE, RMSE, RMSLE, and R 2 .The training database comprised 41 observations, while the testing database had 11 observations.The GP model's MAE was 5.294 for the training dataset and 3.495 for the testing dataset, indicating higher accuracy in predicting shear strength in the testing dataset.Similarly, the RMSE of the GP model was 7.270 in the training dataset and 4.335 in the testing dataset, indicating higher accuracy in the testing dataset.The RMSLE for the training dataset was 0.081, while for the testing dataset, it was 0.110.The R 2 values for the training and testing datasets were 0.971 and 0.987, respectively, indicating high predictive ability of the GP model for the shear strength of the sand and carpet fibre mixture.Overall, the GP model demonstrates good predictive ability for the shear strength of the sand and carpet fibre mixture.Through normalization of the parameter values, the optimal equation (Eq. 3) for the genetic model is obtained.
where r 1 , r 2 , r 3 , and r 4 are constants and CFC is carpet fibre content, e is void ratio, G s is specific gravity, D 50 is mean particle size, and σ n is normal stress.Table 11 shows the values of constants in Equation 3.

GA-ENN (genetic algorithm-emotional neural network)
Lotfi and Akbarzadeh (2014) have recently developed a limbicbased Emotional Neural Network (ENN), which adopts a single-layer structure inspired by the emotional processes of the brain.Unlike Artificial Neural Networks (ANNs) that mimic biological neurons, ENNs are designed to simulate the interactions between four specific neural areas of the emotional brain: the thalamus, sensory cortex, orbitofrontal cortex (OFC), and amygdala.
In this study, a hybrid method (GA-ENN) was used to predict shear strength of sand-carpet fibre mixtures.Following several analyses, the optimal GA-ENN model has been identified and introduced.Table 12 provides an overview of the specifications of the best proposed model.
A comparison between the predicted and actual shear strength values is illustrated in Figure 13 for the best GA-ENN model.This model was determined after conducting extensive GA-ENN modelling analyses.The results obtained demonstrate that the GA-ENN model outperforms all previous models in terms of accuracy.
Table 13 provides an assessment of the GA-ENN model's performance in predicting the shear strength of sand and carpet fibre mixtures, considering both the training and testing databases.Evaluation metrics such as MAE, RMSE, RMSLE, and R 2 were utilized.The training database consisted of 41 observations, while the testing database contained 11 observations.For the training dataset, the GA-ENN model achieved an MAE of 2.639 and for the testing dataset, the MAE was 3.118.These results indicate a higher level of accuracy in predicting shear strength for the training dataset.Similarly, the RMSE values were 3.271 for the training dataset and 3.342 for the testing dataset, again showcasing higher accuracy in the training dataset.The RMSLE value for the training dataset was 0.063, while for the testing dataset, it was 0.101.This indicates the level of accuracy in predicting the logarithmic differences between the predicted and actual shear strength values, with lower values suggesting higher accuracy.Furthermore, the R 2 values for the training and testing datasets were 0.994 and 0.993, respectively.These high R 2 values indicate the strong predictive ability of the GA-ENN model for estimating the shear strength of the sand and carpet fibre mixture.

Discussion
The shear strength of a sand-carpet fibre mixture can be impacted by the mixture's density, and confining pressure, as well as the carpet's fibre content, aspect ratio, and orientation.For a mixture to perform effectively in a variety of situations, it is essential to understand how these factors influence its performance.

Carpet's fibre content, aspect ratio, and orientation
The type and amount of carpet fibre present in the mixtures have a significant influence on its shear strength.The shear strength of the mixture increases as the CFC increases.Fibres, as reinforcement, give the mixture load-bearing capacity as well as resistance to shear stress.However, there is a certain limit to the amount of carpet fibre that can be added to the mixture beyond which no further increase in shear strength will be observed.By neglecting the influence of dry density, and as shown in Figure 14, the gap between loose and dense varies less with increasing carpet fibre content.Furthermore, the figure indicates that CFCs have a greater impact on loose mixtures than on dense mixtures.Shear strength is also affected by the aspect ratio of the fibres (the ratio of the length to the diameter) as well as their orientation.As the aspect ratio of the fibre increases, the surface area of the fibre also increases, and thus, aspect ratio must be considered along with fibre quantity, orientation, and sand particle morphology.Orientation vertically is associated with a higher shear strength than orientation horizontally.It is, however, difficult to manage the carpet fibre orientation in real-life situations.Figure 15 illustrates a schematic diagram of sand reinforced with carpet fibres at different particle shapes, carpet fibre orientations, and aspect ratios.

Active lateral pressure
Lateral earth pressure in a sand-carpet fibre mixture can be affected by carpet fibre content, mixture density, and confining   pressure (Figure 16).A carpet's fibre content can have an impact on the lateral earth pressure in several ways.By increasing the carpet fibre content, the interlocking mechanism between the sand and fibre particles can be improved, resulting in an increase in shear strength.Consequently, the active lateral earth pressure increases.However, the inclusion of carpet fibre to the mixture reduces the density and increases the void ratio between the sand particle which affects and reduces the lateral pressure.As carpet fibre has a low specific gravity in comparison to sand, the density of the mixture decreases as we add more carpet fibre, and therefore the active lateral pressure decreases.The density of the mixture has an impact on the lateral earth pressure.By increasing the total weight of the soil mass, a higher mixture density can increase the lateral earth pressure.In contrast, a lower mixture density may result in a lower lateral earth pressure.Normal stress (confined pressure) can also affect lateral earth pressure.If the confining pressure increases, the lateral earth pressure will increase, whereas if the confining pressure decreases, the lateral earth pressure will decrease.In most cases, confined pressures are caused by the weight of soil and the normal load applied, which increases the shear strength and active lateral pressure.

Comparison of mathematical models
The performance of five distinct mathematical models in predicting the shear strength of sand and carpet fibre mixtures, namely MLR, PLS, CRRF, GP, and GA-ENN, is compared in Table 14.The evaluation was based on various statistical measures computed on both training and testing datasets.The findings showed that GA-ENN outperformed the other models, exhibiting the lowest MAE on the testing database.GP and CRRF closely followed GA-ENN in performance.All AI models including GA-ENN, CRRF, and GP had the lowest MAE values on the training database, while MLR and PLS recorded considerably higher values.Moreover, GA-ENN demonstrated the lowest RMSE values on both the training and testing databases, followed by GP and CRRF.On the other hand, MLR and PLS models recorded relatively higher RMSE values.The logarithmic predictions' accuracy was assessed using RMSLE, where GA-ENN, CRRF and GP models showed the lowest values on both training and testing databases.All models exhibited relatively high R 2 values on both datasets, with GA-ENN, CRRF, and GP showing the highest R 2 values.The findings revealed that the GA-ENN model performed the best on both the training and testing databases.Also, MLR and PLS models may not be the most appropriate models for predicting shear strength of this mixture.

Importance of parameters
The significance of input parameters in mathematical models is a crucial aspect of research in this domain.To assess their    relative importance, a sensitivity analysis was conducted by varying one parameter by ± 100%, while keeping the other parameters constant, and measuring the introduced error.This procedure was applied to all the models developed in this study, and the outcomes are presented in Figure 17, where larger errors indicate greater sensitivity of the model to that particular parameter, and hence higher importance of the parameter.Mean increase error parameter quantifies the average increase in error resulting from sensitivity analysis for each individual parameter.For all models, Table 15 shows the results of ranking the importance of parameters.The results indicate that, generally, normal stress parameter is the most significant, while mean particle size (D 50 ) parameter is the least significant.According to the results, follow normal stress, CFC/G s , and void ratio, respectively, are the most important parameters in four models MLR, PLS, CRRF, and GA-ENN.

Limitations
While the study provides valuable insights into the effects of adding carpet fibre on soil properties and develops predictive models, it is important to acknowledge some limitations.Firstly, the study focused on the effects of carpet fibre on soil properties in direct shear tests, which may not fully capture all potential geotechnical applications or scenarios.Further research is needed to explore a broader range of soil tests and field conditions.Secondly, the study examined the effects of carpet fibre on two specific types of sand, and it is uncertain how applicable the findings are to other soil types.The generalizability of the results to different soil compositions and conditions needs to be further investigated.Thirdly, the study mainly focused on the immediate effects of carpet fibre addition on soil properties.The long-term performance and durability of the modified soil need to be assessed to determine its sustainability and effectiveness over time.Fourthly, while the AI models developed in the study exhibit high predictive accuracy for shear strength, they are based on simplified assumptions and may not fully capture the complexity and variability of real-world geotechnical systems.Further refinement and validation of the models are necessary.Next, although the study examines the importance of input parameters in predicting soil behaviour, it only considers a specific set of parameters.Other factors, such as environmental conditions or additional soil properties, may also influence the effectiveness of carpet fibre and should be considered in future research.Finally, the current study acknowledged a limitation regarding the limited number of datasets utilized.This limitation restricts the scope and applicability of the findings.To enhance the robustness and generalizability of the study, it is recommended to utilize a larger and more comprehensive database for future research.By expanding the dataset, a more comprehensive understanding of the subject matter can be obtained, leading to more reliable and widely applicable conclusions.

Conclusion
The inclusion of carpet fibre in sand is a crucial step in improving the strength and stability of the mixture.The experiments in this study were carried out under various conditions, with varying percentages of carpet fibre (0%, 0.25%, 0.5%, and 1%) mixed with two different types of sand, to assess the effectiveness of using carpet fibre in sand.For the first time, two black-box and one grey-box AI models, namely genetic algorithm-emotional neural network (GA-ENN), classification and regression random forest (CRRF), and genetic programming (GP), and two statistical models, namely multiple linear regression (MLR) and partial least squares (PLS), were used to analyse the data obtained from the experiments and to develop predictive models.The following are the main results in this study: • The shear strength of sand-carpet fibre mixtures is influenced by the mixture's density and confining pressure, as well as the carpet's fibre content, aspect ratio, and orientation.It is important to understand these factors to predict the performance of the mixtures in various applications and to optimize the mixture design for specific requirements.
• The inclusion of carpet fibre to sand decreases the mixture density, and as known, low density is equal to low shear strength.In contrast, the fibre reinforcement increases the shear strength by increasing the interlocking between the sand particles.Overall, the experiment revealed that loose reinforced mixture at specific reinforcement percentage, and in comparison, with dense pure sand, could produce higher shear strength.This depends on the carpet fibre content and instalment orientation.• Equation was created using GP model to predict the shear strength of the mixture using the following inputs: normal stress, void ratio, mean particle size, and CFC/G s .
The accuracy of the model was above 97%.• GA-ENN, CRRF, and GP were the most accurate models in predicting the mixture shear strength, with higher R 2 values and lower RMSE and RMSLE.In contrast, MLR and PLS models are less suitable.• Normal stress is identified as the most significant parameter in all models.Conversely, mean particle size was identified as the least important parameter in all models.

Figure 1 .
Figure 1.Sieve analysis of the used sand.Figure 2. Carpet fibres.

Figure 3 .
Figure 3. Shear strength versus dry density of the mixture at different density and normal stresses.

Figure 4 .
Figure 4. Shear stress and vertical deformation versus horizontal displacement for loose mixture at different normal stresses, a) 200 kPa, b) 100 kPa, and c) 50 kPa.

Figure 7 .
Figure 7. Unit weight and void ratio versus CFC.

Figure 8 .
Figure 8. Coefficient of volume compressibility versus dry density of mixture with different CFC.

Figure 9 .Figure 10 .
Figure 9.The results of MLR for (a) training and (b) testing databases.

Figure 11 .
Figure 11.The results of CRRF model for (a) training and (b) testing databases.

Figure 12 .
Figure 12.The results of GP model for (a) training and (b) testing databases.

Figure 13 .
Figure 13.The results of GA-ENN model for (a) training and (b) testing databases.

Figure 14 .
Figure 14.Shear strength versus carpet fibre content at loose state and dense state.

Figure 16 .
Figure16.Lateral earth pressure versus dry density of sand-carpet fibre mixture at different densities and normal stresses.

Figure 15 .
Figure 15.Schematic diagram of sand reinforced with carpet fibre, showing different particle shapes, carpet fibre orientations, and aspect ratios.

Table 3 .
The statistical information of experimental database.

Table 4 .
Statistical information about the training database.

Table 5 .
Statistical information about the testing database.

Table 6 .
The performance of the MLR model to predict shear strength.

Table 7 .
The performance of the PLS model to predict shear strength.

Table 9 .
The performance of the CRRF model to predict shear strength.

Table 8 .
The specifications of the best CRRF.

Table 12 .
Characteristics of developed GA-ENN model.

Table 10 .
The performance of the GP model to predict shear strength.

Table 13 .
The performance of the GA-ENN model to predict shear strength.

Table 14 .
The results of five mathematical models for training and testing databases.

Table 15 .
The results of variable importance for all AI models.