Research Article  Open Access
Yunfei Xie, Shichun Chi, "Optimization Method of Reducing the Differential Settlements of Piled Raft Foundations Based on PiletoPile Interaction Theory", Advances in Civil Engineering, vol. 2020, Article ID 1521876, 14 pages, 2020. https://doi.org/10.1155/2020/1521876
Optimization Method of Reducing the Differential Settlements of Piled Raft Foundations Based on PiletoPile Interaction Theory
Abstract
In order to reduce the differential settlement of piled raft foundations, an optimization method based on piletopile interaction theory is proposed in this paper, which translates the problem of piletopile interaction (PPI) in pile groups into that in single piles using the interaction factor method. The pile lengths were adjusted via the relationship between load, settlement, and the length of single piles during the optimization design. ANSYS software, in conjunction with nonlinear elastic soil model, is used to analyze piled raft foundation models. Two cases with different safety factors that suffer different kinds of surface loads (uniform load and nonuniform load) are used to verify this method. The differential settlements of the raft in different cases are all reduced by nearly or more than 80% after optimization design. The results show that the optimization method proposed in this paper has high efficiency and stability. This study can help practicing engineers optimize the pile lengths in pile groups to satisfy higher differential settlement requirements.
1. Introduction
Piled raft foundations (PRFs) have been used for many years as a method of transmitting the loads from the superstructure to competent strata at depth in the ground. During the traditional design of pile foundations, a large safety factor is always used to control the settlements and differential settlements. The stiffness distribution of the pile group is always uniform for the convenience of design and construction. However, with the gradual improvement of building requirements for the control of settlement and differential settlement, it is very wasteful to meet the requirements by increasing the safety factor for the reason that piletopile interaction and pilesoilpile interaction always exist in PRFs [1–3]. More recently, the focus has shifted towards more economical and reasonable design, thereby prompting considerable research effort to optimize the piles in the group.
Optimization design of PRFs with dissimilar piles (with different lengths and diameters) [4, 5] and pile arrangement [6] has drawn much attention in recent years, which can minimize the differential settlement on top of the pile group and rearrange the load distribution of a PRF. The critical design criteria are often associated with the differential settlements especially for large pile groups and important buildings [7–11]. To address this issue, Padfield and Sharrock [12] enhanced pile group performance by changing the pile arrangements across the group. Nguyen et al. [13, 14] gave a parametric study for the optimal design of large PRFs on sand, and the results show that differential settlement could be effectively reduced by putting more piles in the center of a raft than in other areas. Leung et al. [15] gave an analysis and optimization method for PRFs which allows the coupling of the stiffness of superstructure with the foundation model. Algin [16] carried out a multiobjective optimization analysis method using the response surface method (RSM) to achieve the most economical design solution that satisfies the present design constraints for JGRs.
These aforementioned studies, however, have some limitations in the convenience of the application of optimization methods and the accuracy of optimization results. Pile layout optimizations proposed via outwardly powerful and inwardly weak or outwardly weak and inwardly powerful are easy to achieve, but the accuracy of optimization results is difficult to control. Optimization designs based on traditional gradient methods or genetic algorithm methods can easily result in accurate optimization results, but they are restricted to design problems that they must be mathematically formulated and differentiable and need expensive computation (especially for largescale pile groups). In the view of construction convenience, variable pile length is undoubtedly the least difficult way compared with the design using variable pile diameter and variable pile spacing because it avoids the increase of positioning complexity caused by different pile spacing and the constructing trouble caused by different pile diameters. In order to achieve a variable stiffness optimization design using different pile lengths, it is necessary to clarify the “loadsettlementpile length” (LSP) relationship of pile groups. It is easy to establish LSP mathematical relationship for a single pile, but difficult for pile groups. Therefore, the key of this paper is how to transform the problems in pile group optimization into the optimization problems among single piles combined with PPI, to achieve the optimization design with dissimilar pile length (PPI optimization method).
An important result of research into the settlement behavior of pile foundations under vertical loads is that predictions only based on elastic theory are excessively conservative, and since the surrounding soil is always assumed to be single phase, the consolidating pilehead settlement is often neglected when a pile foundation is located in deep and saturated soil [17]. Therefore, nonlinear soil behavior is considered to be governed by Duncan–Chang (DC) model [18] in this paper. Within a nonlinear framework, pile settlement is coupled with the ultimate capacity, and therefore, an accurate estimation of a pile capacity is a prerequisite for a rigorous analysis of the serviceability limit state. While simplified design approaches an integral part of preliminary foundation design, threedimensional (3D) nonlinear continuum analyses are becoming more commonplace. Therefore, 3D finite element method (FEM) is used to estimate the bearing capacity of single piles and the loads shared by each pile in the groups. Through the relationship of the settlement of single piles and that in pile group, which is established in the PPI optimization method, each pile length can be adjusted according to the corresponding pile top settlement of the pile group. High calculation cost and too much iteration times are avoided during this process. The optimization results show that the differential settlements of the pile groups in different cases can be highly reduced after optimization design compared with the traditional design no matter what the load is uniform or nonuniform.
2. Method
Based on the verification of the simulation of a single pile load test via FEM, the single pile loadsettlement (LS) curve and the ultimate bearing capacity (UBC) in different pile lengths with the same soil can be estimated. The LSP surface of single piles is obtained via interpolating the PS curves of several single piles with MATLAB. The required pile length can be determined by any given pile top load P_{i} and the corresponding settlement S_{ii} of single pile via LSP. However, LSP is only used for single pile settlement problems because it neglects the complex influence in pile groups. Therefore, LSP cannot be applied to the optimization design of pile group directly, and a relationship between pile group and single piles is necessary to be established.
2.1. Establishment of PiletoPile Influence Coefficient (β_{ij}) in Pile Groups
PPI effects, which especially exist between pile shafts, necessitate rigorous treatment, as settlements are amplified compared with single pile values [19, 20]. This necessitates more reliable modelling of the development of pilesoil interface resistance during loading as well as more realistic treatment of PPI [21, 22].
2.1.1. Introduction of PPI Factors
The interaction factor method (IFM) [23, 24] is used for accounting PPI in the design of pile groups. Interaction factor α_{ij} [25, 26] for pile i due to pile j, corresponding to the spacing between piles i and j, is given as follows:where S_{ii} is the settlement of a single pile i under its own load (P_{i}); ΔS_{ii} = S_{ii}/P_{i} is the settlement of a single pile i under unit load; S_{ij} is the additional settlement of pile i due to adjacent pile j loaded by P_{j}; and ΔS_{ij} = S_{ij}/P_{j} is the additional settlement of pile i due to adjacent pile j loaded by unit load. The total settlement of pile i, which satisfies the principle of superposition, is the sum of its own settlement (S_{ii}) and additional settlement (S_{ij}) due to pile j [27]:
2.1.2. Application of PPI Factor in Pile Groups
Some researchers consider that the use of interaction factors based only on the spacing of any pair of piles in a group may result in ignoring the potential reinforcing effects of center piles [28–31]. From this view, the principle of superposition is not valid in nonlinear engineering problems. Nevertheless, some other research [32] shows that the soil remains in a linear elastic state outside this narrow region near the pilesoil interface. Therefore, PPI is governed by the elastic behavior in this paper.
Superposition of twopile interaction factors is used to predict full pile group behavior. In this method (referring to Figure 1), the settlement S_{i} of pile i within a group of N piles can be obtained by the IFM, which is given as follows [27]:where
β_{i} is defined as the total influence on pile i affected by the other piles in the pile group. Similarly, there is β_{j} = S_{j}/S_{jj}. The influence coefficient of pile j on pile i in the pile group is defined as follows:
Equation (5) shows that the influence coefficient between any two piles in the pile group can be expressed by the single pile settlements of the two piles and the settlements of the two piles in a pile group under the same load condition. This paper assumes that the β_{ij} only depends on the spacing between piles i and j and will not change with the change of pile length.
2.2. PPI Optimization Method
2.2.1. Optimization Design Based on β_{ij}
In order to solve the differential settlement problem in a PRF caused by pilesoil interaction and different loads on top of the piles, variable pile lengths are used in this paper to decrease the differential settlement. The parameters with “′” are defined as the optimization target parameters. For example, S_{ii}′ and S_{i}′ are the target settlement of a single pile i and the target settlement of pile i in a pile group, respectively.
The reference pile (Rpile) needs to be determined at first, and the other piles are optimized with the settlement of Rpile as the reference line. Therefore, if pile j is chosen as the Rpile, there is S_{jj} = S_{jj}′. S_{i}′ = S_{j}′ = S_{j} for any pile i in the pile group. The influence coefficient is unchanged because of the unchanged pile spacing during the optimization procedure, and therefore, β_{ij}′ = β_{ij} is given by
Equation (6) shows that the target single pile settlement of any pile i (S_{ii}′) can be obtained through the single pile settlement (S_{jj}) of Rpile and the influence coefficient (β_{ij}) of pile i due to Rpile. The influence coefficients of Rpile on other piles in the pile group are different and can be expressed as matrix [β]_{R}. The target settlement (S_{ii}′) of each single pile is also different and can be expressed as matrix [S]_{t}. Therefore, equation (6) can be expressed as follows:
The length of each pile ([L]) in the pile group after optimization design can be obtained with [S]_{t}, [P], and LSP surface.
2.2.2. Calculate [β]_{R}
The load distribution on top of piles in a pile group is affected by raft stiffness. Even if the load is uniform, the load can be unevenly distributed on top of the piles due to the pilesoil interaction. Moreover, the load distribution mechanism of pile top under uneven load is more complex. Therefore, it is difficult to calculate [S]_{t} directly. In order to solve this problem, a PRF system is established using 3D FEM, which considers the influence of pile group effect and raft stiffness on load distribution. After analysis of the PRF system, the loads ([P]) and settlements ([S]_{group} = [S_{1},…, S_{i},…, S_{N}]) on top of the piles can be extracted directly, and then, the settlements of single piles ([S]_{single} = [S_{11},…, S_{ii},…, S_{NN}]) under corresponding loads ([P]) can be extracted by LSP surface. Finally, each influence coefficient (β_{ij}) can be obtained by equation (5).
2.3. Optimal Iteration
Optimization efficiency can be improved by recycling the PPI optimization method. After the m_{th} iteration, the load distribution ([P]_{m}), settlements of the piles ([S]_{group,m}), and pile lengths ([L]_{m}) of the optimized model are extracted, and the corresponding [S]_{single,m} can be obtained. Then, renew [β]_{R} and [S]_{t} by equations (5) and (7). Finally, the new pile lengths ([L]_{m+1}) are obtained by LSP surface, where m represents the number of iterations and m = 0 in the initial design.
The different optimized pile lengths need to be merged and combined after optimization design, and the error will always exist. Therefore, this paper agrees that if the increment (ΔFF_{m+1} = FF_{m+1} − FF_{m}) of optimum rate (FF) of differential settlement is less than 5%, the optimization effect tends to be stable and the optimization cycle can be stopped. FF represents the optimization rate and is defined as follows:where R is the variable of the initial design and R′ is that of the optimized design.
2.4. Suggestions and Constraint Condition
The corresponding SF in a different initial design will produce a different differential settlement. In order to reduce the differential settlement more effectively, this paper assumed that there is a linear correlation between the maximum settlement of the pile group and the corresponding safety factor (SF) as shown in Figure 2, and three suggestions about the selection of Rpile are given as follows:(a)When , it belongs to conservative design and means that the optimization space of the PRF is sufficient. It is suggested that the pile with the largest settlement should be taken as Rpile, and the whole optimization design is mainly to reduce the length of piles.(b)When , it belongs to the low safety factor design. The pile with the smallest settlement should be taken as Rpile, and the optimization design is mainly to increase the length of piles.(c)When , the pile corresponding to the moderate settlement is taken as Rpile, and the length of the piles can be increased or reduced at the same time in the whole optimization design. The pile corresponding to this settlement (S in Figure 2) is taken as the Rpile.
In this paper, the SF is defined as the UBC of the pile groups divided by an applied load as shown in equation (9), and it depended on the initial design of the PRF (e.g., the pile length, pile diameter, and pile configuration). In order to facilitate the calculation, the load component of raft [33] is not taken into account within the UBC:where Q_{uk,i} is the UBC of single pile i in the initial design, F_{load} is the total load on top of the raft, L_{i} is the length of pile i, and L_{0i} is the length of pile i in the initial design. Due to the influence of PPI of longshort piles, light microwave vibration of the length of Rpile is unavoidable. The length of Rpile can be adjusted according to equation (10) to ensure that the Rpile is always in a stable and unchanged reference state.
Different pile lengths need to be merged and combined after optimization design in consideration of the feasibility and convenience of construction. The minimum value of any L_{i}−L_{j} needs to be no less than 0.5 m for the reason of too small a difference of pile lengths (e.g., L_{i}−L_{j} <0.5 m) will affect the feasibility of construction and then affect the optimization result.
This paper only studies the optimization design of pile groups in the elastic stage, and then, it is necessary to ensure that the load on top of each pile does not exceed its UBC during the optimization process. The load on top of each pile after optimization design should satisfy the following condition:
If P_{i} does not satisfy equation (11), the safety and affluence of the pile group are insufficient and the selected Rpile in this cycle is unreasonable and needs to be reselected. The largest P_{i,max} of the piles which does not satisfy equation (11) is selected as a new Rpile, and then, the previous cycle should be reoptimized according to SF = 1.
The flowchart of the PPI optimization method is shown in Figure 3.
3. Illustrative Example
3.1. Piled Raft Model (PRF) and Material
Two design cases are given in this paper to investigate the foregoing method. Both of the PRFs in the two cases consist of 121 concrete piles (N = 121) arranged in a 11 × 11 pattern and connected with a raft (49 m × 49 m). The only difference is that the piles are 37 m long in Case 1 and 33 m long in Case 2. Pile diameters are all 1.5 m located with the same centertocenter spacing of 4.5 m. The calculation range of the soil in the horizontal direction is approximately 4 times the length of the raft. Generally, in order to reduce the influence of the boundary effect, the distance from the pile tip to the bottom of the model is 2 times the pile length [34, 35]. In this study, the distance between the pile tip and the top of the bedrock is less than 2 times the pile length while the bedrock can be considered as a fixed end, so the bottom boundary of the model is located 5 m below the top of the bedrock. Table 1 and Figure 4 summarize the pile group parameters used in the analyses. For the calculation of loadsettlement responses in piles, Moayedi and Rezaei [36] and Moayedi and Hayati [37] used neural network, and the software of ANSYS was used in this paper. In the ANSYS analyses, the nodes are shared along the pilesoil, pileraft, and soilraft interfaces, and zero slippage at these interfaces is assumed. The maximum element sizes in the horizontal direction and the vertical direction are 6 m and 3 m, respectively. This results in about 530000 nodes and 840000 elements.

