Abstract

A concave slope is a common type of slope. This paper proposes a simplified method to study the effect of a clay concave slope on laterally loaded piles. The hyperbolic - curve model is selected as the lateral pile-soil interaction model of the concave slope. Considering the two angles of the concave slope, the variation of the ultimate soil resistance with depth is divided into two parts, and the ultimate soil resistance varies nonlinearly with depth. The reduction factor method and normalization method are used to obtain the initial stiffness. The theoretical results will be compared with the calculation results of the 3D FE analysis to prove the rationality of this method. Finally, the simplified method is used to analyze the response of laterally loaded piles under different parameters.

1. Introduction

Pile foundation is one of the most commonly used foundations in bridge engineering, offshore drillings, and offshore wind turbines. These pile foundations are often used on sloping ground, such as river valleys and the seabed [1, 2]. The pile foundation will be subjected to lateral loads caused by traffic loads, lateral wind, and waves.

The bearing capacity of pile foundations depends on the bearing capacity of the rock and soil around the pile. There are three main approaches to study the bearing capacity of rock and soil around the pile: theoretical methods [3, 4], numerical simulations [5], and experimental methods [6]. In the past few decades, the - curve method is often used to study the response of pile foundation bearing lateral load. The main research includes the influence of laterally loaded piles in flat ground and sloping ground. For the flat grounds, many scholars and institutions proposed - curves for different types of soil [79]. For the sloping ground, the soil in front of the pile is weakened, and the damage model of soil is different from that in the horizontal ground [1012]. Therefore, Reese et al. [13] proposed - curves that were suitable for sand and clay sloping ground, respectively. Based on the 3D FE analysis, Georgiadis and Georgiadis [14, 15] obtained the - curves suitable for clay sloping ground. On this basis, the - curves of clay sloping ground were proposed, which considered the distance between the slope and pile. But all the - methods mentioned above only considered level ground and single-angle slope.

However, due to the influence of external factors such as rain erosion and soil accumulation, the slope has more than two angles. Wu et al. [16] and Fan et al. [17] pointed out that slope shapes could be roughly divided into four types, which were the straight type with a single angle, the convex type (the upper slope angle is smaller than the lower slope angle), the concave type (the upper slope angle is greater than the lower slope angle), and a mixed type.

The - method is widely used in engineering because of its simple calculation and short calculation time. However, compared to a slope with a single angle, the distribution law of the concave slope’s ultimate soil resistance and initial stiffness will change. But unfortunately, the existing - curves of sloping ground can only consider the change law of ultimate soil resistance and initial stiffness under a single angle. It leads to errors in the analysis of the horizontal bearing characteristics of piles using the existing - curve. Therefore, it is necessary to carry out further research on the lateral load characteristics and calculation methods of concave slope piles.

This paper focuses on giving a nonlinear analysis method considering the bearing characteristics of the laterally loaded pile in concave sloping ground. In this method, the pile-soil interaction model adopts the hyperbolic - curve model. And the calculation formula of the ultimate soil resistance and the initial stiffness varying with depth in the concave slope is given. Then, the derived - curve is brought into MATLAB for a differential calculation to obtain the nonlinear response of the pile under a lateral load. To prove the correctness, the method is verified by the results of three-dimensional finite element analysis considering the concave clay slope. Furthermore, this article discusses the influence of the different slope angles and upper slope height and obtains the response law of piles under the influence of different parameters.

2. Establish a - Curve of Clay Concave Slope

At present, the - curve method, which regards soil as a nonlinear spring, is widely used to study the lateral load response of the horizontal ground piles. There are many kinds of mathematical models of the - curve. One of the most widely used is the hyperbolic - curve [1821]. The expression is shown as follows: where  = ultimate soil resistance of the soil along the pile,  = initial stiffness of the foundation.

Yang [22] compared different types of - curves and concluded that the hyperbolic - curve has the best fitting effect with the data obtained from field experiments. Therefore, this article adopts a hyperbolic - curve.

2.1. Basic Assumption

To analyze the influence of the concave slope which has two angles on the horizontal bearing characteristics of the pile, an analysis model is established, as shown in Figure 1.

