Evaluation of Induced Settlements of Piled Rafts in the Coupled Static-Dynamic Loads Using Neural Networks and Evolutionary Polynomial Regression

Coupled Piled Raft Foundations (CPRFs) are broadly applied to share heavy loads of superstructures between piles and rafts and reduce total and differential settlements. Settlements induced by static/coupled static-dynamic loads are one of the main concerns of engineers in designingCPRFs. Evaluation of induced settlements of CPRFs has been commonly carried out using threedimensional finite element/finite difference modeling or through expensive real-scale/prototype model tests. Since the analyses, especially in the case of coupled static-dynamic loads, are not simply conducted, this paper presents two practical methods to gain the values of settlement. First, different nonlinear finite differencemodels under different static and coupled static-dynamic loads are developed to calculate exerted settlements. Analyses are performed with respect to different axial loads and pile’s configurations, numbers, lengths, diameters, and spacing for both loading cases. Based on the results of well-validated three-dimensional finite differencemodeling, artificial neural networks and evolutionary polynomial regressions are then applied and introduced as capable methods to accurately present both static and coupled static-dynamic settlements. Also, using a sensitivity analysis based on Cosine Amplitude Method, axial load is introduced as the most influential parameter, while the ratio l/d is reported as the least effective parameter on the settlements of CPRFs.


Introduction
CPRF is representative for combined piled raft foundation, which is commonly applied to suffer the heavy load of skyscrapers through sharing the exerted load between raft and piles.Coefficient of piled raft,  pr , controls load sharing ratio and is defined as the ratio of the sum of the loads carried by piles to the corresponding value of the resistance of the whole system.It varies from 0 for spread footings to 1 for pile foundations.Also, as shown in Figure 1, the settlements of a CPRF ( pr ) can be calculated based on the value of ( pr ) and the settlement of pile foundation ( pf ).
Piles are generally used to decrease the foundation's total settlements.In the last three decades, some researchers like those of [1][2][3] have declared efficient load sharing mechanisms to accurately investigate the piled raft behavior.
New well-developed computational approaches and available finite element (FE) codes such as Defpig and Napra have facilitated modeling of available interactions [4][5][6].Besides, some other small scale model studies [7,8] and numerical and analytical simulations in the literature have made improvements in the design process [9][10][11].In addition, results of monitoring the behavior of piled rafts supporting heavy loads of some skyscrapers are reported [12,13].But there is an obvious need for a comprehensive study about the piled raft settlements in the case of combined vertical static loads and horizontal dynamic ones.Since it is experimentally time-consuming and expensive, numerical modeling is considered as a good alternative to accurately study CPRFs behavior.However, in order to have a comprehensive and accurate model, especially in the case of CPRFs for high rise buildings with large number of piles, there will be a large number of structural elements and the calculation time makes the analysis so difficult.To cope with this difficulty, recently, some research works are focused on the application of soft computing techniques, as the speedy and powerful tool and alternative for other timeconsuming and expensive methods, to the problem of piles and piled raft foundations behaviors [14,15].In the research of Armaghani et al. (2015) [14], the ultimate bearing capacity is investigated, while Baziar et al. (2014) [15] evaluated the static induced settlements.Hence, the dynamically loaded induced settlements of piled rafts and the application of soft computing techniques to this problem have not been completely studied.Then, there is an obvious need for a speedy and accurate tool for predicting piled raft settlements in the combined vertical static-horizontal dynamic loading conditions.Relying on the fact that conducting real-scale or prototype scale field or experimental tests on the dynamic behavior of Coupled Piled Raft Foundations needs a lot of time and money and also considering the process and analysis time of numerical modeling, if an alternative method for predicting the settlements of CPRFs in dynamic loading cases exists, it will be useful.This method should have advantages of both accurate experimental modeling and low calculation time.Hence, the application of neural network modeling and evolutionary polynomial regressions (EPRs), which are believed to be common ways to accurately and timely predict engineering complicated functions, can be examined.Different attempts to apply neural networks and EPRs to model different civil and geotechnical problems are presented in the literature [14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32][33].They are well-applied in a wide range of problems from deep soil stabilizations, concrete, and their related structures, compressive strength of soils, rocks, and stabilized samples, bearing capacity of shallow and deep foundations, lateral spreading, rock mechanics, rock engineering, and soil mechanics [14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32][33].In these methods, a sufficient amount of accurate data is required initially to train artificial neural networks and evolutionary polynomial regression models.In this regard, firstly, based on nonlinear dynamic finite difference method, piled raft foundation models were developed and successfully validated with available field tests.Based on the parametric study on some geometrical parameters of CPRF and loads, a comprehensive sight to the problem has been achieved and a strong database has been prepared and results of induced settlements for both static and coupled static-dynamic conditions are calculated and gathered regarding input concerning parameters.Then, using new powerful forecasting methods, for example, neural network modeling and evolutionary polynomial regression modeling, it has been attempted to present comprehensive, speedy, and accurate models for predicting CPRF's settlements in static and combined staticdynamic loading conditions.Hence, as the novelty of the paper, new neural network and EPR models are presented for modeling settlements of CPRF in the static and coupled static-dynamic loading conditions.As another main result of the paper, the most and the least influential parameters, which, respectively, should be paid more and less attentions in dynamic CPRF modeling, are introduced.

