Numerical Analysis of Bearing Behavior of the Prebored Precast Pile with an Enlarged Base

Prebored precast pile with an enlarged base (PPEB pile) is a new type of green and environmental protection pile foundation developed in China in recent years, which has complex bearing characteristics and many inﬂuencing factors. Based on the static load tests and key parameters’ tests in deep soft soil in Shanghai, a three-dimensional numerical analysis model was established using ABAQUS ﬁnite element software. The transfer law of load among the precast pile, cement soil, and soil around the pile and the action mechanism of the enlarged base were analyzed emphatically, and a sensitivity analysis of the main factors aﬀecting the bearing performance was carried out. The calculation results show that the existence of the enlarged base can greatly improve the compressive bearing capacity, increasing the diameter and height of the enlarged base is beneﬁcial to the bearing capacity, and the inﬂuence of the diameter expansion ratio is more eﬀective. With the increase of the proportion of nodular piles, the ultimate bearing capacity increases slightly, but the deformation increases obviously. Under the condition of cement soil of the test piles, the spacing of the neighboring nodules of nodular piles has no obvious eﬀect on the bearing capacity, and the 1 m spacing commonly used in engineering applications can be optimized. The increase of cement soil thickness is beneﬁcial to the improvement of pile bearing capacity, but the eﬃciency is low. Finally, some improvement measures for the construction technology of the PPEB pile were put forward.


Introduction
Prebored precast pile with enlarged base (PPEB pile) is a new type of green and environmental protection pile foundation developed in China in recent years [1]. is PPEB pile requires a special installation process [2], as shown in Figure 1. First, a special mixing drill is used to drill, stir, and grout the soil to form a borehole with an enlarged base filled with cement soil slurry. en, the precast concrete pile is inserted into the cement soil slurry by its self-weight. After cement soil solidification, the pile structure is composed of rigid pile body wrapped by cement soil. e properties of the cement grout used in the pile shaft and the enlarged base are different. Generally, the premade base-forming cement grout with a water/cement (w/c) ratio of 0.6 is injected into the enlarged base, and side-forming cement grout with a w/c ratio of 1.0 is injected into the borehole. erefore, the pile body will be in the form of cement soil with low strength around the precast pile and cement soil with high strength at the enlarged base. e inner concrete precast piles can be made up of prestressed high-strength concrete (PHC) piles, PHC nodular piles, prestressed reinforced high-strength concrete (PRHC) piles, and so on.
Due to the special construction technology and complex material composition and structure, the bearing characteristics of the PPEB pile are complex and the influencing factors are numerous. On the basis of static load test and model test, many researchers [2][3][4][5][6] preliminarily discussed the load transfer law of the pile under vertical compression and verified its reliable bearing capacity. It should be emphasized that the number of on-site test piles or models is limited, which makes the analysis results have some limitations. Numerical simulation is an effective tool for analyzing complex engineering problems, especially the pile-soil interaction. rough numerical analysis, Ye et al. [7][8][9] discussed load transfer law and the effect of reinforced core on the performance and failure behavior of stiffened deep cement mixing piles. Based on the engineering application of PPEB pile technology in Zhejiang Province, China, Yang et al. [10], Zhou et al. [11], and Xie et al. [12] used the numerical simulation method to analyze the influence of structure size of precast pile, soil properties, cement soil characteristics, and other factors on bearing capacity, which enriched the understanding of the characteristics of PPEB piles. However, the previous numerical research lacks the analysis of the internal stress behavior of pile foundation, which is essential to understand the law of load transfer. rough field tests in Shanghai, the authors preliminarily discussed the bearing behavior of PPEB piles in deep soft soil stratum [2]. In this paper, based on the above static load tests and key parameter tests, the load transfer law between precast pile, cement soil, and soil around pile and the action mechanism of enlarged base were further analyzed by numerical simulation. en, a sensitivity analysis of the main factors affecting the compressive bearing performance was carried out.

