Next Article in Journal
Bayesian Mixture Model to Estimate Freeway Travel Time under Low-Frequency Probe Data
Next Article in Special Issue
Determination of Landslide Displacement Warning Thresholds by Applying DBA-LSTM and Numerical Simulation Algorithms
Previous Article in Journal
Dynamic Hand Gesture Recognition for Smart Lifecare Routines via K-Ary Tree Hashing Classifier
Previous Article in Special Issue
Dynamic Response of a Heavy-Haul Railway Tunnel’s Bottom Structures in Hard Rock
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seismic Stability for 3D Two-Step Slope Governed by Non-Linearity in Soils Using Modified Pseudo-Dynamic Approach

School of Civil Engineering, Central South University, Changsha 410075, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(13), 6482; https://doi.org/10.3390/app12136482
Submission received: 3 June 2022 / Revised: 20 June 2022 / Accepted: 21 June 2022 / Published: 26 June 2022
(This article belongs to the Special Issue Slope Stability and Earth Retaining Structures)

Abstract

:
Seismic events are an active factor in causing slope instability, and the existing pseudo-dynamic method ignores some practical engineering conditions. In this paper, the seismic stability of three-dimensional (3D) two-step slopes is evaluated utilizing the modified pseudo-dynamic method based on the kinematical approach of limit analysis, in which the nonlinear characteristics of soil are considered using the generalized tangent technique. A 3D horn-like rotational failure mechanism is established to solve internal energy dissipation and external force work, in which the seismic work is considered in addition to the soil weight work and obtained by the layer-wise summation method. Based on the force-increase technique, the analytical expression for the safety factor (FS) of 3D two-step slopes is derived more readily. To verify the reliability of the new expression, the present results are compared with already posted solutions and the original pseudo-dynamic solutions. Ultimately, the sensitivity discussions are carried out to investigate the impacts of various factors on slope stability. This has some significance for the design and safety of 3D two-step slopes.

1. Introduction

Stability analysis for slopes has invariably been one of the most classic and hottest subjects in geotechnical engineering. New and numerous innovative approaches have been employed, and feasible engineering conditions are considered in the investigation of slope stability. Most of these research works are directed to slopes in the two-dimensional (2D) plane-strain situation. Nevertheless, slope instability and failure frequently exhibit distinct three-dimensional (3D) properties. Therefore, it is necessary to analyze the 3D stability for slopes in order to obtain more reasonable results. The stability analysis of 3D slopes can normally be carried out by the following three measures: (1) limit equilibrium method (LEM), (2) finite element method (FEM), and (3) limit analysis method (LAM). The existing 3D LEM is generally promoted on the basis of the 2D method; that is, the 2D slices are enlarged to the columns under the 3D space [1,2]. The assumptions of internal forces must be introduced in the analysis so as to ascertain the positive stress distribution on the sliding surface, yet the results obtained from this calculation and the actual stress state often have a large difference [3]. In numerical simulations based on the FEM, the strength reduction method is commonly utilized for calculations [4,5,6]. However, after discounting the slope strength, the stress state inside the slope will be changed so that the location of the most critical sliding surface is somewhat different from the actual situation. Recently, the LAM has been widely implemented in the stability assessment of slopes, which avoids the unreasonable assumptions on internal forces in the LEM and the tedious calculations in the FEM [7,8,9,10].
The primary challenge of applying the LAM to stability evaluation for 3D slopes is to construct a 3D kinematical admissible failure mechanism, and there has been a large amount of literature written on this field of research. The failure mechanisms adopted in the investigation of 3D limit analysis for slope stability can be summarized into two categories: (1) the block sliding mechanism, and (2) the integral sliding mechanism. The multi-block 3D sliding mechanism and the improved multi-block 3D sliding mechanism [11] belong to the first category, and the 3D horn-like rotational mechanism [12] and spherical mechanism belong to the second category. The application of the above two categories of sliding mechanisms in the stability assessment for 3D slopes can obtain more satisfactory results. However, the integral sliding mechanism can ensure calculation accuracy while greatly simplifying the calculation process.
Relevant studies have shown that slope geometry has a significant effect on slope stability. Therefore, in practical excavation engineering, the single-step is often modified into the double-step by adding a step platform to improve slope stability. In addition, earthquakes are also a factor that cannot be ignored in the stability assessment for 3D slopes [13]. The seismic stability of the 3D steep slope was investigated on the basis of the pseudo-static method. Nevertheless, the pseudo-static method evaluates the seismic effect in a highly simplified way, and it does not account for the influence of seismic excitation features on slope stability. A simplified pseudo-dynamic approach was proposed by Steedman and Zeng [14,15], which considers the temporal effects of seismicity but overlooks critical factors such as boundary conditions and material damping. On the basis of the existing methods, Bellezza [16,17] presented a modified pseudo-dynamic method, which has been frequently utilized in seismic stability analysis of geotechnical structures.
Furthermore, the linear Mohr-Coulomb (M-C) criterion, which has previously commonly been used in the stability assessment of soil slopes, cannot reflect the mechanical properties of soil material. Numerous experimental data show that the nonlinear envelope can better depict the shear strength of soils, but it is difficult to be applied in the stability evaluation of soil slopes based on the kinematic method of limit analysis. The generalized tangential technique proposed by Yang and Yin [8] succinctly solves this difficulty and is broadly employed in stability problems of geotechnical structures.
In this paper, the seismic stability of 3D two-step slopes is evaluated using the modified pseudo-dynamic approach based on the kinematical approach of limit analysis, in which the nonlinear characteristic of the soil is considered using the generalized tangent technique. The 3D horn-like rotational failure mechanism is discretized to solve the internal dissipation and the external work, in which the seismic work is considered in addition to the soil weight work and obtained by the layer-wise summation method. Based on the force-increase technique, the safety factor of 3D two-step slopes can be derived more readily. Eventually, the paper verifies the reliability of the results by comparing with already posted solutions and the original pseudo-dynamic solutions, and analyzes and discusses the effects of parameters.

