1. Introduction
As a new type of variable cross-section pile, the squeezed-branch pile was developed from the conventional circular pile. It has a central rounded shaft with at least one branch attached to it. The branch penetrates the relatively hard soil layer using a hydraulic squeezed machine. The soil surrounding each branch is applied to the pile’s top static load. Due to its benefits including a high bearing capacity and low settlement, the squeezed-branch pile is widely used in high-rise buildings. Although the squeezed-branch pile group foundation was constructed with equal stiffness, the differential settlements cause the superstructure to distort and crack, which has a significant impact on a building’s ability to be used normally [
1]. In view of these problems with the traditional design of pile group foundations, Franke et al. [
2] proposed the concept of the variable-stiffness design of pile groups. By adjusting the pile length, the foundation settlement tends to be uniform, reducing the internal force of the foundation and the secondary stress of the superstructure.
Numerous academics have undertaken in-depth research to construct the pile group foundation’s variable-stiffness and leveling design, mostly using the theoretical-analysis method, model-test method, and numerical simulation method. In terms of theoretical analysis, Leung et al. [
3,
4] established the shape function of pile-length optimization in the pile group foundation, optimized the 5 × 5 pile group foundation, and analyzed the overall stiffness and differential settlement of the pile group before and after optimization. Chan et al. [
5] used a genetic algorithm to optimize the vertical bearing capacity and the uneven settlement of the cap. Zhang et al. [
6] combined the pile-group-foundation-displacement function with the pile-length-optimization function, compared different pile-layout schemes, and studied the axial direction of the pile group foundation’s variation law of stiffness ratio and differential settlement ratio. Rasim et al. [
7] used the GWO algorithm to optimize the pile group foundation under vertical load. Han et al. [
8] carried out an improved two-way asymptotic structure optimization method, based on the continuum-topology-optimization theory, to produce the variable stiffness and leveling of a piled raft foundation using friction piles.
In terms of model tests, Zhang et al. [
9] analyzed the relationship between the plane layout, geometric size and settlement of piles by comparing models of variable-stiffness piled-raft foundations. Qian et al. [
10] compared seven sets of physical-variable stiffness and physical-simulation tests and found that changing the pile length has a more significant effect on the variable stiffness and leveling. In terms of numerical simulations, Dang et al. [
11] studied the influences of pile length, pile number, load distribution and raft thickness on the settlement and deformation of piled raft foundations.
In another study using numerical simulation methods, Sheil and McCab [
12] proposed a set of equations which can be used to predict the stiffness efficiency of the pile group based on a finite element analysis (FEA) parametric study of several variables that affect stiffness efficiency. Based on a parametric study of three-dimensional FEA simulations of jet-grouted rafts, Algin [
13] presented a multi-objective-optimization analysis using the surface response method to achieve the most economical design solution.
The majority of the pile group foundation optimizations discussed above use the methods of variable pile length, pile diameter, and variable pile spacing to modify the overall layout of the pile groups to achieve the goal of reducing the differential settlement of pile groups. However, this variable-stiffness-adjustment method is aimed entirely at the constant cross-section pile group foundation. For the squeezed-branch pile group foundation, the variable-stiffness-leveling method is still seldom researched.
Therefore, this paper first analyzes and simplifies the variable-stiffness principle of the pile group in order to address the issue of uneven settlement on the top of the squeezed-branch pile group. It develops an optimization model for the squeezed-branch pile group, with the number of branches, vertical bearing capacity, maximum settlement value and settlement difference as the objective function and constraint conditions. Secondly, the particle swarm algorithm (PSO) was improved to make it more suitable for solving the optimization problem of pile groups with squeezed branches and plates. Finally, the ABAQUS FEA software was redeveloped through the Python language [
14] and an automatic optimization platform for the squeezed-branch pile group was created, to achieve variable rigidity leveling of the squeezed-branch pile group foundation.
2. Optimization Model of Squeezed-Branch Pile Group
As a component of the superstructure–foundation system, the foundation of a high-rise building is constrained by the combined action of three aspects. Liu et al. [
15] presented the equation for the total balance as follows:
where
is the superstructure-stiffness matrix,
is the base-stiffness matrix,
is the support-stiffness matrix of the foundation soil,
is the displacement vector at the bottom of the foundation,
is the load vector of the superstructure, and
is the load vector of the foundation. When the pile group foundation is in the design stage, the stiffness matrices
and
of the superstructure and foundation, respectively, are determined, and the corresponding loads
and
are also determined accordingly. Therefore, only by adjusting the support-stiffness matrix of the foundation soil
can the displacement vector at the bottom of the foundation
be made uniform.
For conventional piles, pile-side resistance and pile-bottom resistance make up the majority of the support stiffness of foundation soil. Because the squeezed-branch pile adds the bearing plate as a member of the pile body, its foundation-soil-supporting stiffness matrix can be expressed as follows:
where
is the pile-bottom stiffness matrix,
is the pile-side stiffness matrix, and
is the branch stiffness matrix.
By adjusting the stiffness matrices of the pile side, pile bottom, and bearing plate, the purpose of reducing the differential settlement of the pile group can be achieved. However, if the three aspects are adjusted at the same point in the optimization process, there are many factors involved and optimization is not easy to achieve. For this reason, it is necessary to grasp the key part of the optimization problem. As the load borne by the branch of the squeezed-branch pile accounts for about 55–65% of the total load [
16], the stiffness matrix of the branch is selected as the key factor in the optimization problem, and the optimization problem is simplified as follows:
- (1)
The pile length L, pile diameter r and pile spacing S of the squeezed-branch piles are not changed during the optimization process.
- (2)
The diameter D of the branch, the height of the branch h, and the ratio t of the branch diameter to the pile diameter are fixed.
- (3)
The spacing of the branch set on the pile body of the squeezed-branch plate meets the specification requirements.
After simplifying the overall balance equation of the squeezed-branch pile group foundation, the soil-support-stiffness matrix is only related to the branch stiffness matrix . For a squeezed-branch pile group, the stiffness matrix of the branch is proportional to the number of branches. Therefore, the stiffness matrix of the branch can be changed by adjusting the number of branches to adjust the settlement of the pile group foundation, in order to achieve the purpose of reducing the differential settlement of pile groups.
Through the principle of variable stiffness adjustment, the optimization problem of a squeezed-branch pile group is transformed into the distribution problem of the pile branch, and the optimization model with the number of branches on each pile as the design variable is established. The total volume of the pile foundation is typically regarded as the objective function of optimization in pile group optimization, because it allows for cost management and straightforward calculation. However, when the total volume is selected as the optimization objective function of a squeezed-branch pile group, a change in the number of pile group branches has little effect on the overall volume; this may lead to difficulties in the optimization process. For this reason, the number of total branches of the pile group is selected as the objective function of the optimization problem. This not only reflects the difficulty of construction, but also makes it possible to comprehensively consider the influence of factors such as time and cost.
A method for combining and reoptimizing pile groups is proposed, due to the significant amount of calculation required for individual optimization of each pile, which is not conducive to actual construction. The pile group foundation is modeled and calculated in ABAQUS FEA software according to the original design scheme, and the axial force
of each pile top is obtained through post-processing operations. As shown in
Figure 1, the maximum value
and the minimum value
of the pile-top axial force are selected as two special values, and the average value
of the pile-top axial force is calculated. Next, the three special values divide the pile-top axial force into two regions. After averaging the special values
times, a total of
special values are generated. The axial force of the pile top is divided into
regions, and the whole pile group is also divided into
groups. The number of pile groups
can be determined according to the number of piles contained in the piled raft foundation and the uneven degree of foundation stress.
The number of squeezed-branch pile-bearing plates in the same group is consistent and replaced by a design variable. After grouping, the objective function of the squeezed-branch pile group optimization problem becomes:
where
is the solution of the squeezed-branch pile group optimization problem,
is the design variable for the
group of piles,
is the objective function,
is the number of pile groups,
is the number of piles contained in group
, and
is the number of branch sets for group
piles.
For the optimization problem of the squeezed-branch pile group, the objective function should be solved under constraint conditions. After considering the vertical bearing capacity, maximum settlement value, and settlement difference of the pile group, the constraint conditions are established as follows:
where
is the vertical-load value of each single pile in the
group of piles,
is the characteristic value of the vertical bearing capacity of a single pile in the pile group,
is the maximum settlement of the raft calculated by the FEA software, and
is the maximum allowable settlement value of the piled raft foundation, which can be determined according to the design level of the building in conjunction with the code.
is the minimum settlement of the raft, while
represents the maximum allowable settlement difference of the piled raft foundation after optimization, which may be determined by the specific project.
3. Improved Particle Swarm Algorithm
To produce the variable stiffness leveling of the pile foundation, the particle swarm algorithm is used. The PSO is a cluster intelligence optimization algorithm for nonlinear function optimization, proposed by Kennedy, J. and Eberhart, R.C. in 1995, which is inspired by the behavioral characteristics of a flock of birds during foraging. The PSO forms an initial population by randomly generating particles of certain sizes and finding the corresponding fitness value according to the objective function of the optimization problem. The initial population flies within the feasible solution to find the maximum or minimum fitness value; the direction and distance of the particle flight are determined by the velocity.
Assuming that the total number of particles in d-dimensional space is N, the displacement of particle
in dimension
at time
is
, and the fitness value of particle
is calculated as follows:
where
is the fitness value of particle
,
is a relatively large number to ensure that the fitness value is positive,
is the degree of violation of each constraint condition by particle
,
is the value of constraint condition
in particle
,
is the constraint value for each constraint condition.
Using the fitness function
to calculate particle-fitness value, the position where the particle is currently searching for the largest fitness value is noted as the optimal position
. The position where all particles in the whole swarm are currently searching for the largest fitness value is denoted as
. The velocity of particle
is
. The velocity and displacement of particle
in dimension
at time
(
and
) are simply derived from Equations (7) and (8).
where
represents the inertia weight coefficient,
and
are the cognitive and social acceleration constants, and
and
are two random numbers that form a uniform distribution [0, 1].
For solving the problem of minimizing the objective function, the PSO updates the
and
using Equations (9) and (10).
To solve the pile group foundation optimization problem using a discretization approach, Equation (8) needs to be improved. Particle-position-update probability is calculated by mapping the particle velocity to [−1, 1] via the Tanh function. The value
denotes the probability that the position update velocity of particle
in dimension d takes 1, −1, 0; the relationship between
and
is shown in Equation (11). Given the large absolute value of
, there is a high probability that
will take on a value of 1 or −1. Meanwhile, the particle-velocity formula remains unchanged; the position is updated according to Equation (12).
where rand() is a function that generates a random number between 0 and 1.
At the initial stage of iteration, the PSO needs to enhance the particles’ global search capability to guarantee that the particles can search the entire solution space at this stage. However, in the late iteration phase, the local search ability of the particles needs to be consolidated, increasing the probability of the algorithm searching for the global optimal solution by boosting the particles’ ability to search in the local solution space.
The particles’ global and local search abilities are related to the inertia weight coefficient , cognitive-acceleration constants and social-acceleration constants in Equation (13). When and take large values, it is beneficial to enhance the particles’ global search ability, while when and take small values, it is beneficial to enhance the particles’ local search ability.
To improve the optimized performance of the PSO, an adaptive adjustment strategy was used to dynamically correct the inertia weight coefficient and acceleration constants
and
, which was calculated as follows:
where
is the initial inertia weight coefficient,
is the number of current iterations of the particle swarm, and
is the total number of iterations specified.
Finally, the PSO is stopped when a specified number of generations is reached and the optimal position of the particle swarm is solved for the optimization problem.
The basic procedure for the overall optimization of a squeezed-branch pile group is as shown in
Figure 2.
FEA is performed on the original design plan, and the axial force of each pile is extracted.
An appropriate number of groups is selected; the pile group foundation is divided into regions.
The single-pile form of squeezed and expanded branch plates that can be selected in each group is determined.
Initial populations are created and fitness value is calculated.
The optimal individual particle position and the optimal particle-population position are updated.
The velocity and position of the particle are updated, and particle quality is evaluated.
Steps e to f are repeated until the specified number of evolutionary generations is reached, and is found in the population of particles corresponding to the optimal design of the squeezed-branch pile group.
5. Optimization Results
For the optimization of the engineering example, the results of the area division were used to arrange the pile groups. The initial population size was set to 50, and the number of evolutionary generations was 50. The inertia weight coefficient
was set to 0.7, the cognitive-acceleration constant
was set to 1.4, and the social-acceleration constant
was set to 0.6. The squeezed-branch pile group automatic optimization platform was used for the calculation, and 10 independent optimizations were taken. The optimized arrangement of the pile groups is shown in
Figure 6, and the distribution of the pile branches is shown in
Table 4.
The optimized pile group foundation increased the number of branches at the center pile and reduced the number of branches at the edge and corner piles. This pile arrangement conformed to the principle of variable stiffness around the weakening of the strengthening center. Reducing the number of branches from 324 to 142 and reducing the number of branches by 56% greatly reduced the cost and accelerated the speed of construction, making the pile group arrangement more reasonable and greatly improving its economic efficiency.
Figure 7 shows the raft settlement before and after the optimization of the pile group foundation. The pile group of the initial design has a large center and a small dish-shaped settlement phenomenon, while the optimized pile group raft weakens this phenomenon. The results before and after the optimization are plotted in
Table 5. It can be concluded from the table that after meeting the maximum settlement value, after optimization, the differential settlement of the pile group raft was reduced from the original 12.5 mm to 7.6 mm, which is a 39% reduction compared with the original design. The maximum axial force of the pile top was reduced from 2377.4 kN to 1569.3 kN, the minimum value of the pile-top axial force was almost unchanged, and the difference value of the pile-top axial force was reduced from 1312.8 kN to 517.1 kN, a 61% reduction.
The results show that the optimized pile group using PSO can greatly reduce the differential settlement of the raft by adjusting the distribution method of the pile-group-bearing plate under the conditions of satisfying the bearing capacity, settlement and other constraints. The maximum value of the pile-top axial force makes the settlement and force of this pile group foundation more uniform, thereby weakening the additional bending moment and shear force generated by the foundation and superstructure, and reducing the risk of cracking or tilting and collapsing of the overall structure. The safety reserve is improved, making it safer during use.
Figure 8 shows the axial force of the piles before and after the optimization of the design. In the initial design scheme, the maximum axial forces of the pile group at
Section 1 and
Section 2 were 1959.2 kN and 2377.4 kN, respectively, and the minimum forces were 1064.6 kN and 1889.6 kN. The axial-force differences between the two sections were 887.6 kN and 487.8 kN, respectively. The maximum and minimum axial forces of the optimized pile group at
Section 1 were 1437.7 kN and 1052.2 kN, respectively. The maximum and minimum values of the axial force at
Section 2 were 1569.3 kN and 1435.8 kN, respectively. The axial-force difference between the two sections decreased from 887.6 kN and 487.8 kN to 382.5 kN and 133.5 kN, respectively. The results show that the settlement curve and axial-force curve tended to be gentle due to the optimization of the extruded and expanded piles, that the dish-shaped distribution of the settlement and axial force of the pile group foundation was weakened, and that the additional shear force and bending moment of the pile group foundation were reduced; these changes were more conducive to the safe use of the superstructure.