三维软硬互层边坡的破坏模式与稳定性研究

工程科学学报，第39卷，第2期：182-189.2017年2月 Chinese Journal of Engineering,Vol.39,No.2:182-189,February 2017 DOI:10.13374/j.issn2095-9389.2017.02.003;http://journals.ustb.edu.cn 三维软硬互层边坡的破坏模式与稳定性研究 姚文敏)，胡斌四，余海兵12)，李华舟)，和大钊 1)中国地质大学工程学院，武汉4300742)武汉市政工程设计研究院有限责任公司.武汉430015 区通信作者，E-mail:40275083@qq.com 摘要采用FLAC3D强度折减法，研究在岩层倾角、岩层与边坡走向夹角变化时三维软硬互层边坡的稳定性状况，并对其破 坏模式进行辨识与归纳分析.结果表明：边坡破坏模式的判别应综合考虑岩层的倾角大小、岩层走向与边坡走向的夹角大小 及坡面上的剪出条件：当岩层与边坡走向夹角B90°时，边坡的破坏模式趋势为塑流-拉裂、滑移-弯曲、弯 曲-拉裂：边坡稳定性系数随走向夹角的增大先增加后减小，B=90°时最大，且α：越大，稳定性系数峰值越大：顺向时随着岩层 倾角的增大，边坡的破坏模式变化趋势为蠕滑-压致拉裂、滑移-拉裂、滑移-弯曲、弯曲一拉裂，稳定性系数变化先减小后增 大，存在一最不利岩层倾角，其对应的稳定性系数最小：反向坡的破坏模式变化趋势为蠕滑-压致拉裂和弯曲-拉裂，稳定性系 数逐渐增加. 关键词边坡：破坏模式：稳定性系数：强度折减法 分类号TU457 Numerical analysis of the failure modes and stability of 3D slopes with interbreeding of soft and hard rocks YAO Wen-min'),HU Bin),YU Hai-bing2),LI Hua-zhou,HE Da-zhao) 1)Faculty of Engineering.China University of Geosciences,Wuhan 430074,China 2)Wuhan Municipal Engineering Design Research Institute Co.Ltd.,Wuhan 430015,China Corresponding author,E-mail:40275083@qq.com ABSTRACT FLAC3D strength reduction method had been adopted to research on the stability of 3D slopes with interbedding of soft and hard rocks,when the dip angle of stratum and the angle between the strike direction of slope and stratum changed,and the failure modes were identified and analyzed.The results show that:to identify the failure modes of slopes,the angle of the rock stratum,the angle between the strike direction of slope and stratum and the cut off condition on the slope surface should be considered.When B (the intersecting angle between the strike direction of slope and that of strata)90,the failure modes change from plastic flowing-ripping to sliding-bending,and then become bending- ripping.With the increase of the angle between the strike direction of slope and stratum,the slope stability coefficient increases at first and then decreases,it reaches the maximum when B=90,and the greater the a is,the greater the coefficient peak is.For the dip slope,with the increase of the angle of stratum,the trend of failure modes is:creeping-pressure induced ripping,sliding-ripping, sliding-bending,and bending-ripping.The stability coefficient decreases first and then increases,and there is a most unfavorable an- gle of the stratum,the corresponding stability coefficient is the least.As for the anti-dip slope,the trend of failure modes is creeping- pressure induced ripping and bending-ripping,and the stability coefficient increases gradually. KEY WORDS slopes;failure modes;stability;strength reduction method 收稿日期：2016-05-16

工程科学学报,第 39 卷,第 2 期:182鄄鄄189,2017 年 2 月 Chinese Journal of Engineering, Vol. 39, No. 2: 182鄄鄄189, February 2017 DOI: 10. 13374 / j. issn2095鄄鄄9389. 2017. 02. 003; http: / / journals. ustb. edu. cn 三维软硬互层边坡的破坏模式与稳定性研究 姚文敏1) , 胡 斌1) 苣 , 余海兵1,2) , 李华舟1) , 和大钊1) 1) 中国地质大学工程学院, 武汉 430074 2) 武汉市政工程设计研究院有限责任公司, 武汉 430015 苣 通信作者, E鄄mail:40275083@ qq. com 收稿日期: 2016鄄鄄05鄄鄄16 摘 要 采用 FLAC 3D强度折减法,研究在岩层倾角、岩层与边坡走向夹角变化时三维软硬互层边坡的稳定性状况,并对其破 坏模式进行辨识与归纳分析. 结果表明:边坡破坏模式的判别应综合考虑岩层的倾角大小、岩层走向与边坡走向的夹角大小 及坡面上的剪出条件;当岩层与边坡走向夹角 茁 90毅时,边坡的破坏模式趋势为塑流鄄鄄拉裂、滑移鄄鄄弯曲、弯 曲鄄鄄拉裂;边坡稳定性系数随走向夹角的增大先增加后减小,茁 = 90毅时最大,且 琢 越大,稳定性系数峰值越大;顺向时随着岩层 倾角的增大,边坡的破坏模式变化趋势为蠕滑鄄鄄压致拉裂、滑移鄄鄄拉裂、滑移鄄鄄弯曲、弯曲鄄鄄拉裂,稳定性系数变化先减小后增 大,存在一最不利岩层倾角,其对应的稳定性系数最小;反向坡的破坏模式变化趋势为蠕滑鄄鄄压致拉裂和弯曲鄄鄄拉裂,稳定性系 数逐渐增加. 关键词 边坡; 破坏模式; 稳定性系数; 强度折减法 分类号 TU457 Numerical analysis of the failure modes and stability of 3D slopes with interbreeding of soft and hard rocks YAO Wen鄄min 1) , HU Bin 1) 苣 , YU Hai鄄bing 1,2) , LI Hua鄄zhou 1) , HE Da鄄zhao 1) 1) Faculty of Engineering, China University of Geosciences, Wuhan 430074, China 2) Wuhan Municipal Engineering Design & Research Institute Co. Ltd. , Wuhan 430015, China 苣 Corresponding author, E鄄mail: 40275083@ qq. com ABSTRACT FLAC 3D strength reduction method had been adopted to research on the stability of 3D slopes with interbedding of soft and hard rocks, when the dip angle of stratum and the angle between the strike direction of slope and stratum changed, and the failure modes were identified and analyzed. The results show that: to identify the failure modes of slopes, the angle of the rock stratum, the angle between the strike direction of slope and stratum and the cut off condition on the slope surface should be considered. When 茁 (the intersecting angle between the strike direction of slope and that of strata) 90毅, the failure modes change from plastic flowing鄄鄄ripping to sliding鄄鄄 bending, and then become bending鄄鄄 ripping. With the increase of the angle between the strike direction of slope and stratum, the slope stability coefficient increases at first and then decreases, it reaches the maximum when 茁 = 90毅, and the greater the 琢 is, the greater the coefficient peak is. For the dip slope, with the increase of the angle of stratum, the trend of failure modes is: creeping鄄鄄 pressure induced ripping, sliding鄄鄄 ripping, sliding鄄鄄bending, and bending鄄鄄ripping. The stability coefficient decreases first and then increases, and there is a most unfavorable an鄄 gle of the stratum, the corresponding stability coefficient is the least. As for the anti鄄dip slope, the trend of failure modes is creeping鄄鄄 pressure induced ripping and bending鄄鄄ripping, and the stability coefficient increases gradually. KEY WORDS slopes; failure modes; stability; strength reduction method

姚文敏等：三维软硬互层边坡的破坏模式与稳定性研究 ·183· 我国西南地区煤系、红层、变质岩等地层分布广泛 连续介质考虑.边坡模型尺寸取为230m×50m×100 发育有大量软硬互层边坡.软硬互层边坡由于其互层岩 m,坡高50m,计算模型的左边界到坡脚的距离为坡高 层强度等力学性质差异显著的特性，差异风化作用更加 的1.5倍，右边界至坡顶的距离为坡高的2.5倍：边界 明显，较均质边坡更容易发生失稳破坏，如宣汉天台乡滑 条件为底部固定约束，模型四周采用法向约束，上部为 坡、重庆武隆鸡尾山滑坡、青海高原龙穆尔沟红层滑坡等. 自由边界；岩层与边坡走向夹角B范围为0°~180° 目前，已有部分学者对二维软硬互层边坡破坏模 (取0°、10°、30°、50°、70°、90°、100°、120°、140°、160°、 式、破坏机理等展开了研究.Hok等山将岩质边坡的180)，岩层倾角a范围为0°~90°（取0°、10°、20°、 破坏类型划分为4类.张倬元等)提出了5种斜坡变 30°、40°45°、60°、70°、90°)，边坡角度为50°.二维软 形破坏模式.林杭等]基于FLAC3D数值模拟软件，运 硬互层边坡破坏模式的研究结果显示]，在坡脚有软 用强度折减法研究层状边坡在不同岩层倾角对应下边 弱岩层剪出的条件下，当软硬岩层厚度在一定范围内 坡的破坏模式和稳定性情况.胡斌等以西南地区某 变化时，软硬岩厚度的不同组合对边坡破坏模式无显 边坡为例，研究了软硬互层边坡的崩塌机理.宋娅芬 著影响，对于岩层倾角较小的边坡稳定性影响也较小， 等[)采用模型试验方法模拟工程开挖及雨水浸润软化 因此本次研究中设定硬岩岩层厚度为8m,软岩岩层厚 过程作用下缓倾软硬互层边坡的变形破坏模式.宋玉 度为2m,软硬岩层间隔分布.采用有限元软件ANSYS 环6]以鲤鱼塘溢洪道边坡为例，对西南地区软硬互层 建立边坡模型，采用自编的ANSYS-LAC"接口程序 边坡变形破坏模式和稳定性分析进行了研究.Hungr 将模型导入到FLAC”软件中进行计算.考虑到单元 等)在滑坡调查统计的基础上，研究了边坡发生滑移、 数量对于计算效率的影响，首先经过试算，保证网格疏 倾倒等破坏的破坏机制.部分学者[s-]采用数值模拟 密程度满足计算精度要求.所采用的模型如图1 及模型试验等方法就岩层与边坡走向夹角对层状边坡 所示. 稳定性情况的影响做了部分研究.随着计算机技术及 数值模拟技术的发展，现在的边坡稳定性分析越来越 多地采用三维模型进行数值模拟研究，Dawson等[s]基 于强度折减法，研究了边界条件、模型形状对三维边坡 稳定性的影响规律.高永涛等采用离散元软件 3DEC研究了结构面倾角与边坡角变化情况下顺层岩 质边坡稳定性的发展规律，解释了顺层岩质边坡的变 形破坏特征.周元辅等]基于强度折减法对三维边 坡的失稳判据进行了研究，并建议在滑带明确的情况 下采用塑性区贯通率增量突变判据，在滑带未知的情 况下采用位移增量突变判据.Wei等[u]采用FLACD 强度折减法和极限平衡法对三维均质边坡的破坏模式 图1 FLAC3D数值模型示意图（单位：m) 进行了对比分析，结果显示两种计算方法的结果相吻 Fig.1 Numerical model in FLAC3D (unit:m) 合，研究表明含软弱夹层边坡夹层部位局部网格的划 分对边坡的失稳规律有很大的影响.刘红帅等)基 稳定性系数是评价边坡稳定性的重要指标，在模 于ABAOUS软件，采用有限元强度折减法研究了边界 型的计算过程中，坡顶临空面小规模软弱夹层是边坡 约束效应对于边坡稳定性的影响.以上诸多学者的研 稳定性最差的部分（如图1中a所示），软件的计算会 究极大地促进了强度折减法及边坡破坏模式方面的研 将该部分的稳定性系数作为整体边坡的稳定性系数， 究进展，但由于三维条件下软硬互层边坡的复杂性，对 但这显然是不合理的].因此，为了减小其对边坡整 于多因素影响下尤其是岩层与边坡走向夹角变化情况 体边坡稳定性系数的影响，在建模的过程中将坡顶小 下软硬互层边坡的变形破坏模式尚缺乏系统的研究. 规模软岩材料也赋值为硬岩从而使边坡的稳定性系数 本文拟采用FLAC”强度折减法，研究在岩层倾 更接近于真实值. 角、岩层与边坡走向夹角变化时三维软硬互层边坡的 本文是在边坡模型尺寸、边坡角度、软硬岩厚度固 变形破坏模式和稳定性状况. 定的条件下，调整岩层倾角及岩层与边坡走向夹角来 研究岩层倾角、岩层与边坡走向夹角变化时三维软硬 1数值模拟方法 互层边坡的稳定性状况，当各项条件改变时，边坡的破 1.1数值模型的建立 坏模式及稳定性系数变化趋势有所不同.在调整岩层 数值模型中，硬岩和软岩均采用实体单元模拟，按 倾角时，以边坡模型左下角点（即坐标原点）为圆心旋