2. Generalized Tangential Technique and Modified Pseudo-Dynamic Method

2.1. Generalized Tangential Technique

In conventional stability analysis for slopes, the linear M-C yield criterion is widely used to describe the shear strength of geotechnical materials. However, the relationship between shear strength and positive stress during actual material failure is often nonlinear [18,19]. The corresponding nonlinear strength envelope is shown in Figure 1, and the analytical equation in the σn-τ stress space is
τ = c 0 ( 1 + σ n / σ t ) 1 / m
where σ n and τ are the normal stress and the shear stress on the failure surface, correspondingly. c 0 is the initial cohesion, σ t is the axial tensile strength at failure, and m is the parameter evaluating the nonlinear feature of the envelope ( m 1 ). It can be seen that for m = 1 , Equation (1) transforms into the classical linear M-C criterion. In addition, φ t can be introduced as an intermediate variable defined as tan φ t = d τ / d σ n , then σ n and τ can be written as the function of φ t
σ n = σ t ( m σ t tan φ t c 0 ) m 1 m σ t
τ = c 0 ( m σ t tan φ t c 0 ) 1 1 m
Nonetheless, the upper limit analysis of the bearing capacity or stability of geotechnical structures on rock under the nonlinear failure criterion has been one of the difficult problems in geotechnical engineering. In this paper, the generalized tangent technique [8] is adopted to obtain the upper limit solution for safety factor FS of 3D slopes subjected to non-linearity in soils. As shown in Figure 1, L is the tangential straight line to any point P on the nonlinear envelope, and its expression can be given as
τ = c t + σ n tan φ t
where c t is the intercept of that tangent line on the axis, bringing Equation (2) and Equation (3) into the Equation (4) whose analytic expression can be deduced
c t = m 1 m c 0 ( m σ t tan φ t c 0 ) 1 1 m + σ t tan φ t
The difficulty of using the nonlinear strength criterion directly in geotechnical engineering problems is avoided by the generalized tangent technique, and instead, the line tangent to the nonlinear envelope is used to describe the shear strength of the soils. Obviously, when the normal stresses are the same, the strength represented by the tangent line is greater than or equal to the strength of the nonlinear envelope. Therefore, the minimum upper limit solution of the safety factor for the 3D slope can be obtained by means of the generalized tangent line technique.

2.2. The Modified Pseudo-Dynamic Method

Seismicity can significantly increase the inducing forces of slope failure and reduce the shear strength of soils, which is one of the main factors inducing slope instability and landslide. Numerous new attempts and advances have been made in the aspect of seismic stability for geotechnical structures. Whereas the pseudo-static method, which is usually used to deliberate the effects of seismicity on 3D slope stability, neglects the fluctuation characteristics of seismic waves that vary with time, the existing pseudo-dynamic method only considers the time effect, ignoring some practical engineering conditions. Therefore, a new pseudo-dynamic method [16,17], which is improved on the basis of the original pseudo-dynamic method, is adopted in this paper. The present study puts forward a kinematic analytical solution for assessing a seismic 3D two-step slope subjected to the non-linearity in soils, which couples the modified pseudo-dynamic method with the upper limit theorem of limit analysis.
Bellezza [16] came up with the modified pseudo-dynamic method based on the assumptions of a visco-elastic medium and Kelvin–Voigt material, which surmounted the imperfections of the existing pseudo-dynamic method, namely: (1) The existing method cannot satisfy the zero-stress boundary condition because it only considers the upward propagation of the incident wave throughout the linear elastic soil medium [16,20]; (2) The acceleration value derived from the existing pseudo-dynamic method is considered to be augmented linearly toward the ground. Hence, an assumed amplification factor is required when using this method [21], which is regarded as equal to 1 in most relevant studies; (3) The existing method does not consider the energy dissipation or damping generated by the real material.
When kv 0.5 kh, Gazetas et al. [22] concluded that the effect of vertical seismic acceleration is almost negligible for the stability of 3D slopes based on the research of horizontal and vertical seismic waves. In order to optimize the program design process and computational speed, the dominant horizontal seismic motivator is taken into account with the least possible influence on the seismic effect in this paper. The horizontal seismic acceleration based on the improved pseudo-dynamic can be given as
a h ( z , t ) = k h g C S 2 + S S 2 [ ( C S C S Z + S S S S Z ) cos ( ω t ) + ( S S C S Z + C S S S Z ) 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 )
where kh and g are horizontal seismic coefficient and gravitational acceleration, respectively, and ys1 and y s 2 are defined as a function of the normalized frequency wH/Vs and damping ratio D, in which H is the height of slope, and Vs is the velocity of S-wave propagating. The remaining parameters: w = 2π/T is angular frequency, where T is the seismic period; t is the arbitrary moment in the period T; z is the value corresponding to the vertical axis.

