#### Abstract

A simplified nonlinear method was proposed to evaluate lateral behavior of a pile located in or nearby a slope, based on the traditional *p*-*y* method. This method was validated with field test results of a steel pipe pile in clay and model tests of piles in sand slopes. The comparison indicated that the calculated horizontal displacement and bending moment of piles agree well with experimental results. Then, parametric studies were performed, and it shows that horizontal displacement, rotation, bending moment, and shear force increase along with increasing slope angles; the depth of maximum moment locates at about 1.6 *D* below ground surface for horizontal ground, while this value turns to be about 3.6 *D* and 5.6 *D* for sloping ground of 30° and 60°, respectively. The study clearly shows that slope angle has a significant effect on the deflection and lateral capacity of piles.

#### 1. Introduction

The urban infrastructure development in China increases the possibility to construct piles in or nearby slopes, to support bridges, high-rise buildings, transmission towers, off-shore structures, retaining walls, etc. [1–5]. The lateral bearing behavior of these piles is extremely complicated due to sloping ground [6–8]. Compared to conventional piles, they may undergo severe reduction in horizontal bearing behavior. Although conventional laterally loaded piles have been studied by many researchers [9–15], limited literature can be found for piles in sloping ground [16–23]. In the past, the lateral behavior of piles could be evaluated with assumed earth pressure distribution, which can also be determined from field or model tests. This method usually assumes linearly increasing subgrade reaction modulus and cannot account for the influence of sloping ground. Thus, further studies have to be carried out to fill this gap.

This paper presented a simplified *p*-*y* method of laterally loaded piles located in or nearby slopes. A field test result of a steel pipe pile in clay and small-scale pile tests in sand slopes were employed to verify the proposed method and to assess the influence of slope angles on the maximum deflection and required length of a pile.

#### 2. The *p*-*y* Analysis for Sloping Ground

For a laterally loaded pile, the *p*-*y* method is a simple and practical way to account for the relationship between soil reaction *p* and pile deflection *y* along pile shaft. The *p*-*y* curve can be measured in field, by loading cells installed on pile-soil interface, or stress meters installed on steel cages as proposed by McClelland and Focht [24].

##### 2.1. The *p*-*y* Curve of Clay Slopes

According to field test results, a *p*-*y* curve of Houston soft clay was proposed by Matlock, as shown in Figure 1, and has been adopted by the American Petroleum Institute (API). It can be expressed as follows [9, 10, 25]:where is the soil reaction; is the pile deflection; is the ultimate soil resistance; is the ratio of soil resistance; is the ratio of pile deflection; and is the pile deflection when soil resistance reaches 50% of its ultimate value, and it can be evaluated by the following equation:where is the strain when the soil resistance reaches a half of its ultimate value and is the pile diameter.

When a laterally loaded pile locates in a clay slope of a slope angle , the ultimate soil resistance in front of a pile can be computed by the following equation [24, 26]:where is the average undrained shear strength; is the average unit weight of soil; is the depth from ground surface to a studied point on the pile; is the slope angle; and can be computed by the following equation:

##### 2.2. The *p*-*y* Curve of Sand Slopes

A hyperbolic model was proposed to best-fit normalized *p*-*y* curves for laterally loaded piles in sand ground, as shown in Figure 2. This *p*-*y* curve is featured with an initial stiffness and the relevant equation is as follows [27–32]:

where is the initial stiffness, which depends on the soil stiffness, the pile stiffness, and the pile diameter. The initial stiffness, , can be assumed to increase linearly with depth in sands as [33–35]:where is the coefficient of horizontal subgrade reaction, which is related to the internal friction angle, the relative density, etc. can be determined according to internal angle and relative density of sands [26].

The ultimate soil resistance in front of a pile located in sand slope can be described by the following equation [24, 26]:wherewhere is the coefficient of static earth pressure; is the internal friction angle; ; and is the angle of the wedge. Bowman [36] suggested that for loose sand and for dense sand; is the coefficient of active earth pressure.

#### 3. Simplified *p-y* Method of Laterally Loaded Pile in Sloping Ground

##### 3.1. Basic Equations

Assuming that the slope is stable and ignoring friction on the pile, a simplified method for laterally loaded piles in sloping ground can be established as shown in Figure 3. This yields to the differential equation as follows [20, 26]:where is the flexural stiffness of a pile and is the soil resistance of sloping ground.

##### 3.2. Boundary Conditions

The boundary condition at pile top can be free, hinged, and partially or fully fixed, while that at pile toe can be fixed or hinged [20, 26]. In this paper, a bending moment and a shear force are considered as external load on the pile head, and the pile toe is fixed, which yields to the following boundary conditions:

Free pile top:

Fixed pile toe:

##### 3.3. Finite Difference Solution

