Stability Analysis of Rock Slope with Multilayer Weak Interlayer

*e existence of the weak interlayer of the rock slope changes its mechanical characteristics. To ensure the safety of the slope, it is necessary to analyze the overall stability of the slope. Taking the double-layer weak interlayer rock slope beside 318 National Road in Qiyue Mountain, Hubei Province, as an example, a slope model with a weak interlayer was established through GTS software, and the model was imported into FLAC3D for calculation, and the deformation of the slope by the double-layer soft interlayer was studied. *e influence of characteristics and safety factors reveals the controlling effect of the double-layer weak interlayer on the stability of the slope and its failure mode. *e potential sliding surface of the slope is determined to be the lower weak interlayer, and the weak interlayer after the anchor cable reinforcement is carried out. Numerical analysis shows that the reasonable application of anchor cables significantly improves slope stability.*e research results can provide reference significance for slope stability analysis of similar projects.


Introduction
With the rapid construction of infrastructure in Western China, a large number of rock slopes are distributed along highways and railways. Some rock slopes contain weak interlayers which play a controlling role in slope stability. e occurrence of slope disasters causes a series of safety and environmental problems. erefore, it is very important to evaluate slope stability with a weak interlayer [1][2][3].
Currently, many scholars have done some research on this issue. Fan et al. [3,4] studied transfer function and dynamic failure energy identification methods for rock slope with a weak interlayer. Wu et al. [5] studied the reinforcement effect of anchor cable-frame support and determined its role in maintaining slope stability. Chen and Yang [6] studied the dynamic response rules of rock slope with antidip weak interlayer under the seismic wave. Qin [7] analyzed the influence of weak interlayer on slope stability on the engineering side. Feng et al. [8] studied the failure mode of gently inclined bedding slope by the bottom friction method. Some scholars [9][10][11] used the numerical simulation methods to study the failure mode and stability of rock slope and have achieved more meaningful results. Fu et al. [12] applied the method of mechanical theory analysis to construct the deformation characteristics and failure mechanism of the slope with a weak interlayer. Li and Yang [13] analyzed the monitoring data after excavation of the bedding slope with weak rock in the Xijiafen area of Tongzi County, Zunyi City, and concluded that the slope is in a stable state. Wu et al. [14] took a typical soft interlayer slope on the Shanghai-Chengdu Expressway as an example and compared the corresponding softening model and the ideal elastoplastic model.
In summary, the stability of rock slope with weak interlayer needs further research. Taking the slope of 318 National Highway in Qiyueshan, Hubei Province as an example, this paper closely combines the engineering practice and does numerical simulation based on the project, analyzes the controlling effect of multilayer weak interlayer on slope stability, the mechanism of slope deformation, and stress and then puts forward slope reinforcement scheme.
is study can provide theoretical and empirical support for similar projects.

General Situation
e slope of 318 National Highway in Qiyue Mountain has eroded hilly landform with undulating hills. e strata in the slope area are a weathered residual layer, slightly weathered limestone, and limestone. Limestone is hard and has strong weathering resistance. Karst is well developed, and karst fissure water is filled between fissures. Affected by karst water, the local rock mass is weak.
(1) Weathered Residual Layer. is layer's soil property is mainly granular clay which has color of yellowbrown and yellow, the granular composition is mainly silicified shale breccia which has grain size of 1∼2 cm, and grain shape is subedge and subcircle. e content of silicified shale breccia is about 20%∼ 30% (wet soft plastic-plastic).
(2) Striated Limestone. Color is yellowish, cyan-grey, and purple-red, the structure is cryptocrystalline tortoise crack shape, and the layer is mediumthick and thick. Two layers of weak interlayer are intercalated in the slope which has main components such as debris, softening mud, and fault gouge.