(a)
(b)
The soil model is assumed as an incremental nonlinear (hyperbolic) stressdependent model (DC model), which is based on Hooke’s law, and the main feature of which is a nonlinear relationship between stress and strain. Before the calculation, the initial Young modulus of each soil layer is calculated using the initial principal stresses (σ_{1} and σ_{3}), where the principal initial stresses of each element are obtained from equation (12). K_{0} is the static lateral pressure coefficient, z is the depth of the soil element center to the ground, and γ is the gravity per unit volume of the soil:
A detailed description about the DC model can be found in Duncan and Chang [18] and Chen et al. [38]. The corresponding 8 DC model parameters (K, n, φ, c, R_{f}, G, F, and D) of the soil adopted in this paper are listed in Table 2 [39]. This paper assumes that the soil with the same property but different thickness has the same DC model parameters.

3.2. Boundary Conditions
All soil nodes at the bottom of the model (with a Z coordinate of −26 m) are considered to be constrained in three directions, i.e., the displacements of the nodes in X, Y, and Z directions are equal to zero. The outboard lateral boundary nodes only have displacement in Z direction, and the motions in X and Y directions are constrained.
3.3. Loading and Analysis Procedure
In order to investigate the universality of the PPI method, two kinds of loads are used in this paper: uniform load and nonuniform load. The distribution of the two loads is shown in Figure 5. The total load (F_{load}) on top of the raft is 900 MN in Case 1 and is 1100 MN in Case 2. There are 21 load steps during the analysis procedure. The first step (i = 1) is the gravity of the PRF system, and the external load on top of the raft is divided into 20 steps (step i from 2 to 21) in order to consider the nonlinear characteristics of soil.
4. Optimization Design via PPI Optimization Method
As mentioned above, the LSP surface of single piles is an important point during the whole optimization procedure. Therefore, the simulation of a single pile load test via FEM is needed to verify the validity of the analysis method.
4.1. Verification
A single pile load test (test pile) from Jiang [39] is simulated in this paper. The pile in the single pile model has a circular cross section with a diameter of 2.8 m and a depth of 125 m (the pile length between +4 m and −46 m is not in the soil). The soil model, which is composed of 9 soil layers, has a diameter of 100 m and a depth of 174 m. The minimum and the maximum sizes of the model are 1 m and 5 m, respectively. This results in about 19000 elements and 22000 nodes. The nodes along the pilesoil interfaces are shared, assuming that there is no slippage at these interfaces. Only a quarter of the actual model was modeled due to its symmetry. The other boundary conditions are the same as the PRF model. The load applied on top of the pile, which is 5 MPa, is divided into 20 steps, increasing from 0 MPa to 5 MPa with an increment load equal to 0.25 MPa.
The concrete pile has Young’s modulus of E = 3.3 × 10^{10} Pa, an effective unit weight of γ = 1.4 MN/m^{3}, and Poisson’s ratio of ν = 0.12. The elastic constitutive model is used to model the concrete pile in this paper. The behavior of soil layers is modeled using the DC model. The corresponding soil parameters of the single pile test are shown in Table 3.

