3. 圆形巷道围岩塑性区经典解析计算

在矿山巷道支护设计中，塑性区半径是最基础、也最常被引用的一个量。它回答的是这样一个问题：开挖之后，围岩中究竟有多大一圈已经进入塑性屈服状态？经典解析解之所以长期被使用，是因为它把复杂的围岩-支护相互作用简化为一个可以手算的力学模型，在方案比选、参数敏感性分析和现场快速估算中依然有不可替代的价值。本节围绕轴对称圆形巷道、Mohr-Coulomb屈服准则这一经典组合，讲清公式从哪里来、什么时候能用、以及实际算的时候每一步该取什么值。

一、力学模型与基本假定

经典解建立在如下理想化条件之上：巷道为圆形，半径记为 a；围岩均质、各向同性、连续；原岩应力为静水压力状态，即水平与垂直地应力相等，记为 p₀；巷道无限长，按平面应变问题处理；围岩服从 Mohr-Coulomb 屈服准则；塑性区内的岩体视为理想弹塑性材料，屈服后强度参数不随塑性应变变化。在这些假定下，问题具有轴对称性，围岩中只有径向应力 σr 和切向应力 σθ 两个主应力，且 σθ 为最大主应力、σr 为最小主应力，剪应力为零。

需要特别强调静水压力这一条。真实地应力往往存在侧压比不等于 1 的情况，此时塑性区不再是同心圆，而呈椭圆形或更复杂的形态，经典公式只能作为量级估计。因此本节的公式适用于侧压比接近 1、围岩相对均质、且以初步设计估算为目的的场景。

二、Mohr-Coulomb 准则的应力表达

Mohr-Coulomb 准则在主应力空间中写作 τ = c + σ·tanφ，转换为最大、最小主应力表达为：

σθ = σc + k·σr
其中  k = (1 + sinφ) / (1 − sinφ)
      σc = 2c·cosφ / (1 − sinφ)

这里 c 为黏聚力，φ 为内摩擦角，k 常被称为塑性系数。k 的物理含义是围压对强度的放大倍率：φ 越大，k 越大，围岩屈服后所能承受的切向应力越高。σc 则是无围压条件下的单轴抗压强度（对符合该准则的材料而言）。这两个量把岩体强度参数直接转成了后续公式的输入。

三、塑性区半径公式的推导思路

推导分两个区域进行。塑性区内，应力必须同时满足 Mohr-Coulomb 准则和轴对称平衡微分方程：

dσr/dr + (σr − σθ)/r = 0

将准则式代入平衡方程，得到一个关于 σr 的一阶常微分方程，结合塑性区内边界条件——在洞壁 r = a 处径向应力等于支护力 pi（无支护时 pi = 0）——积分可得塑性区应力分布：

σr = (pi + σc/(k−1))·(r/a)^(k−1) − σc/(k−1)
σθ = k·σr + σc

弹性区则按厚壁圆筒解处理，其应力随半径增大逐渐趋近原岩应力 p₀。关键在于弹塑性交界面上应力必须连续：交界面处既满足弹性解，又恰好达到屈服。把弹性区在 r = R 处的应力代入屈服准则，并令其与塑性区外侧应力相等，即可解出塑性区半径 R：

R = a · [ (2p₀ + σc·(k−1) − 2pi) / ((k+1)·(pi + σc/(k−1))) ]^(1/(k−1))

当巷道无支护、即 pi = 0 时，公式简化为更常用的形式：

R = a · [ (2p₀ + σc·(k−1)) / (σc·(k+1)/(k−1)) ]^(1/(k−1))

从公式结构可以看出三个直观结论：原岩应力 p₀ 越大，R 越大；支护力 pi 越大，R 越小，这正是支护抑制塑性区扩展的力学依据；强度参数 c、φ 越大，k 与 σc 越大，R 越小。下图把这一从强度参数到塑性区半径的传递关系与计算路径集中呈现出来。

四、适用条件与常见误用

使用该公式前，应逐条核对适用条件，否则算出的数值可能偏离实际很远。第一，断面应接近圆形或已按等效圆处理，等效半径一般取外接圆与内切圆半径的平均值；第二，原岩应力应接近静水压力状态，若侧压比明显偏离 1，应改用非轴对称解或数值方法；第三，围岩应相对均质，不存在控制性软弱夹层或断层带，否则塑性区会沿弱面局部突出；第四，岩体强度参数应取现场尺度下的等效 c、φ，而不是室内完整岩块试验值，通常需要按地质强度指标或经验折减；第五，公式给出的是平面应变条件下的稳定塑性区，未考虑时间效应、流变和开挖扰动引起的强度衰减。

五、计算步骤与示例

实际计算可按以下顺序推进：确定巷道半径 a 与等效半径；根据地应力实测或估算确定 p₀；由岩体分级或试验折减确定 c、φ；计算 k 与 σc；确定对应阶段的支护力 pi；代入公式求 R；最后用 R/a 判断塑性区相对规模并复核参数合理性。

下面给出一个明确属于教学假设的示例：某圆形巷道半径 a = 2.0 m，原岩应力 p₀ = 15 MPa，岩体黏聚力 c = 1.5 MPa，内摩擦角 φ = 25°，无支护，即 pi = 0。先算强度参数：

k  = (1 + sin25°)/(1 − sin25°) = (1+0.4226)/(1−0.4226) ≈ 2.464
σc = 2×1.5×cos25°/(1−sin25°) = 3×0.9063/0.5774 ≈ 4.71 MPa

代入无支护公式：

分子 = 2×15 + 4.71×(2.464−1) = 30 + 6.90 = 36.90
分母 = 4.71×(2.464+1)/(2.464−1) = 4.71×3.464/1.464 ≈ 11.14
R = 2.0 × (36.90/11.14)^(1/1.464)
  = 2.0 × (3.312)^(0.683) ≈ 2.0 × 2.26 ≈ 4.52 m

即塑性区半径约为 4.5 m，塑性区相对半径 R/a ≈ 2.26，说明塑性区已扩展到约 2.3 倍洞径范围，属于需要重视支护的规模。若施加 pi = 0.3 MPa 的支护力，重新代入可得 R 明显减小，这一对比正是支护参数优化的直接依据。需要说明的是，上述数值为教学示例，实际工程应以现场实测地应力和经折减的岩体强度参数为准。

掌握这一经典解的意义不仅在于得到一个半径数值，更在于理解各参数对塑性区的影响方向与敏感程度。它是后续数值模拟结果判读、支护设计取值和现场监测预警的共同参照基准。