Calculation Model
According to the general situation of slope and considering the effect of weak interlayer on slope stability, the slope model with weak interlayer is established by GTS software and imported into FLAC3D for calculation and analysis. e weak interlayer influence on the deformation characteristics and safety coefficient of the slope is analyzed. e calculation gives different physical parameters to the model rock and the weak interlayer, respectively. e parameters are determined according to the large-scale shear test, as shown in Table 1. In the process of calculation, other factors are appropriately simplified, without considering the role of water. e calculation model consists of 5498 nodes and 2683 nodes and used the Mohr-Coulomb material model. In the model, three and five parts of soil are weak interlayer parts, which are given parameters to analyze the overall stability of the slope, respectively. e model is shown in Figure 1. Directional constraints are applied to the bottom, left, and right sides of the slope model. Slope model constraints are shown in Figure 2

Control Analysis of Weak Interlayer on the Upper Slope.
e same parameters of the lower weak interlayer and the surrounding soil are kept, only the upper weak interlayer is retained, and other conditions remain unchanged. Based on these conditions, slope stability is analyzed. erefore, the influence of the upper interlayer on the slope stability is analyzed and determined whether the upper weak interlayer is the control surface of the slope stability. Under these conditions, the horizontal and vertical displacement distributions of the slope are shown in Figures 3 and 4, respectively.
From Figure 3, the maximum horizontal displacement occurs in the upper rock of a weak interlayer, especially on the lower part empty slope of the upper sliding body. It develops along the interlayer toward both ends, which may form a slip zone along the weak interlayer.
From Figure 4, under the influence of the upper weak interlayer, the upper rock has obvious vertical displacement compared with the slope of other parts. e displacement of the negative direction of the Z-axis appears under the influence of self-weight. e maximum vertical displacement occurs at the top of the slope. Under the influence of the weak interlayer, the vertical displacement of the upper rock is larger than that of other parts. Because of the existence of a weak interlayer, the stability of the upper rock becomes worse. e deformation of rock tends to be continuous, and there is no abrupt change in deformation.
e horizontal and vertical stress distributions of the slope are shown in Figures 5 and 6, respectively. e stress distribution of the slope changes with the shape of the slope. From Figure 5, under the natural stress state, the slope is in a local tension state, which occurs in the intersection zone between the weak interlayer and the upper rock, which generally is an unstable region. From Figure 6, the stress distribution of different values is parallel to the slope. Positive value is tension stress, and a negative value is compression stress. Maximum compression stress exists at a certain depth in the slope.
In Figure 7, there is a large shear strain increment along the weak interlayer and surrounding rock, which clearly shows the position of the potential sliding surface of the slope. Figure 8 shows the change of the plastic zone at the weak interlayer. Compared with other parts, there is a more obvious plastic zone along the weak interlayer, most of which is a shear failure. Figure 9 shows that the upper rock displacement of the interlayer is more obvious than that of other parts and also shows the potential displacement slip surface. e principle of the strength reduction method can be summarized as dividing the internal friction angle φ value and the cohesive force c value by the reduction coefficient F si , assuming that the elastic modulus E and Poisson's ratio v remain unchanged, and then input them as new parameters.

Advances in Civil Engineering
In the formula, c is the cohesive force, kPa; φ is the internal friction angle, (°); σ is the total stress, kPa; F si is the strength reduction factor; c i is the cohesive force corresponding to the reduction factor F si , kPa; φ i is the internal friction angle corresponding to the reduction factor F si , (°); and τ i is the shear strength corresponding to the reduction factor F si , kPa.
When the slope is unstable, the strength reduction factor F s is the safety factor. e safety coefficient calculated by the strength reduction method is 1.2, which is higher than the minimum safety coefficient of the slope in the stable state. e slope stability is good, but the influence of a weak interlayer on the slope stability is small. e deformation vector diagram clearly shows the position of the potential sliding surface of the slope. e analysis and simulation of the upper weak interlayer accord with the actual situation of practice.

