Abstract

According to the results of the physical model tests, the failure plane of an anaclinal layered rock slope was a linear-type plane at an angle above the plane normal to the discontinuities, and the failure mode of rock strata was bending tension. However, the shear failure occurred near the slope toe, the effects of the cohesion of the discontinuities on the stability of the slope, and the contribution of tangential force to cross-section axial force were neglected in such studies. Moreover, none of the experts had developed a rigorously theoretical method for determining the angle between the failure plane and the plane normal to the discontinuities. This paper was initiated for the purpose of solving the problems described above. With the cantilever beam model and a step-by-step analytical method, an optimization of the analytical method for determining the flexural toppling failure plane based on the limit equilibrium theory was developed and the corresponding formulations were derived. Based on the present computational framework, comparisons with other studies were carried out by taking a slate slope in South Anhui in China and a rock slope facing the Tehran-Chalus Road near the Amir-Kabir Dam Lake in Iran. Furthermore, the sensitivity analyses of the parameters used in the calculation process of the failure angle of the slate slope in South Anhui in China were performed. The results demonstrated that the failure plane and the safety factor of the stability obtained with the presented method were credible, which verified the proposed method. The dip angle of the slope, the dip angle of the rock stratum, and the friction angle of the discontinuities were the controlling factors for the overall failure of the slate slope in South Anhui in China.

1. Introduction

Toppling failure, being a kind of typical instability mode for rock slopes, widely exists in domestic and foreign water conservancy, hydropower, highways, and open-pit mine slope engineering. The toppling instability of slopes has caused great harm to people’s lives, property safety, and engineering construction [14]. However, due to the cognition limitations of the public with regard to anaclinal layered slopes, these slopes had not been able to arouse people’s attention until anaclinal layered toppling landslides frequently occurred. Currently, more and more experts are devoting themselves to the study of anaclinal layered toppling slopes. In [5], the authors divided the toppling failures into three groups, i.e., block toppling, flexural toppling and block-flexural toppling (block-flexural toppling can be treated as a recombined form of the above two types). Numerous significant results had been found for block toppling, and relatively complete theoretical analysis methods were developed [612], while researchers continued to be preoccupied with flexural toppling [1322].

The limit equilibrium method is one of the most common and effective approaches in studies of flexural toppling failure. However, before using the limit equilibrium method for mechanical analyses, the position and shape of the failure plane of a slope must be determined. Thus, many experts put forward various methods to study this problem, and they carried out physical model tests to verify their methods. In [16], the authors first proposed a stability analysis method for slopes and underground openings under various loading conditions against the flexural toppling failure, derived by employing the equations of the limit equilibrium and the boundary conditions. They verified the method by carrying out base friction tests and found that the failure plane was a linear-type plane, emanating from the toe of the slope and perpendicular to the discontinuities, i.e., the angle between the failure plane and the plane normal to the discontinuities (called the failure angle in this study) was equal to 0°. In [13, 23] and [24], the authors improved Aydan and Kawamoto’s theory of flexural toppling failure in open excavations through centrifuge tests and found that the failure angle was 10°. By using the principles of compatibility and the equations of equilibrium along with the governing equations of elastic deformation for the beams, the authors of [25, 26] derived equations for determining the intercolumn forces in rock masses with the potential for flexural toppling failure. Since this model did not allow for slippage between layers, this method might yield very low safety factors for short layers near the slope toe, which might only be true for small deflections. In [15], the authors carried out a series of model tests on a single column and slopes under dynamic loading and found that the failure angle was in the range of 0° to 15°. In [22], the authors developed a new UDEC Trigon approach for simulating the flexural toppling failure. The simulated results showed that the inclination of the total failure surface was around 13° above the plane normal to the discontinuities.

