阶梯状滑动断续节理顺层边坡稳定性分析

李得建; 贾文韬; 程肖; 赵炼恒; 张迎宾; 余鹏程

doi:10.11779/CJGE202211019

阶梯状滑动断续节理顺层边坡稳定性分析 English Version

1.
西南交通大学土木工程学院, 四川成都 610031
2.
四川大学灾后重建与管理学院深地科学与工程教育部重点实验室, 四川成都 610065
3.
中南大学土木工程学院重载铁路工程结构教育部重点实验室, 湖南长沙 410075

基金项目:

四川省科技厅项目 2021YJ0032

四川省科技厅项目 2021YJ0390

同济大学岩土及地下工程教育部重点实验室开放基金项目 KLE-TJGE-B2004

中央高校基本科研业务费-科技创新项目 2682021CX007

专职博士后研发基金项目 2022SCU12126

详细信息

作者简介:
李得建(1990—)，男，讲师，博士，主要从事岩土构筑物极限状态分析及抗减震消能加固设计等方面的教学和科研工作。E-mail: lidejian@swjtu.edu.cn

通讯作者:
程肖，E-mail: chengxiao@scu.edu.cn

中图分类号: TU43
计量
- 文章访问数: 158
- HTML全文浏览量: 28
- PDF下载量: 51
出版历程
- 收稿日期: 2021-11-15
- 网络出版日期: 2022-12-08
- 刊出日期: 2022-10-31

Stability of stepped sliding of bedding rock slopes with discontinuous joints

1.
School of Civil Engineering, Southwest Jiaotong University, Chengdu 610031, China
2.
Institute for Disaster Management and Reconstruction, MOE Key Laboratory of Deep Earth Science and Engineering, Sichuan University, Chengdu 610065, China
3.
School of Civil Engineering, Key Laboratory of Heavy-Haul Railway Engineering Structure, Ministry of Education, Central South University, Changsha 410075, China

摘要

摘要: 阶梯状滑动顺层岩质滑坡在工程中较为常见，其岩桥和节理特性对阶梯状滑动面的形成至关重要。考虑岩桥和节理面之间强度参数弱化特性，构建阶梯状滑动顺层边坡稳定性分析模型。研究表明：与平面型滑动对比验证了分析模型的准确性。裂隙连通率k_Y、基本摩擦角φ_b、粗糙度系数JRC和弱化系数K_c对安全系数F_s影响显著。岩桥倾角β₃越大，F_s减小越显著，边坡越不稳定。台阶法中开挖次数和坡角角度对边坡稳定性影响显著，通过算例分析验证了该方法与极限平衡解差异率在1%以内。某些情况下当开挖使节理面出露时，F_s显著降低（如β=53°时降低了22.8%），易诱发滑坡的发生。同时，随着拉裂缝倾角β₁和充水高度h_w的增大，F_s显著降低，出流缝被堵塞时水力效应更为显著，此时边坡最不稳定。
- 节理岩质边坡 /
- 阶梯状滑动 /
- 弱化系数 /
- 台阶型坡面 /
- 水力作用
Abstract: The bedding rock landslides with stepped sliding surface are common in engineering, and the characteristics of rock bridges and joints are very important to the formation of stepped sliding surface. The weakening characteristics of the strength parameters of rock bridges and joints are considered. The models for stability of stepped sliding of bedding slopes are established. The results show that the accuracy of the model is verified by comparing with the planar failure model. In addition, the fracture connectivity rate k_Y, basic friction angle φ_b, roughness coefficient JRC and weakening coefficient K_c have significant influences on F_s. In this mode, the greater β₃ is, the more significant the reduction of F_s is, and the more unstable the slope is. In addition, the excavation times and slope angle have significant influences on the slope stability. The difference rate between the proposed method and the limit equilibrium solution is less than 1%. In some cases, when the joint surfaces are exposed by excavation, F_s is significantly reduced (the reduction of 22.8% when β=53°), which easily leads to the occurrence of landslides. Meanwhile, with an increase in angle β₁ of tensile cracks and height h_w of water, F_s significantly decreases. The water pressure is more significant when the outlets are blocked, and the rock slope is most unstable.
- joint rock slope /
- stepped sliding /
- weakening coefficient /
- step-type slope surface /
- water pressure

HTML全文

0. 引言

实际工程中岩质边坡大多是由岩桥和不连续节理面组成的复杂地质体^[1-2]，且大多数岩质边坡节理面相对于工程尺度呈显著的非贯通特征。由于岩桥和不连续节理面大都埋入边坡内部，节理面的发育具有显著的随机性，故对其稳定性分析较为复杂和困难。在断续节理顺层岩质边坡内部存在共面和雁型两种节理面^[3-4]。破坏模式区别在于岩桥贯通模式：若是节理间岩桥贯通则为平面型破坏模式；若是相邻节理间岩桥贯通则为阶梯状（或多阶梯）破坏模式。在强震^[5-6]、降雨^[7-8]、开挖^[9]、长期蠕变时空演化及其他作用^[10-11]下发生的阶梯状滑坡表明：当岩桥发生破坏后，通过与相邻节理面搭接或多层节理面搭接形成阶梯状和多阶梯状破坏模式（图 1）。

图 1 三荔高速顺层岩质阶梯状滑面失稳案例^[8]

Figure 1. Typical case of stepped sliding of jointed rock landslides in Sandu-Libo Expressway^[8]

下载: 全尺寸图片幻灯片

基于平面型破坏模式，Jennings^[12]提出了断续节理强度参数的加权平均计算公式：tanφ = ktanφ_j + (1-k)tanφ_r，c=kc_j+(1-k)c_r（k为裂隙连通率）。该式考虑了岩桥和节理面抗剪强度参数的差异性，却未考虑其相关性（本文研究对象为次生节理，该类节理可视为由完整岩石演化而来）。虽然邹宗兴^[13]探究了顺层岩质边坡滑带抗剪强度参数弱化系数的问题，但是，从规范中相关规定和实际情况可知^{[2, 14]}，内摩擦角正切值与黏聚力两者的弱化程度并不一致。目前研究大都忽略了岩桥和节理面强度参数的相关性和差异性，或对两者仍采用相同的破坏准则，其合理性受到削弱。节理面抗剪性能并不总能满足M-C破坏准则，尤其当 ${\sigma _{\text{n}}}$ 较小时计算出的剪切强度明显偏大^[15]。Barton等^[16-18]通过研究节理面直剪特性，并基于试验结果提出了可用于评估不规则、无充填节理面峰值抗剪强度的非线性JRC-JCS经验公式（即B-B破坏准则），可较好反映节理面的抗剪性能。

综上所述，首先，岩桥和节理面分别服从M-C和B-B破坏准则，岩桥和节理面强度参数采用不同弱化系数K_φ和K_c。然后，开展考虑后缘倾斜拉裂缝存在的阶梯状岩桥剪切破坏模式下边坡稳定性研究，开展对比分析和参数敏感性分析，获取各关键参数对断续节理顺层岩质边坡稳定性的影响规律。最后，基于台阶型边坡外形特征，开展台阶法开挖次数和坡角角度对节理顺层岩质边坡安全系数影响规律的探究。同时考虑水力作用的影响，探究出流缝被堵塞和未堵塞对边坡稳定性的影响，分析拉裂缝倾角和拉裂缝充水高度对边坡稳定性的影响规律。