Subdividing the pile into *N* sections, the length of each section is *h* = *l*/*N*, as shown in Figure 4. According to the principle of the central difference method, two virtual nodes are added at the pile head and toe, respectively. Thus, there are *N* + 1 nodes on the pile shaft (node number: from 0 to *N*), 2 virtual nodes at the top (node number: −2 and −1), and another 2 virtual nodes at the toe (node number: *N* + 1 and *N* + 2).

Let the horizontal displacement at each node be (where ), then (9) can be rewritten as

The slope, , the moment, , and the shear force, , along the pile shaft can be obtained by using the difference method:

The boundary conditions of the pile, namely, (10) and (11), can also be rewritten as follows:

Then, we can establish the basic equation of the pile aswhere is the matrix of load at pile nodes and is the horizontal stiffness matrix. These two matrices are as follows:where

Then, (16) can be rewritten as

The horizontal displacement, , slope, , bending moment, , and shear, , at each pile node can be solved by (19) by using an iteration process as shown in Figure 5. The iteration criterion is indicated as *ε* in the process.

#### 4. Case Studies

##### 4.1. Field Test in Clay

The field test result of a steel pipe pile in clay by Matlock [9] was used to verify the proposed method. Parameters adopted in this analysis are as follows: the diameter, *d* = 0.324 m; the wall thickness of the pipe, *t* = 12.7 mm; the bending rigidity, EI = 31.28 MN·m^{2}; the length, *l* = 12.81 m; the average unit weight of soils, *γ* = 18 kN/m^{3}; the undrained strength of soils, *c*_{u} = 39.1 kPa; *ε*_{50} = 0.012; and the slope angle, *θ* = 0°. The comparison between the calculated and measured pile deflection and the bending moment is plotted in Figure 6.

**(a)**

**(b)**

It is clear in Figure 6 that the calculated horizontal displacement and bending moment agree well with the measured and API method.

##### 4.2. Model Tests of Piles in Sand Slopes

The objective of the model tests is to verify the proposed method. The model piles were made by the PPR (polypropylene random) pipe, which is of 63 mm in outside diameter, 58 mm in inner diameter, and 1680 MPa in elastic modulus. The total length of piles is 1200 mm, and the embedded depth is 900 mm, as shown in Figure 7. The model slope was filled by sand using the rainfall method. The physical and mechanical properties of the sand were tested in laboratory, including the moisture content, = 2.43%; the internal friction angle, *ϕ* = 39°; the unit weight of sand, *γ* = 15.65 kN/m^{3}; the relative density, *D*_{r} = 0.81; and the coefficient of horizontal subgrade reaction, *n*_{h} = 70 MN/m^{3}. The laboratory test carried out and the details of the model test preparation have been presented in the reference [37]. The predicted horizontal displacement and bending moment of piles are compared to the measured in Figure 8 and Table 1.

**(a)**

**(b)**

**(a)**

**(b)**

**(c)**

**(d)**

We can learn from Figure 8 and Table 1 that the pile head deflection and bending moment predicted by the proposed method agree well with the measured in the model tests, and the discrepancy is 2.7% for bending moment and 15.9% for pile head deflection.

As we can learn from Figure 8, the depth of the maximum moment increases from about 10 cm (1.6 *D*, *D* is the pile diameter) below ground level in even ground (*θ* = 0°) to 22.5 cm (3.6 *D*) and 36 cm (5.6 *D*) below ground level in sloping ground of 30° and 60° in slope angle, respectively.

The results also show that the pile head deflection on slope surface rises from 3.3 mm in even ground (*θ* = 0°) to 7.5 mm (127%) and 15.3 mm (364%) in sloping ground of 30° and 60° in angle, respectively; the rotation at the top of the pile rises from −1.15° in even ground (*θ* = 0°) to −1.84° (*θ* = 30°) and −2.74° (*θ* = 60°) in sloping ground; the maximum bending moment increases from 23.5 N·m in even ground (*θ* = 0°) to 31.6 N·m (*θ* = 30°) and 32.2 N·m (*θ* = 60°) in sloping ground, which are 34% and 37% increase, respectively; the maximum shear force increases from −78 N in even ground to −83 N (*θ* = 30°) and −93 N (*θ* = 60°) in sloping ground.

#### 5. Conclusion

A simplified *p*-*y* method of piles located in slopes was proposed and solved using difference method in this paper. The proposed method was verified by the field tests of a steel pipe pile in clay and the model tests of piles in sand slopes. The main influence factor, namely, the slope angle, was discussed by parametric study. The results indicate that the horizontal displacement, rotation, bending moment, and shear force increase with increasing slope angle; the depth of maximum moment is about 1.6 *D* below ground level for even ground and about 3.6 *D* and 5.6 *D* for sloping ground of 30° and 60°, respectively. It is suggested that steep slope should be avoided when designing a laterally loaded pile in sloping ground.

#### 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

The research described in this paper was financially supported by the National Natural Science Foundation of China (Grant nos. 51478051 and 51408066) and Postdoctoral Science Foundation of China (Grant no. 2017M612544).