Next Article in Journal
Influence of Subsoil and Building Material Properties on Mine-Induced Soil–Structure Interaction Effect
Previous Article in Journal
Features of the Stress–Strain State of 3D Metal Objects Produced by Additive Microplasma Deposition of the Powder of a Fe–Cr–Ni–B–Si System
Previous Article in Special Issue
Management System According to ISO/IEC 17025: Method Validation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Solution for Seismic Stability of 3D Rock Slope Reinforced with Prestressed Anchor Cables

1
School of Civil Engineering, Central South University, Changsha 410075, China
2
Department of Civil & Architecture Engineering, Changzhou Institute of Technology, Changzhou 213032, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(10), 4160; https://doi.org/10.3390/app14104160
Submission received: 11 March 2024 / Revised: 21 April 2024 / Accepted: 11 May 2024 / Published: 14 May 2024
(This article belongs to the Special Issue Slope Stability and Earth Retaining Structures)

Abstract

:
Currently, the study of analytical solutions for the seismic stability of slopes under anchorage conditions is one of the hottest subjects in engineering. In this paper, an analytical solution for the seismic safety of the three-dimensional (3D) two-stage rock slope reinforced with prestressed anchor cables governed by the nonlinear Hoek–Brown criterion was deduced, in which the analyses of seismicity were performed by the latest modified pseudo-dynamic method. This method supplements the consideration of the damping effect of the rock medium on seismic waves, which is more in line with the real seismic situation. A mathematical geometric model was developed for calculating the external forces work and internal energy dissipation acting on 3D rotating rock masses reinforced with prestressed anchor cables, in which the seismic work rate was calculated using a new layer-by-layer superposition summation method. The analytical solution of the safety factor could be collated as an explicit function of several variables, and then, the optimal value was obtained by the Particle Swarm Optimization (PSO) algorithm. To corroborate the accuracy of new analytical solutions, the results were contrasted with those of the pertinent literature. The results of the two comparisons were very close. Ultimately, the sensitivity analyses and coupling effects of seismic pseudo-dynamic factors and prestressing anchorage factors were carried out. It was found that even small seismic intensities had a large effect on the stability of rock slopes with developed joints. Increasing the number of steps and prestressing anchors can effectively improve the stability of rock slopes under seismic effects. The conclusions have significant implications for the anchorage design of the 3D two-stage rock slope as seismic events occur.

1. Introduction

The analytical computation for rock slope stability disturbed by large seismicity has traditionally been one of the challenges in the field of engineering. For slopes with poor-quality lithological media, prestressed anchor cables are commonly employed as a means of reinforcement [1,2,3,4]. It can provide the resistance force acting on the structure to bear the external load and creates a compressive stress zone within the anchored rock layer, which allows the weak rock slope to resist heavy seismic effects [5,6]. Moreover, there is a considerable amount of research on two-dimensional (2D) slopes reinforced with prestressed anchor cables in the previous literature, while slopes are indeed frequently in the three-dimensional (3D) status, which has received little attention [7,8]. Currently, the safety assessment approaches of the 3D slope are primarily of the following three types: (1) limit equilibrium (LE), (2) limit analysis (LA), and (3) finite element (FE). Among these, the LA of plasticity ignores the elastic–plastic process during rock mass failure and directly analyzes the limit state of 3D slopes [9]. Since LA can readily obtain the exact analytical solution [10,11] as well as avoid the irrational assumptions of LE and the time-consuming computations of FE [12], this method is utilized in this paper.
Building the kinematically admissible failure mechanism is the main difficulty with utilizing LA to solve the stability solution, and there is a wealth of literature in this area of study [13,14,15]. For the first time, a 3D sliding mechanism composed of numerous logarithmic spirals was created in order to investigate the 3D stability of vertically tangential slopes at different angles [16]. A 3D integral sliding mechanism was constructed, and the effect of slope top foundation load on the overall and local stability of the slope was examined using this sliding mechanism [17]. Considering the volume deformation of sliding bodies for slopes, a horn-shaped rotating sliding mechanism for 3D stability analysis of slopes with soil or rock was proposed [18], and it was adopted in this paper.
For rock media, the Hoek–Brown (H-B) yield criterion is one of the most frequently employed strength criteria [19]; it was proposed based on triaxial tests on hundreds of rock specimens and many statistical analyses for engineering field tests, combined with a large amount of previous rock failure theories in the 1980s [20]. In 1992, Hoek et al. [21] published, the generalized Hoek–Brown yield criterion based on the existing criterion by introducing the Geological Strength Index (GSI), which is applicable to the case of weak rock masses. In 2002, the perturbation factor D was introduced to consider the effects of blasting and stress relief on the stability of the rock masses [22]. Nevertheless, since the H-B yield criterion reflects the relationship between the maximum and minimum principal stresses when destruction occurs in rock structures, this is not easy to utilize directly in the structural stability analysis of rock engineering. In order to address this issue, the generalized tangent technique and the equivalent Mohr–Coulomb strength criterion (EMC) have been proposed by scholars, respectively, for incorporating the H-B yield criterion into rocky stability analysis [23,24], in which the former is employed in the article.
Recently, the impact of seismicity on the safety of slopes has been the subject of numerous study findings [25,26,27]. On the premise of the pseudo-static method, the seismic safety of 3D steep slopes was explored. The seismic safety of 3D steep slopes was investigated based on the pseudo-static approach’s premise [28]. However, the pseudo-static method does not take into consideration the impact of seismic arousal characteristics on the stability of slopes and evaluates seismic effects in a very simple manner [29]. Steedman and Zeng suggested a simplified pseudo-dynamic method that switches the account of the time effect for seismicity while it ignores the significant variables involving boundary state and medium damping [30,31]. The modified pseudo-dynamic method put forward by Bellezza, which was developed from earlier approaches, has been extensively employed in the seismic safety assessment of geotechnical engineering [32,33,34]. Thus, this approach is also utilized in this study.
According to related research, altering the geometry of slopes can also significantly enhance the aseismatic performance of slopes while being more cost-effective than providing additional reinforcement. Thus, the aseismatic effectiveness of coupling the prestressed anchor cable [35] with the two-stage slope, that is, adding a step platform on the slope surface, is investigated in the present article. For the studies of slope anchorage based on the LA method, Li et al. investigated the stability of 2D slopes reinforced with a single-row anchor using the rotating failure mechanism and discovered that there is a risk of local instability of the slope when the anchor is close to the toe of the slope [36]. And Yan et al. proposed a method for assessing the stability of 2D slopes under seismic conditions based on the modified pseudo-dynamic method, which considers the dynamic changes in the axial force of the anchor cable [37]. Based on the foregoing literature, this paper expands the 2D slope condition to 3D.
In this paper, an analytical solution for the seismic safety of the 3D two-stage rock slope reinforced with prestressed anchor cables governed by the nonlinear H-B criterion is presented, in which the analyses of seismicity are performed by the latest modified pseudo-dynamic method. And a 3D horn-shaped rotating failure model is utilized to investigate and determine the external forces and internal energy dissipation. With the prerequisite, FS can be compiled as an explicit function with regard to the five variables. To confirm the correctness for analytical solutions, the results are contrasted with those of the pertinent literature and other approaches, and the errors were all within the acceptable range. Eventually, the sensitivity analyses and coupling effects of the seismic pseudo-dynamic parameters and prestressing anchorage parameters were carried out.

2. Modifications for Hoek–Brown Yield Criterion and the Pseudo-Dynamic Method on Seismicity Analysis

2.1. Nonlinear Hoek–Brown Strength Criterion

At the end of the last century, the nonlinear H-B strength criterion was devised by Hoek and Brown [20], and the characteristic that the failure envelope shows significant nonlinearity during the destruction of rocks was further elucidated. It has been widely used in practical geotechnical engineering, due to its material parameters that can clearly characterize the state or quality of the rock. The nonlinear H-B failure criterion is expressed by the following equation of the major and minor principal (effective) stresses σ 1 and σ 3
σ 1 = σ 3 + σ c m σ 3 σ c + s n
where σ c is the uniaxial compressive strength of intact rock. And the material parameters m, s, and n are defined as
m = m i e G S I 100 28 14 D
s = e G S I 100 9 3 D
and
n = 1 2 + 1 6 e G S I 15 e 20 3
where GSI, mi, and D are the Geological Strength Index, which describes rock quality, the parameter associated with the rock material type, and the rock perturbation coefficient characterizing blast effects and stress relief, respectively. More detailed definitions of these parameters are presented by Hoek et al. [22]. The homogeneity and isotropy of the rock mass are the primary condition for the H-B failure criterion. More precisely, it should be presupposed that the rocky slope is either intact or composed of strongly fractured rock when the H-B strength criterion is employed for stability analysis of the rocky slope.