1. 断续节理顺层边坡阶梯状失稳模式

1.1 B-B破坏准则与M-C破坏准则的参数转换

岩桥服从M-C破坏准则，节理面服从B-B破坏准则。考虑岩桥和节理面抗剪强度参数的差异性，利用不同的弱化系数K_c和K_φ获取岩桥强度参数c_r，φ_r。根据相关学者研究^[19-20]，获得两准则参数具体转换为（图 2）

图 2 B-B与M-C破坏准则相关参数转换示意图

Figure 2. Conversion parameters of B-B and M-C failure criteria

下载: 全尺寸图片幻灯片

$\varphi = \arctan \left( {\frac{{\partial \tau }}{{\partial {\sigma _{\text{n}}}}}} \right)\text{，}$

(1)

$\frac{{\partial \tau }}{{\partial {\sigma _{\text{n}}}}} = tan\left( {{\text{JRC}}\lg \frac{{{\text{JCS}}}}{{{\sigma _{\text{n}}}}} + {\varphi _{\text{b}}}} \right) -$

$\frac{\text{πJRC}}{180\mathrm{ln}10}\left[ta{n}^{2}\left(\text{JRC}\mathrm{lg}\frac{\text{JCS}}{{\sigma }_{\text{n}}}+{\phi }_{\text{b}}\right)+1\right]\text{，}$

(2)

$c = \tau - {\sigma _{\text{n}}}\tan \varphi 。$

(3)

通过式（1）~（3）可获取节理面黏聚力c_j和内摩擦角φ_j，结合K_c，K_φ可获得岩桥黏聚力c_r和内摩擦角φ_r。

1.2 岩桥与节理面强度参数弱化特性

断续节理顺层岩质边坡与贯通节理岩质边坡破坏机理有些不同^[21-22]，因而岩桥和节理面强度参数存在显著的弱化特性，并存在以下两方面的假定：

（1）根据节理接触状态，节理面所受到的法向应力等于总法向应力乘以传压系数。本文假定节理面承担的法向应力与岩桥部分承担的法向应力一致。

（2）实际边坡中岩桥和节理面搭接处是应力集中部位，存在明显弱化现象，且向岩桥内部弱化趋势逐渐减弱。本文假定忽略岩桥局部弱化情况，取完整岩桥相关强度参数。

基于以上两个假定，通过引入弱化系数将弱化特性体现在强度参数上面，可有效体现节理面对断续节理顺层岩质边坡稳定性的影响。

1.3 断续节理顺层岩质边坡阶梯状失稳模式构建

基于上述相关假定，构建阶梯状滑动极限状态失稳模式如图 3所示，为求解方便建立直角坐标系，A点坐标为(0，0)，X轴方向指向向右，Y轴方向指向向上。注意：根据黄达等的研究^{[2, 23]}，当岩桥倾角小于90°时主要发生剪切破坏；岩桥倾角大于90°时主要发生张拉破坏。本文构建的分析模型中岩桥倾角小于90°（中 ${\beta _{\text{3}}}$ ），故分析模型适用于岩桥剪切破坏。

图 3 顺层岩质边坡岩桥剪切贯通阶梯状失稳模式图

Figure 3. Shear cut through mode of rock bridges of stepped sliding of jointed rock slope

下载: 全尺寸图片幻灯片

中，边坡坡角和坡顶角分别为 $\beta$ ， $\alpha$ ，OB水平距离为L，边坡高度为H+Ltan $\alpha$ 。BB₁高度为H₁+Ltan $\alpha$ ，倾角为 ${\beta _{\text{1}}}$ ；B₁B₂，B₃B₄为节理面长度，高度分别为H₂和H₄，倾角分别为 ${\beta _{\text{2}}}$ ， ${\beta _{\text{4}}}$ （岩桥间节理面大都近平行，故 ${\beta _{\text{4}}}$ = ${\beta _{\text{2}}}$ ）。B₂B₃为岩桥剪切破坏长度，高度为H₃，倾角为 ${\beta _{\text{3}}}$ 。BC₂，B₂B₅，B₁B₆与B₄C₂垂直，BC₂₂与B₄C₂₂垂直。坐标参数：B(H/tanβ+L，H+Ltanα)，B₁(H/tanβ+L-(H₁+Ltanα)/tan ${\beta _{\text{1}}}$ , H-H₁)，B₂(H/tanβ+L- (H₁+Ltanα)/tan ${\beta _{\text{1}}}$ -H₂/tan ${\beta _{\text{2}}}$ ，H-H₁-H₂)。

阶梯状滑动裂隙连通率表达式为

${k_Y} = \frac{{\frac{{{H_4} + {H_2}}}{{\sin {\beta _2}}} + \frac{{({H_1} + L\tan \alpha ) \cdot \cos ({\beta _1} - {\beta _2})}}{{\sin {\beta _1}}}}}{{\frac{{{H_4} + {H_2}}}{{\sin {\beta _2}}} + \frac{{{H_3} \cdot {\text{cos(}}{\beta _3} - {\beta _2})}}{{\sin {\beta _3}}} + \frac{{({H_1} + L\tan \alpha ) \cdot \cos ({\beta _1} - {\beta _2})}}{{\sin {\beta _1}}}}} ={k}_{Y1}+{k}_{Y2}\text{，}$

(4)

式中，k_Y为裂隙连通率。k_Y1可提供抗滑力（节理面），k_Y2不可提供抗滑力（后缘倾斜拉裂缝）。

H₄可由H₃，y₄表示：

${H_4}{\text{ = }}H - {H_1} - {H_2} - {y_4} - {H_3} 。$

(5)

通过式（6），y₄可由H₃表示：

${y_4} = {y_{44}} - {H_3}\frac{{\sin ({\beta _3} - {\beta _2}) \cdot \sin \beta }}{{\sin (\beta - {\beta _2}) \cdot \sin {\beta _3}}} 。$

(6)

将式（5），（6）代入式（4），可获得H₃表达式：

${H}_{3}=\frac{(1-{k}_{Y})\frac{({H}_{1}+L\mathrm{tan}\alpha )\cdot \mathrm{cos}({\beta }_{1}-{\beta }_{2})}{\mathrm{sin}{\beta }_{1}}+(1-{k}_{Y})\frac{H-{H}_{1}-{y}_{44}}{\mathrm{sin}{\beta }_{2}}}{({k}_{Y}-1)\frac{\mathrm{sin}({\beta }_{3}-{\beta }_{2})\cdot \mathrm{sin}\beta }{\mathrm{sin}(\beta -{\beta }_{2})\cdot \mathrm{sin}{\beta }_{3}\cdot \mathrm{sin}{\beta }_{2}}+{k}_{Y}\frac{cos({\beta }_{3}-{\beta }_{2})}{\mathrm{sin}{\beta }_{3}}+\frac{(1-{k}_{Y})}{\mathrm{sin}{\beta }_{2}}}。$

(7)

进而获取点B₄的坐标x₄和y₄：