The above theoretical investigations greatly enriched the understanding of the deformation and failure mechanism of such slopes. The significant achievements can be summarized as two viewpoints, i.e., the failure plane of an anaclinal layered rock slope was a linear-type plane at an angle above the plane normal to the discontinuities, and the failure mode of rock strata was bending tension. However, these studies still have some issues needed to be resolved. For instance, the failure of rock strata was not only the result of bending tension as shear failure often occurred near the slope toe, which was supported by field investigations and theoretical analysis. Furthermore, the effects of the cohesion of the discontinuities on the stability of anaclinal layered rock slopes against flexural toppling failure were neglected, which resulted in underestimation of the stability. In addition, the contribution of tangential force to cross-section axial force was neglected, which caused the calculated tensile stress to be overestimated. Last but not least, the failure angle was obtained only from the laboratory data or numerical simulation results, and none of the experts had developed a rigorously theoretical method for determining the failure angle.

In order to solve the above problems, an optimization of the analytical method for determining the flexural toppling failure plane based on the limit equilibrium theory [2729] was developed in this study. The basal failure plane was considered as a plane at which the stress of the slope arrived at the state of limit equilibrium, i.e., the plane, at which the residual sliding force at the toe of the slope was equal to zero, was the basal failure plane. With a cantilever beam model and a step-by-step analytical method, the corresponding formulations were derived in this study. Furthermore, the controlling factors for the overall failure of an anaclinal layered rock slope were obtained after the sensitivity analyses of the parameters used in the calculation process of the slope failure angle were performed.

2. Geological Geometrical Model of an Anaclinal Layered Slope

The geological geometrical model of an anaclinal layered slope is shown in Figure 1, where H is the slope height, β is the dip angle of the slope, η is the dip angle of the rock strata, η0 is the natural slope angle, α is the dip angle of the plane normal to the discontinuities, θ is the dip angle of the failure plane, θr (called failure angle) is the angle between the basal failure plane and the plane normal to the discontinuities, β0 is the difference between β and α, b is the rock thickness, and hi is the contact height between stratum i and stratum i + 1. According to the geometrical conditions, equation (1) can be obtained as follows:

3. Searching Failure Plane

3.1. Searching Principle

According to the results of the physical model tests [13, 15, 16, 2226], the failure plane of the failed strata was a linear-type plane at an angle above the plane normal to the discontinuities. However, none of the experts had developed a rigorously theoretical method for determining the failure angle. In this study, the basal failure plane was considered as a plane at which the stress of the slope arrived at the state of limit equilibrium, i.e., the plane, at which the residual sliding force at the toe of the slope was equal to zero, was the basal failure plane.

Due to the limitation of the computer’s accuracy, it was unable to guarantee the existence of the failure angle that made the residual sliding force at the toe of the slope equal to zero. An attempt was made to control the error of the failure angle to be in an acceptable range by taking a small size of searching step, and thus, the residual sliding force at the toe of the slope was close enough to zero. The searching direction is shown in Figure 2.

Obviously, the above problem was actually an optimization problem with the failure angle θr as a variable parameter, which can be written as follows:

Here, , calculated by the shear failure mode, is the residual sliding force at the toe of the slope, T0, calculated by the tension failure mode, is the residual sliding force at the toe of the slope, β0 is the difference between the dip angle of the slope and the plane normal to the discontinuities, and nn is the searching times. Taking 106 as the value of the searching times, the error of failure angle is less than 10−4, which is in an acceptable range.

Here, Δθr is the size of the searching step. According to the geometrical conditions, the contact height hi between stratum i and stratum i + 1 can be written as follows:

Here, ntp is the number of the first stratum on the top of the slope (numbered from the toe to the top of a slope).

After hi was obtained through equation (5), the mechanical analyses of rock strata can be performed based on the slope geometry and the limit equilibrium theory [2729], which are presented in Sections 3.2 and 3.3.

3.2. Failure Modes