2.3. Force-Increase Technique

The resistance and the induction act together on the slope structure, and the slope structure will be destabilized or collapsed as the induction is critical to the resistance. In the stability analysis of slopes, two methods are frequently used: (1) lowering the resistance to make the slope critical, namely, the strength reduction method (SRM) [22,23,24]; (2) raising the induction to obtain the slope critical as well, namely, the force-increase technique (FIT) [25,26] used in this paper.
The SRM has universally been utilized to measure the safety factor FS of slopes, but what is derived from the SRM is an implicit equation containing FS. Thus, the force-increasing technique is used in this paper to facilitate the program design and subsequent solution. According to the force-increase principle, FS, namely, the factor of safety, can be customarily defined as
F S = γ c r γ
where γ and γ c r are the respective initial unit weight of soil and the unit weight of soil that increases when the slope is in critical failure. It is well established that the external work can be written as the linear expression on the external force. According to this, FS can be correspondingly expressed as the ratio of the critical external work to the real external work
F S = W c r W
where W c r is the critical external work, and W is the real external work. On the basis of the kinematical approach of limit analysis, W c r = D (internal energy dissipation), so FS can be further written as the ratio between internal energy dissipation and real external force work, i.e.,
F S = D W
As we can see from Equation (15), the equation on FS derived by the force increase technique is the explicit function, which greatly facilitates solutions of complex integral expressions obtained from both real external force work and internal energy dissipation in 3D conditions.

3. Kinematical Analysis Model of 3D Soil Slopes with Two Steps

The stability evaluation problems for geological engineering are frequently greatly simplified by adopting the kinematical approach of the LAM [27]. The kinematical approach can be formulated as follows: for any kinematically admissible plastic strain rate fields and velocity fields, if the internal energy dissipation is equated with the external force work, then the calculated solutions must be higher or identical to what actually occurs. Additionally, since the associated flow law must be satisfied in the kinematically admissible velocity field, which corresponds to the M-C criterion, there is a relationship for the principal strain rates, namely,
ε ˙ 1 + ε ˙ 2 + ε ˙ 3 + ( ε ˙ 1 ε ˙ 2 ε ˙ 3 ) sin φ t = 0
where φ t is the slope of the tangent line at any point on the nonlinear envelope aforementioned in the section on generalized tangent technique, namely, the angle of internal friction of the soil. Applying the Equation (16) to the velocity interrupted curved surface formed at the time of slope failure, the following equation can be obtained
v n = v t tan φ t
in which v n is the normal component of the velocity along the interrupted curved surface, and v t is the tangential component of the velocity along the interrupted curved surface. Furthermore, the 3D horn-like rotational mechanism, which has been shown to satisfy this condition well [12], is employed in this paper for the establishment of the failure mechanism for 3D two-step slopes.

3.1. Description of the 3D Failure Mechanism in Two-Step Slopes

Constructing the failure mechanism is one of the most crucial procedures in the stability analysis of geotechnical structures using the kinematic method of limit analysis. The 3D two-step slopes subjected to seismicity are analyzed in the framework of the 3D horn-like rotational failure mechanism in this paper.
In Figure 2, it can be observed that the exterior for the 3D two-step slope consists of the top surface for the slope AB, slope surface above the step BC with inclination β 1 , slope surface below the step CD with inclination β 2 , and step plane DE. The overall slope height is H. The height of slope above the horizontal step surface is ξ 1 H , and the height of slope below the horizontal step surface is ξ 2 H , for which ξ 1 and ξ 2 are the coefficients of depth. Apparently,
ξ 1 + ξ 2 = 1
The 3D horn-like rotational failure mechanism is implemented on the 3D two-step slope by crossing the peak and foot for the slope, and the symmetry plane for this failure mechanism is perfectly decided by two logarithmic spirals, whose analytic expressions can be given separately as
r = r 0 e ( θ θ 0 ) tan φ t
and
r = r 0 e ( θ θ 0 ) tan φ t
in which r is the curve AE, r′ is the curve A′E′, and r0 = OA, r0 = OA, and θ 0 is the angle between the horizontal axis and the line OA, as shown in Figure 2. The 3D horn-like rotational failure mechanism is the curvilinear cone with vertex angle 2 φ t shaped by the circle with progressively increasing radius in continuity around the center of rotation O. The length between point O and the center of the circle, and the radius of the circle can be expressed separately as
r m = ( r r ) / 2 = r 0 f 1
and
R = ( r r ) / 2 = r 0 f 2
in which f 1 and f 2 can be written 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 ]
As illustrated in Figure 3, the 3D slopes can be extended to the case of 2D by embedding the block of width b with the plane-strain property in the 3D failure mechanism for two-step slopes, in accordance with Michalowski and Drescher [12]. The width b of the embedded block can be calculated as
b = B B m a x
where B and Bmax are the width for the entire two-step slope and the maximum width for the 3D section of the two-step slope, correspondingly. The 3D two-step slope is transformed into the 2D condition with the ratio of overall width B to height H of the two-step slope tending to infinity.

