Abstract

Considering that the pile foundation is widely used in soft soil, in order to improve the prediction accuracy of the vertical ultimate compressive bearing capacity of a single pile, this paper puts forward a calculation model of the pile side friction resistance in homogeneous soil based on the previous research and further deduces the calculation method of pile body internal force in layered soil under the influence of gravity field. In this research, the transmission of load along the pile body is related to the previous research on the ultimate end resistance, and an analytical calculation for the vertical ultimate compressive capacity of a single pile is proposed. The novelty of the method is that the ultimate bearing capacity is solved from the perspective of force transfer while considering its own gravity, making the results more accurate. In order to verify the reliability of the method, four engineering test piles in the soft soil area of Suzhou, China, were selected, and it is indicated that the results of the static load test and finite element simulation are consistent with the proposed method, proving that the method proposed in this paper for predicting the vertical ultimate compressive bearing capacity of single pile in soft soil area has an excellent effect.

1. Introduction

As a bearing structure, the pile is widely used in foundation reinforcement, and predicting pile bearing capacity is essential to work. Predecessors have done much research on this [14]. Conte et al. [5] proposed a simple method to evaluate the bearing capacity of piles under oblique loads. Liu et al. [6] deduced the calculation method of bearing capacity of Single Squeezed Branch Pile based on the load transfer method. Guo et al. proposed a method for calculating the ultimate vertical bearing capacity of SDRN piles considering two damage modes [7]. Pu et al. proposed a method for calculating the ultimate bearing capacity of uplifted piles in geotechnical bodies [8]. The above methods promote the progress of the prediction level of ultimate bearing capacity of monopiles, but do not sufficiently consider the influence of the pile’s gravity. In addition, some researchers have also used artificial intelligence techniques to establish the ultimate bearing capacity prediction models [914], and these methods do not require specific formulas and thus have convenient features. However, the establishment and development of prediction models require the collection of a sufficient number of valid test data.

Numerical simulations are also an important way to study the load-carrying capacity of piles [1517]. Kong et al. used FLAC 3D numerical modeling technique to gain insight into the load carrying capacity of wedge piles [18], and the results showed that the expanded bottom wedge pile is a cost-effective pile that can improve the vertical load carrying capacity under compressive and pullout loads. Elsherbiny and Naggar evaluated the compressive capacity and load transfer mechanism of helical piles in sand and clay using a validated finite element model [19]. Saggu and Chakraborty investigated the interaction between energy piles and soil in sandy soils using the finite element software Abaqus to explore the factors influencing the bearing capacity [20]. Nowkandeh and Choobbasti used finite element (FE) models to study the axial compressive performance of helical piles in sandy and cohesive soils, obtained the failure mechanism of pile groups, and provided recommendations for the design of helical pile groups using analytical equations [21]. Although the numerical simulation approach can provide more comprehensive data on the piles when they are subjected to loads, the modeling is generally complex and requires a reasonable intrinsic model to be assigned to each entity [2225].

A static load test is common and effective for predicting bearing capacity at construction sites [26, 27]. Mishra [28] summarized and compared the results of load settlement curves in static load tests in the Bombay area with methods for predicting the ultimate bearing capacity of rock-inlaid piles such as De Beer’s method (1967) [29], Chin-Kondner Extrapolation method (1970) [30], Decourt’s method (1999) [31]. Alielahi and Adampira [32] compared the results of static load test load settlement curves based on the Shahid Rajaee Port Complex Development Project in the south of Iran case study between the API [33] method and AASHTO [34] method for predicting the bearing capacity of monopiles. In general, the static load test is a costly and has a delayed prediction method. The prediction of the ultimate bearing capacity of monopiles based on the load tests of CPT and CPTU is a simple and fast prediction method [3538]. Two parameters, cone tip resistance and sidewall resistance, are obtained from the CPT tests, which can be used to classify the soil layers and estimate the strength and deformation characteristics of the soil at different depths to provide data support for the ultimate load capacity prediction. However, the final ultimate bearing capacity prediction results are derived from empirically generated equations, and empirical values are often a major source of error.