Flexural toppling failure is one of the main types of toppling failures in anaclinal layered rock slopes, and its failure mechanism is due to the bending of layers with or without cross joints. For flexural toppling-type anaclinal layered slopes, although bending tension failure occurred in most of the rock strata, bending tension failure was unlikely to occur in several strata with a small slenderness ratio at the toe of the slopes [22, 30]. Loaded by self-weight and thrust provided by the upper and underlying strata, the strata with a small slenderness ratio were likely subjected to shear failure according to actual field investigations [31, 32] and theoretical analysis [33, 34]. In this study, the failure modes of the rock strata in flexural toppling-type anaclinal layered slopes were shear failure and bending tension failure. The specific failure modes of the rock strata (tension or shear) depended on the types of stress (tension or shear) that reached their critical state earlier.

3.2.1. Tension Failure

If tension failure occurs in stratum i, it should meet the criterion of the maximum tensile stress theory [35], which can be expressed as follows:where σt is the tensile strength of column i, Mi is the bending moment exerted on the center of column i, I is the second moment of inertia, and Ni is the axial force exerted on column i.

As shown in Figure 3, the bending moment exerted on column i can be written as follows:where

Here, is the weight of the part of stratum i above the failure plane, is the equivalent height of stratum i above the failure plane, γ is the unit weight, Ti is the normal force provided by the part of stratum i + 1 above the failure plane, Qi is the tangential force provided by the part of stratum i + 1 above the failure plane, is the normal force provided by the part of stratum i − 1 above the failure plane, Qi−1 is the tangential force provided by the part of stratum i 1 above the failure plane, hi−1 is the contact height of stratum i and stratum i − 1, ϕi is the friction angle of the discontinuities, and χ is the height of the thrust line. As the derivation and verification of χ = 0.6 obtained in [13] were extremely rigorous, the value of χ was equal to 0.6 in this manuscript.

For the calculation of the axial force exerted on column i, the tangential forces (Qi and Qi−1), neglected in previous studies, were reconsidered here, which can be written as follows:

For the calculation of the tangential force (Qi and Qi−1), the cohesion of the discontinuities (ci), neglected in previous studies, was reconsidered here, which can be written as follows:

From equations (6)–(15), the value of can be determined as follows:

3.2.2. Shear Failure

If shear failure occurs in stratum i, it should meet the Mohr–Coulomb criterion, which can be expressed as follows:

Here, τ is the shear stress, σ is the normal stress, ϕ is the friction angle of the rock strata, and c is the cohesion.

Multiplying the bottom length of the stratum on both sides of equation (18), equation (19) can be obtained as follows:

Namely,

Here, S and N represent the tangential force and the normal force of stratum i, respectively.

As shown in Figure 4, the tangential force and the normal force on the bottom of stratum i above the failure plane can be written as follows:

For the calculation of the tangential force (Qi and Qi−1), the cohesion of the discontinuities (ci), neglected in previous studies, was reconsidered in this study, which can be written as follows:

Here, is the normal force provided by the part of stratum i + 1 above the failure plane and is the normal force provided by the part of stratum i − 1 above the failure plane.

From equations (17)–(23), the value of can be determined as follows:

3.2.3. Criteria for Determining Failure Modes

If , it indicated that the thrust causing the shear failure in stratum i was less than that causing bending fracture failure. In this case, we believed that stratum i would be subjected to shear failure rather than bending fracture failure, and vice versa.

3.3. Amount of Failed Strata

During the calculation process, if max , failure would not occur for stratum i, and the interaction force between stratum i and stratum i − 1 would be zero. If max , it indicated that failure occurred for stratum i. Furthermore, , which meant no force provided by the upper stratum was exerted on the last stratum. After the failure modes of the rock strata were determined, the residual sliding force at the toe of the slope, max (P0, T0), can be derived with the step-by-step method.