2.2. Generalized Tangent Technique

Another challenge is that the principal stresses σ 1 σ 3 mode of the H-B criterion presented in Equation (1) cannot be directly enforced in the computations of the rate of energy dissipation along the failure surfaces delineated by the assumed kinematically admissible mechanism. A solution to this problem was to incorporate the H-B criterion into the stability analysis of rock engineering structures using the generalized tangent technique [13,23]. The nonlinear characteristic of the H-B yield criterion is linearized in the technique, i.e., the outer tangent of the H-B criterion envelope is used instead of the envelope curve itself, which overestimates the rock strength but leads to the solution of the upper limit for the ultimate load. The envelope curve of the nonlinear H-B yield criterion and its outer tangent are illustrated in Figure 1, and the outer tangent can be expressed as the following equation.
τ = c t + σ n tan φ t
in which c t is the cut-off distance of the tangential line on the vertical axis, and φ t denotes the inclination of the tangent line about the horizontal axis. It can be observed that the above equation closely resembles the form of the classical Mohr–Coulomb (M-C) criterion expression. Thus, the parameters c t and φ t can be visualized as the cohesive force and the internal friction angle for a particular material, respectively, while the cohesive force c t can be written as a function of the internal friction angle φ t
c t σ c = cos φ t 2 m n ( 1 sin φ t ) 2 sin φ t n 1 n tan φ t m 1 + sin φ t n m n ( 1 sin φ t ) 2 sin φ t 1 1 n + s m tan φ t
As can be seen in Figure 1, the failure strength corresponding to the tangential line is more than or equivalent to the actual H-B criterion strength when the positive stress σ n is determined [23]. Consequently, the solution obtained by the generalized tangent technique is the upper bound of the solution by the nonlinear H-B criterion.

2.3. Modified Pseudo-Dynamic Method

Seismicity could dramatically enhance induced force for slope instability and decrease the shearing intensity for rock mass, which is an essential factor leading to the failure of rocky slopes [38,39]. Within the classical pseudo-static model, the seismic effect is normally simplified to a uniform inertial force, and this approximate treatment ignores the dynamic characteristics of seismicity [40,41]. In the conventional pseudo-dynamic approach, the seismic behavior is typically denoted in the form of a seismic wave motion [42]. Nevertheless, it postulates that the slope medium is a linear elastomer, which causes the seismic acceleration to be infinitely amplified, deviating from the actual reality [43]. To overcome these drawbacks, a pseudo-dynamic approach based on viscoelastic standing waves is proposed [32,33], in which the geotechnical medium is presumed to be the Kelvin–Voigt (K-V) constitutive model; thus, the damping effects of the geotechnical material are considered, such that the seismic acceleration is not infinitely amplified. In accordance with the boundary conditions: (1) the free surface shear stress at the top of the slope is zero; (2) the displacement at the foot of the slope (z = 0) ( u h t = u h 0 cos ( ω t ) ) and the explicit equation of the horizontal displacement of transverse wave within the slope mechanism can be obtained based on the improved pseudo-dynamic method.
u h ( z , t ) = u h 0 C S 2 + S S 2 C S C SZ + S S S SZ cos ( ω t ) + S S C SZ + C S S SZ sin ( ω t )
in which
C s z = cos y S 1 H z H cos h y S 2 H z H
S s z = sin y S 1 H z H sin h y S 2 H z H
C S = cos y S 1 cos h y S 2
S S = sin y S 1 sin h y S 2
y s 1 = ω H V s 1 + 4 D 2 + 1 2 1 + 4 D 2
y s 2 = ω H V s 1 + 4 D 2 1 2 1 + 4 D 2
By second-order differentiation of Equation (7), the horizontal seismic acceleration can be deduced as
a h z , t = k h g C S 2 + S S 2 C S C SZ + S S S SZ cos ω t + S S C SZ + C S S SZ sin ω t
where k h g = ω 2 u h 0 ; k h and g are the level seismicity factor and gravity accelerated speed, respectively. y s 1 and y s 2 are defined as the correlation functions of ω H / V s and D r , where H is the slope’s height and V s is the speed of the transverse wave propagation.
Accordingly, the displacement and acceleration of the longitudinal wave propagation per unit rock layer at a certain moment can be given similarly. However, relevant research has shown that the vertical seismic acceleration during longitudinal wave propagation is almost negligible, with kv ≤ 0.5kh [44]. Hence, the horizontal seismic effect, which has the most impact on slope stability, is taken into consideration in this paper with enough comparisons of the effects of lateral and longitudinal seismic accelerations.

3. Model for 3D Two-Stage Rock Slopes Reinforced with Prestressed Anchor Cables Effected by Seismicity and Its Analytical Solution

The upper limit approach of plastic limit analysis is the application of a set of theorems and assumptions [9], for ideally plastic materials, which permits the mechanism to strictly yield to the minimum upper limit of energy dissipation when failure occurs, and it is consistent with the deformation of the rock mass governed by the associated flow law. As a result, the principal strain rates are required to fulfill such a relationship.
ε ˙ 1 + ε ˙ 2 + ε ˙ 3 + ε ˙ 1 ε ˙ 2 ε ˙ 3 sin φ t = 0
Secondly, there exists the following relationship between the normal and tangential components vn, vt of the velocity on the sliding failure surface St.
v n = v t tan φ t
And in order to gain the safety evaluation index of 3D rock slopes, it is imperative to establish the kinematically admissible failure mechanism.

3.1. The 3D Horn-Shaped Rotating Failure Mechanism

The limit analysis method (LA) requires that the assumed failure mechanism ought to meet the kinematically admissible conditions; specifically, the sliding surface of the mechanism must be continuous, i.e., without interruption points. In 2009, the 3D horn-shaped rotating failure mechanism was designed by Michalowski and Drescher based on the 2D log-helix failure mechanism [18], in which the derived safety factor of 3D soil slopes is quite close to the engineering reality [45,46]. Subsequently, the suitability of the above-mentioned mechanism for rock slopes was extensively discussed and demonstrated when the nonlinear strength criterion was converted to the linear, employing the generalized tangent technique. Therefore, the 3D rotating failure mechanism is adopted in the article for analyzing the safety of 3D rock slopes subjected to seismicity reinforced with prestressed anchor cables.
As illustrated in Figure 2, the 3D rock slope disturbed by the seismicity is reinforced by adding a step at the appropriate vertical height of the slope and installing prestressed anchor cables on the inclined surface of the slope, in which the width CD of the added step is set to a. Additionally, it can be seen that the 3D rotating failure mechanism is a horn-shaped cone with vertex angle 2 φ t , formed by radiating two log-spiral curves at the vertex O and taking the ray segment from the origin O intercepted by these two log-helix curves as the diameter of the cone section, where the lower spiral curve passes precisely through the toe of the two-stage slope. The mathematical expressions for the upper and lower log-spiral curves can be given separately:
r = r 0 e θ θ 0 tan φ t
and
r = r 0 e θ θ 0 tan φ t
where r and r are the polar diameters of the upper and lower log-spiral curves at any polar angle, respectively. As an illustration, r 0 and r 0 are the polar diameters ( O A and O A ) of the upper and lower log-spiral curves at polar angle θ 0 . According to Equations (15) and (16), the expression for the curve formed by the successive centers of the circular section of the horn-shaped cone at any polar angle can be deduced as
r m = r + r / 2 = r 0 f 1
The radius of the circular section for the horn-shaped cone at any polar angle can be similarly shown as
R = r r / 2 = r 0 f 2
where f 1 , f 2 can be expressed as
f 1 = 1 2 e θ θ 0 tan φ t + r 0 r 0 e θ θ 0 tan φ t
f 2 = 1 2 e θ θ 0 tan φ t r 0 r 0 e θ θ 0 tan φ t
In addition, it can be seen that the overall height of the 3D two-stage rock slope is H, and the heights of the upper and lower steps are ξ 1 H and ξ 2 H , respectively, in which ξ 1 and ξ 2 are the height coefficients. Obviously,
ξ 1 + ξ 2 = 1
From the stereoscopic perspective, as shown in Figure 3, the 3D failure mechanism for the two-stage rock slope is inserted in a plane strain block with an arbitrary adjustable length, and with the ever-increasing length b of the block, the 3D slope will progressively change to the 2D case. Apparently, the following correlations exist between the width of the embedded plane strain block b, the maximum width of the 3D portion B max , and the overall slope width B:
b = B B max
The introduction of the plane strain block allows a comparison with the data results for the 2D slope condition in the subsequent sections, providing the design and maintenance basis for diverse slope engineering sites.