Field Test
e field tests in this paper were carried out in the intersection of Jiang Yang South Road and Hong Wan Road of Shanghai. In general, there are eight soil layers within the upper 70 m of soil below the ground surface. e basic physical and mechanical properties of the soils are shown in Table 1. e effective cohesion and friction angle were determined from consolidated undrained tests. e undrained shear strength was obtained from in situ vane shear tests and empirical values. ree test piles, designated piles TP1, TP2, and TP3, were installed at the site. e length and shaft diameter of the test piles were 55 m and 750 mm, respectively. e inner precast concrete piles were assembled by a 40 m long, 600 mm diameter PHC pile with 110 mm wall thickness in the upper part and a 15 m long, 650 mm (500 mm) diameter nodular pile (i.e., 500 mm pile shaft diameter, which increases to 650 mm at the nodules) with 100 mm wall thickness in the lower part. e enlarged bases were 2750 mm long with a diameter of 1200 mm. To ensure the material strength of the precast piles, the test piles adopted C100 high-performance concrete with 100 mm cube compressive strength of 120.9 MPa and elastic modulus of 54.7 GPa.
Static load tests were carried out 43-45 days after pile installation. e maintained load method was adopted in accordance with the Chinese code JGJ106-2014 [13]. e test apparatus mainly consisted of a counterforce device using concrete blocks, a loading device using a hydraulic Jack, and a measuring system using force and displacement sensors. e test piles were instrumented with vibrating-wire strain gauges to obtain the axial force of the pile.

Advances in Civil Engineering
Obvious inflections can be observed on the load-displacement curves of TP1 and TP2. e displacements were small during the initial stages of loading and increased dramatically when reaching the ultimate bearing capacity of 8800 kN. As the pile head of TP3 was inclined during the test, the load test on TP3 was terminated at the applied load of 8000 kN. e load-displacement curve of TP3 had no obvious inflection. e displacements and applied loads of pile heads at several critical stages are summarized in Table 2. Detailed test results were described in a previous paper [2].

Numerical Modeling.
ree-dimensional finite element analysis (FEA) method incorporated in the software ABA-QUS was used to model the load of TP1. As the vertical loading of a single pile was axisymmetric, half of the model was taken for analysis to simplify the calculation. To minimize the influence of the FEA model boundary, the radius of soil around the pile was 20 m, a width more than 20 times the pile diameter in the plane, and the depth was 80 m, as shown in Figure 2. e distance from the pile toe to the bottom of the entire soil model was more than 20 times the pile diameter to ensure that the boundary effect was negligible. e bottom of the soil model was fixed in all three directions, and the other sides were constrained to allow vertical displacement only. e linear elastic model was used for precast pile and the Mohr-Coulomb elastic-plastic model was used for cement soil and soil around the PPEB pile. e initial ground stress of the soil was considered in the model, while the displacement control method was used to apply the vertical load. e parameters of the soil were selected according to Table 1. According to the engineering experience in Shanghai [14] and the results of trial calculation, the elastic modulus (E) of soil can be taken as 5 times the compression modulus (E s ). rough the back analysis of the results of the static load test [2], the elastic modulus of the precast pile was 60 GPa and Poisson's ratio was 0.15.
Cement soil is an important part of PPEB pile. e characteristics and parameters of cement soil, precast pilecement soil interface, and cement soil-soil interface are essential to obtain the correct numerical analysis results. e above parameters related to cement soil would be obtained through comprehensive analysis of field tests and laboratory tests.