Then, the failure angle θr was changed constantly and the above process was repeated until the residual sliding force at the toe of the slope was less than zero for the first time. Comparing the absolute value of the residual sliding force of the last step and that of the penultimate step, the smaller value was the desired result, recorded as F. The failure angle corresponding to F was the final failure angle. In this case, the amount of total failed strata, the amount of shear failed strata and the amount of bending fracture strata, recorded as nm, ns, and nt, respectively, can be obtained.

4. Numerical Example

4.1. Slate Slope in South Anhui

A slate slope in South Anhui in China [31] was taken as an engineering example. Table 1 and Figure 5 show the parameters and the geometry profile of the slope [31], respectively.

Using the proposed method developed in this paper, the final results were presented as follows.

Toppling failure occurred in the slate slope in South Anhui in China, which was consistent with field observations and the results reported in [30, 31]. The angle between the failure plane and plane normal to discontinuities, θr, was 7.9766°, and the amount of total failed strata was 38. Shear failure occurred in strata 1–3, while bending tension failure occurred in the remaining failed strata. The safety factor of stability was 0.7573, and the residual sliding force was 26.2239 MN. The failure angle obtained in [30] was 13°, and the amount of failed strata was 37. Shear failure occurred in strata 1–4, and bending tension failure occurred in the remaining failed strata. The safety factor of stability calculated with Aydan and Kawamoto’s method [16] was 0.9343, and the residual sliding force was 4.6781 MN (see Table 2). The safety factor of stability calculated with Majdi and Amini’s method [26] was 0.5430. The result comparisons of the failure planes are shown in Figure 6, and the convergence curves of the safety factor of the stability are shown in Figure 7.

According to the data presented in Section 4.1, it was found that the safety factor of the stability of the slate slope in South Anhui in China obtained with the proposed method was smaller than that calculated with Aydan and Kawamoto’s method [16]. For Aydan and Kawamoto’s method [16], the plane normal to the discontinuities was considered to be the failure plane. However, the overall failure of the slate slope in South Anhui in China occurred before the damage developed to the plane normal to the discontinuities. Therefore, the safety factor of stability of the slate slope in South Anhui in China calculated with Aydan and Kawamoto’s method [16] was overestimated. As the effects of the cohesion of the discontinuities on the stability of the slate slope in South Anhui in China were neglected for the method developed in [30], the failure plane obtained in [30] was above the actual failure plane. Thereupon, the failure plane and the safety factor of the stability of the slate slope in South Anhui in China obtained with the presented method were credible.

4.2. Rock Slope Facing the Tehran-Chalus Road

A rock slope facing the Tehran-Chalus Road near the Amir-Kabir Dam Lake in Iran [36] was taken as another example. Table 3 and Figure 8 show the parameters and the geometry profile of the slope [36], respectively.

Using the proposed method developed in this paper, the final results were presented as follows.

The rock slope facing the Tehran-Chalus Road near the Amir-Kabir Dam Lake in Iran was stable, which was consistent with field observations and the results reported in [36]. The safety factor of stability obtained with the proposed method was 3.1666, and the residual sliding force was −6.1078 MN. The safety factor of the stability calculated with Majdi and Amini’s method [26] was 3.5990. The safety factor of the stability calculated with Aydan and Kawamoto’s method [16] was 3.1082, and the residual sliding force was −7.3278 MN. The safety factor of the stability obtained in [36] was 2.6037, and the residual sliding force was −3.1683 MN. The result comparisons of the safety factor of the stability and the residual sliding force are shown in Table 4, and the convergence curves of the safety factor of the stability are shown in Figure 9.

