# Buckling Stability Analysis for Piles in the Slope Foundation Based on Cusp Catastrophe Theory

Authors: Rui Zhang; Chuanxun Li; Dandan Jin
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The aim of this paper is to analyze the buckling stability problem for piles in the slope foundation based on cusp catastrophe theory. Formulation of critical buckling load of piles in the slope foundation is obtained. The influential factors of slope angle, distribution of landslide thrust behind the pile, pile-embedded ratio, pile constraints, pile-side friction, pile-side soil resistance, and pile socketed ratio upon buckling stability characteristic for piles in the slope foundation are examined. The results reveal that when pile diameter remains unchanged, critical buckling load increases with the increase of pile length when pile-embedded ratio reaches 60%. When pile length remains unchanged, critical buckling load increases with the increase of pile diameter. Critical buckling load with the assumption of nonlinear horizontal elastic resistance of pile-side soil in the paper is more close to the value based on horizontal elastic resistance of pile-side soil suggested in the code. When slope angle increases, decreased extent of buckling critical load for piles 30–60 m in length is more obvious than the piles which are 10–30 m in length. Strengthening of pile constraints and increase of pile-embedded ratio and socketed ratio are helpful to pile critical buckling load increase. The influential factors of pile-side friction and landslide thrust behind the pile upon pile critical buckling load are tiny and can be neglected.


---

## Body

## 1. Introduction

During the 14th Five-Year Plan Period, the pace of construction and upgrading will be promoted continuously in China. Meanwhile, as a mountainous country, mountainous and hilly areas cover two-thirds of the total land area. Therefore, with the continuous impetus of national infrastructure, as well as increasing piles of civil buildings, roads and bridges will be located in slope foundation. Piles in slope foundation are influenced by asymmetric pile-soil interaction, slope sliding force, and nonlinear pile-side soil resistance, whose buckling stability research is very vital. It is of important theoretical and practical significance to make deep research on buckling stability for piles in slope foundation [1].Extensive studies with different approaches have been performed by many scholars on piles buckling stability analysis. Based on three-dimensional finite element method, Kurian and Srinivas [2] discussed bearing capacity for tapered piles and found that ultimate load-carrying capacities of tapered piles were higher than those of uniform piles by 10%. Budkowska and Szymczak [3] utilized energy method to analyze buckling stability for partially embedded pile and revealed that the skin friction could be neglected if the ratio of the embedded length of the pile to its total length was less than 0.6. With dimensionless analytical method, Heelis et al. [4] made research into buckling stability for fully and partially embedded pile. Lee et al. [5] designed model experiment to estimate axial bearing capacity for tapered piles. The proposed method was valid for piles in sands of medium to dense conditions. Based on variation theory, Ren and Wang [6] discussed vertical bearing capacity for cement soil reinforced piles and predicted the performance of single JPP pile. Zou and Zhao [7] employed difference method to analyze axial loading performance for superlong piles on bedrock in soft soil and got the conclusion that when l/d reaches a certain value, regardless of whether the pile toe is socketed in rock or not, the tip resistance will not be exerted. With cavity expansion theory and model experiment, Manandhar and Yasufuku [8] discussed vertical bearing capacity for tapered piles in sand and discovered that the pile load-settlement curves can easily be determined with the aid of the simple fundamental properties of soils with the closed-form solution. By means of numerical analysis method, Kong et al. [9] made research of bearing capacity for belled wedge pile and drew the conclusion that belled wedge pile can improve vertical load-carrying capacities under both compression and uplift loads and reduce the amount of negative skin friction developed. Bhattacharya and Goda [10] formulated probabilistic method to analyze buckling stability characteristic for axially loaded piles in liquefiable soils and found that boundary condition of the pile was essential to pile buckling stability. Nadeem et al. [11] adopted finite element method on the research of nonlinear stability for slender defect piles. The analyses’ results showed that the load-carrying capacity of the piles decreased with increasing pile curvature and increasing slenderness ratio of the initially bent piles. With the help of energy method, Zhao et al. [12] carried out buckling stability research on bridge piles in steep slopes and put forward the following conclusions: (1) When pile-embedded ratio was <0.61, buckling point was in the free section of the pile. (2) Pile critical buckling load increased with increasing pile socketed ratio. When socketed ratio reached a certain value, its increasing effect to critical buckling load was exerted. Yin et al. [13] employed similarity theory in the indoor model experiment to discuss vertical bearing capacity and buckling stability for piles in slope foundation and concluded that the bigger the slope angle or free section of the pile, the more decrease of critical buckling load. Lu and Zhao. [14] adopted energy method to analyze calculated length for buckling stability of steel pipe pile and obtained the equations on calculated length of piles in uniform foundation and layered foundation. Lee et al. [15] applied Runge–Kutta method and Regula-Falsi method on the buckling stability research for tapered piles. The results showed that the proposed model provided good predictions of pile buckling response and can determine an optimal taper ratio such that the pile is capable of carrying a maximum flexural buckling load.Cusp catastrophe theory [16] is a branch on discontinuity, which has developed greatly during the last 30 years and had a wide range of applications in pile ultimate bearing capacity prediction [17]. Some scholars have studied the buckling stability of piles based on cusp catastrophe theory. With static load test, Zhang and Cui [18] derived computing formula on ultimate vertical bearing capacity for a single end bearing pile and found that pile displacement value of instability point could be calculated out with parabolic load-displacement curve. Chen et al. [19] discussed buckling stability for Y-shape piles and established cusp catastrophe model for piles whose tips were socketed and tops were free; critical buckling loads of Y-shape piles were higher than those of uniform piles by 30%. Li et al. [20] performed buckling stability analysis for tapered piles considering self-weight influence and concluded that self-weight of piles had a great influence on critical buckling load for tapered piles. When embedded ratio was less than 40% of pile length, larger critical buckling load would be obtained, which would cause unsafe bearing capacity designed according to the strength of pile body.Based on the analysis of the existing literature, there are still some challenges on buckling stability analysis. (1) Cusp catastrophe theory is seldom adopted in the existing articles for in-slope pile buckling stability analysis. (2) Potential energy function is not complete, which may lead to incomprehensive influential factors discussion on pile buckling stability analysis.The objective of this study is to analyze buckling stability for piles in slope foundation based on cusp catastrophe theory. The paper establishes total potential energy equation with assumptions of nonlinear distributed horizontal elastic resistance of pile-side soil and landslide thrust behind pile. Furthermore, the paper discusses influences of slope angle, landslide thrust distribution behind pile, pile-embedded ratio, pile constraints, pile-side soil friction, pile-side soil resistance, and pile socketed ration upon buckling stability for piles in slope foundation.