Nonlinear Finite Difference Modeling of CPRF
In order to develop CPRF models, FLAC 3D software was used and CPRFs with regular 3 × 3, 4 × 4, and 5 × 5 arrangements were modeled in both static and combined static-dynamic conditions.
In order to achieve a comprehensive model, three different arrangements and topologies were considered for piles (pattern N for near piles, pattern M for mediate piles, and pattern F for far piles).Also, 4, 5, and 6 were the ratios of / (spacing of piles to the piles' diameters).Also, assuming four (/) conditions (16, 20, 24, and 32) and two different static vertical pressure states (60 and 90 kPa) static and combined static-dynamic settlements were calculated.
In order to model three-dimensional soil and pile geometry, 8-node brick elements and 6-node cylinder elements are, respectively, used.Names, shapes and other specifications of the used mesh shapes are presented in Figure 2.
Also, there are some well-developed constitutive models presented in FLAC as built in material behaviors.Among them, Mohr-Coulomb's model applies more simple constitutive parameters (in view of characterization), internal friction angle, cohesion, and dilation.Hence, it is considerably applied to model the behavior of surrounding soil.Also, in order to investigate the interaction between the soil and piles, shear and normal coupling springs are used.In this paper, using nonassociative Mohr-Coulomb model (the yield surface is presented in Figure 3), the soil nonlinear behavior is modeled, while slippage and separation during the motion are described using coupling springs.The normal behavior of the pile/grid interface is represented by a spring with a limiting normal force, which is dependent on the direction of piles' nodal movement.In addition, the shear behavior of the pile/grid interface is represented as a spring-slider system at the nodes of modeled piles.The shear behavior of the interface during relative displacement between the nodal points of pile and the grid is numerically attained by adjusting the proper shear and normal stiff-nesses.Elements of the interface and parameters of the used constitutive model in the interface along with the interface elements along the bottom and the side are presented in Figure 4. : Specifications of the used elements [16].For the analysis of the piled raft foundation, the soil as well as the structure is discretized into elements.Usual practice is to divide the soil mass into rectangular zones of aspect ratios less than 4 : 1.For the region of greater interest, finer discretization is used.The software internally breaks each of the rectangular elements into four overlapping triangular elements.The limits of the discretized grid should be properly planned to suit the geometry of the problem.To avoid boundary effects, the mesh should extend sufficiently beyond the region of interest.Experiences of the previous researchers are taken as guidelines to decide the extent of the grid.The software also allows graded discretization, which is more efficient than abrupt change in zone sizes.In this regard, a sensitivity analysis on the number of grids (mesh refinement) is conducted to obtain the proper mesh size, in which the size of elements does not meaningfully change the accuracy of results.Figure 5 presents the adopted mesh for piles cap and the soil and the pile's length in the optimum refined mesh view.A typical shape for one of the developed models has been shown in Figure 6, where elements near the piles are finer than others to accurately track the sensitivities in this zone (with higher stress intensity).In addition, in order to neglect the boundary effects, the lateral boundaries are in the distances equal to 50 from the edge of raft.Moreover, the distance between lower boundaries with the tip of piles is assumed to be 3 (length of piles).It should be noted that interactions between soil and piles have been accurately considered using interface elements and [34] where  is the soil's bulk modulus,  is the shear modulus of the soil, Δ min is the minimum width of neighbor elements at the interface, and   and   , respectively, present shear and normal stiffness of elements at interface nodes.
2.1.Dynamic Analysis.Critical time step is calculated using [34] where   is the speed of  wave,  represents the area of triangular element, and Δ max is the maximum dimension of the area.In the case of damping proportional to the stiffness, the critical time step is calculated using [36]  where  max is maximum predominant frequency of system and  represents a ratio of critical damping at this frequency and are calculated using [34]  In these equations,  min and  min are Rayleigh's damping and angular frequency.

Boundary Conditions. Quiet boundaries introduced by
Lysmer and Kuhlemeyer [34] are used to neglect the effects of wave reflections in the model.In these boundaries, normal and shear quiet tensions are modeled using [34] where   and   are, respectively, normal and shear components of wave velocity at boundary. is specific gravity and   and   represent  and  wave's speed, respectively.Also,   and   , respectively, show the normal and shear quiet tensions.Figure 7 shows quiet boundaries used in FLAC for rigid and flexible beds.

Loading Condition.
As mentioned, in order to evaluate the effect of coupled vertical static loading and horizontal dynamic loading, 60 and 90 kPa vertical loads were separately exerted to the structure and the model was analyzed.Also, dynamic loading was modeled using horizontal stress wave as shown in Figure 8.Since quiet boundaries were used, acceleration/speed time histories at boundaries cannot be exerted to the models.Hence, based on (6) [34], the wave velocities were changed to the waves' stress terms and shown in Figure 8.
where   and   are, respectively, normal and shear components of the wave velocity at boundary. is specific gravity and   and   represent  and  wave's velocity, respectively.

Validation of Numerical Model
In order to further study and investigate the role of each of the concerning parameters, the previously well-validated threedimensional finite difference modeling procedure (previously applied by the corresponding author of the paper) is used [6].This model was validated using a 1 g physical model test on medium dense sand with the soil and pile geometry shown in Figure 9.In addition, the adopted mesh along with the comparison of physical modeling results with the results of numerical modeling are also shown in Figure 9.As it can be seen, there is an acceptable agreement between the results of three-dimensional nonlinear finite difference modeling and 1 g physical modeling [6].

Parametric Study on Finite Difference Model
As it was mentioned, the effect of some of the geometrical pile parameters and the loads on the static and combined staticdynamic CPRF's settlements has been studied.Considering three pile patterns (different configurations), two different pile diameters, three ratios for /, four ratios for /, and two different static load states, a databank with 144 static data series and 144 dynamic ones has been achieved.Figures 10 and 11 show the results of static and combined staticdynamic settlements for the 3 × 3 CPRF as a sample of the results.As shown, F-architecture is the best architecture, in which the minimum settlements are observed in both static and coupled static-dynamic models.Besides, increasing / ratio will result in an increase in the calculated settlements.However, piles with larger lengths cause smaller induced static and coupled static-dynamic settlements.The combined effects of / and / can be better tracked using a sensitivity analysis.The following sections present the sensitivity analysis, which shows the relative importance of increasing / and decreasing / in the reduction of induced settlements.
As shown and comparing to the effect of increasing / ratio, decreasing s/d ratio plays more important role in the reduction of the settlements in both cases.Furthermore, Tables 1 and 2, respectively, present the results of parametric study on the whole affecting parameters for static loading condition and combined static-dynamic loading condition.
The following subsections present parametric studies conducted on different concerning parameters (e.g., piles' lengths, diameters, spacing, architecture, and axial loads).In this regard and as sample cases, settlements of CPRFs under coupled static-dynamic loads are presented and discussed for different parameters.

Effect of Piles' Architecture (Pattern) on the Settlements of CPRF.
To investigate the effect of piles' arrangement and their architecture on the static and coupled static-dynamic settlements, three different patterns are considered, where N, F, and M symbols stand for the near-, far-, and mediumdistance piles architectures, respectively.As described, 0.3 and 0.5 meters are the values of piles' diameter used to calculate the induced settlements.In this section, effects of different piles' architectures on the coupled static-dynamic settlements of 0.5 meters in diameter piles under 60 kPa static axial load are evaluated.Figure 12  should be noted that the reported value for the settlement of each CPRF is the average of pile tip's settlement since a rigid behavior is assumed for the piles' cap.
As shown, N-architecture CPRFs experience the largest settlements, while F-architecture shows a better resistance against settlements.This is because of the piles group's performance in the CPRF.Indeed, in far-piles architecture (F), with decreasing the relative effect of piles, each individual pile behaves as a single pile and enhances the resistance of the whole CPRF.Nevertheless, in the N-architecture, overlapping stress bubbles of individual piles, the maximum resistance of the whole CPRF against the settlement decreases.It is seen that, for a constant architecture, with increasing the length of piles, value of coupled static-dynamic settlements decreases.

Effect of Axial Pile Loads on the Settlements of CPRF.
In this section, as a sample case, settlements of a 3 * 3architecture CPRF (with constant diameter piles,  = 0.5 m) under 60 kPa and 90 kPa static axial force subjected to a horizontal dynamic load are evaluated.As shown in Figure 14, the value of axial force affects the induced settlements.It is observed that, with increasing the static axial force, the induced settlements increase.

Artificial Neural Networks.
Artificial Intelligence has different branches.Among these branches, Artificial Neural Networks (ANNs) are known as the black box, while Evolutionary Polynomial regression (EPR) is referred to as the grey box method.In ANNs, generally, experimental datasets are used to attain relationships between input and output parameters.In this method, the large number of databases is one of the requirements to get the results with less error [17][18][19][20][21][22]37].
The relationship between input and output parameters is gained using the learning rules.In spite of classical statistical methods, neural networks do not need any previous knowledge about the quality and mechanics of the problem and their concerning parameters.One of the most in-demand kinds of neural networks is multilayer perceptron (MLP) networks, which is formed using correct definition of its constructing layers, input, and hidden and output layers [23].Regarding the type of problem complexity and its nonlinearity, the number of MLP layers is defined [24].The constructing elements of each layer are neurons, which are connected to the neurons of other layers, but cannot make connection with the neurons of the same layer.The relationship between neurons of a layer to the neurons of the next neighboring layer is carried out using some connections.The analyzed information will be multiplied to the assigned weight of each node and then will be transformed using the activation function [38].In an MLP, as shown in Figure 15, independent input variables are connected to the neurons of hidden layers and can predict the multiple dependent output variables based on network training (by detecting the similarities between input and output parameters and minimizing the prediction errors) [25].Hence, architecture of the network, learning rule, and the transfer function are key parameters, which should be correctly defined to build a proper neural network [35].

Neural Network Training.
Feed-Forward Backpropagation, FFBP, is the most applied algorithm used to train the neural networks for a broad range of engineering applications [26].In this training algorithm, firstly and in forward pass, assuming a primary value for connection between neurons, outputs are forecasted and then the computational error is calculated.This error will be backcalculated to update the  neuron's weights and obtain accurate outputs.This second phase is called backpropagation phase [27].
In order to make all the input and output parameters dimensionless, (7) [28,29] is used and all the parameters are normalized to a 0-1 scale.
In this paper, training and cross-validation datasets were made through random selection of 85% of the whole data.The remaining 15% are the testing data series.Different transfer functions, for example, TANSIG, LOGSIG, and PURELIN, are available to be used as the transfer function to mathematically represent, in terms of spatial or temporal frequency, the relation between the input and output.The transfer functions usually have a sigmoid shape, but they may also take the form of other nonlinear functions, piecewise linear functions, or step functions.Hence, to evaluate the ability of different transfer functions in the prediction of coupled static-dynamic CPRF's settlements, the mentioned transfer functions were applied and TANSIG function showing the best performance was selected for the modeling purposes.Figure 16 shows the TANSIG transfer function.Also, (8) shows the formula used to calculate TANSIG function [39].
where   represents the weighted sum for each of the concerning input parameters [39].
In addition, underfitting and overfitting should be avoided during constructing the networks.Indeed, using incorrect number of solution epochs leads to these two phenomena.Overfitting corresponds to the case of too many training cycles, while underfitting declares that insufficient epochs are used to train the network [30].squares method to estimate the constants of a previously evolutionary developed model [40].The general steps applied in EPR are depicted in Figure 17.

Evolutionary Polynomial Regression
In this method, type of adopted functions, number of terms exponents' range, and the number of generations are the parameters which affect the output relationship [31].Besides, similar to the neural network modeling, the coefficient of determination evaluates the degree of accuracy of the proposed equation.If the proposed model satisfied both termination criteria (maximum number of generations and sentences) and required prediction accuracy, it will terminate.Otherwise, it goes through another evolution [32].

Results
In order to evaluate and compare the developed models and techniques, their performances have been investigated in different various terms.
where   and   are calculated and forecasted outputs, respectively.Also,  is the number of datasets.

Neural Network Performances.
To evaluate the performance of the optimum models, the normalized predicted values of static and combined static-dynamic settlements (for testing data series) are predicted using both conventional and genetic algorithm based neural networks.
(1) Performance of Neural Networks in the Static Loading Case.
The value of RMSE and coefficient of correlation along with the slope of fitting line, , for different network architectures are shown in Table 3, where different terms, digits, of each of the network architectures represent the number of neurons in its corresponding layers.Figure 18 presents correlation between predicted settlements versus calculated ones for the best network, 5-3-12-1.
Moreover, the predicted and calculated datasets are shown in Figure 19.developed models in the coupled static-dynamic loading cases are presented in Table 4.
Figure 20 presents correlation between predicted settlements versus calculated ones for the best network, 5-8-8-1.Moreover, the predicted and calculated datasets are shown in Figure 21.

EPR Models.
Using EPR modeling, different relations have been developed for both static and coupled staticdynamic settlements based on input parameters.Equations (10) show the best optimum relations for settlement prediction in static and coupled static-dynamic loading conditions, respectively.
As it is shown in Figures 22 and 23, the coefficients of determination for both static and static-dynamic coupled models presented using EPR are acceptable values of 97.49% and 97.26%, respectively, which mean that the proposed models are capable for predicting the settlements based on the studied input parameters.Moreover, the values of RMSE for static and coupled static-dynamic settlement models are 0.271 and 0.364, respectively.

Sensitivity Analysis.
The strength of relationship between affecting parameters (number of piles, radius, ratio of /, ratio of /, and axial force) and the values of settlements of CPRF in both static and coupled static-dynamic conditions can be found using sensitivity analysis based on Cosine Amplitude Method (CAM) [33] of CPRF in both static and coupled static-dynamic conditions.In this method, the strength of relationship between target parameter and inputs is found using [33,41]: where each of the input parameters is expressed as one of  array's elements shown in [33]  = {  where each of its elements is a vector with the length of  and is presented in [33]   = { 1 ,  2 ,  3 , . . .  } .
Figures 24 and 25 show the strength of relationship between input parameters and CPRF settlements in static and coupled static-dynamic case, respectively.As it can be seen, axial load is the most important parameter in both conditions.Also, / is the least effective parameter on the static and coupled static-dynamic CPRF's settlements.Besides, effects of other concerning parameters on the settlements are quiet the same.

Conclusion
Knowing the piled raft settlements in static and coupled static-dynamic loading conditions plays an important role in preliminary design phase of these combined foundations.
Regarding the high expenses of full-scale model tests and experimental tests in determining the earthquake induced settlements of the combined piled raft foundation and, also, crucial need for quite a speedy and inexpensive predictive method, the applicability of neural network modeling has been investigated and compared.Moreover, the large calculation times available in commercial finite element or finite difference software programs for dynamic complicated problems makes the problem time-consuming.Hence, a speedy and accurate tool for predicting these values will significantly help the engineers to facilitate the design procedure.In this regard, the highly in-demand predictive tools, ANN, and EPR modeling have been considered to forecast the settlements in both static and coupled dynamic-static loading conditions.In this regard, dynamic nonlinear finite difference modeling has been used to prepare the required datasets.After wellvalidating the used numerical procedure and applying a parametric study, a databank consisting of 144 datasets has been gained.Using the developed databank, different ANNs and polynomial regression models have been constructed, which can be efficiently used to predict the values of settlements for all the cases that there are similar available input parameters.Moreover, using the sensitivity analysis, the most and the least effective parameters on the static and coupled static-dynamic settlements have been determined.Based on the obtained results, the following can be concluded: (i) Neural networks are introduced as capable effective tools for predicting static and coupled static-dynamic settlements of CPRFs.
(iii) In the case of coupled static-dynamic loads, the network with the architecture of 5-8-8-1, RMSE of 0.039,  2 of 0.97, and  of 0.983 is the best obtained neural network.
(iv) Based on EPR modeling, a new relationship with coefficient of determination of 97.49% and RMSE of 0.271 and the line slope of 0.975 for predicting the static settlements of combined piled raft foundations has been proposed.
(v) A new relationship with coefficient of determination of 97.26% and RMSE of 0.137 and the line slope of 0.973 for predicting the coupled static-dynamic settlements of combined piled raft foundations has been presented.
(vi) Sensitivity analysis emphasized that axial load plays an important role in calculating static/coupled staticdynamic settlements and neglecting its effects will lead to significant errors in calculations.It was also shown that the effect of / on the settlements is less than other concerning parameters.
Based on the results, some efficient and highly capable tools for predicting static and coupled static-dynamic settlements of CPRFs are suggested, which can be successfully used to predict the values of settlements in preliminary step of CPRF design.

Figure 5 : 2 SFigure 6 :
Figure 5: Two different grid sizes used to carry out the sensitivity analysis on the mesh sizes.

Figure 9 :
Figure 9: (a) Geometry of conducted physical model.(b) Three-dimensional finite difference mesh.(c) Comparison of the results of finite difference and physical model [6].

4. 2 .
Effect of Pile's Length on the Settlements of CPRF.To show the effect of piles' length, as sample results, settlements of different CPRFs are shown in Figure 13.Figures 13(a), 13(b), and 13(c), respectively, describe the effect of piles' length on the coupled static-dynamic settlements of the studied 3 * 3, 4 * 4, and 5 * 5 CPRFs under 60 kPa static axial load for 0.5 meters in diameter piles.

( 2 )
Performance of Neural Networks in the Combined Static-Dynamic Loading Case.The results of evaluations of the

Figure 22 :
Figure 22: The calculated values of static settlement versus its corresponding EPR prediction.

Figure 23 :
Figure 23: The calculated values of coupled static-dynamic settlements versus its corresponding EPR prediction.

Figure 25 :
Figure 25: Strength of relationship between input parameters and CPRF settlements in the coupled static-dynamic case.

Table 1 :
Parametric study on the whole affecting parameters for static loading condition.

Table 2 :
Parametric study on the whole affecting parameters for combined static-dynamic loading condition.

Table 3 :
Performance of neural networks in the static case.

Table 4 :
Performance of neural networks in the coupled static-dynamic loading case.