Field and Laboratory Tests of Cement Soil.
Consistent with the requirements of the pile hole in the field test, the cement soil mixing pile K1 without precast pile was constructed. It could be considered that the cement soil property of K1 was the same as that of the test piles. erefore, the core samples were obtained from the K1 for testing. After 48 days of curing, the unconfined compressive strength (UCS) of cement soil was 1.7∼4.0 MPa, the cohesion was 342∼731 kPa, and the internal friction angle was 37∼45°. e deformation parameters of cement soil are usually measured by the deformation modulus (E 50 is the secant modulus of cement soil at 50% compressive strength) [15]. Based on the stress-strain relationship of unconfined compression specimens, the deformation modulus E 50 of cement soil was estimated to be 102∼222 MPa.
Due to the difficulty in drilling construction, the cement soil at the enlarged base could not be sampled and tested. erefore, according to the composition of the actual enlarged base, the laboratory proportioning tests of cement soil were carried out. e typical silty clay and silty sand in Shanghai were selected as the soil materials, and P·O 42.5 Portland cement was selected as cementitious material. In the actual construction, the amount of cement grout at the pile toe was 100% of the enlarged base volume, and then the drill pipe was stirred up and down evenly for 1 to 5 times. After the completion of grouting, the volume ratio (cement grout/slurry) was about 1.0 to 1.5. In the laboratory test, the specific gravity of the slurry was 1.5, while the water-cement ratio of cement grout was 0.6. e volume ratio (cement grout/slurry) was set to 1, 1.5, and 2 to prepare the cement soil slurry of enlarged base. e rough samples were made by the hanging bag method (see Figure 3) and then cut into cylinders with a height of 100 mm and a diameter of 50 mm.
After 28 days of curing in the standard curing room, the unconfined compressive strength test and shear test were carried out with the loading speed of 1 mm/min. e laboratory test results are shown in Table 3. According to the field and laboratory test results, in the FEA model, the values of cement soil parameters are shown in Table 4. e Poisson's ratios were assumed upon the empirical values.

Interface Definition.
In this FEA model, three interface elements were defined: precast pile-cement soil interface, precast pile-soil interface, and cement soil-soil interface. e Table 1: Summary of soil properties.

Soil layer
Depth contact behavior of all three interfaces was defined using an isotropic Coulomb friction model, which allowed for slippage and separation between the master and slave surfaces. e feasibility and accuracy of using the Coulomb friction model to simulate the pile-soil interaction had been well verified by many researchers [16][17][18]. eir ultimate friction resistance is expressed as where μ is the friction coefficient of the interface, σ y is the normal stress, and τ max is the control value of ultimate friction, which was obtained by back analysis of the axial force of pile shaft in field test (as shown in Table 1).
To understand the characteristics of concrete-cement soil interface more clearly, the model tests were carried out using a cylindrical interface model device, as shown in Figure 4. e strength grade of concrete material was C40, and the volume ratio of cement grout and slurry was set to 0.15 and 0.6 in two groups. e test results showed that when the UCS of cement soil was 522∼4572 kPa, the interfacial shear strength increased with the increase of UCS; the value was 0.167 to 0.224 (mean value 0.193). e UCS of cement soil of pile shaft was conservatively taken as 1.0∼2.0 MPa. According to the results of this interface model test, it was conservatively estimated that the interface friction force provided by concrete-cement soil was 170∼340 kPa. However, the ultimate local shaft resistance of cement soil-soil in PPEB pile in Shanghai was about 100 kPa [2]. erefore, it could be seen that the friction performance of precast pile-cement soil interface was obviously better than that of cement soil-soil interface. erefore, the friction coefficient of the precast pile-cemented soil interface was set to 0.80 to ensure high friction performance between precast pile and cement soil. e friction coefficient (μ) of the pile-soil interface can be calculated by equations (2) and (3) [19]:

Advances in Civil Engineering
where φ is the friction angle of soil and δ is the internal friction angle of the pile-soil interface. When the friction angle φ is equal to 15°to 30°, the friction coefficient μ is equal to 0.23 to 0.35 from (2) and (3). e internal friction angle of the pile-soil interface can be also estimated according to (0.7∼1.0) φ [20]. When φ � 15∼30°, δ � 0.8φ, the friction coefficient μ is equal to 0.21 to 0.45. Considering the field geological conditions, the friction coefficient between the precast pile and soil in the model was taken as 0.35. Part of the cement soil slurry infiltrates into the soil around the PPEB pile, so the friction performance of the cement soil-soil interface is generally better than that of the concrete-soil interface. Combined with the results of the back analysis of the field test, the friction coefficient of the cement soil-soil interface was set to 0.53.

Validation of the Numerical Results.
e load-displacement curve and axial force distribution of the test pile calculated by ABAQUS program are shown in Figure 5. e FEA modeling well simulated the bearing and deformation behavior of the field test pile, which showed that the established three-dimensional model and the selected parameters were reasonable. erefore, based on this three-dimensional numerical model, it was reliable to analyze the main factors affecting the compressive bearing performance.

