Analytical Solution for Estimating Bearing Capacity of a Closed Soil Plug: Veriﬁcation Using An On-Site Static Pile Test

: When the open-ended pile penetrates the soil layer, the resistance generated by the soil plug cannot be ignored. A pile with a full-size pressure sensor installed at pile tip can detect resistance more accurately than a microsensor when the pile penetrates into the soil. In this paper, the pile installed full-size pressure sensor was used for penetration test and the relationship between formation parameters and pile tip force is obtained. Using the solution of the Kelvin problem in inﬁnite space and the plane stress distribution function, the analytical solution of the bearing capacity of the soil plug is derived under the condition that the displacements of the bottom of the pile and the soil plug are consistent. The results show that the ultimate stress of the soil plug is closely related to the pile diameter and pipe thickness. The bearing capacity of the soil plug is closely related to the properties of the soil layer. The analytical solution of the bearing capacity of the soil plug has a linear relationship with the formation parameters SPT and CPT. The analytical solution of the ultimate bearing capacity of the soil plug has been veriﬁed by ﬁeld test data and has a good match with the geometric dimensions of the pile tip and the formation parameters. the discrete element soils multi-size piles pile ﬃ This study presents a new analytical solution for exploring relationship of bearing capacity of plug with that of properties of soil layers. Using the solution of the Kelvin problem in inﬁnite space and the analytical solution obtained by the plane stress distribution function, the proportion of the penetration force of the soil plug is directly related to the diameter of the pile tip and the pipe thickness. When the diameter is ﬁxed, the greater the pipe thickness, the smaller the proportion of the ultimate load-bearing stress of the soil plug. Through the static pressure penetration test of the pile foundation installed with a full-section pressure sensor, the obtained formula can be well used in the bearing capacity of the closed soil plug. The results show that the bearing capacity of the soil plug is closely related to the properties of the soil layer. The analytical solution of the ultimate bearing capacity of the soil plug has a linear relationship with the formation parameters SPT, CPT. The conclusion of this article can provide a reference for the design.


Introduction
As one of the deep foundations, piles have the advantages of high bearing capacity [1,2], good stability, small settlement and small environmental impact, and are widely used in construction engineering and marine engineering. Among them, the length of the offshore pile foundation often exceeds 100 m, and the diameter can be larger 2 m. Such a super-long and large-diameter pile foundation needs sufficient penetration force to reach the design depth. Piles usually have two forms: open-ended and closed-ended. Open-ended piles are easier to drive to the design depth than closed-ended piles. However, when the open-ended penetrates, a soil plug will be formed inside the pile [3]. When the soil plug is closed, the open-ended pile behaves similarly to the closed-ended pile [4]. The presence of soil plugs influence the form of pile load transfer [5,6].
Numerous research have been conducted related to formation of soil plugs in open-ended piles from various perspectives such as, (1) Parameters in the dynamic process of forming soil plugs, such as the diameter of the pile [7], the penetration method [8], the penetration speed [9,10] and the condition of the soil layer [11,12]; (2) IFR (incremental filling ratio) [13] and (3) Dynamic effects [14] through the analysis of parameters to more accurately calculate the effects of soil plugs. At the same time, more experimental methods such as centrifuge test [15], model test research [16][17][18], and field test [19][20][21][22] are also adopted to study the soil plug effect. The scale model test and centrifugal test are affected by the particle size, and there are two obvious defects. (1) The number of soil particle samples on the micro sensor is obviously too small, and the data is prone to extremely fluctuating, (2) The test soil samples are mostly remodeled soil, it is difficult to restore the state of the on-site soil samples.
Some scholars [23][24][25][26] explored the analytical solution of pile-soil interaction during penetration process, and analyzed the stress in the process of half-space and full-space penetration, and deduced the mechanical formulas.Few research [27][28][29][30] have utilized model piles for load tests and field load tests. They believe that SPT and CPT are closely related to soil plugs. Many studies [31][32][33][34][35][36][37][38][39][40] utilized the finite element and discrete element methods to numerically evaluate the soil plug effect and to further analyze the plugging mechanism. In the finite element analysis, the Euler-Lagrange model is mostly used to simulate the penetration and plugging mechanism of the pile. The discrete element can better simulate the displacement movement of particles.
The traditional measurement method is to arrange a miniature pressure sensor on the inner or outer surface of the pile [4,13,17,41]. Due to the small size of the sensor, the small number of particle samples on the contact surface of the sensor is easy to cause violent data fluctuations. Past studies have focused on the formation of soil plugs and the stress distribution on it. The distribution of stress at the pile tip is evenly distributed stress, which does not satisfy the condition of consistent displacement at each point of the pile tip. Therefore, to accurately obtain the bearing capacity of the soil plug, the key is to accurately measure the bearing capacity of the pile tip and then rationally divide the bearing capacity of the soil plugs and the pile pipe. The objective of this study is to develop a new analytical solution for exploring relationship between bearing capacity of the soil plug under closed condition with that of properties of soil layers. In order to achieve this, the full-section pressure sensor was used to accurately measure the penetration force of the pile. The relationship between the bearing capacity of the soil plug and the properties of soil layer is studied based on the stress function satisfying the displacement condition of the pile tip. An on-site static pressure pile test was also conducted in order to verify such relationship.