3.2. Internal Energy Dissipation

When the 3D two-stage rock slope is in the critical failure state, the deformation of the rock body DV and the sliding along the failure surface St between the slope body and the rock base Dt are the primary contributors of the energy dissipation inside the slope to retain its original condition. According to the theory of plastic limit analysis, DV and Dt can be deduced, individually, as
D V = V c t cos φ t ε ˙ 1 ε ˙ 2 ε ˙ 3 d V
and
D t = S t c t v cos φ t d S t
v and V are velocity at the arbitrary micro-element and the deformation volume generated inside the 3D stepped slope failure mechanism, respectively, in which v can be calculated as
v = r m + y ω
In the above equation, y represents the value that corresponds to the coordinates of the Y-axis in the regional Cartesian system, which is established on any circular slice of the horn-shaped cone mechanism and originates at the center of the cone section, as depicted in Figure 4. The total internal energy dissipation could be stated as follows, based on the divergence theorem of Gauss:
D = D t + D V = c t cot φ t S r v n d S r
in which S r represents the residual surface section for the kinematically admissible model after reducing the 3D sliding surface. v denotes the speed vector on S r , n represents the unit normal vector extrorsely, and v n is their dot product.
According to Figure 3, the internal energy dissipation of the two-stage rock slope is separated into two parts, namely the 3D section and the 2D insert section. Likewise, the internal energy dissipation along the outer contour for the two-stage rock slope could be separated into four parts, namely DAB, DBC, DCD, and DDE. In accordance with Equation (28), the 3D section can be sequentially derived as
D AB 3 D = 2 ω c t cot φ t θ 0 θ B 0 x 1 * r 0 2 sin 2 θ 0 cos θ sin 3 θ d x d θ
D BC 3 D = 2 ω c t cot φ t θ B θ C 0 x 2 * r 0 2 sin 2 θ B + β 1 sin 2 θ 0 sin 2 θ B cos θ + β 1 sin 3 θ + β 1 d x d θ
D CD 3 D = 2 ω c t cot φ t θ C θ D 0 x 3 * r 0 sin θ 0 + ξ 1 H 2 cos θ sin 3 θ d x d θ
D DE 3 D = 2 ω c t cot φ t θ D θ h 0 x 4 * r 0 2 e 2 θ h θ 0 tan φ t sin 2 θ h + β 2 cos θ + β 2 sin 3 θ + β 2 d x d θ
The second integral’s upper bound for integration x i ( i = 1 , 2 , 3 , 4 ) in the dual integration formulas above can be expressed as x i * = R 2 d i 2 . In accordance with Figure 4, the expressions for d i ( i = 1 , 2 , 3 , 4 ) are
d 1 = r 0 sin θ 0 sin θ r m = r 0 f 3
d 2 = r 0 sin θ 0 sin θ B + β 1 sin θ B sin θ + β 1 r m = r 0 f 4
d 3 = r 0 sin θ 0 + ξ 1 H sin θ r m = r 0 f 5
d 4 = r 0 e θ h θ 0 tan φ t sin θ h + β 2 sin θ + β 2 r m = r 0 f 6
where f 3 , f 4 , f 5 , and f 6 are formulations without dimensions and are represented as
f 3 = sin θ 0 sin θ f 1
f 4 = sin θ 0 sin θ B + β 1 sin θ B sin θ + β 1 f 1
f 5 = sin θ 0 + ξ 1 η 1 sin θ f 1
f 6 = e θ h θ 0 tan φ t sin θ h + β 2 sin θ + β 2 f 1
In the light of the geometric qualities of the log-helix mechanism in Figure 4, the angles θ B , θ C , and θ D can be derived as
θ B = arccos cos θ 0 η 2 sin 2 θ 0 + cos θ 0 η 2 2
θ C = arccos cos θ 0 η 2 ξ 1 cot β 1 η 1 sin θ 0 + ξ 1 η 1 2 + cos θ 0 η 2 ξ 1 cot β 1 η 1 2
θ D = arccos e θ h θ 0 tan φ t cos θ h + ξ 2 cot β 2 η 1 e θ h θ 0 tan φ t sin θ h ξ 2 η 1 2 + e θ h θ 0 tan φ t cos θ h + ξ 2 cot β 2 η 1 2
in which η 1 = H / r 0 and η 2 = L / r 0 are two functions without dimensions, denoted as
η 1 = e θ h θ 0 tan φ t sin θ h sin θ 0
η 2 = cos θ 0 e θ h θ 0 tan φ t cos θ h ξ 1 cot β 1 + ξ 2 cot β 2 + λ a η 1
L is the length of the line segment AB at the peak of the slope. Furthermore, the internal energy dissipation of the embedded section with the plane strain characteristic can be represented separately as
D AB embed = 2 ω c t cot φ t θ 0 θ B 0 b / 2 r 0 2 sin 2 θ 0 cos θ sin 3 θ d x d θ
D BC embed = 2 ω c t cot φ t θ B θ C 0 b / 2 r 0 2 sin 2 θ B + β 1 sin 2 θ 0 sin 2 θ B cos θ + β 1 sin 3 θ + β 1 d x d θ
D CD embed   = 2 ω c t cot φ t θ C θ D 0 b / 2 r 0 sin θ 0 + η 1 H 2 cos θ sin 3 θ d x d θ
D DE embed   = 2 ω c t cot φ t θ D θ h 0 b / 2 r 0 2 e 2 θ h θ 0 tan φ t sin 2 θ h + β 2 cos θ + β 2 sin 3 θ + β 2 d x d θ
Resultantly, the whole internal energy dissipation of the 3D horn-shaped mechanism for the two-stage slope could be written as
D = D AB 3 D + D BC 3 D + D CD 3 D + D DE 3 D + D AB   embed + D BC   embed + D CD   embed + D DE   embed = ω c t cot φ t r 0 3 g 1 + g 2
where g1 and g2 are non-dimensional formulas of the internal dissipation for the three-dimensional section and embedded section, individually, and are represented as
g 1 = 2 sin 2 θ 0 θ 0 θ B cos θ sin 3 θ f 2 2 f 3 2 d θ           2 sin 2 θ B + β 1 sin 2 θ 0 sin 2 θ B θ B θ C cos θ + β 1 sin 3 θ + β 1 f 2 2 f 4 2 d θ           2 sin θ 0 + ξ 1 η 1 2 θ C θ D cos θ sin 3 θ f 2 2 f 5 2 d θ           2 e 2 θ h θ 0 tan φ t sin 2 θ h + β 2 θ D θ h cos θ + β 2 sin 3 θ + β 2 f 2 2 f 6 2 d θ
g 2 = 1 2 b H η 1 sin 2 θ 0 sin 2 θ B 1 + sin 2 θ 0 sin 2 θ B sin 2 θ B + β 1 sin 2 θ C + β 1 1                 + sin θ 0 + ξ 1 η 1 2 1 sin 2 θ D 1 sin 2 θ C                 + e 2 θ h θ 0 tan φ t 1 sin 2 θ h + β 2 sin 2 θ D + β 2

3.3. External Force Work

With respect to the external forces, the safety of the 3D two-stage rock slope is influenced by gravity, seismicity, and prestressed anchor cables, in which gravity and seismicity are negative effects, while the prestressed anchor cables contribute positive effects. In the following, the derivation processes of external work for gravity, seismicity, and prestressed anchor cables are introduced sequentially in the framework of the horn-shaped rotating mechanism of the 3D two-stage rock slope.

3.3.1. Gravitational Work

