SciELO - Scientific Electronic Library Online

 
vol.40 issue3Adsorption Thermodynamics of Cr(VI) Removal by Using Agro-Industrial Waste of Oil Palm Bagasse and Plantain Peels author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Ingeniería e Investigación

Print version ISSN 0120-5609

Ing. Investig. vol.40 no.3 Bogotá Sep./Dec. 2020  Epub Jan 25, 2021

https://doi.org/10.15446/ing.investig.v40n3.83170 

Original articles

Bearing Capacity and Settlement Prediction of Multi-Edge Skirted Footings Resting on Sand

Capacidad de Carga y Predicción de Asentamiento de Zapatas Bordeadas de Bordes Múltiples que Descansan Sobre Arena

Tammineni Gnananandarao1 

Vishwas Nandkishor Khatri2 

Rakesh Kumar Dutta3 
http://orcid.org/0000-0002-4611-9950

1Civil Engineering Andhra University, M.E. Andhra University, Ph.D. National Institute of Technology Hamirpur, India. Affiliation: Research Scholar, Civil Engineering Department, NIT Hamirpur, India. E-mail: anandrcwing@gmail.com

2Civil Engineering Shivaji University, M.Tech. Indian Institute of Science, Bangaluru, Ph.D. Indian Institute of Science, Bengaluru, India. Affiliation: Assistant professor, Department of Civil Engineering, Indian Institute of Technology, Dhanbad, India. Email: vishuiisc@gmail.com

3Civil Engineering K.N.I.T. Sultanpur, M.Tech. UniversityofRookee, Ph.D. Indian Institute of Technology Delhi, India. Affiliation: Professor, Department of Civil Engineering, NIT Hamirpur, India. E-mail: rakeshkdutta@gmail.com


ABSTRACT

This paper presents the application of artificial neural networks (ANN) and multivariable regression analysis (MRA) to predict the bearing capacity and the settlement of multi-edge skirted footings on sand. Respectively, these parameters are defined in terms of the bearing capacity ratio (BCR) of skirted to unskirted footing and the settlement reduction factor (SRF), the ratio of the difference in settlement of unskirted and skirted footing to the settlement of unskirted footing at a given pressure. The model equations for the prediction of the BCR and the SRF of the regular shaped footing were first developed using the available data collected from the literature. These equations were later modified to predict the BCR and the SRF of the multi-edge skirted footing, for which the data were generated by conducting a small scale laboratory test. The input parameters chosen to develop ANN models were the angle of internal friction (𝜙) and skirt depth (D s ) to the width of the footing (B) ratio for the prediction of the BCR; as for the SRF one additional input parameter was considered: normal stress (σ). The architecture for the developed ANN models was 2-2-1 and 3-2-1 for the BCR and the SRF, respectively. The R for the multi-edge skirted footings was in the range of 0,940-0,977 for the ANN model and 0,827-0,934 for the regression analysis. Similarly, the R for the SRF prediction might have been 0,913-0,985 for the ANN model and 0,739-0,932 for the regression analysis. It was revealed that the predicted BCR and SRF for the multi-edge skirted footings with the use of ANN is superior to MRA. Furthermore, the results of the sensitivity analysis indicate that both the BCR and the SRF of the multi-edge skirted footings are mostly affected by skirt depth, followed by the friction angle of the sand.

Keywords: square/circular skirted footings; multi-edged skirted footings; bearing capacity ratio; settlement reduction factor; artificial neural networks; multivariable regression analysis

RESUMEN

Este documento presenta la aplicación de redes neuronales artificiales (ANN) y el análisis de regresión multivariable (MRA) para predecir la capacidad de carga y el asentamiento de las zapatas bordeadas de bordes múltiples en arena. Estos parámetros se definen, respectivamente, en términos de la relación de capacidad de carga (BCR) de carga de la zapata con zócalo y sin zócalo y el factor de reducción de asentamiento (SRF), la razón de la diferencia en la solución de zócalo sin zócalo y zapatas bordeadas para el asentamiento de zapatas sin falda a una presión determinada. Las ecuaciones modelo para predecir la BCR y el SRF de la zapata de forma regular se desarrollaron primero utilizando los datos disponibles recopilados de la literatura. Estas ecuaciones se modificaron posteriormente para predecir la BCR y el SRF de la zapata bordeada de bordes multiples, para la cual se generaron los datos mediante la realización de una prueba de laboratorio a pequeña escala. Los parámetros de entrada elegidos para desarrollar modelos ANN fueron el ángulo de fricción interna (𝜙), la profundidad del faldón (D s ) al ancho de la relación de zapata (B) para la predicción del BCR; en cuanto al SRF, se consideró un parámetro de entrada adicional: la tensión normal (σ). La arquitectura para los modelos ANN desarrollados fue 2-2-1 y 3-2-1 para la BCR y el SRF, respectivamente. El R para las zapatas bordeadas de bordes múltiples estuvo en el rango de 0,940-0,977 para el modelo ANN y 0,827-0,934 para el análisis de regresión. De manera similar, el R para la predicción del SRF pudo haber sido de 0,913-0,985 para el modelo ANN y 0,739-0,932 para el análisis de regresión. Se reveló que la BCR predicha y el SRF para las zapatas con borde de múltiples bordes con el uso de ANN es superior al MRA. Además, los resultados del análisis de sensibilidad indican que tanto el BCR como el SRF de las zapatas bordeadas de bordes múltiples se ven más afectados por la profundidad de la falda, seguida del ángulo de fricción de la arena.

Palabras clave: zapatas de zócalo cuadradas/circulares; zapatas bordeadas de bordes múltiples; relación de capacidad de carga; factor de reducción de liquidación; redes neuronales artificiales; análisis de regresión multivariable

Introduction