For simplicity, the following assumptions are made: (1)The slope is stable without a sliding surface, and slope failure and instability are not considered in the calculation(2)The soil resistance along the pile changes nonlinearly with the increase in lateral displacement. When the ultimate soil resistance is achieved, its value remains constant as the increase in lateral displacement(3)This article only considers concave slopes, the upper slope angle is larger than the lower slope angle,

2.2. Ultimate Soil Resistance of Concave Slope Varying with Depth

For the ultimate soil resistance , its expression is as follows: where is the ultimate soil-resistance-bearing factor, is the undrained shear strength of soil, is the pile diameter.

Equation (3) shows that the value of changes nonlinearly with the increase in depth [14]. The initial value of is at the ground surface, and the maximum value is :

The value of , , and are related to the adhesion factor of the pile-soil interface [23, 24]).

The adhesion factor is the ratio of the interface shear strength to the undrained shear strength of the soil. In this article, can be obtained from Figure 2 [14].

When , we have

When , we have

When , we have

is a dimensionless parameter that changes with the adhesion factor. For smooth piles, is 0.55 (); for rough piles, is 0.4 ():

is the ultimate soil-resistance-bearing factor at the ground surface:

is the ultimate soil-resistance-bearing factor based on the deep soil flow failure model by Randolph and Houlsby [25]: where .

To obtain the ultimate soil resistance along the pile in concave sloping ground, the ultimate soil-resistance-bearing factor is the key.

Jiang et al. [26] found that with the increase in the lateral displacement of the pile in sloping ground, the stress and strain of the soil show a wedge-shaped distribution in the three-dimensional numerical simulation. Similarly, Georgiadis and Georgiadis [15] found the same rule when studying the influence of the distance between the pile and slope on laterally loaded piles, as shown in Figure 3.

It shows in Figure 3(a) that when the pile head deflection is small, the soil has a damaged wedge in a shallow soil. The sliding surface of the damaged wedge only intersects with the level ground. Therefore, the slope will not affect the ultimate soil resistance near the mudline. As the pile head deflection augments, the depth of the damaged wedge increases. When the depth of the damaged wedge reaches , the sliding surface of the wedge spreads to the slope. As shown in Figure 3(b), when the depth of the damaged wedge is less than , the expression of the ultimate soil-resistance-bearing factor is the same as that of the flat ground since the pile is only affected by the level ground. When the sliding surface of the damaged wedge intersects with the slope surface, the ultimate soil resistance along the pile is weakened compared to the level ground. The expression of needs to consider the influence of the slope.

Georgiadis and Georgiadis [15] believed that the depth of the damaged wedge would change with the increase of lateral load, and it had a nonlinear relationship with the distance of pile from the crest of the slope. When the sliding surface intersects with the slope crest, the critical depth is expressed as follows: where indicates the distance from the pile core to the crest of the slope. When , the pile is at the crest of the slope.

For concave slopes, this article will use the same method to get the expression of , as shown in Figure 4.

Figure 4(a) shows that the sliding surface of the damaged wedge intersects the slope with the angle of when , and the ultimate soil resistance is only affected by the upper slope. When the depth reaches , the sliding surface of the wedge intersects with the intersection of the upper slope and lower slope. When the depth is greater than , the sliding surface intersects the slope with the angle of , and the lower slope begins to affect the ultimate soil resistance. It is pointed out in Figure 4(b) the changing law of in the concave slope with two angles. The expression of is only controlled by when the depth of the damaged wedge is less than . When the depth is greater than , the expression of is affected by . It is reasonable that the value of increases at relative to that of the single-angle slope () because the lower slope increases the volume of soil in front of the pile. To determine the critical depth , this paper adopts the modified critical depth expression:

Because the slope is a homogeneous medium, the expression of the ultimate soil resistance is continuous and derivable when the depth increases, and there will be no sudden change point. Therefore, when the depth reaches the critical depth (), the piecewise function of should be the same value , as shown in Figure 4(b). When the value of is , the corresponding depths are for the slope with angle and for the slope with angle . And the difference betweenandis defined as, as shown in Figure 5.