Soft soil is an inevitable form of soil. The poor bearing capacity of foundation in soft soil areas is an important reason for some engineering accidents. Whether some of the above methods can meet the accuracy requirements for predicting the monopile bearing capacity in soft soil areas needs to be further verified. The purpose of this paper is to propose a method for predicting the ultimate bearing capacity of monopiles suitable for the prediction accuracy in soft soil areas from the perspective of force transfer and considering the influence of the pile’s own gravity.

The research line of this paper starts from the characteristics of load transfer mode and pile end damage. The method of Randolph [39] can well express the transfer law of pile top load along the pile body. Based on this, the author considered the limit pile lateral friction resistance limit, considered the effect of gravity, and proposed a pile lateral friction resistance calculation model for pile sections in homogeneous soil. In this paper, this calculation model is introduced into the laminated soil, thus further revealing the transfer law of pile top load along the pile body in the laminated soil. When piles in soft soil areas are subjected to ultimate loads, the soil at the pile end undergoes a severe squeezing tendency, and Janbu [40] proposed pile end deformation characteristics that match this. Its ultimate end resistance calculation method can be well coupled with the force characteristics of the pile in the upper soft soil region when the overall shear deformation occurs. The author establishes a unified connection between the above-mentioned transfer law of pile top load and the study of ultimate end resistance and proposes an analytical equation to calculate the vertical ultimate compressive bearing capacity of monopiles in soft soil areas.

The advantage of the analytical equation proposed in this paper for calculating the vertical ultimate compressive bearing capacity of the monopile in soft soil area is that it can start from the perspective of transferring the load at the top of the pile along the pile body, consider the limitation of the ultimate pile lateral friction resistance, consider the influence of the pile’s own gravity, coupled with the ultimate end resistance calculation method proposed by previous authors according to the deformation characteristics of the pile end, which makes the final monopile vertical ultimate bearing capacity prediction results more accurate.

As shown in Figure 1, in order to test the reliability of the prediction method in this paper, four engineering test piles in Suzhou, China, were chosen to carry out static load tests on-site. As a reliable bearing capacity detection method, the static load test can get the bearing capacity of single piles according to the characteristics of the load settlement curve, and the result is consistent with the predicted value. Finally, based on the field engineering situation, the finite element method simulates the loading process of four test piles, and the result is further compared with the data output by this prediction method from the aspects of deformation characteristics and structural stress. It is found that the deformation characteristics of the soil around the pile and the variation trend of the axial force along the pile shaft obtained by the simulation are consistent with the corresponding predicted results. In general, two research methods, static load test and numerical simulation, are chosen to corroborate the reliability of the newly proposed analytical equations for calculating the vertical ultimate compressive load capacity of monopiles in soft soil areas. The corroboration shows that the method proposed in this paper is a simple, economical, and accurate method for predicting the vertical ultimate compressive load capacity of monopiles in soft soil areas, which provides valuable experience for similar projects.

2. Calculation of Vertical Ultimate Compressive Bearing Capacity of Single Pile in Soft Soil Area

The approximate derivation line of the method is given as shown in Figure 2. The vertical ultimate compressive load capacity prediction method for monopiles proposed in this paper consists of two parts, one of which is the further derivation of the load transfer law in morphogenic soils based on Randolph’s study on the internal forces of pile bodies in homogeneous soils, and the other part is Janbu’s study on the ultimate end resistance of pile ends. In this paper, the two parts were linked and the analytical equations for the vertical ultimate compressive load capacity of a single pile were constructed.

2.1. Calculation of Pile Side Friction Resistance and Pile End Resistance in Layered Soils

In homogeneous soils, Randolph [41] proposed a method for calculating pile end resistance.

Here, where is the load on the pile top, is the pile length, is the influence radius of the pile, is the radius of the pile, is the radius of the pile end, is the Poisson’s ratio of the soil layer, is the soil inhomogeneity coefficient around the pile, is the average shear modulus within the pile length, is the shear modulus of the soil at depth , is Young’s modulus of the soil layer, is the shear modulus of the soil at the pile end, and is the elastic modulus of the pile.

For cylindrical piles, , , and . In this paper, considering the influence of the pile gravity , Equation (1) can be further transformed into

Thus, the pile side frictional resistance is calculated as follows.