$\left.\begin{array}{l}{x}_{4}=\frac{H-{H}_{1}-\mathrm{tan}{\beta }_{2}\left(\frac{H}{\mathrm{tan}\beta }+L-\frac{{H}_{1}+L\text{tan}\alpha }{\mathrm{tan}{\beta }_{1}}\right)}{\mathrm{tan}\beta -\mathrm{tan}{\beta }_{2}}-\\ {H}_{3}\frac{\mathrm{sin}({\beta }_{3}-{\beta }_{2})\cdot \mathrm{cos}\beta }{\mathrm{sin}(\beta -{\beta }_{2})\cdot \mathrm{sin}{\beta }_{3}}\text{，}\\ {y}_{4}={x}_{4}\cdot \mathrm{tan}\beta 。\end{array}\right\}$

(8)

将式（8）中y₄代入式（5），可获得H₄表达式。将H₃，H₄及其他相关参数代入阶梯状滑动面岩桥剪切破坏稳定性分析中，可获取其安全系数表达式。

2. 断续节理顺层岩质边坡阶梯状失稳模式的构建

2.1 断续节理阶梯状滑动岩桥剪切破坏稳定性分析

边坡沿缓倾裂隙面滑移导致岩桥发生剪切破坏贯通。因此，可沿缓倾裂隙方向（图 3中B₄C₂方向）计算抗滑力F_resist和下滑力F_induce从而得到该模式下安全系数表达式^[12]。

（1）边坡体沿B₄C₂方向抗滑力的计算

节理面和岩桥承担的法向应力相同，投影在B₄C₂上坡体法向应力 ${\sigma _{{\text{n}}Y}}$ 为

${\sigma _{{\text{n}}Y}} = \frac{{W\cos {\beta _4}}}{{{B_4}{B_6}}} \text{，}$

(9)

式中，W为滑坡体自重，W=γS_Y滑，S_Y滑为滑坡体面积。S_Y滑的计算过程如下：

a）四边形OBC₂₂B₄面积S_Y1的求解

${S_{Y1}} = \frac{{L \cdot L\tan \alpha }}{2} + L \cdot ({H_1} + {H_2} + {H_3} + {H_4}) +$

$\text{ }\frac{{({H}_{1}+{H}_{2}+{H}_{3}+{H}_{4})}^{2}}{2\mathrm{tan}\beta }\text{ }。$

(10)

b）多边形BC₂₂B₄B₃B₂B₁面积S_Y2的求解

$\begin{array}{c}{S}_{Y2}=\frac{{(L\mathrm{tan}\alpha +{H}_{1})}^{2}}{2\mathrm{tan}{\beta }_{1}}+\left(\frac{L\mathrm{tan}\alpha +{H}_{1}}{\mathrm{tan}{\beta }_{1}}+\frac{{H}_{2}}{2\mathrm{tan}{\beta }_{2}}\right)\cdot {H}_{2}+\\ \left(\frac{L\mathrm{tan}\alpha +{H}_{1}}{\mathrm{tan}{\beta }_{1}}+\frac{{H}_{2}}{\mathrm{tan}{\beta }_{2}}+\frac{{H}_{3}}{2\mathrm{tan}{\beta }_{3}}\right)\cdot {H}_{3}+\\ \left(\frac{L\mathrm{tan}\alpha +{H}_{1}}{\mathrm{tan}{\beta }_{1}}+\frac{{H}_{2}}{\mathrm{tan}{\beta }_{2}}+\frac{{H}_{3}}{\mathrm{tan}{\beta }_{3}}+\frac{{H}_{4}}{2\mathrm{tan}{\beta }_{4}}\right)\cdot {H}_{4}\text{，}\end{array}$

(11)

${S}_{Y滑}={S}_{Y1}-{S}_{Y2} 。$

(12)

用H₅≥0对滑动面在坡表出露进行约束：

${H_5} = H - {H_1} - {H_2} - {H_3} - {H_4} (H_{5}≥0) 。$

(13)

c）B₄B₆长度

${B_4}{B_6} = \frac{{{H_4}}}{{\sin {\beta _4}}} + \frac{{{H_3} \cdot {\text{cos(}}{\beta _3} - {\beta _4})}}{{\sin {\beta _3}}} + \frac{{{H_2}}}{{\sin {\beta _2}}} 。$

(14)

阶梯状滑动面为岩桥剪切破坏与相邻节理面连接贯通而成，岩桥倾角大于节理面倾角，故β₃ > β₂=β₄。

由于岩桥和节理面对边坡稳定性的贡献不同，故B₄B₆方向抗滑力为

${F_{{\text{resist}}}} = ({c_r} + {\sigma _{{\text{n}}Y}} \cdot \tan {\varphi _r}) \cdot {B_3}{B_5} + ({B_5}{B_6} + {B_4}{B_3}) \cdot$

$\left\{{\sigma }_{\text{n}Y}\mathrm{tan}\left[{\phi }_{\text{b}}+\text{JRC}\mathrm{lg}\left(\frac{\text{JCS}}{{\sigma }_{\text{n}Y}}\right)\right]\right\}\text{ }。$

(15)

（2）边坡体沿B₄B₆方向下滑力的计算

${F_{{\text{induce}}}} = W \cdot \sin {\beta _4} 。$

(16)

（3）阶梯状岩桥剪切破坏模式安全系数的计算

$\begin{gathered} {F_{\text{s}}} = {\text{\{ }}({c_r} + {\sigma _{{\text{n}}Y}} \cdot \tan {\varphi _{\text{r}}}) \cdot {B_4}{C_2} \cdot (1 - {k_Y}) + \hfill \\ \left\{ {{\sigma _{{\text{n}}Y}}\tan \left[ {{\varphi _{\text{b}}} + {\text{JRC}}\lg \left( {\frac{{{\text{JCS}}}}{{{\sigma _{{\text{n}}Y}}}}} \right)} \right]} \right\} \cdot \hfill \\ \end{gathered}$

${B}_{4}{C}_{2}\cdot {k}_{Y1}\text{}\}/W\cdot \mathrm{sin}{\beta }_{4}\text{ }。$

(17)

将相关参数代入式（17），可求得阶梯状滑动岩桥剪切破坏时边坡安全系数F_s。

2.2 各关键参数对阶梯状滑动岩桥剪切破坏影响分析

（1）对比分析

开展平面型与阶梯状失稳模式下安全系数的对比分析对比结果如图 4所示。参数取值：① k_Y=0.5~1.0， ${\beta _3}$ =35°~80°，K_φ=0.6，K_c=0.2；②k_Y=0.5~1.0， ${\beta _3}$ =35°~80°，K_φ=0.8，K_c=0.4。

图 4 弱化系数K_φ和K_c下安全系数F_S随岩桥倾角β₃和裂隙连通率k_Y的变化曲线

Figure 4. Safety factor F_S versus rock bridge angle β₃ and fracture connectivity rate k_Y under weakening coefficients K_φ and K_c

下载: 全尺寸图片幻灯片

其他参数取值：H=100 m，L=25 m， $\alpha$ =10°， $\beta$ =75°， $\gamma$ =25 kN/m³，H₁=2 m，H₂=3 m， ${\beta _1}$ =75°， ${\beta _2}$ =35°，φ_b=20°，JRC=6，JCS=40 MPa。

由可知，当 ${\beta _{\text{3}}}$ = ${\beta _{\text{2}}}$ =35°时为平面型失稳模式。随 ${\beta _{\text{3}}}$ 的增大，平面型失稳模式转化为阶梯状失稳模式，F_s呈减小趋势，且 ${\beta _{\text{3}}}$ 越大，减小得越快，这与文献[4]中数值模拟的结果规律相一致。由图 4（a）可知，k_Y越小岩桥对边坡稳定性的贡献越大，故F_s越高。此时随着 ${\beta _{\text{3}}}$ 的变化，对F_s的影响越来越显著。由图 4对比可知，随K_φ和K_c的增大，岩桥抗剪强度参数相应减小，故F_s相对降低，这是由于以节理面参数为基础，K_φ和K_c越大，岩桥强度参数越小，故F_s越低。