In Figure 6, the settlement results of the single pile load test are compared with those of the analytical prediction. The results indicate that the simulation results are in good agreement with experimental results, and the settlement error is only 3% under the same ultimate load, as shown in Table 4. That means the FEM adopted in this paper can simulate the stress state of the pile before its UBC. Before the load reaches the ultimate load, a sharp subsidence of the pile will not appear; i.e., the slippage of the interface between pile and soil is not severe. Therefore, unless the load on top of the pile exceeds its UBC, the assumption of rigid contact between pile and soil is valid.

4.2. Getting LSP Surface
The single pile loadsettlement with different pile lengths in the soil layers of PRF can be simulated via the same simulation method in Section 4.1, and the corresponding loadsettlement curves can be obtained (see Figure 7(a)). These different loadsettlement curves of different piles length can be translated into LSP surface via MATLAB software, as is shown in Figure 7(b).
(a)
(b)
The UBCs of the piles in different lengths need to be calculated by other methods because the simulation method can only simulate the state of pile top load before reaching the ultimate load. This paper uses empirical formulas (equation (13)) from the Chinese code [40] to estimate the UBC:where Q_{uk} is the UBC of a single pile; u is the pile perimeter; A_{p} is the section area at the tip of the pile; l_{i} is the i_{th} soil layer’s thickness; q_{sik} and q_{pk} are the ultimate shaft and tip resistance of i_{th} soil layer, respectively; and and are the effect coefficients of the shaft and tip resistance for piles in large diameter (diameter greater than 0.8 m), respectively. The corresponding empirical parameters used in equation (13) of different soil are shown in Table 5.
 