The pile side frictional resistance increases with the load on the pile top but is limited by the ultimate pile side frictional resistance. In Equation (4), the limit of pile side ultimate frictional resistance is not taken into account, which will be considered and corrected in the next part of this paper.

In this paper, the calculation model of the ideal elastoplastic double fold line is selected to describe the relationship between the unit pile lateral frictional resistance and the relative displacement of the pile and soil, as shown below. where is the pile side friction resistance per unit length at the relative displacement, is the body shear deformation parameter, is the pile-soil relative displacement, is the ultimate elastic displacement of the soil around the pile, is the ultimate pile side friction resistance. The calculation model is shown in Figure 3.

According to the previous research [31], the unit limit pile lateral friction resistance can be calculated using the following equation. where is the lateral pressure coefficient; is the in situ earth pressure coefficient, and it can be estimated through the equation ; is the friction angle at the pile-soil interface; is the effective internal friction angle of the soil on the pile side; is the effective overburden weight at a given depth; is the soil capacity. As for the value of , Kulhawy [42] suggested that the value of should be in the range of 0.7 to 1.2 for smooth H-shaped piles, steel pipe piles, and small displacement concrete piles. For smooth H-shaped piles, steel pipe piles, and large displacement concrete piles, the value of should be in the range of 1 to 2.

Some researchers have summarized the suggested values of [43], as shown in Table 1.

Regarding , another calculation method is given by researchers as follows [47]:

In homogeneous soils, when the pile segment is subjected to the ultimate pile side frictional resistance, as shown in Figure 4, the pile side frictional resistance per unit length of the pile segment has also reached the ultimate value.

Therefore, for a pile segment, the ultimate pile side friction resistance is where is the starting height of the pile section; is the ending height of the pile section. Let

After considering the limit pile side friction resistance limit in the pile section, then the modified pile side friction of this paper is calculated as follows. where is the load on the pile top corresponding to the limit of pile side friction resistance, and its value can be back-calculated from Equation (4), i.e.: . The model for calculating the pile side friction resistance of the pile section is shown in Figure 5.

Therefore, the pile side friction resistance ratio of this pile section is

The pile end resistance of the pile section after correction in this paper is

Thus, the pile end resistance ratio of this pile section is

In this paper, the method is introduced to the layered soil, and the pile in each layer of soil is regarded as a new pile segment, so the pile end resistance of the pile segment in the previous layer is regarded as the axial force there and also the load on the pile top of the next layer. The load transfer schematic is shown in Figure 6.

As shown in Figure 6, is the load on pile section ; is the pile side friction of pile section ; is the pile end force of pile section . Assuming that the pile corresponds to a soil layer with layers, after transfers of the pile top load along the pile body, the pile end resistance obtained is where is the pile end resistance ratio of pile section .

2.2. Ultimate Pile End Resistance

Janbu [40] proposed a calculation method for unit ultimate end resistance based on the ball hole expansion theory. where is the cohesion force of the soil at the pile end; is the average effective pressure on the lateral side at the pile end plane. As described by researchers [48], and are the dimensionless bearing capacity coefficients that reflect the cohesive force of the soil and the effect of lateral earth pressure at the pile end plane, respectively. They can be calculated by the following method. where is the effective vertical stress at the pile end; is the angle between the pile end compaction core boundary and the horizontal plane in the failure mode of pile end soil. The compressibility of soil can be determined by in situ testing methods such as penetration tests. Generally, it is about for soft clay and about for dense sandy soil [48].

Combined with the above discussion, the ultimate pile end resistance calculation method derived in this paper is as follows:

2.3. Analytical Equation of Vertical Ultimate Compressive Load Capacity of Single Pile

The load on the top of the pile is transmitted downward from the pile body to the pile end. When the pile end resistance reaches the ultimate value, the pile end will occur sliding shear failure, and its shape is shown in Figure 7. At this time, the load on the pile top is the ultimate bearing capacity of a single pile.

Combining Equation (14) and Equation (17), the analytical equation for the vertical ultimate compressive bearing capacity of a single pile in a layered soil is established as follows.

For load , if there is a corresponding pile end resistance ratio , ..., of each pile section, making Equation (18) tenable, the solution of the equation is the vertical ultimate compressive bearing capacity of a single pile.

