3. 岩体强度参数选取与折减

边坡稳定性计算的第一步不是套公式,而是把岩体强度“翻译”成计算模型能接受的参数。工程中最常用的两条路径是:以黏聚力 c、内摩擦角 φ 为核心的 Mohr-Coulomb(M-C)准则,以及以地质强度指标 GSI、岩石单轴抗压强度 σci、扰动系数 D 为核心的 Hoek-Brown(H-B)准则。二者并非互斥,而是同一岩体强度的两种表达,关键在于参数来源可靠、换算路径清晰、折减方式与准则匹配。

一、Mohr-Coulomb 参数的获取

M-C 参数获取有三条典型途径,可靠性依次降低,工程中应优先采用高等级途径并交叉验证。

  • 试验直取:对含结构面较少的均质岩体,可通过三轴压缩试验获得一组(σ3,σ1)数据,在 τ-σ 平面上拟合强度包络线,截距为 c,倾角为 φ。完整岩块试验值往往偏高,需按岩体尺度折减。
  • 规范/经验取值:依据岩体质量分级(RMR、BQ、GSI)查表取 c、φ 范围,再结合边坡等级、破坏后果调整。此类取值应注明依据与折减系数,避免“拍脑袋”。
  • 反演校核:对已有变形或失稳迹象的边坡,用极限平衡反算滑面强度,作为正演取值的上限约束。

需要强调:c、φ 是“等效”参数,对节理岩体而言,它们综合反映了岩块与结构面的共同作用,不能简单等同于室内岩块试验值。对于层状、片理发育的边坡,还应区分沿结构面滑动与切岩滑动两种模式,分别给出参数。

二、Hoek-Brown 参数换算为 M-C 参数

H-B 准则的广义形式为:

σ1 = σ3 + σci · [ mb·(σ3/σci) + s ]^a

其中 mb = mi·exp[(GSI−100)/(28−14D)],s = exp[(GSI−100)/(9−3D)],a = 0.5 + (e−GSI/15 − e−20/3)/6。换算为等效 M-C 参数时,关键是选定一个“等效围压范围” σ3max,因为 H-B 包络线是曲线,等效 c、φ 随围压变化。常用做法是令 σ3max = 0.72·(σcm/γH)−0.91·σcm(σcm 为岩体强度,γ 为重度,H 为坡高),再在该区间内对 H-B 曲线做线性拟合,得到:

φ = arcsin[ (6·a·mb·(s + mb·σ3n)^(a−1)) / (2(1+a)(2+a) + 6·a·mb·(s + mb·σ3n)^(a−1)) ]
c = σci·[ (1+2a)s + (1−a)mb·σ3n ]·(s + mb·σ3n)^(a−1) / [ (1+a)(2+a)·√(1 + 6·a·mb·(s+mb·σ3n)^(a−1)/((1+a)(2+a))) ]

式中 σ3n = σ3max/σci。工程上更常用的是查表或软件内置换算,但必须核对所取 σ3max 是否与坡高、应力水平匹配——坡越高,等效围压区间越大,等效 φ 通常越低、c 相应变化,直接套用浅坡参数会高估深部岩体强度。

示例(假设值,仅用于说明流程):某边坡 H=120 m,γ=26 kN/m³,σci=60 MPa,GSI=45,mi=12,D=0.7。计算得 mb≈1.19、s≈0.0016、a≈0.51,取 σ3max≈2.1 MPa,换算得等效 φ≈31°、c≈0.42 MPa。该结果应再与现场结构面调查、反演结果比对后使用。

三、强度折减法中的折减路径

强度折减法(SRM)的核心是把强度参数按同一比例逐步下调,直至边坡达到极限平衡,此时的折减系数即安全系数 F。折减路径必须与所用准则一致:

  • M-C 模型:同时折减 c 与 tanφ,即 cF=c/F,tanφF=tanφ/F。注意是折减 tanφ 而非 φ 本身,二者在数值上不同,误用会系统性偏差。
  • H-B 模型:不能直接折减 c、φ,而应折减 σci 与 mb(或等效折减 GSI),保持 s、a 的物理含义,或按软件约定同步折减。不同软件对 H-B 折减的实现不同,建模前必须确认。

折减路径还涉及“是否同步折减变形参数”。SRM 判定失稳依靠强度储备耗尽后的位移突变或计算不收敛,弹性模量、泊松比一般不参与折减;若同步折减,会改变变形场,影响失稳判据。工程取值时,F 的接受标准应结合边坡等级:一般临时边坡可取 1.1~1.2,永久边坡 1.25~1.35,重要边坡或存在软弱结构面时适当提高,并以多种判据(位移突变、塑性区贯通、计算不收敛)相互印证。

岩体强度参数获取—换算—折减路径 M-C 参数来源 • 三轴试验拟合 c、φ • 岩体分级查表取值 • 失稳边坡反演校核 注意:等效参数≠岩块值 H-B 参数输入 σci、GSI、mi、D → mb、s、a 选定 σ3max 等效围压 坡越高,σ3max 越大 H-B → M-C 等效换算 在 σ3max 区间线性拟合 得到等效 c、φ σ3max 与坡高 H 匹配 结果需与反演/调查比对 曲线包络→分段等效 M-C 折减路径 cF = c / F tanφF = tanφ / F ✗ 不可直接折减 φ 角值 H-B 折减路径 折减 σci 与 mb 或等效折减 GSI 保持 s、a 物理含义 须核对软件实现约定 失稳判据与工程取值 位移突变 · 塑性区贯通 · 计算不收敛 → F;临时坡 1.1~1.2,永久坡 1.25~1.35

四、工程取值的关键控制点

第一,参数来源要“可追溯”。每一组 c、φ 或 GSI、σci 都应记录试验编号、分级依据、折减系数和采用理由,便于复核与反演修正。第二,注意尺度效应与结构面控制。室内岩块强度通常高于岩体强度,GSI 与 D 的取值直接决定 mb、s,扰动系数 D 对爆破、开挖卸荷区应取 0.7~1.0,对机械开挖可取 0.3~0.5,取值不当会使强度高估数倍。第三,折减路径与准则必须匹配,M-C 折 tanφ、H-B 折 σci 与 mb,混用会导致 F 失真。第四,安全系数不是唯一判据,应结合滑面位置、塑性区形态和位移量级综合判断;当 F 接近临界值时,参数微小变化即可改变滑面出口,需做敏感性分析。

常见误区:把室内岩块 c、φ 直接用于节理岩体边坡;用 H-B 换算时随意取 σ3max;在 M-C 模型中折减 φ 角值而非 tanφ;以及忽略 D 值对 mb、s 的指数影响。这些都会使安全系数偏离真实值,且偏差方向往往偏于不安全。

综上,参数选取与折减是一条“来源—换算—折减—判据”的完整链条:先明确岩体质量与结构面控制模式,再选择准则并获取参数,按准则匹配的路径折减,最后用多判据确认安全系数。只有链条上每一环都可追溯、可校核,边坡稳定性计算结果才具备工程决策价值。