（2）参数敏感性分析

探究裂隙连通率k_Y和基本摩擦角φ_b、粗糙度系数JRC和壁面有效抗压强度JCS、弱化系数K_φ和K_c、岩桥倾角 ${\beta _{\text{3}}}$ 和节理面倾角 ${\beta _{\text{2}}}$ 对阶梯状滑动边坡稳定性的影响，具体如图 5所示。参数取值如表 1所示，其他参数取值：图 5（a）~中H=50 m， $\alpha$ =10°， $\gamma$ =25 kN/m³，H₁=2 m，H₂=3 m。中H=100 m， $\alpha$ =10°， $\gamma$ =25 kN/m³，H₁=1 m，H₂=3 m。节理面参数取值时，应保证φ_b+JRClg(JCS/ ${\sigma _{\text{n}}}$ )最大值不大于70°。

图 5 阶梯状失稳模式各关键参数对边坡稳定性的影响规律

Figure 5. Influences of key parameters on stability of stepped failure mode of slope

下载: 全尺寸图片幻灯片

表 1 参数分析图 5参数取值

Table 1. Values of parameters in Fig. 5

参数	k_Y	φ_b/(°)	K_φ	K_c	JRC	JCS /MPa	β₃/(°)	β₂ /(°)	β/(°)	L/m	β₁/(°)
图 5(a)	0.2~1.0(0.1)	20~36(4)	0.8	0.4	6	50	50	40	60	16	60
图 5(b)	0.8	28	0.650.6~0.95 (0.05)	0.05~0.40 (0.05)	6	30	60	45	60	15	60
图 5(c)	0.65	20	0.6	0.3	0~18(3)	20~100(20)	50	40	60	22	60
图 5(d)	0.6	20	0.8	0.4	10	40	30~75(5)	25~50(5)	58	10	60

下载: 导出CSV

| 显示表格

由图 5可知，随k_Y，K_φ和K_c的增大，F_S显著减小。φ_b增大同样提高了岩桥和节理面抗剪强度参数，Fs呈非线性增大趋势。F_S随JRC和JCS的增大均呈非线性增大趋势，只是后者增加得越来越慢。F_S随β₃的增大而非线性减小，且减小趋势非常显著，即岩桥倾角越大，减小越快。这是由于岩桥倾角越大，构成的阶梯状滑动面更容易下滑，边坡越不稳定。故所呈现的效果与滑动面的倾角越大边坡越不稳定的情况类似，因此，F_S越小。F_S随 ${\beta _{\text{2}}}$ 的增大而非线性减小，且减小趋势较为显著。

3. 阶梯状滑动台阶型坡面岩质边坡稳定性分析

3.1 台阶型坡面边坡阶梯状失稳模式构建

典型顺层岩质边坡表现为节理面在坡面出露或者随着开挖进行在开挖面出露。目前针对边坡的开挖处置大都基于台阶法开展，通常情况下坡率设置为1︰0.75（53.13°）或1︰1（45°），台阶之间设置平台（通常设置为2 m，可根据现场实际情况相应调整），达到对边坡卸载和使之更加稳定的目的。台阶型坡面顺层岩质边坡阶梯状岩桥剪切贯通模式如图 6所示。

图 6 台阶型坡面顺层岩质边坡阶梯状岩桥剪切贯通模式图

Figure 6. Shear cut through mode of rock bridges of stepped sliding of jointed rock slope based on step-type slope surface

下载: 全尺寸图片幻灯片

基于图 3，6分别设置若干级边坡和台阶，自上而下为A_iD_i (i=1, 2, …，n-1)和A_nA，倾角为α_i (i=1, 2, …, n-1)，高度为h₁+L_A1tanα₁、h_i (i=2, 3, …, n)，平台宽度为d_i (i=1, 2, …, n-1)。考虑边坡开挖顺序均是自上而下进行，故着重探究开挖次数和各级边坡坡率对阶梯状滑动岩质边坡稳定性的影响分析。

（1）台阶型坡面引起边坡几何尺寸的变化

基于2.1节，通过边坡裂隙连通率、几何尺寸和岩土体参数获得阶梯状滑动面位置和安全系数表达式，具体依照2.1节相关内容求解。

台阶法使阶梯状滑动面上部岩体重力荷载减小，若按照2.1节计算公式，滑动面上的正应力、岩桥和节理面参数均随之变化，故为获得开挖次数和坡角角度的影响，设定每组情况下岩桥和节理面抗剪强度参数均一致，而正应力因为上覆荷载的变化而随之变化，进而影响抗滑力计算。开挖效应主要引起边坡体几何尺寸和潜在滑动面上覆岩土体重度的变化，不考虑开挖时滑坡体面积为S_Y滑（式（12）），开挖掉的滑坡体面积为S_Y挖，具体相关几何尺寸计算如下：

台阶型坡面点A₁至坡顶点O的水平距离L_A1为

${L_{A1}} = \frac{{\sum\limits_{i = 1}^{n - 1} {\frac{{{h_i}}}{{\tan {\alpha _i}}}} + \frac{{{h_n}}}{{\tan \beta }} + \sum\limits_{i = 1}^{n - 1} {{d_i}} - \frac{H}{{\tan \beta }}}}{{1 - \frac{{\tan \alpha }}{{\tan {\alpha _1}}}}} 。$

(18)

每级开挖体面积计算为