3. Application of Engineering Cases

In this paper, four engineering test piles with different lengths, TS1, TS2, TS3, and TS4, were selected for the application of the method in the soft soil area of Suzhou, China; among them, TS1, TS2, and TS4 were nondestructible test piles,while TS3 was a destructible test pile. The type of pile is the bored pile. The design strength grade of concrete is submerged C30. the diameter of all four test piles is 0.6 m. As shown in Table 2, the calculated parameters of each soil layer in this paper are extracted from the engineering geological survey report and summarized as follows: among them, Poisson’s ratio , cohesion , and internal friction angle of the soil can be measured by triaxial tests; the density and elastic modulus of the soil can be measured by indoor geotechnical tests.

The calculated parameters of the pile end for the four test piles are summarized in Table 3. Among them, the value of , since the soil at the pile end consists mainly of clay, as suggested by Zhang [48], is 70 degrees as taken in this paper.

The data in the table were substituted into the proposed method to predict the bearing capacity of four engineering test piles. Among them, the ultimate pile side friction resistance of each pile section is shown in Figure 8. The magnitude of the vertical line represents the limit of frictional resistance that can be achieved for each pile segment at each depth. From the trend of the graph line, it is indicated that the ultimate pile side friction resistance will increase with the increase of depth, and the changing trend is the same for each test pile.

The final bearing capacities of test piles TS1, TS2, TS3, and TS4 calculated by this method are 1856.337 kN, 2179.249 kN, 2362.212 kN, and 3102.119 kN, respectively. When the test piles are subjected to ultimate compressive bearing capacity, the bottom axial force and pile end resistance of each soil layer from top to bottom are shown in Table 4.

4. Single Pile Static Load Test

In order to test the reliability of the method in this paper, static load tests were conducted in the field on the four engineering test piles mentioned above. The installation of the field static load test is shown in Figure 9.

4.1. Testing Method

The test method of the vertical ultimate compressive bearing capacity of the test piles in this project adopts the slow maintenance load method. Testing was based on the technical code for testing of building foundation piles [49]. The reaction force of the load adopted the reaction force device of the deadweight platform, and the stack load required to be applied to the deadweight platform is not less than 1.2 times the maximum load value of the pile top load. The graded loading value of the pile top load of each test pile is shown in Table 5.

4.1.1. Reading Method

After each level of load was applied to the pile top, readings of pile top settlement were taken at the 5th, 15th, 30th, 45th, and 60th minutes. Thereafter, pile top settlement readings were taken at 30-minute intervals.

4.1.2. Judgment Criteria for Settlement Stability

The settlement stability is judged by the following criterion: under each level of load, the settlement of the pile top is found to be no more than 0.1 mm per hour in two consecutive tests. When the settlement rate of the pile top reaches the relative stability criterion, then the next level of load is applied.

4.1.3. Conditions for Termination of Loading

The technical code for testing of building foundation piles stipulates the following five conditions for termination of loading [49]: (1)Under a certain level of loading, if the settlement at the top of the pile is greater than 5 times the settlement under the previous level of loading and the total settlement at the top of the pile exceeds 40 mm, the loading should be terminated(2)Under a certain level of loading, if the settlement at the top of the pile is greater than two times the settlement under the previous level of loading, and the standard of settlement stability has not been reached after 24 hours, the loading should be terminated(3)If the pile top load loading value has reached the maximum loading value required by design and the pile top settlement has reached the standard of settlement stability(4)If the anchor pile is used as an anchor pile, the loading should be terminated if the anchor pile has reached the permitted value(5)When the load-settlement curve is slowly changing, it can be loaded until the total settlement at the top of the pile is 60 mm~80 mm; when the resistance at the end of the pile is not fully developed, it can be loaded until the cumulative settlement at the top of the pile is more than 80 mm

4.2. Analysis of Results