## 2. Cusp Catastrophe Theory

Potential function of cusp catastrophe model is defined as [16](1)Vx=14x4+12ux2+vx,where x is the system state variable and u and v are system control variables.Making the derivation of potential functionV (x) zero, we can get equilibrium surface of cusp catastrophe model:(2)V′x=x3+ux+v=0.This surface is called catastrophe manifold, which indicates condition of system equilibrium state. The corresponding point is called critical or extreme point in mathematics.System singular point sets satisfy equation (2) and the second derivative of potential function V″x=0; that is,(3)V″x=3x2+u=0.Eliminatex through the combination of equations (2) and (3), and parametric plane equation is obtained as(4)D=4u3+27v2=0.Equation (4) is called system bifurcation sets, which denotes projection on the plane determined by control variables u and v of singular point sets. When system control variables satisfy equation (4), the system is in critical equilibrium state and will jump to a stable equilibrium state finally to achieve system catastrophe.CurveD is called bifurcation sets, whose discriminant is written as D=4u3+27v2.WhenD > 0, the system is in a stable state. When D < 0, the system is in an unstable state. Catastrophe occurs when u < 0 and D = 0. In this state, there are three real roots for equation (4).


## 3. Computation Model and Basic Assumptions

Assuming that L is the total pile length, it can be divided into three sections: socketed sectionl1, slope influenced section l2, and free section l; that is, L = l1 + l2 + l. Vertical load on pile top is P. Pile embedded depth is h: h = l1 + l2. The coordinate origin o is at the center of pile bottom. Simplified computation model for buckling stability analysis of an in-slope pile is shown in Figure 1.Figure 1
Simplified computation model for buckling stability analysis of an in-slope pile.Pile bottom is a fixed constraint, pile top is free,α is the slope angle, d is the pile diameter, Q(x) is the landslide thrust behind the pile, l2 = 4dtanα [13], and τ is the ultimate friction strength between pile and soil.Basic assumptions in the computation are listed as follows:(1)
Ultimate friction strength between pile and soilτ is a constant.(2)
Horizontal elastic resistance of pile-side soilq (x, y) is nonlinearly distributed within the range of embedded depth of pile, whose expression is regarded as follows [13]:(5)qx,y=mh−l2−xb0y0≤x≤l1,mh−x2l2b0yl1≤x≤h,where m is proportional coefficient of horizontal elastic resistance of pile-side soil. b0 is calculated width of piles. For circular piles, b0 = 0.9 (1.5 d + 0.5), when d < 1.0 m. b0 = 0.9 (d + 1), when d ≥ 1.0 m [20]. Generalized expression of landslide thrust behind pile Q (x) is regarded as [21](6)Qx=a1x2+b1x+c1,where a1, b1, and c1 are undetermined parameters related to the resultant force of landslide thrust F.


## 4. Solution of Critical Buckling Load

### 4.1. Establishment of Deflection Functions

According to the boundary conditions shown in Figure1, the deflection function for a pile whose tip is socketed and top is free is established as [12, 20](7)y=∑n=1∞Cn1−cos2n−12Lπx,where n is half-wave number of the deflection function. Cn is undetermined constant term. Under the small deformation assumption of pile-soil system, take n = 1. Therefore, pile deflection function is expressed as(8)y=C1−cosπx2L.From equation (8), the first and second derivatives of pile deflection function are obtained as(9)y′=πC2Lsinπx2L,(10)y″=π2C4L2cosπx2L.


### 4.2. Solution Based on Cusp Catastrophe Theory

