Pile-Soil Stress Ratio and Settlement of Composite Foundation Bidirectionally Reinforced by Piles and Geosynthetics under Embankment Load

'e settlement calculation of composite foundation bidirectionally reinforced by piles and geosynthetics is always a difficult problem. 'e key to its accuracy lies in the determination of pile-soil stress ratio. Based on the theory of double parameters of the elastic foundation plate, the horizontal geosynthetics of composite foundation are regarded as the elastic thin plate, and the vertical piles and surrounding soil are regarded as a series of springs with different stiffness. 'e deflection equation of horizontal geosynthetics considering its bending and pulling action is obtained according to the static equilibrium conditions. 'e equation is solved by using Bessel function of complex variable, and the corresponding deflection function of horizontal geosynthetics is deduced. 'en, the calculation formula of pile-soil stress ratio and settlement of composite foundation is derived by considering the deformation coordination of pile and soil. 'e results of engineering case analysis show that the theoretical calculation results are in good agreement with the measured values, which indicates that the proposedmethod is feasible and the calculation accuracy is good. Finally, the influence of composite modulus of horizontal geosynthetics, tensile force of geosynthetics, and pile-soil stiffness ratio on pile-soil stress ratio and settlement is further analyzed. 'e results show that the pile-soil stress ratio increases with the increase of the composite modulus of the horizontal geosynthetics, the tensile force of geosynthetics, and the pile-soil stiffness ratio, and the settlement decreases with the increase of the composite modulus of the horizontal geosynthetics, the tensile force of geosynthetics, and the pile-soil stiffness ratio. When the flexural stiffness of the horizontal geosynthetics is small, the influence of the tensile action of the geosynthetics on the pile-soil stress ratio and settlement of the composite foundation cannot be ignored.