The prediction of the bearing capacity and the footing settlement with a reasonable accuracy is required in the field of foundation design and affects the overall economy of a project. For this purpose, the bearing capacity equations for conventional footings such as strip, square, circular, and rectangular were well established and are available in the literature. However, different, unconventional geometries for shallow footings were sometimes required for the economy and architectural matters. Numerical analysis using FLAC 3D software was carried out by Ghazavi and Mokhtari (2008) to observe the failure behavior of the sand beneath multi-edge shallow footings. The laboratory tests were performed by Davarci, B., Ornek, M., and Turedi, Y. (2014) on multi-edge footings and reported that the performance of such footings was better than that of square ones of the same width. Furthermore, for these types of footings, no bearing capacity and settlement equations were available in the literature. Hence, the researchers had to resort to experimental research, but generating this kind of data is very expensive. An alternative approach could be building a mathematical model to understand the relationships between the various parameters by calibrating and fitting the generated experimental data. To this effect, the power of ANN to store, learn and capture the complicated relationships between multiple parameters without any prior assumptions is the best choice in determining the bearing capacity ratio and the settlement reduction factor. Hence, the paper presents the application of ANN to model such parameters in unconventional multi-edge skirted footings on the sand.

Background

There have been several studies with ANN in the geotechnical engineering, such as those related to bearing capacity and the settlement of regular shaped footings in different mediums (Kalinli, Acar, and Gunduz, 2011; Marto, Hajihassani, and Momeni, 2014; Ziaee, Sadrossadat, Alavi, and Shadmehri, 2015; Nazir et al., 2015a; Rezaei, Nazir, and Momeni, 2016; Momeni, Armaghani, Fatemi, and Nazir, 2017; Khudier, 2018). The bearing capacity of the footing on sand was predicted by Kalinli et al. (2011) using ANN based on 97 datasets by varying the footing width, embedment depth, geometry, unit weight of sand, and the friction angle of the cohesionless soil. On similar lines with Kalinli et al. (2011), other researchers have explored the application of ANN on the prediction of the bearing capacity of footings resting on sand and rock. The list of these references with the scope of their work, type of material, type of data collected, settlement to width ratio, data set, and input parameter is shown in Table 1.

Table 1 Dataset and parameters varied in the literature for the development of ANN model 

Source: Authors

Note that, for the references listed in Table 1, the dataset in various cases ranged between 75 and 150. Similarly, some papers predicted the ultimate bearing capacity and settlement of the piles using experimental data with help from different soft computing techniques (Nazir, Momeni, Marsono, and Sohaie, 2013; Nazir et al. 2015b; Momeni et al. 2015a; Harandizadeh, Armaghani, and Khari, 2019; Chen etal. 2020; Khari et al. 2020; Yong et al. 2020). Furthermore, the results obtained from all the studies above indicated that ANN-based predictive models could be satisfactorily used in predicting the bearing capacity and the settlement of regular shaped footings. In the present paper, the data for regular shaped (square and circular) skirted footings were collected from the published literature, whereas the data for the multi-edge (T, Plus, Double box) skirted footings were generated through experimentation in the laboratory. The data was used during training, testing, and the predictive phase of ANN models. Input variables for the bearing capacity ratio were skirt depth to width ratio of the footing and friction angle of the sand for ANN modeling. Similarly, to model the settlement reduction factors, the considered input variables were skirt depth to width ratio of the footing, friction angle of the sand, and normal stress. The outputs for these ANN models were bearing capacity ratio and the settlement reduction factor, respectively.

Experimental Materials and Methods

As mentioned earlier, the data related to the bearing capacity and the settlement of multi-edge skirted footings on the sand was generated by performing experiments in the laboratory. All tests were performed in s tank (700 mm length x 450 mm width x 600 mm depth), which was prepared with a 15 mm thick perspex sheet, stiffened by the mild steel plate as shown in Figure 1.

Source: Authors

Figure 1 Complete setup of the testing program. 

The 10 mm thick steel footings with plan dimensions 80 mm x 80 mm were prepared in the shape of a square, T, plus, and double box, as shown in Figure 2. The skirts with 5 mm thickness and depths between 0,25 B and 1,5 B were firmly welded to these footings to simulate skirted footings, as shown in Figure 3.

Source: Authors

Figure 2 Plan and a sectional view of various shapes of unskirted and skirted footings: (a) square, (b) Plus, (c) T shape, (d) Double box. 

Source: Authors

Figure 3 Photograph of model footings. 