New Analytical Solution of Closed Soil Plugs
When the open-ended pile penetrates the soil layer, the soil is squeezed into the pile and a soil plug is formed inside the pile. When the soil plug is closed, the bearing capacity of the open-ended pile tip is the same as that of the closed-ended pile tip, and the internal height of the soil plug does not change with the penetration depth. The bearing capacity of the closed soil plug can be derived with the analytical solution of the bearing capacity of the closed pile tip. The soil plug is often formed in deep soil. The penetration force at the pile tip can be deduced by the solution of the Kelvin problem in infinite space. If the section of the soil plug is to be studied in detail, a stress function is needed, because the penetration force of the Kelvin problem is a concentrated load.

Analytical Solution of Penetration Force and Displacement in Infinite Space
According to the solution of Kelvin problem, the relationship between the vertical displacement and the vertical force is shown in Equation (1) and Figure 1: among them, Where u z is the vertical displacement, r = x 2 + y 2 + z 2 , E is the compression modulus of soil, υ is the Poisson's ratio of the soil, and P z is the penetration force of the pile tip. Taking the center of the coordinate axis as the center of the pile tip and in the infinite body, the vertical displacement is simplified as: Taking the center of the coordinate axis as the center of the pile tip and in the infinite body, the vertical displacement is simplified as:

Analytical Solution of Pile Tip Stress in Infinite Space
The force of the Kelvin problem is a vertical concentrated force, and the vertical displacement caused by multiple vertical concentrated forces can be obtained by superposition. To study the penetrating force at the pile tip, it is not possible to use the method of superposition of a limited number of vertical concentrated forces, it is necessary to use the stress distribution function, assuming the stress distribution function:

Analytical Solution of Pile Tip Stress in Infinite Space
The force of the Kelvin problem is a vertical concentrated force, and the vertical displacement caused by multiple vertical concentrated forces can be obtained by superposition. To study the penetrating force at the pile tip, it is not possible to use the method of superposition of a limited number of vertical concentrated forces, it is necessary to use the stress distribution function, assuming the stress distribution function: Among them, p 0 is the compressive stress amplitude, a is the radius of the pile tip, r is the distance from any point to the axis, and D is the area of the pile tip with radius a.
Due to the rotational symmetry of the pile tip, the vertical displacement of a point is only related to the distance r from the center point. For this purpose, only the displacement of the point on the x axis needs to be calculated. Take point A coordinate (x, y) on the x axis and the B coordinate (x , y ) at any point on the pile tip, according to the positional relationship shown in Figure 2.
Stress at point B: The vertical displacement of point A is equal to the superposition of the influence of any point on the pile tip, according to Equation (6): Equation (7) calculates: The vertical displacement of each point in the contact area of the closed-ended pile tip is consistent with the displacement coordination condition.
Calculated by Equation (4) Exchange points: The simultaneous Equations (8) and (9) yield: The vertical displacement of point A is equal to the superposition of the influence of any point on the pile tip, according to Equation (6): Position relations of various points in the pile tip area. Among them, A is any point on the x axis in the circle, B is any point in the circle, a is the diameter, s is the distance between the two points of A and B, s 1 is the distance between the extension line of AB and the circle at a point, and t is the distance from B to the center O, ϕ is the angle between AB and the x axis, and θ is the angle between OB and the x axis.