Introduction
In recent years, with the rise of highway and railway construction in China, the problems of soft soil foundation have become increasingly prominent, the bearing capacity of the foundation is insufficient, the settlement is too large, and uneven settlement is particularly serious. e composite foundation reinforced by piles and geosynthetics, a new type of soft foundation treatment consisting of vertical piles and horizontal geosynthetics, has a good effect on the treatment of the abovementioned foundation problems, and the foundation treatment method has simple construction and rapid progress. In recent years, it has been widely used in railway soft foundation treatment [1].
Research on the working mechanism of composite foundation has achieved many results so far [2][3][4]. Among them, the pile-soil stress ratio and settlement are important parameters reflecting the working behavior of composite foundation reinforced by piles and geosynthetics, and the research on its calculation method has also been a hot spot in this field. Among them, Yang et al. [5] uses the helix to simulate the deformation shape of the cushion, taking into account the net effect of the cushion; Chen et al. [6] considered the tensile effect of the geosynthetics, assumed that the shape of the cushion after deformation is a rotating paraboloid, and deduced the calculation formula of pile-soil stress ratio based on Winkler foundation; Abusharar et al. [7] assumed that the stress of soil between piles is uniformly distributed, and a simple method for calculating the stress ratio of pile to soil is put forward by using the arc assumption of cushion deformation. For a considerable part of the project, single-layer grille cannot meet the needs of the project, geogrid, or multilayer grille should be installed in the cushion. At the same time, reference [2] points out that when this kind of structure is different from the single-layer grille, the geosynthetics and the contained bulk fillers can form a "flexible raft" with strong bending resistance. Obviously, the above method based on the thin film theory is not suitable for geogrid or multilayer grille because it cannot take into account the bending effect. erefore, Rao and Zhao [8] assume that the geosynthetics is a thin plate, and the shape of the thin plate is a parabola in order to consider the "flexible raft effect." e formula for calculating the stress ratio of pile to soil is derived; Zheng et al. [9] think that this kind of cushion not only has a "flexible raft effect" but also has "pulling membrane effect," so it is assumed that the deformation of geosynthetics is in the form of heavy trigonometric function, and the geometric nonlinearity in the deformation process is considered. As a result, the pile-soil stress ratio and settlement calculation method considering both "flexible raft effect" and "pulling membrane effect" of cushion are obtained. Ma et al. [10,11] have done some work in the calculation of soft soil settlement. He et al. [12] made some attempts to chemically improve special soils. In addition, it can be seen from the above literature that the deformation shape of geosynthetics is the premise and key to calculating the pile-soil stress ratio and settlement, and the pile-soil stress ratio and settlement obtained by different deformation shapes are also different. erefore, the closer it is to the actual deformation of the geosynthetics, the more accurate the pile-soil stress ratio and settlement will be. However, the above literature all assume the deformation of the geosynthetics in advance, which cannot really reflect the stress state of the geosynthetics, which can easily lead to a large error between the calculated results and the measured results. For this reason, Tan et al. [13] derived the flexure function of horizontal stiffened body in the threedimensional case based on Winkler elastic foundation rectangular plate theory, but due to the complexity of rectangular plate boundary conditions, it is difficult to consider the influence of pile settlement on geosynthetics. After that, Zhao et al. [14] took the composite foundation within the influence area of a single pile as the research object, discretized the pile and soil into a series of springs with different stiffness, and with the aid of the circular plate theory of elastic foundation and the deformation coordination of cushion at the pile-soil interface, the deflection function of geosynthetics considering the common deformation of pile and soil is derived, and on this basis, the calculation method of pile-soil stress ratio is put forward. However, the above two analysis methods still only consider the "flexible raft effect" of geosynthetics, but not the "pulling membrane effect" of the reinforcement.
From the above analysis, we can know how to comprehensively consider the flexible raft effect, net effect, and pile-soil joint deformation of cushion without assuming the deformation of geosynthetics is the key to accurately solve the pile-soil stress ratio and settlement of composite foundation reinforced by piles and geosynthetics. erefore, on the basis of the above research work, based on the Filonenko-Borodich two-parameter foundation model [15], this paper comprehensively considers the "flexible raft effect" and net effect of the geosynthetics, considers the pilesoil deformation at the same time, deduces the analytical function of the flexural deformation of the geosynthetics according to the Bessel complex function, and then obtains the calculation expression of pile-soil stress ratio and settlement of composite foundation. In order to further improve the calculation method of pile-soil stress ratio and settlement of composite foundation.

Calculation Model
Take the typical unit within the influence range of the single pile as shown in Figure 1 for analysis. e diameter of the pile is d, the center distance is s d , the diameter of single pile reinforcement range is d e , square pile d e � 1.13s d , and piles are arranged as the shape of plum blossoms: d e � 1.05s d .
In order to further simplify the calculation, the following assumptions are made: (1) As shown in Figure 2, the geosynthetics within the range of single pile reinforcement is regarded as an elastic circular thin plate placed on the Filonenko-Borodich two-parameter foundation model [16], and among them, the tension T of the reinforcement can be obtained by the following formula: where ε g is the average strain of the reinforcement and E g is the tensile stiffness of the reinforcement. (2) It is assumed that the vertical supporting force P p of the pile to the thin plate is uniformly distributed, and according to references [6,13,14], the deformation of pile and soil between piles accords with Winkler foundation model. (3) e distribution of embankment load caused by differential settlement is not considered, that is, the embankment load is uniformly distributed [7,8,14].

Analysis of Geosynthetics
Let w 1 (r) represent the flexure function of the pile top reinforcement. According to the Filonenko-Borodich twoparameter foundation model, under the joint action of the embankment load q and the pile top reaction force p p , the control differential equation of the horizontal geosynthetics is as follows: where D is the bending stiffness of the thin plate.D � (Eδ 3 /12(1 − ] 2 )), where E is the elastic modulus of the plate, ] is Poisson's ratio of the plate, and δ is the thickness of the plate: . en, the homogeneous equation corresponding to equation (2) can be transformed into e general solution of the homogeneous equation corresponding to equation (2) is where J N and Y N are the first and second kinds of N-order Bessel functions and C 1 , C 2 , C 3 , and C 4 are undetermined constants.
Under the action of uniformly distributed load q and supporting force p p , the special solution of equation (2) is Combining formula (4), the flexural deformation function of the thin plate at the top of the pile is