The total potential energy equation of pile-soil system is established as follows [12, 20]:(11)Π=UP+US+VQ+VP+VG+Wf,where UP is the bending strain energy of the pile, US is the elastic deformation energy of soil beside pile, VQ is the potential energy of landslide thrust behind the pile, VP is the potential energy of pile top load, VG is the potential energy of pile self-weight, and Wf is the potential energy of pile-side friction embedded in soil. Substituting equations (8)∼(10), deduction of each potential energy item is listed as follows.Bending strain energy of pileUP is as follows:(12)UP=EI2∫0Ld2y/dx221+dy/dx23dx=EI2∫0Ld2ydx221−3dydx2dx=π4EIC232L4∫0Lcos2πx2L1−3π2C24L2sin2πx2Ldx=π4EI64L3C2−3π6EI1024L5C4,where E is the elastic modulus of the pile and I is the moment of inertia of the pile section.Deformation energy of soil beside pileUS is as follows:(13)US=12∫0hqydx=mb02∫0l1h−l2−xC21−cosπx2L2dx+mb02l2∫l1hh−x2C21−cosπx2L2dx=−mb0C28π230L2−3π2l12+2L2−16cosπl12L+cosπl1L−mb0C24π3l2−π3h3+π3l13hl2+l12+32πL2l2cosπl12L−116cosπl1L−64L3sinπh2L−132sinπhL−8L−8L2+π2l22sinπl12L+L−2L2+π2l22sinπl1L=−mb0C28π2AS1+32L2AS2−128L3πl2AS3−16πLl2AS4=−mb0AS8π2C2,where(14)AS1=30L2−3π2l12−2π2h3l2+2π2l13h+l12l2,AS2=cosπl12L−116cosπl1L,AS3=sinπh2L−132sinπhL−sinπl12L−132sinπl1L,AS4=sinπl12L−18sinπl1L.Potential energy of landslide thrust behind pileVQ is as follows:(15)VQ=−12∫l1ha1x2+b1x+c1l2y2dx=−12∫l1hC2l2a1x2+b1x+c11−cosπx2L2dx=−C28π3l2π3h6c1+h3b1+2a1h−π3l16c1+l13b1+2a1l1−32πL2b1+2a1hcosπh2L−116cosπhL−16L−8a1L2+π2a1h2+b1h+c1sinπh2L+2L−2a1L2+π2a1h2+b1h+c1sinπhL+32πL2b1+2a1l1cosπl12L−116cosπl1L+16L−8a1L2+π2a1l12+b1l1+c1sinπl12L−2L−2a1L2+π2a1l12+b1l1+c1sinπl1L=−C28π3l2BQ1+BQ2cosπh2L−116cosπhL+BQ3cosπl12L−116cosπl1L+BQ4sinπh2L−18sinπhL+BQ5sinπl12L−18sinπl1L+128a1L3BQ6=−BQ8π3l2C2,where(16)BQ1=π3h6c1+h3b1+2a1h−π3l16c1+l13b1+2a1l1,BQ2=−32πL2b1+2a1h,BQ3=32πL2b1+2a1l1,BQ4=−16π2La1h2+b1h+c1,BQ5=16π2La1l12+b1l1+c1,BQ6=sinπh2L−132sinπhL−sinπl12L−132sinπl1L.Potential energy of pile top loadVP is as follows:(17)VP=−P∫0L1+dydx2−1dx=−P2∫0Ldydx2−14dydx4dx=−P2∫0Lπ2C24L2sin2πx2L−π4C464L4sin4πx2Ldx=−π2P16LC2+3π4P1024L3C4.Potential energy of pile self-weightVG is as follows:(18)VG=−γA∫0LL−x1+dydx2−1dx=−γA2∫0LL−xdydx2−14dydx4dx=−γA2∫0LL−xπ2C24L2sin2πx2L−π4C464L4sin4πx2Ldx=−γA32π2−4C2+π2γA3π2−162048L2C4,where γ is the pile unit weight and A is the pile section area.Potential energy of pile-side frictionWf is as follows:(19)Wf=τU∫0hh−x1+dydx2−1dx=τU2∫0hh−xdydx2−14dydx4dx=τU2∫0hh−xπ2C24L2sin2πx2L−π4C464L4sin4πx2Ldx=τUC232L2−2L2+π2h2+2L2cosπhL+π2τUC44096L415L2−6π2h2+L2−16cosπhL+cos2πhL=τUCf32L2C2+π2τUC¯f4096L4C4,where U is the pile section circumference. The expressions of Cf and C—f are written as(20)Cf=π2h2−2L21−cosπhL,C¯f=15L2−6π2h2−16L2cosπhL−116cos2πhL.Substituting equations (12)–(19) into equation (11), the total potential energy equation of pile-soil system is expressed as(21)∏=a¯0+a¯1C+a¯2C2+a¯3C3+a¯4C4.From equations (12)–(19), we can see that each expression of partial potential energy is the function of C2 and C4, so a¯0=a¯1=a¯3=0. The expressions of a¯2anda¯4 can be written as(22)a¯2=π4EI64L3−mb0AS8π2−BQ8π3l2−π2P16L−γA32π2−4+τUCf32L2,a¯4=−3π6EI1024L5+3π4P1024L3+π2γA3π2−162048L2+π2τUC¯f4096L4.Substituting variables ofC, let C=t1/4a¯41/4, a¯4>0, and the total potential energy equation (21) is expressed as(23)∏=14t4+12ut2+vt,where u=a¯2/a¯4 and v=a¯1/4a¯41/4=0.According to cusp catastrophe theory, the discriminant is written as(24)D=4u3+27v2=4a¯2a¯43.WhenD = 0, system catastrophe occurs. Therefore, it is assumed that(25)D=4u3+27v2=4a¯2a¯43=0.From equation (25), the expression of pile critical buckling load Pcr is obtained as(26)Pcr=π2EI4L2−2mb0LASπ4−2LBQπ5l2−γAL2π2π2−4+τUCf2π2L.


## 4.1. Establishment of Deflection Functions

According to the boundary conditions shown in Figure1, the deflection function for a pile whose tip is socketed and top is free is established as [12, 20](7)y=∑n=1∞Cn1−cos2n−12Lπx,where n is half-wave number of the deflection function. Cn is undetermined constant term. Under the small deformation assumption of pile-soil system, take n = 1. Therefore, pile deflection function is expressed as(8)y=C1−cosπx2L.From equation (8), the first and second derivatives of pile deflection function are obtained as(9)y′=πC2Lsinπx2L,(10)y″=π2C4L2cosπx2L.


## 4.2. Solution Based on Cusp Catastrophe Theory

