4. 极限平衡法安全系数计算原理

极限平衡法是矿山边坡稳定性分析中应用最广、工程经验最丰富的一类定量方法。它的核心思想并不复杂:把可能滑动的岩土体划分成若干竖向条块,对每个条块建立静力平衡(或力矩平衡)方程,再通过引入安全系数这一折减系数,使滑面上的抗剪强度与剪应力之间达到临界平衡状态。真正需要吃透的,是条分法的基本假定如何影响方程的可解性、安全系数如何定义、以及为什么必须迭代求解。

一、条分法为什么要“分条”

均质土坡的滑动面通常是曲线(圆弧或近似圆弧),滑面上的正应力和剪应力沿弧长连续变化,直接对整个滑体列平衡方程会得到含未知应力分布的积分方程,无法手算求解。条分法把滑体沿竖直方向切成 n 个条块,假定条块底面为滑面的一段,条块之间通过条间力相互作用。这样,每个条块底面的法向力 N 和剪切力 T 可以用条块自重、条间力和底面倾角表示,问题转化为有限个未知量的代数方程组。

以瑞典条分法(Fellenius 法)为例,其关键假定是:忽略条块之间的相互作用力,即条间力为零。此时每个条块独立平衡,底面法向力 Ni = Wicosαi(Wi 为条块自重,αi 为底面倾角),安全系数由抗滑力矩与滑动力矩之比给出:

F = Σ[c'·l_i + (W_i·cosα_i - u_i·l_i)·tanφ'] / Σ(W_i·sinα_i)

式中 c'、φ' 为有效抗剪强度指标,li 为条块底弧长,ui 为底面孔隙水压力。瑞典法满足整体力矩平衡,但不满足条块间的力平衡,对条间力较大的边坡偏于保守,安全系数常偏低 5%~20%。

二、从瑞典法到通用条分:假定决定方程可解性

更严格的条分法(毕肖普法、简布法、摩根斯顿-普赖斯法等)都保留条间力,但条间力的方向、大小和作用点是未知的。对 n 个条块,未知量包括各条块底面的 Ni、Ti,条间法向力 Ei、条间剪力 Xi,以及安全系数 F,未知量个数远多于可列的平衡方程数,属于超静定问题。因此必须补充假定:

  • 毕肖普法:假定条间力合力水平,即条间剪力 Xi = 0,只保留水平条间力 Ei,满足整体力矩平衡,得到含 F 的隐式方程;
  • 简布法:假定条间力作用点位置(推力线),通过迭代同时满足力和力矩平衡;
  • 摩根斯顿-普赖斯法:假定条间力方向用一个可调函数 X/E = λ·f(x) 描述,λ 为待定比例系数,严格满足全部平衡条件。

这些假定的本质,是在“方程数不足”与“工程可算”之间取得平衡。假定越少、越接近真实受力,计算越复杂,但对安全系数的修正越合理。工程实践中,圆弧滑面常用毕肖普法,任意形状滑面常用摩根斯顿-普赖斯法或简布法。

单条块受力(第 i 条) 滑面(圆弧段) W_i 自重 N_i T_i E_i X_i E_i+1 X_i+1 底面:N_i、T_i;条间:E_i、X_i;自重 W_i 瑞典法令 E_i = X_i = 0;毕肖普法令 X_i = 0 安全系数迭代求解逻辑 ① 假定初始 F₀(如 1.00) 划分条块,计算 W_i、l_i、α_i、u_i ② 由 F 反算条间力 / 底面力 T_i = [c'l_i + (N_i-u_i l_i)tanφ'] / F ③ 代入平衡方程求新 F′ 力矩平衡 / 力平衡方程解出 F′ ④ 收敛判断 |F′ − F| ≤ ε ε 取 0.001~0.005,满足则输出 F ⑤ 输出安全系数 F 与对应滑面 不收敛则回代 注:F 同时出现在抗剪强度折减与平衡方程两侧,故为隐式方程,需迭代;不同条分假定对应不同补充方程。

三、安全系数的定义与隐式方程