Control Analysis of Weak Interlayer on the Lower Slope.
e same parameters of the upper weak interlayer and the surrounding soil are kept, only the lower weak interlayer is retained, and other conditions remain unchanged. Based on these conditions, the influence of the lower interlayer on the slope stability is analyzed and determined whether the lower weak interlayer is the control surface of the slope stability. Under these conditions, the horizontal and vertical displacement distributions of the slope are shown in Figures 10  and 11, respectively.
From Figure 10, the maximum horizontal displacement occurs at the slope foot position. e upper rock of the weak interlayer has a larger displacement than other parts; the closer to the slope foot, the greater the displacement; it develops along the interlayer toward both ends, which may form a slip zone along the weak interlayer.
From Figure 11, the displacement of the negative direction of the Z-axis appears under the influence of selfweight. e maximum vertical displacement occurs at the top of the slope. Under the influence of the weak interlayer, the vertical displacement of the upper rock is larger than that of other parts. Because of the existence of a weak interlayer, the stability of the upper rock becomes worse. e deformation diagram can clearly show the position of the slope potential sliding surface.
Slope horizontal and vertical stress distributions are shown in Figures 12 and 13, respectively. From Figure 12, we can see the horizontal stress distribution of the slope. ere is larger tension stress along the weak interlayer upper part, while there is a larger compression stress concentration at the slope foot position. Different colors in the figure indicate the change of stress. A positive value is tension stress, and a negative value is compression stress. e maximum compression stress in this area exists at the slope foot position. From Figure 13, the same stress distribution is parallel to the slope. e different colors in the figure indicate the change of stress.
e maximum vertical compression stress in this region exists at a certain depth in the slope. Figure 14 shows that there is a large shear strain increment along the weak interlayer and surrounding rock. It also has a large shear strain increment at the slope foot position and a large shear increment along the weak interlayer lower part. e potential slope sliding surface is clearly shown. As we can see from Figure 15, compared with

Block State
None shear-n shear-p shear-n shear-p tension-p shear-n tension-n shear-p tension-p shear-p tension-p  Advances in Civil Engineering slope other parts, the plastic zone along the weak interlayer is more obvious, and most of which are a shear failure. It shows that the existence of the weak interlayer makes the surrounding rock become deformed obviously. Figure 16 shows that the interlayer displacement of the upper rock is more obvious than that of other parts and also shows the potential displacement slip surface. e safety coefficient calculated by the strength reduction method is 1.05, which is close to the minimum safety coefficient in a stable state and is on the verge of stability. e deformation vector diagram clearly shows the position of the potential sliding surface of the slope. e analysis and simulation of the lower weak interlayer accord with the actual situation of practice. Because the Contour of X-Displacement Interval = 1.0e-003 -4.9816e-004 to 0.0000e+000 0.0000e+000 to 1.0000e-003 1.0000e-003 to 2.0000e-003 2.0000e-003 to 3.0000e-003 3.0000e-003 to 4.0000e-003 4.0000e-003 to 5.0000e-003 5.0000e-003 to 6.0000e-003 6.0000e-003 to 7.0000e-003 7.0000e-003 to 8.0000e-003 8.0000e-002 to 8.8999e-003   Contour of SXX Interval = 1.0e+005 -5.0000e+005 to -4.0000e+005 -5.4610e+005 to -5.0000e+005 -4.0000e+005 to -3.0000e+005 -3.0000e+005 to -2.0000e+005 -2.0000e+005 to -1.0000e+005 -1.0000e+005 to 0.0000e+000 0.0000e+000 to 5.8990e+004   slope is on the verge of stability, it is necessary to reinforce the slope.