The total potential energy equation of pile-soil system is established as follows [12, 20]:(11)Π=UP+US+VQ+VP+VG+Wf,where UP is the bending strain energy of the pile, US is the elastic deformation energy of soil beside pile, VQ is the potential energy of landslide thrust behind the pile, VP is the potential energy of pile top load, VG is the potential energy of pile self-weight, and Wf is the potential energy of pile-side friction embedded in soil. Substituting equations (8)∼(10), deduction of each potential energy item is listed as follows.Bending strain energy of pileUP is as follows:(12)UP=EI2∫0Ld2y/dx221+dy/dx23dx=EI2∫0Ld2ydx221−3dydx2dx=π4EIC232L4∫0Lcos2πx2L1−3π2C24L2sin2πx2Ldx=π4EI64L3C2−3π6EI1024L5C4,where E is the elastic modulus of the pile and I is the moment of inertia of the pile section.Deformation energy of soil beside pileUS is as follows:(13)US=12∫0hqydx=mb02∫0l1h−l2−xC21−cosπx2L2dx+mb02l2∫l1hh−x2C21−cosπx2L2dx=−mb0C28π230L2−3π2l12+2L2−16cosπl12L+cosπl1L−mb0C24π3l2−π3h3+π3l13hl2+l12+32πL2l2cosπl12L−116cosπl1L−64L3sinπh2L−132sinπhL−8L−8L2+π2l22sinπl12L+L−2L2+π2l22sinπl1L=−mb0C28π2AS1+32L2AS2−128L3πl2AS3−16πLl2AS4=−mb0AS8π2C2,where(14)AS1=30L2−3π2l12−2π2h3l2+2π2l13h+l12l2,AS2=cosπl12L−116cosπl1L,AS3=sinπh2L−132sinπhL−sinπl12L−132sinπl1L,AS4=sinπl12L−18sinπl1L.Potential energy of landslide thrust behind pileVQ is as follows:(15)VQ=−12∫l1ha1x2+b1x+c1l2y2dx=−12∫l1hC2l2a1x2+b1x+c11−cosπx2L2dx=−C28π3l2π3h6c1+h3b1+2a1h−π3l16c1+l13b1+2a1l1−32πL2b1+2a1hcosπh2L−116cosπhL−16L−8a1L2+π2a1h2+b1h+c1sinπh2L+2L−2a1L2+π2a1h2+b1h+c1sinπhL+32πL2b1+2a1l1cosπl12L−116cosπl1L+16L−8a1L2+π2a1l12+b1l1+c1sinπl12L−2L−2a1L2+π2a1l12+b1l1+c1sinπl1L=−C28π3l2BQ1+BQ2cosπh2L−116cosπhL+BQ3cosπl12L−116cosπl1L+BQ4sinπh2L−18sinπhL+BQ5sinπl12L−18sinπl1L+128a1L3BQ6=−BQ8π3l2C2,where(16)BQ1=π3h6c1+h3b1+2a1h−π3l16c1+l13b1+2a1l1,BQ2=−32πL2b1+2a1h,BQ3=32πL2b1+2a1l1,BQ4=−16π2La1h2+b1h+c1,BQ5=16π2La1l12+b1l1+c1,BQ6=sinπh2L−132sinπhL−sinπl12L−132sinπl1L.Potential energy of pile top loadVP is as follows:(17)VP=−P∫0L1+dydx2−1dx=−P2∫0Ldydx2−14dydx4dx=−P2∫0Lπ2C24L2sin2πx2L−π4C464L4sin4πx2Ldx=−π2P16LC2+3π4P1024L3C4.Potential energy of pile self-weightVG is as follows:(18)VG=−γA∫0LL−x1+dydx2−1dx=−γA2∫0LL−xdydx2−14dydx4dx=−γA2∫0LL−xπ2C24L2sin2πx2L−π4C464L4sin4πx2Ldx=−γA32π2−4C2+π2γA3π2−162048L2C4,where γ is the pile unit weight and A is the pile section area.Potential energy of pile-side frictionWf is as follows:(19)Wf=τU∫0hh−x1+dydx2−1dx=τU2∫0hh−xdydx2−14dydx4dx=τU2∫0hh−xπ2C24L2sin2πx2L−π4C464L4sin4πx2Ldx=τUC232L2−2L2+π2h2+2L2cosπhL+π2τUC44096L415L2−6π2h2+L2−16cosπhL+cos2πhL=τUCf32L2C2+π2τUC¯f4096L4C4,where U is the pile section circumference. The expressions of Cf and C—f are written as(20)Cf=π2h2−2L21−cosπhL,C¯f=15L2−6π2h2−16L2cosπhL−116cos2πhL.Substituting equations (12)–(19) into equation (11), the total potential energy equation of pile-soil system is expressed as(21)∏=a¯0+a¯1C+a¯2C2+a¯3C3+a¯4C4.From equations (12)–(19), we can see that each expression of partial potential energy is the function of C2 and C4, so a¯0=a¯1=a¯3=0. The expressions of a¯2anda¯4 can be written as(22)a¯2=π4EI64L3−mb0AS8π2−BQ8π3l2−π2P16L−γA32π2−4+τUCf32L2,a¯4=−3π6EI1024L5+3π4P1024L3+π2γA3π2−162048L2+π2τUC¯f4096L4.Substituting variables ofC, let C=t1/4a¯41/4, a¯4>0, and the total potential energy equation (21) is expressed as(23)∏=14t4+12ut2+vt,where u=a¯2/a¯4 and v=a¯1/4a¯41/4=0.According to cusp catastrophe theory, the discriminant is written as(24)D=4u3+27v2=4a¯2a¯43.WhenD = 0, system catastrophe occurs. Therefore, it is assumed that(25)D=4u3+27v2=4a¯2a¯43=0.From equation (25), the expression of pile critical buckling load Pcr is obtained as(26)Pcr=π2EI4L2−2mb0LASπ4−2LBQπ5l2−γAL2π2π2−4+τUCf2π2L.