As for the 3D part of the two-stage slope, the work done by gravity on the rock can be stated by Equation (53):
W γ 3 D = 2 ω γ θ 0 θ B d 1 R 0 x * r m + y 2 cos θ d x d y d θ + θ B θ C d 2 R 0 x * r m + y 2 cos θ d x d y d θ + θ C θ D d 3 R 0 x * r m + y 2 cos θ d x d y d θ + θ D θ h d 4 R 0 x * r m + y 2 cos θ d x d y d θ
in which γ is the bulk density of the rock. And the upper bound on the first integral of each triple integral in Equation (53) can be inferred as x * = R 2 y 2 . Ultimately, with the addition of the intervening 2D plane strain component, the complete rock weight project can be achieved as
W γ = W γ 3 D + W γ i n s e r t             = ω γ r 0 4 g 3 + g 4
where g3 and g4 are non-dimensional formulas of internal dissipation for the 3D portion and the embedded portion, individually, and are shown as
g 3 =           2 θ 0 θ B 1 8 f 2 2 f 3 1 4 f 3 3 2 3 f 1 f 3 2 1 2 f 1 2 f 3 + 2 3 f 1 f 2 2 f 2 2 f 3 2 + 1 8 f 2 4 + 1 2 f 1 2 f 2 2 arccos f 3 f 2 cos θ d θ                     + 2 θ B θ C 1 8 f 2 2 f 4 1 4 f 4 3 2 3 f 1 f 4 2 1 2 f 1 2 f 4 + 2 3 f 1 f 2 2 f 2 2 f 4 2 + 1 8 f 2 4 + 1 2 f 1 2 f 2 2 arccos f 4 f 2 cos θ d θ                     + 2 θ C θ D 1 8 f 2 2 f 5 1 4 f 5 3 2 3 f 1 f 5 2 1 2 f 1 2 f 5 + 2 3 f 1 f 2 2 f 2 2 f 5 2 + 1 8 f 2 4 + 1 2 f 1 2 f 2 2 arccos f 5 f 2 cos θ d θ                     + 2 θ D θ h 1 8 f 2 2 f 6 1 4 f 6 3 2 3 f 1 f 6 2 1 2 f 1 2 f 6 + 2 3 f 1 f 2 2 f 2 2 f 6 2 + 1 8 f 2 4 + 1 2 f 1 2 f 2 2 arccos f 6 f 2 cos θ d θ
g 4 = b H η 1 f 7 f 8 f 9 f 10 f 11
in which f 7 , f 8 , f 9 , f 10 , and f 11 are dimensionless expressions that are written as
f 7 = 3 tan φ t cos θ h + sin θ h e 3 θ h θ 0 tan φ t 3 tan φ t cos θ 0 sin θ 0 3 1 + 9 tan 2 φ t
f 8 = 1 3 η 2 sin θ 0 cos θ 0 1 2 η 2
f 9 = 1 3 ξ 1 η 1 cos θ 0 η 2 + sin θ 0 cot β 1 cos θ 0 η 2 1 2 ξ 1 η 1 cot β 1
f 10 = 1 3 λ a η 1 sin θ 0 + ξ 1 η 1 cos θ 0 η 2 ξ 1 η 1 cot β 1 1 2 λ a η 1
f 11 = 1 3 ξ 2 η 1 e θ h θ 0 tan φ t sin θ h + β 2 sin β 2 e θ h θ 0 tan φ t cos θ h + 1 2 ξ 2 η 1 cot β 2

3.3.2. Seismic Force Work

Seismicity is the other critical factor of slope instability. The seismic force can be decomposed into a horizontal seismic force and a vertical seismic force, where the vertical seismic force has much less impact on slope stability than the horizontal seismic force at kv ≤ 0.5kh, according to the results of Gazetas et al. [44], which can be neglected to simplify the computational model. In the later paper, the horizontal stratification method is employed to analyze the effect of the horizontal seismic force on slope stability, since the seismic acceleration stays invariant in the arbitrary horizontal plane.
As depicted in Figure 5a, a horizontal plane, denoted as PP’, is taken at random height Z of the slope of height H, and the polar angles corresponding to the two intersection points P and P’ of this horizontal plane and the slope profile are denoted as α and α , separately. Therefore, an arbitrary height Z within the range of slope height H in the Cartesian coordinate system x-E-z, as originated by the foot of the slope, can be stated as a function of the polar angle α .
z α = r 0 e θ h θ 0 tan φ t sin θ h r p sin α
where r P is the polar diameter corresponding to point P, which could be shown as
r p = r 0 e α θ 0 tan φ t
By putting Equation (62) into Equation (14), we can obtain seismic accelerated speed a h on the level plane of any height Z within the height H of the slope. And furthermore, as shown in Figure 5b, the work exerted by a level seismic force within volume micro-element dV of height dz can be expressed as
d W s = a h g γ v sin θ d V
Of these, the geometric equation for the volume micro-element in the 3D portion of the slope, as shown in Figure 5b, can be derived as
d V 3 D = l θ sin α sin 2 θ r p d θ d z
Likewise, as illustrated in Figure 5b, the volume micro-element in the embedded plane strain block of the slope can be similarly shown as
d V embed   = b sin α sin 2 θ r p d θ d z
In Equation (65), l ( θ ) is the breadth for the 3D section of the two-stage slope at an arbitrary polar angle θ   θ 0 θ θ h , which is shown as
l θ = 2 R 2 d 2
d = r p sin α sin θ r m = r 0 f 12
where f 12 is a non-dimensional equation, listed as
f 12 = e α θ 0 tan φ t sin α sin θ f 1
Also, the micro-element dz in Equations (65) and (66) can be stated as
d z = cos α φ t cos φ t r p d α
Thereby, with the combination of Equations (62)–(70), we can obtain the seismic force work within the 3D portion for the slope and embedded plane strain block, individually.
W s 3 D = ω γ θ 0 θ h α α l ( θ ) r p 3 sin 2 α sin 2 θ cos α φ t cos φ t a h g d α d θ
W s e m b e d = ω γ θ 0 θ h α α b r p 3 sin 2 α sin 2 θ cos α φ t cos φ t a h g d α d θ
The integration upper limit α of the dual integrals in Equations (71) and (72) can be derived from α , whose segmental function expression is represented as
α = arccos r 0 sin θ 0 tan θ B r p sin α r 0 sin θ 0 tan β 1 r 0 sin θ 0 tan θ B r p sin α r 0 sin θ 0 tan β 1 2 + r p 2 sin 2 α θ 0 α θ F     arccos r 0 sin θ 0 + ξ 1 H tan θ D r p sin α r 0 sin θ 0 ξ 1 H tan β 2 r 0 sin θ 0 + ξ 1 H tan θ D r p sin α r 0 sin θ 0 ξ 1 H tan β 2 2 + r p 2 sin 2 α       θ F < α θ h
Resultantly, the overall seismic force work of the 3D horn-shaped mechanism for the two-stage slope could be further written as
W s = W s 3 D + W s e m b e d         = ω γ k h r 0 4 g 5 + g 6
in which g5 and g6 are non-dimensional formulas for the seismic force work, individually, and are listed as
g 5 = 2 θ 0 θ h α α e 3 α θ 0 tan φ t sin 2 α sin 2 θ cos α φ t cos φ t a h k h g f 2 2 f 12 2 d α d θ
g 6 = b H η 1 θ 0 θ h α α e 3 α θ 0 tan φ t sin 2 α sin 2 θ cos α φ t cos φ t a h k h g d α d θ

3.3.3. Anchoring Force Work