${S}_{Y挖i}=\left({L}_{A1}\cdot \left(1-\frac{\mathrm{tan}\alpha }{\mathrm{tan}{\alpha }_{1}}\right)-\left({\displaystyle \sum _{1}^{i-1}(\frac{{h}_{j}}{\mathrm{tan}{\alpha }_{j}}+{d}_{j}-}\frac{{h}_{j}}{\mathrm{tan}\beta }\right)\right){h}_{i}+\\ \frac{{h}_{i}{}^{2}}{2\mathrm{tan}\beta }-\frac{{h}_{i}{}^{2}}{2\mathrm{tan}{\alpha }_{i}}\text{ (}i\text{=1,2}\text{，}\cdots \text{，}n-\text{1)}\text{，}$

(19)

故开挖体总面积为

${S}_{Y挖}={\displaystyle \sum _{i}^{n-1}{S}_{Y挖i}} 。$

(20)

此时阶梯状滑动台阶型坡面滑坡体区域面积为

${S}_{Y台}={S}_{Y滑}-{S}_{Y挖}\text{ }。$

(21)

有些情况下滑坡体潜在滑动面从台阶型边坡坡面或者平台间出露，即此时未完全包含挖坡体，故在计算开挖后滑坡体面积时需将多减去的部分重新加回。

（2）边坡体沿B₄C₂方向抗滑力的计算

a）若B₃B₄未因开挖而出露，可参照式（9）进行法向应力计算。故投影在B₄C₂上坡体法向应力 ${\sigma }_{\text{n}Y台}$ 为

${\sigma }_{\text{n}Y台}=\frac{{W}_{台}\cdot \mathrm{cos}{\beta }_{4}}{{B}_{4}{B}_{6}} \text{，}$

(22)

式中，W_台为滑坡体自重，W_台= $\gamma {S}_{Y台}$ _， ${S}_{Y台}$ 为滑坡体面积，具体计算见式（21）。

由于开挖时潜在滑动面并未出露，故B₄B₆的长度和抗滑力计算如式（14），（15）不变。

b）若B₃B₄因开挖而出露，此时节理面长度变短，即采用新的B₃B_4台代入式（14），（15）求解因节理面出露而造成法向应力、B₄B₆的长度和抗滑力改变后的数值。

（3）边坡体沿B₄B₆方向下滑力为

${F}_{\text{induce}}={W}_{台}\cdot \mathrm{sin}{\beta }_{4} 。$

(23)

（4）台阶型坡面阶梯状失稳模式安全系数为

${F}_{\text{s}}=\{({c}_{r}+{\sigma }_{\text{n}Y台}\cdot \mathrm{tan}{\phi }_{\text{r}})\cdot {B}_{3}{B}_{5}+({B}_{5}{B}_{6}+{B}_{4}{B}_{3})\cdot$

$\left\{{\sigma }_{\text{n}Y台}\mathrm{tan}\left[{\phi }_{\text{b}}+\text{JRC}\mathrm{lg}\left(\frac{\text{JCS}}{{\sigma }_{\text{n}Y}}\right)\right]\right\}\}{W}_{台}\cdot \mathrm{sin}{\beta }_{4}\text{，}$

(24)

式中，若因开挖节理面出露，采用新的 ${\sigma }_{\text{n}Y台}$ 和B₃B_4台代入相应的位置进行替代。

将相关参数代入式（24）中，可以求得台阶型坡面阶梯状失稳模式岩桥剪切破坏时安全系数F_s。

3.2 算例分析

（1）对比分析

锦屏一级水电站右岸泄洪洞引渠内侧岩质边坡共形成六级台阶，属于典型多级台阶状坡面岩质高边坡。由于文献[24]中采用的各节理面裂隙连通率为100%，故此时发生平面型滑动破坏模式。相关参数如下：后缘拉裂缝处内摩擦角为16.70°，黏聚力为20 kPa，抗拉强度为200 kPa。滑动面g_8sz-1错动带处内摩擦角为30.96°，黏聚力为170 kPa。滑动面g_yj7错动带处内摩擦角为28.81°，黏聚力为150 kPa^[24]。基于本节方法开展与文献[24，25]中计算结果的对比分析（表 2）。

表 2 本文计算结果与文献[24，25]计算结果的对比分析

Table 2. Comparison between calculated results and those obtained in References [24, 25]

滑动面组合	不考虑拉裂缝的贡献^[24]
g_8sz-1错动带	文献[24]极限分析上限解	文献[24]极限平衡解	本文
g_8sz-1错动带	1.058657	1.058652	1.047927
g_yj7错动带	文献[24]极限分析上限解	文献[24]极限平衡解	本文
g_yj7错动带	1.018695	1.018693	1.007019
滑动面组合	考虑拉裂缝的贡献^[25]
g_8sz-1错动带	文献[25]极限平衡解	本文
g_8sz-1错动带	1.289	1.216197
g_yj7错动带	文献[25]极限平衡解	本文
g_yj7错动带	1.256	1.214398

下载: 导出CSV

| 显示表格

由表 2可知，与文献[24]计算结果十分接近，差异率在1%以内，与文献[25]计算结果也较为接近，差异率在5%以内，从而验证了本文方法的准确性。结果表明：伴随着后缘拉裂缝的开展，F_S降低达17.08%，表明后缘拉裂缝开展对边坡稳定性的显著影响。

（2）算例分析

台阶型开挖次数和坡角角度对节理岩质边坡稳定性影响显著，故探究它们对台阶型岩质边坡稳定性的影响规律。参数取值为：k_Y=0.75， ${\beta _1}$ =60°， ${\beta _2}$ =30°， ${\beta _3}$ =50°， $\beta$ =45°~53°（坡率1︰1~1︰0.75）， $\alpha$ =0°，H=50 m，H₁=4 m，H₂=3 m，L=15 m，γ=25 kN/m³，K_φ=0.6，K_c=0.2，φ_b=20°，JRC=3，JCS=20 MPa。本节以五级边坡为例，开挖次数和坡角角度对岩质边坡阶梯状失稳模式安全系数的影响规律具体如图 7所示。图 7（a）中虚线表示未考虑节理面出露对抗滑力贡献的丧失，图 7（b），7（c）中红色虚线表示出露的节理面长度。

图 7 开挖次数和坡角角度对阶梯状滑动台阶型坡面边坡安全系数和破坏模式的影响规律

Figure 7. Influences of excavation times and slope angle on safety factor and failure mode of stepped sliding of slope based on step-type slope surface

下载: 全尺寸图片幻灯片

由图 7可知，开挖坡率较大时（倾角较小），节理面不易随开挖而出露。此时，随着开挖次数的进行，达到对边坡安全卸载的目的，故F_S相应提高（图 7（a））。而随开挖坡率的减小，节理面在边坡下部出露。由于保持裂隙连通率不变，使得阶梯状滑动面剪出口位置相应提高，故其节理面出露长度相应减小（图 7（b），7（c））。由于节理面出露部分对抗滑力贡献的丧失，使F_S显著降低，开挖坡率越小，降低程度越显著，降低程度分别从 $\beta$ =47°的7.2%增加至 $\beta$ =53°的22.8%，说明节理面对边坡抗滑力的影响。

4. 水力作用下阶梯状滑动岩质边坡稳定性分析

4.1 水力作用下阶梯状滑动岩质边坡失稳模式构建

由于边坡后缘拉裂缝的存在，为水力作用提供了过水通道，因而降雨或地下水会对顺层岩质边坡的稳定性产生更为不利影响。基于图 3，考虑出流缝被堵塞和出流缝未堵塞两种工况^[26]，其中，出流缝未堵塞时采用地下水位中点处水压力最大的水力分布情况^[26-27]，构建了水力作用下阶梯状滑动岩质边坡失稳模式如图 8所示。

图 8 水力作用下阶梯状岩质边坡失稳模式图

Figure 8. Failure mode of stepped sliding of jointed rock slope under water pressure

下载: 全尺寸图片幻灯片

基于图 3相关参数，设地下水的重度为γ_w，拉裂缝充水高度为h_w；出流缝被堵塞时节理面、岩桥和拉裂缝所受的水压力分别为U₁，T₁和V₁；出流缝未堵塞时各部分所受的水压力分别为U₂，T₂和V₂。

（1）出流缝被堵塞时边坡安全系数F_S1的计算

根据图 8（a），出流缝被堵塞时节理面、岩桥和拉裂缝所受的水压力分别为

$\left.\begin{array}{l}{U}_{1}=\frac{{\gamma }_{\text{w}}{H}_{2}(2{h}_{\text{w}}+{H}_{2})}{2\mathrm{sin}{\beta }_{2}}+\frac{{\gamma }_{\text{w}}{H}_{4}(2{h}_{\text{w}}+2{H}_{2}+2{H}_{3}+{H}_{4})}{2\mathrm{sin}{\beta }_{4}}\text{，}\\ {T}_{1}=\frac{{\gamma }_{\text{w}}{H}_{3}(2{h}_{\text{w}}+2{H}_{2}+{H}_{3})}{2\mathrm{sin}{\beta }_{3}}\text{，}\\ {V}_{1}=\frac{{\gamma }_{\text{w}}{h}_{\text{w}}^{2}}{2\mathrm{sin}{\beta }_{1}}。\end{array}\right\}$