## 5. Correctness Verification of Critical Buckling Load Expression

According to critical buckling load expression of equal section pile in level ground from [20], suppose that slope angle α = 0, and horizontal elastic resistance of pile-side soil q (x, y) is adopted as q (x, y) = mb0 (h − x) y. In the meantime, the influences of ultimate friction strength between pile and soil τ and landslide thrust behind pile are neglected, and formula (26) can be written as follows:(27)Pcr=π2EI4L2+2mb0L3π4−30+3π2hL2+216cosπh2L−cosπhL−γAL2π2π2−4.The derived critical buckling load expression of equal section pile in level ground from formula (26) is the same as that in [20].Furthermore, considering pile self-weight zero and pile-embedded depthh = 0, expression (27) can be simplified as(28)Pcr=π2EI4L2.According to expression (28), we obtain critical buckling load of ordinary compressed bar whose tip is socketed and top is free. Therefore, critical buckling load expression (26) for piles in slope foundation derived in the paper is reliable.


## 6. Case Study

The pile is in slope foundation, pile concrete strength grade isC25, slope soil is soft plastic clay, and slope angle α is 30°. Pile-embedded ratio h/L is 0.6. Landslide thrust behind pile Q (x) is assumed as triangular distribution. Composition of landslide thrust behind pile is 80 kN/m. Other pile and soil parameters are obtained from the work of Zhao et al. [12] and listed in Table 1. Physical meanings of the input parameters are explained in the previous section.Table 1
Material parameters for the pile and soil.
E (MPa)γ (kN∙m−3)m (MN∙m−4)τ (kPa)2.38 × 104255.020


### 6.1. Variation Analysis of Pile Buckling Load with Different Pile Lengths

Critical buckling load changes for piles in slope foundation with varying pile diameter and length are shown in Figure2. Pile length ranges from 10 m to 60 m. Pile diameter d is 0.5 m, 0.75 m, 1.0 m, 1.5 m, and 2.0 m. Pile-embedded ratio h/L is 0.6.Figure 2
Pcr-L curve with different pile lengths.It can be seen from Figure2 that when pile diameter d remains unchanged, with the increase of pile length L, pile socketed length increases, and critical buckling load of piles in slope foundation increases accordingly. When pile length L remains unchanged, with the increase of pile diameter d, the pile slenderness ratio L/d decreases, and critical buckling load of piles in slope foundation increases accordingly.


### 6.2. Influence of Slope Angle

With slope angleα being 0°, 10°, 20°, 30°, 40°, 50°, and 60°, the influence of slope angle on critical buckling load of piles in slope foundation is discussed. Critical buckling load change curves of piles in slope foundation with varying slope angle are shown in Figure 3.Figure 3
Critical buckling load with varying slope angle.It can be seen from Figure3 that, compared with pile in level ground (α = 0°), in-slope pile critical buckling load decreases. For medium-long piles in the range of 10–20 m, critical buckling load reduction is small when slope angle increases. For long piles in the range of 30–40 m and superlong piles over 50 m, critical buckling load reduction is obvious when slope angle increases. Regardless of the pile length, critical buckling load decreases proportionally with the angle approximately. That is because when pile-embedded ratio is constant, socketed length l1 decreases proportionally with slope angle approximately.


### 6.3. Influence of Landslide Thrust behind the Pile

According to landslide thrust distribution function table in [21], parabola distribution, linear distribution, trapezoidal distribution, triangular distribution, and no thrust distribution of landslide thrust are employed, and the influence of landslide thrust behind pile upon critical buckling load is discussed. Composition of landslide thrust behind pile is 80 kN/m. Critical buckling load change curves of piles in the slope foundation with varying landslide thrust are listed in Table 2.Table 2
Critical buckling load for piles in the slope foundation with varying landslide thrust distribution (MN).
L (m)Parabola distributionLinear distributionTrapezoidal distributionTriangular distributionNo thrust distribution1036.10836.09936.11236.11436.1532068.22168.20868.22768.22968.36330240.959240.946240.965240.968241.20640632.638632.627632.646632.648632.996501336.5871336.5781336.6001336.6021336.991602447.8792447.8772447.8992447.9012448.464It can be seen from Table2 that when landslide thrust behind pile is considered, critical buckling load of piles in slope foundation decreases with tiny reduction. Compared with critical buckling load without thrust distribution, the reduction is less than 0.5%. Critical buckling load varies with tiny reduction with parabola, linear, trapezoidal, and triangular distribution. Therefore, for in-slope pile buckling stability analysis, influence of landslide thrust behind pile upon critical buckling load is tiny and can be neglected.


### 6.4. Influence of Pile-Embedded Ratio

Pile-embedded ratioh/L is 0.15, 0.20, 0.30, 0.40, 0.50, 0.60, 0.70, and 0.75, and the influence of pile-embedded ratio upon critical buckling load is discussed. Critical buckling load change curves of piles in slope foundation with varying pile-embedded ratio are shown in Figure 4.Figure 4
Critical buckling load with varying pile-embedded ratio. (a) Global graph. (b) Partial enlarged graph.
(a)(b)It can be seen from Figure4 that when pile-embedded ratio h/L ≤ 0.20, critical buckling load of piles in slope foundation decreases gradually with the increase of pile length. At this moment, critical buckling point is in the free section of the pile. When pile-embedded ratio h/L is in the range of 0.30 to 0.50, critical buckling load of piles in slope foundation decreases first and increases afterwards with the increase of pile length. Critical buckling point transits from free section to influenced section of piles. When pile-embedded ratio h/L ≥ 0.60, critical buckling load of piles in slope foundation increases with the increase of pile length. Therefore, pile-embedded ratio h/L should reach 60% for piles in slope foundation.