3.2. Internal Energy Dissipation

The internal energy dissipation inside the 3D two-step slope can be divided into two portions: (1) the dissipation generated along the non-continuous surface St where sliding is about to occur Dt; (2) the dissipation generated by the internal volume deformation within the slope DV. In accordance with the upper limit theorem of limit analysis, Dt and DV can be, separately, given as
D t = S t c t v cos φ t d S t
and
D V = V c t cos φ t ( ε ˙ 1 ε ˙ 2 ε ˙ 3 ) d V
where St is the velocity non-continuous surface, and V is the volume of deformation occurring within the slope. v, namely, the velocity in the failure mechanism of 3D stepped slope, can be calculated by the following Equation (28) deduced from the local Cartesian coordinate systems founded on any circular cross-section of the 3D horn-like rotational failure mechanism with the center of the circular cross-section as the origin, as shown in Figure 2.
v = ( r m + y ) ω
in which y is the coordinate value corresponding to the y-axis in the local coordinate system. The rest of the parameters involved in Equations (26) and (27) are explicitly stated in the previous sections and will not be repeated here. According to the Gaussian dispersion theorem and after the rigorous derivation, the total internal energy dissipation for the sliding mechanism can be ultimately written as
D = D t + D V = c t cot φ t S r v n d S r
where Sr is the remaining surface of the kinematically admissible failure mechanism after removing the velocity discontinuity surface. n is the outward unit normal vector of the surface Sr, and v is the corresponding velocity vector.
In general, the overall internal energy dissipation of the failure mechanism relatively easily deduced from the Equation (29) can be classified as four portions along the profile of the two-step slope: DAB, DBC, DCD, and DDE are the corresponding energy dissipation rates on the planes represented by AB, BC, CD, and DE, respectively, and their 3D portions can be shown 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 upper limit of integration x i * (i = 1, 2, 3, 4) of the second integral of the above double integration Equations can be written as x i * = R 2 d i 2 , in which, according to Figure 2, the expressions of d i can be given as, separately
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 dimensionless expressions are written 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 accordance with the geometric features, the angles θ B , θ C , and θ D in Figure 2 can be deduced 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 dimensionless expressions, shown 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 horizontal line AB’s length at the top of the slope, as illustrated in Figure 2. Consequently, the overall internal dissipation for two-step slopes after adding the plane insert section can be expressed in the following form
D = D AB 3 D + D BC 3 D + D CD 3 D + D DE 3 D + D insert = ω c t cot φ t r 0 3 ( g 1 + g 2 )
where g1 and g2 are dimensionless expressions for the internal dissipation of 3D section and 2D plane-strain section, respectively, given 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

The stability for two-step slopes is affected by the gravitational force of the soil and the horizontal seismic force, in which the effect of vertical seismic force is almost negligible at kv ≤ 0.5 kh, based on the conclusion of Gazetas et al. The soil weight work in the two wings of the two-step slope, namely, the 3D section, can be shown as
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 θ ]
where γ is the unit weight of soil. The upper limit of integration x * of the third integral for the above triple integration Equations can be written as x * = R 2 y 2 . Hence, the overall soil weight work after adding the intermediate 2D plane-strain section can be given as
W γ = W γ 3 D + W γ i n s e r t = ω γ r 0 4 ( g 3 + g 4 )
where g3 and g4 are dimensionless expressions of the soil weight work for 3D section and 2D plane-strain section, respectively, 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 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 ]
In this paper, the effects of seismic force on the stability of two-step slopes are evaluated based on the modified pseudo-dynamic approach. Since the seismic acceleration keeps consistent on the same horizontal plane, the seismic force work can be solved using the layer-wise summation method.
As illustrated in Figure 4a, the planar right-angle coordinate system x-E-z is founded with E as the origin. P is the intersection of velocity discontinuity surface and symmetry plane for the 3D failure mechanism at arbitrary height, while P′ is the intersection of slope surface and symmetry plane for the 3D failure mechanism at the corresponding level. The polar angles of P and P′ are denoted as α and α , respectively. Hence, the value on the z-axis corresponding to the horizontal surface PP’ can be derived as
z ( α ) = r 0 e ( θ h θ 0 ) tan φ t sin θ h r p sin α
in which rP is the distance from point P to point O, and can be written as
r p = r 0 e ( α θ 0 ) tan φ t
Bringing Equation (59) into Equation (6), the horizontal seismic acceleration ah of arbitrary height in the two-step slope can be obtained.
In accordance with the layer-wise summation method, the seismic force work in the volume unit dV formed by the vertical increment dz and the horizontal surface PP′ within the slope failure mechanism, as illustrated in Figure 4b, can be given as
d W s = a h g γ v sin θ d V
in which v is the velocity vector, and dV is the volume unit. According to the geometric relations, the volume units corresponding to the 3D section of two wings and the intermediate 2D plane-strain section inserted in the two-step slope failure mechanism can be deduced as, separately,
d V 3 D = l ( θ ) sin α sin 2 θ r p d θ d z
d V i n s e r t = b sin α sin 2 θ r p d θ d z
l ( θ ) in Equation (62) represents the width of the 3D portion for the stepped slope at any polar angle θ   ( θ 0 θ θ h ) , which can be derived as
l ( θ ) = 2 R 2 d 2
d = r p sin α sin θ r m = r 0 f 12
where f 12 is the dimensionless expression, shown as
f 12 = e ( α θ 0 ) tan φ t sin α sin θ f 1
In addition, the vertical increment dz can be further written as
d z = cos ( α φ t ) cos φ t r p d α  
Hence, the seismic force work in the 3D section of two wings and the intermediate 2D plane-strain section can be listed as the following double integral expressions, respectively
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 i n s e r t = ω γ θ 0 θ h α α b r p 3 sin 2 α sin 2 θ cos ( α φ t ) cos φ t a h g d α d θ
The upper limit of integration α for the second integral in the above two equations is the function of α , whose segmental function expressions are deduced 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
After simplification to Equations (68) and (69), the total seismic force work in the two-step slope failure mechanism can be shown as
W s = W s 3 D + W s insert = ω γ k h r 0 4 ( g 5 + g 6 )
in which g5 and g6 are dimensionless expressions of the seismic force work for the 3D section and the 2D plane-strain section, respectively, given 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 θ
The overall external force work is equal to the sum of the seismic force work and the soil weight work, which is
W = W γ + W s

