两种边坡稳定性分析方法比较研究.pdfVIP

第10卷 第10期 中 国 水 运 Vol.10 No.10 2010年 10月 China Water Transport October 2010 摘 要：以某水电工程岩质高边坡做为实例，将强度折减理论与 FLAC3D 软件相结合，通过有限差分程序 FLAC3D 软件来模拟分析其稳定性。并与极限平衡方法的分析结果对比，探索两种方法的差异性与结果的可靠性，为确定适 合工程建设实际的岩质边坡稳定分析方法提出了有益的参考。 关键词：强度折减法；极限平衡法；边坡稳定性 中图分类号：P642.1 文献标识码：A 文章编号：1006-7973（2010）10-0197-03 一、引言 目前，国内在建和待建的大型水电工程大多坐落在西南、 西北高山峡谷地区。我国的水电建设面临着一系列高边坡稳 定问题。在现代岩土工程和科学技术的新成就的支持下，确 定适合工程建设实际的岩质边坡稳定分析方法，是摆在水利 [1] 水电工程技术人员面前的任务 。 目前工程实践中岩质边坡稳定性定量分析主要有三种 方法：解析法（最常用的是极限平衡法）、数值方法和概率 法。极限平衡法是最常用的解析法，它是在边坡滑动面确 定的情况下，根据滑裂面上抗滑力和滑动力比值直接计算 安全系数，此外，关键块理论也属于这样的确定性分析方 法。数值方法则是借助计算机进行数值分析（例如有限元、 快速拉格朗日分析法、离散元、块体元和 DDA 等）从而确 图1 边坡剖面 定边坡的位移场和应力场，再用超载法、强度折减法等使 三、强度折减法 边坡处于极限状态，从而间接得到安全系数。这种方法同 强度折减系数法的基本原理是将坡体强度参数凝聚力 c 时可以考虑位移协调条件和岩体本构关系等。概率法是将 和内摩擦角 f 值同时除以一个安全系数 K，得到一组新的 c 、 k 概率统计理论被引用到边坡岩体的稳定性分析中来，它通 fk 值，然后作为新的资料参数输入，再进行试算，当计算不 过现场调查，以获得影响边坡稳性影响因素的多个样本， 收敛时，对应的 K 被称为坡体的最小稳定安全系数，此时坡 然后进行统计分析，求出它们各自的概率分布及其特征参 体达到极限状态，发生剪切破坏，同时可得到坡体的破坏滑 数，再利用某种可靠性分析方法，来求解边坡岩体的破坏 动面。 [2] 概率即可靠度 。 FLAC3D （ Three Dimensional Fast Lagrangian 文中选用某水电工程岩质高边坡做为实例，采用强度折 Analysis of Continua）是美国 Itasca Consulting Goup 减法和极限平衡法对岩质高边坡的稳定性进行对比分析。 lnc 开发的三维快速拉格朗日分析程序。该程序能较好地模 二、边坡工程地质条件 拟地质材料在达到强度极限或屈服极限时发生的破坏的力学 模型宽约为 700m，高约为 700m。 行为，特别适用于分析渐进破坏和失稳。 基岩以中粒结构的灰白色、微红色黑云二长花岗岩为主， 文中利用 FLAC3D，采用“二分法”[3]实现强度折减法， 并有辉绿岩脉（β）、花岗细晶岩脉、闪长岩脉等各类脉岩穿 求解安全系数。 插发育于花岗岩中，尤以辉绿岩脉分布较多。建模过程中考 所建计算模型节点为 29,646 个，单元为 24,005 个。 虑了岩体中对边坡稳定影响较大的几个岩脉。 模型的边界条件：模型四周法向约束，底部固定约束，顶部 根据岩体风化特点，岸坡岩体由表向内可划分为全风化 自由，仅受重力作用。 [4] 带、强风化带、弱风化带、微风化—新鲜岩体。岩体风化的 研究表明，随着剪胀角的增大，安全系数也逐渐增大 。 水平、垂直分带性明显。 不过，Vermeer 和 de Borst（1984 年）研究证明，一般土 边坡内无地下水分布。 体、岩石和混凝土的剪胀角要比它们的摩擦角小得多，且通 边坡剖面如图 1 所示。 常在 0°～20°内变化[5] 。因此，剪胀角对强度折减法计算 收稿日期：2010-06-11 作者简介：马玉岩（1987-），男，黑龙江绥化人，武汉大学水资源与水电工程科学国家重点实验室水利水电工程施工与 管理专业硕士研究生，主要研究方向为岩土边坡工程研究以及结构设计。 198 中 国 水 运 第10卷 结果的影响是有一定范围的。 文中在分析过程中选用 Mohr-Coulomb 本构模型，采 用非关联流动法则（剪胀角ψ 0）求解，并使用剪胀角 ψ 20°的情况加以校核。 计算模型所需的材料参数如表 1 所示，岩脉选用全风化 带参数。 表1 岩体物理力学参数 （a）瑞典圆弧法 （b）简化Bishop 法 3 对模型进行计算，当数值计算达稳定后，可获得剪切应 变增量云图，如图 2 所示。 （c