### 6.5. Influence of Pile Constraints

Constraint types of pile tip and top include free, elastic embedded, articulated, and fixed. Suppose that pile tip constraint is fixed and pile top constraint is one of the above four constraints. Combine the constraint conditions of pile tip and pile top afterwards. With the assumption of small deformation of pile-soil system, influence of pile constraints upon critical buckling load of piles in slope foundation is discussed. The corresponding pile deflection functions are listed in Table3 [22].Table 3
Deflection functions under different pile constraints.
Pile top constraintDeflection functionFreey=C1−cosπx/2LElastic embeddedy=C1−cosπx/LArticulatedy=Ccos3πx/2L−cosπx/2LFixedy=C1−cos2πx/LCritical buckling load change curves of piles in slope foundation with varying pile top constraints are shown in Figure5.Figure 5
Critical buckling load with varying pile top constraints.It can be seen from Figure5 that critical buckling load is of the minimum value when pile top is free. Piles are most prone to buckling instability. Critical buckling load is of the maximum value when pile top is fixed. Piles are the least prone to buckling instability. Critical buckling load is between the previous values when pile top is elastic embedded or articulated. Critical buckling load with articulated pile top is a little more than the value with elastic embedded pile top. Therefore, when pile tip is fixed, critical buckling load of piles in slope foundation can be improved effectively with strengthening of pile top constraints.


### 6.6. Influence of Pile-Side Friction

Let ultimate friction strength between pile and soilτ be 0, 30, and 60 kPa. Pile-embedded ratio h/L is 0.6. The influence of pile-side friction upon critical buckling load is discussed. Critical buckling loads with different ultimate friction strength between pile and soil and pile length are shown in Table 4.Table 4
Critical buckling load for piles in the slope foundation with varying limit friction strength (MN).
L (m)τ = 0 kPaτ = 30 kPaτ = 60 kPa1036.08436.12936.1742068.16968.25968.34830240.879241.013241.14740632.531632.710632.888501336.451336.671336.90602447.712447.982448.25It can be seen from Table4 that when pile-side friction is considered, critical buckling load of piles in slope foundation increases. But the extent of increase is less than 0.5%. Therefore, the influence of pile-side friction upon critical buckling load can be neglected.


### 6.7. Influence of Horizontal Elastic Resistance of Pile-Side Soil

The typical methods to calculate horizontal elastic resistance of pile-side soil include Constant method,m method, and c method. Among them, m method is the method adopted in the current national code [23] to calculate foundation counterforce. At present, some scholars combine m method and c method with Constant method to calculate pile-side soil resistance. In this section, critical buckling load of piles in slope foundation is obtained with different pile-side soil resistance computation methods. The expressions assumed in the above methods are listed in Table 5 (the expressions assumed in the paper are given in equation (5)).Table 5
Pile-side soil resistance expressions.
Method nameq(x, y)0 ≤x < l1l1≤x ≤ hm methodmb0h−xyc methodmb0h−x0.5ym-constant methodmb0h−l1ymb0h−xyc-constant methodmb0h−l10.5ymb0h−x0.5yq(x, y) can be expressed as(29)qx,y=kxy,where K (x) is pile-side soil resistance coefficient. Therefore, pile-side soil resistance distribution along pile length is different for the above methods. Pile-side soil resistance is linear distribution along pile length for m method. Pile-side soil resistance is 0.5 power exponential distribution for c method. For m-Constant method and c-Constant method, pile-side soil resistance is constant when 0 ≤ x < l1 and linear or 0.5 power exponential distribution when l1 ≤ x ≤ h. As shown in expression (5), for pile-side soil resistance suggested in the paper, it is linear distribution when 0 ≤ x < l1 and square exponential distribution when l1 ≤ x ≤ h.Critical buckling load change curves of piles with varying pile-side soil resistance computation methods are shown in Figure6.Figure 6
Critical buckling load with varying pile-side soil resistance computation methods.It can be seen from Figure6 that, for medium-long piles in the range of 10 m–20 m, difference of critical buckling load with varying pile-side soil resistance computation methods is not wide. For long piles in the range of 30–40 m and superlong piles over 50 m, difference of critical buckling load with varying pile-side soil resistance computation methods appears gradually. Critical buckling load obtained with pile-side soil resistance computation method suggested in the paper is close to the value with m method adopted in the current national code. Therefore, pile-side soil resistance computation method suggested in the paper is of certain usability.


### 6.8. Influence of Pile Socketed Ratio

For piles in slope foundation with different lengths, the influence of pile socketed ratiol1/L upon critical buckling load is discussed through increase of pile socketed ratio. Critical buckling load change curves of piles in slope foundation with different lengths and pile socketed ratio are shown in Figure 7.Figure 7
Critical buckling load with varying pile socketed ratio.It can be seen from Figure7 that critical buckling load of piles in slope foundation increases with the increase of pile socketed ratio. That is to say, appropriate increase of pile socketed length is helpful to increase pile critical buckling load. Furthermore, the more of pile length, the more of increase latitude of pile critical buckling load. There is especially obvious increase latitude for superlong piles over 50 m. Critical buckling load change curves will become flat when pile socketed ratio reaches a certain value. That is to say, there is a critical fixed depth [12]. Critical buckling load change curves of medium-long piles in the range of 10 m–20 m tend to be flat earlier.


## 6.1. Variation Analysis of Pile Buckling Load with Different Pile Lengths

