2. Hoek-Brown准则的数学形式与参数体系

Hoek-Brown准则之所以能在岩体工程中被广泛采用,关键在于它用一组物理意义明确的参数,把“完整岩石强度”与“节理岩体强度”之间的折减关系表达了出来。本节从广义方程出发,逐一辨析 σci、mi、GSI、D、Ei 等参数的含义,并说明它们如何共同决定岩体的强度包络线。

一、从狭义到广义:方程的推导脉络

最初的Hoek-Brown准则(1980年)针对完整岩石提出,形式为:

σ1 = σ3 + σci · (mi · σ3 / σci + 1)^0.5

其中 σ1、σ3 分别为破坏时的最大、最小主应力(压为正),σci 为完整岩石单轴抗压强度,mi 为完整岩石的Hoek-Brown常数。该式只适用于完整岩石,无法描述节理岩体的强度折减。

为把节理、层理等结构面的影响纳入,Hoek与Brown引入地质强度指标 GSI 和扰动系数 D,得到广义Hoek-Brown准则:

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

其中三个岩体参数由 GSI、D、mi 换算得到:

mb = mi · exp[(GSI - 100) / (28 - 14D)]
s  = exp[(GSI - 100) / (9 - 3D)]
a  = 0.5 + (1/6) · [exp(-GSI/15) - exp(-20/3)]

这组换算式是广义准则的核心:GSI 越低、D 越大,mb 与 s 衰减越剧烈,岩体强度包络线整体下移。当 GSI=100、D=0 时,mb=mi、s=1、a=0.5,广义式退化为完整岩石的狭义式,这一“退化检验”是判断参数输入是否合理的第一道关卡。

广义Hoek-Brown准则参数体系与强度包络线 输入参数 σci 完整岩石单轴抗压强度 mi 完整岩石摩擦常数 GSI 地质强度指标 0~100 D 扰动系数 0~1 Ei:完整岩石弹性模量 岩体参数换算 mb = mi·exp[(GSI-100)/(28-14D)] 岩体Hoek-Brown常数 s = exp[(GSI-100)/(9-3D)] 岩体完整程度系数 a = 0.5 + (1/6)(e^(-GSI/15) - e^(-20/3)) 包络线曲率指数 广义Hoek-Brown准则 σ1 = σ3 + σci·(mb·σ3/σci + s)^a 压为正;σ3 为最小主应力 GSI=100, D=0 时退化为狭义式 σ1–σ3 强度包络线(示意) σ3 σ1 GSI高、D小 GSI中、D中 GSI低、D大 参数影响方向:GSI↓ 或 D↑ → mb、s 减小、a 增大 → 包络线下移、曲率增大 → 岩体强度降低 Ei 不进入强度方程,但用于估算岩体变形模量 Em = Ei·[0.02 + 1/(1+e^((60-15D-GSI)/11))]

二、各参数的物理意义辨析

σci(完整岩石单轴抗压强度)是准则的“强度标尺”,量纲为应力。它应由不含明显结构面的完整岩芯试样,按ISRM建议方法做单轴压缩试验获得,通常取5个以上有效试样的平均值。σci 决定了包络线的整体高度:σci 越大,同样应力状态下岩体可承受的 σ1 越高。工程中若缺乏试验数据,可依据点荷载试验换算或参考同类岩性经验值,但必须在报告中注明来源。

mi(完整岩石Hoek-Brown常数)反映完整岩石的摩擦特性与脆性程度,无量纲。它由一组三轴压缩试验结果,以 (σ1−σ3) 对 σ3 做非线性回归拟合得到。经验上,结晶质硬岩 mi 偏高(如花岗岩约25~33),软岩偏低(如页岩约4~10)。mi 越大,包络线越陡,岩石摩擦分量越强。注意 mi 是“完整岩石”参数,不能直接用于节理岩体,必须经 GSI、D 换算为 mb。

GSI(地质强度指标)取值0~100,是连接地质描述与力学参数的桥梁。它由岩体结构(块状、层状、碎裂等)与结构面表面条件(粗糙度、风化程度、充填物)两个维度,对照Hoek提出的GSI图表查得。GSI 不是试验量测值,而是地质观察的经验评分,因此必须由有经验的地质人员现场判定,并记录判定依据。GSI 降低意味着结构面增多、咬合力下降,直接导致 mb、s 减小。

D(扰动系数)取值0~1,描述爆破损伤与应力释放对岩体的扰动程度。D=0 表示未扰动(如TBM掘进、控制爆破效果良好),D=1 表示严重扰动(如普通爆破、边坡开挖无控制)。D 通过分母 (28−14D) 与 (9−3D) 放大 GSI 的折减效应:同样 GSI 下,D 从0增至1,mb 与 s 显著下降。实际取值应结合开挖方式、单响药量与岩体质量综合判断,并在报告中说明依据。

Ei(完整岩石弹性模量)不进入强度方程,但它是估算岩体变形模量 Em 的基准。常用经验式为 Em = Ei·[0.02 + 1/(1+e^((60−15D−GSI)/11))]。Ei 通常由单轴压缩试验的应力—应变曲线线性段斜率确定,单位与应力相同(GPa)。当 GSI 较低时,Em/Ei 比值急剧下降,这正是节理岩体“强度与刚度同步折减”的体现。

三、参数输入的合理性检验

在把参数代入方程前,建议做三项检查。第一,退化检验:令 GSI=100、D=0,确认 mb=mi、s=1、a=0.5,说明换算式实现无误。第二,量级检验:用示例参数 σci=80 MPa、mi=12、GSI=55、D=0.7,可算得 mb≈12·exp(−45/18.2)≈1.01,s≈exp(−45/6.9)≈0.0015,a≈0.5+(1/6)(e^(−3.67)−e^(−6.67))≈0.504。此时岩体单轴抗压强度 σcm=σci·s^a≈80×0.0015^0.504≈3.1 MPa,与“GSI中等、强扰动”的定性判断一致。第三,边界检验:σ3=0 时 σ1=σci·s^a,应等于岩体单轴抗压强度;σ3 增大时包络线应单调上升且斜率递减。若计算结果违反这些基本性质,应回查参数单位(σci 与 σ3 是否同为 MPa)与 GSI、D 的取值区间。

教学示例中的 σci=80 MPa、mi=12、GSI=55、D=0.7 为合理假设值,用于演示换算流程与量级判断,并非某工程实测结果。实际项目应以试验与现场地质编录为准。

理解这套参数体系的关键,是分清“哪些是试验量测的(σci、mi、Ei)、哪些是地质判定的(GSI、D)”。前者决定完整岩石的强度基准,后者决定从完整岩石到节理岩体的折减路径。只有两类参数都来源可靠、取值有据,广义Hoek-Brown准则给出的强度包络线才具备工程可信度。