Note. The tip resistances (q_{pk}) of the soil take the lower limit of the given range in the Chinese code [35]. 
The empirical value of UBC of the test pile is 43.6 MN, with only 3% error compared with the field test value (45 MN). This indicates that the empirical parameters of the soil layers can be used to estimate the UBC of the other piles which have different pile lengths but in the same soil types (the thickness of the soil layer can be different). The empirical values of UBCs of the piles in different lengths in Figure 7(a) are shown in Figure 8. The turning point of the straight line in Figure 8 is caused by soil stratification.
4.3. Optimization Design of the Cases
The UBC of a single pile in 37 m long is 17.2 MN, which can be determined by equation (13) or Figure 8. The safety factor SF of the PRF in Case 1 is more than that in Case 2 (SF is 2.3 under uniform load and 2.2 under nonuniform load). Therefore, the pile with the maximum settlement in the initial design is selected as Rpile (the pile in the circle), as is shown in Figure 9. Finally, the optimization process of Case 1 converges after 4 iterations under uniform load and 3 iterations under nonuniform load. The final pile lengths and the corresponding arrangements are shown in Table 6 and Figure 9. The pile length varies from 33.5 m to 37 m under the uniform load and from 31.5 m to 37 m under nonuniform load. The purpose of this paper is to reduce the differential settlement, and the total volume of material will change according to the optimization results under different SF, rather than remain constant.
(a)
(b)