In the case of rough footings, the sand was pasted to the base of the footing and also to the inner surfaces of the skirt; while, as for partly rough footings, a machined surface was used as it was. Beas river sand was used in this research and had a specific gravity (2,67), coefficient of uniformity (C u = 1,46), coefficient of curvature (C c = 0,98), and minimum and maximum dry unit weights of 13,06 kN/m3 and 15,97 kN/m3, respectively. As per IS 1498, the sand was classified as poorly graded (SP). The consolidated drained triaxial friction angle of sand at relative densities of 30%, 40%, 50%, and 60% was measured as 36,06°, 38,64°, 39,86°, and 41,72°, respectively. It is pertinent to note that due to the dilatancy effect, the friction angle of sand depends on the prevalent stress level or confining pressure in the test. Hence, it is likely that the friction angle developed in the small scale laboratory test, such as this one, will be much greater than the mobilized friction angle for the field size footing at failure with similar loading and soil conditions. Therefore, the predicted bearing capacity for the field size footing will be generally higher if extrapolated linearly from the results of the small-scale laboratory test data. This aspect is referred to as the size or scale effect on the bearing capacity (Chakraborty and Kumar 2013; Tang et al. 2014). To study this effect, it is necessary to carry out the load test with different footing sizes, which is a cumbersome task, especially with respect to multi-edge skirted footings. Hence, in the present study, no attempt has been made to study it. This implies that care is paramount when extrapolating the results of the present study for the field size footings. The tank was filled with sand to achieve the targeted relative density, which was varied from 30% to 60% in this investigation. The sand bed was placed in the tank in 8 equal 60 mm thick layers. The weight of the sand in each layer was calculated corresponding to the required relative density by using the unit weight of the sand and the volume of the layer. The weighed sand was then poured from a constant height and was compacted using a wooden 6 N rammer by giving several blows (obtained using the trial and error method) to reach the required relative density. It was ensured that the difference in measured relative densities was within ±1%. This was achieved by placing four steel bowls of known volume in each of the layers and measuring the achieved relative density of sand in them. As shown in Figure 1, the test on the prepared sand bed was performed with a strain-controlled loading frame of 50 kN and the employment of a load cell of 5 kN capacity. Note that, rather than measuring the pressure below the footing, the load applied on the top of the footing was recorded with a data acquisition system. This implies that this measured load, divided by the plan area of the footing, represents the magnitude of the average uniform pressure below the footing. Furthermore, any non-linear variation of normal stress was not considered (Kumar, 2009). This assumption is in linewith the literature, wherein the magnitude of the bearing capacity is defined based on the average uniform pressure. All the tests were performed using a 0,24 mm/min strain rate. It was intended that at this low strain rate, it would be possible to capture the non-linear pressure-settlement behavior. It was assumed that the measured penetration depth of the applied load would be a reasonable representation of the immediate settlement of the footing. However, further study is required to confirm this, since it is beyond the scope of the present work. After the completion of each test, the pressure-settlement curve was plotted, and, subsequently, the bearing capacity and the settlement at a given pressure were calculated. It was anticipated that the applied load on the skirted footing would be resisted by the shear resistance mobilized along the skirt-soil interface, along with the shear strength of the sand. The bearing capacity of the skirted footing was denoted in terms of the BCR. Further, the reduction in settlement of the footing, due to the provision of the skirts, at a given pressure, was expressed in terms of settlement reduction factor. For the development of the models (ANN and MRA), for the prediction of the BCR and the SRF, the data from the literature on square and circular footings was collected. After the development of the models for regular-shaped skirted footings, the models were extended for the prediction of the BCR and the SRF of T, Plus, and double box shapes by means of simple multiplication factors. Finally, the ANN and MRA model predictions were compared to bring out the superiority of the ANN over the MRA model.

Test Results

Pressure settlement curve and Bearing capacity variation of footing with skirt depth

The pressure-settlement behavior of the footings (skirted and unskirted) with different plan shapes corresponding to a relative density of 30% is shown in Figure 4. It is pertinent to mention here that the ultimate bearing capacity was taken with respect to the minimum of the peak pressure or the pressure corresponding to an s/B ratio of 10% on the pressure settlement curve. The settlement ratio was used for the calculation of the settlement of the footing. Furthermore, if a clear peak in the pressure-settlement ratio curve was not obtained, the ultimate bearing capacity was calculated by using a double tangent method.

Source: Authors

Figure 4 Pressure settlement behavior for footings with partly rough (a, c, e, g) and rough (b, d, f, h) interfaces for square (a, b), plus (c, d), double box (e, f), T (g, h) corresponding to a relative density of 30%. 

The results obtained in the current study related to the bearing capacity of the unskirted square footings (partly rough and completely rough) and were compared with the classical bearing capacity formula reported by Terzaghi (1943). To calculate the bearing capacity of the footing, the friction angle obtained from the triaxial test for different relative densities was used. This comparison is shown in Table 2.

Table 2 Comparison of bearing capacity values of the square footing with literature 

R d (%) Present work Terzaghi (1943)
Partly rough Rough
30 65,50 73,30 40,68
40 120,39 153,24 45,09
50 158,53 207,81 53,73
60 228,31 268,93 91,09

Source: Authors

From this table, it can be seen that the observed bearing capacity in the present case was higher in comparison to the one obtained by using Terzaghi's formula in all circumstances. This is perhaps due to a higher mobilized friction in the test on account of the dilatancy of the sand, which is more significant at a low-stress level; and due to a slight localized densification of the sand nearby the footing given the applied load. The study of Figure 4 shows that, irrespective of the plan shape of the footing, the bearing capacity increases with the increase in skirt depth. The obtained bearing capacity for unskirted footings is shown in Table 3.

Table 3 Variation of bearing capacity with relative density for the different plan shape of footing 

𝑹𝒅 (%) Interface condition Ultimate bearing capacity (kPa)
Square Plus Double box T
30 Partly rough 65,5 74,29 65,5 77,14
Rough 73,3 91,43 73,3 100
40 Partly rough 120,39 128,57 120,39 131,43
Rough 153,24 160 153,24 165,71
50 Partly rough 158,53 165,71 158,53 168,57
Rough 207,81 210,57 207,81 211,43
60 Partly rough 228,31 231,43 228,31 234,29
Rough 268,93 270 268,93 271,42

Source: Authors

A careful study of this table indicates that, for a given relative density, irrespective of the interface condition of the footing, T-shaped footings is provide the highest bearing capacity, followed by Plus, double box, and square shapes. The difference in the bearing capacity of the footing is quite substantial for the lower relative density, i.e., at Rd = 30%. It should be noted that, without any skirts, double box and square-shaped footings are both the same, hence the identical value of the bearing capacity.

The variation of the bearing capacity of the skirted footing (expressed in the form of BCR) with the skirt depth for the different relative density considering various plan shapes is shown in Figure 5. In all cases, the BCR for the square and the multi-edge skirted footings increases with the increase in skirt depth. The BCR at a relative density of 30% for a partly rough footing with D s /B = 0,25 was observed to be 1,62, 1,67, 1,91, and 2,03 for square, T, plus shape, and double box shapes, respectively. When the D s /B increased to 1,5, the BCR increased to 3,51, 3,86, 4,43, and 4,64, respectively. For a given skirt depth, the BCR decreased with the increase in relative density.

Source: Authors

Figure 5 BCR variation for partly (a, c, e, g) and completely (b, d, f, h) rough footing with skirt depth for plan square (a, b), plus (c, d), double box (e, f), and T (g, h) for relative density of 30%. 