3.4. Factor of Safety

In order to make the expression of the final derived safety factor more concise, two dimensionless parameters are introduced here, namely, the cohesion coefficient λ c and the step-width coefficient λ a , which are defined as follows
λ c t = c t / ( γ H )
λ a = a / H
In accordance with Equation (3), and combining with the external force work and the internal energy dissipation, the factor of safety FS can be derived as
F S = λ c t κ 1 cot φ t ( g 1 + g 2 ) g 3 + g 4 + k h ( g 5 + g 6 )
Furthermore, the geometric parameters must meet the following conditions in order for the failure mechanism to be kinematically admissible.
{ 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
After specifying the remaining parameters, an optimization program can be utilized to discover the minimum over total feasible kinematical solutions by optimizing the objective value FS in Equation (77) regarding the five variables θ 0 , θ h , φ t , r 0 / r 0 , and t .

4. Results and Discussions

4.1. Comparisons

Michalowski and Drescher [12] investigated the stability for 3D single-step slopes without considering the soil’s nonlinear features and seismic effects, and came up with the stability number N s ( N s = γ H c r / c t ), from which the critical height for 3D single-step slopes can be derived. Where the slope with the critical height is about to collapse, the safety factor is FS = 1 at this time. To demonstrate the applicability of the current approach, the cohesion coefficient λ c deduced from the stability number N s ( N s = 1 / λ c ) is introduced in Equation (77) to calculate the safety factor FS, corresponding to λ a = 0.0 , m = 1 , k h = 0.0 , where the double steps effect, non-linearity effect, and seismicity effect are eliminated. The detailed data are shown in Table 1, and it is evident that the result is pretty near to one, which verifies the effectiveness of the present method in this scenario.
In the research of seismic effects, the fundamental distinction between the modified pseudo-dynamic method and the existing pseudo-dynamic method lies in the various assumptions on seismic acceleration. The horizontal seismic acceleration in the original pseudo-dynamic analysis is written as
a h = [ 1 + z H ( f 1 ) ] k h g sin ω ( t z V s )
where f is the amplification factor of the seismic wave. According to the above equation, after eliminating the time and position variables and setting f = 1.0 , the horizontal seismic acceleration will be simplified to k h g , namely, the seismic acceleration in the pseudo-static analysis. Similarly, the seismic acceleration in the modified pseudo-dynamic analysis represented by Equation (6) will also be simplified to k h g by removing the time and position variables and setting the damping ratio D = 0.0 . Thus, with the same input parameters set, the results of the modified pseudo-dynamic method and the original pseudo-dynamic method should be theoretically consistent, which can demonstrate the modified pseudo-dynamic method’s efficiency in the seismic analysis of 3D slopes. Table 2 shows the outcomes of the two pseudo-dynamic approaches, and they are strikingly similar, which verifies the applicability of the proposed method once more. The fundamental input parameters are as follows: H = 15   m , β 1 = β 2 = 45 ° , λ a = 0.0 , m = 1.6 , c 0 = 25   kPa , σ t = 30   kPa , T = 0.1   s , V s = .

4.2. Parametric Effects Discussions

In accordance with the Equation (77), nonlinear properties of soil medium, geometric parameters of two-step slopes, and seismic parameters in the modified pseudo-dynamic analysis are all factors that affect the safety factor FS for 3D two-step slopes, the effects of which are depicted in Figure 5, Figure 6, Figure 7 and Figure 8, respectively, with the corresponding basic parameters: B / H = 3 , H = 15   m , ξ 1 = 0.4 , β 1 = 45 ° , β 2 = 60 ° , λ a = 0.1 , m = 1.6 , c 0 = 25   kPa , σ t = 30   kPa , T = 0.1   s , D = 0.0 , V s = 150   m / s , k h = 0.2 .
As shown in Figure 5, the nonlinear coefficient m and axial tensile strength σ t are inversely correlated with the safety factor FS; i.e., the larger the nonlinear coefficient and axial tensile strength lead to the smaller FS. Nonetheless, the initial cohesion c0 is positively connected with FS, and the magnitude of the effect is stronger, as the FS varies from 0.4 to 2.3 corresponding to the initial cohesion value ranging from 12.5 to 32.5 kPa.
Then, Figure 6 illustrates the interactions between the nonlinear parameters. It can be observed that the trend of FS remains constant with σ t and c 0 , but when m is small, the slope stability changes dramatically. Moreover, when nonlinear parameters m , σ t and c 0 are fixed, FS decreases as B/H rises. It means that the stability of 3D slopes is better compared with 2D slopes at the same other parameters. This also demonstrates that treating the 3D slope simply as the 2D case in real engineering is overly conservative.
The effects of geometric parameters on slope stability are shown in Figure 7. It is obvious that the depth coefficient ξ 1 and step-width coefficient λ a are positively correlated with FS, while the inclined angles β 1 and β 2 are negatively correlated with FS. It is worth noting that the FS variation range corresponding to β 2 is larger compared with β 1 in the same range of 45~75°, indicating that β 2 has a stronger effect on slope stability. This has some significance for the design of two-step slopes. Furthermore, for the geometric parameters ξ 1 , λ a , β 1 , and β 2 , when the value of B/H is changed from 2 to 4, there is a significant drop in FS. With a rising B/H, the FS falls slowly and gravitates to 2D solutions.
In Figure 8, the influence of pseudo-dynamic parameters varying with B/H on slope stability is described. The variation pattern of FS with B/H is consistent with the foregoing conclusions; although, FS variations with the modified pseudo-dynamic parameters are all minor, not surpassing 5%.
By comparing the results obtained by the original pseudo-dynamic method and the previous relevant research literature, the reliability of the proposed method has been verified, repeatedly. The results obtained from the sensitivity tests of geometric parameters, pseudo-dynamic parameters, and nonlinear parameters of soil have considerable significance for the design and safety evaluation of three-dimensional stepped slopes in current practical projects.
The design mode of a two-order slope can resist the negative impacts of seismic events on slope stability to a certain extent. In addition, the greater the step width, the greater the depth of the step level, and the smaller the inclination angle β 1 , the better the slope stability.

5. Conclusions

The present paper puts forward kinematic analytical solutions for assessing seismic 3D two-step slopes subjected to non-linearity in soils that couples the modified pseudo-dynamic method with the upper bound theorem. The 3D horn-like rotational failure mechanism is employed to solve internal energy dissipation and external force work, in which seismic work is considered in addition to soil weight work and obtained by the layer-wise summation method. Based on the force-increase technique, the safety factor of 3D two-step slopes can be derived more readily, and the results of this derivation can be applied to geotechnical engineering.
To verify the reliability of the approach in this paper, the results of Michalowski and Drescher [12] are extrapolated for the calculation of FS, which are extremely close to one. Additionally, the modified pseudo-dynamic solutions are compared to the original pseudo-dynamic solutions in order to evaluate the modified pseudo-dynamic method’s efficacy in the seismic stability analysis of 3D two-step slopes.
In this paper, the sensitivity tests are performed to investigate the impacts of various factors on slope stability. Among nonlinear parameters, the initial cohesion c 0 has a larger magnitude of influence on FS. When the nonlinear coefficient m is small, FS varies more with the initial cohesion c 0 and axial tensile strength σ t . Among the geometric parameters, β 2 has a stronger effect on slope stability compared with β 1 . This has some significance for the design of two-step slopes. When the value of B/H is changed from 2 to 4, there is a significant drop in FS. With a rising B/H, the FS falls slowly and gravitate to 2D solutions. It is significant to mention that FS varies slightly with the pseudo-dynamic parameters. The above conclusions have a certain meaning for the design, construction, and stability evaluation for 3D two-step slopes subjected the seismic effects.
In the future, the authors will further study the stability of the 3D two-step slope in the following aspects: (1) The support methods of the 3D slope affected by seismic effects, such as prestressed anchor cable, anti-sinking pile reinforcements, etc.; (2) Special conditions of slope medium and complex engineering conditions, for instance, non-uniformity of soils, infiltration of surface water, etc.

Author Contributions

Conceptualization, methodology, software, formal analysis, investigation, resources, data curation, Y.Y.; writing—original draft preparation, writing—review and editing, visualization, D.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Guizhou Province Science and Technology Plan Project under grant no. 2013-122-001.

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.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hovland, H.J. Three-dimensional slope stability analysis method. J. Geotech. Eng. Div. 1977, 103, 971–986. [Google Scholar] [CrossRef]
  2. Su, L.J.; Sun, C.N.; Yu, F.W.; Sarfraz, A. Seismic stability analysis of slopes with pre-existing slip surfaces. J. Mt. Sci. 2018, 15, 1331–1341. [Google Scholar] [CrossRef]
  3. Krahn, J. The 2001 RM Hardy Lecture: The limits of limit equilibrium analyses. Can. Geotech. J. 2003, 40, 643–660. [Google Scholar] [CrossRef]
  4. Luo, W.J.; Yang, X.L. 3D stability of shallow cavity roof with arbitrary profile under influence of pore water pressure. Geomech. Eng. 2018, 16, 569–575. [Google Scholar]
  5. Hou, C.T.; Yang, X.L. 3D stability analysis of tunnel face with influence of unsaturated transient flow. Tunn. Undergr. Space Technol. 2022, 123, 104414. [Google Scholar] [CrossRef]
  6. Griffiths, D.V.; Lane, P.A. Slope stability analysis by finite elements. Geotechnique 1999, 49, 387–403. [Google Scholar] [CrossRef]
  7. 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]
  8. Yang, X.L.; Yin, J.H. Slope stability analysis with nonlinear failure criterion. J. Eng. Mech. 2004, 130, 267–273. [Google Scholar] [CrossRef]
  9. Yang, X.L.; Li, L.; Yin, J.H. Seismic and static stability analysis for rock slopes by a kinematical approach. Geotechnique 2004, 54, 543–549. [Google Scholar] [CrossRef]
  10. 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]
  11. Farzaneh, O.; Askari, F. Three-dimensional analysis of nonhomogeneous slopes. J. Geotech. Geoenvironmental Eng. 2003, 129, 137–145. [Google Scholar] [CrossRef]
  12. Michalowski, R.L.; Drescher, A. Three-dimensional stability of slopes and excavations. Géotechnique 2009, 59, 839–850. [Google Scholar] [CrossRef]
  13. 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]
  14. Steedman, R.S.; Zeng, X. The influence of phase on the calculation of pseudo-static earth pressure on a retaining wall. Geotechnique 1990, 40, 103–112. [Google Scholar] [CrossRef]
  15. Zeng, X.; Steedman, R.S. On the behaviour of quay walls in earthquakes. Geotechnique 1993, 43, 417–431. [Google Scholar] [CrossRef]
  16. Bellezza, I. A new pseudo-dynamic approach for seismic active soil thrust. Geotech. Geol. Eng. 2014, 32, 561–576. [Google Scholar] [CrossRef]
  17. Bellezza, I. Seismic active earth pressure on walls using a new pseudo-dynamic approach. Geotech. Geol. Eng. 2015, 33, 795–812. [Google Scholar] [CrossRef]
  18. Yang, X.L.; Yin, J.H. Slope equivalent Mohr-Coulomb strength parameters for rock masses satisfying the Hoek-Brown criterion. Rock Mech. Rock Eng. 2010, 39, 505–511. [Google Scholar] [CrossRef]
  19. 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]
  20. Choudhury, D.; Katdare, A.D.; Pain, A. New method to compute seismic active earth pressure on retaining wall considering seismic waves. Geotech. Geol. Eng. 2014, 32, 391–402. [Google Scholar] [CrossRef]
  21. Choudhury, D.; Bhattacharyya, S.K.; Pain, A. Seismic stability of retaining wall–soil sliding interaction using modified pseudo-dynamic method. Géotechnique Lett. 2015, 5, 56–61. [Google Scholar]
  22. Gazetas, G.; Garini, E.; Anastasopoulos, I.; Georgarakos, T. Effects of near-fault ground shaking on sliding systems. J. Geotech. Geoenvironmental Eng. 2009, 135, 1906–1921. [Google Scholar] [CrossRef] [Green Version]
  23. Dawson, E.M.; Roth, W.H.; Drescher, A. Slope stability analysis by strength reduction. Geotechnique 1999, 49, 835–840. [Google Scholar] [CrossRef]
  24. Shen, J.; Karakus, M. Three-dimensional numerical analysis for rock slope stability using shear strength reduction method. Can. Geotech. J. 2014, 51, 164–172. [Google Scholar] [CrossRef]
  25. Giasian, H.; Rajabian, A.; Jafari, S. Assessment of Safety Factor of Reinforced Slopes with Respect to Interaction between Soil and Geotextile. J. Transp. Res. 2016, 13, 93–100. [Google Scholar]
  26. Rao, P.; Wu, J.; Mo, Z. 3D Limit analysis of the transient stability of slope during pile driving in nonhomogeneous and anisotropic soil. Adv. Civ. Eng. 2020, 2020, 7560219. [Google Scholar] [CrossRef]
  27. Li, T.Z.; Yang, X.L. An efficient uniform design for Kriging-based response surface method and its application. Comput. Geotech. 2019, 109, 12–22. [Google Scholar] [CrossRef]