Simulation Analysis of Entire
Stability. Different parameters are assigned to the upper and lower weak interlayer to do the model calculation. e corresponding parameters are given to the upper and lower weak interlayer, respectively, and the influence of slope stability control surface and potential sliding surface on the slope is studied. e horizontal and vertical displacement distribution of the slope is shown in Figures 17 and 18, respectively. From Figure 17, under the influence of the upper and lower weak interlayers, along the lower part of the two weak interlayers, the displacement of rock is relatively large. e nearer to the slope foot position, the greater the displacement. Developing along the interlayer toward both ends, a slip zone along the weak interlayer may be formed.
From Figure 18, the maximum vertical displacement occurs at the top of the slope, mainly in the upper rock of the weak interlayer. Due to the low strength and relatively high compression stress of the weak interlayer, the large vertical displacement of the upper part of the weak interlayer occurs. e upper rock of slope is affected by the weak interlayer, especially affected by two weak interlayers; the vertical displacement is larger than that of other parts. Because of the existence of a weak interlayer, the stability of the upper rock becomes worse. e deformation diagram can clearly show the position of the potential sliding surface of the slope. e horizontal and vertical stress distributions of the slope are shown in Figures 19 and 20, respectively. From Figure 19, there is large tensile stress along the upper part of the weak interlayer. Affected by two weak interlayers, the area is larger than that of one weak interlayer, while there is a large concentration of compression stress at the slope foot position. e maximum compression stress in this region exists at the slope foot position. From Figure 20, the same stress distribution is parallel to the slope, and the different colors showed the change of stress. e maximum vertical compression stress in this region exists at a certain depth in the slope. Figure 21 shows that there are large shear strain increments along with two weak interlayers and surrounding rock. e shear strain increment along the upper weak interlayer mainly concentrates on the lower part of the interlayer. e shear strain increment along the lower weak interlayer almost runs through the whole weak interlayer. At the same time, there is a large shear strain increment at the slope foot position. e position of the potential sliding surface of the slope is shown. ere are obvious plastic zones along with the two weak interlayers in Figure 22, most of which are shear failures. From Figure 23, the displacement of the two interlayers is obvious. e safety coefficient calculated by the strength reduction method is 1.05, which is close to the minimum safety factor of the slope in the stable state. It is concluded that the slope is on the verge of stability. e deformation vector diagram clearly shows the position of the potential sliding surface of the slope. e analysis and simulation of the lower weak interlayer accord with the actual situation of practice. Because of the importance of slope stability, measurement should be taken to improve slope stability.

Numerical Simulation of Prestressed Anchor
Cable. e anchor cable used at the project site and its layout method is taken as an example. FLAC3D software is used to simulate the reinforcement measures of prestressed anchor cables to   (Figure 24). e reinforcement effect of the anchor cable is evaluated by calculation. en, the corresponding reinforcement measurement is put forward.

Analysis of Reinforcement Effect of Prestressed Anchor
Cable. From Figure 25, the displacement from weak interlayer to the slope is reduced after the anchor cable is applied, and there is only a small maximum displacement zone. From Figure 26, the vertical displacement of the weak interlayer decreases significantly when the anchor cable is applied. From Figure 27, a plastic zone exists in the weak interlayer and slope area before and after anchor cable reinforcement, and the shear strain area decreases, and the maximum shear strain increment decreases.
Contour of X-Displacement Interval = 5.0e-004 0.0000e+000 to 5.0000e-004 -2.5829e-004 to 0.0000e+000 5.0000e-004 to 1.0000e-003 1.0000e-003 to 1.5000e-003 1.5000e-003 to 2.0000e-003 2.0000e-003 to 2.5000e-003 2.5000e-003 to 3.0000e-003 3.0000e-003 to 3.4816e-003 Figure 25: Horizontal displacement diagram after application of the anchor cable. e safety coefficient calculated by the strength reduction method after the anchor cable is applied is 2.2, which can effectively guarantee the safety requirements. e safety coefficient of the weak interlayer slope without anchor cable is 1.05. e reinforcement of the anchor cable significantly improves the stability of the slope.

Conclusion
e numerical model is established based on the rock slope with a multilayer weak interlayer. e influence of upper and lower weak interlayer on slope stability is considered comprehensively. Conclusions are as follows: (1) e weak interlayer plays a controlling role in slope stability. When the slope is damaged, the position of the weak interlayer has been sliding or has a large displacement, and the shear yield stress usually occurs at the position of the weak interlayer. (2) e stability of slope with multilayer weak interlayer is controlled by the state of the main weak interlayer. e lower weak interlayer is the main controlling factor, and the slope is on the edge of stability and lacks sufficient safety reservation.
(3) Reinforcing rock slope with weak interlayer by prestressed anchor cable can effectively improve the integrity of the slope and reduce the occurrence of geological disasters such as landslides.

Data Availability
e data used to support the findings of this study are included within the article.