Deformation Analysis of Thin Plate
Let w 2 (r) denote the flexural function of the thin plate at the top of the pile, under the action of uniformly distributed load q and foundation soil support, the governing differential equation of thin plate at the top of interpile soil is as follows: where k s is the spring coefficient of soil foundation between piles.
To simplify the calculation, e homogeneous equation corresponding to equation (6) can be transformed into the following form: It can be seen from formula (7) that the values of a 21 and a 22 are divided into real and imaginary numbers according to the positive and negative of T 2 − 4Dk s . e solution of equation (8) is different in the two cases, which are discussed as follows.
(1) When T 2 ≧ 4Dk s , order: e general solution of equation (8) is where N n is the second kind of N-order virtual variable Bessel function and C 5 , C 6 , C 7 , and C 8 are undetermined constants.
Under the action of uniformly distributed load q, the special solution of equation (6) is w 2 * � q/k s , and the combination formula (11) can be used to obtain the flexural deformation function of thin plate at the top of soil between piles at T 2 ≧ 4Dk s : When T 2 < 4Dk s , the general solution of equation (6) can be constructed as follows:

Advances in Civil Engineering
where u n and v n are the real parts of the first kind of Bessel function of order N and the real parts of Hankel function and f n and g n are the imaginary part of Bessel function of order N and the imaginary part of Hankel function, respectively. C 5 , C 6 , C 7 , and C 8 are undetermined constants.

Solution of Undetermined Coefficient
In order to determine the bending function of the thin plate, it is also necessary to compare the unknown parameters C 1 , C 2 , . . ., C 8 and p p in equations (5), (11), and (12), which are solved by boundary and continuous conditions.
As the turning angle of the thin plate at the center of the circular plate is 0, C 2 � C 3 � 0 can be obtained. At the pile-soil junction, the continuity conditions of thin plate deflection, rotation angle, radial moment, and shear force are as follows: In addition, according to hypothesis (2), the deflection w 1 (d/2) of the circular plate at the pile-soil interface is e value of pile top reaction p p is the sum of the upper embankment load and the load transferred from the cushion to the pile, that is, According to reference [2], the corner and shear boundary conditions at r � d e /2 of the element are as follows: At the same time, according to the bending functions (5), (11), and (12) of the thin plate, the expressions of the rotation angle, radial bending moment, and shear force of the thin plate can be obtained.
Within the range of the top of the pile, that is, Within the range of the top of the pile, that is,
us, the deflection expression of the plate is obtained.

Pile-Soil Stress Ratio and Settlement
From the above analysis, if the strain of the steel bar is known before the calculation, the tension of the steel bar can be calculated directly by using formula (1), and then the undetermined coefficient can be solved according to the steps in Section 5. However, in most projects, the deformation of reinforcement is not measured beforehand. In view of this situation [4] through field test and study, it is found that the relationship between average strain and settlement of reinforcement under embankment load is approximately as follows: where α and β are the relevant fitting parameters, respectively, which can be selected according to the settlement calculation position. According to the table listed in reference [4], 50 kN/m is the tensile stiffness of the grid used in the test. According to the above research results, the maximum deflection w 2 (d e /2) of geosynthetics is selected as the settlement S in formula (22). Combined with the method in Section 5, when the tension of the reinforcement is unknown, the steps of solving the undetermined coefficient are as follows: (1) First of all, assuming an initial tension value T 0 (T 0 can be a very small value), it is substituted into equations (2) and (6)  e flowchart of the above steps is shown in Figure 3 as follows.
After obtaining the flexure function and p p of the geosynthetics, it can be seen from the calculation model shown in Figure 2 that, according to the stress balance in the z direction, the average vertical stress p s at the top of the soil between piles can be obtained as follows: en, the pile-soil stress ratio n can be calculated according to the following formula [16]:

Advances in Civil Engineering
Reference [8] thinks that, in the case of large area overloading, the settlement at the top of the pile is relatively small, and the settlement of soil between piles can be regarded as the settlement of the foundation. (25)

Stiffness Coefficient of Pile k p .
In order to comprehensively consider the nonlinearity in the process of pile deformation, the deformation stiffness of pile is taken as the secant slope of the load test q-s curve of pile foundation.

Spring Coefficient of Surrounding Soil k s .
In order to comprehensively consider the nonlinearity in the process of soil deformation between piles, the coefficient of soil foundation between piles k s takes the secant slope on the load test q-s curve. If there is no measured data, it can also be taken in accordance with regional experience or determined according to the following formula [13]: where E s is the deformation modulus of soil, the weighted average of multilayer soil is depth, and H s is the thickness of soil layer.

Calculated ickness of Cushion δ.
From the point of view of the joint action of the geosynthetics and the bulk material pile, the thickness of the composite structure formed by the geosynthetics and its wrapped filler is taken as the thin plate thickness if it is a multilayer grille, that is, the distance between the top layer and the bottom grille. If it is a geotechnical cell, the thickness of the geotechnical cell can be taken directly.

Composite Elastic Modulus of Cushion E.
For the grid cushion, Zheng et al. [9] gave a method to determine its composite elastic modulus, that is, the weighted average value of the grid elastic modulus and the deformation modulus of the cushion filler. In the case of geotechnical cell, Yang et al. [17] gave the average elastic modulus of different types of geotechnical cell combined with various common fillers through stacked beam test.  [4], and the pile is powder injection pile, the diameter of the pile is 0.50 m, the center of the pile is 1.0 m, piles are arranged as the shape of plum blossoms and de � 1.05 m, and powder injection pile is laid with a double-layer geogrid with 50 kN/m tensile stiffness and a distance of 30 cm. e foundation is quaternary alluvial soft soil, the bedrock is mudstone, the upper filling load is 20 kN/m 3 , the filling height of roadbed is 10 m, and the crushed stone cushion is 25 kN/m 3 . According to the static load test, the stiffness coefficient of pile is k p � 2000 kN/ m. e thickness of the thin plate is 0.3 m, the composite elastic modulus E is 52 MPa, and Poisson's ratio v is 0.3. Other calculation parameters are shown in Table 1. e comparison between the calculated and measured values of the central settlement of the embankment and the pile-soil stress ratio under the geogrid geosynthetics is shown in Table 2.

Example and Parameter Analysis
e development trend of foundation settlement S of DK10 + 336 test section with the increase of embankment filling height H is shown in Figure 4.
As can be seen from Tables 2 and 3, the pile-soil stress ratio and settlement of double-layer geogrid-reinforced composite foundation obtained by this method are close to the measured values.

Example 2.
e test section of an expressway lying on soft soil foundation in Hunan is treated with geocell + mixing piles. e pile is arranged in the form of plum blossoms, and its diameter is 0.50 m with its spacing 1.2 m. e pile top is filled with a thick 30 cm sand cushion, and the center is equipped with a geotechnical cell with the thickness of 10 cm. e foundation is muddy clay, the load of the upper fill is 20 kN/m 3 , and the filling height of the roadbed in the test section is 4 m. e measured settlement value S is 5.3 cm, and the pile-soil stress ratio is n. e stiffness  coefficient of pile k p is 2355 kN/m. After treatment, the coefficient of soil foundation between piles k s is 1024 kN/m 3 , the thickness of thin plate δ is 0.10 m, the composite elastic modulus E is 40 MPa, and the composite Poisson's ratio v is 0.3. e pile-soil stress ratio and settlement are calculated as shown in Table 3.
It can be seen from Table 3 that the pile-soil stress calculated in this paper is closer to the measured value than the method in reference [14]. In addition, the settlement is also close to the measured value. is is due to the fact that compared with the method in reference [14], the tensile effect of the cell body is considered in this paper.