Variation of SRF with skirt depth

The provision of the skirt along the periphery of the footing increased the bearing capacity and decreased the settlement below the footing. The settlement reduction in skirted footings was defined in a quantitative manner by using a settlement reduction factor:

where, s and s sk are the settlement of the unskirted and skirted footings, respectively, corresponding to a given pressure σ.

In the present study, the SRF was calculated at pressures of 25 kPa, 50 kPa, 100 kPa, and 200 kPa. Its variation with the skirt depth, for the different relative densities, at a pressure of 100 kPa, and for the various plan shape of the footing, is presented in Figure 6. As anticipated, for a given relative density, the SRF increased with the increase in skirt depth. In contrast, for the constant skirt depth value, the SRF decreased with the increase in relative density for all the plan-shaped footings. The considerable reduction in the settlement was obtained for the footing with a rough interface and at D s /B = 1,5 and R d = 30% as the SRF was observed to be 0,87, 0,83, 0,89, and 0,86 for square, Plus, double box, and T plan-shaped footings, respectively. This implies that, at a relative density of 30%, the reduction in the settlement due to the provision of the skirts could range between 80% and 90%, which is quite significant. Thus, skirted footings can be considered as an alternative ground improvement technique in loose sand wherein large settlements are often anticipated.

Source: Authors

Figure 6 SRF Variation for partly (a, c, e, g) and completely (b, d, f, h) rough footing with skirt depth with different relative densities of sand for cr = 100 kPa with a plan square (a, b), plus (c, d), double box (e, f), and T (g, h) shapes respectively. 

Artificial Neural Network

This work aims to model the neural network architecture to predict the bearing capacity ratio and settlement reduction factor for multi-edge footings on sand. However, considering the data availability on the bearing capacity and the settlement of the square and circular skirted footings in literature, an ANN model was first developed, which was applicable for these shapes. This model was later modified to suit the aim of this research. Generally, the bearing capacity of non-cohesive soil is dependent on the friction angle (Meyerhof, 1963; Vesic, 1973). In this study, along with the friction angle, skirt depth to width of the footing ratio was also used, since it is considered to affect the bearing capacity. Skirts play a significant role in improving the bearing capacity, as reported by some studies (Khatri, Debbarma, Dutta, and Mohanty, 2017; Khatri and Kumar, 2019). Regarding settlement prediction, normal stress was also plays a role, along with the friction angle and skirt depth (Al-Aghbari and Dutta, 2008; Al-Aghbari and Mohamedzein, 2018; Gnananandarao, Dutta, and Khatri, 2020). Hence, to predict BCR, the nondimensional skirt depth (D s /B) and friction angle (ф) were considered as input variables. In contrast, the prediction of SRF requires an additional input of pressure, along with the variables above.

The difficult job in the ANN model development is to determine the number of hidden layers and their neurons. The accuracy of the network model is dependent on the initially assigned weights and other associated parameters, as well as its architecture. However, until now, there is no defined technique to achieve the optimal architecture and parameter settings for an ANN model. Therefore, researchers have to follow the time-consuming trial and error approach. However, ANN have a disadvantage, such as being stuck in local minima and slow learning rates (Marto et al. 2014). Still, this technique is popular in Geotechnical applications.

Network Structure Preparation and Data Set Used

The performance of an artificial neural network entirely depends on its structure. The first step in its creation is to fix the number of hidden layers and hidden layer nodes. Generally, a rule of thumb has been adopted for choosing the hidden layer (one in this case) and the neurons in the hidden layer, which is 2/3 of the size of the input variables (Shahin, Maier, and Jaksa, 2002; Dutta etal. 2015b; Rezaeietal. 2016; Dutta, Rani, and Gnananandarao, 2018; Gnananandarao, Dutta, and Khatri, 2018). Considering this, the number of hidden layers and the hidden layer nodes were 1 and 2, respectively. The chosen ANN architecture for the prediction of the BCR and the SRF is shown in Figure 7.

Source: Authors

Figure 7 Artificial neural network diagram for (a) BCR and (b) SRF. 

The next difficult task is to fix the number of epochs. An excess number leads to overfitting of the data, whereas a low number leads to poor prediction. The statistical parameter (mean square error) was calculated between the actual and the predicted value corresponding to different epochs. The lowest mean squared error corresponding to an epoch is chosen to develop the neural network model. Byfollowing the procedure above, 100 and 450 were selected as the number of epochs for the BCR and the SRF, respectively. Finally, in this study, 2-2-1 and 3-2-1 structures for developing the ANN model were chosen for the BCR and the SRF, respectively. A total of 43 and 131 model experimental data points were collected from the literature (Prasanth and Kumar, 2017; Momeni etal. 2015a; Al-Aghbari, 2002; Al-Aghbari, 2007; Eid, Alansari, Odeh, Nasr, and Sadek, 2009) to predict these parameters in regular shaped skirted footings (square and circular).

Furthermore, an additional 94 and 266 data points were collected from the author's experimental work reported elsewhere (Gnananandarao, Khatri, and Dutta, 2018). Hence the total data points used for modeling were 137 and 397 for the BCR and the SRF, respectively. The data for the other regular shaped footings (rectangular and strip) were excluded due to scarcity in the literature. The range of the various input parameters for which the developed model is applicable is also presented in Table 4.

Table 4 Details of parameters used in ANN and regression models for BCR and SRF 

Parameters Minimum value Maximum value
The angle of internal friction (σ) (Deg.) 32,50 42,9
Skirt depth to width of the footing ratios (D s IB) 0 2
Bearing capacity of the skirted footing to bearing capacity of unskirted footing ratio (BCR) 1 7,8
Normal stress (σ) (kPa) 25 200
Settlement of skirted footing to settlement of unskirted footing ratio (SRF) 0,06 1

Source: Authors

