摘要
地表超载是威胁盾构隧道结构安全运营的重要因素,且随着隧道服役时间的增长,钢筋锈蚀导致的盾构隧道性能劣化可能进一步加剧地表超载作用下结构损伤破坏的风险.鉴于此,提出了地表超载作用下锈蚀盾构隧道结构易损性分析方法,以软土地区某典型浅埋盾构隧道为例,建立了地表超载作用下锈蚀盾构隧道管片三维精细化有限元分析模型,考虑地层抗力系数与土体重度的不确定性开展了数值分析;在此基础上,以接头张开量为损伤指标,采用双参数对数正态分布模型建立了地表超载作用下锈蚀盾构隧道结构易损性曲线,揭示了超载和钢筋锈蚀率等因素对隧道结构易损性的影响规律.研究结果表明,隧道易损性概率随地表超载水平与钢筋锈蚀率的增大而增大;当地表超载较小时,钢筋锈蚀率对隧道易损性的影响更为突出;同一工况下,隧道腰部接头相较于拱顶或拱底接头具有更高的损伤概率.
近年来,随着我国城市建设的蓬勃发展,各大城市纷纷加快地铁建设以应对日益严重的交通拥堵问题.截至2022年底,我国已有53个城市开通城市轨道交通,总运营里程达到9 584 km,地铁已成为沟通城市各大功能分区的重要交通方式.然而,随着城市地下空间开发的不断深入,不可避免地需要在既有地铁隧道附近开展各项工程建设活动.其中,地表超载作为一种典型的外部扰动对既有盾构隧道的运营安全带来了严重威胁.如2008年上海某地铁区间盾构隧道上方突发大量堆土,最大堆土高度约7 m,事故进一步引发了隧道结构横向收敛变形、接缝渗漏水、衬砌管片开裂等病害,严重威胁隧道结构与运营安
目前,国内外众多学者已通过解析分析、数值模拟和试验研究等方法对地表超载下隧道受力与变形特性展开了大量研究.在解析方面,学者们基于经典理论模型对地表超载下隧道横向与纵向受力与变形展开分析.如房营光
上述工作丰富了对超载作用下隧道结构受力与变形规律的认识.然而,随着盾构隧道服役时间的增加,钢筋锈蚀等因素将导致隧道力学性能不可避免地产生不同程度的劣
另外,易损性分析是合理评价隧道结构安全状态的重要手段之一,其中,易损性的概念源于军事领域,意指系统对外部灾害作用的承受能力.目前,易损性分析方法广泛应用于地震、滑坡、恐怖袭击和泥石流等灾害作用下建筑、桥梁、隧道等结构的安全评
因此,本文首先提出了地表超载作用下锈蚀盾构隧道结构易损性分析流程;以软土地区典型盾构隧道为例,建立了超载作用下锈蚀盾构隧道三维有限元分析模型;通过考虑地层抗力系数与土体重度的不确定性,采用蒙特卡洛方法开展大量数值计算;基于数值计算结果,进一步以隧道接头张开量为损伤指标,建立了地表超载作用下锈蚀盾构隧道结构易损性曲线,揭示了超载作用和钢筋锈蚀等因素对盾构隧道结构易损性的影响规律.
1 盾构隧道结构易损性分析方法
1.1 盾构隧道结构损伤指标及其等级
开展地表超载作用下盾构隧道结构易损性分析的重要工作之一是确定评价结构损伤的指标及其分级标准.根据柳献、吴庆、刘四进等学
在盾构隧道结构易损性研究中通常将损伤指标(DI)分为5个等级,分别对应无破坏D0,轻微破坏D1,中等破坏D2,严重破坏D3和完全破坏D4.本文结合沈贤
损伤等级 | D0 | D1 | D2 | D3 | D4 |
---|---|---|---|---|---|
接头张开量/mm | 0~2 | >2~4 | >4~6 | >6~8 | >8 |
1.2 盾构隧道结构易损性曲线的定义
本文采用双参数对数正态分布模
(1) |
式中:P(Damage≥Di|P0)表示地表超载大小为P0时隧道结构损伤等级达到Di的概率;结构损伤等级通过结构损伤指标来表征,在本文的研究中,选取隧道接头张开量作为结构损伤指标,详见1.1节;Φ为标准正态分布的累计分布函数;x为意外荷载的强度水平,即易损性曲线的横坐标,在本文的研究中对应地表超载P0;mxi为当结构损伤等级为i时,对应的地表超载水平P0的对数均值;βxi为第i级损伤状态对应的总对数标准差,描述易损性曲线的总体变异性.本研究中考虑易损性曲线的总体变异性主要来源于三个方面,即隧道损伤状态确定及分级范围的不确定性βDi、隧道结构响应计算结果不确定性βC以及地表超载的不确定性βl,假设三者相互独立,则βxi可表示为
(2) |
根据美国联邦紧急事务管理局与国家建筑科学研究所的相关研
(3) |
式中:k为回归曲线的斜率,b为回归曲线的截距.基于