Analysis of Compressive Bearing
Characteristics. e development of the plastic zone of the PPEB pile-soil system during loading is shown in Figure 6. AC YIELD is the yield mark [21], which is equal to 1 when yielding occurs and 0 in the case of unyielding. Since the value of the node is determined by interpolation, it may be more than 1 or less than 0 when displayed, but it does not affect the reliability of the calculation results. e representative loading conditions of the test pile are also shown in Figure 6, in which Q 1 is the elastic limit load, Q 3 is the ultimate load, Q 2 is the intermediate transition load, and Q 4 is the maximum load.
Before the initial loading stage Q 1 , the shaft resistance of PPEB pile in the upper soil was gradually mobilized; at the same time, the toe resistance had also played a role and increased with the loading. After loading Q 2 , the soil at the outer edge of the enlarged base was the first to enter the plastic state; there was a trend of separation between the shoulder of the enlarged base and the soil, so the stress relaxation of the soil led to yield. As the loading continued, the shaft resistance reached the limit, and more loads would be borne by the soil around the enlarged base. e plastic zone of the soil at the outer edge of the enlarged base gradually expanded and ran through the bottom of the enlarged base. When loaded to the ultimate load Q 3 , the soil stress relaxation zone whose length was about 1 time the diameter of the enlarged base was formed above the variable cross section of the enlarged base. After the ultimate load Q 3 , the continued loading was mainly borne by the soil at the end of the pile, and the plastic zone developed to the upper part of soil. e cement soil under the last nodule of the nodular pile also formed a yield zone due to the increase of stress and gradually developed from the bottom to the top of the enlarged base. Finally, the PPEB pile reached the maximum load Q 4 . During the whole loading process, there was no plastic zone in the cement soil around the pile shaft, which indicated that the strength of cement soil could ensure the bond strength of the interface between the precast pile and the cement soil.

Stress Analysis of the Pile Shaft.
e comparison between the total axial force of PPEB pile and the axial force of the precast pile under representative loading conditions is shown in Figure 7. For the nonenlarged base part, the stiffness of the cement soil around the pile was lower, so the axial force of the PPEB pile was mainly controlled by the precast pile. As the property of the cement soil at the enlarged base was much better than that around the precast pile, the load sharing of the enlarged base could not be ignored. △Q in Figure 7 is the load shared by the cement soil in the PPEB pile. e axial force of the cement soil at the enlarged base increased sharply with the depth and reached the maximum below the last nodule. On the contrary, the axial force of the precast pile decreased sharply within the range of the enlarged base. is reflected the role of the enlarged base. e force was transferred to the cement soil of the enlarged base through the nodules and then to the soil around the enlarged base. It was the reaction of a larger range of soils that enhanced the end bearing capacity of the PPEB pile. Under the ultimate load Q 3 , the enlarged base resistance (axial force of 52 m section) accounted for about 26% of the applied load, which was close to the field test results [5]. Under the loading of Q 1 , Q 2 , Q 3 , and Q 4 , the toe resistance (axial force of 55 m section) was 830, 965, 1415, and 2765 kN, respectively, accounting for 14%, 13%, 16%, and 27% of the total load. e stress of the cement soil was much less than that of the precast pile, so that the load was first transferred from the precast pile to the cement soil and then to the soil around the pile. e vertical stress of cement soil in the upper pipe pile increased gradually along the depth. e vertical stress distribution of the lower nodular pile section was more complex. In the process of loading, each nodule was squeezed with its lower cement soil and had a tension trend with the upper cement soil. erefore, within the range of each nodule, the vertical stress of the cement soil under the nodule was the largest, while the vertical stress of the cement soil above the next nodule was the smallest. e vertical stress of the cement soil of the enlarged base was also clearly different above and below each nodule, as shown in Figure 8.    Advances in Civil Engineering under the last nodule was the largest and developed upward. After the elastic load Q 1 , the tensile stress zone occurred above the last nodule. With the ultimate load Q3 applied, the tensile stress zone appeared above the three nodules in the enlarged base. After Q 3 , the extrusion of nodules was strengthened. However, the extrusion deformation was restricted by the displacement of the pile toe; the tensile stress zone was gradually closed from bottom to top. e cement soil of the enlarged base was also subjected to a higher shear stress. Figure 9 shows that, under the ultimate load Q 3 , the shear stress (1260 kPa) of the cement soil under the last nodule was obviously higher than that of the soil at the end of the pile (265 kPa). erefore, when the strength of the cement soil was not enough, the compression failure or shear failure of the cement soil might occur before the soil failure at the end of the pile.