Parameter Analysis.
In order to further discuss and analyze the effects of geosynthetics composite elastic modulus, grille tension and pile-soil stiffness ratio on the settlement and pile-soil stress ratio of composite foundation reinforced by piles and geosynthetics, based on the parameters in example 2, the corresponding parameters are analyzed according to the abovementioned factors.

Influence of Pile-Soil Stiffness Ratio on Settlement and
Pile-Soil Stress Ratio. Without considering the influence of tension of geosynthetics, k p � 4K p /(πd 2 ) is introduced to characterize the comprehensive influence of pile    Advances in Civil Engineering 7   Advances in Civil Engineering deformation stiffness and replacement ratio, and other parameters remain unchanged. As shown in Figure 5, the influence of pile-soil stiffness ratio on settlement increases with the increase of elastic modulus E. At the same elastic modulus, the settlement of composite foundation decreases with the increase of pile-soil stiffness ratio. However, when the pile-soil stiffness ratio exceeds a certain value, continuing to increase the pile-soil stiffness ratio, it has little effect on reducing settlement.
As can be seen from Figure 6, the relationship between pile-soil stiffness ratio and pile-soil stress ratio is nonlinear under different elastic modulus, and the smaller the elastic modulus is, the more obvious the nonlinear phenomenon is. In addition, the larger the elastic modulus is, the more obvious the influence of pile-soil stiffness ratio on pile-soil stress ratio is.

Effect of Reinforcement Tension on Settlement and
Pile-Soil Stress Ratio. As can be seen from Figures 7 and 8, for the same elastic modulus E, the settlement decreases with the increase of reinforcement tension T, while the pile-soil stress ratio n increases with the increase of tension T.
In addition, the effect of Ton both pile-soil stress ratio and settlement decreases with the increase of E. It can be seen that when the elastic modulus of geosynthetics-reinforced cushion is small, using graded material and properly increasing the interface friction can effectively increase the pile-soil stress ratio and reduce the settlement.

Conclusions
Based on theory of the Filonenko-Borodich two-parameter elastic foundation model, the horizontal geosynthetics of composite foundation are regarded as the elastic thin plate, and the vertical piles and surrounding soil are regarded as a series of springs with different stiffness. e deflection equation of horizontal geosynthetics considering its bending and pulling action is obtained according to the static equilibrium conditions. e equation is solved by using Bessel function of complex variable, and the corresponding deflection function of horizontal geosynthetics is deduced. e calculation method for pile-soil stress ratio and settlement of composite foundation is derived by considering the deformation coordination of pile and soil.
(1) e settlement of composite foundation decreases with the increase of pile-soil stiffness ratio, geosynthetics tension, and composite elastic modulus of geosynthetics. e pile-soil stress ratio of composite foundation increases with the increase of pile-soil stiffness ratio, elastic modulus of geosynthetics, and its tension. When the bending stiffness of geosynthetics is small, the influence of tensile action of reinforcement on pile-soil stress ratio and settlement cannot be ignored. (2) e key to calculating with this method is to accurately measure the composite elastic modulus of geosynthetics-reinforced cushion, the strain of geosynthetics, and the stiffness of pile and soil. (3) Although this presented method does not rely on preassumed deformation when calculating the pilesoil stress ratio and settlement of composite foundation, the net effect of geosynthetics-reinforced cushion can be considered only when the strain of geosynthetics is known, which brings defects to the convenience of calculation. e proposed method still cannot consider the coupling relationship between deformation of geosynthetics-reinforced cushion and tension of geosynthetics, which can be carried out in the next work.

Abbreviations
T: e tension of geosynthetics ε g : e average strain of geosynthetics E g : e tensile stiffness of geosynthetics D: e bending stiffness of the thin plate

Conflicts of Interest
e authors declare that they have no conflicts of interest.