(25)

基于式（25），坡体法向应力 ${\sigma _{{\text{n}}1}}$ 为

${\sigma _{{\text{n}}1}} = \frac{{W\cos {\beta _4} - {U_1} - {T_1}\cos ({\beta _3} - {\beta _4}) - {V_1}\cos ({\beta _1} - {\beta _4})}}{{{B_4}{B_6}}} 。$

(26)

坡体抗滑力F_resist1为

${F_{{\text{resist1}}}} = ({c_r} + {\sigma _{{\text{n}}1}} \cdot \tan {\varphi _r}) \cdot {B_3}{B_5} + ({B_5}{B_6} + {B_4}{B_3}) \cdot {\text{ }}$

$\left\{{\sigma }_{\text{n1}}\mathrm{tan}\left[{\phi }_{\text{b}}+\text{JRC}\mathrm{lg}\left(\frac{\text{JCS}}{{\sigma }_{\text{n1}}}\right)\right]\right\}\text{ }。$

(27)

坡体下滑力F_induce1为

${F}_{\text{induce1}}=W\mathrm{sin}{\beta }_{4}+{T}_{1}\mathrm{sin}({\beta }_{3}-{\beta }_{4})+{V}_{1}\mathrm{sin}({\beta }_{3}-{\beta }_{4})。$

(28)

出流缝被堵塞时边坡安全系数F_S1为

$\begin{gathered} {F_{{\text{S1}}}} = {\text{\{ }}({c_r} + {\sigma _{{\text{n1}}}} \cdot \tan {\varphi _{\text{r}}}) \cdot {B_4}{C_2} \cdot (1 - {k_Y}) + \hfill \\ \left\{ {{\sigma _{{\text{n}}1}}\tan \left[ {{\phi _{\text{b}}} + {\text{JRC}}\lg \left( {\frac{{{\text{JCS}}}}{{{\sigma _{{\text{n1}}}}}}} \right)} \right]} \right\} \cdot {B_4}{C_2} \cdot {k_{Y1}}{\text{\} /}} \hfill \\ \end{gathered}$

$\text{(}W\cdot \mathrm{sin}{\beta }_{4}+{T}_{1}\mathrm{sin}({\beta }_{3}-{\beta }_{4})+{V}_{1}\mathrm{sin}({\beta }_{1}-{\beta }_{4})\text{)}。$

(29)

（2）出流缝未堵塞时边坡安全系数F_S2的计算

根据图 8（b），出流缝未堵塞时节理面、岩桥和拉裂缝所受的水压力分别为

$\left.\begin{array}{l}{U}_{2}=\frac{{\gamma }_{\text{w}}{H}_{2}(2{h}_{\text{w}}+{H}_{2})}{2\mathrm{sin}{\beta }_{2}}+\frac{{\gamma }_{\text{w}}{H}_{4}^{2}}{2\mathrm{sin}{\beta }_{4}}\text{，}\\ {T}_{2}={\gamma }_{\text{w}}[({h}_{\text{w}}+{H}_{2}+{H}_{3}-{H}_{4})({h}_{\text{w}}+{H}_{2}+{H}_{3}-3{H}_{4})+\\ (-{h}_{\text{w}}-{H}_{2}+{H}_{3}+{H}_{4})(3{h}_{\text{w}}+3{H}_{2}+{H}_{3}+{H}_{4})]/8\mathrm{sin}{\beta }_{3}\text{，}\\ {V}_{2}=\frac{{\gamma }_{\text{w}}{h}_{\text{w}}^{2}}{2\mathrm{sin}{\beta }_{1}}。\end{array}\right\}$

(30)

基于式（30），坡体法向应力 ${\sigma _{{\text{n2}}}}$ 为

${\sigma _{{\text{n2}}}} = \frac{{W\cos {\beta _4} - {U_2} - {T_2}\cos ({\beta _3} - {\beta _4}) - {V_2}\cos ({\beta _1} - {\beta _4})}}{{{B_4}{B_6}}} 。$

(31)

坡体抗滑力F_resist2为

${F_{{\text{resist2}}}} = ({c_{\text{r}}} + {\sigma _{{\text{n}}2}} \cdot \tan {\varphi _{\text{r}}}) \cdot {B_3}{B_5} + ({B_5}{B_6} + {B_4}{B_3}) \cdot$

$\left\{{\sigma }_{\text{n}2}\mathrm{tan}\left[{\phi }_{\text{b}}+\text{JRC}\mathrm{lg}\left(\frac{\text{JCS}}{{\sigma }_{\text{n}2}}\right)\right]\right\}\text{ }。$

(32)

坡体下滑力F_induce2为

${F}_{\text{induce2}}=W\cdot \mathrm{sin}{\beta }_{4}+{T}_{2}\mathrm{sin}({\beta }_{3}-{\beta }_{4})+{V}_{2}\mathrm{sin}({\beta }_{1}-{\beta }_{4})。$

(33)

出流缝被堵塞时边坡安全系数F_S2为

$\begin{gathered} {F_{{\text{S2}}}} = {\text{\{ }}({c_r} + {\sigma _{{\text{n2}}}} \cdot \tan {\varphi _r}) \cdot {B_4}{C_2} \cdot (1 - {k_Y}) + \hfill \\ \left\{ {{\sigma _{{\text{n2}}}}\tan \left[ {{\varphi _{\text{b}}} + {\text{JRC}}{{\lg }_{10}}\left( {\frac{{{\text{JCS}}}}{{{\sigma _{{\text{n2}}}}}}} \right)} \right]} \right\} \cdot {B_4}{C_2} \cdot {k_{Y1}}{\text{\} }}/ \hfill \\ \end{gathered}$

$\text{(}W\cdot \mathrm{sin}{\beta }_{4}+{T}_{2}\mathrm{sin}\text{(}{\beta }_{3}-{\beta }_{4}\text{)}+{V}_{2}\mathrm{sin}({\beta }_{1}-{\beta }_{4})\text{)}。$

(34)

将相关参数代入式（29），（34），可求得出流缝被堵塞和未堵塞时边坡的安全系数F_S1，F_S2。

4.2 拉裂缝相关参数和拉裂缝充水高度对边坡稳定性的影响

裂缝倾角 ${\beta _1}$ 和拉裂缝充水高度h_w对边坡稳定性影响显著。基于此，开展出流缝被堵塞与未堵塞时的对比分析与 ${\beta _1}$ 和h_w对边坡稳定性的参数敏感性分析，如所示。参数取值： ${\beta _1}$ =60°~90°，h_w=0~2.0 m。其他参数取值：H=100 m，L=25 m， $\alpha$ =10°， $\beta$ =75°， $\gamma$ =25 kN/m³，H₁=2 m，H₂=3 m，β₂=35°，φ_b=20°，JRC=6，JCS=40 MPa， ${\gamma _{\text{w}}}$ =10 kN/m³。

图 9 F_s随拉裂缝倾角β₁和拉裂缝充水高度h_w的变化曲线

Figure 9. F_S versus angle β₁ and depth h_w of water in tensile cracks