To check the generalization capability of the model, about 70% and 30% of the data from the total data were respectively selected randomly for training and testing purposes. After its development, the model for these regular-shaped skirted footings was modified to account for the shape of multi-edge skirted footings with the introduction of multiplication factors. These factors were obtained by dividing the BCR/SRF of multi-edge footings with those of the predicted BCR/SRF in regular-shaped footings from the model. The multiplication factor for the BCR and the SRF was termed as F mbcr and F msrf, respectively. A total of 120 and 450 data on the BCR and the SRF of T, Plus, and double-box-shaped skirted footings were taken from the present experimentation. These data points were used to calculate the average multiplication factors to predict the BCR and the SRF of multi-edge skirted footings.

Activation Function Selection, Performance Measure, and Sensitivity Analysis

Artificial neural networks make use of various activation functions to create a relationship between the input and the output variables at each neuron layer. These functions are generally mathematical expressions and are used to produce the outputs. The different activation functions used were linear, threshold, threshold symmetric, sigmoid, sigmoid stepwise, sigmoid symmetric, sigmoid symmetric stepwise, gaussian, gaussian symmetric, gaussian stepwise, elliot, elliot symmetric, linear piece, linear piece symmetric, sin symmetric, cos symmetric, sin, and cos. All of these are available in the open-source Agiel neural network software. To choose the best activation function among these 18 functions, a comparison was made with various performance measures such as the coefficient of determination (R2), variance account for (VAF), mean square error (MSE), root mean square error (RMSE), mean absolute error (MAE), and mean absolute percentage error (MAPE).

The predictive models with an R2 of 1, VAF of 100%, and MSE, RMSE, and MAE of zero indicate a perfect model (Hajihassani, Abdullah, Asteris, and Armaghani, 2019; Huang, Asteris, Koopialipoor, Armaghani, and Tahir, 2019; Xu, Zhou, Asteris, Armaghani, and Tahir, 2019). The mathematical formulas for these performance measures (reported by Dutta et al. 2015b; Dutta et al. 2018; Gnananandarao et al. 2019) are provided in Table 5 for the BCR. In the same table, the performance measures for the SRF can be written just by replacing the BCR with the SRF. Based on the best statistical results, sigmoid symmetric and sigmoid were selected as the activation functions for the BCR and the SRF, respectively. The performance measures for these functions in the training and testing phase are shown in Table 6.

Table 5 Mathematical expressions for the performance measures for the BCR 

Source: Authors

Table 6 Performance measures for the training and testing of data for BCR and SRF for the best activation function 

Parameter Training Testing
BCR SRF BCR SRF
Activation function Sigmoid symmetric Sigmoid Sigmoid Symmetric Sigmoid
R2 0,96 0,91 0,94 0,93
VAF (%) 91,21 87,62 92,84 89,28
MSE 0,30 0,01 0,36 0,01
RMSE 0,55 0,09 0,60 0,09
MAE 0,36 0,07 0,42 0,06
MAPE 15,97 19,25 17,61 18,90

Source: Authors

The next important step is to choose a learning rate that analyses the performance of the activation function both for the training and the testing datasets after fixing the optimal epochs. The learning rate is vital for understanding any neural network structure, as it explains the influence of the errors on the hidden weights and biases. If the learning rate is low, the values will take a longer time to converge; otherwise, the model may be overfit to the target value. Thus, a default value of 0,7 was considered in this study. After getting the predicted BCR and the SRF of the regular shaped footings, using the procedure mentioned above, the next step was to check the accuracy of the predicted parameters regarding the target values. It was achieved after minimizing the errors. The comparison of predicted and the target BCR and the SRF values during training and testing step is shown in Figure 8a and b, respectively. This figure indicates that the coefficient of determination ranges between 0,91 and 0,96, which implies a good fit.

Furthermore, the sensitivity analysis was performed to know the direct or the indirect relation between the input and the output parameter saccording to the method used by Erzin and Gul (2014). The connection weights and biases obtained in the neural network were examined to understand the relative importance of the input and output parameters.

Source: Authors

Figure 8 Plot of predicted versus targeted (a) BCR (b) SRF values using ANN for training and testing. 

The result of the sensitivity analysis is presented in Figure 9, which suggests that for the prediction of the BCR as well as SRF, the input variable D s /B affects the most that is about 64% and 57%, respectively.

Source: Authors

Figure 9 Sensitivity analysis for (a) BCR and (b) SRF. 

Comparison with multivariable regression analysis

In the present study, apart from the development of ANN models, additional multivariable regression analyses were carried out to develop expressions for the predictions of the BCR and the SRF. The forecasts from these expressions were compared with ANN predictions. On similar lines with the ANN models, the MRA was first performed on the data related to regular shaped footings, and later, the developed expressions were modified for prediction of the BCR and the SRF of multi-edge footings with the introduction of multiplication factors such as Fmbcr and Fmsrf. To decide the form of the expression, the BCR of the skirted footing was compared with the ratio of the bearing capacity of the embedded footing to the surface footing:

therefore,

where, s q - s γ - shape factors, and d q - d γ - depth factors, according to Meyerhof's bearing capacity theory (1951,1963).

A careful study of Equation (3) suggests that the form of the expression for the regression analysis can be chosen as

In the expression above, P and Q are constants that will be obtained through regression analysis, since the N q /N γ ratio is dependent on the friction angle, which is replaced with the function of 𝜙 i.e. f(𝜙). The depth factor d γ is determined as follows:

In the present study, the non-linear regression analysis was carried out by using the Datafit software, version 9,1 (trial version). The final equation after regression analysis for the prediction of the BCR is:

Finally, the BCR of the multi-edge skirted footing can be related to the BCR given by Equation (6) as:

where F mbcr can be calculated by a procedure similar to the one described in the previous sections. The SRF equation for the regular shaped footing was obtained after the regression analysis as:

Furthermore, the SRF of multi-edge footings is related regular shaped footings as:

The R2 for Equations (6) and (8) is about 0,87 and 0,84, respectively, which is acceptable. A comparison between the BCR and the SRF predictions for the square and the multi-edge footings with partly rough and completely rough interfaces obtained from the regression analysis and the ANN is shown in Figures 10 and 11, respectively. These figures indicate that the predicted BCR and SRF are almost within the ±20% of the line of equality. The R2 ranges between 0,940 and 0,977 for the ANN model and 0,827 and 0,934 for the regression analysis. Furthermore, the study of Figure 11 suggests that the R2 for the SRF prediction may lie between 0,913 and 0,985 for the ANN model and 0,739 and 0,932 for the regression analysis. A comparison of Fmbcr and Fmsrf obtained from the ANN and the regression analysis for the parameters in question is provided in Table 7. It is quite noteworthy that these factors are comparable, even though they were obtained with different methodologies.

Source: Authors

Figure 10 Comparison of ANN with MRA models after the prediction of BCR for partly rough (a, c, e, g) and rough (b, d, f, h) interfaces of square (a, b), plus (c, d), double box (e, f), T (g, h) plan shapes. 

Source: Authors

Figure 11 Comparison of ANN model with MRA after the prediction of SRF for partly rough (a, c, e, g) and rough (b, d, f, h) interfaces of square (a, b), plus (c, d), double box (e, f), T (g, h) footings. 

Table 7 Multiplication factors Fmbcr and Fmsrf for prediction of BCR and SRF for multi-edge footings in ANN and regression models 

Plan shape of footing Multiplication factor (F mbcr and F msrf )
Artificial neural networks Regression Analysis
BCR SRF BCR SRF
Plus 1,03 0,98 1,02 0,99
Double box 1,11 0,89 1,09 0,92
T 1,04 0,97 1,03 0,98

Source: Authors

Equation for the BCR and the SRF from ANN

The goal of the present study was to develop an equation for the reliable prediction of the BCR and the SRF of the regular shaped and multi-edge skirted footings. From the comparison of the predictions obtained with the ANN and the regression analyses, it is quite clear that ANN predictions were always superior. Therefore, the equations for predicting the BCR and the SRF are also presented here.

The ANN model was developed with the optimum number of epochs (100) using the open-source Agiel software to obtain the weights and the biases introduced between the hidden layer and the output layer. The generalized function of the ANN model for the output (BCR) of the regular shaped skirted footing is defined by Equation 10:

The BCRSquare/circle resulting from Equation 14 is between - 1 and 1 for the activation function (sigmoid symmetric). Hence, there is a need for the denormalization of the output to get the actual value. The denormalized equation is as follows:

where BCRSquare/Circle(max) and BCRSquare.Circle(min) are the maximum and the minimum predicted bearing capacity ratios of regular shaped skirted footings, respectively.

The BCR of multi-edge footings can be obtained by following Equation (7), which was defined earlier. On similar lines with the BCR, the SRF can be predicted by following Equations 16-19 as given below:

The SRFSquare/Circle predicted by Equation (19) ranged between -1 and 1 for the activation function (sigmoid symmetric). Hence, there is a need for the denormalization of output (SRF) to get the actual value. The denormalized equation (20) is shown below:

Where SRFSquare/Circle(max) and SRF Square/Circle(min) are the maximum and the minimum predicted settlement reduction factors of regular shaped skirted footings. Furthermore, the SRF of multi-edge footings is obtained following Equation (9) and Table 7. The weights and biases between the hidden layer neurons with input and the output nodes for the BCRSquare/Circle and the SRFSquare/Circle prediction are provided in Tables 8 and 9, respectively.

Table 8 Weights and biases between hidden layer neurons with input and output nodes for BCR prediction 

Source: Authors

Table 9 Weights and biases between hidden layer neurons with input and output nodes for SRF prediction 

Source: Authors

Conclusions

The study carried out in this work showed the feasibility of using a simple ANN and multivariable regression analysis to predict the bearing capacity ratio and settlement reduction factor of multi-edge skirted footings resting on sand. Based on the obtained results, the following conclusions can be drawn:

  1. For a given relative density, the non-dimensional bearing capacity of square/circular and multi-edge footings was defined in terms of BCR increases with an increase in depth of the skirt. The BCR of multi-edge footings was marginally higher than that of its square/circular counterpart, which is desirable.

  2. For a given constant pressure, the reduction in settlement of skirted footings, expressed in the form of SRF, increased with the increase in skirt depth. However, contrary to the observation of BCR, the SRF of multi-edge footings was marginally smaller in comparison to regular-shaped skirted footings.

  3. The BCR and SRF predicted by using ANN or regression analyses are in agreement with experimental values, although ANN provides superior predictions in all cases.

  4. The developed equations can be used to predict the BCR and SRF of partly rough and completely rough skirted footings as the predictions are not much affected by interface roughness.

It is anticipated that the outcome of this study will help in emphasizing the use of multi-edge skirted footings wherever possible. Furthermore, the developed equations can be used to predict the bearing capacity and reduction in settlement of such footings at a given pressure without conducting the field tests. Neural network models, in general, cannot provide reasons and reasoning beyond the model (black boxes) thus obtained; once a network knows one set of weights, any new learning results in disastrous forgetting. The suitability of alternative techniques such as vector supporting machines, particle swarm optimization, or genetic programming may also be explored in the future.

References

Al-Aghbari, M. Y., and Khan, A. J. (2002). Behaviour of shallow strip foundations with structural skirts resting on dense sand. Proceedings of challenges of concrete constructions, 6, 737-746. 10.1680/cfec.31784.0072 [ Links ]

Al-Aghbari, M. Y. (2007). Settlement of shallow circular foundations with structural skirts resting on sand. The Journal of Engineering Research, 4(1), 11-16. 10.24200/tjer.vol4iss1pp11-16 [ Links ]

Al-Aghbari, M. Y. and Mohamedzein, Y. E.-A. (2020). The use of skirts to improve the performance of a footing in sand. International Journal of Geotechnical Engineering, 14(2), 134-141. 10.1080/19386362.2018.1429702 [ Links ]