Parameter Analysis
Horizontally, the structure of the PPEB pile is nonuniform, and longitudinally the size of the pile shaft is different (with enlarged base), as shown in Figure 10, so there are many Precast pile PPEB pile

e Height of the Enlarged Base.
e setting of the enlarged base is an important feature that distinguishes the PPEB piles from other bored precast piles and concretecored deep cement mixing piles. e structure of the enlarged base has a great influence on the bearing performance.
In the modeling calculation, the height of the enlarged base (L b ) was taken as 0, 1.0, 2.0, 2.7, 4, and 5 m, respectively (i.e., 0D, 1.3D, 2.7D, 3.6D, 5.3D, and 6.7D; D is the diameter of borehole), including 0, 1, 2, 3, 4, and 5 nodules, respectively. e calculation results are shown in Figure 11. Combined with the characteristics of Q-s curve and the development of soil plastic zone, the ultimate bearing capacity was defined, as shown in Figure 11(a). Comparative analysis of the ultimate bearing capacity of PPEB piles under different enlarged base heights is shown in Figure 11(b). Based on the case without enlarged base (L b � 0), when L b � 1.3D, 2.7D, 3.6D, 5.3D, 6.7D, the ultimate bearing capacity increased by 11.7%, 13.0%, 13.9%, 15.5%, and 16.7%, respectively. It could be shown that the existence of the enlarged base (i.e., including one nodule) could dramatically improve the bearing capacity, but with the increase of the height of the enlarged base, the increase of bearing capacity was limited. erefore, in actual projects, the height of the enlarged base is 3 times the borehole diameter, which contains 2∼3 nodules, and the force of the enlarged base can be guaranteed.

e Diameter Expansion Ratio of the Enlarged Base.
e diameter expansion ratio (λ) was defined as the ratio of the enlarged base diameter to the borehole diameter. In the modeling calculation, the diameter expansion ratio λ was set to 1.0, 1.3, 1.6, and 2.0, and the end bearing layer was set to medium-dense silt mixed silty clay (⑦) and dense silty sand (⑧). e numerical results under different conditions of diameter expansion ratio and pile end bearing layer are shown in Figure 12. e better the end bearing layer was, the greater the ultimate bearing capacity was. Similarly, the larger the diameter expansion ratio was, the greater the ultimate bearing capacity was; in addition, the proportion of toe resistance increased linearly. e choice of expansion ratio is more controlled by whether the strength of the enlarged base meets the requirements of punching or compression resistance and the limitation of the construction equipment capacity. In practical engineering, the diameter expansion ratio is generally controlled at 1.6.

Proportion of the Nodular Pile.
PHC nodular pile is generally used in the lower part of PPEB pile, and the nodules can effectively enhance the bonding force between the precast pile and cement soil. Especially at the end of the pile shaft, the PHC nodular pile and the enlarged base form a whole and bear the load together. In the numerical model, the length of the nodular pile segment was set to 0, 4, 15, 28, 40, and 55 m from the pile toe, accounting for 0%, 7%, 27%, 50%, 73%, and 100% of the total pile length, respectively.
With the proportion of nodular piles increased from 7% to 100%, the pile head displacement increased significantly while the ultimate bearing capacity increased only by about 3.5%, as shown in Figure 13. Since the stiffness of the nodular piles was smaller than that of PHC piles with the same diameter, the deformation increased with the increase of the proportion of nodular piles under the same loads. In particular, the PHC nodular pile used in the upper part of the pile shaft (such as 100%) was subject to a large load, which would share a larger pile deformation. On the other hand, the larger pile head displacement led to the larger pile-soil relative displacement, which improved shaft resistance to a certain extent. erefore, in projects with high requirements for deformation, it is a more reasonable form to set up 1 ∼2 section nodular piles in the lower part of the PPEB pile.