When the depth is less than , the expression of is as follows:

When the depth is greater than , the expression of is as follows:

2.3. Initial Stiffness of Concave Slope Varying with Depth

Carter [27] proposed the expression of the initial stiffness of the clay level ground, and believed that the initial stiffness is related to the initial elastic modulus of the soil, Poisson’s ratio, and other factors: where is the pile diameter; is the pile diameter reduction factor, usually taken as 1; is the bending stiffness of the pile; is the Poisson’s ratio of the soil; is the initial elastic modulus of the soil.

Kondner and Robertson et al. [28, 29] proposed the equation of elasticity modulus and believed that the initial elastic modulus can be related to the elastic modulus at fifty percent of the failure stress : where is deviatoric stress; is elasticity modulus; is deviatoric failure stress; is the ratio of deviatoric failure stress over deviatoric ultimate stress, usually taken equal to 0.8.

Setting , , and , equation (14) becomes

Equation (16) reveals that the initial stiffness of the foundation is only related to the characteristics of the soil and the pile. The normalization method is commonly used to obtain the initial stiffness expression of the sloping ground. On the basis of the initial stiffness expression of the level ground, many scholars adopt the reduction factor to establish the initial stiffness expression of the sloping ground [14, 30]: where is the initial stiffness of the sloping ground; is the initial stiffness of the level ground.

When studying the influence of the initial stiffness of the sloping ground, Georgiadis and Georgiadis [14] proposed that the initial stiffness reduction factor is related to the depth of the slope, as shown in equation (18). The reduction factor changes linearly with the increase of depth, as shown in Figure 6. For slopes at an arbitrary angle, when the depth is less than , the slope weakens the initial stiffness, and the degree of weakening decreases when the depth becomes greater. When the depth is greater than , the initial stiffness of the sloping ground is the same as that in level ground, and remains constant with the increase of depth:

As shown in Figure 7, for a concave slope, when the value of exceeds (red line (1) in Figure 7(a)), the initial stiffness of the concave sloping ground at reaches the value of that in level ground. The lower slope does not affect the initial stiffness, and the reduction factor is only controlled by the upper slope (blue line in Figure 7(b)). When is 0, the concave slope becomes a single-angle slope with (red line (2) in Figure 7(a)). The reduction factor is controlled by the slope with angle (red line in Figure 7(b)). When ranges from 0 to (red line (3) in Figure 7(a)), the reduction factor varies with depth in the range between the two limit conditions (shaded part in Figure 7(b)). Georgiadis and Georgiadis [14] point out that the reduction factor has a small effect on both pile head deflection and maximum bending moment. This paper assumes that the initial value of the reduction factor changes uniformly from to as increases from 0 to .

When , we have

When , we have

And the initial stiffness of concave sloping foundation can be expressed as follows:

To verify the rationality of the above assumptions, three cases of reduction factors are used to obtain the pile head load-displacement curve. is calculated by the method in this paper, is obtained by the calculation method considering the single-angle slope with angle , and is obtained from the calculation method considering the single-angle slope with angle . Undrained shear strength  kPa, and elastic modulus of soil is 14 MPa. Pile diameter  m; pile length  m, 12 m, and 20 m; and elastic modulus of piles  kPa. ; and upper slope height  m. The load-displacement curve of the pile head under different conditions is shown in Figure 8. It indicates that the load-displacement curve of the pile head is almost unchanged even if is taken as two limit conditions (, ). For the pile length  m, the maximum discrepancies are only 16%. The above simplifying method is sufficient for determining the reduction factor.

3. Approach Verification

A series of three-dimensional finite element analysis models of the laterally loaded pile in a concave slope is established to verify the correctness of the calculation method in this paper.

3.1. Establishment and Verification of 3D FEA Model