Stress Distribution of Closed Soil Plug
Since, the displacement of the pile tip is the same under closed soil plug, the stress distribution can be calculated using Equation (1). However, the action range needs to be calculated separately. The section of the soil plug is shown in Figure 3.

Stress Distribution of Closed Soil Plug
Since, the displacement of the pile tip is the same under closed soil plug, the stress distribution can be calculated using Equation (1). However, the action range needs to be calculated separately. The section of the soil plug is shown in Figure 3. Figure 3. Soil plug section.Among them, a is the outer diameter of the pile, b is the inner diameter of the pile, and c is the length obtained by ab and the Pythagorean theorem, θ is the angle between a and b. When the soil plug is closed, the proportion of the penetration force of the soil plug is directly related to the diameter of the pile tip and the pipe thickness. With a fixed diameter, the greater the pipe thickness and the larger is θ and smaller is the proportion of ultimate load stress of the soil plug. The greater is the proportion of the bearing stress of the pipe (as shown in Figure 4).
When the soil plug is closed, the proportion of the penetration force of the soil plug is directly related to the diameter of the pile tip and the pipe thickness. With a fixed diameter, the greater the pipe thickness, and the larger is θ, and the smaller is the proportion of the ultimate load stress of the soil plug; then the greater is the proportion of the bearing stress of the pipe (as shown in Figure 4).

Test Overview
An on-site static pressure pile test was carried out to verify the relationship between the bearing

Test Overview
To verify the relationship between the bearing capacity and displacement of Closed soil plugs, an on-site static pressure pile test was carried out. The test site is located in Dongying, Shandong Province, about 20 km away from the Yellow River. It belongs to the Yellow River impact plain and has clear soil layers. It is mainly composed of silt and silty clay. The soil layer is shown in Table 1. The pile diameter is 400 mm and the pile length is 12 m. A 100t full-section pressure sensor is placed at a pile tip. In order to minimize error due to signal attenuation caused by long cable (60 m), a RS-485 signal transmission is used. Static pressure penetration was applied with the average penetration speed is 1.2 m/min. The test position and sensor are shown in Figure 5.

Test Results
The penetration test of the static pressure pile is divided into 9 strokes. The penetration depth of each stroke is about 1.4 m. The relationship between the resistance of the pile tip and the penetration depth is shown in Figure 6. After the start of static pressure penetration, there is a stage of stored force pressurization at each stroke before the pile tip pierces the soil layer. During this period, the deformation of the soil body is mainly compression, and the settlement is less. When the pile tip force breaks through the critical value, the pile foundation quickly penetrates the soil layer. Table 2 shows the bearing capacity of pile tips at each stroke.

Test Results
The penetration test of the static pressure pile is divided into 9 strokes. The penetration depth of each stroke is about 1.4 m. The relationship between the resistance of the pile tip and the penetration depth is shown in Figure 6. After the start of static pressure penetration, there is a stage of stored force pressurization at each stroke before the pile tip pierces the soil layer. During this period, the deformation of the soil body is mainly compression, and the settlement is less. When the pile tip force breaks through the critical value, the pile foundation quickly penetrates the soil layer. Table 2 shows the bearing capacity of pile tips at each stroke.

Relationship Between Soil Plug Resistance and Soil Layer
According to the relationship (Figure 6), between the penetration resistance at the pile tip and the soil layer, the initial penetration resistance in the soil layer is greater than the secondary penetration in the soil layer. After the initial penetration into the soil layer, the soil around the pile undergoes large deformation, and the part of soil layer, where it penetrates, again has been affected by the previous penetration. A similar phenomenon occurs in the strokes 1~2, 6~7, 8~9. It is proved that the maximum resistance of the layered soil pile tip occurs at the peak value generated, when the layered soil penetrates initially.  According to Equation (12), the analytical solution of closed soil plug is calculated and listed in Table 2. The Equation (15) obtained by fitting the analytical solution of the closed soil plug stress in Table 2 to CPT has a good correlation, R 2 = 0.96502.
The Equation (16) obtained by fitting the analytical solution of the closed soil plug stress in Table 2 to SPT has a good correlation, R 2 = 0.96176, as shown in Figure 7. This test pile passes through the plain fill, silt, and silty clay, and it can be considered that Equations (15) and (16) have good applicability to these three soil layers. According to Equation (12), the analytical solution of closed soil plug is calculated and listed in Table 2. The Equation (15) obtained by fitting the analytical solution of the closed soil plug stress in Table 2 to CPT has a good correlation, R 2 = 0.96502.
The Equation (16) obtained by fitting the analytical solution of the closed soil plug stress in Table  2 to SPT has a good correlation, R 2 = 0.96176, as shown in Figure 7. This test pile passes through the plain fill, silt, and silty clay, and it can be considered that Equations (15) and (16) have good applicability to these three soil layers.