According to the data presented in Section 4.2, it was found that the safety factor of the stability of the rock slope facing the Tehran-Chalus Road, obtained with the proposed method, was larger than that, obtained in [36]. In [36], the authors assumed that the joints completely cut some of the rock strata and no tensile stress existed in these rock strata, which was relatively rare in the actual slope. Therefore, the safety factor of stability of the rock slope facing the Tehran-Chalus Road obtained in [36] was underestimated. On the one hand, the effects of the cohesion of the discontinuities and the cohesion of rock strata on the stability of a slope were neglected for Aydan and Kawamoto’s method [16], resulting in the underestimation of the safety factor of stability of the slope. On the other hand, the potentially shear failure of stratum 1 at the toe of the slope was neglected for Aydan and Kawamoto’s method [16], resulting in the overestimation of the safety factor of stability of the slope. The safety factor of stability of the rock slope facing the Tehran-Chalus Road obtained with Aydan and Kawamoto’s method was almost equal to that obtained with the presented method due to the combined effects of the two aspects. For the Majdi and Amini’s method [26], the effects of the cohesion of rock strata, the cohesion of the discontinuities, the friction angle of the rock strata, and the friction angle of the discontinuities on the stability of anaclinal layered rock slopes against flexural toppling failure were neglected, leading to a decrease of the accuracy of the safety factor of the stability. Thereupon, the safety factor of the stability of the rock slope facing the Tehran-Chalus Road obtained with the presented method was credible.

5. Parametric Sensibility Analyses for the Failure Angle

The failure angle θr is an important indicator for measuring the degree of instability of anaclinal layered slopes. The position of the failure plane becomes higher as the failure angle θr increases, which indicates that it takes less time to arrive at the state of limit equilibrium, and thus, the slope is more unstable. During the computational process, many parameters may influence the failure angle θr. In this section, the parametric sensibility analyses for the failure angle θr are discussed in detail. All these parameters are classified into two groups, i.e., the geometrical parameters of a slope and the mechanical parameters of rock mass. The geometrical parameters of a slope contain the height of the thrust line χ, the rock thickness b, the slope height H, the dip angle of the slope β, the dip angle of the rock stratum η, and the natural slope angle η0. The mechanical parameters of the rock mass contain the tensile strength σt, the cohesion of rock strata c, the cohesion of the discontinuities ci, the friction angle of the rock strata ϕ, the friction angle of the discontinuities ϕi, and the unit weight γ. When the value of one parameter changed, the others remained the same as in Table 1.

5.1. Sensitivity Analyses of Geometrical Parameters of the Slate Slope in South Anhui

The value of one geometrical parameter of the slate slope in South Anhui was changed constantly, and the corresponding failure angles calculated through the above process are listed in Tables 5 and 6. The failure angles changed with the changes of the geometrical parameters of the slope, and the changing laws are shown in Figure 10.

From the data in Tables 5 and 6 and Figure 10, some conclusions can be drawn as follows:(1)The relatively sensitive parameters for the failure angle θr in the geometrical parameters of the slate slope in South Anhui were the dip angle of the slope β and the dip angle of the rock stratum η, which was consistent with the result obtained in [30]. The effect degree of these two sensitive parameters for the failure angle θr was β>η. For the anaclinal layered rock slopes, their failure was mainly associated with the upper parts of the rock strata above the plane normal to the discontinuities. When the dip angle of the slope or the dip angle of the rock stratum became larger, the possible failure zones of the slope increased, and thus, the slope was more unstable.(2)When β = 45° or η = 55°, the slope was stable. This result further verified Goodman and Bray’s viewpoint [5] that only if inequality (25) was satisfied, the overall failure of an anaclinal layered rock slope would occur.(3)The thickness b is an important indicator for measuring the ability to resist the bending deformation [37]. However, the thickness b was not the most sensitive parameter for the failure angle θr in the geometrical parameters of the slate slope in South Anhui according to the results obtained in this study, i.e., it was not the controlling factor for the overall failure of the slate slope, which was consistent with the opinion of [30].(4)When the geometrical parameters of the slope, except for the rock thickness b, became larger, the failure angle θr increased, which indicated that the slope was more unstable.

5.2. Sensitivity Analyses of Mechanical Parameters of Rock Mass of the Slate Slope in South Anhui