The load-settlement curves obtained from the static load test are shown in Figure 10. It can be seen that the curves corresponding to TS1, TS2, and TS4 belong to slow deformation, while that of TS3 belongs to steep settlement type. According to the technical code for testing of building foundation piles [49], the loading load on TS1, TS2, and TS4 test piles should be terminated when the designed maximum loading is reached and the settlement of the pile top meets the standard of relative stability. For TS3, when the load reaches the maximum loading value of 2400 kN, the load-settlement curve shows obvious steep drop characteristics, which means that the test pile has been damaged. Therefore, the vertical ultimate compressive bearing capacity of the TS3 test pile should be taken as the loading value of the previous test pile without damage, i.e., the load value corresponding to the starting point of an obvious steep drop, which is 2160 kN.

Table 6 summarizes the vertical ultimate compressive load bearing capacity obtained from the static load test and the prediction method of this paper. From the comparison of the results, the predicted ultimate bearing capacity of the static load test is conservative, while the pile actually does not reach the ultimate bearing state at this time. Therefore, in theory, the results calculated by the method in this paper should be greater than the results of the static load test, which is also fully reflected in the table. Further analysis of the data shows that the difference between the vertical ultimate compressive load capacity obtained from the static load test and the method in this paper is not much, which fully demonstrates the high accuracy of the calculation method. In addition, it is also obtained that the longer the length of the pile, the stronger the law of its vertical ultimate compressive load capacity. Therefore, the results of the field static load tests of TS1, TS2, TS3, and TS4 can coincide with the results of the method in this paper.

As shown in Figure 11, this paper also incorporated the comparison of the results of the method without considering the effect of the pile’s own gravity. It can be clearly observed from the Figure 11 that the vertical ultimate compressive load bearing capacity results obtained by the method without considering the effect of the pile’s self-gravity are generally larger than considering its self-gravity, which is more obvious compared with the results of static load test detection, which greatly increases the error of prediction results. Therefore, it is a necessity to consider the pile’s self-gravity to calculate the vertical ultimate compressive load-bearing capacity.

In addition, the static load test results contain both the damaged engineering test piles, including undamaged test piles, which are more valuable for comparison. If the results of the damaged test piles are used to evaluate the accuracy of the calculation results in the static load test, the loading value of the TS3 test pile is 2400 kN, and the vertical ultimate compressive bearing capacity calculated by the proposed is 2362.212 kN, the error rate between them is very small, roughly estimated to be only 1.57%. Therefore, the prediction method of the vertical ultimate compressive load capacity for single pile in soft soil area proposed in this paper is an excellent method.

5. Finite Element Simulation of Test Pile Loading Process

In order to further test the reliability of the prediction method, the finite element simulation of the loading process was carried out for four test piles with the same dimensions of 20 m in length, 20 m in width, and 50 m in height. The pile material ontology model adopts the elastic model, considering the pile concrete material is submerged C30, its elastic modulus is taken as 30 Gpa, its capacitance is 27 kN/m3 and Poisson’s ratio is 0.16. The soil ontology relationship adopts the Mohr-Coulomb model. The soil material parameters defined by the model include soil thickness, soil Poisson’s ratio, soil cohesion, angle of internal friction, soil gravity, and soil modulus of elasticity from Table 2. The model should be empirically set up for the contact interface between pile and soil, and the interface can well simulate the relationship between pile lateral resistance and relative displacement through the embedded displacement-friction function. The model has fixed X-directional constraints on the left and right, fixed Y-directional constraints on the front and back, and three fixed XYZ constraints on the bottom. The finite element simulation considered the effect of the gravity field, and the simulation process was divided into three stages, including the initial stress stage, piling stage, and loading stage. In order to fit the loading process of the static load test to the greatest extent, the displacements of the first two stages are cleared to zero. The finite element model is shown in Figure 12.

As can be seen from Figure 13, the settlements of TS1, TS2, TS3, and TS4 obtained from the finite element simulation under all levels of pile top loading do not differ much from the settlements obtained from the static load test. On the whole, the results of the finite element simulation are slightly larger. This indicates that the finite element simulation is a reliable and safe method for settlement prediction. In addition, the fact that the load settlement curves obtained from the static load test and the finite element simulation can be fitted also indicates that the finite element simulation of these four test piles is a valid simulation, which can provide more comprehensive information on the internal force and pile end deformation as a reference to verify the reliability of the method in this paper.