All the basic parameters of the 3D finite element analysis model in this paper are obtained from the literature [5, 14]. The basic parameters of the pile are set as follows: pile length L = 20 m; diameter ; elastic modulus of pile  kPa; Poisson’s ratio of pile is 0.1; and the density of pile is 2500 kg/m3. Piles are all embedded in the soil, the load is applied at the pile head, and the pile head is free. Slope angles are 20° and 40°. The undrained shear strength of the soil is 70 kPa, elastic modulus of the soil at 50% ultimate stress is 14 MPa, Poisson’s ratio of soil is 0.49, and bulk unit weight is 18 kN/m3. C3D8R grids are used for piles and soil. In addition, there is more detailed meshing around the pile. The number of meshes in all models is approximately 25,000. The bottom boundary of the model is fixed in all directions, and the other boundaries, except for the top boundary, are only fixed in the normal direction. The soil is established based on the Mohr-Coulomb model. The contact surface between pile and soil adopts normal behavior and tangential behavior. Normal behavior is set as a “hard” contact mode, and tangential behavior is set as “ penalty” function. The “penalty” factor of the pile side is 0.5, and that between pile tip and pile-tip soil is considered to approach 1 for the slender pile, and assumed to be 0.5 for rigid pile. The friction angle is taken as a smaller value. In this paper, the friction angle is taken as 10° and the results obtained from the 3D FEA model for different working conditions are fitted to the data from Georgiadis and Georgiadis, as shown in Figure 9. In general, the fitting is well for any working conditions, proving that the modeling method and the selection of parameters in this paper are reasonable and correct.

3.2. Verification of the Calculation Method

In this section, the proposed method is validated by comparing the load-displacement curve of the pile head and the ultimate soil resistance with the 3D FE analysis results of three cases, respectively. The three cases are shown in Table 1. The modeling method and correlation parameters are the same as those in Section 3.1. Figures 10 and 11 show the pile and soil three-dimensional model of case 2 (; ; and  m).

Case 1 considers two upper slope heights ( m, 3 m) as variates under the same slope angle (; ). The pile head displacement of Case 1 is shown in Figure 12. The theoretical calculation has a good agreement with the result of the 3D FE analysis. And Figure 13 shows the load-displacement curve of the pile head of Case 2. In Case 2, two upper slope angles (, ), one lower slope angle (), and one upper slope height ( m) are considered. The verification result is good.

Case 3 carries out a verification between the theoretical ultimate soil resistance and the results of 3D FE analysis under different angles and heights. A higher degree polynomial is used to accurately fit the shearing force curve of the pile which is obtained from the three-dimensional model [26, 31]. The curve of soil resistance versus pile depth under different lateral load is obtained by differentiating the fitted shearing force curve of the pile. Through combining - curves and the pile displacement curves under different lateral loads , the ultimate soil resistance of each point of the pile is obtained. The verification result is displayed in Figure 14. The theoretical result is in good agreement with the data of 3D FE analysis.

4. Parameter Analysis

4.1. The Effect of Upper Slope Angle

In this section, the effect of on the pile under lateral load is discussed for the case of , and the height of upper slope  m and 2 m by selecting 5 values of , 50°, 40°, 30°, 20° as variables. And calculation parameters are pile length , diameter , elastic modulus of pile  kPa, the undrained shear strength of the soil  kPa, and the elastic modulus  MPa.

Figures 15(a) and 15(b) show the influence of different on the load-displacement curve for the case of  m and 2 m, respectively. Under the same lateral load, the deflection of the pile head increases with the increment of . And the deflection growth rate of the pile head is greater when the lateral load becomes larger. In addition, with the augmentation of, the influence of the upper slope is enhanced, and the dispersion of load-displacement curves under different slope conditions becomes larger. For example, considering that goes up from 0° to 60°, the pile head deflection increases by 34.3% under the condition that the lateral load is 1500 kN and  m. And the maximum growth rate can reach 51% for the case of  m.

Figure 16 investigates the curve of maximum bending moment varying with the pile head displacement under the effect of . Compared with the load-displacement curve, has a similar but smaller influence on the maximum bending moment of the pile. The maximum bending moment of the pile becomes larger when the angle decreases under the same displacement of the pile head. And the dispersion of curves is pronounced for larger . At a pile head displacement , the maximum bending moment of the pile for is higher than that for by 4.4% and 7.2% for the  m and  m, respectively.