下载: 全尺寸图片幻灯片

由图 9可知，当出流缝被堵塞或未堵塞时，F_S均随着h_w的增加而降低，且均随着 ${\beta _1}$ 的增加而降低。由与（b）对比可知，当 ${\beta _1}$ =60°时，出流缝被堵塞时F_S1随着h_w的增加降低了22.0%，出流缝未堵塞时F_S2降低了3.1%，出流缝被堵塞时F_S对h_w的变化较敏感。在不考虑水力作用的情况下，当 ${\beta _1}$ =60°时，F_S=1.41；当出流缝被堵塞且h_w=2.0 m时，与不考虑水力作用时相比，F_S下降了62.8%；当出流缝未堵塞且h_w=2.0 m时，与不考虑水力作用时相比，F_S下降了20.3%。因此，水力作用对带裂缝边坡的稳定性有重要的影响。

5. 结论

本文开展了阶梯状岩桥剪切贯通失稳模式下各关键参数对断续节理顺层岩质边坡稳定性影响规律的研究，获得4点结论。

（1）假定岩桥服从M-C破坏准则，节理面服从B-B破坏准则，考虑岩桥和节理面强度参数的弱化系数K_φ，K_c，依据节理面强度参数和弱化系数K_φ，K_c获得了岩桥内摩擦角φ_r和黏聚力c_r。

（2）基于断续节理顺层岩质边坡阶梯状岩桥剪切失稳模式，通过对比分析验证了该模式的准确性。研究表明：裂隙连通率k_G(k_Y)、基本摩擦角φ_b、粗糙度系数JRC、弱化系数K_c和岩桥倾角 ${\beta _3}$ 对F_S影响显著。尤其β₃越大，F_S减小越显著，边坡越不稳定。

（3）通过对比分析验证了台阶型坡面岩质边坡阶梯状失稳模式稳定性计算结果的准确性。研究表明：开挖次数和坡率对断续节理顺层岩质边坡稳定性影响显著。开挖范围不涉及潜在滑动面时边坡更为稳定，而开挖一旦使节理面出露便使得边坡稳定性大幅下降，最高下降达22.8%，因此当开挖范围涉及到潜在滑动面时，此时使得边坡稳定性显著下降。故在实际工程中，不能盲目放坡，应先对目标边坡地质情况进行细致调研，才能真正达到放坡促使边坡更为稳定的目的。

（4）通过对比分析以及参数分析探究了水力作用对边坡稳定性的影响规律。研究表明：水力作用对边坡稳定性的不利影响十分显著。其中，出流缝被堵塞时边坡F_S下降非常显著，h_w=2.0 m时F_S下降达62.8%，出流缝未堵塞时安全系数下降较为显著；拉裂缝倾角和拉裂缝充水高度对边坡稳定性有一定的不利影响。在实际工程中对于带裂缝边坡，不能忽视水力作用的影响，在实际工程中，对易滑边坡及时做出排水措施。

图 1 三荔高速顺层岩质阶梯状滑面失稳案例^[8]

Figure 1. Typical case of stepped sliding of jointed rock landslides in Sandu-Libo Expressway^[8]

下载: 全尺寸图片幻灯片

图 2 B-B与M-C破坏准则相关参数转换示意图

Figure 2. Conversion parameters of B-B and M-C failure criteria

下载: 全尺寸图片幻灯片

图 3 顺层岩质边坡岩桥剪切贯通阶梯状失稳模式图

Figure 3. Shear cut through mode of rock bridges of stepped sliding of jointed rock slope

下载: 全尺寸图片幻灯片

图 4 弱化系数K_φ和K_c下安全系数F_S随岩桥倾角β₃和裂隙连通率k_Y的变化曲线

Figure 4. Safety factor F_S versus rock bridge angle β₃ and fracture connectivity rate k_Y under weakening coefficients K_φ and K_c

下载: 全尺寸图片幻灯片

图 5 阶梯状失稳模式各关键参数对边坡稳定性的影响规律

Figure 5. Influences of key parameters on stability of stepped failure mode of slope

下载: 全尺寸图片幻灯片

图 6 台阶型坡面顺层岩质边坡阶梯状岩桥剪切贯通模式图

Figure 6. Shear cut through mode of rock bridges of stepped sliding of jointed rock slope based on step-type slope surface

下载: 全尺寸图片幻灯片

图 7 开挖次数和坡角角度对阶梯状滑动台阶型坡面边坡安全系数和破坏模式的影响规律

Figure 7. Influences of excavation times and slope angle on safety factor and failure mode of stepped sliding of slope based on step-type slope surface

下载: 全尺寸图片幻灯片

图 8 水力作用下阶梯状岩质边坡失稳模式图

Figure 8. Failure mode of stepped sliding of jointed rock slope under water pressure

下载: 全尺寸图片幻灯片

图 9 F_s随拉裂缝倾角β₁和拉裂缝充水高度h_w的变化曲线

Figure 9. F_S versus angle β₁ and depth h_w of water in tensile cracks

下载: 全尺寸图片幻灯片

表 1 参数分析图 5参数取值

Table 1 Values of parameters in Fig. 5

参数	k_Y	φ_b/(°)	K_φ	K_c	JRC	JCS /MPa	β₃/(°)	β₂ /(°)	β/(°)	L/m	β₁/(°)
图 5(a)	0.2~1.0(0.1)	20~36(4)	0.8	0.4	6	50	50	40	60	16	60
图 5(b)	0.8	28	0.650.6~0.95 (0.05)	0.05~0.40 (0.05)	6	30	60	45	60	15	60
图 5(c)	0.65	20	0.6	0.3	0~18(3)	20~100(20)	50	40	60	22	60
图 5(d)	0.6	20	0.8	0.4	10	40	30~75(5)	25~50(5)	58	10	60

下载: 导出CSV

表 2 本文计算结果与文献[24，25]计算结果的对比分析

Table 2 Comparison between calculated results and those obtained in References [24, 25]

滑动面组合	不考虑拉裂缝的贡献^[24]
g_8sz-1错动带	文献[24]极限分析上限解	文献[24]极限平衡解	本文
g_8sz-1错动带	1.058657	1.058652	1.047927
g_yj7错动带	文献[24]极限分析上限解	文献[24]极限平衡解	本文
g_yj7错动带	1.018695	1.018693	1.007019
滑动面组合	考虑拉裂缝的贡献^[25]
g_8sz-1错动带	文献[25]极限平衡解	本文
g_8sz-1错动带	1.289	1.216197
g_yj7错动带	文献[25]极限平衡解	本文
g_yj7错动带	1.256	1.214398

下载: 导出CSV

参考文献(27)