图1 mxi与βl参数确定方法示意图
Fig.1 Method of determining the parameter mxi and βl
1.3 地表超载作用下锈蚀盾构隧道易损性分析流程
本文提出的地表超载作用下锈蚀盾构隧道的易损性曲线的分析流程如

图2 锈蚀盾构隧道易损性分析流程
Fig.2 Fragility analysis of corroded shield tunnel
1) 选取典型的盾构隧道结构形式与材料参数,确定钢筋锈蚀率与地表超载水平的取值;
2) 确定不确定性参数及其概率密度分布形式,如本文选取了地层抗力系数与土体重度作为不确定参数,并对参数进行抽样;
3) 根据步骤1)与步骤2)中的相关参数建立地表超载作用下锈蚀盾构隧道的有限元计算模型,确定典型工况开展大量数值模拟计算;
4) 确定地表超载作用下盾构隧道结构损伤指标及结构损伤等级的阈值范围,本文选取盾构隧道接头张开量作为结构损伤指标,结构损伤等级的划分方法见
5) 基于数值计算结果与步骤1)中的结构损伤指标及其阈值,采用
6) 根据步骤5)中对两个关键参数的计算结果,结合
2 锈蚀盾构隧道数值分析模型
2.1 锈蚀盾构隧道有限元分析模型
本文采用有限元软件ABAQUS建立盾构隧道计算模型来模拟地表超载作用下管片的力学性能.由于后续需进行大量参数抽样与数值计算,因此采用荷载结构法模拟管片受到的外力作用.荷载结构模型如