姚文敏等: 三维软硬互层边坡的破坏模式与稳定性研究 我国西南地区煤系、红层、变质岩等地层分布广泛, 发育有大量软硬互层边坡. 软硬互层边坡由于其互层岩 层强度等力学性质差异显著的特性,差异风化作用更加 明显,较均质边坡更容易发生失稳破坏,如宣汉天台乡滑 坡、重庆武隆鸡尾山滑坡、青海高原龙穆尔沟红层滑坡等. 目前,已有部分学者对二维软硬互层边坡破坏模 式、破坏机理等展开了研究. Hoek 等[1] 将岩质边坡的 破坏类型划分为 4 类. 张倬元等[2]提出了 5 种斜坡变 形破坏模式. 林杭等[3] 基于 FLAC 3D数值模拟软件,运 用强度折减法研究层状边坡在不同岩层倾角对应下边 坡的破坏模式和稳定性情况. 胡斌等[4]以西南地区某 边坡为例,研究了软硬互层边坡的崩塌机理. 宋娅芬 等[5]采用模型试验方法模拟工程开挖及雨水浸润软化 过程作用下缓倾软硬互层边坡的变形破坏模式. 宋玉 环[6]以鲤鱼塘溢洪道边坡为例,对西南地区软硬互层 边坡变形破坏模式和稳定性分析进行了研究. Hungr 等[7]在滑坡调查统计的基础上,研究了边坡发生滑移、 倾倒等破坏的破坏机制. 部分学者[8鄄鄄12] 采用数值模拟 及模型试验等方法就岩层与边坡走向夹角对层状边坡 稳定性情况的影响做了部分研究. 随着计算机技术及 数值模拟技术的发展,现在的边坡稳定性分析越来越 多地采用三维模型进行数值模拟研究,Dawson 等[13]基 于强度折减法,研究了边界条件、模型形状对三维边坡 稳定性的影响规律. 高永涛等[14] 采用离散元软件 3DEC 研究了结构面倾角与边坡角变化情况下顺层岩 质边坡稳定性的发展规律,解释了顺层岩质边坡的变 形破坏特征. 周元辅等[15] 基于强度折减法对三维边 坡的失稳判据进行了研究,并建议在滑带明确的情况 下采用塑性区贯通率增量突变判据,在滑带未知的情 况下采用位移增量突变判据. Wei 等[16] 采用 FLAC 3D 强度折减法和极限平衡法对三维均质边坡的破坏模式 进行了对比分析,结果显示两种计算方法的结果相吻 合,研究表明含软弱夹层边坡夹层部位局部网格的划 分对边坡的失稳规律有很大的影响. 刘红帅等[17] 基 于 ABAQUS 软件,采用有限元强度折减法研究了边界 约束效应对于边坡稳定性的影响. 以上诸多学者的研 究极大地促进了强度折减法及边坡破坏模式方面的研 究进展,但由于三维条件下软硬互层边坡的复杂性,对 于多因素影响下尤其是岩层与边坡走向夹角变化情况 下软硬互层边坡的变形破坏模式尚缺乏系统的研究. 本文拟采用 FLAC 3D 强度折减法,研究在岩层倾 角、岩层与边坡走向夹角变化时三维软硬互层边坡的 变形破坏模式和稳定性状况. 1 数值模拟方法 1郾 1 数值模型的建立 数值模型中,硬岩和软岩均采用实体单元模拟,按 连续介质考虑. 边坡模型尺寸取为 230 m 伊 50 m 伊 100 m,坡高 50 m,计算模型的左边界到坡脚的距离为坡高 的 1郾 5 倍,右边界至坡顶的距离为坡高的 2郾 5 倍;边界 条件为底部固定约束,模型四周采用法向约束,上部为 自由边界;岩层与边坡走向夹角 茁 范围为 0毅 ~ 180毅 (取 0毅、10毅、30毅、50毅、70毅、90毅、100毅、120毅、140毅、160毅、 180毅),岩层倾角 琢 范围为 0毅 ~ 90毅 (取 0毅、10毅、20毅、 30毅、40毅、45毅、60毅、70毅、90毅),边坡角度为 50毅. 二维软 硬互层边坡破坏模式的研究结果显示[18] ,在坡脚有软 弱岩层剪出的条件下,当软硬岩层厚度在一定范围内 变化时,软硬岩厚度的不同组合对边坡破坏模式无显 著影响,对于岩层倾角较小的边坡稳定性影响也较小, 因此本次研究中设定硬岩岩层厚度为8 m,软岩岩层厚 度为 2 m,软硬岩层间隔分布. 采用有限元软件 ANSYS 建立边坡模型,采用自编的 ANSYS鄄鄄 FLAC 3D接口程序 将模型导入到 FLAC 3D软件中进行计算. 考虑到单元 数量对于计算效率的影响,首先经过试算,保证网格疏 密程度 满 足 计 算 精 度 要 求. 所 采 用 的 模 型 如 图 1 所示. 图 1 FLAC 3D数值模型示意图(单位:m) Fig. 1 Numerical model in FLAC 3D (unit: m) 稳定性系数是评价边坡稳定性的重要指标,在模 型的计算过程中,坡顶临空面小规模软弱夹层是边坡 稳定性最差的部分(如图 1 中 a 所示),软件的计算会 将该部分的稳定性系数作为整体边坡的稳定性系数, 但这显然是不合理的[8] . 因此,为了减小其对边坡整 体边坡稳定性系数的影响,在建模的过程中将坡顶小 规模软岩材料也赋值为硬岩从而使边坡的稳定性系数 更接近于真实值. 本文是在边坡模型尺寸、边坡角度、软硬岩厚度固 定的条件下,调整岩层倾角及岩层与边坡走向夹角来 研究岩层倾角、岩层与边坡走向夹角变化时三维软硬 互层边坡的稳定性状况,当各项条件改变时,边坡的破 坏模式及稳定性系数变化趋势有所不同. 在调整岩层 倾角时,以边坡模型左下角点(即坐标原点)为圆心旋 ·183·