Critical buckling load changes for piles in slope foundation with varying pile diameter and length are shown in Figure2. Pile length ranges from 10 m to 60 m. Pile diameter d is 0.5 m, 0.75 m, 1.0 m, 1.5 m, and 2.0 m. Pile-embedded ratio h/L is 0.6.Figure 2
Pcr-L curve with different pile lengths.It can be seen from Figure2 that when pile diameter d remains unchanged, with the increase of pile length L, pile socketed length increases, and critical buckling load of piles in slope foundation increases accordingly. When pile length L remains unchanged, with the increase of pile diameter d, the pile slenderness ratio L/d decreases, and critical buckling load of piles in slope foundation increases accordingly.


## 6.2. Influence of Slope Angle

With slope angleα being 0°, 10°, 20°, 30°, 40°, 50°, and 60°, the influence of slope angle on critical buckling load of piles in slope foundation is discussed. Critical buckling load change curves of piles in slope foundation with varying slope angle are shown in Figure 3.Figure 3
Critical buckling load with varying slope angle.It can be seen from Figure3 that, compared with pile in level ground (α = 0°), in-slope pile critical buckling load decreases. For medium-long piles in the range of 10–20 m, critical buckling load reduction is small when slope angle increases. For long piles in the range of 30–40 m and superlong piles over 50 m, critical buckling load reduction is obvious when slope angle increases. Regardless of the pile length, critical buckling load decreases proportionally with the angle approximately. That is because when pile-embedded ratio is constant, socketed length l1 decreases proportionally with slope angle approximately.


## 6.3. Influence of Landslide Thrust behind the Pile

According to landslide thrust distribution function table in [21], parabola distribution, linear distribution, trapezoidal distribution, triangular distribution, and no thrust distribution of landslide thrust are employed, and the influence of landslide thrust behind pile upon critical buckling load is discussed. Composition of landslide thrust behind pile is 80 kN/m. Critical buckling load change curves of piles in the slope foundation with varying landslide thrust are listed in Table 2.Table 2
Critical buckling load for piles in the slope foundation with varying landslide thrust distribution (MN).
L (m)Parabola distributionLinear distributionTrapezoidal distributionTriangular distributionNo thrust distribution1036.10836.09936.11236.11436.1532068.22168.20868.22768.22968.36330240.959240.946240.965240.968241.20640632.638632.627632.646632.648632.996501336.5871336.5781336.6001336.6021336.991602447.8792447.8772447.8992447.9012448.464It can be seen from Table2 that when landslide thrust behind pile is considered, critical buckling load of piles in slope foundation decreases with tiny reduction. Compared with critical buckling load without thrust distribution, the reduction is less than 0.5%. Critical buckling load varies with tiny reduction with parabola, linear, trapezoidal, and triangular distribution. Therefore, for in-slope pile buckling stability analysis, influence of landslide thrust behind pile upon critical buckling load is tiny and can be neglected.


## 6.4. Influence of Pile-Embedded Ratio

Pile-embedded ratioh/L is 0.15, 0.20, 0.30, 0.40, 0.50, 0.60, 0.70, and 0.75, and the influence of pile-embedded ratio upon critical buckling load is discussed. Critical buckling load change curves of piles in slope foundation with varying pile-embedded ratio are shown in Figure 4.Figure 4
Critical buckling load with varying pile-embedded ratio. (a) Global graph. (b) Partial enlarged graph.
(a)(b)It can be seen from Figure4 that when pile-embedded ratio h/L ≤ 0.20, critical buckling load of piles in slope foundation decreases gradually with the increase of pile length. At this moment, critical buckling point is in the free section of the pile. When pile-embedded ratio h/L is in the range of 0.30 to 0.50, critical buckling load of piles in slope foundation decreases first and increases afterwards with the increase of pile length. Critical buckling point transits from free section to influenced section of piles. When pile-embedded ratio h/L ≥ 0.60, critical buckling load of piles in slope foundation increases with the increase of pile length. Therefore, pile-embedded ratio h/L should reach 60% for piles in slope foundation.


## 6.5. Influence of Pile Constraints

Constraint types of pile tip and top include free, elastic embedded, articulated, and fixed. Suppose that pile tip constraint is fixed and pile top constraint is one of the above four constraints. Combine the constraint conditions of pile tip and pile top afterwards. With the assumption of small deformation of pile-soil system, influence of pile constraints upon critical buckling load of piles in slope foundation is discussed. The corresponding pile deflection functions are listed in Table3 [22].Table 3
Deflection functions under different pile constraints.
Pile top constraintDeflection functionFreey=C1−cosπx/2LElastic embeddedy=C1−cosπx/LArticulatedy=Ccos3πx/2L−cosπx/2LFixedy=C1−cos2πx/LCritical buckling load change curves of piles in slope foundation with varying pile top constraints are shown in Figure5.Figure 5
Critical buckling load with varying pile top constraints.It can be seen from Figure5 that critical buckling load is of the minimum value when pile top is free. Piles are most prone to buckling instability. Critical buckling load is of the maximum value when pile top is fixed. Piles are the least prone to buckling instability. Critical buckling load is between the previous values when pile top is elastic embedded or articulated. Critical buckling load with articulated pile top is a little more than the value with elastic embedded pile top. Therefore, when pile tip is fixed, critical buckling load of piles in slope foundation can be improved effectively with strengthening of pile top constraints.


## 6.6. Influence of Pile-Side Friction

