Bearing Capacity of Ring Foundations on Anisotropic and Heterogenous Clays: FEA, NGI-ADP, and MARS

Axisymmetric solutions for the bearing capacity of ring foundation resting on anisotropic and heterogenous clays are presented in this paper using finite element analysis (FEA). The NGI-ADP model in PLAXIS FEA, a widely used anisotropic soil model, is adopted to study the stability responses of ring foundations, with special consideration given to the effects of increasing undrained shear strength with the depth. Numerical results are formulated in terms of a dimensionless stability number (bearing capacity ratio) that is a function of three dimensionless input parameters: namely, the ratio of inner and outer radius, the increasing strength gradient ratio, and the anisotropic shear strength ratio. The influence of each dimensionless input parameter on the bearing capacity ratio is investigated using design charts and failure mechanisms, and they are scored by relative importance indexes in multivariate adaptive regression splines (MARS) model—a machine learning approach. A highly accurate equation generated from the MARS model is proposed as an effective tool for engineering practitioners.


Introduction
The use of shallow circular and spudcan footings were popular in the past few decades in supporting axisymmetric structures such as transmission towers, water towers, annular platforms, silos, storage tanks, and chimneys. Nowadays, the more economical ring foundations are often being considered, with the increasing numbers of recent research in this direction.
Several experimental studies were reported in relation to the performance of ring footing in sand (e.g., Ohri et al. 1997;Hataf and Razav 2003) as well as in clay (e.g., Demir et al. 2012;Shalaby 2017). Calculation approaches have been constantly developed, Abstract Axisymmetric solutions for the bearing capacity of ring foundation resting on anisotropic and heterogenous clays are presented in this paper using finite element analysis (FEA). The NGI-ADP model in PLAXIS FEA, a widely used anisotropic soil model, is adopted to study the stability responses of ring foundations, with special consideration given to the effects of increasing undrained shear strength with the depth. Numerical results are formulated in terms of a dimensionless stability number (bearing capacity ratio) that is a function of three dimensionless especially in the area of numerical simulation. For instance, Zhao and Wang (2008) performed the bearing capacity of ring foundations in sandy soil using FLAC; whilst Benmebarek et al. (2012) adopted FLAC and Chavda and Dodagoudar (2019) to estimate N c , N q , N γ considering the effects of the dilation angle, smooth and rough ring footings. Choobbasti et al. (2010) assessed the bearing capacity and settlement of ring footings on homogeneous soils using the finite element analysis. Kumar and Chakraborty (2015) adopted the finite element limit analysis (FELA), Chavda and Dodagoudar (2021) applied FEA, and Gholami and Hosseininia (2017) used method of characteristics to estimate the bearing capacity factors for ring foundations on cohesive-frictional soils. The bearing capacity solution of a ring footing rest on clay was also investigated in Lee et al. (2016), Benmebarek et al. (2017), Khatri and Kumar (2009), Yang et al. (2020), Birid and Choudhury (2021), Keawsawasvong and Lai (2021), and Papadopoulou and Gazetas (2020). Recently, the bearing capacity of a ring foundation on rock mass was studied by Yodsomjai et al. (2021a).
It was generally recognized that undrained shear strength of soils not only increases with depth but also it is directly dependent on the orientation of the major principal stress to the vertical axis or depositional direction (see e.g., Ladd and Degroot 2003;Ladd 1991;Yu and Sloan 1994). In other words, the undrained shear strength of clay exhibits somewhat anisotropic behavior in nature. Recently, some failure criteria were proposed to consider the complex inherent of anisotropic clay such as the Anisotropic Undrained Shear (AUS) model in Krabbenhoft et al. (2019). In addition, the NGI-ADP constitutive model was presented by Grimstad et al. (2012). Although the recent finite element limit analysis is considered as a powerful technique to solve various stability problems in the geotechnical engineering field (e.g., Shiau and Al-Asadi 2020a, b, c, d;Keawsawasvong et al., 2022;Keawsawasvong and Ukritchon, 2017a, 2021Lai et al. 2021a;Yodsomjai et al. 2021b, c;Beygi et al. 2020;Ukritchon and Keawsawasvong, 2017, 2019Tho et al. 2014;Bhattacharya and Sahoo 2017;Bhattacharya 2016Bhattacharya , 2018, the displacement-based finite element analysis is also useful in solving the complex boundary conditions (e.g. Shiau and Yu 2000;Shiau and Smith 2006;Shiau et al. 2006a;b, 2018;Shiau and Watson 2008;Halder and Manna 2020;Pirastehfar et al. 2020;Chakraborty and Goswami, 2021;Chatterjee and Murali Krishna 2021;Keawsawasvong and Ukritchon, 2017b;Hamouma et al. 2021;Huynh et al. 2022a, b;Ukritchon et al. 2017a, b).
The application of the NGI-ADP model has been used to analyze several geotechnical problems such as in Zhang et al. (2021), Li et al. (2021), andLangford et al. (2021) for the deep excavation problem; Li and Zhang (2020) for the pile response problem; Aamodt et al. (2021) for the slope stability problem; and Zhang et al. (2020) for the tunnel stability problem. Moreover, the anisotropy effect of soils was investigated for the works related to passive failure modes of plate anchors as well as foundations by , Keawsawasvong and Lawongkerd (2021), Nguyen et al. (2021), and Lai et al. (2022). Nevertheless, very few studies were reported in the literature in relation to the bearing capacity of ring footings with the combined effects of clay inhomogeneity and strength anisotropy using the NGI-ADP model.
In this paper, the NGI-ADP model is adopted to investigate the bearing capacity factor and failure mechanism of ring foundations resting on anisotropic and heterogeneous clays. The effects of the three dimensionless input parameters, namely the ratio of the inner to outer radius of the ring foundation, the increasing shear strength gradient ratio, and the anisotropic strength ratio of the cohesive clay on the performance of the ring foundation are investigated. The final outcomes are presented using design tables, charts, and an empirical equation. Due to the complex effect of each parameter on the bearing capacity factor, the multivariate adaptive regression splines (MARS) model is used to examine the sensitivity of each parameter using the output results, as well as to build a correlation equation between the multi-input parameters and output results. The tools provided in this paper would contribute to practical designs of ring footings resting on anisotropic and heterogeneous clays. Fig. 1 is the problem definition of a rigid ring footing resting on anisotropic and heterogeneous clay. The soil anisotropy is simulated using the NGI-ADP constitutive model, whilst the heterogeneous clay is studied with the undrained shear strength linearly increasing with depth. The rigid ring footing is subjected to a uniform pressure q and has an external radius r o and an internal radius r i .

Shown in
The main soil parameters of the NGI-ADP model are divided into 2 groups: stiffness and strength. The stiffness parameters include G ur ∕s A u (ratio of unloading/reloading shear modulus over the "active" undrained shear strength), u (Poisson's ratio), and ∕ 1 + s P u ∕s A u , is assumed to be a harmonic mean (Krabbenhoft et al. 2019;Keawsawasvong et al. 2022). Figure 2 shows the failure criterion of the NGI-ADP model in the π-plane. Detailed descriptions of the model are stated in Grimstad et al. (2012), Ukritchon and Boonyatee (2015) and Brinkgreve and Vermeer (2019) and will not be repeated here.
Following the discussion of the anisotropic soil model, heterogeneous clay is expressed by using the three anisotropic undrained shear strengths that are linearly increasing with depth. They are presented in Eqs. (1)-(3).
(1)  2 Failure criterion of the NGI-ADP model in the πplane (After Grimstad et al. 2012) where s A u0 , s DSS u0 , s P u0 are the anisotropic undrained shear strengths at the ground level, ρ is the linear-gradient with a unit of kPa/m (per meter depth), z indicates the depth determined from the ground surface (see Fig. 1). Note that these undrained shear strengths s A u , s DSS u , s P u can be obtained from three modes of shearing, including undrained tests of triaxial compression (for s A u ), triaxial extension (for s P u ), and direct simple shear (for s DSS u ). For typical cohesive soils, undrained shearing strength of triaxial compression, s A u is the largest, followed by direct simple shear, s DSS u , and triaxial extension, s P u , which is the lowest. A dimensionless bearing capacity factor N is defined as the ratio of the uniform pressure q over s A u0 , and it is a function of the three dimensionless design parameters, namely the ratio of the inner and outer radii r i ∕r 0 , the strength gradient ratio m = r 0 ∕s A u0 , and the anisotropic ratio r e = s P u ∕s A u . This is shown in Eq. (4) The finite element model with axial symmetry is adopted to investigate the problem, as shown in Fig. 3, by using Plaxis2D v20 (Brinkgreve and Vermeer 2019). The ring footing is modelled as rigid plates subjecting to a uniform pressure q. The soils are simulated by using 15-noded triangular elements with the NGI-ADP soil model. The model size was carefully chosen so that the overall velocity field would not interfere the boundary and the effect on the produced results can be minimized. This is shown in Fig. 3. The bottom boundary is fixed in x, y-directions, the right-hand boundary is fixed only in the x-direction, the left-hand boundary is an axisymmetric line, and the ground surface is free. The ranges of parameters chosen for the study are: (1) r i ∕r 0 = 0, 0.25, 0.33, 0.5, and 0.75; (2) m = r 0 ∕s A u0 = 0, 1, 2.5, 5, and 15; (3) s P u ∕s A u = 0.4, 0.5, 0.6, 0.7, 0.8, and 1.0. It is to be noted that s P u ∕s A u = 1 indicates an isotropic clay. The range of andm was selected based on work of Lee et al. (2016), Benmebarek et al. (2017). According to D'Ignazio et al. (2017) and Hansen and Clough (1981), typical values of s P u ∕s A u ratios vary from 0.3 to 0.81. Following the recommendation in the FE Plaxis code, the values of s C,TX u ∕s A u = 0.99, 0 ∕s A u = 0.7, u = 0.495 are adopted for all analyses in the paper. The selected E u ∕s A u ratio follows a previous stability analysis using FEA , such that there is very small to none effect of this parameter on the limit load of this stability problem. It should be noted that 0 is initial in situ maximum shear stress and E u is undrained young modulus.

Validations, Results, and Discussions
To improve the confidence in later parametric analyses, numerical results from finite element analysis (FEA) must be compared with published results (Shiau et al. 2014(Shiau et al. , 2016a(Shiau et al. , b, 2017. This is done by comparing with those in Lee et al. (2016) and is presented in Fig. 4 for cases of ring footing on isotropic clay (where the anisotropic ratio s P u ∕s A u = 1) and heterogeneous clay (the increase shear gradient ratio m = 0, 1, 2.5, 5, and 15). The comparison study has shown that the present FEA results of bearing capacity factor N are in good agreement with those in Lee et al. (2016) and that the current model can be further be used to study the anisotropic effect with reasonable confidence.
With the success of model validation, the next task is to study the effects of the ratio of inner and outer radius ( r i ∕r 0 ), the anisotropic ratio ( s P u ∕s A u ), and the shear strength gradient ratio (m) on the bearing capacity factor (N). Figure 5 shows the linear relationships of the anisotropic ratio ( s P u ∕s A u ) on the bearing capacity factor N for the various cases of m and r i /r o . Noting that ( s P u ∕s A u = 1) indicates an isotropic clay, a decrease in anisotropic ratio ( s P u ∕s A u ) means an increased effect of soil anisotropy. For the various m considered in Fig. 5 (m = 0, 1, 2.5, 15), numerical results have shown that a decrease in ( s P u ∕s A u ) results in a decrease in the bearing capacity factor N. They have also shown that the bearing capacity factor N decreases with the increasing r i ∕r 0 . Noting that the larger the r i ∕r 0 , the smaller the footing contact area, it is therefore not surprised to see the decreased bearing capacity owing to the reduced footing area (increasing r i ∕r 0 ). Besides, it was found that the rate of decrease is more pronounced for larger values of the shear strength gradient ratio m.
Using the same results, Fig. 6 shows a nonlinear relationship of the shear strength gradient ratio (m) on the bearing capacity factor N . The greater the value of m, the larger the N is. The rate of increase (the gradient) in N decreases as the ratio of inner and outer radius r i /r o increases (less footing area).
The potential failure mechanisms of ring footings in anisotropic and heterogeneous clay are investigated in Figs. 7, 8 and 9. The influence of the anisotropic ratio ( s P u ∕s A u ) on the potential failure mechanisms of ring footings is presented in Fig. 7 for the case of (m = 5 and r i /r o = 0.5). The results of the failure zone and shear band have indicated an unchanged mechanism despite the increase of the (s P u ∕s A u ). It means that the anisotropic ratio ( s P u ∕s A u ) has a little effect on the failure mechanism. More specifically speaking, it is true for the undrained clay only where the soil frictional angle is zero. The same observation was made by the recent studies of undrained anisotropic clay in Nguyen et al. (2021) and Keawsawasvong (2021). Despite this, it is not yet to be concluded for drained soils with non-zero frictional angles owing to the lack of published literatures.
A further failure mechanism study on the influence of shear gradient ratio (m) is shown in Fig. 8 for the case of ( s P u ∕s A u = 0.5 and r i /r o = 0.5). The plots of failure mechanisms have shown that the overall failure zone is reduced in both horizontal and vertical directions as the value of shear strength gradient ratio (m) increases. The interference effect (overlapping) of the failure zone diminishes as m increases and the smallest failure zone occurs at the largest value of m = 15. For the more interesting study, the effects of inner and outer radius ratio ( r i ∕r 0 ) on the failure mechanism of ring footing are presented in Fig. 9 for the case ( s P u ∕s A u = 0.5 and m = 5). The failure mechanism for a solid circular footing is firstly presented for r i ∕r 0 = 0. A ring foundation forms as r i ∕r 0 > 0, and the overlapping of failure zone diminishes as the value of r i ∕r 0 increases, resulting in a Prandtl type of general failure.

Sensitive Analysis and MARS
The multivariate adaptive regression splines (MARS) model is an automated regression modelling tool. Using several piecewise linear segments (splines) with differing gradients, the relationship between the input variable and output results can be established in multi-dimensions. Recently, the use of the MARS model as the machine learning method in analyzing geotechnical data has become more and more common. Lai et al. (2021a, b) adopted the MARS model to assess the impacts of input design parameters on the output ground movements due to the effects of installing twin caisson foundations. Zhang et al. (2017) proposed an empirical equation to determine In the MARS model, the different splines are connected using a knot representing the end of one spline and the beginning of another. The fitted basic functions (BFs) have better flexibility to the studied model where the bends, thresholds, and other derivations from linear functions are allowed (Zhang, 2019). The basic function can be written as: where x is an input variable, and t is a threshold value. MARS model produces BFs by searching in a stepwise process, of which the knot locations are automatically determined using the adaptive regression algorithm. MARS model is presented by a two-step procedure. The first (or called "forward") step is to provide BFs as well as to find their potential knots to optimize the model performance and fitting accuracy. The second (or "backward") step uses a pruning algorithm to delete the least effective terms, resulting in the generation of an optimal model used for the problem prediction. To build the correlation equation between the input and output variables, the MARS model combines all linear basic functions (BFs) which are described in Eq. (6), where a 0 is the constant, N is the number of BFs, g n is the nth BF, a n is the coefficient of g n . Note that increasing the number of basic functions can increase the accuracy of the MARS model.
The current study uses the MARS model to investigate sensitivity analyses of each input variables (i.e., r i ∕r 0 , s P u ∕s A u and m). The aim is to build a mathematical equation for predicting the N value, considering all effects of input variables. A training data set in the MARS model is selected from all FEA results of bearing capacity factors N. The corresponding 150 design combinations and input parameters are shown in Table 1.
To achieve the best accuracy, the chosen number of basic functions is varied to check the performance of MARS using two criteria of statistical analyses named Mean Squared Error (MSE) and the coefficient of determination (R 2 value). MSE represents the mean square error between the predicted output variables and the real output results. The lower MSE value, the better model can be obtained. The closer the value of R 2 is to 1, the better the linear regression fits the data. R 2 value equals 0.0 means that the model fails to predict real value, whilst 1.0 means that the forecast model is highly reliable.
Shown in Fig. 10, is the variation of MSE and R 2 due to the changes in the number of basic functions. Noting that the MSE decreases sharply as the The sensitivity of each input variable is described by the relative important index (RII), as shown in Fig. 11. The value of RII indicates the weight of impaction. RII of 100% means that the respective input variable has the most important influence on the output results. This is shown for the increasing shear strength gradient ratio (m) in this study. The increasing shear gradient ratio has the most significant influence on the bearing capacity factor N. This is followed by the anisotropic ratio ( s P u ∕s A u ), and the ratio of inner and outer radius ( r i ∕r 0 ) with RII of 50.30% and 50.13% respectively. This finding has suggested that all the three investigated parameters play an important role in the design of ring foundations considering the anisotropic and heterogeneous behavior of cohesive soil. Table 2 shows the basic functions and the mathematical equation generated by the MARS model. They can be written as Eq. (7).
To verify the proposed Eq. (7), a comparison between the bearing capacity factors N from FEA results and those from equation prediction is presented in Fig. 12. Numerical results have shown that the predicted N has an excellent fit to those from FEA-with the high value of R 2 = 99.99%. It can be concluded that Eq. (7) is an effective tool to estimate the bearing capacity of ring footings rest on anisotropic and heterogeneous clay.

Conclusion
Axisymmetric solutions for the bearing capacity of ring foundation resting on anisotropic and heterogenous clays have been successfully investigated in this paper using Plaxis finite element analysis and the NGI-ADP soil model-a widely used anisotropic soil model. The influences of inner and outer radius ratio ( r i ∕r 0 ), anisotropic ratio ( s P u ∕s A u ) and the shear strength gradient ratio (m) on the bearing capacity factor (N) and the failure mechanism of ring footings resting on anisotropic and heterogeneous clay are determined. The following conclusions are drawn based on the study results. • FEA results have shown that the bearing capacity factor N decreases with the increasing inner and outer radius ratio ( r i ∕r 0 ) and the decreasing anisotropic ratio ( s P u ∕s A u ). In contrast, the bearing capacity factor N increases with the increasing shear strength gradient ratio (m).
• Using the MARS model for sensitive analysis of each parameter, it was concluded that the increasing shear gradient ratio (m) is the most influential parameter with a relative importance index RII = 100%. This is followed by the anisotropic ratio ( s P u ∕s A u ) and ratio of inner and outer radius ( r i ∕r 0 ) with RII = 50.30 and 50.13%, respectively. • An accurate equation is proposed with R 2 = 99.99%, which is considered as an effective tool for engineering practitioners to evaluate the bearing capacity of ring foundations in anisotropic and heterogenous clays. Note that the value of is the highest accuracy MARS models can provide according to some previous works , b) Zhang et al. 2017, 2018, 2019Zhou et al. 2021).
Funding Open Access funding enabled and organized by CAUL and its Member Institutions.  Fig. 12 Comparison of results-the finite element analysis results (Plaxis) and the proposed equation Data Availability Statement All data, models, or code that support the findings of this study are available from the corresponding author upon reasonable request.

Conflict of Interest
The authors declare no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.