Prestressed anchor cables are widely employed for slope reinforcement [47,48]. When seismic events occur, the slope is in a state of critical instability. At this time, the anchor cables apply the anchoring force to the sliding block of the slope through axial tensile deformation. In the following section, the negative work done by the anchoring force on the slope, that is, keeping the slope in the steady state, is analyzed and derived in detail.
As illustrated in Figure 6, the anchor cable numbers installed on the upper step slope are noted as T u 1 , T u 2 , …, T u i , …, T u n , with a total of n rows of anchor ropes; the anchor ropes on the lower step slope are similarly recorded as T l 1 , T l 2 , …, T l j , …, T l m , with m rows of anchor cables overall. The above numbers also represent the anchoring force values of the corresponding anchor ropes. In addition, δ u i and h u i are the anchorage angle and anchor head height of the i-th row of anchor cables on the upper slope of the step, respectively; likewise, δ lj and h l j are the anchorage angle and anchor head height of the j-th row of the anchor ropes installed on the lower step slope. For the failure mechanism depicted in Figure 6, the external force work on the sliding block of the slope done by the i-th row of anchors on the upper step slope and the j-th row of anchors on the lower step slope can be expressed, respectively, as
W ˙ u i = ( b + B u i ) T u i ω r u i sin ( θ u i δ u i )
W ˙ l j = ( b + B l j ) T l j ω r l j sin ( θ l j δ l j )
where B u i , r u i , and θ u i are the breadth for the 3D section of the two-stage slope at the altitude of the i-th row of anchor heads on the upper step slope as well as the polar diameter and polar angle of this row of anchor heads; B l j , r l j , and θ l j are the breadth for the 3D section of the two-stage slope at the altitude of the j-th row of anchor heads on the lower step slope, the polar diameter, and the polar angle of this row of anchor heads in the polar coordinate system, separately. The negative value of the anchoring force work indicates that the anchor cables do negative work on the sliding block of the slope and are the stabilizing factors of the slope. Moreover, θ u i in Equation (77) and θ l j in Equation (78) can be obtained by the geometric relationship shown in Figure 6.
θ u i = arctan e ( θ h θ 0 ) tan φ t sin θ h h u i / r 0 cos θ 0 η 2 ( η 1 h u i / r 0 ) cot β 1
θ l j = arctan e ( θ h θ 0 ) tan φ t sin θ h h l j / r 0 e ( θ h θ 0 ) tan φ t cos θ h + h l j / r 0 cot β 2
In combination with the dimensionless quantities f 1 , f 2 , f 4 and f 6 appearing and employed in the preceding chapters, the entire anchorage force work can be further written as
W T = i = 1 n W ˙ u i + j = 1 m W ˙ l j             = ω r 0   4 ( T u i i = 1 n g 7 i + T l j j = 1 m g 8 j )
in which g 7 i   ( i = 1 , 2 , , n ) and g 8 j   ( j = 1 , 2 , , m ) are non-dimensional formulas for the anchorage force work and are listed as
g 7 i = b H η 1   3 + 2 η 1   2 f 2 ( θ u i ) 2 f 4 ( θ u i ) 2 H 2 f 1 ( θ u i ) + f 4 ( θ u i ) sin ( θ u i δ u i )
g 8 j = b H η 1   3 + 2 η 1   2 f 2 ( θ l j ) 2 f 6 ( θ l j ) 2 H 2 f 1 ( θ l j ) + f 6 ( θ l j ) sin ( θ l j δ l j )
Hence, the entire external force work would be equated to the totality for the gravitational work, seismic work, and negative work contributed by the anchor cable groups, which is
W = W γ + W s + W T

3.4. Analytical Solution for the Factor of Safety

According to the upper bound method for LA, the factor of safety (FS) for the 3D two-stage slope is the rate for internal energy dissipation to external force work, which can be deduced as
F S = D W γ + W s + W T
Taking Equations (50), (54), (74), and (81) into the above formula and introducing the stability coefficient λ c = c / ( γ H ) , the analytical solution for FS can ultimately be organized as
F S = η 1 λ c cos φ t ( g 1 + g 2 ) g 3 + g 4 + k h ( g 5 + g 6 ) γ 1 ( T u i i = 1 n g 7 i + T l j j = 1 m g 8 j )
From the above formula, it can be observed that the slope is in the critical state of instability when FS = 1. Thus, FS > 1 is the necessary condition for slope stability. And it is known from the limit analysis theorem that FS in Equation (86) is the upper bound of the real FS. Additionally, the function on FS given in Equation (86) involves five variables altogether: θ 0 , θ h , φ t , r 0 / r 0 and t , whose definition domains are represented by the following set of inequalities.
0 < θ 0 < θ B < θ C θ D < θ h < π 0 < φ t < 2 / π 0 < r 0 / r 0 < 1 0 < ξ 1 < 1 0   <   t   <   T 0 < b / H + B max / H B / H
where B max is the maximum width of the three-dimensional part of the two-stage slope. With θ 0 , θ h , φ t , r 0 / r 0 , and t as optimization variables and FS as the target variable, a global optimization program based on the Particle Swarm Optimization (PSO) was designed and executed to search for the minimum value of the feasible solution for FS.

4. Results and Discussions

In the following, the results are compared with those presented in the previous literature to demonstrate the validity of the 3D two-stage rock slope model reinforced with prestressed anchor cables governed by the Hoek–Brown criterion and the accuracy of the corresponding analytical solutions. Afterwards, parametric analyses of seismic pseudo-dynamic, prestressing anchorage, rock mass characterization, and geometric features for the two-stage slope are given to explore their respective independent effects on slope safety and their coupling effects on each other.

4.1. Comparisons

Another assessment measure for slope stability, the stability number N s ( N s = γ H c r / c t ), is utilized to analyze the stability of the slopes governed by the linear M-C yield criterion in the 3D horn-shaped rotating failure mechanism [9], in which H c r is the critical height of the slope in the case of instability and the corresponding safety factor FS = 1 . Meanwhile, the safety number N s can be represented by the reciprocal of λ c t ( N s = 1 / λ c t ). In Equation (1), if the parameter n = 1, the nonlinear H-B criterion will be transformed into a linear M-C criterion. Hence, to verify the correctness of the rock slope failure mechanism, the results from Michalowski and Drescher can be converted into λ c t and brought into Equation (1) to calculate FS, corresponding to T u i = 0.0   ( i = 1 ,   2 ,     ,   n ) , T l j = 0.0   ( j = 1 ,   2 ,     ,   m ) , k h = 0.0 , and λ a = 0.0 ( λ a = a / H , the step width coefficient). The comparison results are displayed in Table 1, which shows that the safety factor is close to 1, indicating that the failure model proposed in this paper is effective.
In order to examine the accuracy of the modified pseudo-dynamic approach used in the article, it is possible to compare the results with those obtained from the pseudo-static method, whose horizontal acceleration can be expressed as: a h = k h g . After eliminating the effect of the remaining parameters, the comparison outcomes employing the two methods are listed in Table 2, corresponding to f = 1.0 and D = 0.0 . It can be observed that the FS derived from the two pseudo-dynamic methods are remarkably analogous. The rest of the input parameter values are: H = 15   m , β 1 = β 2 = 45 ° , λ a = 0.0 , σ c = 10 Mpa , γ = 25 kPa , GSI = 20 , m i = 10 , D = 0.0 , T = 0.2   s , V s = 300   m / s , T u i = 0.0   ( i = 1 ,   2 ,     ,   n ) , and T l j = 0.0   ( j = 1 ,   2 ,     ,   m ) .
On the basis of converting the nonlinear H-B criterion to the linear M-C criterion, set λ a = 0.0 , B / H = , and the number of anchor rows n = 1, m = 0 for the upper and lower step of the slope, respectively, and the 3D two-stage slope reinforced by multiple rows of prestressed anchors is reduced to the 2D single-stage slope reinforced by a single row of anchors, which will facilitate a comparison with the data in the previous literature [36,49]. From the comparison results in Table 3, it is clear that the data derived in this paper employing the LA in the 3D horn-shaped rotating failure mechanism are similar to the results obtained using other methods, which correspond to the rest of the geometric parameters: H = 8   m , h u i = 4   m , δ u i = 15 ° .

4.2. Parametric Effects Discussions

According to the previous sections of this paper, the following aspects affect the safety of 3D rock slopes: nonlinear H-B criterion parameters, seismic pseudo-dynamic parameters, anchorage parameters, and slope geometry parameters, whose effects on rock slopes are plotted in Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11, respectively. The rows of upper and lower step anchor cables are: n = 1 , m = 1 . And the corresponding fundamental parameters are: h u i = 9 , h l j = 3 , δ u i = δ l j = 15 , T u i = T l j = 40 kN , H = 12   m , σ c = 10 Mpa , and γ = 25 kPa .

4.2.1. Nonlinear H-B Criterion Parametric Analysis

As shown in Figure 7, the Geological Strength Index G S I and the parameter for the class of rock substance m i are proportional to FS, i.e., the larger that G S I and m i are, the higher the stability of the slope, where the FS rise is more noticeable in Figure 7a. In the G S I rating system, the value of G S I is determined by the structure of the rock mass and the quality condition of the structural surface,
Under the rating system, the quality conditions of the structural surface influence the value of G S I , from which it can be observed that the overall quality of the rock mass has more of an impact on the stability of rock slopes than the quality characteristics of the rock in the stability analysis of the rock slopes. Thus, special attention should be paid to both the structure of the rock mass and the quality condition of the structural face in the engineering design.
Furthermore, the fall in FS with an increase in B/H occurs when nonlinear H-B criteria factors are maintained. This indicates that, given the same other factors, 3D slopes are more stable than 2D slopes. This shows that it is overly cautious to simply handle the 3D slope in actual engineering as if it were just a 2D case. Additionally, Figure 7c shows the trend of FS with m i for different values of G S I . And it can be noticed that the FS growth amplitude is also small when the G S I is small, implying that the rock type and quality do not significantly enhance the stability of the rock slopes when the joint fractures in the rock mass are excessively developed.
Figure 7. Impacts of nonlinear H-B criterion parameters on the safety of 3D two-stage rock slopes: (a) the Geological Strength Index G S I ; (b) the parameter for the class of rock substance m i ; (c) the interactions between G S I and m i .
Figure 7. Impacts of nonlinear H-B criterion parameters on the safety of 3D two-stage rock slopes: (a) the Geological Strength Index G S I ; (b) the parameter for the class of rock substance m i ; (c) the interactions between G S I and m i .
Applsci 14 04160 g007