Let ultimate friction strength between pile and soilτ be 0, 30, and 60 kPa. Pile-embedded ratio h/L is 0.6. The influence of pile-side friction upon critical buckling load is discussed. Critical buckling loads with different ultimate friction strength between pile and soil and pile length are shown in Table 4.Table 4
Critical buckling load for piles in the slope foundation with varying limit friction strength (MN).
L (m)τ = 0 kPaτ = 30 kPaτ = 60 kPa1036.08436.12936.1742068.16968.25968.34830240.879241.013241.14740632.531632.710632.888501336.451336.671336.90602447.712447.982448.25It can be seen from Table4 that when pile-side friction is considered, critical buckling load of piles in slope foundation increases. But the extent of increase is less than 0.5%. Therefore, the influence of pile-side friction upon critical buckling load can be neglected.


## 6.7. Influence of Horizontal Elastic Resistance of Pile-Side Soil

The typical methods to calculate horizontal elastic resistance of pile-side soil include Constant method,m method, and c method. Among them, m method is the method adopted in the current national code [23] to calculate foundation counterforce. At present, some scholars combine m method and c method with Constant method to calculate pile-side soil resistance. In this section, critical buckling load of piles in slope foundation is obtained with different pile-side soil resistance computation methods. The expressions assumed in the above methods are listed in Table 5 (the expressions assumed in the paper are given in equation (5)).Table 5
Pile-side soil resistance expressions.
Method nameq(x, y)0 ≤x < l1l1≤x ≤ hm methodmb0h−xyc methodmb0h−x0.5ym-constant methodmb0h−l1ymb0h−xyc-constant methodmb0h−l10.5ymb0h−x0.5yq(x, y) can be expressed as(29)qx,y=kxy,where K (x) is pile-side soil resistance coefficient. Therefore, pile-side soil resistance distribution along pile length is different for the above methods. Pile-side soil resistance is linear distribution along pile length for m method. Pile-side soil resistance is 0.5 power exponential distribution for c method. For m-Constant method and c-Constant method, pile-side soil resistance is constant when 0 ≤ x < l1 and linear or 0.5 power exponential distribution when l1 ≤ x ≤ h. As shown in expression (5), for pile-side soil resistance suggested in the paper, it is linear distribution when 0 ≤ x < l1 and square exponential distribution when l1 ≤ x ≤ h.Critical buckling load change curves of piles with varying pile-side soil resistance computation methods are shown in Figure6.Figure 6
Critical buckling load with varying pile-side soil resistance computation methods.It can be seen from Figure6 that, for medium-long piles in the range of 10 m–20 m, difference of critical buckling load with varying pile-side soil resistance computation methods is not wide. For long piles in the range of 30–40 m and superlong piles over 50 m, difference of critical buckling load with varying pile-side soil resistance computation methods appears gradually. Critical buckling load obtained with pile-side soil resistance computation method suggested in the paper is close to the value with m method adopted in the current national code. Therefore, pile-side soil resistance computation method suggested in the paper is of certain usability.


## 6.8. Influence of Pile Socketed Ratio

For piles in slope foundation with different lengths, the influence of pile socketed ratiol1/L upon critical buckling load is discussed through increase of pile socketed ratio. Critical buckling load change curves of piles in slope foundation with different lengths and pile socketed ratio are shown in Figure 7.Figure 7
Critical buckling load with varying pile socketed ratio.It can be seen from Figure7 that critical buckling load of piles in slope foundation increases with the increase of pile socketed ratio. That is to say, appropriate increase of pile socketed length is helpful to increase pile critical buckling load. Furthermore, the more of pile length, the more of increase latitude of pile critical buckling load. There is especially obvious increase latitude for superlong piles over 50 m. Critical buckling load change curves will become flat when pile socketed ratio reaches a certain value. That is to say, there is a critical fixed depth [12]. Critical buckling load change curves of medium-long piles in the range of 10 m–20 m tend to be flat earlier.


## 7. Conclusions

Based on cusp catastrophe theory, buckling stability parametric study for piles in slope foundation was made. The expression of pile critical buckling load is derived. The influences of slope angle, landslide thrust distribution behind pile, pile-embedded ratio, pile constraints, pile-side soil friction, pile-side soil resistance, and pile socketed ration upon buckling stability for piles in slope foundation are discussed. The key observations of the present work are as follows:(1)
When pile diameter remains unchanged, critical buckling load of piles in slope foundation increases with the increase of pile length if theh/L ratio is constant. When pile length remains unchanged, with the increase of pile diameter, the pile slenderness ratio decreases, and critical buckling load of piles in slope foundation increases. We get the result making pile-embedded ratio 60%.(2)
Compared with pile in level ground, critical buckling load of piles in slope foundation decreases. For long piles in the range of 30–40 m and superlong piles over 50 m, the extent of reduction of critical buckling load is obvious when the slope angle increases.(3)
Critical buckling load of piles in slope foundation is influenced by pile constraints. Suppose that pile tip constraint is fixed. Critical buckling load is of the minimum value when the pile top is free. Critical buckling load is of the maximum value when the pile top is fixed. Piles are the least prone to buckling instability. Therefore, critical buckling load of piles in slope foundation can be improved effectively with strengthening of pile constraints.(4)
Critical buckling load of piles in slope foundation is greatly influenced by pile-embedded ratio and pile socketed ratio. Pile-embedded ratio should reach 60% for piles in slope foundation. Critical buckling load of piles in slope foundation increases with the increase of pile socketed ratio. When socketed ratio reached certain value, its increasing effect to critical buckling load was exerted and critical buckling load change curves will become flat. Medium-long piles in the range of 10 m–20 m tend to reach critical socketed depth earlier.(5)
Critical buckling load of piles over 30 m is greatly influenced by pile-side soil resistance. Critical buckling load obtained with the pile-side soil resistance computation method suggested in the paper is close to the value withm method adopted in the current national code. The influence of pile-side friction and landslide thrust behind pile upon critical buckling load of piles in slope foundation is tiny, and both of them can be neglected.



---

*Source: 5560790-2021-06-21.xml*