For Case 2, the UBC of a single pile in 33 m long is 14.1 MN. The safety factor SF is 1.55 under uniform load and 1.5 under nonuniform load. 1< SF <2, therefore, the pile with settlement close to 0.118 m is selected as Rpile under uniform load and the pile with settlement close to 0.121 m is selected as Rpile under nonuniform load, as is shown in Figure 10 (Rpile is the pile in the circle). The optimization process converges after 2 iterations under both uniform load and nonuniform load. The pile length and the corresponding arrangements are shown in Table 7 and Figure 10. The pile length changes from 28 m to 35 m under uniform load and from 27 m to 35 m under nonuniform load.
(a)
(b)

5. Results and Discussion
Differential settlements can result in an additional bending moment of the raft and then increase the cost of the raft. Therefore, the settlement of the raft is chosen as the main analysis objective in this paper to investigate the effectiveness of the optimization method.
5.1. Result of Case 1
5.1.1. Settlement of the Raft
The settlements on top of the raft in Case 1 under different loads before and after optimization design are compared in Figures 11 and 12. The result shows that the settlement surface of the raft becomes more uniform after optimization design compared with the initial design, and the maximum settlement of the raft is moved out from the center of the raft.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
Differential settlement is defined as the difference between the maximum and minimum settlement in this paper. Table 8 gives a detailed comparison of the differential settlement of the raft (DS raft) under different conditions. The DS raft of Case 1 is reduced by 80% from 0.015 m to 0.003 m under uniform load and by 87% from 0.02 m to 0.0025 m under nonuniform load.
 
Note. IU = initial design under uniform load; IN = initial design under nonuniform load. 
5.1.2. Verification of UBC of Piles in Different Lengths
The comparison of the maximum pile top loads of different pile lengths in different models and their corresponding UBCs is shown in Tables 9 and 10. The result indicates that no pile exceeds the corresponding UBC during the optimization procedure, which can also be verified in Figure 8.
 
Note. P_{max} = the maximum load on top of the piles with the same pile length. 
 
Note. P_{max} = the maximum load on top of the piles with the same pile length. 
5.2. Result of Case 2
5.2.1. Settlement of the Raft
Figures 13 and 14 give the comparison of the settlement on top of the raft in Case 2 under different loads before and after optimization design. The result is similar to that in Case 1, where the settlement surface of the raft is more uniform by using different pile lengths compared with the initial design, and the maximum settlement of the raft is moved out from the center of the raft. A detailed comparison of the DS rafts in Case 2 under different loads is shown in Table 11. The DS raft is reduced from 0.047 m to 0.01 m under uniform load, by 79%, and from 0.051 m to 0.0076 m under nonuniform load, by 85%.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
 
Note. IU = initial design under uniform load; IN = initial design under nonuniform load. 
5.2.2. Verification of Ultimate Bearing Capacity of Piles in Different Lengths
The maximum pile top loads of different pile lengths are compared with their UBCs in Tables 12 and 13. The result indicates that no pile exceeds the corresponding UBC during the optimization procedure.
 