·184· 工程科学学报，第39卷，第2期 转岩层：调整岩层与边坡走向夹角时，以纵坐标轴(z 表1岩体物理力学参数表 轴)为轴旋转岩层 Table I Physio-mechanical parameters of rock mass 1.2计算方法 弹性模量/泊松内聚力/内摩擦角/重度/抗拉强度/ 岩体 FLAC”软件基于三维显式有限差分原理，采用 GPa比MPa ((kN-m-3)MPa “混合离散法”来模拟材料的塑性破坏和塑性流动，可 硬岩7.00 0.21.00 35 25 0.500 采用非线性本构模型求解大变形问题，因此在岩土工 软岩。 0.68 0.30.02 20 20 0.016 程领域应用广泛. 强度折减法于20世纪70年代由Zienkiewiczl]首 构判断边坡的剪出条件，以边坡临界破坏（即计算不 收敛)时的位移云图、速度矢量图进行分析，图中不同 先提出，其计算原理为对岩土体抗剪强度参数c和P 同时除以折减系数K进行折减，然后进行计算，通过 颜色代表不同的水平位移值，颜色越深（棕色一橙 不断调整折减系数直至边坡处于临界破坏状态，此时 色一黄色一绿色一蓝色)，位移值越大：箭头表示速度 的折减系数K即为边的稳定性系数： 矢量，主要代表岩土体运动方向，参考刘汉超对于坡体 结构的分类方法[]，对不同坡体结构进行分类，依据 'tan g'tang K (1) 张倬元提出的5种层状边坡变形破坏模式[)，对计算 式中，黏聚力c和内摩擦角p为折减前的岩土体抗剪 所得边坡变形破坏模式进辨识（目前三维软硬岩层状 强度参数，'、p'为折减后的抗剪强度参数，K为折减后 边坡的破坏模式尚缺乏统一的定量评价标准，本文在 边坡处于临界状态的稳定性系数 数值模拟计算结果的基础上，采用定性的方法对边坡 随着计算机技术的发展，目前强度折减法可应用 的破坏模式进行辨识)，将边坡结构与破坏模式归纳 于各类数值模拟软件中，在破坏判据与计算精度方面 如表2所示，并选取各破坏模式中典型示例分述 也有了长足的进步[20].结合数值模拟软件的强度折 2.1蠕滑-压致拉裂 减法不需要事先假定滑动面，根据数值模拟软件的不 如图2(a)所示，对于近水平层状边坡，受到自重 同，可对边坡的应力-应变、位移发展特征及其渐进性 影响，边坡坡顶临空面附近首先产生向下的滑移，由于 破坏过程进行研究 软岩的蠕滑特性，使硬岩产生沿软岩的卸荷回弹，在坡 本次计算选用能同时考虑关联流动拉伸屈服和非 脚及坡面上岩层面产生张拉裂缝，由于上部岩体的蠕 关联流动剪切屈服的Mohr-Coulomb屈服准则，其表达 滑及压应力作用，下部岩体顺软岩层面向临空面方向 式如下。 产生滑移，坡脚及附近岩体被挤出.这类滑动具有压 剪切屈服准则： 剪性质，该类边坡的破坏模式属于蠕滑-压致拉裂. f=I1sinp/3-ccos p+√J2(cos0。+sin0,/5)=0. 如图2(b)所示，对于近直立近正交斜向倾内层状 (2) 边坡，由于岩层走向与边坡走向大于90°，岩体不具有 拉伸屈服准则： 向临空面方向滑移的岩层结构，同近水平层状边坡变 f=03-0 (3) 形破坏机理基本一致，边坡发生蠕滑-压致拉裂破坏. 式中：f(分别表示剪切和拉伸屈服函数；山，、J2分别为 2.2滑移-拉裂 应力张量第一不变量和应力偏量第二不变量：。为应 如图3(a)所示，对于缓倾角、中倾角顺向层状边 力罗德角：c、P分别为黏聚力和内摩擦角；σ为抗拉强 坡、缓倾角近顺向层状边坡及中倾角、陡倾角近正交斜 度：03为第三主应力. 向倾外层状边坡，坡脚或坡面上有软岩出露，具备剪出 强度折减法的失稳判别标准目前主要有数值计算 条件，岩体在重力的作用下沿软硬岩层界面发生滑移， 的不收敛、塑性区的贯通、特征部位位移的突变)]. 坡顶产生张裂缝，滑体各部分位移及速度方向基本一 本次以数值计算不收敛（不平衡力比率小于105）作 致.随着走向夹角的增大，边坡软岩在重力的作用下 为失稳判据. 沿软硬岩层界面向岩层倾向方向发生滑移，由于侧限 1.3模型参数 作用的原因，滑动方向有发生变化的趋势，但因为走向 本次计算以沪蓉高速公路西昌一恩施段出露的三 夹角较小，滑体滑动速度方向并无明显的改变，因此剪 叠系巴东组粉砂质泥岩和灰岩为软、硬岩层代表，结合 出口位置、堆积体大小沿边坡走向方向无明显差异，边 勘察报告和室内试验得到模型计算所需岩体物理力学 坡发生顺软硬岩层界面的滑移-拉裂破坏. 参数如表1所示. 如图3(b)所示，对于缓倾角、中倾角斜向倾外层 状边坡，随着走向夹角的增大，滑体在滑动时速度方向 2边坡变形破坏模式分析 发生较明显的改变，沿软硬岩层界面向边坡临空面滑 对于数值计算所得结果，综合边坡的岩层分布结 出，剪出口沿边坡走向方向不在同一高度，坡脚一侧由

工程科学学报,第 39 卷,第 2 期 转岩层;调整岩层与边坡走向夹角时,以纵坐标轴( z 轴)为轴旋转岩层. 1郾 2 计算方法 FLAC 3D软件基于三维显式有限差分原理,采用 “混合离散法冶来模拟材料的塑性破坏和塑性流动,可 采用非线性本构模型求解大变形问题,因此在岩土工 程领域应用广泛. 强度折减法于 20 世纪 70 年代由 Zienkiewicz [19]首 先提出,其计算原理为对岩土体抗剪强度参数 c 和 渍 同时除以折减系数 K 进行折减,然后进行计算,通过 不断调整折减系数直至边坡处于临界破坏状态,此时 的折减系数 K 即为边的稳定性系数: c忆 = c K ,tan 渍忆 = tan 渍 K . (1) 式中,黏聚力 c 和内摩擦角 渍 为折减前的岩土体抗剪 强度参数,c忆、渍忆为折减后的抗剪强度参数,K 为折减后 边坡处于临界状态的稳定性系数. 随着计算机技术的发展,目前强度折减法可应用 于各类数值模拟软件中,在破坏判据与计算精度方面 也有了长足的进步[20鄄鄄22] . 结合数值模拟软件的强度折 减法不需要事先假定滑动面,根据数值模拟软件的不 同,可对边坡的应力鄄鄄应变、位移发展特征及其渐进性 破坏过程进行研究. 本次计算选用能同时考虑关联流动拉伸屈服和非 关联流动剪切屈服的 Mohr鄄鄄Coulomb 屈服准则,其表达 式如下。 剪切屈服准则: f = I1 sin 渍/ 3 - ccos 渍 + J2 (cos 兹滓 + sin 兹滓 / 3) = 0. (2) 拉伸屈服准则: f t = 滓3 - 滓 t . (3) 式中:f、f t 分别表示剪切和拉伸屈服函数;I1 、J2 分别为 应力张量第一不变量和应力偏量第二不变量;兹滓 为应 力罗德角;c、渍 分别为黏聚力和内摩擦角;滓 t 为抗拉强 度;滓3 为第三主应力. 强度折减法的失稳判别标准目前主要有数值计算 的不收敛、塑性区的贯通、特征部位位移的突变[15] . 本次以数值计算不收敛(不平衡力比率小于 10 - 5 ) 作 为失稳判据. 1郾 3 模型参数 本次计算以沪蓉高速公路西昌—恩施段出露的三 叠系巴东组粉砂质泥岩和灰岩为软、硬岩层代表,结合 勘察报告和室内试验得到模型计算所需岩体物理力学 参数如表 1 所示. 2 边坡变形破坏模式分析 对于数值计算所得结果,综合边坡的岩层分布结 表 1 岩体物理力学参数表 Table 1 Physio鄄mechanical parameters of rock mass 岩体 弹性模量/ GPa 泊松 比 内聚力/ MPa 内摩擦角/ (毅) 重度/ (kN·m -3 ) 抗拉强度/ MPa 硬岩 7郾 00 0郾 2 1郾 00 35 25 0郾 500 软岩 0郾 68 0郾 3 0郾 02 20 20 0郾 016 构判断边坡的剪出条件,以边坡临界破坏(即计算不 收敛)时的位移云图、速度矢量图进行分析,图中不同 颜色代表不同的水平位移值,颜色越深 ( 棕色—橙 色—黄色—绿色—蓝色),位移值越大;箭头表示速度 矢量,主要代表岩土体运动方向,参考刘汉超对于坡体 结构的分类方法[23] ,对不同坡体结构进行分类,依据 张倬元提出的 5 种层状边坡变形破坏模式[2] ,对计算 所得边坡变形破坏模式进辨识(目前三维软硬岩层状 边坡的破坏模式尚缺乏统一的定量评价标准,本文在 数值模拟计算结果的基础上,采用定性的方法对边坡 的破坏模式进行辨识),将边坡结构与破坏模式归纳 如表 2 所示,并选取各破坏模式中典型示例分述. 2郾 1 蠕滑鄄鄄压致拉裂 如图 2(a)所示,对于近水平层状边坡,受到自重 影响,边坡坡顶临空面附近首先产生向下的滑移,由于 软岩的蠕滑特性,使硬岩产生沿软岩的卸荷回弹,在坡 脚及坡面上岩层面产生张拉裂缝,由于上部岩体的蠕 滑及压应力作用,下部岩体顺软岩层面向临空面方向 产生滑移,坡脚及附近岩体被挤出. 这类滑动具有压 剪性质,该类边坡的破坏模式属于蠕滑鄄鄄压致拉裂. 如图 2(b)所示,对于近直立近正交斜向倾内层状 边坡,由于岩层走向与边坡走向大于 90毅,岩体不具有 向临空面方向滑移的岩层结构,同近水平层状边坡变 形破坏机理基本一致,边坡发生蠕滑鄄鄄压致拉裂破坏. 2郾 2 滑移鄄鄄拉裂 如图 3(a)所示,对于缓倾角、中倾角顺向层状边 坡、缓倾角近顺向层状边坡及中倾角、陡倾角近正交斜 向倾外层状边坡,坡脚或坡面上有软岩出露,具备剪出 条件,岩体在重力的作用下沿软硬岩层界面发生滑移, 坡顶产生张裂缝,滑体各部分位移及速度方向基本一 致. 随着走向夹角的增大,边坡软岩在重力的作用下 沿软硬岩层界面向岩层倾向方向发生滑移,由于侧限 作用的原因,滑动方向有发生变化的趋势,但因为走向 夹角较小,滑体滑动速度方向并无明显的改变,因此剪 出口位置、堆积体大小沿边坡走向方向无明显差异,边 坡发生顺软硬岩层界面的滑移鄄鄄拉裂破坏. 如图 3( b)所示,对于缓倾角、中倾角斜向倾外层 状边坡,随着走向夹角的增大,滑体在滑动时速度方向 发生较明显的改变,沿软硬岩层界面向边坡临空面滑 出,剪出口沿边坡走向方向不在同一高度,坡脚一侧由 ·184·