Spacing of the Neighboring Nodules.
e spacing of the neighboring nodules (L n ) of nodular piles was set to 1, 2, 3, and 4 m for analysis. Under the condition of cement soil (UCS � 1.7∼4.0 MPa) of the field test piles, the load-displacement curves were almost the same when using the above different spacing, as shown in Figure 14. e curves showed that the change of spacing of the neighboring nodules from 1 m to 4 m had no obvious effect on the bearing capacity. ere might be two main reasons for this result. First, the cement soil around the precast pile had good properties, which could ensure the high bond strength of the interface between the precast pile and the cement soil; high interface strength resisted the relative displacement between nodular pile and cement soil, thus weakening the influence of nodules and spacing of the Note: the solid is marked as ultimate load. neighboring nodules on the bearing performance. Second, in the numerical model, the nodular piles were in the lower part and accounted for only 27% of the pile length, so that the overall influence of the nodular pile segment was not significant. Hence, when the property of cement soil around the precast pile is good and the proportion of nodular piles is relatively small, the 1 m spacing commonly used in the project can be optimized.

ickness of Cement
Soil. e thickness of cement soil in PPEB pile varies with the diameter of the borehole. In the FEA model, the 650 mm (500 mm) diameter nodular pile was used in the lower part. When the diameter of cement soil borehole (D) was set to 750, 850, 950, and 1050 mm, the thickness of cement soil was 125, 175, 225, and 275 mm, respectively, in nonnodules.
Under the condition that the size of the enlarged base was constant, the bearing capacity increased nonlinearly with the increase of the thickness of cement soil, as shown in Figure 15(a). e efficiency analysis was based on the borehole diameter of 750 mm of the test pile, as shown in Figure 15(b). Although the increase of shaft resistance was similar to the development of bearing capacity, the increase rate of shaft resistance per unit area continued to decrease. It can be inferred that, for deep soft soil areas, the efficiency of improving the bearing capacity by increasing the thickness of cement soil is low. Advances in Civil Engineering 11

Conclusions
Based on the static load test and key parameters test in deep soft soil in Shanghai, a three-dimensional FEA model was established using ABAQUS. en, the bearing characteristics of the PPEB pile were analyzed. e main conclusions are as follows: (1) e Mohr-Coulomb elastic-plastic model was adopted for the cement soil around the pile shaft and at the enlarged base, and the elastic modulus was 150 MPa and 1500 MPa, respectively. e friction coefficients of precast pile-cement soil interface and cement soil-soil interface were 0.80 and 0.53, respectively. e numerical results well simulated the bearing and deformation behavior of the field test pile.
(2) At the initial stage of loading, the load had been transferred to the pile toe, and with the increase of the applied load, the shaft resistance was gradually mobilized from top to bottom. In the nonenlarged base part, the axial force of the PPEB pile was mainly controlled by the precast pile. In the enlarged base part, the load was transferred to the cement soil of the enlarged base through nodules and then to the soil around the enlarged base. e cement soil of the enlarged base was subject to higher vertical compressive stress and shear stress. (3) e existence of an enlarged base could greatly improve the compressive bearing capacity of the PPEB pile. Increasing the diameter and height of the enlarged base was beneficial to the bearing capacity, and the effect of diameter expansion ratio was more effective. In practical engineering, it is reliable to control the enlarged base diameter at 1.6 times the borehole diameter and adopt the enlarged base height of 3 times the borehole diameter. (4) With the increase of the proportion of nodular piles, the ultimate bearing capacity increased slightly, and the deformation increased obviously. Setting up 1 ∼2 section nodular piles in the lower part of the PPEB pile is a more reasonable measure for projects with higher deformation requirements. When the property of cement soil around the precast pile is good and the proportion of nodular piles is relatively small, the 1 m spacing commonly used in projects can be optimized. (5) e bearing capacity increased nonlinearly with the increase of the thickness of cement soil. For deep soft soil areas, the efficiency of improving bearing capacity by increasing the thickness of cement soil is low.

Data Availability
e data used to support this study are available from the corresponding author upon request.

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