Figure 1. Nonlinear strength envelope and corresponding tangential straight line.
Figure 1. Nonlinear strength envelope and corresponding tangential straight line.
Applsci 12 06482 g001
Figure 2. The 3D horn-like rotational failure mechanism for the 3D two-step slope.
Figure 2. The 3D horn-like rotational failure mechanism for the 3D two-step slope.
Applsci 12 06482 g002
Figure 3. Stereoscopic diagram for 3D horn-like rotational mechanism: (a) without plane insert; (b) with plane insert.
Figure 3. Stereoscopic diagram for 3D horn-like rotational mechanism: (a) without plane insert; (b) with plane insert.
Applsci 12 06482 g003
Figure 4. Schematic diagram of the layer-wise summation method: (a) entire graph; (b) the fine-tuned layer.
Figure 4. Schematic diagram of the layer-wise summation method: (a) entire graph; (b) the fine-tuned layer.
Applsci 12 06482 g004
Figure 5. Impacts of nonlinear factors on the stability for 3D two-step slopes: (a) nonlinear coefficient m ; (b) initial cohesion c 0 ; (c) axial tensile strength σ t .
Figure 5. Impacts of nonlinear factors on the stability for 3D two-step slopes: (a) nonlinear coefficient m ; (b) initial cohesion c 0 ; (c) axial tensile strength σ t .
Applsci 12 06482 g005
Figure 6. The interactions between the nonlinear factors on the stability for 3D two-step slopes: (a) nonlinear coefficient m and initial cohesion c 0 ; (b) nonlinear coefficient m and axial tensile strength σ t .
Figure 6. The interactions between the nonlinear factors on the stability for 3D two-step slopes: (a) nonlinear coefficient m and initial cohesion c 0 ; (b) nonlinear coefficient m and axial tensile strength σ t .
Applsci 12 06482 g006
Figure 7. Impacts of geometric factors on the stability for 3D two-step slopes: (a) step-width coefficient λ a ; (b) depth coefficient ξ 1 ; (c) inclined angle β 1 ; (d) inclined angle β 2 .
Figure 7. Impacts of geometric factors on the stability for 3D two-step slopes: (a) step-width coefficient λ a ; (b) depth coefficient ξ 1 ; (c) inclined angle β 1 ; (d) inclined angle β 2 .
Applsci 12 06482 g007
Figure 8. Impacts of seismic factors on the stability for 3D two-step slopes: (a) horizontal seismic coefficient kh; (b) damping ratio D; (c) period T; (d) S-wave propagating velocity Vs.
Figure 8. Impacts of seismic factors on the stability for 3D two-step slopes: (a) horizontal seismic coefficient kh; (b) damping ratio D; (c) period T; (d) S-wave propagating velocity Vs.
Applsci 12 06482 g008
Table 1. Comparisons of both the present outcomes and the solutions of Michalowski and Drescher [12], corresponding to λ a = 0.0 , m = 1 , k h = 0.0 .
Table 1. Comparisons of both the present outcomes and the solutions of Michalowski and Drescher [12], corresponding to λ a = 0.0 , m = 1 , k h = 0.0 .
B/HResults β 1 = β 2  
45 ° 60 ° 75 ° 90 °
1.0Ns by Michalowski and Drescher54.85023.83514.70111.028
FS by the present method0.98930.98230.99930.9849
1.5Ns by Michalowski and Drescher46.84520.77312.9768.935
FS by the present method0.99220.97430.97540.9924
2.0Ns by Michalowski and Drescher42.73219.10312.1098.604
FS by the present method0.99940.99190.97960.9598
3.0Ns by Michalowski and Drescher39.95617.87311.1847.974
FS by the present method1.00191.00100.99550.9925
5.0Ns by Michalowski and Drescher37.99417.06310.6287.266
FS by the present method1.00631.00330.99860.9939
10.0Ns by Michalowski and Drescher36.70316.52710.2656.944
FS by the present method1.00851.00261.00121.0004
Table 2. Comparisons between the modified pseudo-dynamic solutions and the original pseudo-dynamic solutions, corresponding to H = 15   m , β 1 = β 2 = 45 ° , λ a = 0.0 , m = 1.6 , c 0 = 25   kPa , σ t = 30   kPa , T = 0.1   s , V s = .
Table 2. Comparisons between the modified pseudo-dynamic solutions and the original pseudo-dynamic solutions, corresponding to H = 15   m , β 1 = β 2 = 45 ° , λ a = 0.0 , m = 1.6 , c 0 = 25   kPa , σ t = 30   kPa , T = 0.1   s , V s = .
B/HModified Pseudo-Dynamic MethodOriginal Pseudo-Dynamic Method
kh = 0.2kh = 0.5kh = 0.8kh = 0.2kh = 0.5kh = 0.8
1.51.8441.7471.6931.8551.7571.688
3.01.5941.5271.4371.5891.5171.429
5.01.5011.4401.3531.5021.4281.358
10.01.4451.3781.3081.4451.3811.322
2D1.4121.3301.2761.4231.3331.277
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/app12136482

AMA Style

Yang Y, Zhou D. Seismic Stability for 3D Two-Step Slope Governed by Non-Linearity in Soils Using Modified Pseudo-Dynamic Approach. Applied Sciences. 2022; 12(13):6482. https://doi.org/10.3390/app12136482

Chicago/Turabian Style

Yang, Yushan, and De Zhou. 2022. "Seismic Stability for 3D Two-Step Slope Governed by Non-Linearity in Soils Using Modified Pseudo-Dynamic Approach" Applied Sciences 12, no. 13: 6482. https://doi.org/10.3390/app12136482

APA Style

Yang, Y., & Zhou, D. (2022). Seismic Stability for 3D Two-Step Slope Governed by Non-Linearity in Soils Using Modified Pseudo-Dynamic Approach. Applied Sciences, 12(13), 6482. https://doi.org/10.3390/app12136482

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