The value of one mechanical parameter of the rock mass of the slate slope in South Anhui was changed constantly, and the corresponding failure angles calculated through the above process are listed in Tables 7 and 8. The failure angles changed with the changes of the mechanical parameter of the rock mass, and the changing laws are shown in Figure 11.

From the data in Tables 7 and 8 and Figure 11, some conclusions can be drawn as follows:(1)The most sensitive parameter for the failure angle θr in the mechanical parameters of the rock mass of the slate slope in South Anhui was the friction angle of the discontinuities, ϕi, which verified Goodman and Bray’s viewpoint [5] that the intense interlayer movement was the precondition for the toppling failure.(2)The tensile strength σt is an important indicator for measuring the ability to resist the bending tension failure [38]. The cohesion c and the friction angle of the rock strata ϕ are the major parameters of shear failure [39]. However, these three parameters were proven not to be the most sensitive parameters for the failure angle θr according to the results obtained in this study, i.e., they were not the controlling factors for the overall failure of the slate slope, which was consistent with the opinion of [40].(3)When the mechanical parameters of the rock mass, except for the unit weight γ, became larger, the failure angle θr became smaller, which indicated that the slope was more stable.

6. Discussion

Based on the above analyses, it was clear that the dip angle of the slope, β, the dip angle of the rock stratum, η, and the friction angles of the discontinuities, ϕi, were the most sensitive parameters for the failure angle θr, which indicated that these three parameters were the controlling factors for the overall failure of the slate slope in South Anhui in China. The effect degree of these three controlling factors for the overall failure of the slate slope was β>η>ϕi.

As we all know, the failure of a slope depends on the stress distribution and the strength of the rock mass. The stress distribution of a slope is mainly affected by the geometrical morphology of the slope, and the strength of the rock mass is mainly controlled by the characteristics of the rock mass discontinuities. For an anaclinal layered slope, the dip angle of the slope, β, and the dip angle of the stratum, η, are the two main factors of the geometrical morphology. The friction angle of the discontinuities, ϕi, is a significant parameter for measuring the ability to resist movement between layers. Hence, these three parameters are the controlling factors for the overall failure of the slate slope in South Anhui in China.

7. Conclusions

As the shear failure occurred near the slope toe, the effects of the cohesion of the discontinuities on the stability of anaclinal layered rock slopes against flexural toppling failure and the contribution of tangential force to cross-section axial force were neglected in previous studies, and none of the experts had developed a rigorously theoretical method for determining the angle between the failure plane and the plane normal to the discontinuities for a slope, an optimization of the analytical method for determining the flexural toppling failure plane based on the limit equilibrium theory was developed and the corresponding formulations were derived in this study. After that, the effectiveness of the proposed method was verified by taking a slate slope in South Anhui in China and a rock slope facing the Tehran-Chalus Road near the Amir-Kabir Dam Lake in Iran, and the sensitive analyses of the parameters used in the calculation process of the slope failure angle were performed as well. Eventually, some conclusions can be drawn as follows.

The proposed method, developed based on the limit equilibrium theory, is feasible to determine the flexural toppling failure plane, which was verified by the case history study.

If the geometrical parameters of a slope, except for the rock thickness, increase, the slope will be more unstable. If the mechanical parameters of the rock mass, except for the unit weight, increase, the slope will be more stable. The dip angle of the slope, the dip angle of the rock stratum, and the friction angle of the discontinuities are the controlling factors for the overall failure of the slate slope in South Anhui in China.

It should be noted that the results obtained in this study only applied to the ultimate failure of anaclinal layered slopes and can be used for preliminary evaluation of the stability.

Data Availability

All data generated or analyzed during this study are included in this article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This study has been financially supported by the Doctoral Scientific Research Fund of the Anyang Institute of Technology (Grant no. BSJ2018009). The authors would like to thank LetPub (http://www.letpub.com) for providing linguistic assistance during the preparation of this manuscript.