Analytical Verification
When the soil plug is closed, the displacement of each point at pile tip is the same. Equation (8) derived from Equation (3) can ensure that the displacement of each point at pile tip is consistent. It can be seen from Equation (3) that the force of each point at pile tip is related to the distance from the center point. Many scholars have conducted theoretical derivation and experimental research on the formation of soil plugs. However, it is more difficult to separate the pile pipe stress and the soil plug stress at pile tip. Only a few scholars have successfully distinguished it, the experimental data of De Nicola [13] and Paik Kyuho [2] are shown in table 3-5, but it is still difficult to test the values of soil plug stress and plie pipe stress at pile tip when the soil plug is closed. During the penetration of the pile, the experimental value of the stress distribution between the soil plug and the pipe still has a certain reference.

Analytical Verification
When the soil plug is closed, the displacement of each point at pile tip is the same. Equation (8) derived from Equation (3) can ensure that the displacement of each point at pile tip is consistent. It can be seen from Equation (3) that the force of each point at pile tip is related to the distance from the center point. Many scholars have conducted theoretical derivation and experimental research on the formation of soil plugs. However, it is more difficult to separate the pile pipe stress and the soil plug stress at pile tip. Only a few scholars have successfully distinguished it, the experimental data of De Nicola [13] and Paik Kyuho [2] are shown in Tables 3-5, but it is still difficult to test the values of soil plug stress and plie pipe stress at pile tip when the soil plug is closed. During the penetration of the pile, the experimental value of the stress distribution between the soil plug and the pipe still has a certain reference.
Tables 3-5 lists the parameters of the open-ended pile and the results calculated by Equations (13) and (14). The results are made into Figure 8. According to the results shown in Figure 8, more than 90% of the test data are within the scope of the analytical solution proposed in this paper, the formulas in this paper can be used as a reference for the calculation of soil plug stress.  (13) and (14). The results are made into Figure 8. According to the results shown in Figure 8, more than 90% of the test data are within the scope of the analytical solution proposed in this paper, the formulas in this paper can be used as a reference for the calculation of soil plug stress.

Conclusions
Due to the difficulty of the test, only a few tests measured the pipe stress and the soil plug stress at the same time, which brought great difficulty to the verification work of this paper. If the discrete element method can be used to simulate the soils and multi-size piles better, the force at pile tip will no longer be difficult to distribute when the pipe pile is jacked or the soil plug is closed, which is also a new research direction.
This study presents a new analytical solution for exploring relationship of bearing capacity of plug with that of properties of soil layers. Using the solution of the Kelvin problem in infinite space and the analytical solution obtained by the plane stress distribution function, the proportion of the penetration force of the soil plug is directly related to the diameter of the pile tip and the pipe thickness. When the diameter is fixed, the greater the pipe thickness, the smaller the proportion of

Conclusions
Due to the difficulty of the test, only a few tests measured the pipe stress and the soil plug stress at the same time, which brought great difficulty to the verification work of this paper. If the discrete element method can be used to simulate the soils and multi-size piles better, the force at pile tip will no longer be difficult to distribute when the pipe pile is jacked or the soil plug is closed, which is also a new research direction.
This study presents a new analytical solution for exploring relationship of bearing capacity of plug with that of properties of soil layers. Using the solution of the Kelvin problem in infinite space and the analytical solution obtained by the plane stress distribution function, the proportion of the penetration force of the soil plug is directly related to the diameter of the pile tip and the pipe thickness. When the diameter is fixed, the greater the pipe thickness, the smaller the proportion of the ultimate load-bearing stress of the soil plug. Through the static pressure penetration test of the pile foundation installed with a full-section pressure sensor, the obtained formula can be well used in the bearing capacity of the closed soil plug. The results show that the bearing capacity of the soil plug is closely related to the properties of the soil layer. The analytical solution of the ultimate bearing capacity of the soil plug has a linear relationship with the formation parameters SPT, CPT. The conclusion of this article can provide a reference for the design.