姚文敏等：三维软硬互层边坡的破坏模式与稳定性研究 ·185· 表2三维软硬互层边坡变形破坏模式 Table 2 Failure modes of 3D slopes with interbedding of soft and hard rocks 岩层倾角，a 岩层走向与边坡走向的夹角，B 边坡类型 破坏模式 0°≤a<10° 0°≤B<180 近水平层状边坡 蠕滑-压致拉裂 0°≤β<30° 缓倾角顺向、近顺向层状边坡 滑移-拉裂 30°≤B<70° 缓倾角斜向倾外层状边坡 转动型滑移-拉裂 10°≤a≤20 70°≤B≤90° 缓倾角近正交斜向倾外层状边坡 塑流-拉裂 90°<B<180° 缓倾角斜向倾内层状边坡 塑流-拉裂 B=180° 缓倾角反向层状边坡 蠕滑-压致拉裂 B=0° 中倾角顺向层状边坡 滑移-拉裂 0°<B<70° 中倾角近顺向、斜向倾外层状边坡 转动型滑移-拉裂 20°<a≤50° 70°≤B≤90° 中倾角近正交斜向倾外层状边坡 滑移-拉裂 90°<B<160 中倾角斜向倾内层状边坡 滑移-弯曲 160°<B≤180° 中倾角近反向、反向层状边坡 弯曲-拉裂 B=0° 陡倾角顺向层状边坡 滑移-弯曲 0°<B<70° 陡倾角近顺向、斜向倾外层状边坡 转动型滑移-弯曲 50°<a≤70° 70°≤B<90° 陡倾角近正交斜向倾外层状边坡 滑移-拉裂 90°≤B≤180° 陡倾角斜向倾内、近反向、反向层状边坡 弯曲-拉裂 B=0° 近直立顺向层状边坡 弯曲-拉裂 0°<β<70° 近直立近顺向、斜向倾外层状边坡 转动型弯曲-拉裂 70°<a≤90° 70°≤B<90° 近直立近正交斜向倾外层状边坡 弯曲-拉裂 90°≤B≤120° 近直立近正交斜向倾内层状边坡 蠕滑-压致拉裂 120°<B≤180° 近直立斜向倾内、近反向、反向层状边坡 弯曲-拉裂 (b) 图2蠕滑-压致拉裂破坏模式.(a)a=0°；(b)B=120°，a= 90°：(c)二维示意图(B=120°，a=90) 图3滑移-拉裂破坏模式.（a)B=10°，a=20°；(b)B=50°， Fig.2 Failure modes of creeping-pressure induced ripping:(a)a= a=20°；(c)B=70°，a=45°；(d)二维示意图(B=70°，a= 0°；(b)B=120°，a=90°；(c）two-dimensional sketch map(B= 45°) 120°，a=90°) Fig.3 Failure mode of sliding-ripping:(a）B=l0°，d=20°；(b) 于滑体的压应力作用被挤出，边坡发生顺软硬岩层界 B=50°，a=20°；(c)B=70°，a=45:(d)two-dimensional 面的转动型滑移-拉裂破坏. sketch map(B=70°，a=45) 如图3()所示，对于缓倾角、中倾角近正交斜向 2.3塑流-拉裂 倾外层状边坡，边坡坡面上有软岩出露，具备剪出条 如图4(a)所示，对于缓倾角近正交斜向倾外层状 件，在重力和侧限作用下，坡顶产生张裂缝，岩体沿软 边坡，上部岩体在重力的作用下向下产生变形，软岩在 硬岩层界面向边坡临空面方向发生滑移，滑体为以软 坡面上出露，成为良好的剪出条件，下部岩体在上覆岩 岩层面为底的楔形体，滑体各部分位移及速度方向基 体的压力作用下，产生塑性流动并向临空面方向挤出， 本一致，边坡产生顺软硬岩界面的滑移-拉裂破坏. 从而使上覆岩体及上部岩体被拉裂缝切割后，临空侧

姚文敏等: 三维软硬互层边坡的破坏模式与稳定性研究 表 2 三维软硬互层边坡变形破坏模式 Table 2 Failure modes of 3D slopes with interbedding of soft and hard rocks 岩层倾角,琢 岩层走向与边坡走向的夹角,茁 边坡类型 破坏模式 0毅臆琢 < 10毅 0毅臆茁 < 180毅 近水平层状边坡 蠕滑鄄鄄压致拉裂 0毅臆茁 < 30毅 缓倾角顺向、近顺向层状边坡 滑移鄄鄄拉裂 30毅臆茁 < 70毅 缓倾角斜向倾外层状边坡 转动型滑移鄄鄄拉裂 10毅臆琢臆20毅 70毅臆茁臆90毅 缓倾角近正交斜向倾外层状边坡 塑流鄄鄄拉裂 90毅 < 茁 < 180毅 缓倾角斜向倾内层状边坡 塑流鄄鄄拉裂 茁 = 180毅 缓倾角反向层状边坡 蠕滑鄄鄄压致拉裂 茁 = 0毅 中倾角顺向层状边坡 滑移鄄鄄拉裂 0毅 < 茁 < 70毅 中倾角近顺向、斜向倾外层状边坡 转动型滑移鄄鄄拉裂 20毅 < 琢臆50毅 70毅臆茁臆90毅 中倾角近正交斜向倾外层状边坡 滑移鄄鄄拉裂 90毅 < 茁 < 160毅 中倾角斜向倾内层状边坡 滑移鄄鄄弯曲 160毅 < 茁臆180毅 中倾角近反向、反向层状边坡 弯曲鄄鄄拉裂 茁 = 0毅 陡倾角顺向层状边坡 滑移鄄鄄弯曲 50毅 < 琢臆70毅 0毅 < 茁 < 70毅 陡倾角近顺向、斜向倾外层状边坡 转动型滑移鄄鄄弯曲 70毅臆茁 < 90毅 陡倾角近正交斜向倾外层状边坡 滑移鄄鄄拉裂 90毅臆茁臆180毅 陡倾角斜向倾内、近反向、反向层状边坡 弯曲鄄鄄拉裂 茁 = 0毅 近直立顺向层状边坡 弯曲鄄鄄拉裂 0毅 < 茁 < 70毅 近直立近顺向、斜向倾外层状边坡 转动型弯曲鄄鄄拉裂 70毅 < 琢臆90毅 70毅臆茁 < 90毅 近直立近正交斜向倾外层状边坡 弯曲鄄鄄拉裂 90毅臆茁臆120毅 近直立近正交斜向倾内层状边坡 蠕滑鄄鄄压致拉裂 120毅 < 茁臆180毅 近直立斜向倾内、近反向、反向层状边坡 弯曲鄄鄄拉裂 图 2 蠕滑鄄鄄压致拉裂破坏模式 郾 (a) 琢 = 0毅; (b) 茁 = 120毅, 琢 = 90毅; (c) 二维示意图(茁 = 120毅, 琢 = 90毅) Fig. 2 Failure modes of creeping鄄鄄pressure induced ripping: (a) 琢 = 0毅; (b) 茁 = 120毅, 琢 = 90毅; (c) two鄄dimensional sketch map ( 茁 = 120毅, 琢 = 90毅) 于滑体的压应力作用被挤出,边坡发生顺软硬岩层界 面的转动型滑移鄄鄄拉裂破坏. 如图 3(c)所示,对于缓倾角、中倾角近正交斜向 倾外层状边坡,边坡坡面上有软岩出露,具备剪出条 件,在重力和侧限作用下,坡顶产生张裂缝,岩体沿软 硬岩层界面向边坡临空面方向发生滑移,滑体为以软 岩层面为底的楔形体,滑体各部分位移及速度方向基 本一致,边坡产生顺软硬岩界面的滑移鄄鄄拉裂破坏. 图 3 滑移鄄鄄拉裂破坏模式 郾 (a) 茁 = 10毅, 琢 = 20毅; (b) 茁 = 50毅, 琢 = 20毅; (c) 茁 = 70毅, 琢 = 45毅; ( d) 二维示意图( 茁 = 70毅, 琢 = 45毅) Fig. 3 Failure mode of sliding鄄鄄ripping: (a) 茁 = 10毅, 琢 = 20毅; (b) 茁 = 50毅, 琢 = 20毅; ( c) 茁 = 70毅, 琢 = 45毅; ( d) two鄄dimensional sketch map (茁 = 70毅, 琢 = 45毅) 2郾 3 塑流鄄鄄拉裂 如图 4(a)所示,对于缓倾角近正交斜向倾外层状 边坡,上部岩体在重力的作用下向下产生变形,软岩在 坡面上出露,成为良好的剪出条件,下部岩体在上覆岩 体的压力作用下,产生塑性流动并向临空面方向挤出, 从而使上覆岩体及上部岩体被拉裂缝切割后,临空侧 ·185·

·186· 工程科学学报，第39卷，第2期 的危岩体由于下伏岩体的挤出、风化、剥蚀和淘空导致 危岩体重心向外偏移，坡脚及其前方岩体发生大范围 的挤出与位移，当其倾倒力矩大于抗倾力矩时，块体就 发生倾倒破坏24]，边坡发生塑流-拉裂的破坏. 图5滑移-弯曲破坏模式.(a)B=0°，a=70°：(b)B=30° a=70°：(c)B=120°，a=45°：(d)二维示意图(B=30°，a= 70°) 图4塑流-拉裂破坏模式.(a)B=70°，a=20°：(b)B=140°， Fig.5 Failure mode of sliding-bending:(a)B=0°，a=70°：(b) a=20°：(c）二维示意图(B=140°，a=20°) B=30°，a=70°；(c）B=120°，a=45°；(d)tbwo-dimensional Fig.4 Failure mode of flowing-ripping:(a）B=70°，a=20°； sketch map(B=30°，a=70°) (b)B=140°，a=20°；(c）two-dimensional sketch map(B= 破坏 140°，a=20°） 2.5弯曲-拉裂 如图4(b)所示，对于缓倾角斜向倾内层状边坡， 如图6(a)所示，对于近直立顺向层状边坡，边坡 上部岩体在重力的作用下向下产生变形，软岩在坡面 坡面及坡脚不具备剪出条件，边坡坡顶临空面附近软 上出露，虽然岩层倾向坡体内部，但由于岩层夹角较 岩首先发生上部向临空面的弯曲变形，并进一步影响 小，仍能成为良好的剪出条件，同缓倾角近正交斜向倾 坡面和坡顶岩体，软硬岩界面产生较大的裂缝（尤其 外层状边坡变形破坏机理一致，边坡发生塑流-拉裂 是坡面上)，边坡发生弯曲-拉裂破坏 的破坏 (a 2.4滑移-弯曲 如图5(a)所示，对于陡倾角顺向层状边坡，边坡 坡面及坡脚不具备剪出条件，上部岩体在重力的作用 下沿软岩层面发生向下的滑移，坡顶产生张裂缝，但由 于坡脚及以下岩体的阻滑作用，岩体发生由坡脚向上， 沿软岩层面的弯曲，坡面上软硬岩界面有较大的张裂 缝，边坡发生滑移-弯曲破坏. 如图5(b)所示，对于陡倾角近顺向、斜向倾外层 状边坡，边坡坡面及坡脚不具备剪出条件，上部岩体在 重力的作用下沿软岩层面先岩层倾向方向发生向下的 图6弯曲-拉裂破坏模式.(a)B=0°，a=70°：(b)B=30° 滑移，坡顶产生张裂缝，但由于侧边界、坡脚及以下岩 a=70°：(c)B=120°，a=45°：(d)二维示意图(B=140°，a= 体的阻滑作用，岩体发生由坡脚向上、沿软岩层面的弯 70°) 曲，坡面上软硬岩层界面有较大的张裂缝，边坡发生转 Fig.6 Failure mode of bending-ripping:(a）B=0°，a=70°：(b) 动型滑移-弯曲破坏，同时由于侧边界的侧限作用，滑 B=30°，a=70°；(c)B=120°，a=45°：(d)two-dimensional sketch map(β=140°，a=70) 动体岩层倾向一侧弯曲程度更大. 如图5(c)所示，对于中倾角近反向层状边坡，上 如图6(b)所示，对于近直立近顺向、斜向倾外层 部岩体在重力的作用下，发生向斜坡临空面方向的滑 状边坡，边坡坡面及坡脚不具备剪出条件，边坡坡顶临 移，并沿软硬岩层界面挤出，从而使岩层产生弯曲，伴 空面附近软岩首先发生上部向临空面的弯曲变形，并 随着上部岩体向临空面产生弯曲，坡顶产生张裂缝，越 进一步影响坡面和坡顶岩体，软硬岩界面产生较大的 靠近下部的岩体在上覆岩体压力的作用下滑移距离越 裂缝（尤其是坡面上），边坡发生转动型弯曲-拉裂破 长，下部坡脚及其下部岩体由于滑动体的冲击产生弯 坏.由于侧边界的侧限作用，滑动体岩层倾向一侧弯 曲，坡脚前岩体大范围被挤出，边坡发生滑移-弯曲的 曲程度更大.进一步发展下部岩体将发生断裂，从而

工程科学学报,第 39 卷,第 2 期 的危岩体由于下伏岩体的挤出、风化、剥蚀和淘空导致 危岩体重心向外偏移,坡脚及其前方岩体发生大范围 的挤出与位移,当其倾倒力矩大于抗倾力矩时,块体就 发生倾倒破坏[24] ,边坡发生塑流鄄鄄拉裂的破坏. 图4 塑流鄄鄄拉裂破坏模式郾 (a) 茁 = 70毅, 琢 = 20毅; (b) 茁 = 140毅, 琢 = 20毅; (c) 二维示意图(茁 = 140毅, 琢 = 20毅) Fig. 4 Failure mode of flowing鄄鄄 ripping: ( a) 茁 = 70毅, 琢 = 20毅; (b) 茁 = 140毅, 琢 = 20毅; ( c) two鄄dimensional sketch map ( 茁 = 140毅, 琢 = 20毅) 如图 4( b)所示,对于缓倾角斜向倾内层状边坡, 上部岩体在重力的作用下向下产生变形,软岩在坡面 上出露,虽然岩层倾向坡体内部,但由于岩层夹角较 小,仍能成为良好的剪出条件,同缓倾角近正交斜向倾 外层状边坡变形破坏机理一致,边坡发生塑流鄄鄄 拉裂 的破坏. 2郾 4 滑移鄄鄄弯曲 如图 5(a)所示,对于陡倾角顺向层状边坡,边坡 坡面及坡脚不具备剪出条件,上部岩体在重力的作用 下沿软岩层面发生向下的滑移,坡顶产生张裂缝,但由 于坡脚及以下岩体的阻滑作用,岩体发生由坡脚向上, 沿软岩层面的弯曲,坡面上软硬岩界面有较大的张裂 缝,边坡发生滑移鄄鄄弯曲破坏. 如图 5( b)所示,对于陡倾角近顺向、斜向倾外层 状边坡,边坡坡面及坡脚不具备剪出条件,上部岩体在 重力的作用下沿软岩层面先岩层倾向方向发生向下的 滑移,坡顶产生张裂缝,但由于侧边界、坡脚及以下岩 体的阻滑作用,岩体发生由坡脚向上、沿软岩层面的弯 曲,坡面上软硬岩层界面有较大的张裂缝,边坡发生转 动型滑移鄄鄄弯曲破坏,同时由于侧边界的侧限作用,滑 动体岩层倾向一侧弯曲程度更大. 如图 5(c)所示,对于中倾角近反向层状边坡,上 部岩体在重力的作用下,发生向斜坡临空面方向的滑 移,并沿软硬岩层界面挤出,从而使岩层产生弯曲,伴 随着上部岩体向临空面产生弯曲,坡顶产生张裂缝,越 靠近下部的岩体在上覆岩体压力的作用下滑移距离越 长,下部坡脚及其下部岩体由于滑动体的冲击产生弯 曲,坡脚前岩体大范围被挤出,边坡发生滑移鄄鄄 弯曲的 图 5 滑移鄄鄄弯曲破坏模式 郾 ( a) 茁 = 0毅, 琢 = 70毅; ( b) 茁 = 30毅, 琢 = 70毅; (c) 茁 = 120毅, 琢 = 45毅; ( d) 二维示意图( 茁 = 30毅, 琢 = 70毅) Fig. 5 Failure mode of sliding鄄鄄bending: (a) 茁 = 0毅, 琢 = 70毅; (b) 茁 = 30毅, 琢 = 70毅; ( c) 茁 = 120毅, 琢 = 45毅; ( d) two鄄dimensional sketch map (茁 = 30毅, 琢 = 70毅) 破坏. 2郾 5 弯曲鄄鄄拉裂 如图 6(a)所示,对于近直立顺向层状边坡,边坡 坡面及坡脚不具备剪出条件,边坡坡顶临空面附近软 岩首先发生上部向临空面的弯曲变形,并进一步影响 坡面和坡顶岩体,软硬岩界面产生较大的裂缝(尤其 是坡面上),边坡发生弯曲鄄鄄拉裂破坏. 图 6 弯曲鄄鄄拉裂破坏模式 郾 ( a) 茁 = 0毅, 琢 = 70毅; ( b) 茁 = 30毅, 琢 = 70毅; (c) 茁 = 120毅, 琢 = 45毅; (d) 二维示意图( 茁 = 140毅, 琢 = 70毅) Fig. 6 Failure mode of bending鄄ripping: (a) 茁 = 0毅, 琢 = 70毅; (b) 茁 = 30毅, 琢 = 70毅; ( c) 茁 = 120毅, 琢 = 45毅; ( d) two鄄dimensional sketch map (茁 = 140毅, 琢 = 70毅) 如图 6( b)所示,对于近直立近顺向、斜向倾外层 状边坡,边坡坡面及坡脚不具备剪出条件,边坡坡顶临 空面附近软岩首先发生上部向临空面的弯曲变形,并 进一步影响坡面和坡顶岩体,软硬岩界面产生较大的 裂缝(尤其是坡面上),边坡发生转动型弯曲鄄鄄 拉裂破 坏. 由于侧边界的侧限作用,滑动体岩层倾向一侧弯 曲程度更大. 进一步发展下部岩体将发生断裂,从而 ·186·

姚文敏等：三维软硬互层边坡的破坏模式与稳定性研究 ·187· 使上部岩体发生倾倒破坏现象. 稳定性系数均在走向夹角α=90时达到最大值，且随 如图6（℃）所示，对于中倾角、陡倾角、近直立斜向 着岩层倾角的增大，边坡稳定性系数峰值逐渐增大 倾内层状边坡（包括反向层状边坡），上部岩体在重力 对于近水平边坡，随着走向夹角的增大，稳定性系 的作用下向临空面方向发生弯曲，坡顶及软硬岩层面 数逐渐增加，当走向夹角为70°时与水平层状边坡稳 产生张裂缝，越靠近下部的岩体在上覆岩体压力的作 定性系数基本相等.当走向夹角大于90°后随着其继 用下滑移距离越长，上部岩体的弯曲程度最大，坡脚前岩 续增大，稳定性系数小幅增加. 体大范围被挤出，边坡发生弯曲-拉裂的破坏.进一步发 对于缓倾角边坡，随着走向夹角的增大，稳定性系 展当坡顶岩体首先发生断裂，滑体可能产生倾倒破坏 数先有一小幅的下降段，之后稳定性系数逐渐增加. 对于变形破坏模式的分析总结如下： 当走向夹角大于90°后随着其继续增大，稳定性系数 (1)当岩层走向与边坡走向夹角小于90时，随着 小幅减小 岩层倾角的增大，边坡的破坏模式由蠕滑-压致拉裂、 对于中倾角边坡，当为斜向倾外边坡时，随着走向 塑流拉裂、滑移-拉裂逐渐转变为滑移-弯曲、弯曲-拉 夹角的增大，稳定性系数有一波动段，稳定性系数先增 裂.当岩层倾角小于坡角时，边坡坡面上有软弱层出 加后减小再增加.当走向夹角继续增大，稳定性系数 露，为潜在的剪出条件.边坡上部岩体产生变形，变形 逐渐减小 逐渐积累，由于坡面剪出条件的不同，产生蠕滑-压致 对于陡倾角边坡，当走向夹角小于90°时，随着走 拉裂、塑流-拉裂、滑移-拉裂的破坏.当岩层倾角大于 向夹角的增大，稳定性系数曲线呈先增大、后减小再增 坡角时，边坡坡面上没有软弱层出露，不具备剪出条 大的趋势：当走向夹角大于90°时，随着走向夹角的增 件，边坡上部岩体的变形产生圆弧形的滑移面，岩体产 大，稳定性系数逐渐降低. 生向临空面方向的弯曲，边坡发生滑移-弯曲的破坏； 对于近直立边坡，边坡以弯曲-拉裂的破坏模式 当岩层倾角进一步增大接近直立，边坡坡顶岩体的变 为主，随着走向夹角的增大，稳定性系数曲线呈两边 形转变为向临空面方向的弯曲，并进一步影响坡面和 低、中间高的近似抛物线形. 坡顶岩体，软硬岩层界面产生较大的裂缝（尤其是坡 面上)，边坡发生弯曲-拉裂破坏 4.0 (2)当岩层走向与边坡走向及夹角大于90°时，随 着岩层倾角的增大，边坡的破坏模式由塑流-拉裂向 滑移-弯曲、弯曲-拉裂转变.当岩层倾角较小时，上部 2.5 岩体在重力的作用下向下产生变形，软岩在坡面上出 露，虽然岩层倾向坡体内部，但由于岩层夹角较小，仍 2.0 一■一地层02。一地层10 能形成良好的剪出口条件，下部岩体在上覆岩体的压 15 一◆一地层20°，-地层40° 力作用下，产生塑性流动并向临空面方向挤出，继续变 一◆一地层60°一←一地层90° 0 形，上部临空侧危岩体将发生倾倒变形，边坡产生塑 20 406080100120140160.180 流-拉裂的破坏.当岩层倾角较大时，上部岩体发生向 岩层与边坡走向夹角) 斜坡临空面方向的滑移，并沿软硬岩层界面挤出，从而 图7斜交坡稳定性系数统变化图 使岩层产生弯曲，坡顶产生张裂缝，越靠近下部的岩体 Fig.7 Stability coefficient of diagonal slopes 在上覆岩体压力的作用下滑移距离越长，下部坡脚及 3.2非斜交坡稳定性分析 其下部岩体由于滑动体的冲击产生弯曲，坡脚前岩体 非斜交坡（即岩层走向与边坡走向夹角为0°或 大范围被挤出，边坡发生滑移-弯曲的破坏.当岩层倾 180°)的稳定性系数变化如图8所示，其中对于顺向坡 角进一步增大，边坡上部岩体的变形不再以滑移为主， 来说（图8中红色曲线），稳定性系数呈现先减小、后 而转变为向临空面方向的弯曲，坡顶及软硬岩层面产 增大的变化趋势，存在一最不利岩层倾角，其对应的稳 生张裂缝，越靠近下部的岩体在上覆岩体压力的作用 定性系数最小，本次研究中对应的最不利岩层倾角 下滑移距离越长，上部岩体的弯曲程度最大，坡脚前岩 α=40°时，这与二维软硬互层边坡稳定性的研究结论 体大范围被挤出，边坡发生弯曲-拉裂的破坏 相一致[1.此外，林杭等)、高永涛等、黄洪波 等[]的研究也认为，对于顺层层状边坡，岩层倾角对 3稳定性分析 边坡稳定性影响较大，存在一最不利岩层倾角，该最不 3.1斜交坡稳定性分析 利岩层倾角的大小受岩体抗剪强度指标、坡角、坡体结 从图7中可以直观地看出，各岩层倾角的边坡的 构等因素的影响.当α=0°时，边坡发生蠕滑-压致拉

姚文敏等: 三维软硬互层边坡的破坏模式与稳定性研究 使上部岩体发生倾倒破坏现象. 如图 6(c)所示,对于中倾角、陡倾角、近直立斜向 倾内层状边坡(包括反向层状边坡),上部岩体在重力 的作用下向临空面方向发生弯曲,坡顶及软硬岩层面 产生张裂缝,越靠近下部的岩体在上覆岩体压力的作 用下滑移距离越长,上部岩体的弯曲程度最大,坡脚前岩 体大范围被挤出,边坡发生弯曲鄄鄄拉裂的破坏. 进一步发 展当坡顶岩体首先发生断裂,滑体可能产生倾倒破坏. 对于变形破坏模式的分析总结如下: (1)当岩层走向与边坡走向夹角小于 90毅时,随着 岩层倾角的增大,边坡的破坏模式由蠕滑鄄鄄 压致拉裂、 塑流拉裂、滑移鄄鄄拉裂逐渐转变为滑移鄄鄄弯曲、弯曲鄄鄄拉 裂. 当岩层倾角小于坡角时,边坡坡面上有软弱层出 露,为潜在的剪出条件. 边坡上部岩体产生变形,变形 逐渐积累,由于坡面剪出条件的不同,产生蠕滑鄄鄄压致 拉裂、塑流鄄鄄拉裂、滑移鄄鄄拉裂的破坏. 当岩层倾角大于 坡角时,边坡坡面上没有软弱层出露,不具备剪出条 件,边坡上部岩体的变形产生圆弧形的滑移面,岩体产 生向临空面方向的弯曲,边坡发生滑移鄄鄄 弯曲的破坏; 当岩层倾角进一步增大接近直立,边坡坡顶岩体的变 形转变为向临空面方向的弯曲,并进一步影响坡面和 坡顶岩体,软硬岩层界面产生较大的裂缝(尤其是坡 面上),边坡发生弯曲鄄鄄拉裂破坏. (2)当岩层走向与边坡走向及夹角大于 90毅时,随 着岩层倾角的增大,边坡的破坏模式由塑流鄄鄄 拉裂向 滑移鄄鄄弯曲、弯曲鄄鄄拉裂转变. 当岩层倾角较小时,上部 岩体在重力的作用下向下产生变形,软岩在坡面上出 露,虽然岩层倾向坡体内部,但由于岩层夹角较小,仍 能形成良好的剪出口条件,下部岩体在上覆岩体的压 力作用下,产生塑性流动并向临空面方向挤出,继续变 形,上部临空侧危岩体将发生倾倒变形,边坡产生塑 流鄄鄄拉裂的破坏. 当岩层倾角较大时,上部岩体发生向 斜坡临空面方向的滑移,并沿软硬岩层界面挤出,从而 使岩层产生弯曲,坡顶产生张裂缝,越靠近下部的岩体 在上覆岩体压力的作用下滑移距离越长,下部坡脚及 其下部岩体由于滑动体的冲击产生弯曲,坡脚前岩体 大范围被挤出,边坡发生滑移鄄鄄弯曲的破坏. 当岩层倾 角进一步增大,边坡上部岩体的变形不再以滑移为主, 而转变为向临空面方向的弯曲,坡顶及软硬岩层面产 生张裂缝,越靠近下部的岩体在上覆岩体压力的作用 下滑移距离越长,上部岩体的弯曲程度最大,坡脚前岩 体大范围被挤出,边坡发生弯曲鄄鄄拉裂的破坏. 3 稳定性分析 3郾 1 斜交坡稳定性分析 从图 7 中可以直观地看出,各岩层倾角的边坡的 稳定性系数均在走向夹角 琢 = 90毅时达到最大值,且随 着岩层倾角的增大,边坡稳定性系数峰值逐渐增大. 对于近水平边坡,随着走向夹角的增大,稳定性系 数逐渐增加,当走向夹角为 70毅时与水平层状边坡稳 定性系数基本相等. 当走向夹角大于 90毅后随着其继 续增大,稳定性系数小幅增加. 对于缓倾角边坡,随着走向夹角的增大,稳定性系 数先有一小幅的下降段,之后稳定性系数逐渐增加. 当走向夹角大于 90毅后随着其继续增大,稳定性系数 小幅减小. 对于中倾角边坡,当为斜向倾外边坡时,随着走向 夹角的增大,稳定性系数有一波动段,稳定性系数先增 加后减小再增加. 当走向夹角继续增大,稳定性系数 逐渐减小. 对于陡倾角边坡,当走向夹角小于 90毅时,随着走 向夹角的增大,稳定性系数曲线呈先增大、后减小再增 大的趋势;当走向夹角大于 90毅时,随着走向夹角的增 大,稳定性系数逐渐降低. 对于近直立边坡,边坡以弯曲鄄鄄 拉裂的破坏模式 为主,随着走向夹角的增大,稳定性系数曲线呈两边 低、中间高的近似抛物线形. 图 7 斜交坡稳定性系数统变化图 Fig. 7 Stability coefficient of diagonal slopes 3郾 2 非斜交坡稳定性分析 非斜交坡(即岩层走向与边坡走向夹角为 0毅或 180毅)的稳定性系数变化如图 8 所示,其中对于顺向坡 来说(图 8 中红色曲线),稳定性系数呈现先减小、后 增大的变化趋势,存在一最不利岩层倾角,其对应的稳 定性系数最小,本次研究中对应的最不利岩层倾角 琢 = 40毅时,这与二维软硬互层边坡稳定性的研究结论 相一 致[18] . 此 外, 林 杭 等[3] 、 高 永 涛 等[14] 、 黄 洪 波 等[25]的研究也认为,对于顺层层状边坡,岩层倾角对 边坡稳定性影响较大,存在一最不利岩层倾角,该最不 利岩层倾角的大小受岩体抗剪强度指标、坡角、坡体结 构等因素的影响. 当 琢 = 0毅时,边坡发生蠕滑鄄鄄压致拉 ·187·

·188· 工程科学学报，第39卷，第2期 裂的破坏.随着岩层倾角的增大，边坡上部岩体下滑 可成为良好的剪出条件，随着岩层倾角的增大，边坡的 力增大，抗滑力减小：岩层在坡面上的出露也为边坡岩 破坏模式变化趋势为塑流-拉裂、滑移-弯曲、弯曲- 体的滑动提供了良好的剪出条件.因此，边坡的稳定 拉裂. 性系数逐渐减小.随着岩层倾角继续增大，边坡岩体 (3)各岩层倾角的边坡稳定性系数随着岩层走向 产生更大的位移量，导致岩层产生弯曲：但一方面坡面 与边坡走向夹角的增大基本呈先增加后减小的趋势， 上没有软硬岩层界面出露，岩体的抗弯强度较大，另一 在走向夹角α=90°时稳定性系数达到最大值，且随着 方面岩体较大的抗弯强度和坡脚岩体挤出堆积提供了 岩层倾角的增大，边坡稳定性系数峰值逐渐增大.对 一定的抗滑力，从而使得边坡的稳定性增强.当岩层 于中倾角斜向倾外边坡，随着走向夹角的增大，稳定性 近直立时，边坡岩体的变形方式由滑移向弯曲转变，岩 系数有一波动段，稳定性系数先增加后减小再增加 体的抗弯强度与抗拉强度均较大，因此边坡的稳定性 (4)顺向时随着岩层倾角的增大，坡面剪出条件 进一步增强 逐渐变差，边坡的破坏模式变化趋势为蠕滑-压致拉 4.0 裂、滑移-拉裂、滑移-弯曲、弯曲-拉裂，边坡的稳定性 3.5 系数变化呈现先减小，后增大的变化趋势，存在一最不 利岩层倾角，其对应的稳定性系数最小 (5)反向时由于岩层倾向坡体内部，随着岩层倾 2.5 角的增大，边坡的破坏模式变化趋势为蠕滑-压致拉 裂、弯曲-拉裂.边坡岩体主要在重力作用下产生弯 2.0 。一反向坡 曲，岩层倾角增大，边坡岩体弯曲方向由斜向下方逐渐 一·一顺向坡 转变为水平方向，弯曲力逐渐减小，因此边坡的稳定性 系数变化整体呈现逐渐增加的变化趋势. 0 40 60 80 岩层倾角) 参考文献 图8非斜交坡稳定性系数变化图 [1]Hoek E,Bray J W.Rock Slope Engineering.Washington DC:In- Fig.8 Stability coefficient of non diagonal slopes stitution of Mining Metallurgy,1977 反向坡的稳定性系数变化如图8中黑色曲线所 [2]Zhang Z Y,Wang S T,Wang L S.Principle of Engineering Geolo- 示，呈现先增加，后减小再增加的变化趋势，但整体呈 gy Analysis.Beijing:Geological Publishing House,1994 (张倬元，王士天，王兰生.工程地质分析原理.北京：地质 现增加的趋势.由于岩层倾向坡体内部，随着岩层倾 出版社，1994) 角的增大，边坡岩体的变形方式由蠕滑向弯曲转变 [3]Lin H,Cao P,Li J T,et al.Numerical analysis of failure modes 边坡岩体主要在重力作用下产生弯曲，岩层倾角增大， and stability of stratified rock slopes.Rock Soil Mech,2010,31 边坡岩体弯曲方向由斜向下方逐渐转变为水平方向， (10):3300 及弯曲力逐渐减小，因而边坡的稳定性逐渐增强 (林杭，曹平，李江腾，等.层状岩质边坡破坏模式及稳定性 的数值分析.岩土力学，2010,31(10)：3300) 4结论 [4]Hu B,Huang R Q.Collapse mechanism and treatment measures of slopes with interbeddings of soft and hard rocks.I Eng Geol, (1)边坡破坏模式的判别应综合考虑岩层的倾角 2009,17(2):200 大小、岩层走向与边坡走向的夹角大小及其组合形成 (胡斌，黄润秋.软硬岩互层边坡崩塌机理及治理对策研究 的坡面上的剪出条件.随着岩层倾角的变化，边坡岩 工程地质学报，2009,17(2)：200) 体的变形方式发生变化：随着坡面上剪出条件的变化， [5] Song Y F,Chen C X,Zheng Y,et al.Model experimental study 边坡岩体的破坏模式发生变化. of deformation failure mechanism of low-angled slopes with inter- bedding of soft and hard rocks.Rock Soil Mech,2015,36(2): (2)以张倬元提出的5种斜坡变形破坏模式为依 487 据，对不同类型边坡的变形破坏模式进行归纳总结，并 (宋娅芬，陈从新，郑允，等.缓倾软硬岩互层边坡变形破坏 对不同类型边坡的变形破坏机理进行了分析：当B90时，岩层倾角较小时，坡面上出露的软硬岩层仍 (宋玉环.西南地区软硬互层岩质边坡变形破坏模式及稳定

工程科学学报,第 39 卷,第 2 期 裂的破坏. 随着岩层倾角的增大,边坡上部岩体下滑 力增大,抗滑力减小;岩层在坡面上的出露也为边坡岩 体的滑动提供了良好的剪出条件. 因此,边坡的稳定 性系数逐渐减小. 随着岩层倾角继续增大,边坡岩体 产生更大的位移量,导致岩层产生弯曲;但一方面坡面 上没有软硬岩层界面出露,岩体的抗弯强度较大,另一 方面岩体较大的抗弯强度和坡脚岩体挤出堆积提供了 一定的抗滑力,从而使得边坡的稳定性增强. 当岩层 近直立时,边坡岩体的变形方式由滑移向弯曲转变,岩 体的抗弯强度与抗拉强度均较大,因此边坡的稳定性 进一步增强. 图 8 非斜交坡稳定性系数变化图 Fig. 8 Stability coefficient of non diagonal slopes 反向坡的稳定性系数变化如图 8 中黑色曲线所 示,呈现先增加,后减小再增加的变化趋势,但整体呈 现增加的趋势. 由于岩层倾向坡体内部,随着岩层倾 角的增大,边坡岩体的变形方式由蠕滑向弯曲转变. 边坡岩体主要在重力作用下产生弯曲,岩层倾角增大, 边坡岩体弯曲方向由斜向下方逐渐转变为水平方向, 及弯曲力逐渐减小,因而边坡的稳定性逐渐增强. 4 结论 (1)边坡破坏模式的判别应综合考虑岩层的倾角 大小、岩层走向与边坡走向的夹角大小及其组合形成 的坡面上的剪出条件. 随着岩层倾角的变化,边坡岩 体的变形方式发生变化;随着坡面上剪出条件的变化, 边坡岩体的破坏模式发生变化. (2)以张倬元提出的 5 种斜坡变形破坏模式为依 据,对不同类型边坡的变形破坏模式进行归纳总结,并 对不同类型边坡的变形破坏机理进行了分析:当 茁 90毅时,岩层倾角较小时,坡面上出露的软硬岩层仍 可成为良好的剪出条件,随着岩层倾角的增大,边坡的 破坏模式变化趋势为塑流鄄鄄 拉裂、滑移鄄鄄 弯曲、弯曲鄄鄄 拉裂. (3)各岩层倾角的边坡稳定性系数随着岩层走向 与边坡走向夹角的增大基本呈先增加后减小的趋势, 在走向夹角 琢 = 90毅时稳定性系数达到最大值,且随着 岩层倾角的增大,边坡稳定性系数峰值逐渐增大. 对 于中倾角斜向倾外边坡,随着走向夹角的增大,稳定性 系数有一波动段,稳定性系数先增加后减小再增加. (4)顺向时随着岩层倾角的增大,坡面剪出条件 逐渐变差,边坡的破坏模式变化趋势为蠕滑鄄鄄 压致拉 裂、滑移鄄鄄拉裂、滑移鄄鄄弯曲、弯曲鄄鄄拉裂,边坡的稳定性 系数变化呈现先减小,后增大的变化趋势,存在一最不 利岩层倾角,其对应的稳定性系数最小. (5)反向时由于岩层倾向坡体内部,随着岩层倾 角的增大,边坡的破坏模式变化趋势为蠕滑鄄鄄 压致拉 裂、弯曲鄄鄄拉裂. 边坡岩体主要在重力作用下产生弯 曲,岩层倾角增大,边坡岩体弯曲方向由斜向下方逐渐 转变为水平方向,弯曲力逐渐减小,因此边坡的稳定性 系数变化整体呈现逐渐增加的变化趋势. 参 考 文 献 [1] Hoek E, Bray J W. Rock Slope Engineering. Washington DC: In鄄 stitution of Mining & Metallurgy, 1977 [2] Zhang Z Y, Wang S T, Wang L S. Principle of Engineering Geolo鄄 gy Analysis. Beijing: Geological Publishing House, 1994 (张倬元, 王士天, 王兰生. 工程地质分析原理. 北京: 地质 出版社, 1994) [3] Lin H,Cao P, Li J T, et al. Numerical analysis of failure modes and stability of stratified rock slopes. Rock Soil Mech, 2010, 31 (10): 3300 (林杭, 曹平, 李江腾, 等. 层状岩质边坡破坏模式及稳定性 的数值分析. 岩土力学, 2010, 31(10): 3300) [4] Hu B, Huang R Q. Collapse mechanism and treatment measures of slopes with interbeddings of soft and hard rocks. J Eng Geol, 2009, 17(2): 200 (胡斌, 黄润秋. 软硬岩互层边坡崩塌机理及治理对策研究. 工程地质学报, 2009, 17(2): 200) [5] Song Y F, Chen C X, Zheng Y, et al. Model experimental study of deformation failure mechanism of low鄄angled slopes with inter鄄 bedding of soft and hard rocks. Rock Soil Mech, 2015, 36 (2): 487 (宋娅芬, 陈从新, 郑允, 等. 缓倾软硬岩互层边坡变形破坏 机制模型试验研究. 岩土力学, 2015, 36(2): 487) [6] Song Y H. Failure Models and Instability of Slopes with Soft and Hard Rock Interbedded Laminated Structure in Southwest China [Dissertation ]. Chengdu: Chengdu University of Technology, 2011 (宋玉环. 西南地区软硬互层岩质边坡变形破坏模式及稳定 ·188·

姚文敏等：三维软硬互层边坡的破坏模式与稳定性研究 ·189· 性研究一以鲤鱼塘水库溢洪道边坡为例[学位论文].成 analysis by the strength reduction and limit equilibrium methods. 都：成都理工大学，2011) Comput Geotech,2009,36(1-2):70 [7]Hungr O,Evans SG.The occurrence and classification of massive [17]Liu HS,Nian T K,Wan SS.Effect of boundary constraint con- rock slope failure.Felsbau,2004,22(2):16 dition on the stability analysis of 3D slope.J Jilin Univ Earth Sci [8]Lu Z M,Chen C X,Zuo B C,et al.Experimental research on Ed.2010,40(3):638 factors influencing stability of anti-dip layered slope.Rock Soil (刘红帅，年廷凯，万少石.三维边坡稳定性分析中的边界 Mech.2006.27(4):629 约束效应.吉林大学学报（地球科学版），2010,40(3)： (卢增木，陈从新，左保成，等.对影响逆倾层状边坡稳定性 638) 因素的模型试验研究.岩土力学，2006,27(4)：629) [18]Zheng Z Y,Yu H B,Xu H Q.Numerical analysis of failure [9]Yang Z W.A simplified algorithm of dip slope considers the angle modes and stability of soft and hard rock interbedded slope. between strike of slope and strata /The Academic Essays of Rock Yangtze River Sci Res Inst,2016,33(9):102 Mechanics and Engineering Society in Guizhou in 2013.Guiyang, (郑志勇，余海兵，徐海清.软硬岩互层边坡的破坏模式及 2013:234 稳定性研究.长江科学院院报，2016,33(9)：102) (杨泽伟.考虑岩层与坡面走向夹角的顺层岩坡稳定性简化 [19]Zienkiewicz O C.Humpheson C,Lewis R W.Associated and 算法/贵州省岩石力学与工程学会2013年学术年会论文 non-associated visco-plasticity and plasticity in soil mechanics. 集.贵阳，2013：234) Geotechnique,1975,25(4):671 [10]Feng J,Jiang N,Zhou D P.Influence of angle between strike of [20]Shen J Y,Karakus M.Three-dimensional numerical analysis for slope and strata on slope stability.Ind Constr.2008,38(10): rock slope stability using shear strength reduction method.Can 76 Geotech J.2014,51(2):164 (冯君，江南，周德培.顺层边坡岩层走向与边坡走向夹角 [21]Kelesoglu M K.The evaluation of three-dimensional effects on 对其稳定性的影响.工业建筑，2008,38(10)：76) slope stability by the strength reduction method.KSCE Civ [11]Bai Y F.Zhou D P.Feng J.Upper limit of angle between strikes Eng,2016,20(1):229 of slope and strata for dip slope.J Southucest Jiaotong Univ, [22]Chen G Q,Huang R Q,Shi Y C.et al.Stability analysis of 2005,40(3):326 slope based on dynamic and whole strength reduction methods. (白云峰，周德培，冯君.顺向坡岩层走向与边坡走向夹角 Chin J Rock Mech Eng,2014,33(2):243 的上限值.西南交通大学学报，2005,40(3)：326) (陈国庆，黄润秋，石豫川.等.基于动态和整体强度折减法 [12]Chai B,Yin K L.Influence of intersection angle between trend 的边坡稳定性分析.岩石力学与工程学报，2014,33(2)： of slope and strata on stability of bedding slope.Chin J Rock 243) Mech Eng,2009.28(3):628 [23]Liu HC.Chen M D,Li K W.The Geological Enrironment Eral- (柴波，殷坤龙.顺向坡岩层倾向与坡向夹角对斜坡稳定性 uation Research.Chengdu:Chengdu University of Science and 的影响.岩石力学与工程学报，2009,28(3)：628) Technology Press,1993 [13]Dawson E M,Roth W H,Drescher A.Slope stability analysis by (刘汉超，陈明东，黎克武，等.库区环境地质评价研究.成 strength reduction.Geotechnique,1999,49(6):835 都：成都科技大学出版社，1993) [14]Gao Y T,Xiao S,Wu S C,et al.Numerical simulation of the [24]Chen X T,Huang B L,Liu G N,et al.Analysis of deforming deformation and failure characteristics of consequent rock slopes mode transforming of gentle soft-hard interbed slopes in three gor- and their stability.Chin J Eng,2015,37(11):1403 ges reservoir.J Geol Hazards Environ Preserv,2009,20(2):57 (高永涛，肖术，吴顺川，等。顺层岩质边坡变形破坏特征及 (陈小婷，黄波林，刘广宁，等.三峡库区平缓层状软硬相间 稳定性数值模拟.工程科学学报，2015,37(11)：1403) 斜坡变形模式变化分析.地质灾害与环境保护，2009,20 [15]Zhou Y F,Deng J H,Cui Y L,et al.Instability criterion of (2):57) three-dimensional slope based on strength reduction method. [25]Huang H B.Stability Analysis of Stratified Rock Slope Disserta- Rock Soil Mech,2014,35(5):1430 tion].Hangzhou:Zhejiang University,2003 (周元轴，邓建辉，崔玉龙，等.基于强度折减法的三维边 (黄洪波.层状岩质边坡的稳定性分析[学位论文].杭州： 坡失稳判据.岩土力学，2014,35(5)：1430) 浙江大学，2003) [16]Wei W B,Cheng Y M,Li L.Three-dimensional slope failure

姚文敏等: 三维软硬互层边坡的破坏模式与稳定性研究 性研究———以鲤鱼塘水库溢洪道边坡为例[学位论文]. 成 都: 成都理工大学, 2011) [7] Hungr O, Evans S G. The occurrence and classification of massive rock slope failure. Felsbau, 2004, 22(2): 16 [8] Lu Z M, Chen C X, Zuo B C, et al. Experimental research on factors influencing stability of anti鄄dip layered slope. Rock Soil Mech, 2006, 27(4): 629 (卢增木, 陈从新, 左保成, 等. 对影响逆倾层状边坡稳定性 因素的模型试验研究. 岩土力学, 2006, 27(4): 629) [9] Yang Z W. A simplified algorithm of dip slope considers the angle between strike of slope and strata / / The Academic Essays of Rock Mechanics and Engineering Society in Guizhou in 2013. Guiyang, 2013: 234 (杨泽伟. 考虑岩层与坡面走向夹角的顺层岩坡稳定性简化 算法 / / 贵州省岩石力学与工程学会 2013 年学术年会论文 集. 贵阳, 2013: 234) [10] Feng J, Jiang N, Zhou D P. Influence of angle between strike of slope and strata on slope stability. Ind Constr, 2008, 38(10): 76 (冯君, 江南, 周德培. 顺层边坡岩层走向与边坡走向夹角 对其稳定性的影响. 工业建筑, 2008, 38(10): 76) [11] Bai Y F, Zhou D P, Feng J. Upper limit of angle between strikes of slope and strata for dip slope. J Southwest Jiaotong Univ, 2005, 40(3): 326 (白云峰, 周德培, 冯君. 顺向坡岩层走向与边坡走向夹角 的上限值. 西南交通大学学报, 2005, 40(3): 326) [12] Chai B, Yin K L. Influence of intersection angle between trend of slope and strata on stability of bedding slope. Chin J Rock Mech Eng, 2009, 28(3): 628 (柴波, 殷坤龙. 顺向坡岩层倾向与坡向夹角对斜坡稳定性 的影响. 岩石力学与工程学报, 2009, 28(3): 628) [13] Dawson E M, Roth W H, Drescher A. Slope stability analysis by strength reduction. G佴otechnique, 1999, 49(6):835 [14] Gao Y T, Xiao S, Wu S C, et al. Numerical simulation of the deformation and failure characteristics of consequent rock slopes and their stability. Chin J Eng, 2015, 37(11): 1403 (高永涛, 肖术, 吴顺川, 等. 顺层岩质边坡变形破坏特征及 稳定性数值模拟. 工程科学学报, 2015, 37(11): 1403) [15] Zhou Y F, Deng J H, Cui Y L, et al. Instability criterion of three鄄dimensional slope based on strength reduction method. Rock Soil Mech, 2014, 35(5): 1430 (周元辅, 邓建辉, 崔玉龙, 等. 基于强度折减法的三维边 坡失稳判据. 岩土力学, 2014, 35(5): 1430) [16] Wei W B, Cheng Y M, Li L. Three鄄dimensional slope failure analysis by the strength reduction and limit equilibrium methods. Comput Geotech, 2009, 36(1鄄2): 70 [17] Liu H S, Nian T K, Wan S S. Effect of boundary constraint con鄄 dition on the stability analysis of 3D slope. J Jilin Univ Earth Sci Ed, 2010, 40(3): 638 (刘红帅, 年廷凯, 万少石. 三维边坡稳定性分析中的边界 约束效应. 吉林大学学报( 地球科学版), 2010, 40 (3 ): 638) [18] Zheng Z Y, Yu H B, Xu H Q. Numerical analysis of failure modes and stability of soft and hard rock interbedded slope. J Yangtze River Sci Res Inst, 2016, 33(9): 102 (郑志勇, 余海兵, 徐海清. 软硬岩互层边坡的破坏模式及 稳定性研究. 长江科学院院报, 2016, 33(9): 102) [19] Zienkiewicz O C, Humpheson C, Lewis R W. Associated and non鄄associated visco鄄plasticity and plasticity in soil mechanics. G佴otechnique, 1975, 25(4): 671 [20] Shen J Y, Karakus M. Three鄄dimensional numerical analysis for rock slope stability using shear strength reduction method. Can Geotech J, 2014, 51(2): 164 [21] Kelesoglu M K. The evaluation of three鄄dimensional effects on slope stability by the strength reduction method. KSCE J Civ Eng, 2016, 20(1): 229 [22] Chen G Q, Huang R Q, Shi Y C, et al. Stability analysis of slope based on dynamic and whole strength reduction methods. Chin J Rock Mech Eng, 2014, 33(2): 243 (陈国庆, 黄润秋, 石豫川, 等. 基于动态和整体强度折减法 的边坡稳定性分析. 岩石力学与工程学报, 2014, 33 (2): 243) [23] Liu H C, Chen M D, Li K W. The Geological Environment Eval鄄 uation Research. Chengdu: Chengdu University of Science and Technology Press, 1993 (刘汉超, 陈明东, 黎克武, 等. 库区环境地质评价研究. 成 都: 成都科技大学出版社, 1993) [24] Chen X T, Huang B L, Liu G N, et al. Analysis of deforming mode transforming of gentle soft鄄hard interbed slopes in three gor鄄 ges reservoir. J Geol Hazards Environ Preserv, 2009, 20(2): 57 (陈小婷, 黄波林, 刘广宁,等. 三峡库区平缓层状软硬相间 斜坡变形模式变化分析. 地质灾害与环境保护, 2009, 20 (2): 57) [25] Huang H B. Stability Analysis of Stratified Rock Slope [Disserta鄄 tion]. Hangzhou: Zhejiang University, 2003 (黄洪波. 层状岩质边坡的稳定性分析[学位论文]. 杭州: 浙江大学, 2003) ·189·