The settlement deformation in Z-direction obtained for each engineering test pile under the corresponding load is shown in Figure 14. The results show that the settlement of the soil near the pile perimeter is the largest, and a new compacted dense core area will be formed at the pile end, and a part of the soil will be settled below the compacted dense core area is relatively small. Combined with Figure 7, the soil settlement deformation trend near each pile end is consistent with that revealed by Janbu. Therefore, the ultimate pile end resistance calculation method referenced by the analytical equations proposed in this paper is appropriate.

Figure 15 shows the variation of axial force along the pile body at loading values of 1750 kN, 2100 kN, 2160 kN, and 2700 kN for test piles TS1, TS2, TS3, and TS4 extracted from the finite element model, respectively, and also shows the variation of axial force along the pile body at the vertical ultimate compressive load capacity for each test pile obtained by the calculation method in this paper. From a single test pile, the axial force of each part of the pile calculated by the method in this paper is generally larger compared with the results of finite element simulation, which is consistent with the actual situation, because the top load of the pile corresponding to the finite element simulation at this time does not reach the ultimate bearing capacity state that the pile can withstand. However, the trend of the axial force along the pile body is the same, and the axial force decreases gradually along the pile body. The reason is that the pile side frictional resistance of the shallow soil to the pile body is smaller, and a larger axial force is required to keep the pile body in a state of force equilibrium. The trend of the axial force along the pile body further illustrates that the vertical ultimate compressive load capacity calculation method proposed in this paper is a practical method in soft soil areas from the structural force characteristics.

6. Limitations and Future Works

The reliability of the analytical equations proposed in this paper for calculating the vertical ultimate compressive load capacity of monopiles in soft soil areas is verified from both numerical simulations and finite elements. The results show that the method in this paper is a simple and accurate method. This means that the method can be widely used for the prediction of ultimate bearing capacity of monopiles in soft soil areas, thus making it possible to reduce the number of static load tests, which can reduce certain costs for the project. However, the method also has some limitations. (1)From the basic principle of the equation, the method can be extended to other types of soil areas. However, it is also necessary to consider whether the damage deformation characteristics of pile ends in other types of soil areas are similar to those in soft soil areas, and sufficient test data are needed to verify the reliability of the method(2)In this paper, the ideal elastoplastic bifold model was chosen to describe the relationship between unit pile lateral friction resistance and pile-soil relative displacement, which is a basic and common model, but the model does not consider the softening effect between pile and soil. In future research, we can try to express the relationship between unit pile lateral frictional resistance and pile-soil relative displacement with a model that considers more influencing factors

Therefore, the new analytical equation proposed in this paper is used to calculate the vertical ultimate compressive bearing capacity of monopiles in soft soil areas with certain engineering value, but there is a need to do further research to make it produce more application value.

7. Conclusion

Based on all the research contents and results obtained in this study, the following conclusions can be drawn: (1)Based on Randolph’s study on the internal force of piles, a model for calculating the frictional resistance of pile segments in homogeneous soils considering the effect of self-gravity is proposed in this paper, and a method for calculating the internal force of piles in laminated soils is further derived. Combined with Janbu’s research on ultimate end resistance, an analytical equation is proposed for calculating the vertical ultimate compressive bearing capacity of monopiles in soft soil areas. In order to verify the reliability of the method, static load tests and finite element simulations were carried out on four engineering test piles in the soft soil area of Suzhou, China(2)The results show that the results of TS1, TS2, TS3, and TS4 for predicting the ultimate bearing capacity of monopiles in soft soil areas are in good agreement with the static load test results, and the error rate is only 1.57% if the results of the damaged test pile TS3 are taken to measure the deviation of the results. Finite element simulation was revealed from the two perspectives of pile end deformation characteristics and load transfer. the ultimate bearing capacity prediction method proposed in this paper based on Randolp’s study of pile internal forces and Janbu’s study of ultimate end resistance fits the actual situation and is an accurate and reliable method, which has important reference significance for other similar projects

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the Ministry of Industry and Information Technology of the People's Republic of China IoT Integrated Innovation and Convergence Application Project Fund (Grant Nos. 2018470) and Key R&D Projects in Shaanxi Province (No. 2020SF-373).