Note. P_{max} = the maximum load on top of the piles with the same pile length. 
 
Note. P_{max} = the maximum load on top of the piles with the same pile length. 
5.3. Discussion
Pilesoil interaction is the main reason leading to different pile lengths when the load is uniform, and the difference between different pile lengths is usually not too large. Therefore, the effect of the error due to the merging and combining of pile lengths is more prominent, and then, more iterations are needed to make up for it. When the load is nonuniform, the difference of different pile lengths under nonuniform load is more obvious after optimization compared with that under uniform load. Even though some differences of the pile lengths are eliminated by merging and combining, the influence on the optimization result by this error is not as obvious as that under uniform load, and then fewer iterations are needed, which can be seen from Tables 8 and 11.
The safety factor in Case 1 is more than 2, and it is between 1 and 2 in Case 2. Increasing and reducing pile lengths can be adopted at the same time in Case 2 during the optimization process as is suggested before, while Case 1 can only choose reducing the pile lengths. Therefore, the models in Case 1 need more iterations compared with those in Case 2 (see Figure 9 and Table 6). The DS raft always shows a higher optimization efficiency under nonuniform load because of the merging and combing of pile lengths compared with that under uniform load.
Compared with the optimization method [41] proposed by the authors, the optimization method proposed in this paper is more effective: (1) the FF of differential settlement of piled raft foundation is significantly improved; (2) the realization of optimization method is more automatic, avoiding the error caused by manual participation; and (3) it is easier to change the pile length for construction compared with change pile diameter and pile spacing.
6. Conclusions
This paper proposes an optimization method based on PPI theory aiming to remove the differential settlement of the PRFs. ANSYS package software is used in conjunction with DC nonlinear elastic soil model to simulate PRFs under different surface loads to investigate the optimization efficiency of the PPI optimization method. After the verification of a single pile load test simulation, the settlement load curves of single piles with different lengths via the same simulation method are obtained to establish LSP surface of single piles. The results emphasize the superior advantages (the optimization rate FF is nearly more than 80% without too many iterations) of the PPI optimization method over the initial design by the traditional method in reducing the differential settlements of PRFs. Based on the results, other conclusions are proposed:(1)The problem of unknown pile top loads of pile groups in the application of conventional IFM is solved via FEM. The simplified method, which defined raft as rigid and flexible, is avoided.(2)According to Poulos’s formula [27] for calculating the settlement of a pile in pile group, which is based on the superposition method, the influence coefficient (β_{i}) of other piles in pile group on pile i is obtained. The relationship between the settlement of a pile in a single pile system and the settlement of the pile in a pile group system under the same pile top loads is established.(3)The results before and after optimization show that the PPI optimization method has high optimization efficiency. The DS raft is reduced by nearly 80% under uniform load and more than 85% under nonuniform load, which results in higher optimization rate compared with published works [5, 41]. It can be seen that the uneven distribution of the load does not affect the efficiency of the optimization method proposed in this paper.(4)The traditional uniform pile layout design is selected as the initial design to facilitate comparison in this paper. However, the optimization method does not limit the load type and the initial design of the pile arrangement, as is shown in Figure 3. The reason is the optimization of each pile is independent in the optimization process. Therefore, this method has strong applicability, especially in the optimization design with the differential settlement of the pile foundation as the control objective.
The selection of the Rpile is a key point of this research method, which depends on the SF of the initial design of the pile group and determines the bearing capacity level of PRF after the optimized design. The SF is directly proportional to the economic investment of pile foundation; therefore, the proposed method does not directly incorporate economic optimization into the optimization objective, but gives different optimization suggestions for PRFs under different safety factors. In order to further study the optimization method in the near future, the centrifuge test or smallscale test will be studied to more accurately assess the performance of the model for pile groups and piled raft, and the parametric analysis will be studied to research the importance of parameters. Furthermore, with the continuous development of artificial intelligence or artificial neural network, future research efforts also can be focused on the combination of artificial intelligence and optimization design.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
There are no conflicts of interest regarding the theoretical and numerical observations introduced in this paper.
Acknowledgments
This work was supported by the National Key Research and Development Project in China (grant no. 2016YFB0201001).
Supplementary Materials
The supplementary material is the database for the Figures 11–14 within the manuscript. (Supplementary Materials)
References
 A. T. Ghalesari, A. Barari, P. F. Amini, and L. B. Ibsen, “Development of optimum design from static response of pileraft interaction,” Journal of Marine Science and Technology, vol. 20, no. 2, pp. 331–343, 2015. View at: Google Scholar
 H. G. Poulos, “Pile behaviourtheory and application,” Géotechnique, vol. 39, no. 3, pp. 365–415, 1989. View at: Publisher Site  Google Scholar
 A. T. Ghalesari and H. Rasouli, “The effect of gravel layer on the behavior of piled raft foundations,” in Proceedings of the Advances in Soil Dynamics and Foundation Engineering, pp. 373–382, ASCE, Shanghai, China, May 2014. View at: Publisher Site  Google Scholar
 F. Liang, L. Chen, and J. Han, “Integral equation method for analysis of piled rafts with dissimilar piles under vertical loading,” Computers and Geotechnics, vol. 36, no. 3, pp. 419–426, 2009. View at: Publisher Site  Google Scholar
 Y. F. Leung, A. Klar, and K. Soga, “Theoretical study on pile length optimization of pile groups and piled rafts,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 136, no. 2, pp. 319–330, 2010. View at: Publisher Site  Google Scholar
 K. N. Kim, S.H. Lee, K.S. Kim, C.K. Chung, M. M. Kim, and H. S. Lee, “Optimal pile arrangement for minimizing differential settlements in piled raft foundations,” Computers and Geotechnics, vol. 28, no. 4, pp. 235–253, 2001. View at: Publisher Site  Google Scholar
 L. D. Sanctis, A. Mandolini, G. Russo, and C. Viggiani, “Some remarks on the optimum design of piled rafts,” in Proceedings of the International Deep Foundations Congress, pp. 405–425, ASCE, Orlando, FL, USA, February 2002. View at: Google Scholar
 Y. L. Li, J. F. Wang, Q. Wang, and K. Yang, “Optimization design study on variable stiffness of the pile raft foundation in highrise building,” Advanced Materials Research, vol. 243249, pp. 2498–2502, 2011. View at: Publisher Site  Google Scholar
 M. M. Sales, M. Prezzi, R. Salgado, Y. S. Choi, and J. Lee, “Loadsettlement behaviour of model pile groups in sand under vertical load,” Journal of Civil Engineering and Management, vol. 23, no. 8, pp. 1148–1163, 2017. View at: Publisher Site  Google Scholar
 S. Shrestha, N. Ravichandran, and P. Rahbari, “Geotechnical design and design optimization of a pileraft foundation for tall onshore wind turbines in multilayered clay,” International Journal of Geomechanics, vol. 18, no. 2, 2017. View at: Publisher Site  Google Scholar
 A. T. Ghalesari and A. J. Choobbasti, “Numerical analysis of settlement and bearing behaviour of piled raft in Babol clay,” European Journal of Environmental & Civil Engineering, vol. 22, no. 8, pp. 978–1003, 2016. View at: Publisher Site  Google Scholar
 C. J. Padfield and M. Sharrock, Settlement of Structures on Clay Soils, Construction Industry Research & Information Institute, London, UK, 1983.
 D. D. C. Nguyen, D. S. Kim, and S. B. Jo, “Settlement of piled rafts with different pile arrangement schemes via centrifuge tests,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 139, no. 10, pp. 1690–1698, 2013. View at: Publisher Site  Google Scholar
 D. D. C. Nguyen, D.S. Kim, and S.B. Jo, “Parametric study for optimal design of large piled raft foundations on sand,” Computers and Geotechnics, vol. 55, pp. 14–26, 2014. View at: Publisher Site  Google Scholar
 Y. F. Leung, A. Klar, K. Soga, and N. A. Hoult, “Superstructurefoundation interaction in multiobjective pile group optimization considering settlement response,” Canadian Geotechnical Journal, vol. 54, no. 10, pp. 1408–1420, 2017. View at: Publisher Site  Google Scholar
 H. M. Algin, “Optimised design of jetgrouted rafts subjected to nonuniform vertical loading,” KSCE Journal of Civil Engineering, vol. 4, pp. 1–15, 2017. View at: Google Scholar
 T. Senjuntichai, N. Sornpakdee, J. Teerawong, and R. K. N. D. Rajapakse, “Timedependent response of an axially loaded elastic bar in a multilayered poroelastic medium,” Journal of Engineering Mechanics, vol. 133, no. 5, pp. 578–587, 2007. View at: Publisher Site  Google Scholar
 J. M. Duncan and C. Y. Chang, “Nonlinear analysis of stress and strain in soils,” Journal of the Soil Mechanics and Foundations Division, vol. 96, no. 5, pp. 1629–1653, 1970. View at: Google Scholar
 M. Modarresi, H. Rasouli, A. T. Ghalesari, and M. H. Baziar, “Experimental and numerical study of piletopile interaction factor in sandy soil,” Procedia Engineering, vol. 161, pp. 1030–1036, 2016. View at: Publisher Site  Google Scholar
 A. T. Ghalesari, A. Barari, P. F. Amini, and L. B. Ibsen, “The settlement behavior of piled raft interaction in undrained soil,” in Proceedings of the IACGE 2013: Challenges and Recent Advances in Geotechnical and Seismic Research and Practices, pp. 605–612, ASCE, Chengdu, China, October 2013. View at: Publisher Site  Google Scholar
 S. C. Wong and H. G. Poulos, “Approximate piletopile interaction factors between two dissimilar piles,” Computers and Geotechnics, vol. 32, no. 8, pp. 613–618, 2005. View at: Publisher Site  Google Scholar
 Q. Zhang and Z. Zhang, “Study on interaction between dissimilar piles in layered soils,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 35, no. 1, pp. 67–81, 2011. View at: Publisher Site  Google Scholar
 H. G. Poulos, “Analysis of the settlement of pile groups,” Géotechnique, vol. 18, no. 4, pp. 449–471, 1968. View at: Publisher Site  Google Scholar
 H. G. Poulos, “Modified calculation of pile‐group settlement interaction,” Journal of Geotechnical Engineering, vol. 114, no. 6, pp. 697–706, 1988. View at: Publisher Site  Google Scholar
 B. A. McCabe and B. B. Sheil, “Pile group settlement estimation: suitability of nonlinear interaction factors,” International Journal of Geomechanics, vol. 15, no. 3, 2015. View at: Publisher Site  Google Scholar
 B. B. Sheil, B. A. Mccabe, E. M. Comodromos, and B. M. Lehane, “Pile groups under axial loading: An appraisal of simplified nonlinear prediction models,” Géotechnique, vol. 69, no. 7, pp. 565–579, 2019. View at: Publisher Site  Google Scholar
 H. G. Poulos, “Pile group settlement estimation—research to practice,” in Proceedings of the Foundation Analysis and Design, pp. 1–22, ASCE, Shanghai, China, June 2006. View at: Publisher Site  Google Scholar
 B. E. Sharnouby and M. Novak, “Stiffness constants and interaction factors for vertical response of pile groups,” Canadian Geotechnical Journal, vol. 27, no. 6, pp. 813–822, 1990. View at: Publisher Site  Google Scholar
 P. H. Southcott and J. C. Small, “Finite layer analysis of vertically loaded piles and pile groups,” Computers and Geotechnics, vol. 18, no. 1, pp. 47–63, 1996. View at: Publisher Site  Google Scholar
 G. Mylonakis and G. Gazetas, “Settlement and additional internal forces of grouped piles in layered soil,” Géotechnique, vol. 48, no. 1, pp. 55–72, 1998. View at: Publisher Site  Google Scholar
 F. Basile, “Analysis and design of pile groups,” in Numerical Analysis and Modelling in Geomechanics, pp. 278–315, CRC Press, Boca Raton, FL, USA, 2003. View at: Google Scholar
 Y. F. Leung, K. Soga, B. M. Lehane, and A. Klar, “Role of linear elasticity in pile group analysis and load test interpretation,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 136, no. 12, pp. 1686–1694, 2010. View at: Publisher Site  Google Scholar
 J. H. Lee, Y. Kim, and S. Jeong, “Threedimensional analysis of bearing behavior of piled raft on soft clay,” Computers and Geotechnics, vol. 37, no. 12, pp. 103–114, 2010. View at: Publisher Site  Google Scholar
 M. Bryan and S. Bria, “Pile group settlement estimation: suitability of nonlinear interaction factors,” International Journal of Geomechanics, vol. 15, no. 3, 2015. View at: Publisher Site  Google Scholar
 H. M. Algin, “Optimised design of jetgrouted rafts subjected to nonuniform vertical loading,” KSCE Journal of Civil Engineering, vol. 22, no. 2, pp. 494–508, 2018. View at: Publisher Site  Google Scholar
 H. Moayedi and A. Rezaei, “An artificial neural network approach for underreamed piles subjected to uplift forces in dry sand,” Neural Computing and Applications, vol. 31, no. 2, pp. 327–336, 2019. View at: Publisher Site  Google Scholar
 H. Moayedi and S. Hayati, “Applicability of a CPTbased neural network solution in predicting loadsettlement responses of bored pile,” International Journal of Geomechanics, vol. 18, no. 6, 2018. View at: Publisher Site  Google Scholar
 S. P. Chen, W. J. Yao, and S. Q. Zhu, “Nonlinear finite element analysis of superlong pile and soil interaction in soft soil,” Advanced Materials Research, vol. 201203, pp. 1601–1605, 2011. View at: Publisher Site  Google Scholar
 T. Jiang, Three Dimensional Nonlinear Finite Element Analysis of the PileGroup Effect for LargeDiameter and Superlong PileGroup Foundation, HoHai University, Nanjing, China, 2007, in Chinese.
 Ministry of Housing and UrbanRural Development, JGJ94, Technical Code for Building Pile Foundations, Ministry of Housing and UrbanRural Development, Beijing, China, 2008.
 Y. F. XIE and S. C. CHI, “Optimization method for irregular piled raft foundation on layered soil media,” Advances in Civil Engineering, vol. 2019, Article ID 5713492, 15 pages, 2019. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Yunfei Xie and Shichun Chi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.