4.2. The Effect of Lower Slope Angle

In order to investigate the influence of on the laterally loaded pile, , 15°, 30°, 45°, and 60° are selected as variables in this section. When  m, 2 m, and is 60°, the variation rules of pile head deflection and maximum bending moment are analyzed. The calculation parameters of pile and soil are the same as those in Section 4.1.

Figures 17(a) and 17(b) show the influence of different on the load-displacement curve under the different upper slope heights , respectively. The figures represent that as the angle increases, the displacement of the pile head under the same load raises nonlinearly. And the growth rate of deflection is positively correlated with . Besides, the influence of the lower slope can be weakened, and the dispersion of load-displacement curves under different slope conditions becomes smaller when the upper slope height increases. Considering that goes up from 0° to 60°, the pile head deflection increases by 40% under the condition that the lateral load is 1500 kN and  m. When is 2 m, the maximum growth rate is 20%.

Figure 18 shows the effect of in the maximum bending moment. The maximum bending moment of the pile becomes larger when the angle decreases under the same displacement of the pile head. And the dispersion of curves is smaller for larger . At a pile head displacement , the maximum bending moment of the pile for is higher than that for by 11.6% and 7.1% for the  m and  m, respectively.

4.3. The Effect of the Normalized Height

To study the influence of the normalized height , this section considers as 0, 0.5, 1.0, 2.0, 4.0, and 5.0. And upper slope angle , lower slope angle . Other calculation parameters of pile and soil are the same as in Section 4.1.

Figure 19 shows the load-displacement curves of the pile head considering different sloping conditions. It indicates that the augmentation in increases the deflection of the pile head. When the normalized height changes from 0 to 2, the pile head deflection increases rapidly. For the lateral load , the growth rate of pile head deflection is 32%. When the normalized height is greater than 2, the growth rate of deflection gradually becomes flat. For and applied load , the deflection of the pile head is only 10% higher than that for  m. It also can be seen that the deflection of the pile head is similar to that for the case of when the normalized height exceeds 5 in Figure 19.

The influence of the normalized height on the maximum bending moment of the pile is illustrated in Figure 20. The maximum bending moment is smaller for larger under the same deflection of the pile head. The max of a pile decreases steeply first and then gently with the increase of the normalized height . When the value of exceeds 5, the normalized height has no effect on the maximum bending moment. However, the effect of the increase in on max is quite small. When the displacement of the pile head , the growth rate of maximum bending moment of the pile is only 9.2% with the normalized height range 0 to ∞.

5. Conclusion

This paper has proposed the - curve which is suitable for a concave slope, and the lateral response of a pile has been studied. Equations of initial stiffness varying with depth were obtained through the reduction coefficient method and the normalization method. The nonlinear formulas of the ultimate soil resistance with depth were obtained by using the soil damaged model in front of piles. The rationality of the theory in this paper was verified by comparing the result of three-dimensional finite element analysis. The upper slope angle , the lower slope angle , and the normalized height were discussed. The following conclusions can be drawn: (1)Both ultimate soil resistance and deflection of pile head were predicted using a new method in this paper, which is in good agreement with those calculated by 3D FE analysis(2)The slope angle (, ) has a significant effect on the pile head deflection and a moderate effect on the maximum bending moment. Both deflection and maximum bending moment increase with the increase in slope angle. In addition, the effects of the angle and the height of the upper slope are mutual. As increases, the influence of the lower slope angle weakens and that of the upper slope angle enhances(3)The normalized height is also a remarkable factor for the concave slope. The increase of the deflection and maximum bending moment is greater for larger under the same load. And the influence scope of is 0 to 5(4)The response of lateral load piles on concave slopes differs markedly from the response of lateral load piles on single-angle slopes. There is a large error using the existing - curves to predict the response of the laterally loaded pile on the concave slope. Thus, the proposed - curve of a concave slope in this paper has practical value

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 there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is supported by National Natural Science Foundation of China (Grant Nos. 51678570 and 51978665).