Deterioration Characteristics of Tunnel Support Structures under Surrounding Rock Creep
-
摘要:
为探究蠕变效应下支护结构的长期劣化特性,针对隧道支护结构体系,建立锚杆破断、钢拱架屈服和混凝土塑性损伤力学模型,采用数值算例验证支护结构劣化力学模型的有效性,并探究竖向应力为主、静水压力和水平应力为主条件下锚杆破断、钢拱架屈服和衬砌损伤的劣化特性. 研究表明:破断首先发生在隧道边墙中部位置的锚杆,再沿隧道环向向两侧发展;钢拱架轴力呈先加速增长、后缓慢发展、最后急速降低的特征,轴力快速降低的同时伴随着弯矩的剧烈变化,部分测点弯矩出现由负变正的特征;受压损伤破坏区主要分布在隧道边墙和墙脚位置,受拉损伤首先出现在二次衬砌边墙中部的表面;侧压力系数越大,锚杆破断、钢拱架屈服、衬砌形成贯通的受压破坏区、受拉损伤达到最大值的时间越早.
Abstract:To investigate the long-term deterioration characteristics of tunnel support structures under creep effect, mechanical models of anchor bolt fracture, steel arch frame yielding, and concrete plastic damage were established for the tunnel support structure system. Numerical examples were used to verify the validity of the mechanical models for support structure deterioration. The deterioration characteristics of anchor bolt fracture, steel arch frame yielding, and lining damage were explored under conditions dominated by vertical stress, hydrostatic pressure, and horizontal stress. The results show that the fracture first occurs at the anchor bolt at mid-height of the tunnel sidewall and then develops circumferentially towards both sides. The axial force of the steel arch frame first increases rapidly, then develops slowly, and finally decreases significantly. The rapid decrease of axial force is accompanied by drastic changes in the bending moment, with some measuring points appearing a change of bending moment from negative to positive. The compressive damage zones are mainly distributed at the sidewall and wall foot positions of the tunnel, while tensile damage first appears on the surface of the secondary lining at the mid-height of the sidewall. As the lateral pressure coefficient increases, the anchor bolt fracture, the steel arch frame yielding, formation of a continuous compressive damage zone in the lining, and the maximum tensile damage appear earlier.
-
Key words:
- creep /
- anchor bolt fracture /
- steel arch frame yielding /
- lining damage /
- deterioration characteristic
-
矿山法隧道多采用复合式衬砌作为支护结构的主体,构成“围岩-初期支护-二次衬砌”共同作用的复合承载结构[1-3]. 围岩条件较好时,初期支护承担主要围岩荷载,二次衬砌主要作为安全储备[4, 5];围岩条件较差时,二次衬砌承担的荷载比例大为增加[6, 7]. 在复合承载体系中,初期支护主要保证施工期隧道围岩的稳定性,二次衬砌主要保证运营期隧道结构的长期安全性. 《公路隧道设计细则》[8]规定隧道衬砌须按照满足100年正常使用的永久性结构设计,《公路工程技术标准》[8]根据公路等级不同将隧道的设计使用年限规定为50~100年不等,《铁路隧道设计规范》[5]明确提出隧道衬砌的设计使用年限为100年,这些规范对隧道支护结构的耐久性提出了严格的要求. 然而,软岩隧道岩体中存在大量节理裂隙和渗水通道,由系统锚杆、钢拱架和喷射混凝土组成的支护体系长期处于恶劣的工作状态,易造成锚杆破断、钢拱架屈服、混凝土性能退化[9-12];同时,围岩在高地应力条件下产生明显的蠕变效应,导致围岩荷载随着使用年限的增加而不断增长[13-19]. 因此,在围岩蠕变作用下,支护结构发生腐蚀劣化,围岩荷载逐步由初期支护向二次衬砌转移,作用于二次衬砌的形变压力不断增长,甚至引起二次衬砌的开裂、破坏和失稳,致使混凝土开裂成为二次衬砌的常见病害.
针对隧道支护结构的长期稳定性,学者们展开了大量的研究工作. 戴永浩等[20]通过将钢支撑的力学参数等效到二次衬砌中,采用弹塑性模型来评估大坂膨胀性泥岩引水隧洞长期稳定性,指出围岩蠕变是引起隧洞底部隆起的主要原因之一. 李建军等[21]采用Burgers蠕变模型,通过二次衬砌应力、轴力、弯矩和安全系数等指标,分析了蠕变效应下二次衬砌力学性能的演变特征. 徐国文等[22]构建了蠕变作用下壳-接头的二次衬砌开裂模型,评估了二次衬砌裂损程度、开裂位置和裂纹数量对结构长期安全性的影响. 王迎超等[23]采用广义Kelvin模型,研究了蠕变荷载下支护时机、支护刚度对二次衬砌承载性能的影响. 师亚龙等[24]研究了木寨岭隧道在不同流变周期内支护结构受力随时间变化的规律,指出围岩流变性对衬砌结构长期稳定性影响极大,衬砌结构在数年后发生压溃开裂的风险较大. 钱文喜等[25]拟合得到泥质粉砂岩的Cvisc模型参数,重点研究了隧道埋深对运营隧道衬砌结构稳定性的影响. 上述的研究工作都考虑了蠕变荷载对衬砌长期稳定性与安全性的影响,但忽略了支护体系中锚杆、钢拱架可能发生的锈蚀断裂行为,也忽视了衬砌在达到峰值强度后承载性能的降低.
本文以宗思隧道为工程背景,针对隧道支护结构体系的长期劣化行为,建立锚杆破断、钢拱架屈服、混凝土塑性损伤力学模型,分析蠕变效应下支护结构长期劣化特性,为蠕变围岩隧道的设计提供参考,为隧道长期安全运营和维护提供依据.
1. 围岩蠕变与支护结构劣化力学模型
1.1 围岩蠕变力学模型
隧道开挖后,岩体的应力重新分布导致隧道围岩发生变形. 如果忽略岩体的时效特征,则开挖完成后隧道围岩完成全部变形. 但是,大量的现场监控量测数据表明隧道开挖后围岩的变形随时间推移不断增长,表现出明显的时间效应,蠕变作用产生的附加荷载使隧道变形不断增长,导致支护结构受力增大,进而导致结构发生劣化,甚至会引起隧道的失稳坍塌. 因此,在研究隧道开挖后支护体系的力学行为时,需要考虑围岩的蠕变特征,才能获得支护结构的真实力学响应[26]. 因同时考虑隧道开挖后围岩的塑性和蠕变行为,围岩本构模型采用Burgers粘塑性蠕变模型(CVISC模型),由Kelvin体、Maxwell体和塑性体串联而成,见图1. 图中:σ为应力,EM、ηM分别为Maxwell体的弹性模量、黏滞系数;EK、ηK分别为Kelvin体的弹性模量、黏滞系数. 模型中塑性体采用Mohr-Coulomb屈服准则,其屈服函数f如式(1)所示.
f={σ1−σ3Nφ+2c√Nφ,剪切破坏,σ3−σt,拉伸破坏, (1) 式中:σ1、σ3分别为第一、第三主应力;Nφ=(1 + sinφ)/(1-sinφ),φ为内摩擦角;c为内聚力;σt为抗拉强度.
当主应力满足屈服准则时,围岩既发生塑性变形也发生蠕变变形. 当主应力未达到塑性屈服条件时,模型退化为黏弹性的Burgers模型,其蠕变方程为
ε=[1EM+1ηMt+1EK(1−e−EKηKt)]σ, (2) 式中,ε为应变,t为时间.
1.2 锚杆破断力学模型
在隧道支护体系中,锚杆作为轴向受力构件,主要承受拉力,在达到一定的伸长率后,锚杆发生断裂. 对于端固式锚杆,杆体断裂后轴力下降为0,锚杆对围岩变形的控制作用完全丧失;对于全长粘结式锚杆,断裂处轴力为0,其余部位发生卸荷回弹,锚杆对围岩形变的约束效应大大减弱. 在包括FLAC3D(fast lagrangian analyses of continua in three-dimensions)[27]、Abaqus[28]、Ansys[29]等数值模拟软件中,锚杆多采用弹性或弹塑性本构模型,弹性模型的轴力随应变的增加而线性增大,不存在上限;弹塑性模型的轴力存在最大值,但进入塑性状态后,无论应变多大,锚杆轴力保持不变(图2(a)),即锚杆不会发生断裂. 因此,弹性模型和弹塑性模型都不能真实体现锚杆支护破断后的力学效应. 为分析锚杆在围岩发生大变形条件下的工作性能,需将锚杆的破断效应(图2(b))引入,其本构关系为
F={EAu,F⩽Ft,u<umax,Ft,F>Ft,u<umax,0,u⩾umax, (3) 式中,F为锚杆单元轴力,u为锚杆单元变形量,E为锚杆弹性模量,A为锚杆横截面积,EA为锚杆轴向刚度,Ft为锚杆抗拉强度,umax为锚杆极限变形量.
锚杆破断力学模型通过Fish语言编程实现. 当锚杆单元轴力F达到抗拉强度Ft且变形量u大于极限变形量umax时,将锚杆单元的横截面积降低为初始面积的万分之一,即可实现锚杆破断轴力降低的力学效果. 验证锚杆破断行为的数值模型见图3(a),采用长度为3 m全长黏结式锚杆,将锚杆划分为10个单元,编号为C01~C10,每个单元长度为300 mm. 锚杆物理力学参数[30]见表1,锚杆破断荷载约为127 kN. 将节点N02固定,在节点N01施加速度1 × 10−5 m/步的位移边界条件,计算步设为
50000 步. 锚杆单元的轴力随计算步的变化曲线见图3(b). 可以看出,随着计算步的增加,锚杆单元的轴力快速增长至127 KN,之后轴力保持不变;当计算步达到24800 步时,单元C01发生破断,轴力快速降低至0,其余单元的轴力随后也逐步降低至0. 锚杆单元轴力与长度的关系见图3(c). 当单元长度达到约348 mm时,锚杆延伸率为16%,锚杆发生破断,轴力降低为0,而弹塑性模型的轴力始终保持在127 kN. 由此可见,破断模型能够实现锚杆真实的支护力学行为.表 1 锚杆物理力学参数Table 1. Physical and mechanical parameters of anchor bolt弹性模量E/GPa 抗拉强度Ft/MPa 延伸率/% 直径d/mm 200 335 16 22 1.3 钢拱架屈服力学模型
钢拱架既能承受轴向拉压荷载,也能发挥抵抗弯曲变形的作用. 当轴力或弯矩超过钢拱架的极限承载能力时,钢拱架必然发生屈服,轴力和弯矩迅速显著降低. 钢拱架的力学行为常采用线弹性特征的梁单元来模拟,其轴力和弯矩会随变形无限增大,夸大钢拱架的支护效果. 因此,引入参数极限轴力Nu和极限弯矩Mu来表征钢拱架的屈服力学行为,通过Fish语言在梁单元节点之间增加拉压屈服弹簧和扭转屈服弹簧以实现钢拱架的屈服. 当钢拱架轴力达到极限轴力时,将拉压屈服弹簧和扭转弹簧的最大承载力都降低并接近0;当弯矩达到极限弯矩时,将扭转弹簧和拉压屈服弹簧的最大承载力降低至趋于0,实现钢拱架承载能力达到峰值后降低的力学效果.
以I20a型工字钢为例,钢拱架屈服的数值模型见图(4a)、(4b),钢拱架被划分为10个梁单元,编号为B01~B10,钢拱架的物理力学参数[30]见表2. 模型共有20个节点,在梁单元连接处存在2个坐标相同的节点,例如N10和N11的坐标相同. 相同坐标节点之间通过屈服弹簧相连接. 验证钢拱架轴向屈服的特性时,将节点N20固定,在节点N01处施加速度为1 × 10−7 m/步的水平位移边界条件,计算
50000 步;验证钢拱架的极限弯矩时,固定节点N01和N20的竖向位移,在节点N10处向下施加速度为1 × 10−8 m/步的位移边界条件,计算15000 步. 图(4c)为钢拱架轴力随时间步变化曲线,初始模型的轴力随时间步线性增长,最大值达到1460 kN,改进模型的轴力达到设定的849 kN后迅速降低至0. 图(4d)为钢拱架弯矩随时间步变化曲线,初始模型的弯矩达到100 kN•m,而修正模型的最大弯矩在达到最大值71 kN•m后迅速跌落至0. 由此可见,改进模型能反映钢拱架的屈服力学特性.表 2 钢拱架物理力学参数Table 2. Physical and parameters of steel arch frame钢拱架型号 弹性模量E/GPa 泊松比v 极限轴力Nu/kN 极限弯矩Mu/kN•m I20a 200 0.3 849 71 1.4 混凝土塑性损伤力学模型
混凝土是一种组分变化显著的人工复合材料,其材料配比和浇筑养护环境的差异,导致力学性能复杂多变. 混凝土在单轴压缩和单轴拉伸条件下,表现出高度的非线性特征. 在隧道施工力学行为中,衬砌多采用线弹性本构关系,缺乏专门的本构模型,忽略了衬砌开裂后力学性能的降低,无法准确获得衬砌的支护力学响应特征. 因此,针对受压和受拉状态下不同的破坏方式,以Mohr-Coulomb破坏准则为基础,建立混凝土的损伤力学模型.
混凝土的受压塑性损伤可以分为峰前和峰后2个阶段,见图5(a). 在峰前损伤区,受压损伤Ds以割线模量的降低来定义,则有
Ds=1−E1E0 = 1−σE0ε, (4) 式中:E0为混凝土初始弹性模量,E1为混凝土割线弹性模量,σ为应力.
峰值应力σp对应的损伤值Dp为
Dp=1−σpE0εp. (5) 在峰后损伤区,损伤以应力的降低来定义,则有
Ds=Dp+(1−Dp)(1−σσp) = 1−(1−Dp)σσp. (6) 混凝土在受拉状态下,假设应力达到抗拉强度σt时进入损伤状态(图5(b)),则受拉损伤Dt为
Dt=1−σσt. (7) 为实现塑性损伤模型的二次开发,需要将损伤关系转化为三维差分形式. 在主应力空间,三维差分形式的应力增量表达式为
{Δσ1 = α1Δεe1+α2(Δεe2+Δεe3)Δσ2 = α1Δεe2+α2(Δεe1+Δεe3)Δσ3 = α1Δεe3+α2(Δεe1+Δεe2) (8) {α1 = (1−D)(K+4G/4G33)α2 = (1−D)(K−2G/2G33) (9) 式中:Δσ1、Δσ2和Δσ3分别表示第一、第二和第三主应力的增量,Δεe 1、Δεe 2和Δεe 3分别表示第一、第二和第三主应变的弹性增量部分,α1和α2均为材料常数;K为体积模量,K=E/[3(1−2v)];G为剪切模量,G=E/[2(1 + v)].
当应力水平达到峰值强度的约30%~80%时,混凝土由弹性状态过渡到塑性状态[31]. 本研究中,取其为0.7,则受压状态下屈服准则为
fs={σ1−σ3Nϕ+2c√Nϕ(0.7+0.3DsDp),Ds⩽Dp,σ1−σ3Nϕ+2c√Nϕ(1−Ds)1−Dp,Ds>Dp. (10) 式中:ϕ为内摩擦角,Nϕ=(1 + sinϕ)/(1-sinϕ).
受拉状态下屈服准则为
ft=σ3−(1−Dt)σt. (11) 求解过程中,根据主应力空间中的判断函数h(σ1,σ3)=0来判断单元进入塑性状态后的屈服类型,其表达式为
h=σ3−σt+(√1+N2ϕ+Nϕ)(σ1−σtNϕ+2c√Nϕ). (12) 单元的屈服类型根据式(13)进行判断.
屈服类型={受压屈服,fs≤0且h≤0受拉屈服,ft≤0且h>0 (13) 当单元发生受压损伤时,主应力为
{σN1=σI1−λs(α1−α2Nψ),σN2=σI2−λsα2(1−Nψ),σN3=σI3−λs(−α1Nψ+α2), (14) λs=fs(σI1,σI3)(α1−α2Nψ)+(−α1Nψ+α2)Nϕ, (15) 式中:σN 1、σN 2和σN 3分别表示更新的第一、第二和第三主应力,σI 1、σI 2和σI 3分别表示以弹性状态计算的第一、第二和第三主应力,λs是与受压屈服相关的正值常数,Ψ为剪胀角,NΨ=(1 + sinΨ)/(1−sinΨ).
当单元发生受拉损伤时,主应力为
{σN1=σI1−λtα2,σN2=σI2−λtα2,σN3=σI3−λtα1, (16) λt = σI3−σtα1. (17) 式中,λt是与受拉屈服相关的正值常数.
通过建立塑性应变与损伤变量之间的函数映射关系,即可得到混凝土的损伤力学行为,这种映射关系可通过对材料的塑性损伤行为反分析来确定. 受压状态下塑性应变采用塑性剪切应变κs,受拉状态下塑性应变采用塑性拉伸应变κt,则对应的函数关系F如式(18)所示.
{Ds=F(κs),Dt=F(κt). (18) 塑性剪切应变κs为塑性剪切应变增量Δκs之和,如式(19)所示.
{κs=∑Δκs=√(Δεps1−Δεpsm)2+Δεpsm2+(Δεps3−Δεpsm)22,Δεpsm=Δεps1+Δεps33,Δεps1=λs,Δεps3=−λsNψ. (19) 塑性拉伸应变κt为塑性拉伸应变增量Δκt之和,如式(20)所示.
{κt=∑Δκt=∑|Δεpt3|,Δεpt3=λt. (20) 根据上述的损伤力学模型,推导出的本构方程三维差分形式,即可在C + + 语言环境下完成对模型进行程序的二次开发,编译成动态链库文件,通过命令config cppudm和model load实现本构模型的注册及调用.
以混凝土结构设计规范GB50010—2002[32]中的C25混凝土曲线为参照,用单轴压缩试验和单轴拉伸试验对提出的损伤模型进行验证. 选用的混凝土基本力学参数见表3,损伤参数见表4. 建立边长为150 mm的立方体数值模型. 固定试件的底部,从顶部向下压缩,记录顶部中心的竖向位移和顶面反力.
表 3 混凝土力学参数Table 3. Mechanical parameters of concrete强度
等级弹性模量E /GPa 泊松比ν 内聚力c/MPa 内摩擦角φ/(°) 抗压强度σc /MP 抗拉强度σt /MPa C25 28.0 0.2 3.0 50 16.7 1.78 C35 31.5 0.2 3.7 55 23.4 2.20 表 4 混凝土弹塑性损伤参数Table 4. Elasto-plastic damage parameters of concreteC25混凝土 C35混凝土 κs/
( × 10−3)Ds κt/
( × 10−3)Dt κs/
( × 10−3)Ds κt/
( × 10−3)Dt 0 0 0 0 0 0 0 0 0.06 0.01 0.42 0.55 0.07 0.01 0.45 0.55 0.20 0.02 1.02 0.75 0.23 0.02 1.13 0.75 0.45 0.03 3.60 0.99 0.56 0.04 3.92 0.99 1.20 0.11 / / 1.40 0.12 / / 3.00 0.40 / / 2.90 0.39 / / 4.50 0.70 / / 4.10 0.67 / / 7.00 0.99 / / 6.00 0.99 / / 单轴压缩下,模拟应力-应变曲线在峰值应力的70%之前呈现弹性应力状态,见图5(c). 当荷载超过峰值的70%后,混凝土进入剪切破坏状态,曲线斜率逐渐减小,直至应力达到峰值. 在峰值应力后,曲线表现出明显的应变软化行为. 此时,数值模拟得到的峰值应力非常接近标准值. 在单轴拉伸下,模拟应力-应变曲线在峰值应力前也表现出明显的弹性,峰后阶段表现出应变软化行为,如图5(d)所示. 曲线与规范中的曲线吻合较好,开发的损伤模型能够有效反映混凝土在单轴受拉时的力学行为. 因此,开发的混凝土弹塑性损伤模型能较好地描述混凝土在受压状态、受拉状态下的非线性力学行为.
2. 工程背景及数值模型
2.1 工程背景
宗思隧道是丽江至香格里拉铁路的控制性工程之一. 隧道全长2 205 m,大围岩1 250 m,V级围岩955 m,最大埋深190 m,隧道纵、横断面见图6. 隧道位于青藏高原东南缘的滇西地震带,区内地质构造复杂、构造运动强烈、地震活动频繁. 受褶皱、断裂和地震的影响,区内岩层扭曲变形严重,节理裂隙发育、层间结合性差,岩体结构松散破碎,手掰即碎. 岩体的强度低,属极软质岩石,浸水后立即出现软化、泥化现象,具有塑性蠕变形特性. 施工初期支护开裂剥落、侵限情况严重,隧道收敛变形量大,水平收敛最大为738.9 mm,收敛速率达最大为151.8 mm/d,致使隧道大范围换拱,换拱比例极高,隧道大变形已成为影响施工安全、制约工程进度、降低工程经济效益的关键问题. 受复杂的地质条件及严重的大变形困扰,可以预见在隧道运营期的必将面对支护结构劣化问题.
2.2 数值模型的建立
以宗思隧道为研究背景建立数值模型,模型宽80 m,高80 m,走向长10 m,见图7(a). 锚杆长度为3.0 m,直径为22 mm,采用cable单元;钢拱架型号为I20a,间距为1.0 m,采用beam单元;初期支护采用C25混凝土,二次衬砌采用C35混凝土,围岩、初期支护和二次衬砌采用实体单元,见图7(b). 模型共有
38720 实体单元、800个beam单元、1900 个cable单元. 锚杆、钢拱架和衬砌的物理力学参数分别见表1 ~ 4. 结合炭质页岩力学特性试验成果[33-34],综合宗思隧道施工过程中的围岩位移表现,选取炭质页岩的蠕变力学参数取值见表5. 模型的前、后、左、右和底面采用位移边界条件,控制其法向位移为0;模型的顶面采用应力边界条件,以隧道埋深100 m处为研究对象,隧道洞身处的竖向主应力约为2.7 MPa,因此,模型顶部施加大小为1.485 MPa的面力,设定围岩蠕变时间为100 a. 宗思隧道现场未开展地应力测试,地应力分布情况不明,因此考虑侧压力系数λ的取值0.5、1.0和1.5共3种工况,分别代表竖向应力为主、静水压力和水平应力为主3种条件. 假设隧道一次开挖完成,开挖后及时施作初期支护,当初期支护位移稳定后施作二次衬砌,再进行蠕变计算. 在数值模型中,对顺时针方向编号的锚杆轴力进行监测,同时对B1 ~ B4点位的钢拱架轴力与弯矩进行监测,监测点布置见图7(c).表 5 围岩蠕变力学参数Table 5. Creep mechanical parameters of surrounding rockK/
GPaEM/
GPaEK/
GPaηM/
GPa.aηK/
(GPa.a)c/
MPaφ/(°) σt/
MPa2 2.5 4 3 1.2 7 26 1 3. 围岩蠕变作用下隧道支护结构劣化特征
3.1 锚杆破断行为分析
3.1.1 锚杆破断时间演化特征
侧压力系数λ为0.5的轴力监测曲线见图8(a). 当蠕变时间达到54 a时,锚杆C3首先发生破断,其他锚杆随时间的增加而逐渐破断;当达到85 a时,每根锚杆都发生断裂破坏. 锚杆破断顺序为C3→ C17→C18→C2→C16→C4→C1→C19→C15→C5→C13→C14→C12→C6→C7→C8→C10→C9→C11. 侧压力系数λ为1.0的轴力监测曲线见图8(b). 首先发生破断的锚杆为C17,破断时间为26 a;当达到68 a时,所有锚杆均发生破断,最后发生破断的为锚杆C11. 锚杆破断顺序为C3→C17→C16→C4→C2→C18→C1→C19→C5→C15→C6→C14→C7→C13→C8→C12→C10→C9→C11. 侧压力系数λ为1.5的轴力监测曲线见图8(c). 蠕变时间达到14 a时,锚杆C3首先发生破断;当达到47 a时,所有锚杆都发生破断. 锚杆的破断顺序为C3→C17→C4 →C16→C2→C18→C5→C1→C15→C6→C19→C14→C7→C8→C13→C9→C12→C10→C11.
3.1.2 锚杆破断空间分布特征
蠕变作用下锚杆破断部位的空间分布情况见图9. 图中可以看出,无论地应力以竖向应力为主,或者以水平应力为主,还是静水压力条件,每根锚杆都发生破断,隧道两侧的锚杆几乎在全长范围发生破断,拱部锚杆的破断长度有所减小、破断部位分布在锚杆中部. 此外,在侧压力系数λ增大的同时,位于拱部锚杆的破断长度逐步增加. 锚杆破断时间分布见图10. 结果表明,3种侧压力条件下,破断首先出现在隧道边墙中部的锚杆C3和锚杆C17,然后沿隧道环向向两侧发展;靠近拱部和位于最下侧的锚杆破断时间较晚. 此外,侧压力系数越大,锚杆发生破断时间越早.
3.2 钢拱架屈服行为分析
3.2.1 钢拱架轴力演化特征
钢拱架轴力监测曲线见图11,弯矩监测曲线见图12. 当侧压力系数λ为0.5时,轴力在初期的2年内迅速上升,而后呈缓慢增长趋势,轴力最大达到约550 kN;当蠕变达到63 a时,拱顶B1处钢拱架轴力为360 kN、弯矩为71 kN•m,因达到极限弯矩发生屈服,此后轴力迅速降低;当达到75a时,测点B2、B3和B4的轴力快速降低至200 kN以内. 当侧压力系数λ为1.0时,轴力同样在初期的2年内快速上升,之后进入一个稳定发展时期,轴力变化缓慢;当达到26 a时,拱顶B1处钢拱架轴力为410 kN、弯矩为71 kN•m,达到极限弯矩而未达到极限轴力,钢拱架屈服、轴力迅速降低;当达到62 a后,测点B2的轴力快速降低至75 kN以内,测点B3和B4的轴力短暂上扬后也快速降低至150 kN以内,其中测点B4的轴力最大达到600 kN. 当侧压力系数λ为1.5时,钢拱架轴力整体变化规律与前面相似,蠕变时间约为17 a时,拱顶B1处钢拱架轴力因达到极限弯矩而屈服,轴力迅速降低;B3测点的轴力最大接近750 kN;当达到40 a时,测点B2、B3和B4的轴力开始下降,其中测点B3的轴力保持在440 kN左右,而B2和B4测点的轴力降低至225 kN之内. 因此,钢拱架的轴力呈现先快速增长、再缓慢发展、后期快速降低的特征;侧压力系数越大,蠕变后期轴力降低的时间越早.
3.2.2 钢拱架弯矩演化特征
从弯矩上看,侧压力系数λ为0.5时,B1测点弯矩在达到极限承载力之前近似线性增长,之后降低为0;B2测点始终为正弯矩,B3测点始终为负弯矩,B3测点的弯矩在80 a左右由负弯矩变为正弯矩;在75 a左右,弯矩发生剧烈变化,轴力在相应时间同样发生快速变化. 当侧压力系数λ为1.0时,拱顶B1测点轴力在3年内快速上升,此后到26 a之前平稳上升,直至达到极限弯矩屈服;B2测点在前95 a全部为正弯矩,在最后5 a时处于正负间波动状态;B3和B4测点在前65 a为负弯矩,此后逐步转变为正弯矩,B4测点的弯矩绝对值始终大于B3测点. 当侧压力系数λ为1.5时,拱顶B1测点同样因达到极限弯矩而在17 a时发生屈服,B2、B3测点弯矩在40 a左右开始逐步由负弯矩转变为正弯矩,且B3测点的弯矩绝对值始终大于B2测点的弯矩绝对值;B4测点的弯矩始终为正值,前40 a在20 kN•m左右波动,后面则超过30 kN•m. 整体上看,拱顶测点B1都因达到极限弯矩而屈服,且侧压力系数越大,钢拱架发生屈服的时间越早;在蠕变后期,B2、B3和B4测点的钢拱架轴力都降低,但部分测点弯矩则出现由负变正的特征.
3.3 衬砌损伤行为分析
3.3.1 衬砌受压损伤演化特征
蠕变作用下衬砌受压损伤分布见图13,其中Contour of Property damagec表示受压损伤值. 在侧压力系数λ为0.5的条件下,蠕变时间为10 a时,损伤区分布于边墙位置的初期支护、墙脚位置的初期支护与二次衬砌中,最大受压损伤值为0.02;此后直至70 a,边墙位置和墙脚位置的损伤都快速发展,但边墙位置的损伤发展比墙脚位置更显著,70 a时的最大受压损伤值为0.46;当达到80 a时,损伤达0.99的区域在衬砌厚度方向贯通,形成贯通的受压破坏区,其中,边墙位置受压破坏区贯通初期支护和二次衬砌,墙脚位置的受压破坏区仅贯通二次衬砌;当达到90 a时,墙脚处初期支护形成贯通的受压破坏区;达到100 a时,形成了边墙至墙脚范围的受压破坏区,拱部区域损伤开始逐步发展. 在侧压力系数λ为1.0的
条件下,受压损伤区的发展规律与侧压力系数λ为0.5时的基本一致,但形成贯通受压破坏区的时间提前. 当蠕变时间达到70 a时,在边墙和墙脚位置形成了贯通的受压破坏区;当达到80 a时,在拱部的二衬形成损伤达到0.99的受压破坏区,此后受压破坏区沿衬砌厚度方向和环向继续发展. 在侧压力系数λ为1.5的条件下,蠕变时间为10 a时,在拱部、边墙和墙脚出现受压损伤,最大损伤值为0.07;此后,边墙部位的损伤区域快速发展;当达到50 a时,在边墙和墙脚出现贯通的受压破坏区;蠕变时间达到60 a时,在拱部出现对称的损伤区;时间达到90 a时,拱部对称的受压破坏区已经贯穿二次衬砌,逐步向初期支护扩展.
通过上述分析可知,侧压力系数越大,衬砌中形成贯通的受压破坏区的时间越早;侧压力系数对隧道衬砌受压损伤的空间分布影响不明显,损伤破坏区主要分布在隧道边墙和墙脚位置;此外,拱部也会形成受压损伤区,且侧压力系数越大,拱部的受压损伤越严重.
3.3.2 衬砌受拉损伤演化特征
蠕变作用下衬砌受拉伤分布情况见图14,其中Contour of Property damaget表示受拉损伤值. 侧压力系数λ为0.5的条件下,蠕变时间为10 a时,在仰拱中心和墙脚位置的二次衬砌出现受拉损伤区,最大受拉损伤值为0.08;当达到20 a时,边墙位置的二次衬砌在表面出现受拉损伤区;随着时间的发展,边墙位置的二次衬砌受拉损伤快速发展;当达到70 a时,边墙位置二次衬砌表面的最大受拉损伤达到0.58;当达到80 a时,边墙位置二次受拉损伤达到0.99;达到90 a时,在拱部的衬砌出现对称分布的受拉损伤区,此后损伤区继续扩展. 在侧压力系数λ为1.0的条件下,衬砌受拉损伤演化特征与前面基本一致,但达到70 a时,最大受拉损伤已经达到0.99,同时在拱部的衬砌出现对称分布的受拉损伤区. 在侧压力系数λ为1.5的条件下,蠕变时间达到50 a时,衬砌最大受拉损伤已经达到0.99. 整体上看,隧道衬砌受拉损伤区域和损伤严重程度小于受压损伤,即以受压损伤破坏为主.
综上所述,随着侧压力系数增大,衬砌受拉损伤达到最大值0.99的时间提前,但对受拉损伤区的空间分布影响不明显;受拉损伤首先出现在边墙中部位置的二次衬砌表面,且发展速度快.
3.4 验证与讨论
在宗思隧道现场,对DK124 + 850 ~ 920段落二次衬砌裂纹进行了监测,裂纹分布见图15. 图中可知,二次衬砌中裂纹主要沿隧道轴线方向分布,少量裂纹沿隧道环向或与轴线斜交方向分布. 沿轴线方向分布的裂纹又主要分布于拱墙分界线附近,或者靠近隧道矮边墙. 沿隧道轴线方向裂纹的分布特征与衬砌损伤分布特征相似,这在一定程度上验证了隧道支护结构劣化力学模型的正确性.
根据蠕变作用下隧道支护体系中锚杆、钢拱架和衬砌的劣化特性分析可知,侧压力系数对锚杆破断、钢拱架屈服和衬砌损伤的空间分布特征影响不明显,但对它们随时间的发展变化影响显著. 侧压力系数越大,锚杆发生破断、钢拱架屈服和衬砌形成贯通的损伤破坏区的时间都会越早,这是因为3种侧压力系数下,竖向应力的大小始终保持不变,侧压力系数增大的同时,水平应力随之增大. 在蠕变条件下,更大的水平应力产生更显著的蠕变荷载,导致作用在支护结构上的作用力增加,相同时间条件下支护结构必然发生更大的变形,致使劣化效应发生的时间提前.
4. 结 论
1) 建立了锚杆破断、钢拱架屈服和混凝土塑性损伤力学模型,实现了支护结构劣化力学模型的二次开发,采用数值算例验证了3种劣化力学模型的有效性.
2) 无论侧压力系数的大小,隧道边墙中部位置的锚杆都先发生破断,然后沿隧道环向向两侧发展,靠近拱部和最下侧锚杆的破断时间较晚;侧压力系数越大,锚杆发生破断时间越早.
3) 无论侧压力系数大小,拱顶处钢拱架均因达到极限弯矩而屈服;钢拱架轴力呈先加速增长、后缓慢发展、最后急速降低的特征;钢拱架轴力快速降低时,弯矩剧烈变化,部分监测点出现由负变正的特征;侧压力系数越大,钢拱架屈服时间越早,轴力降低和弯矩剧烈变化时间也越早.
4) 侧压力系数对衬砌受压损伤、受拉损伤区的空间分布影响不明显,受压损伤破坏区主要分布在隧道边墙和墙脚位置,受拉损伤首先出现在边墙位置的二次衬砌表面;拱部也会形成受压损伤区,且侧压力系数越大损伤越严重;侧压力系数越大,衬砌中形成贯通的受压破坏区、受拉损伤达到最大值的时间越早.
致谢:西华大学人才引进项目(Z231016)资助.
-
表 1 锚杆物理力学参数
Table 1. Physical and mechanical parameters of anchor bolt
弹性模量E/GPa 抗拉强度Ft/MPa 延伸率/% 直径d/mm 200 335 16 22 表 2 钢拱架物理力学参数
Table 2. Physical and parameters of steel arch frame
钢拱架型号 弹性模量E/GPa 泊松比v 极限轴力Nu/kN 极限弯矩Mu/kN•m I20a 200 0.3 849 71 表 3 混凝土力学参数
Table 3. Mechanical parameters of concrete
强度
等级弹性模量E /GPa 泊松比ν 内聚力c/MPa 内摩擦角φ/(°) 抗压强度σc /MP 抗拉强度σt /MPa C25 28.0 0.2 3.0 50 16.7 1.78 C35 31.5 0.2 3.7 55 23.4 2.20 表 4 混凝土弹塑性损伤参数
Table 4. Elasto-plastic damage parameters of concrete
C25混凝土 C35混凝土 κs/
( × 10−3)Ds κt/
( × 10−3)Dt κs/
( × 10−3)Ds κt/
( × 10−3)Dt 0 0 0 0 0 0 0 0 0.06 0.01 0.42 0.55 0.07 0.01 0.45 0.55 0.20 0.02 1.02 0.75 0.23 0.02 1.13 0.75 0.45 0.03 3.60 0.99 0.56 0.04 3.92 0.99 1.20 0.11 / / 1.40 0.12 / / 3.00 0.40 / / 2.90 0.39 / / 4.50 0.70 / / 4.10 0.67 / / 7.00 0.99 / / 6.00 0.99 / / 表 5 围岩蠕变力学参数
Table 5. Creep mechanical parameters of surrounding rock
K/
GPaEM/
GPaEK/
GPaηM/
GPa.aηK/
(GPa.a)c/
MPaφ/(°) σt/
MPa2 2.5 4 3 1.2 7 26 1 -
[1] 关宝树. 漫谈矿山法隧道技术第四讲:钢架[J]. 隧道建设,2016,36(2): 123-130. doi: 10.3973/j.issn.1672-741X.2016.02.001GUAN Baoshu. Tunneling by mining method: lecture Ⅳ: steel arch[J]. Tunnel Construction, 2016, 36(2): 123-130. doi: 10.3973/j.issn.1672-741X.2016.02.001 [2] 李天胜,何川,方砚兵,等. 基于围岩变形失效的隧道结构可靠度设计方法[J]. 西南交通大学学报,2023,58(3): 613-621. doi: 10.3969/j.issn.0258-2724.20220137LI Tiansheng, HE Chuan, FANG Yanbing, et al. Reliability-based design method of tunnel structures based on deformation failure of surrounding rock[J]. Journal of Southwest Jiaotong University, 2023, 58(3): 613-621. doi: 10.3969/j.issn.0258-2724.20220137 [3] 潘钦锋,张丙强,黄志斌. 隧道下穿诱发既有管道-土体非协调变形解析研究[J]. 西南交通大学学报,2024,59(3): 637-645. doi: 10.3969/j.issn.0258-2724.20230334PAN Qinfeng, ZHANG Bingqiang, HUANG Zhibin. Analytical study for uncoordinated deformation of existing pipeline and soil induced by tunnel undercrossing[J]. Journal of Southwest Jiaotong University, 2024, 59(3): 637-645. doi: 10.3969/j.issn.0258-2724.20230334 [4] 中华人民共和国交通运输部. 公路隧道设计细则:JTG/T D70—2010[S]. 北京:人民交通出版社,2010. [5] 国家铁路局. 铁路隧道设计规范:TB 10003—2016[S]. 北京:中国铁道出版社,2017. [6] 周建,杨新安,蔡键,等. 深埋复合式衬砌隧道二衬分担比研究[J]. 公路交通科技,2020,37(5): 92-99. [7] 王耀东,朱力强,余祖俊,等. 基于样本自动标注的隧道裂缝病害智能识别[J]. 西南交通大学学报,2023,58(5): 1001-1008,1036. doi: 10.3969/j.issn.0258-2724.20210092WANG Yandong, ZHU Liqiang, YU Zujun, et al. Intelligent tunnel crack recognition based on automatic sample labeling[J]. Journal of Southwest Jiaotong University, 2023, 58(5): 1001-1008,1036. doi: 10.3969/j.issn.0258-2724.20210092 [8] 中华人民共和国交通运输部. 公路工程技术标准:JTG B01—2014[S]. 北京:人民交通出版社,2015. [9] 陈志豪. 锈蚀加筋喷砼拱肋支护承载性能劣化与损伤评价[D]. 长沙:长沙理工大学,2021. [10] 丁万涛,李术才. 锈蚀对海底隧道锚固支护结构加固性能的影响[J]. 华南理工大学学报(自然科学版),2013,41(6): 100-107,139. doi: 10.3969/j.issn.1000-565X.2013.06.016DING Wantao, LI Shucai. Effect of corrosion on reinforcement performance of anchorage support structure of subsea tunnel[J]. Journal of South China University of Technology (Natural Science Edition), 2013, 41(6): 100-107,139. doi: 10.3969/j.issn.1000-565X.2013.06.016 [11] 杜建明,房倩,程荔琼,等. 围岩-支护结构时效性共同作用下运营隧道服役特性[J]. 中南大学学报(自然科学版),2022,53(10): 4001-4011. doi: 10.11817/j.issn.1672-7207.2022.10.020DU Jianming, FANG Qian, CHENG Liqiong, et al. Service characteristics of operational tunnel under combined action of time effects of surrounding rock and supporting structure[J]. Journal of Central South University (Science and Technology), 2022, 53(10): 4001-4011. doi: 10.11817/j.issn.1672-7207.2022.10.020 [12] 席红兵,李柏生. 硫酸盐-冻融共同作用下隧道衬砌支护喷射混凝土劣化性能研究[J]. 隧道建设(中英文),2022,42(7): 1219-1226.XI Hongbing, LI Bo Sheng. Deterioration performance of tunnel lining support shotcrete under sulfate-freeze-thaw interaction[J]. Tunnel Construction, 2022, 42(7): 1219-1226. [13] 杜学才,周辉,郭鹏云,等. 香炉山引水隧洞围岩的蠕变特性及模型研究[J]. 长江科学院院报,2022,39(12): 75-81. doi: 10.11988/ckyyb.20221011DU Xuecai, ZHOU Hui, GUO Pengyun, et al. Creep characteristics and modelling of surrounding rock of Xianglushan water diversion tunnel[J]. Journal of Yangtze River Scientific Research Institute, 2022, 39(12): 75-81. doi: 10.11988/ckyyb.20221011 [14] 封坤,王胤丞,马文帅,等. 围岩蠕变对盾构隧道受荷特征影响研究[J]. 铁道标准设计,2022,66(8): 117-124,131.FENG Kun, WANG Yincheng, MA Wenshuai, et al. Influence of surrounding rock creep on loading characteristics of shield tunnel[J]. Railway Standard Design, 2022, 66(8): 117-124,131. [15] 经纬,陈洪恩,杨仁树,等. 基于岩石蠕变及D-P准则的深部巷道围岩弹塑性分析[J]. 力学学报,2022,54(7): 1982-1993. doi: 10.6052/0459-1879-21-619JING Wei, CHEN Hongen, YANG Renshu, et al. Elastoplastic analysis of surrounding rock of deep roadway based on rock creep and d-p criterion[J]. Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(7): 1982-1993. doi: 10.6052/0459-1879-21-619 [16] 路军富,肖铮,喻渝,等. 铁路隧道底鼓段围岩蠕变参数反演方法研究[J]. 铁道工程学报,2021,38(1): 66-71. doi: 10.3969/j.issn.1006-2106.2021.01.012LU Junfu, XIAO Zheng, YU Yu, et al. Research on the back analysis method of creep parameters of surrounding rock at floor heave section of railway tunnel[J]. Journal of Railway Engineering Society, 2021, 38(1): 66-71. doi: 10.3969/j.issn.1006-2106.2021.01.012 [17] 张海洋,徐文杰,王永刚. 各向异性软岩隧道围岩蠕变特征及支护方案对比分析[J]. 铁道建筑,2018,58(8): 71-74. doi: 10.3969/j.issn.1003-1995.2018.08.18ZHANG Haiyang, XU Wenjie, WANG Yonggang. Creep characteristics of surrounding rock in anisotropic soft rock tunnel and comparative analysis of supporting schemes[J]. Railway Engineering, 2018, 58(8): 71-74. doi: 10.3969/j.issn.1003-1995.2018.08.18 [18] YANG K, YAN Q X, SHI Z D, et al. Numerical study on the mechanical behavior of shotcrete lining with yielding support in large deformation tunnel[J]. Rock Mechanics and Rock Engineering, 2023, 56(2): 1563-1584. doi: 10.1007/s00603-022-03126-w [19] XU G W, GUTIERREZ M. Study on the damage evolution in secondary tunnel lining under the combined actions of corrosion degradation of preliminary support and creep deformation of surrounding rock[J]. Transportation Geotechnics, 2021, 27: 100501. doi: 10.1016/j.trgeo.2020.100501 [20] 戴永浩,陈卫忠,于洪丹,等. 大坂膨胀性泥岩引水隧洞长期稳定性分析[J]. 岩石力学与工程学报,2010,29(增1): 3227-3234.DAI Yonghao, CHEN Weizhong, YU Hongdan, et al. Long-term stability analysis of daban diversion tunnel buried in swelling mudstone[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(S1): 3227-3234. [21] 李建军,张志强. 岩石蠕变对隧道二次衬砌结构影响的研究[J]. 现代隧道技术,2011,48(6): 58-64,81. doi: 10.3969/j.issn.1009-6582.2011.06.012LI Jianjun, ZHANG Zhiqiang. Study of the influence of rock creep on the structure of a tunnel’s secondary lining[J]. Modern Tunnelling Technology, 2011, 48(6): 58-64,81. doi: 10.3969/j.issn.1009-6582.2011.06.012 [22] 徐国文,何川,汪耀,等. 流变荷载作用下隧道裂损二次衬砌结构安全性能研究[J]. 土木工程学报,2016,49(12): 114-123.XU Guowen, HE Chuan, WANG Yao, et al. Study on the safety performance of cracked secondary lining under action of rheological load[J]. China Civil Engineering Journal, 2016, 49(12): 114-123. [23] 王迎超,尚岳全,孙红月,等. 复合式衬砌在围岩蠕变过程中的受力规律研究[J]. 水文地质工程地质,2010,37(2): 49-54. doi: 10.3969/j.issn.1000-3665.2010.02.011WANG Yingchao, SHANG Yuequan, SUN Hongyue, et al. Study on mechanical rules of double-lining in creeping surrounding rock[J]. Hydrogeology & Engineering Geology, 2010, 37(2): 49-54. doi: 10.3969/j.issn.1000-3665.2010.02.011 [24] 师亚龙,郑波,刘志强,等. 流变效应作用下高地应力铁路软岩隧道衬砌结构长期稳定性研究[J]. 铁道标准设计,2017,61(10): 95-99.SHI Yalong, ZHENG Bo, LIU Zhiqiang, et al. Study on long term stability of tunnel lining structure in high stress soft rock tunnel under the effect of rheology[J]. Railway Standard Design, 2017, 61(10): 95-99. [25] 钱文喜,耿大新,梁国卿. 围岩蠕变对运营隧道衬砌安全性的影响[J]. 公路交通科技,2020,37(9): 90-96. doi: 10.3969/j.issn.1002-0268.2020.09.012 [26] 汪耀. 高地应力下围岩蠕变对隧道结构长期安全性的影响分析[D]. 成都:西南交通大学,2017. [27] Itasca Consulting Group Inc. Fast lagrangian analyses of continua in three-dimensions[CP/OL]. (2012-06-23)[2023-08-12]. https://www.itascacg.com/software/FLAC3D. [28] Hibbitt, Karlsson and Sorensen, Inc. Abaqus analysis user’s guide[CP/OL]. (2014-05-06)[2023-08-12]. https://www.3ds.com/products/simulia/abaqus. [29] Ansys. Mechanical 12.0 user’s guide & theory guide [CP/OL]. (2009-10-14)[2023-08-12]. https://www.ansys.com/. [30] 徐国文. 层状千枚岩地层隧道稳定性分析[D]. 成都:西南交通大学,2017. [31] 过镇海,时旭东. 钢筋混凝土原理和分析[M]. 北京:清华大学出版社,2003. [32] 中华人民共和国建设部. 混凝土结构设计规范:GB 50010—2002[S]. 北京:中国建筑工业出版社,2004. [33] YANG K, YAN Q X, ZHANG C, et al. Study on mechanical properties and damage evolution of carbonaceous shale under triaxial compression with acoustic emission[J]. International Journal of Damage Mechanics, 2021, 30(6): 899-922. doi: 10.1177/1056789521991193 [34] 刘新喜,李盛南,徐泽沛,等. 冻融循环作用下炭质页岩蠕变模型研究[J]. 中国公路学报,2019,32(11): 137-145.LIU Xinxi, LI Shengnan, XU Zepei, et al. Research on creep model of carbonaceous shale under freeze-thaw cycle[J]. China Journal of Highway and Transport, 2019, 32(11): 137-145. -