Armaghani, D. J., Faradonbeh R. S., Rezaei, H., Rashid A. S. A., and Amnieh, H. B. (2018). Settlement prediction of the rock-socketed piles through a new technique based on gene expression programming. Neural Computing and Applications, 29, 1115-1125. 10.1007/s00521-016-26188 [ Links ]

Armaghani, D. J., Shoib, R. S. N. S. B. R., Faizi, K., and Rashid, A. S. A. (2017). Developing a hybrid PSO-ANN model for estimating the ultimate bearing capacity of rock-socketed piles. Neural Computing and Applications, 28, 391-405. 10.1007/s00521-015-2072-z [ Links ]

Chakraborty, D. and Kumar, J. (2013). Dependency of Ny on footing diameter for circular footings. Soils and Foundations, 53(1), 173-180. 10.1016/j.sandf.2012.12.013 [ Links ]

Chen, W., Sarir, P., Bui, X-N., Nguyen, H., Tahir, M. M., and Armaghani D. J. (2019). Neuro -genetic, neuro -imperialism and genetic programing modelsin predicting ultimate bearing capacity of pile. Engineering with Computers, 36, 1101-1115. 10.1007/s00366-019-00752-x [ Links ]

Dawarci, B., Ornek, M., and Turedi, Y. (2014). Analysis of multi-edge footings rested on loose and dense sand. Periodica Polytechnica Civil Engineering, 58(4), 355-370. 10.3311/PPci.2101 [ Links ]

Dutta, R. K., Dutta K., Jeevanandham, S. (2015). Prediction of deviator stress of sand reinforced with waste plastic strips using neural network. International Journal of Geosynthetics and Ground Engineering. 1(2), 1-12. 10.1007/s40891-015-0013-7 [ Links ]

Dutta, R. K., Rani, R., and Gnananandarao, T. (2018). Prediction of ultimate bearing capacity of skirted footing resting on sand using artificial neural networks. Journal of Soft Computing in Civil Engineering, 2(4), 34-46. 10.22115/SCCE.2018.133742.1066 [ Links ]

Eid, H. T., Alansari, O. A., Odeh, A. M., Nasr, M. N., and Sadek, H. A. (2009). Comparative study on the behavior of square foundations resting on confined sand. Canadian Geotechnical Journal, 46, 438-453. 10.1139/T08-134 [ Links ]

Erzin, Y. and Gul, T. (2014). The use of neural networks for the prediction of the settlement of one-way footings on cohesionless soils based on standard penetration test. Neural Computing and Applications, 24, 891-900. 10.1007/s00521-012-1302-x [ Links ]

Ghazavi, M. and Mokhtari, S. (2008). Numerical investigation of load-settlement characteristics of multi-edge shallow foundations. In ed. Jadhav, M. N. (Ed.) Proceedings of The 12th International Conference of International Association for Computer Methods and Advances in Geomechanics (lACMAG), Red Hook, NY: Curran. pp. 3344-3351. [ Links ]

Gnananandarao T., Khatri V. N., and Dutta R. K. (2018). Performance of multi-edge skirted footings resting on sand. Indian Geotechnical Journal, 48(3), 510-519. 10.1007/s40098-017-0270-6 [ Links ]

Gnananandarao, T., Dutta, R. K. and Khatri, V. N. (2019). Application of artificial neural network to predict the settlement of shallow foundations on cohesionless soils. Geotechnical Applications, 13, 51-58. 10.1007/978-98113-0368-5_6 [ Links ]

Gnananandarao, T., Dutta, R. K. and Khatri, V. N. (2020). Model studies of plus and double box shaped skirted footings resting on sand. International Journal of Geo-engineering, 11 (2), 1-17. 10.1186/s40703-020-00109-0 [ Links ]

Hajihassani, M., Abdullah, S. S., Asteris, P. G., and Armaghani, D. J. (2019). A Gene Expression Programming Model for Predicting Tunnel Convergence. Applied Sciences, 9, 4650. 10.3390/app9214650 [ Links ]

Harandizadeh, H., Armaghani, D. J., and Khari, M. (2019). A new development of ANFIS-GMDH optimized by PSO to predict pile bearing capacity based on experimental datasets. Engineering with Computers, 1-16. 10.1007/s00366-019-00849-3 [ Links ]

Huang, L., Asteris, P. G., Koopialipoor, M., Armaghani, D. J., and Tahir, M. M. (2019). Invasive Weed Optimization Technique-Based ANN to the Prediction of Rock Tensile Strength. Applied Sciences, 9, 5372. 10.3390/app9245372 [ Links ]

IS 1498 (1970). Classification and identification of soils for general engineering purposes. Delhi, India: Bureau of Indian Standards. [ Links ]

IS 6403 (1981). Determination of bearing capacity of shallow foundation. Delhi, India: Bureau of Indian Standards . [ Links ]

Kalinli, A., Acar, M. C., and Gunduz, Z. (2011).New approaches to determine the ultimate bearing capacity of shallow foundations based on artificial neural networks and ant colony optimization. Engineering Geology, 117, 29-38. 10.1016/j.enggeo.2010.10.002 [ Links ]

Khari, M., Armaghani, D. J., and Dehghanbanadaki, A. (2020). Prediction of Lateral Defection of Small-Scale Piles Using Hybrid PSO-ANN Model. Arabian Journal for Science and Engineering, 45, 3499-3509. 10.1007/s13369-019-04134-9 [ Links ]

Khatri, V. N. and Kumar, J. (2019). Finite-Element Limit Analysis of Strip and Circular Skirted Footings on Sand. International Journal of Geomechanics, 19(3), 06019001. 10.1061/(ASCE)GM.1943-5622.0001370 [ Links ]

Khatri, V. N., Debbarma, S. P., Dutta, R. K., and Mohanty, B. (2017). Pressure-settlement behavior of square and rectangular skirted footings resting on sand. Geomechanical Engineering 12(4),689-705. 10.12989/gae.2017.12.4.68 [ Links ]