4.2.2. Improved Pseudo-Dynamic Parametric Analysis

In Figure 8, the effects of several typically improved pseudo-dynamic factors on the safety for two-stage rock slopes as B/H varies are depicted. It can be seen that FS gradually decreases as the horizontal seismicity factor k h , velocity of transverse wave propagating V s , seismic period T , and damping ratio D r increase, respectively. But their reductions are not linear, i.e., it is as the value of the seismic pseudo-dynamic parameters change at smaller values that FS decreases more. For example, fixing B/H = 3, FS reduces by 35% as k h increases from 0.1 to 0.2, while FS decreases by less than 10% as k h increases from 0.4 to 0.5, as shown in Figure 8a. The findings demonstrate that the smaller magnitude of aftershocks is additionally significant for the slope reinforcement engineering. Then, Figure 9 exhibits the interactions of the horizontal seismicity coefficient k h with the rest of the modified pseudo-dynamic parameters. At arbitrary k h ( 0.1 k h 0.5 ), the variations of FS with the unit change for V s , T , and D r are all minor and do not exceed 10%.
Figure 8. Effects of pseudo-dynamic parameters: (a) the horizontal seismicity factor k h ; (b) velocity of transverse wave propagating V s ; (c) seismic period T ; (d) damping ratio D r .
Figure 8. Effects of pseudo-dynamic parameters: (a) the horizontal seismicity factor k h ; (b) velocity of transverse wave propagating V s ; (c) seismic period T ; (d) damping ratio D r .
Applsci 14 04160 g008
Figure 9. The interactions between the seismic pseudo-dynamic parameters of the safety of 3D two-stage rock slopes: (a) the reciprocal actions between V s and k h ; (b) the reciprocal actions between T and k h ; (c) the reciprocal actions between D r and k h .
Figure 9. The interactions between the seismic pseudo-dynamic parameters of the safety of 3D two-stage rock slopes: (a) the reciprocal actions between V s and k h ; (b) the reciprocal actions between T and k h ; (c) the reciprocal actions between D r and k h .
Applsci 14 04160 g009aApplsci 14 04160 g009b

4.2.3. Slope Geometry Parametric Analysis

The influence of the slope geometry parameters on slope safety at different seismic classes is illustrated in Figure 10. It can be observed that the height coefficient ξ 1 and the factor of stage width λ a are favorably proportional to FS, while β 1 and β 2 are inversely proportional to FS. And it is noteworthy that the design pattern for the two-stage slope could significantly withstand the adverse effects of seismicity on slope safety, as shown in Figure 10a, and its enhancement on FS is more obvious when the seismicity is at a smaller level. Moreover, the location of the additional step on the middle and lower portions of the slope can better improve the stability of the slope, as demonstrated in Figure 10b. It is important to note that, at the same range of 35°–60°, the FS fluctuation magnitude relating to β 2 is greater than that of β 1 , which implies that β 2 has a stronger impact on slope safety. And the above findings have some reference to the design for the two-stage rock slope.
Figure 10. Effects of slope geometry parameters: (a) the factor of stage width λ a ; (b) the height factor ξ 1 ; (c) inclined angle of upper stage β 1 ; (d) inclined angle of lower step β 2 .
Figure 10. Effects of slope geometry parameters: (a) the factor of stage width λ a ; (b) the height factor ξ 1 ; (c) inclined angle of upper stage β 1 ; (d) inclined angle of lower step β 2 .
Applsci 14 04160 g010aApplsci 14 04160 g010b

4.2.4. Anchorage Parametric Analysis

When the anchorage force for each row of anchor cables is consistent, the critical anchorage force of the prestressed anchor cable can be derived by the variation of Equation (86):
T c r = γ g 3 + g 4 + k h ( g 5 + g 6 ) η 1 λ c cos φ t ( g 1 + g 2 ) i = 1 n g 7 i + j = 1 m g 8 j
The above equation can provide the reference or verification basis for the anchorage design of rock slopes. Moreover, Figure 11 depicts the effects of prestressing anchorage parameters on the reinforcement of two-stage rock slopes with the variation in seismicity classes. It can be found that FS gradually reduces with the increase in the anchorage angles δ u i and δ l j and the anchor head heights h u i and h l j , while the slope stability is progressively enhanced with the growth in the anchoring forces T u i and T l j and the number of anchor rows n and m. It is worth noting that the anchor head heights have little effect on the stability of the slope at small seismicity magnitudes, whereas the effect of the anchor cable arrangement at different heights on FS varies greatly at larger seismicity magnitudes, as illustrated in Figure 11b. In addition, Figure 11a shows that the slippage of the sliding body for the slope can be more effectively controlled with the anchorage angles close to horizontal. These indicate that the anchorage design is essential for the seismic stability of slopes.
Figure 11. Impacts of anchorage parameters on the safety of 3D two-stage rock slopes: (a) the anchorage angles δ u i and δ l j ; (b) the anchor head heights h u i and h l j ; (c) the number of anchor rows n and m, the anchoring forces T u i and T l j .
Figure 11. Impacts of anchorage parameters on the safety of 3D two-stage rock slopes: (a) the anchorage angles δ u i and δ l j ; (b) the anchor head heights h u i and h l j ; (c) the number of anchor rows n and m, the anchoring forces T u i and T l j .
Applsci 14 04160 g011

5. Conclusions

According to the upper bound theorem of LA, this paper presents the new mathematical formula for the analytical solution of the seismic safety of the 3D two-stage rock slope reinforced with prestressed anchor cables governed by the nonlinear Hoek–Brown criterion, in which the analyses of seismicity are performed by the latest modified pseudo-dynamic method. And the 3D failure mechanism is employed to analyze and derive internal energy dissipation and external forces work within the failure model of two-stage rock slopes, where the power of the seismic force is calculated by the horizontal stratification method, and the work exerted by the anchor groups on the upper and lower stage surface is derived and calculated, separately. With these prerequisites, the factor of safety (FS) defined by the rate for internal energy dissipation and external force work can be collated as an explicit function of the five variables. And the precise analytical solution of FS can be found by combining it with the PSO algorithm.
To verify the correctness of the assessment model for the seismic safety for the 3D two-stage rock slope reinforced with prestressing anchor cables and the accuracy of the corresponding analytical solutions presented in this paper, the results were compared with those of the pertinent literature and other approaches, and the errors were all within the acceptable range. Subsequently, the sensitivity analyses and the coupling effects of the seismic pseudo-dynamic parameters and prestressing anchorage parameters within the kinematic model of the 3D two-stage rock slopes governed by the Hoek–Brown criterion were carried out. And the key findings are shown below.
(1)
For the rock mass characteristics, the analyses indicate that the overall quality of the rock mass has more of an impact on the stability of rock slopes than the quality characteristics of the rock in the stability analysis of the rock slopes. And the rock type and quality do not significantly enhance the stability of the rock slopes when the joint fractures in the rock mass are excessively developed.
(2)
For the seismicity factors, their variations with FS are not linear. And it is as the value of the seismic pseudo-dynamic parameters change at smaller values that FS decreases more, which demonstrates that the smaller magnitude of seismicity is also additionally significant for the slope reinforcement engineering.
(3)
For the slope geometry features, the design pattern for the two-stage slope could significantly withstand the adverse effects of seismicity on slope safety, and its enhancement on FS is more obvious when the seismicity is at a smaller level. Moreover, the location of the additional step on the middle and lower portions of the slope can better improve the stability of the slope.
(4)
Equation (87) can provide the reference or verification basis for the anchorage design of rock slopes. For the prestressing anchorage parameters, the slope stability is progressively enhanced with the growth of the anchoring forces T u i and T l j and the number of anchor rows n and m. Additionally, the slippage of the sliding body for the slope can be more effectively controlled with the anchorage angles close to horizontal. These indicate that the anchorage design is essential for the seismic stability of slopes.