[1]	EBERHARDT E, STEAD D, COGGAN J S. Numerical analysis of initiation and progressive failure in natural rock slopes—the 1991 Randa rockslide[J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(1): 69–87. doi: 10.1016/S1365-1609(03)00076-5
[2]	DA HUANG, CEN D F, MA G W, et al. Step-path failure of rock slopes with intermittent joints[J]. Landslides, 2015, 12(5): 911–926. doi: 10.1007/s10346-014-0517-6
[3]	朱雷, 黄润秋, 严明, 等. 基于裂纹扩展模式的岩质斜坡阶梯状滑移破裂机制研究[J]. 岩土工程学报, 2017, 39(7): 1216–1224. http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract16953.shtml ZHU Lei, HUANG Run-qiu, YAN Ming, et al. Step-path failure mechanism of rock slopes based on crack coalescence modes in rock mass[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(7): 1216–1224. (in Chinese) http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract16953.shtml
[4]	郭牡丹. 基于岩体结构面特征的三维网络模拟研究[D]. 沈阳: 东北大学, 2014. GUO Mu-dan. Study on 3-Dimensional Network Simulation Based on the Character of Structural Planes in Rock Mass[D]. Shenyang: Northeastern University, 2014. (in Chinese)
[5]	黄润秋, 许强. 中国典型灾难性滑坡[M]. 北京: . 科学出版社, 2008. HUANG Run-qiu, XU Qiang. Catastrophic landslides in China[M]. Beijing: . China Science Press, 2008. (in Chinese)
[6]	GARCÍA M, PASTÉN C, SEPÚLVEDA S A, et al. Dynamic numerical investigation of a stepped-planar rockslide in the Central Andes, Chile[J]. Engineering Geology, 2018, 237: 64–75. doi: 10.1016/j.enggeo.2018.02.001
[7]	TANNANT D D, GIORDAN D, MORGENROTH J. Characterization and analysis of a translational rockslide on a stepped-planar slip surface[J]. Engineering Geology, 2017, 220: 144–151. doi: 10.1016/j.enggeo.2017.02.004
[8]	ZHAO L H, LI D J, TAN H H, et al. Characteristics of failure area and failure mechanism of a bedding rockslide in Libo County, Guizhou, China[J]. Landslides, 2019, 16(7): 1367–1374. doi: 10.1007/s10346-019-01188-6
[9]	杨绪波. 大型岩质边坡开挖的岩体结构效应研究-以云南澜沧江小湾水电站^#4山梁边坡为例[D]. 成都: 成都理工大学,, 2005. YANG Xu-bo. Study on the rock mass structure effect of large rock slope during excavation-taking slope NO. 4 in Xiaowan hydropower project on Lancang river of Yunnan as an example[D]. Chengdu: Chengdu University of Technology, 2005. (in Chinese)
[10]	BRIDEAU M A, YAN M, STEAD D. The role of tectonic damage and brittle rock fracture in the development of large rock slope failures[J]. Geomorphology, 2009, 103(1): 30–49. doi: 10.1016/j.geomorph.2008.04.010
[11]	黄润秋, 陈国庆, 唐鹏. 基于动态演化特征的锁固段型岩质滑坡前兆信息研究[J]. 岩石力学与工程学报, 2017, 36(3): 521–533. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201703001.htm HUANG Run-qiu, CHEN Guo-qing, TANG Peng. Precursor information of locking segment landslides based on transient characteristics[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(3): 521–533. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201703001.htm
[12]	JENNINGS J E. A mathematical theory for the calculation of the stability of slopes in open cast mine[C]// Proceeding of the Symposium on the Planning Open Pit Mines. Johannesburg, 1970.
[13]	邹宗兴. 顺层岩质滑坡演化动力学研究[D]. 武汉: 中国地质大学, 2014. ZOU Zong-xing. Research on the Evolution Dynamics of the Consequent Bedding Rockslides[D]. Wuhan: China University of Geosciences, 2014. (in Chinese)
[14]	重庆市设计院. 工程地质勘察规范DBJ50—043—2005[S]. 2005. Chongqing Design Institute. Code for Engineering Geological Investigation DBJ50—043—2005[S]. 2005. (in Chinese)
[15]	陈祖煜, 汪小刚, 杨健, 等. 岩质边坡稳定分析—原理、方法、程序[M]. 北京: 中国水利水电出版社, 2005. CHEN Zu-yu, WANG Xiao-gang, YANG Jian, et al. Rock Slope Stability Analysis-Theory Methods and Programs[M]. Beijing: China WaterPower Press, 2005. (in Chinese)
[16]	BARTON N. Review of a new shear-strength criterion for rock joints[J]. Engineering Geology, 1973, 7(4): 287–332. doi: 10.1016/0013-7952(73)90013-6
[17]	BARTON N. The shear strength of rock and rock joints[J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1976, 13(9): 255–279. https://www.sciencedirect.com/science/article/pii/0148906276900036
[18]	BANDIS S, LUMSDEN A C, BARTON N R. Experimental studies of scale effects on the shear behaviour of rock joints[J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1981, 18(1): 1–21. https://www.sciencedirect.com/science/article/pii/014890628190262X
[19]	HOEK E. Rock Engineering[M]. Netherlands: A A Balkema, 1995: 70–72.
[20]	PRASSETYO S H, GUTIERREZ M, BARTON N. Nonlinear shear behavior of rock joints using a linearized implementation of the Barton-Bandis model[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2017, 9(4): 671–682. doi: 10.1016/j.jrmge.2017.01.006
[21]	唐志成, 夏才初, 刘远明. 岩桥渐进弱化的Jennings抗剪强度准则[J]. 岩土工程学报, 2012, 34(11): 2093–2099. http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract14891.shtml TANG Zhi-cheng, XIA Cai-chu, LIU Yuan-ming. Modified Jennings shear strength criterion based on mechanical weakening model of rock bridges[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(11): 2093–2099. (in Chinese) http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract14891.shtml
[22]	邓华锋, 齐豫, 李建林, 等. 水–岩作用下断续节理砂岩力学特性劣化机理[J]. 岩土工程学报, 2021, 43(4): 634–643. http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract18591.shtml DENG Hua-feng, QI Yu, LI Jian-lin, et al. Degradation mechanism of intermittent jointed sandstone under water-rock interaction[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(4): 634–643. (in Chinese) http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract18591.shtml
[23]	ZHANG X P, WONG L N Y. Crack initiation, propagation and coalescence in rock-like material containing two flaws: a numerical study based on bonded-particle model approach[J]. Rock Mechanics and Rock Engineering, 2013, 46(5): 1001–1021.
[24]	王根龙, 伍法权, 张茂省. 平面滑动型岩质边坡稳定性极限分析上限法[J]. 工程地质学报, 2011, 19(2): 176–180. https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201102005.htm WANG Gen-long, WU Fa-quan, ZHANG Mao-sheng. Method of upper bound limit analysis for plane sliding of rock slopes[J]. Journal of Engineering Geology, 2011, 19(2): 176–180. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201102005.htm
[25]	梁文辉. 锦屏一级水电站引渠边坡稳定性研究[D]. 成都: 西南交通大学, 2009. LIANG Wen-hui. Study on Stability of Approach Channel Slope for Jinping 1 Hydropower Station[D]. Chengdu: Southwest Jiaotong University, 2009. (in Chinese)
[26]	舒继森, 王兴中, 周毅勇. 岩石边坡中滑动面水压分布假设的改进[J]. 中国矿业大学学报, 2004, 33(5): 509–512. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD200405003.htm SHU Ji-sen, WANG Xing-zhong, ZHOU Yi-yong. Improving on assumption for water pressure distributing on failure surface in rock slope[J]. Journal of China University of Mining & Technology, 2004, 33(5): 509–512. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD200405003.htm
[27]	ZHAO L H, ZUO S, LI L, et al. System reliability analysis of plane slide rock slope using Barton-Bandis failure criterion[J]. International Journal of Rock Mechanics and Mining Sciences, 2016, 88: 1–11.