The presence of a weak layer in a slope requires special attention because it has a negative impact on slope stability. However, limited insight into the seismic stability of slopes with a weak layer exists. In this study, the seismic stability of a pile-reinforced slope with a weak thin layer is investigated. Based on the limit analysis theory, a translational failure mechanism for an earth slope is developed. The rotational rigid blocks in the previous rotational-translational failure mechanism are replaced by continuous deformation regions, which consist of a sequence of n rigid triangles. The predicted static factor of safety and collapse mechanism in two typical examples of slopes with a weak layer compare well with the results obtained from the available literature and by using the Discontinuity Layout Optimization (DLO) technique. The lateral forces provided by the stabilizing piles are evaluated using the theory of plastic deformation. An analytical solution for estimating the critical yield acceleration coefficient for the pile-reinforced slopes is derived. Based on the proposed translational failure mechanism and the corresponding critical yield acceleration coefficient, Newmark’s analytical procedure is employed to evaluate the cumulative displacement. Considering different real earthquake acceleration records as input motion, the effect of stabilizing piles and varying the spacing of piles on the cumulative displacement of slopes with a weak layer is investigated.
National Natural Science Foundation of China513783451. Introduction
Many catastrophic slope failures have been reported in the past due to earthquakes. Estimation of the stability of slopes subjected to seismic loads is a very important task in geotechnical engineering. There are two practical estimation methods that exist for this problem: the first one is calculating the factor of safety of slopes by considering pseudostatic earthquake body forces within a soil mass (e.g., Seed et al. [1]; Seed [2]; Chen [3]). The concept of factor of safety has been widely used because it is simple and straightforward extension of static considerations. However, it provides no details regarding the shaking process. On the other hand, the technique of slopes stabilized by reinforcement is widely employed by geotechnical engineers to enhance the stability of slopes. The pseudostatic approach may underestimate the stability of reinforced slopes for large earthquake acceleration in engineering design (Ling et al. [4]; Michalowski [5]).
The second one calculates the cumulative displacement subjected to seismic loads. The most common approach is Newmark’s [6] sliding block method, which calculates the cumulative displacement of slopes by integrating earthquake acceleration in a one-block translational mechanism. This approach has the advantage of providing information during an earthquake and being less time consuming; it has been further extended to the rotational mechanism (e.g., Chang et al. [7]; Li et al. [8]) and multiblock mechanism (e.g., Michalowski [9]) of slope through limit analysis. Moreover, this approach has been employed for the seismic displacements of reinforced slopes by other researchers (e.g., Ling et al. [4]; Ling and Leshchinsky [10]; Michalowski and You [11]; He et al. [12]).
In engineering practice, the existence of a weak layer in slopes requires special attention because the low shear strength of a weak layer has an adverse effect on the performance of a slope. The previous studies on slopes with a weak layer were performed under the static condition. Fredlund and Krahn [13] compared several methods for analyzing the stability of a nonhomogeneous soil slope with a weak layer. Moreover, finite element (FE) analysis, incorporated with the shear strength reduction technique, was applied to investigate the stability of a slope with a weak layer (Griffiths and Marquez [14]; Ho [15]). Based on the upper-bound method, Huang et al. [16] proposed a rotational-translational collapse mechanism to assess the factor of safety of slopes with a weak layer. Their analytical results were verified using the finite element method. However, few studies have been conducted regarding evaluating the seismic performance of slopes with a weak layer.
In this study, a translational failure mechanism is developed to evaluate the stability of slopes with a weak layer, which is validated using other solutions with respect to the static factor of safety. Furthermore, a pseudostatic method is employed within the limit analysis framework to calculate the critical yield acceleration coefficient of a reinforced slope. Finally, considering different real earthquake acceleration records as input motion, Newmark’s analytical approach is used to assess the cumulative displacement of two typical cases of slopes with a weak thin layer.
2. Critical Yield Acceleration Coefficient for Pile-Reinforced Slopes with a Weak Layer
In limit analysis theory, soil is assumed to deform plastically according to the normality rule associated with the Mohr-Coulomb yield criterion (e.g., Lu et al. [17]). For homogeneous slopes, the rotational log-spiral failure mechanism has been found to be the most adverse for the stability of slopes (e.g., Chen [3]; Chang et al. [7]; Li et al. [8]; Gao et al. [18]; Gao et al. [19]). When a weak layer exists (the strength of a layer is relatively weak) in a slope, the slip surface of the soil slides along the weak layer (Griffiths and Marquez [14]; Huang et al. [16]). Therefore, the weak layer governs the failure mechanism and the conventional rotational log-spiral failure mechanism is not suitable. Farzaneh et al. [20] developed a rotational-translational mechanism to solve the bearing capacity problems. Recently, Huang et al. [16] proposed a rotational-translational mechanism, which contains three rigid blocks, for a slope with a weak layer, as shown in Figure 1. The velocity of block b is vb, while the angular velocities of block a and block c are ωa and ωc, respectively. Every point on discontinuous surfaces FN and O’M should satisfy the kinematically admissible velocity. Thus, computational effort is needed to determine those discontinuous surfaces, which might become an issue when deriving upper-bound limit solution in seismic analysis.
Rotational-translational mechanism (after Huang et al. 2013 [16]).
In this study, a translational failure mechanism was developed. As shown in Figure 2, the rigid rotational blocks a and c (see Figure 1) were replaced by continuous deformation regions a and c (log-spiral shear zones), consisting of a sequence of n rigid triangles. These triangles satisfy the admissible velocity field of every discontinuous surface. The number of triangles in those regions can be increased without extra analytical effort. Those continuous deformation regions a and c are separated by the translational rigid block b. Due to velocity discontinuities between adjacent blocks, energy is dissipated at their interfaces. The velocity hodographs are shown in Figure 3, in which vb is the velocity of rigid block b; vah is the velocity at the end point N of the log-spiral shear zone a; vc0 is the velocity at the start point M of the log-spiral shear zone c. The directions of the velocities incline at an angle ψ=φ with the corresponding discontinuous surfaces. These velocities constitute a kinematically admissible velocity field.
Translational collapse mechanism for pile-reinforced slopes with a weak layer.
Velocity hodographs at points (a) M and (b) N.
According to the theory of limit analysis, the slopes will collapse if the rate of work performed by external forces and soil weight exceeds the energy dissipation rate for any assumed kinematically admissible failure mechanism. The rate of external work due to the soil weight in this translational failure mechanism is the sum of different elements: (1)Wγ=Wa+Wb+Wc.The external work due to the self-weight of the log-spiral shear zones a and c can be expressed as (2)Wa=∫θa0θahma·g·va·cosθdθ=∫θa0θah12γra2-la2cosθ·va0expθ-θa0·tanφ1dθ,Wc=∫θc0θchmc·g·vc·cosθdθ=∫θc0θch12γrc2-lc2cosθ·vc0expθ-θc0·tanφ1dθ.The external work due to the self-weight of block b can be written as(3)Wb=mbgsinδ-φ2vb=γSbsinδ-φ2vb,where ra and rc are the radius of the log-spiral with respect to angles θa0 and θc0, respectively; γ is the unit weight of the soils; φ1 and φ2 are the internal friction angles of the soils, respectively. la and lc are shown in Figure 2; Sb is the area of block b.
Once the slope with a weak layer is subjected to seismic loads, the rate of the inertial force should be considered in the energy balance equation. The collapse mechanism is considered as an incipient mechanism and does not influence the shaking mode (Li et al. [8]; He et al. [12]). The effects of the seismic loads acting on the potential sliding mass are represented by horizontal forces applied to the center of soil mass. The rate of external work due to the inertial force is the sum of different elements (4)Ws=Was+Wbs+Wcs.The external work due to the inertial force of log-spiral shear zones a and c can be expressed as (5)Was=k∫θa0θahma·g·va·sinθdθ=k∫θa0θah12γra2-la2sinθ·va0expθ-θa0·tanφ1dθ,Wcs=k∫θc0θchmc·g·vc·sinθdθ=k∫θc0θch12γrc2-lc2sinθ·vc0expθ-θc0·tanφ1dθ.The external work due to the inertial force of block b can be written as(6)Wbs=kmbgcosδ-φ2vb=kγSbcosδ-φ2vb.
The technique of stabilizing piles, which has been proven as an efficient approach to enhancing the stability of a slope, is widely adopted in geotechnical engineering practice. It provides lateral force for resistance to the sliding of a slope (He et al. [21]). The total energy dissipation of pile-reinforced slopes is the sum of Dc and Dp, which can be expressed as (7)D=Dp+Dc,where Dp is the rate of work dissipation caused by the piles; Dc is the rate of energy dissipation of soil cohesion. To calculate the effect of the piles, a lateral force is assumed to act at the potential failing surface. The energy dissipation rate can be written as(8)Dp=Fpvicosδ-φ2,where Fp is the lateral force exerted by the piles per unit thickness of the sliding mass; vi is the velocity of the sliding soil.
The theoretical method proposed by Ito and Matsui [22] was adopted in this analysis to calculate lateral force Fp acting on the piles. The lateral force Fp can be expressed as follows:(9)Fp=AcH1NφtanφexpD1-D2D2Nφtanφtanπ8+φ4-2Nφ1/2tanφ-1+2tanφ+2Nφ1/2+Nφ-1/2Nφ1/2tanφ+Nφ-1-cHD12tanφ+2Nφ1/2+Nφ-1/2Nφ1/2tanφ+Nφ-1-2D2Nφ-1/2+γH22NφAexpD1-D2D2Nφtanφtanπ8+φ4-D2,where D1 is the center-to-center spacing between neighboring piles; D2 is the clear spacing between neighboring piles; H is the height of the unstable soil layer from the ground surface to the slip surface; Nφ=tan2(π/4+φ/2); A=D1(D1/D2)(Nφ1/2tanφ+Nφ-1). It can be observed that the lateral force increases as the ratio D2/D1 decreases. The method assumes that the piles are rigid and the soil surrounding the piles deforms plastically. Li et al. [8] reported that when the spacing D2/D1 is below approximately 0.5–0.6, piles act like a continuous wall and soil around the piles cannot reach the plastic state. Thus, values of D2/D1 larger than 0.5 were used in this study.
The rate of internal energy dissipation Dc induced by soil cohesion in the collapse mechanism can be calculated as follows:(10)Dc=Da+DFN+Db+DEM+Dc,in which(11)Db=lMNc2vbcosφ2,DFN=c1lFNsin2φ2vb,DEM=c1lEMsinφ1-φ2vb.
Energy dissipation in the radial shear regions a and c is the sum of the rate of dissipation along the discontinuity surfaces and the rate of dissipation within the continually deforming log spiral shear regions. As shown in Figure 2, energy dissipation takes place along the radial lines OA0, OA1, OA2, and so forth because all the triangles move as rigid blocks in directions that make an angle φ1 with the discontinuity surfaces. Energy is also dissipated on the discontinuity surfaces A0, A1, A2, A3, A4, A5, and A6 between the material at rest and the material in motion. The energy dissipation along the discontinuity surfaces is the cohesion multiplied by the relative velocities and the length of discontinuity lines in the soil mass(12)Da=c1∫θa0θah2ra0-la·va0exp2θ-θa0·tanφ1dθ,Dc=c1∫θc0θch2rc0-lc·vc0exp2θ-θc0·tanφ1dθ.
To determine the critical yield acceleration coefficient kc, at which the slope will be in the limit state, kc can be calculated by equating the rate of internal energy dissipation to the external rate of work; that is, (13)Wγ+Ws=Dp+Dc.
The critical yield acceleration coefficient kc can be calculated as(14)kc=Dp+Dc-Wγ∫θa0θahma·g·va·sinθdθ+mbgcosδ-φ2vb+∫θc0θchmc·g·vc·sinθdθ,where ma, mb, and mc are the masses of different elements; g is the gravitational acceleration. The minimum critical yield acceleration coefficient was obtained with respect to the parameters of θa0, θah, θc0, θch, vah/vb, and vc0/vb.
3. Assessment of Cumulative Displacement of Pile-Reinforced Slopes with a Weak Layer Subjected to Seismic Loads
Once the horizontal earthquake acceleration exceeds the critical value, the displacement of the soil mass begins to accumulate. Newmark’s sliding block technique was adopted to evaluate the cumulative displacement of the slope subjected to seismic loads. The classical sliding block technique involves only one block; this method was employed by other researchers to evaluate the stability of a slope under the rotational failure mechanism in both 2D and 3D problems (e.g., Chang et al. [7]; Li et al. [8]; He et al. [12]). This technique was adopted in this study for the translational mechanism of a pile-reinforced slope with a weak layer. When the earthquake acceleration is exceeded, the sliding block starts to accelerate with acceleration items u¨a, u¨b, and u¨c; thus, balancing the work rate equation yields (15)Wγ+Ws=Dc+Dp+∫θa0θahmavau¨adθ+∫θc0θchmcvcu¨cdθ+mbvbu¨b.
According to Newmark’s analytical approach [6], the sliding mass moves only downhill, irrespective of the direction of the earthquake acceleration. Earthquake acceleration varies during the process of shaking. The slope movement will occur when the seismic coefficient k exceeds its critical yield acceleration coefficient kc. The velocity of the failing mass increases from zero and reaches its peak value when the earthquake acceleration decreases to a certain critical value. After that, the failing mass will still move until the acceleration drops back down to some level below the critical yield acceleration coefficient kc (He et al. [12]). The acceleration can be integrated to calculate the velocity of the failing soil mass, while cumulative displacement can be estimated by integrating velocity (Newmark [6]). Subtracting (14) from (15), the acceleration can be expressed as(16)gk-kc∫θa0θahma·vavb·sinθdθ+mbcosδ-φ2+∫θc0θchmc·vcvb·sinθdθ=u¨b∫θa0θahma·vavb·u¨au¨bdθ+∫θc0θchmc·vcvb·u¨cu¨bdθ+mb.
To satisfy the kinematical admissibility, the ratios u¨ah/u¨b and u¨b/u¨c0 are the same as vah/vb and vb/vc0 (Michalowski and You [11]), which can be expressed as(17)vahvb=u¨ahu¨b=Aa,vc0vb=u¨c0u¨b=Ac.Thus, the acceleration u¨b becomes(18)u¨b=gk-kc∫θa0θahma·Aa·sinθdθ+mbcosδ-φ2+∫θc0θchmc·Ac·sinθdθ∫θa0θahmaAa2dθ+mb+∫θc0θchmcAc2dθ=gk-kcC.
This expression indicates that soil acceleration is not related to the piles directly but the piles affect the failure mechanism and the critical yield acceleration coefficient kc. The failure mechanism tends to reach below the toe of slopes due to the location of the weak layer. The horizontal displacement in the toe of slopes must be determined, though the analysis allows for calculation of the displacement of any point in the failing mass; that is, (19)ubx=cosδ-φ2∬tu¨bdtdt=Ccosδ-φ2∬tgk-kcdtdt.
Coefficient C was calculated based on the proposed translational mechanism. Additionally, the geometry of the slope and soil properties are contained in coefficient C, which is involved in the optimization of the collapse mechanism. The double integral in (19) is dependent on the earthquake acceleration record and the critical yield acceleration coefficient of the pile-reinforced slope.
4. Comparison with Other Results
Before the proposed translational mechanism was used for the analysis of cumulative displacement, it was validated using other solutions with respect to the static factor of safety. The Monte Carlo technique (Hammersley and Handscomb [23]) was employed for determining the critical collapse mechanism and the critical yield acceleration coefficient kc. This technique is efficient in the search for the critical slip surface of slope stability analysis (e.g., He et al. [12]). To obtain a more reliable solution, efficient searches are needed (500,000 trials).
For verification of the proposed solution, the factor of safety predicted by the proposed limit analysis was compared with the results obtained by using different computational approaches. The solutions from the available literature and the numerical program LimitState: GEO (LimitState [24]) were selected for comparison. LimitState: GEO, which has been successfully adopted in evaluating slope stability and other geotechnical problems (e.g., Smith and Gilbert [25, 26]; Leshchinsky and Ambauen [27]), employs the Discontinuity Layout Optimization (DLO) technique to evaluate the critical failure mechanisms and associated upper-bound solution. As shown in Figure 4, it automatically searches for the critical failure surface in a soil mass that is at a limit state and has no need to preassume a critical collapse mechanism. Thus, it can be considered as a useful tool for ensuring the accuracy of the failure mechanism assumed in this study. Two typical cases of slopes with a weak thin layer were selected. The comparison is based on natural slopes without reinforcement.
Procedure of DLO to determine critical layout of slip-lines: (a) starting problem; (b) discretization of soil by nodes; (c) interconnection of nodes with potential discontinuities; (d) identification of critical subset of potential discontinuities by optimization (adapted from Smith and Gilbert [25]).
Slope A, shown in Figure 5, is taken from Huang et al. [16]. A nonhomogeneous slope for β=45° has a height of 6 m and is built from a soil characterized by φ1=10°, c1=20 kPa, and γ1=20 kN/m3, with a 0.5 m thick weak layer characterized by φ2=5°, c2=12 kPa, and γ2=20 kN/m3. To validate the derived equation and the random search approach, the factor of safety computed by the translational mechanism was compared with the results obtained by Huang et al. [16] and the DLO. The comparison is shown in Table 1. It can be observed that the calculated values from the proposed translational mechanism and the existing rotational-translational mechanism (Huang et al. [16]) were slightly higher than the result from DLO (≈4% and 1%, resp.). The computed critical failure mechanism, that is, panhandling failure mechanism that partially penetrates the weak layer, can be observed in Figure 7. While the values of the factor of safety are quite close, the slip surface from Huang et al. [16] is not similar to that of DLO. In their study, the velocities of point N in block a and point M in the block c were fixed. It means that the ratios vah/vb and vc0/vb were not considered in the optimization, which contributes to the deviation of the slip surface. On the other hand, the predicted critical failure surface from the translational mechanism was similar to that of DLO, which proves the reliability of the random search approach.
Comparison of the results of factor of safety (slope A).
Factor of safety
Present study
LA (Huang et al. [16])
FEM (Huang et al. [16])
DLO
1.31
1.27
1.28
1.25
Schematic diagram of slope A.
Slope B is a nonhomogeneous slope underlain by a thin weak layer, as shown in Figure 6 (Fredlund and Krahn [13]). A 1 : 2 slope has a height of 10 m, and the soil properties are φ1=20°, c1=29 kPa, and γ1=18.8 kN/m3, with a weak layer (φ2=10°, c2=0 kPa, and γ2=18.8 kN/m3) in the foundation of the slope. Fredlund and Krahn [13] compared several methods for calculating this slope. The limit equilibrium solutions (Spencer’s method and Janbu’s rigorous method) published by Fredlund and Krahn [13] were selected in this section, whereas the results of DLO were also compared. As shown in Table 2, the present upper-bound solution yields a factor of safety of 1.36, which is in the vicinity of that of other methods. Figure 8 shows the comparison of the predicted failure mechanism between the upper-bound solution and limit equilibrium method for slope B. Furthermore, the critical slip surface from the limit analysis is reasonably consistent with that of DLO but quite different from the one of Fredlund and Krahn [13]. In the limit equilibrium method (Fredlund and Krahn [13]), the analytical parameters (e.g., failure radii, number of slices, and interslice forces) require high computational effort. The discrepancy of the slip surfaces may be attributed to the limitation of computational constraints in the limit equilibrium. Thus, the factor of safety is sensitive to the slip surface.
Comparison of the results of factor of safety (slope B).
Factor of safety
Present study
Spencer’s method (Fredlund and Krahn [13])
Janbu’s rigorous method (Fredlund and Krahn [13])
DLO
1.36
1.377
1.432
1.32
Schematic diagram of slope B.
Comparison of failure mechanisms of slope A.
Comparison of failure mechanisms of slope B.
The agreement is satisfactory, in terms of both factor of safety and the location of its associated critical failure surface. The computed results show good agreement with the values obtained by other approaches. With the confidence gained through the validation, the next section examines the cumulative displacement of pile-reinforced slopes with a weak thin layer.
5. Example
The cumulative displacements of slopes with a weak layer subjected to seismic loads were evaluated by using Newmark’s analytical approach. Moreover, the effects of the stabilizing piles on the cumulative displacement were also compared. The properties and geometries of the slope are shown in Figures 5 and 6. Ho [15] employed the finite element method (FEM) to draw the conclusion that the middle portion of the slope is the optimal pile position for slopes with a weak layer. Thus, the pile was assumed to be in the middle of the slope to assess the seismic displacement in this study.
The displacement of a structure due to seismic loads depends on the different real earthquake acceleration time histories. To make the results useful, the cases selected in this study were evaluated for different earthquake acceleration records. As shown in Figures 9–11, three typical earthquake records (He et al. [12]) were used to estimate the seismic displacements. Figure 9 is the Imperial Valley-06 Earthquake records, with peak earthquake acceleration (PGA) of 0.307 g and a record duration of 14.76 s. Figure 10 is the Kobe Earthquake record, with PGA = 0.345 g and a record duration of 20 s. Figure 11 is the Parkfield-02 Earthquake record, with PGA = 0.373 g and a duration of 21.07 s. The constant time interval of the Imperial Valley-06 and Kobe Earthquake records is 0.01 s, while the constant time interval of the Parkfield-02 Earthquake record is 0.005 s. These earthquake records were used one after another to assess the performance of the slopes (with/without reinforcement) subjected to seismic loads. The results are listed in Table 3. Calculations of the slope displacements are illustrated for the above two example slopes.
Results of cumulative displacement for slopes with a weak layer with and without piles.
Without piles
With piles
Slope A
Slope B
Slope A
Slope B
Slope A
Slope B
Slope A
Slope B
D2/D1
—
—
0.9
0.9
0.85
0.85
0.8
0.8
kc
0.145
0.122
0.257
0.172
0.273
0.190
0.294
0.211
ux (cm)
Imperial Valley-06
16.89
29.49
0.324
8.54
0.068
5.422
0.006
2.907
Kobe
42.43
67.59
1.673
23.20
1.002
15.17
0.415
8.428
Parkfield-02
49.13
68.45
8.845
30.99
6.490
23.06
3.955
16.75
Acceleration records of Imperial Valley-06 Earthquake (PGA = 0.307 g, duration = 14.76 s).
Acceleration records of Kobe Earthquake (PGA = 0.345 g, duration = 20 s).
Acceleration records of Parkfield-02 (CA) Earthquake (PGA = 0.373 g, duration = 21.07 s).
The seismic loads of the Imperial Valley-06 Earthquake, Kobe Earthquake, and Parkfield-02 Earthquake were exerted on the slopes, resulting in displacements of 16.89 cm, 42.43 cm, and 49.13 cm for slope A and 29.49 cm, 67.59 cm, and 68.45 cm for slope B, respectively. The maximum horizontal displacements of the slopes increase with the peak earthquake acceleration of the earthquake experiences. Michalowski and You [11] introduced a displacement criterion for a slope subjected to seismic loads. In the criterion, a slope is allowed to move with a small displacement. The maximum horizontal displacement at the toe ux should be smaller than 0.005H. According to the results, the displacements of those two cases exceed this criterion. The piles must be embedded to improve the performance of the slope, and the different ratios of spacing between neighboring piles D2/D1 of 0.80, 0.85, and 0.90 were compared in this study. The ratio D2/D1 of 0.80 should be used for the Parkfield-02 Earthquake with PGA = 0.373 g, while the ratio D2/D1 of 0.90 is sufficient for slope A under the Imperial Valley-06 Earthquake. Generally, the results show the efficiency of piles in limiting the deformation of the slope under earthquake loading, especially for slope A. This may be because the location of the weak layer in slope A is deeper. Thus, the height of the portion of piles above the sliding surface was relatively higher, resulting in larger lateral force Fp of the piles. For Parkfield-02 Earthquake record (PGA = 0.373 g), the reduced displacement ratios for slope B (H=10 m) are 45%, 66%, and 75.5% with different D2/D1 ratios, respectively. Li et al. [8] evaluated the effect of piles (D2/D1=0.6) on the displacement of homogeneous slopes with a height of 13.7 m. Under a seismic load with PGA = 0.4 g, the reduced ratios of displacement ranging from 68% to 86% were reported in their study. For Parkfield-02 Earthquake record, He et al. [12] reported that the piles (D2/D1=0.85) reduced the cumulative displacement ratio of a homogeneous slope (H=12 m) by about 50%. This comparison further proves that the D2/D1 ratio significantly influences the deformation of a slope. Comparing with the homogeneous slope, the piles contribute to the larger benefit in slopes with a weak layer. This is because the pile reduces the depth of slip surface in homogeneous slopes (Ausilio et al. [28]), leading to the limiting effect of pile lateral force.
6. Conclusion
Previous studies of slopes with a weak layer have focused on the static factor of safety. The seismic performance of pile-reinforced slopes with a weak layer is evaluated in this study. Based on the upper-bound theorem of limit analysis, a translational mechanism is developed to determine the critical yield acceleration coefficient and the accumulative displacement. The rigid rotational blocks in the rotational-translational failure mechanism proposed by Huang et al. [16] were replaced by continuous deformation regions a and c, with a sequence of n rigid triangles. The factor of safety and failure mechanism, evaluated by using the proposed translational failure mechanism, exhibit excellent agreement with other solutions of the available literature and the DLO technique, which proves the accuracy of the proposed failure mechanism.
To account for the presence of piles, the plastic theory is employed to calculate the lateral force and energy dissipation of the piles. Newmark’s analytical procedure is adopted to assess the cumulative displacement of pile-reinforced slopes with a weak layer subjected to seismic loads. Two typical cases existing in the available literature are shown to illustrate the influence of an earthquake on the displacements of slopes. The earthquake experiences and the peak earthquake acceleration have a significant influence on the displacement of the slope. Stabilizing piles can be regarded as an efficient approach in reducing the cumulative displacement of slopes with a weak thin layer, and reducing the spacing of piles is needed for large earthquake acceleration.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This research was funded by the National Natural Science Foundation of China (Grant no. 51378345). The authors would also like to thank the China Scholarship Council for their support.
SeedH. B.LeeK. L.IdrissI. M.Analysis of the sheffield dam failure196995614531490SeedH. B.Considerations in the earthquake-resistant design of earth and rockfill dams197929321526310.1680/geot.1979.29.3.2152-s2.0-0018655834ChenW. F.Plasticity in soil mechanics and landslides19801063443464LingH. I.LeshchinskyD.PerryE. B.Seismic design and performance of geosynthetic-reinforced soil structures199747593395210.1680/geot.1997.47.5.9332-s2.0-0031468596MichalowskiR. L.Soil reinforcement for seismic design of geotechnical structures1998231-211710.1016/S0266-352X(98)00016-02-s2.0-0032464864NewmarkN. M.Effects of earthquakes on dams and embankments196515213916010.1680/geot.1965.15.2.139ChangC.-J.ChenW. F.YaoJ. T. P.Seismic displacements in slopes by limit analysis1984110786087410.1061/(asce)0733-9410(1984)110:7(860)LiX.HeS.WuY.Seismic displacement of slopes reinforced with piles2010136688088410.1061/(asce)gt.1943-5606.00002962-s2.0-77953345157MichalowskiR. L.Displacements of multiblock geotechnical structures subjected to seismic excitation2007133111432143910.1061/(ASCE)1090-0241(2007)133:11(1432)2-s2.0-36349011798LingH. I.LeshchinskyD.Effects of vertical acceleration on seismic design of geosynthetic-reinforced soil structures199848334737310.1680/geot.1998.48.3.3472-s2.0-0032099877MichalowskiR. L.YouL.Displacements of reinforced slopes subjected to seismic loads2000126868569410.1061/(ASCE)1090-0241(2000)126:8(685)2-s2.0-0034255302HeY.HazarikaH.YasufukuN.HanZ.LiY.Three-dimensional limit analysis of seismic displacement of slope reinforced with piles20157744645210.1016/j.soildyn.2015.06.0152-s2.0-84938076841FredlundD. G.KrahnJ.Comparison of slope stability methods of analysis197714342943910.1139/t77-0452-s2.0-0017526261GriffithsD. V.MarquezR. M.Three-dimensional slope stability analysis by elasto-plastic finite elements200757653754610.1680/geot.2007.57.6.5372-s2.0-34548219020HoI.-H.Numerical study of slope-stabilizing piles in undrained clayey slopes with a weak thin layer201415510.1061/(asce)gm.1943-5622.00004452-s2.0-84941702007HuangM.WangH.ShengD.LiuY.Rotational-translational mechanism for the upper bound stability analysis of slopes with weak interlayer20135313314110.1016/j.compgeo.2013.05.0072-s2.0-84879767538LuX.WangH.HuangM.Upper bound solution for the face stability of shield tunnel below the water table201420141172796410.1155/2014/7279642-s2.0-84908440036GaoY.ZhangF.LeiG. H.LiD.WuY.ZhangN.Stability charts for 3D failures of homogeneous slopes201313991528153810.1061/(ASCE)GT.1943-5606.00008662-s2.0-84883341540GaoY.YangS.ZhangF.LeshchinskyB.Three-dimensional reinforced slopes: evaluation of required reinforcement strength and embedment length using limit analysis201644213314210.1016/j.geotexmem.2015.07.0072-s2.0-84939426245FarzanehO.GanjianN.AskariF.Rotation-translation mechanisms for upper-bound solution of bearing capacity problems201037444945510.1016/j.compgeo.2010.01.0042-s2.0-77952546851HeY.HazarikaH.YasufukuN.TengJ.JiangZ.HanZ.Estimation of lateral force acting on piles to stabilize landslides20157931981200310.1007/s11069-015-1942-02-s2.0-84947019440ItoT.MatsuiT.Methods to estimate lateral force acting on stabilizing piles1975154435910.3208/sandf1972.15.4_432-s2.0-0016644650HammersleyJ. M.HandscombD. C.1964London, UKMethuenLimitState2013Sheffield, UKLimitStateSmithC. C.GilbertM.Identification of rotational failure mechanisms in cohesive media using discontinuity layout optimisation201363141194120810.1680/geot.12.P.0822-s2.0-84893005010SmithC. C.GilbertM.Application of discontinuity layout optimization to plane plasticity problems200746320862461248410.1098/rspa.2006.1788MR23528682-s2.0-36348967811LeshchinskyB.AmbauenS.Limit equilibrium and limit analysis: comparison of benchmark slope stability problems2015141100401504310.1061/(ASCE)GT.1943-5606.0001347AusilioE.ConteE.DenteG.Stability analysis of slopes reinforced with piles200128859161110.1016/S0266-352X(01)00013-12-s2.0-0034875979