Author Contributions

Methodology, Y.Y.; resources, D.Z. and H.L.; data curation, Y.Y.; writing—original draft, Y.Y.; writing—review and editing, D.Z. and H.L.; supervision, D.Z. and J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

The Changzhou Science and Technology Support Programme: Study on Key Technologies for the Prevention and Control of Surface Deformation Induced by the Construction and Operation of the Changzhou Metro Tunnel (CE20235048) funded this study.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The analysis data are available from the corresponding author on reasonable request. The data are not publicly available due to privacy.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Jiang, S.H.; Li, D.Q.; Zhang, L.M.; Zhou, C.B. Time-dependent system reliability of anchored rock slopes considering rock bolt corrosion effect. Eng. Geol. 2004, 175, 1–8. [Google Scholar] [CrossRef]
  2. Blanco-Fernandez, E.; Castro-Fresno, D.; Díaz, J.D.C.; Lopez-Quijada, L. Flexible systems anchored to the ground for slope stabilisation: Critical review of existing design methods. Eng. Geol. 2011, 122, 129–145. [Google Scholar] [CrossRef]
  3. Eskandarinejad, A. Seismic vertical uplift capacity of horizontal strip anchors embedded in sand adjacent to slopes using finite element limit analysis. Arab. J. Geosci. 2022, 15, 16. [Google Scholar] [CrossRef]
  4. Sun, J.; Yu, T.; Dong, P. Pseudo-dynamic analysis of reinforced slope with anchor cables. Soil Dyn. Earthq. Eng. 2022, 162, 107514. [Google Scholar] [CrossRef]
  5. Tao, Z.; Zhu, C.; He, M.; Karakus, M. A physical modeling-based study on the control mechanisms of Negative Poisson’s ratio anchor cable on the stratified toppling deformation of anti-inclined slopes. Int. J. Rock Mech. Min. Sci. 2021, 138, 104632. [Google Scholar] [CrossRef]
  6. Chen, J.F.; Du, C.C.; Peng, M.; Sun, R.; Zhao, F.; Shi, Z.M. System reliability analysis of a slope stabilized with anchor cables and piles under seismic loading. Acta Geotech. 2023, 18, 4493–4514. [Google Scholar] [CrossRef]
  7. Zhang, Z.L.; Yang, X.L. Unified solution of safety factors for three-dimensional compound slopes considering local and global instability. Comput. Geotech. 2023, 155, 105227. [Google Scholar] [CrossRef]
  8. Zhong, J.H.; Yang, X.L. Seismic stability of three-dimensional slopes considering the nonlinearity of soils. Soil Dyn. Earthq. Eng. 2021, 140, 106334. [Google Scholar] [CrossRef]
  9. Zhang, S.; Xu, Q.; Peng, D.; Zhu, Z.; Li, W.; Wong, H.; Shen, P. Stability analysis of rock wedge slide subjected to ground water dynamic evolution. Eng. Geol. 2020, 270, 105528. [Google Scholar] [CrossRef]
  10. Xu, S.; Zhou, D. An Analytical Framework for Assessing the Unsaturated Bearing Capacity of Strip Footings under Transient Infiltration. Mathematics 2023, 11, 3480. [Google Scholar] [CrossRef]
  11. Kang, X.D.; Zhou, D. Analytical Solution for Bearing Capacity of Reinforced Strip Footings on Unsaturated Soils under Steady Flow. Mathematics 2023, 11, 3746. [Google Scholar] [CrossRef]
  12. Griffiths, D.V.; Marquez, R.M. Three-dimensional slope stability analysis by elasto-plastic finite elements. Géotechnique 2007, 57, 537–546. [Google Scholar] [CrossRef]
  13. Yang, X.L.; Yin, J.H. Slope stability analysis with nonlinear failure criterion. J. Eng. Mech. 2004, 130, 267–273. [Google Scholar] [CrossRef]
  14. Yang, X.L.; Wang, J.M. Ground movement prediction for tunnels using simplified procedure. Tunn. Undergr. Space Technol. 2011, 26, 462–471. [Google Scholar] [CrossRef]
  15. Yang, X.L.; Huang, F. Collapse mechanism of shallow tunnel based on nonlinear Hoek-Brown failure criterion. Tunn. Undergr. Space Technol. 2011, 26, 686–691. [Google Scholar] [CrossRef]
  16. Giger, M.W.; Krizek, R.J. Stability analysis of vertical cut with variable corner angle. Soils Found. 1975, 15, 63–71. [Google Scholar] [CrossRef]
  17. De Buhan, P.; Garnier, D. Three dimensional bearing capacity analysis of a foundation near a slope. Soils Found. 1998, 38, 153–163. [Google Scholar] [CrossRef] [PubMed]
  18. Michalowski, R.L.; Drescher, A. Three-dimensional stability of slopes and excavations. Géotechnique 2009, 59, 839–850. [Google Scholar] [CrossRef]
  19. Zhong, J.H.; Yang, X.L. Pseudo-dynamic stability of rock slope considering Hoek-Brown strength criterion. Acta Geotech. 2021, 17, 2481–2494. [Google Scholar] [CrossRef]
  20. Hoek, E.; Brown, E.T. Empirical strength criterion for rock masses. J. Geotech. Eng. Div. 1980, 106, 1013–1035. [Google Scholar] [CrossRef]
  21. Hoek, E.; Wood, D.; Shah, S. A modified Hoek–Brown failure criterion for jointed rock masses. In Rock Characterization: ISRM Symposium, Eurock’92, Chester, UK, 14–17 September 1992; Thomas Telford Publishing: London, UK, 1992; pp. 209–214. [Google Scholar]
  22. Hoek, E.; Carranza-Torres, C.; Corkum, B. Hoek-Brown failure criterion-2002 edition. Proc. NARMS-Tac 2002, 1, 267–273. [Google Scholar]
  23. Yang, X.L.; Yin, J.H. Upper bound solution for ultimate bearing capacity with a modified Hoek-Brown failure criterion. Int. J. Rock Mech. Min. Sci. 2005, 42, 550–560. [Google Scholar] [CrossRef]
  24. Yang, X.L.; Li, L.; Yin, J.H. Seismic and static stability analysis for rock slopes by a kinematical approach. Géotechnique 2004, 54, 543–549. [Google Scholar] [CrossRef]
  25. Loukidis, D.; Bandini, P.; Salgado, R. Stability of seismically loaded slopes using limit analysis. Géotechnique 2003, 53, 463–479. [Google Scholar] [CrossRef]
  26. Brink, U.T.; Andrews, B.; Miller, N. Seismicity and sedimentation rate effects on submarine slope stability. Geology 2016, 44, 563–566. [Google Scholar] [CrossRef]
  27. Macedo, J.; Candia, G.; Lacour, M.; Liu, C. New developments for the performance-based assessment of seismically-induced slope displacements. Eng. Geol. 2020, 277, 105786. [Google Scholar] [CrossRef]
  28. Michalowski, R.L.; Martel, T. Stability charts for 3D failures of steep slopes subjected to seismic excitation. J. Geotech. Geoenviron. Eng. 2011, 137, 183–189. [Google Scholar] [CrossRef]
  29. Li, Z.W.; Yang, X.L.; Li, T.Z. Static and seismic stability assessment of 3D slopes with cracks. Eng. Geol. 2020, 265, 105450. [Google Scholar] [CrossRef]
  30. Steedman, R.S.; Zeng, X. The influence of phase on the calculation of pseudo-static earth pressure on a retaining wall. Géotechnique 1990, 40, 103–112. [Google Scholar] [CrossRef]
  31. Zeng, X.; Steedman, R.S. On the behaviour of quay walls in earthquakes. Géotechnique 1993, 43, 417–431. [Google Scholar] [CrossRef]
  32. Bellezza, I. A new pseudo-dynamic approach for seismic active soil thrust. Geotech. Geol. Eng. 2014, 32, 561–576. [Google Scholar] [CrossRef]
  33. Bellezza, I. Seismic active earth pressure on walls using a new pseudo-dynamic approach. Geotech. Geol. Eng. 2015, 33, 795–812. [Google Scholar] [CrossRef]
  34. Chanda, N.; Ghosh, S.; Pal, M. Analysis of slope using modified pseudo-dynamic method. Int. J. Geotech. Eng. 2018, 12, 337–346. [Google Scholar] [CrossRef]
  35. Shi, K.; Wu, X.; Liu, Z.; Dai, S. Coupled calculation model for anchoring force loss in a slope reinforced by a frame beam and anchor cables. Eng. Geol. 2019, 260, 105245. [Google Scholar] [CrossRef]
  36. Li, X.; He, S.; Wu, Y. Limit analysis of the stability of slopes reinforced with anchors. Int. J. Numer. Anal. Methods Geomech. 2012, 36, 1898–1908. [Google Scholar] [CrossRef]
  37. Yan, M.; Xia, Y.; Liu, T.; Bowa, V.M. Limit analysis under seismic conditions of a slope reinforced with prestressed anchor cables. Comput. Geotech. 2019, 108, 226–233. [Google Scholar] [CrossRef]
  38. Li, A.J.; Lyamin, A.V.; Merifield, R.S. Seismic rock slope stability charts based on limit analysis methods. Comput. Geotech. 2009, 36, 135–148. [Google Scholar] [CrossRef]
  39. Kokane, A.K.; Sawant, V.A.; Sahoo, J.P. Seismic stability analysis of nailed vertical cut using modified pseudo-dynamic method. Soil Dyn. Earthq. Eng. 2020, 137, 106294. [Google Scholar] [CrossRef]
  40. Agarwal, E.; Pain, A. Probabilistic stability analysis of geosynthetic-reinforced slopes under pseudo-static and modified pseudo-dynamic conditions. Geotext. Geomembr. 2021, 49, 1565–1584. [Google Scholar] [CrossRef]
  41. Xu, S.; Zhou, D. Seismic Bearing Capacity Solution for Strip Footings in Unsaturated Soils with Modified Pseudo-Dynamic Approach. Mathematics 2023, 11, 2692. [Google Scholar] [CrossRef]
  42. Zhang, Z.L.; Yang, X.L. Seismic stability analysis of slopes with cracks in unsaturated soils using pseudo-dynamic approach. Transp. Geotech. 2021, 19, 100583. [Google Scholar] [CrossRef]
  43. Qin, C.B.; Chian, S.C. Kinematic analysis of seismic slope stability with a discretisation technique and pseudo-dynamic approach: A new perspective. Géotechnique 2018, 68, 492–503. [Google Scholar] [CrossRef]
  44. Gazetas, G.; Garini, E.; Anastasopoulos, I.; Georgarakos, T. Effects of near-fault ground shaking on sliding systems. J. Geotech. Geoenviron. Eng. 2009, 135, 1906–1921. [Google Scholar] [CrossRef]
  45. Yang, Y.; Zhou, D. Seismic Stability for 3D Two-Step Slope Governed by Non-Linearity in Soils Using Modified Pseudo-Dynamic Approach. Appl. Sci. 2022, 12, 6482. [Google Scholar] [CrossRef]
  46. Yang, Y.; Liao, H.; Zhu, J. Stability Analysis of Three-Dimensional Tunnel Face Considering Linear and Nonlinear Strength in Unsaturated Soil. Appl. Sci. 2024, 14, 2080. [Google Scholar] [CrossRef]
  47. Chen, G.; Chen, T.; Chen, Y.; Huang, R.; Liu, M. A new method of predicting the prestress variations in anchored cables with excavation unloading destruction. Eng. Geol. 2018, 241, 109–120. [Google Scholar] [CrossRef]
  48. Cheng, Y.; Xu, Y.; Wang, L.; Wang, L. Stability of expansive soil slopes reinforced with anchor cables based on rotational-translational mechanisms. Comput. Geotech. 2022, 146, 104747. [Google Scholar] [CrossRef]
  49. Qiang, L. Stability Study for Rock Slope and Anchoring Parametric Analysis. Ph.D. Thesis, Central South University, Changsha, China, 2010. [Google Scholar]