Khudier, A. S. (2018). Prediction of bearing capacity for soils in basrah city using artificial neural network (ANN) and multilinear regression (MLR) models. International Journal of Civil Engineering and Technology, 9(4), 853-864. [ Links ]

Kumar, J. (2009). The variation of Ny with footing roughness using the method of characteristics. International Journal for Numerical and Analytical Methods in Geomechanics, 33(2), 275-284. 10.1002/nag.716 [ Links ]

Marto, A., Hajihassani, M., and Momeni, E. (2014). Bearing Capacity of Shallow Foundation's Prediction through Hybrid Artificial Neural Networks. Applied Mechanics and Materials, 567, 681-686. 10.4028/www.scientific.net/AMM.567.681 [ Links ]

Meyerhof, G. G. (1951). The ultimate bearing capacity of foundations. Geotechnique, 2(4) 301-332. 10.1680/geot.1951.2.4.301 [ Links ]

Meyerhof, G. G. (1963). Shallow foundations. Journal of Soil Mechanics and Foundation Division. ASCE, 91(SM2), 21-31. [ Links ]

Meyerhof, G. G. (1963). Some recent research on bearing capacity of foundations. Canadian Geotechnical Journal, 1, 16-26. 10.1139/t63-003 [ Links ]

Momeni, E., Nazir, R., Armaghani, D. J., and Maizir, H. (2015a). Application of Artifcial Neural Network for Predicting Shaft and Tip Resistances of Concrete Piles. Earth Sciences Research Journal, 19(1), 85-93. 10.15446/esrj.v19n1.38712 [ Links ]

Momeni, E., Armaghani, D. J., Nazir, R., and Sohaie, H. (2015b). Bearing capacity of precast thin-walled foundation in sand. Proceedings of the Institution of Civil Engineers Geotechnical Engineering, 168(GE6), 539-550. 10.1680/jgeen.14.00177 [ Links ]

Momeni, E., Armaghani, D. J., Fatemi, S. A., and Nazir, R. (2017). Prediction of bearing capacity of thin-walled foundation: a simulation approach. Engineering with Computers, 3(2), 319-327. 10.1007/s00366-017-0542-x [ Links ]

Nazir, R., Momeni, E., Marsono, K., and Sohaie, H. (2013). Precast spread foundation in industrialized building system. In Hossain, M. Z. and Hossain, M. S. (Eds.) Proceedings of the 3rd International Conference on Geotechnique, Construction Materials and Environment - GEO-MATE 2013 (pp. 47-52). Nagoya, Japan: Nagoya Institute of Technology. [ Links ]

Nazir, R., Momeni, E., Marsono, K., and Maizir, H. (2015a). An artificial neural network approach for prediction of bearing capacity of spread foundations in sand. Journal Teknologi, 72(3), 9-14.10.11113/jt.v72.4004 [ Links ]

Nazir, R., Momeni, E., and Marsono, K. (2015b). Prediction of bearing capacity for thin-wall spread foundations using ICA-ANN predictive model. In: Proceedings of the International Conference on Civil, Structural, and Transportation Engineering, Ottawa, Ontario-May, 4-5. [ Links ]

Prasanth, T. and Kumar, P.R. (2017). A study on load carrying capacity of skirted foundation on sand. International Journal of Science and Research, 6(6), 2231-2235. [ Links ]

Rezaei, H., Nazir, R., and Momeni, E. (2016). Bearing capacity of thin-walled shallow foundations: an experimental and artificial intelligence-based study. Journal of Zhejiang University Science A: Applied Physics and Engineering, 17(4), 273-285. 10.1631/jzus.A1500033 [ Links ]

Shahin, M. A., Maier, H. R., and Jaksa, M. B. (2002). Predicting settlement of shallow foundations using neural networks. Journal of Geotechnical and Geoenvironmental Engineering, ASCE, 128(9), 785-793.10.1061/(ASCE)1090-0241(2002)128:9(785) [ Links ]

Tang, C., Phoon, K.K., and Toh, K.C. (2014). Effect of footing width on Nγ and failure envelope of eccentrically and obli(6), 694-707. 10.1139/cgj-2013-0378 [ Links ]

Terzaghi, K. (1943). Theoretical soil mechanics. New York: John Wiley and Sons. 10.1002/9780470172766 [ Links ]

Vesic, A.S. (1973) Analysis of Ultimate Loads of Shallow Foundations. Journal of the Soil Mechanics and Foundations Division, 99, 45-73. https://cedb.asce.org/CEDBsearch/record.jsp?dockey=0020165Links ]

Xu, H., Zhou, J., Asteris, P. G., Armaghani, D. J., and Tahir, M. Md. (2019). Supervised Machine Learning Techniques to the Prediction of Tunnel Boring Machine Penetration Rate. Applied sciences, 9, 1-19. 10.3390/app9183715 [ Links ]

Yong, W., Zhou, J., Armaghani, D. J., Tahir, M. M., Tarinejad, R., Pham, B. T., and Huynh, V. V. (2020). A new hybrid simulated annealing-based genetic programming technique to predict the ultimate bearing capacity of piles. Engineering with Computers, 1-17. 10.1007/s00366-019-00932-9 [ Links ]

Ziaee, S. A., Sadrossadat, E., Alavi, A. H., and Shadmehri, D. M. (2015). Explicit formulation of bearing capacity of shallow foundations on rock masses using artificial neural networks: application and supplementary studies. Environmental Earth Science, 73(7), 3417-3431. 10.1007/s12665-014-3630-x [ Links ]

How to cite: Gnananandarao, T., Khatri, V. N., and Dutta, R. K. (2020). Bearing capacity and settlement prediction of multi-edge skirted footings resting on sand. Ingeniería e Investigación, 40(3), 9-21. 10.15446/ing.investig.v40n3.83170

Received: October 26, 2019; Accepted: July 30, 2020

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License