极限平衡法中安全系数的标准定义是抗剪强度折减系数:将滑面上的实际抗剪强度指标除以 F,使折减后的强度恰好与滑面上的剪应力平衡。即令 cm = c'/F、tanφm = tanφ'/F,代入莫尔-库仑准则,条块底面可提供的剪力为:

T_i = (c'·l_i + (N_i - u_i·l_i)·tanφ') / F

把 Ti 代入条块或整体平衡方程后,F 会同时出现在方程两端——因为 Ni 本身又通过条间力和 Ti 与 F 耦合。以简化毕肖普法为例,其安全系数表达式为:

F = Σ{[c'l_i + (W_i - u_i l_i)tanφ'] / m_αi} / Σ W_i sinα_i
m_αi = cosα_i + (sinα_i·tanφ') / F

可以看到,F 出现在 mαi 中,方程无法直接解出 F,只能迭代。这正是极限平衡法“隐式求解”的根源,也是很多初学者困惑的地方:不是公式写错,而是安全系数定义方式决定了它必须循环逼近。

四、迭代求解的实际操作与收敛判断

以某矿山排土场边坡为例(示例参数,仅用于说明计算流程):坡高 40 m,坡比 1:2.0,采用简化毕肖普法,划分 20 个条块,土体重度 γ = 19 kN/m³,有效黏聚力 c' = 12 kPa,有效内摩擦角 φ' = 22°,地下水位以下条块考虑孔隙水压力。计算步骤如下:

  1. 假定 F0 = 1.00,按条块几何计算各条块 Wi、li、αi、ui;
  2. 用当前 F 计算每个条块的 mαi,进而求得新的 F′;
  3. 比较 |F′ − F|,若大于收敛容差 ε(工程上常取 0.001~0.005),令 F = F′ 返回第 2 步;
  4. 收敛后输出 F,并检查是否满足 F ≥ 设计安全系数(如《非煤露天矿边坡工程技术规范》对一般边坡常要求 F ≥ 1.15~1.25,具体按边坡等级确定)。

实际迭代通常 4~8 次即可收敛。若出现振荡不收敛,常见原因是条块划分过粗、底面倾角跨越滑面拐点导致 αi 突变,或孔隙水压力取值与水位线不一致。此时应加密条块、检查滑面几何与水位线的一致性,而不是盲目减小容差。

条块划分建议:滑面曲率大的区段条块宜加密(宽度 1~3 m),曲率平缓段可适当放宽;条块底面不应跨越地层分界或滑面拐点,否则 αi 与强度参数的代表性会失真。

五、不同条分方法的选用与结果判读

方法核心假定满足的平衡条件适用场景
瑞典条分法忽略条间力整体力矩平衡均质圆弧滑面、初步估算
简化毕肖普法条间力水平(X=0)整体力矩平衡圆弧滑面、常规矿山边坡
简布法假定推力线位置力与力矩平衡任意滑面、需校核条间力
摩根斯顿-普赖斯法条间力方向函数 X/E=λf(x)严格全部平衡复杂地层、任意滑面、复核

判读结果时要注意:同一滑面用不同方法算出的 F 存在差异是正常的,瑞典法通常最低,摩根斯顿-普赖斯法更接近真实。工程上应以规范指定的方法为准,并用另一种方法做交叉校核;若两者差值超过 10%,应检查条间力假定、强度参数取值和孔隙水压力分布是否合理,而不是简单取平均。安全系数的意义在于相对比较和规范对标,脱离方法与参数谈 F 值高低没有工程意义。

安全系数小于 1 并不等于边坡“一定已经滑动”,它表示按所取强度参数和滑面,抗剪强度已不足以维持极限平衡;反之 F 大于规范值也不代表绝对安全,因为滑面搜索是否遗漏、强度参数是否代表现场、地下水与爆破振动等工况是否覆盖,都会改变结论。计算结果是判据之一,不是唯一判据。

掌握条分法的关键,是理解“假定—方程—安全系数定义—迭代”这条逻辑链:假定决定方程是否可解,安全系数的折减定义决定方程是否隐式,隐式方程决定必须迭代。把这三层关系理清,再结合具体滑面几何与参数,才能把极限平衡法用对、用好。