Figure 1. The envelope curve of the nonlinear H-B yield criterion and its outer tangent.
Figure 1. The envelope curve of the nonlinear H-B yield criterion and its outer tangent.
Applsci 14 04160 g001
Figure 2. Three-dimensional horn-shaped rotating model for stepped rock slope reinforced with anchor groups.
Figure 2. Three-dimensional horn-shaped rotating model for stepped rock slope reinforced with anchor groups.
Applsci 14 04160 g002
Figure 3. Stereoscopic diagram for 3D failure mechanism with the plane strain block and anchor group reinforcement.
Figure 3. Stereoscopic diagram for 3D failure mechanism with the plane strain block and anchor group reinforcement.
Applsci 14 04160 g003
Figure 4. Schematic of the derivation of internal energy dissipation with local coordinate systems.
Figure 4. Schematic of the derivation of internal energy dissipation with local coordinate systems.
Applsci 14 04160 g004
Figure 5. Diagram of the horizontal stratification method: (a) full graph; (b) detail zoom.
Figure 5. Diagram of the horizontal stratification method: (a) full graph; (b) detail zoom.
Applsci 14 04160 g005
Figure 6. Schematic layout of anchor groups on the upper step slope and lower step slope.
Figure 6. Schematic layout of anchor groups on the upper step slope and lower step slope.
Applsci 14 04160 g006
Table 1. Comparisons between this paper’s findings and Michalowski and Drescher’s solutions.
Table 1. Comparisons between this paper’s findings and Michalowski and Drescher’s solutions.
B/HSolutions β 1 = β 2
45 ° 60 ° 75 ° 90 °
1.0Ns form Drescher and Michalowski54.85023.83514.70111.028
FS using the current approach0.98940.98240.99920.9848
1.5Ns form Drescher and Michalowski46.84520.77312.9768.935
FS using the current approach0.99210.97420.97550.9927
2.0Ns form Drescher and Michalowski42.73219.10312.1098.604
FS using the current approach0.99950.99180.97950.9599
3.0Ns form Drescher and Michalowski39.95617.87311.1847.974
FS using the current approach1.00101.00110.99560.9924
5.0Ns form Drescher and Michalowski37.99417.06310.6287.266
FS using the current approach1.00621.00320.99850.9938
10.0Ns form Drescher and Michalowski36.70316.52710.2656.944
FS using the current approach1.00861.00251.00131.0003
Table 2. Comparisons of the original pseudo-static results with the amended pseudo-dynamic results.
Table 2. Comparisons of the original pseudo-static results with the amended pseudo-dynamic results.
B/HAmended Pseudo-Dynamic ResultsOriginal Pseudo-Static Results
kh = 0.2kh = 0.5kh = 0.8kh = 0.2kh = 0.5kh = 0.8
1.51.9811.2150.9261.9321.1350.920
3.01.2760.8790.6761.1850.8440.634
5.01.1360.8350.6601.0910.7950.577
10.01.0310.7710.5980.9920.7550.545
2D1.0160.7350.5790.9620.7180.503
Table 3. Comparisons of solutions with other methods for prestressing anchorage.
Table 3. Comparisons of solutions with other methods for prestressing anchorage.
Methods T   ( kN )
040
Finite Element Method/Cai and Ugai [36]1.091.21
Bishop’s Simplified Method/SLOPE/W [36]1.0951.207
Finite Difference Method/FLAC [49]1.111.22
Limit Analysis Method (2D)/Li et al. [49]1.0961.213
Limit Analysis Method (3D)/This paper1.1061.251
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, Y.; Liao, H.; Zhou, D.; Zhu, J. Analytical Solution for Seismic Stability of 3D Rock Slope Reinforced with Prestressed Anchor Cables. Appl. Sci. 2024, 14, 4160. https://doi.org/10.3390/app14104160

AMA Style

Yang Y, Liao H, Zhou D, Zhu J. Analytical Solution for Seismic Stability of 3D Rock Slope Reinforced with Prestressed Anchor Cables. Applied Sciences. 2024; 14(10):4160. https://doi.org/10.3390/app14104160

Chicago/Turabian Style

Yang, Yushan, Hong Liao, De Zhou, and Jianqun Zhu. 2024. "Analytical Solution for Seismic Stability of 3D Rock Slope Reinforced with Prestressed Anchor Cables" Applied Sciences 14, no. 10: 4160. https://doi.org/10.3390/app14104160

APA Style

Yang, Y., Liao, H., Zhou, D., & Zhu, J. (2024). Analytical Solution for Seismic Stability of 3D Rock Slope Reinforced with Prestressed Anchor Cables. Applied Sciences, 14(10), 4160. https://doi.org/10.3390/app14104160

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop