隧道开挖破坏了隧址区的应力平衡,导致围岩应力及位移场发生复杂变化,引发围岩变形过大及支护体系开裂失稳等问题,是地下工程研究的重点问题^[1-2].

针对围岩应力及位移场复杂计算问题,朱艳峰等^[3]提出将Hoek-Brown强度准则的应变软化模型应用于隧道塑性区分析中,推导了隧道塑性区半径表达式.陆晓清等^[4]针对岩体开挖导致应力重分布问题,以深埋轴对称圆形硐室为研究对象,推导了理想弹塑性解析解.贺耕夫等^[5]依据Hoek-Brown破坏准则,推导了寒区圆形隧洞理想弹塑性围岩的塑性区半径解析解.夏才初等^[6]基于三维非线性Hoek-Brown强度准则,构建了考虑应变软化特性的圆形隧道开挖后围岩非线性力学响应的求解方法.文献[7-8]基于Hoek-Brown本构模型,推导了圆形隧道的弹塑性解,探讨了其适用性.Massinas等^[9]在双极坐标的基础上,利用Mohr-Coulomb强度准则推导了圆形隧道围岩应力和塑性区的解析解.Exadaktylos等^[10-11]推导了半圆形隧道的应力及位移计算表达式.李培楠等^[12]和施有志等^[13]分析了工程中常用的单心圆仰拱马蹄形隧道的解析解.程长清等^[14]根据水工隧洞运行期的受力情况,运用Mohr-Coulomb与Hoek-Brown破坏准则,对无限大均质体中圆形水工隧洞在有无衬砌两种情况下的最小覆盖层厚度进行了弹塑性力学分析.

综上所述,针对服从Mohr-Coulomb准则时围岩的应力及位移解研究较多,对Hoek-Brown准则下应力及位移的求解研究较为有限,但在高围压条件下岩体节理裂隙发育,此时Hoek-Brown强度准则更合适.基于此,本文利用非线性Hoek-Brown强度准则对深埋隧道开挖产生的弹塑性区进行求解,探讨塑性区影响半径的影响因素,并对特征参数进行敏感性分析,以期为隧道开挖考虑非强度准则的应力及位移求解提供理论基础.

构建无限大均质弹性体内圆形隧道开挖问题的简化计算模型,如图1所示.

图1 地下隧道开挖计算分析模型
Fig.1 Calculation and analysis model of underground tunnel excavation

图中:σ_θ、σ_r分别为环向应力及径向应力; σ₀为作用在岩土体上初始静置土压力; p_i为隧道开挖边界上支护压力; r₀为隧道开挖的半径; r_p为隧道开挖后围岩的塑性区半径.

Hoek-Brown强度准则的屈服条件为

σ_θ-σ_r=σ_c(m(σ_r)/(σ_c)+S)^a　(1)

式中:系数m、s及a为表征岩石的基本特征的半经验参数,与参数GSI关系为

m=m_iexp((GSI-100)/(28-14D))　(2)

s=exp((GSI-100)/(9-3D))　(3)

a=1/2+1/6[exp(-(GSI)/(15))-exp(-(20)/3)]　(4)

式(3)和(4)中D是岩石应力松弛程度或破坏程的参数,取值在0～1之间; m_i取值分别为5,10,13,16和24.

轴对称问题中,岩体平衡微分方程为

(dσ_r)/(dr)+k(σ_r-σ_θ)/r=0　(5)

式中,k=1,2,分别为洞身和掌子面附近区域,即图1中的(a)与(b).

边界条件为:

为验证本文理论的有效性与可靠性,取文献[15]的实验参数进行分析.具体为:r₀=5 m,D=0.参数a、m_b和S的值则根据式(2)～(4)求得,结果为:a=0.011,m_b=5.128,S=0.037,其他参数如表1所示.

表1 计算参数取值
Tab.1 Calculation parameter values

求解参数k为不同参数的应力值,绘制应力与塑性区半径的关系曲线,如图2及图3所示.

图2 应力与rp/r0的关系图
Fig.2 The relationship between stress and rp/r0

图3 应力与开挖半径的关系
Fig.3 The relationship between stress and excavation radius

图2及图3表明,应力随着r_p/r₀的增加呈先增后减的趋势,环向应力远小于径向应力,开挖半径和r_p/r₀对环向应力的影响较径向应力显著.

开挖半径对塑性区半径的影响如图4所示,不同m_i值时,r_p/r₀与隧道开挖半径的关系如图5所示.

图4表明,塑性区半径随着开挖半径的增大而增大.图5显示,当m_i≤20时,r_p/r₀随着开挖半径的增大非线性增大; 当m_i=25,r_p/r₀随着开挖半径的增大呈现先增后减趋势,r_p/r₀与GSI的关系如图6所示.

图4 塑性区半径与隧道开挖半径的关系
Fig.4 The relationship between the radius of the plastic zone and the radius of tunnel excavation

图5 rp/r0与隧道开挖半径的关系图
Fig.5 The relationship between rp/r0 and the radius of the tunnel excavation

图6 rp/r0与GSI的关系图
Fig.6 The relationship between u/r0 and GSI

图6表明,不同m_i值时,r_p/r₀随GSI增大而非线性增大,其敏感性逐渐减弱.

利用本文推导公式,绘制u/r₀随着剪胀和开挖半径的分布规律如图7所示.

图7 u/r0与剪胀系数的关系图
Fig.7 The relationship between u/r0 and the dilatancy coefficient

由图7可知,剪胀和开挖半径对u/r₀的影响中,开挖半径对位移的影响更加显著.

(1)本文基于非线性Hoek-Brown强度准则,求解了隧道开挖后塑性区的应力和位移解析解;

(2)利用本文推导公式,对影响应力和位移解的多个因素进行了特征参数敏感性分析,探讨了参数间的耦合影响及作用机制;

(3)塑性区半径随开挖半径的增大而增大;

(4)应力随塑性半径的增加呈现先增后减的演化趋势,同时环向应力远小于径向应力.