图3 荷载结构法计算模型
Fig.3 Calculation model of load structure method
该模型是国际隧道协会推荐的隧道设计方法,在隧道设计与建造中被广泛使用.模型考虑隧道的自重应力、土体应力与水压力,并采用水土分算的形式.其荷载可以用p1~p5表示,其计算方法见式(4)~
(4) |
(5) |
(6) |
(7) |
(8) |
式中:P0为地表超载水平(kPa);γ为土体饱和重度(kN/
本文基于国内较为常用的地铁盾构隧道结构建立了一环衬砌精细化数值模型来开展后续分析.隧道外径为6.2 m,内径为5.5 m.衬砌环宽1 200 mm,由一块封顶块F(16°)、两块标准块B(65°)、两块邻接块L(65°)和一块封底块D(84°)共6片管片组成.管片间建立螺栓孔与手孔模型,并通过2根6.8级的M30直螺栓连接.计算模型如

图4 盾构隧道有限元模型
Fig.4 Finite element model of shield tunnel
所有螺栓长度均为400 mm,螺栓孔直径为 42 mm.模型管片间采用直接接触,设置接触面法向为硬接触关系,切向为摩擦接触关系,摩擦系数为0.4.螺母内侧与手孔断面设置Tie接触以保证螺栓与管片保持共同变形.
基于学者们对软土地区广泛分布的典型土层的研究,模型中土体的侧压力系数λ设为0.6.同时,考虑浅埋隧道较深埋隧道更易受到地表超载作用的影响,本文的研究只考虑隧道的浅埋工况,选取典型的浅埋隧道埋深(顶板覆土厚度) H=10 m.
为确保模型计算精度与计算效率,对模型不同的几何部分分别进行网格划分.管片总体设置为六面体结构化网格单元,局部手孔与螺栓部位划分为扫掠网格.为提高模型的计算效率,模型中的材料均采用线弹性本构,将钢筋混凝土管片简化为均质的混凝土材料,混凝土强度等级为C55,其弹性模量为Ec=3.55×1
既有研究表明,锈蚀造成的管片内钢筋有效截面面积减小或锈胀引起的混凝土与钢筋的黏结锚固变化是钢筋锈蚀诱发盾构隧道结构性能劣化的两个主要因素.上述两个因素都会直接导致管片整体抗弯刚度低于设计值.本文假设钢筋发生均匀锈蚀,通过降低管片的等效抗弯刚度反映钢筋锈蚀对管片性能的影响.根据Cabrer
(9) |
(10) |
式中:Ds为钢筋的设计直径;Dc为锈蚀后钢筋剩余的有效直径;η为钢筋的面积锈蚀率;EI为管片的设计等效刚度;EIcorr为锈蚀后管片的剩余等效刚度.基于上述理论,本文通过改变有限元模型中管片的弹性模量来表达不同钢筋锈蚀率下管片的劣化性能.
2.2 模型验证
由于开展整环锈蚀隧道模型试验需要耗费大量的经济与时间成本,因此在相关领域鲜有可以参考的研究成果,这也为模型的直接验证带来了困难.研究表明,管片的接头性能是决定隧道结构力学特性的关键因素.因此,本文通过开展管片的接头试验,并以接头张开量为指标对比了试验数据与数值计算结果来验证上述有限元模型的合理性,试验模型如

(a) 管片接头力学性能试验

(b) 接头张开量有限元与试验结果对比
图5 接头力学性能的试验验证
Fig.5 Experimental verification on mechanical performance of joint
2.3 不确定性参数与地表超载水平
本研究考虑地层抗力系数与土体重度的变异性.假定地层抗力系数k服从对数正态分布,土体重度γ服从正态分布,其特征参数如
参数 | 分布形式 | 均值 | 标准差 |
---|---|---|---|
地层抗力系数k/(kN· | 对数正态分布 | 0.6 | 0.1 |
土体重度γ /(kN· | 正态分布 | 18 | 1.8 |
盾构隧道的钢筋锈蚀程度往往随着隧道运营时间的增加逐渐加深.现有数据表明,上海地区盾构隧道钢筋的锈蚀速率可按0.02 mm/年计算.以锈蚀诱导期10年计算,我国大多数盾构隧道在达到30年的运营使用期时对应的钢筋锈蚀率为20%,作为轻度锈蚀状态;100年的设计使用寿命时对应的钢筋锈蚀率为70%,作为重度锈蚀状态;同时取运营时间为70年对应的钢筋锈蚀率50%表征中等锈蚀状态.此外,考虑地表超载最大为200 kPa,模型的计算工况如
工况 | 钢筋锈蚀率/% | 地表超载水平/kPa |
---|---|---|
1 | 20 | 0~200 |
2 | 50 | 0~200 |
3 | 70 | 0~200 |
3 超载作用下锈蚀隧道易损性评价
基于第2节提出的地表超载作用下锈蚀盾构隧道结构易损性曲线分析方法,对

(a) 接头1张开量
(b) 接头2张开量

(c) 接头3张开量
图6 钢筋锈蚀率为20%时盾构隧道结构易损性曲线
Fig.6 Fragility curve of shield tunnel when the corrosion rate of steel bar is 20%

(a) 接头1张开量
(b) 接头2张开量

(c) 接头3张开量
图7 钢筋锈蚀率为50%时盾构隧道结构易损性曲线
Fig.7 Fragility curve of shield tunnel when the corrosion rate of steel bar is 50%

(a) 接头1张开量
(b) 接头2张开量

(c) 接头3张开量
图8 钢筋锈蚀率为70%时盾构隧道结构易损性曲线
Fig.8 Fragility curve of shield tunnel when the corrosion rate of steel bar is 70%

(a) 接头1张开量

(b) 接头2张开量

(c) 接头3张开量
图9 不同钢筋锈蚀率隧道结构损伤指标超越D0的概率
Fig.9 Probability of structure damage index reaching D0 level under different reinforcement corrosion rate

(a) 接头1张开量

(b) 接头2张开量

(c) 接头3张开量
图10 不同钢筋锈蚀率隧道结构损伤指标超越D1的概率
Fig.10 Probability of structure damage index reaching D1 level under different reinforcement corrosion rate

(a) 接头1张开量

(b) 接头2张开量

(c) 接头3张开量
图11 不同钢筋锈蚀率隧道结构损伤指标超越D2的概率
Fig.11 Probability of structure damage index reaching D2 level under different reinforcement corrosion rate
3.1 钢筋锈蚀率为20%时隧道结构易损性曲线
3.2 钢筋锈蚀率为50%时隧道结构易损性曲线
3.3 钢筋锈蚀率为70%时隧道结构易损性曲线
上述研究表明,在同一钢筋锈蚀率下,隧道内拱顶接头、腰部接头与拱底接头超越某一特定损伤等级的概率均随着地表超载水平的增加而增加.同时,对比同一工况下不同位置接头的损伤概率,腰部接头比拱顶或拱底接头具有更高的损伤概率.
3.4 不同钢筋锈蚀率下隧道结构易损性对比
本节分别选取了50 kPa、100 kPa、150 kPa与200 kPa的地表超载水平,对不同钢筋锈蚀率下结构损伤状态超越D0、D1、D2损伤等级上限的概率进行了对比分析,揭示了不同钢筋锈蚀率对盾构隧道易损性的影响规律,对比结果如
当以接头1张开量为损伤指标,地表超载水平达到50 kPa时,钢筋锈蚀率为20%时结构损伤概率为37%,较钢筋锈蚀率为50%与70%时低18%,地表超载水平达到200 kPa,钢筋锈蚀率为20%时结构损伤概率为85%,较钢筋锈蚀率为50%与70%时低10%;当以接头2张开量为损伤指标时,在所有情况下,结构损伤指标超越D0损伤等级的概率均为100%;当以接头3张开量为损伤指标时,钢筋锈蚀率为20%时,结构损伤概率在地表超载水平为50 kPa时达到77%,地表超载水平为200 kPa时达到99%,钢筋锈蚀率超过50%时,隧道结构损伤状态超越D0损伤等级的概率始终为100%.
4 结 论
本文构建了超载条件下锈蚀盾构隧道易损性分析方法,并基于大量数值计算研究了超载条件下不同锈蚀盾构隧道结构易损性发展规律.得出的主要结论如下:
1)锈蚀盾构隧道的易损性随着地表超载水平的增大而增大,在同一地表超载水平下,随着钢筋锈蚀率的增加,隧道结构易损性逐渐增加.
2)当地表超载较小时,钢筋锈蚀对隧道易损性的影响更为突出.例如,当地表超载水平为200 kPa时,钢筋锈蚀率为70%接头1发生超越D0损伤等级的概率较钢筋锈蚀率为20%时高10%,而当地表超载水平为50 kPa时,对应的差值增加为18%.
3)在相同工况下,隧道腰部接头相对于拱顶与拱底上限接头具有更高的损伤概率.例如,当地表超载水平为200 kPa时,20%钢筋锈蚀率下腰部接头张开量超越D1损伤等级上限的概率为94%,而位于拱顶与拱底接头超越概率分别为63%与78%.因此,在盾构隧道日常运营维护中应加强对腰部接头的监测与保护.
参考文献
HUANG H W, ZHANG D M. Resilience analysis of shield tunnel lining under extreme surcharge: Characterization and field application[J]. Tunnelling and Underground Space Technology, 2016, 51(1): 301-312. [百度学术]
邵华,黄宏伟,张东明,等.突发堆载引起软土地铁盾构隧道大变形整治研究[J].岩土工程学报,2016,38(6):1036-1043. [百度学术]
SHAO H,HUANG H W, ZHANG D M, et al. Case study on repair work for excessively deformed shield tunnel under accidental surface surcharge in soft clay[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(6): 1036-1043.(in Chinese) [百度学术]
房营光, 孙钧.地面荷载下浅埋隧道围岩的粘弹性应力和变形分析[J].岩石力学与工程学报,1998,17(3):239-247. [百度学术]
FANG Y G,SUN J.Analysis of visco-elastic stress and deformation of surrounding rocks for shallow tunnel under ground load[J].Chinese Journal of Rock Mechanics and Engineering,1998, 17(3): 239-247.(in Chinese) [百度学术]
姜启元, 管攀峰, 叶蓉. 软土盾构隧道的纵向变形分析[J]. 地下工程与隧道, 1999, 2(4): 21-46. [百度学术]
JIANG Q Y, GUAN P F, YE R. The analysis of longitudinal deformation of shield tunnel in soft soil[J]. Underground Engineering and Tunnels, 1999, 2(4): 21-46. (in Chinese) [百度学术]
康成, 梅国雄, 梁荣柱, 等.地表临时堆载诱发下既有盾构隧道纵向变形分析[J].岩土力学, 2018, 39(12): 4605-4616. [百度学术]
KANG C,MEI G X,LIANG R Z,et al.Analysis of the longitudinal deformation of existing shield tunnel induced by temporary surface surcharge[J].Rock and Soil Mechanics, 2018, 39(12): 4605-4616.(in Chinese) [百度学术]
魏纲, 张书鸣, 余剑英, 等.地面堆载对盾构隧道围压影响的模型试验与理论分析[J].岩土工程学报, 2022, 44(10):1789-1798. [百度学术]
WEI G,ZHANG S M,YU J Y,et al.Model tests and theoretical analyses of influences of surface surcharge on confining pressure of shield tunnels[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(10): 1789-1798.(in Chinese) [百度学术]
YAMAMOTO K, LYAMIN A V, WILSON D W, et al. Stability of a circular tunnel in cohesive-frictional soil subjected to surcharge loading[J]. Computers and Geotechnics,2011, 38(4): 504-514. [百度学术]
叶耀东. 软土地区运营地铁盾构隧道结构变形及健康诊断方法研究[D]. 上海: 同济大学, 2007. [百度学术]
YE Y D. Study on structural deformation and health diagnosis of shield tunnel in soft soil area[D]. Shanghai: Tongji University, 2007. (in Chinese) [百度学术]
张东明, 代鉷锋, 王慧, 等.考虑地层变异的浅基础承载力分析[J].地下空间与工程学报,2020,16(5):1412-1419. [百度学术]
ZHANG D M,DAI H F,WANG H,et al.Analysis on shallow foundation capacity considering geological heterogeneity[J].Chinese Journal of Underground Space and Engineering, 2020, 16(5):1412-1419.(in Chinese) [百度学术]
王如路, 张冬梅.超载作用下软土盾构隧道横向变形机理及控制指标研究[J].岩土工程学报,2013, 35(6):1092-1101. [百度学术]
WANG R L,ZHANG D M.Mechanism of transverse deformation and assessment index for shield tunnels in soft clay under surface surcharge[J]. Chinese Journal of Geotechnical Engineering,2013,35(6): 1092-1101.(in Chinese) [百度学术]
柳献, 张浩立, 鲁亮, 等. 超载工况下盾构隧道结构承载能力的试验研究[J]. 地下工程与隧道, 2013, 4(1): 10-15. [百度学术]
LIU X, ZHANG H L, LU L, et al. Experimental study on bearing capacity of shield tunnel structure under overload condition[J]. Underground Engineering and Tunnels, 2013, 4(1): 10-15. (in Chinese) [百度学术]
吴庆,杜守继.地面堆载对既有盾构隧道结构影响的试验研究[J].地下空间与工程学报,2014,10(1): 57-66. [百度学术]
WU Q, DU S J.Model test on influence of ground heaped load on existing shield tunnel structure[J]. Chinese Journal of Underground Space and Engineering, 2014, 10(1): 57-66.(in Chinese) [百度学术]
黄大维,周顺华,赖国泉,等.地表超载作用下盾构隧道劣化机理与特性[J].岩土工程学报,2017,39(7):1173-1181. [百度学术]
HUANG D W,ZHOU S H,LAI G Q,et al.Mechanisms and characteristics for deterioration of shield tunnels under surface surcharge[J].Chinese Journal of Geotechnical Engineering,2017,39(7): 1173-1181.(in Chinese) [百度学术]
刘四进,何川,封坤,等.受荷状态下盾构隧道管片锈蚀劣化破坏过程研究[J].土木工程学报, 2018,5 1(6): 120-128. [百度学术]
LIU S J,HE C,FENG K,et al.Research on corrosion deterioration and failure process of shield tunnel segments under loads[J].China Civil Engineering Journal,2018,51(6):120-128.(in Chinese) [百度学术]
HE Z S,AKIYAMA M,HE C,et al.Life-cycle reliability analysis of shield tunnels in coastal regions:emphasis on flexural performance of deteriorating segmental linings[J].Structure and Infrastructure Engineering, 2019, 15(7): 851-871. [百度学术]
NIELSON B G,DESROCHES R.Analytical seismic fragility curves for typical bridges in the central and southeastern United States[J].Earthquake Spectra, 2007, 23(3): 615-633. [百度学术]
吴俊陶. 基于易损性的某连体复杂超高层结构强震损伤控 制[D]. 广州: 广州大学, 2021. [百度学术]
WU J T. Study on seismic damage control of connected super high-rise structure base on seismic fragility analysis[D]. Guangzhou: Guangzhou University, 2021. (in Chinese) [百度学术]
沈贤达. 软土盾构隧道结构易损性分析[D]. 上海: 同济大学, 2015. [百度学术]
SHEN X D. Vulnerability analysis of shield tunnel in soft soil[D]. Shanghai: Tongji University, 2015. (in Chinese) [百度学术]
禹海涛,李心熙,袁勇,等.沉管隧道纵向地震易损性分析方法[J].中国公路学报,2022,35(10):13-22. [百度学术]
YU H T, LI X X, YUAN Y, et al.Seismic vulnerability analysis method for longitudinal response of immersed tunnels[J]. China Journal of Highway and Transport, 2022, 35(10): 13-22.(in Chinese) [百度学术]
NIBS. HAZUS-MH: Technical manuals [M]. Washington D C:Federal Emergency Management Agency and National Institute of Building Science, 2004. [百度学术]
CABRERA J G. Deterioration of concrete due to reinforcement steel corrosion[J]. Cement and Concrete Composites, 1996, 18(1): 47-59. [百度学术]
BALLIM Y,REID J C.Reinforcement corrosion and the deflection of RC beams—an experimental critique of current test methods[J]. Cement and Concrete Composites, 2003, 25(6): 625-632. [百度学术]