刘学增, 赖浩然, 桑运龙,3, 张 强, 张竹清, 刘金栋
(1. 同济大学 土木信息技术教育部工程研究中心, 上海 200092; 2. 同济大学 土木工程学院, 上海 200092;3. 上海同岩土木工程科技股份有限公司, 上海 200092; 4. 上海地下基础设施安全检测与养护装备工程技术研究中心,上海 200092; 5. 厦门轨道交通集团有限公司, 福建 厦门 361008;6. 青岛银盛泰集团有限公司, 山东 济南 266000)
随着城市地铁路网的不断发展,地铁隧道保护区基坑施工也逐渐增多,基坑施工对周边地层的应力释放作用,不可避免地会对地铁隧道造成一定的影响,严重时造成地铁隧道结构损伤[1-3].
目前,针对盾构隧道结构力学性能已开展了较多研究[4-8],但针对卸载工况下盾构隧道结构受力变形特征的研究尚少.柳献等[9]针对卸载工况下盾构隧道结构承载能力进行了足尺试验,记录了试验现象及相关破坏形式,并得到了盾构隧道在卸载工况下的受力变形机理及关键性能点,但由于足尺试验规模大、周期长,地层难以模拟,因此没有考虑地层与结构之间的相互作用.徐秀峰等[10]针对上海地铁16号线外径为11.6 m的大直径盾构隧道,建立包含内部结构的错缝拼装隧道三维模型,对隧道两侧土体对称开挖卸载工况下管片环的直径变形量、接缝张开量、中隔墙应力及石岩棉压缩量与土体开挖深度的关系进行了分析,主要以变形预测为主,缺乏对管片结构承载性能演化规律的研究.
基于几何相似比1∶5的室内模型试验以及三维有限元数值建模分析,研究淤泥地层双侧卸载工况下盾构隧道错缝拼装结构受力变形特征、损伤过程以及破坏模式,寻求隧道管片应力与变形、构件承载力与变形之间的对应规律,并提出双侧卸载工况下的隧道变形控制指标,为今后盾构隧道保护区内基坑施工变形控制提供依据.
模型试验以厦门地铁2号线厦门岛至海沧段为背景,选取地铁隧道所穿越的陆域段淤泥质地层作为模型试验条件,地质横断面如图1所示.每环管片由1个封顶块(K)、2个邻接块(B1、B2)和3个标准块(A1、A2、A3)组成,按45°错缝拼装.混凝土强度等级为C55,管片环间由16个环缝连接螺栓(M30)连接,管片间由12个纵缝连接螺栓(M30)连接.隧道外径为6 700 mm,内径为6 000 mm,环宽为1 500 mm,管片厚度为350 mm,分块如图2所示.
图1 厦门地铁2号线陆域段工程横断面
取几何相似比CL=5、容重相似比Cγ=1作为基础相似比,根据相似理论,其他物理力学参数相似比如表1所示.
管片由水、中砂和水泥按1.0∶6.0∶1.5的质量比配置而成,并通过模具浇筑成型,如图3所示.管片钢筋采用直径为2 mm、网格大小为40 mm×15 mm的镀锌铁丝网模拟.试验中不考虑管片手孔构造,螺栓采用M6 AZ91d镁合金模拟.每环管片外侧等距布设36个径向弹簧模拟地层.管片采用ABA型错缝拼装模式,相邻环错动角度为45°.
图2 厦门地铁2号线管片衬砌分块图
物理量相似常数几何尺寸/m5位移/m5弹性模量/(N·m-2)5应力/(N·m-2)5弹性抗力系数/(N·m-3)1应变1泊松比1均布荷载/(N·m-2)5轴力/N125弯矩/(N·m)625
图3 拼装后的盾构隧道试验模型
1.3.1试验装置
试验采用三环拼装结构,如图4所示.卧式加载,可实现12个方向荷载的独立施加,每个方向千斤顶均采用“一拖三”的形式,保持压力一致.通过地层弹簧及荷载分配梁(曲板)以均布荷载的形式施加于管片模型结构.
a 实物图
b 示意图
1.3.2荷载施加过程
依据埋深及地下水位信息,采用水土合算法按照全覆土柱理论计算原型和模型隧道的外部水土压力,地层参数如表2所示,水土压力如表3所示.
表2 地层参数
试验加载方式为静力全周加载[11],通过12个方向的弹簧-弧形钢板装置来模拟均布荷载.均布荷载分成3组,分别为q1、q2与q3,加载时完全同步,如图5所示.卸载工况按以下2个阶段进行:
表3 水土压力
模拟实际埋深的加载阶段:由零应力初始状态,分5步对称加载至实际埋深,每步荷载增量Δq1=7.20 kPa,Δq2=5.48 kPa,Δq3=5.52 kPa.
模拟基坑开挖的卸载阶段:保持q1不变,由实际埋深状态分步对称卸载,调节q2和q3,每步卸载增量Δq2=-3.91 kPa,Δq3=-1.96 kPa,直至q2=0 kPa或结构破坏.
图5 环向荷载分组示意图
1.3.3量测方案
试验过程中,为了能有效揭示盾构隧道管片衬砌结构渐进性破坏规律,测试管片应变和位移.各物理量测点布置如图6所示.
1.4.1管片变形及破坏过程
图7反映了侧向荷载与结构变形(中环管片)的对应关系.总体上,结构变形表现为“横鸭蛋”,即横向扩张、竖向收敛.加载阶段,变形近似线性增加,并且竖向收敛值大于横向扩张值.实际埋深条件下,中环管片竖向收敛7.60 mm,横向扩张4.20 mm,拱底内侧、拱顶内侧和拱腰外侧出现微裂缝,裂缝未张开,其他部位未见裂损现象.
图6 中环管片测点布置示意图
图7 中环管片断面收敛与侧向荷载关系(试验)
开始卸载后,由于水平约束的减小,横向扩张量非线性增加,快速向竖向收敛值逼近,如图7所示.卸荷比超过72%后,横向扩张与竖向收敛接近同步发展;侧向荷载降低为零时,中环管片竖向收敛22.20 mm,横向扩张21.46 mm.随着变形的快速发展,拱底内侧、拱腰外侧、拱顶内侧相继出现新的裂缝,拱顶、拱底最为密集;横向扩张与竖向收敛一致时,上环与中环出现明显的错台,导致上环拱部保护层开裂,并沿环向迅速贯通.中环管片裂缝展布如图8所示,隧道最终破坏形态如图9所示.
图8 中环管片裂缝展布
图9 隧道最终破坏形态
卸载工况下管片结构破坏过程分为以下几个阶段:①设计荷载作用下管片受力状态→②中环拱底内侧开裂→③中环左拱腰外侧开裂→④上、中环拱顶错台加深,螺栓变形增加→⑤上环拱顶保护层裂缝贯通,结构变形过大并失去承载性能.破坏过程与中环管片断面收敛关系如表4所示.
表4 破坏过程与中环管片断面收敛关系
1.4.2管片应力
按应力-应变关系,通过公式σm=εEm(σm为混凝土应力,ε为混凝土应变,Em为混凝土弹性模量)将隧道模型混凝土应变值转换为应力,中环管片混凝土拉应力随侧向荷载变化如图10所示.
由图10可知,在设计荷载条件下,中环管片拱顶内侧、拱底内侧和拱腰外侧的混凝土应力较小,约为0.5 MPa,超过其抗拉强度,与微裂缝产生现象相符.随着侧向荷载的不断减小,中环拱底内侧应力增长较为迅速,为管片受拉最不利位置,混凝土拉应力为0.89 MPa时产生一批新裂缝.两侧卸载作用下,中环左、右拱腰外侧混凝土拉应力增长趋势基本一致,侧向卸载29%时,左拱腰外侧产生一批新裂缝,混凝土拉应力为0.86 MPa.管片变形所产生的内力除了自身承受之外,还通过接头向两侧传递,由于顶部封顶块接缝距拱顶较近,有利于释放管片自身内力,拱顶内侧混凝土拉应力相较拱顶增长较为缓慢.
图10 中环管片混凝土拉应力随侧向荷载变化
以模型试验为对象,建立基于荷载-结构法的三维有限元计算模型,分析管片应力、螺栓应力、接缝张开等随卸荷的发展过程,从细观角度揭示淤泥层错缝45°角拼装的盾构隧道在双侧卸载条件下隧道结构的承载性能演化规律.
采用六面体实体单元模拟衬砌,采用杆单元模拟钢筋与接头螺栓,如图11所示.管片衬砌采用混凝土损伤塑性模型,允许结构产生裂缝,最大裂缝允许值为0.15 mm,混凝土裂缝宽度按照《混凝土结构设计规范》[12]提供的公式对裂缝宽度进行验算.螺栓与钢筋采用理想弹塑性模型模拟,当应力达到屈服强度后钢筋、螺栓的承载性能不再增长.混凝土、钢筋及螺栓的物理力学参数如表5和表6所示.
图11 隧道模型
管片之间设置摩擦接触,并通过螺栓连接.管片结构与地层的相关作用通过设置地层弹簧模拟,同时设置剪切弹簧以约束整环的旋转运动.荷载条件、加载方式和卸荷方式均与模型试验保持一致,考虑重力方向与隧道纵向平行,底面设置约束,顶面为自由端.
表5 混凝土材料参数
表6 螺栓和钢筋材料参数
图12反映了侧向荷载与结构变形(中环管片)的对应关系,结构总体呈现“横鸭蛋”变形模式.由初始状态加载至实际埋深的过程中,管片变形表现为线性增加.实际埋深下,产生的竖向收敛为5.41 mm,横向扩张为5.01 mm,左拱腰接缝张开0.12 mm.管片拱底内侧混凝土拉应力最大,为2.213 MPa,拱顶螺栓应力最大,为37.1 MPa,钢筋拉应力最大位置在拱底内侧,为10.5 MPa.设计荷载作用下结构内力如图13所示.
图12 中环管片断面收敛与侧向荷载关系(数值计算)
Fig.12 Relation between cross section convergence and lateral load of middle ring segment (numerical calculation)
开始卸载后,由于淤泥地层地基抗力小、侧压力系数较大,侧向约束逐渐削弱,拱底内侧和左拱腰外侧混凝土拉应力相继达到抗拉强度2.74 MPa,并产生微裂缝.计算中以钢筋拉应力验算裂缝宽度,侧向卸载31%时,拱底内侧钢筋应力达136 MPa,混凝土裂缝宽度首先达设计允许最大值0.15 mm,随后在侧向卸载62%时,左拱腰外侧混凝土裂缝宽度达0.15 mm.由于三环错缝拼装在隧道纵向的不对称性,因此左侧接缝张开量大于右侧,侧向卸载77%时,左拱腰78.75°接缝向外张开5.34 mm,环向螺栓首先屈服.拱底A2块收敛大于顶部,并且两侧接缝处于压紧状态,使得底部刚度较弱,钢筋最不利位置位于拱底内侧.侧向卸载100%时,产生的竖向收敛为103.17 mm,横向扩张为106.50 mm,11.25°接缝张开5.79 mm,78.75°接缝张开8.56 mm,拱底内侧钢筋应力为346 MPa.卸载工况下管片达到极限承载状态时变形及混凝土最大主应力如图14所示.
图13 设计荷载作用下中环管片混凝土最大主应力云图
图14 双侧卸载工况下中环管片达到极限承载状态时混凝土最大主应力云图
Fig.14 Maximum principal stress nephogram of concrete under bilateral unloading when reaching ultimate bearing state
淤泥地层双侧卸载作用下管片损伤过程为:ⓐ实际埋深下管片变形受力状态→ⓑ拱底内侧裂缝宽度达到设计允许值0.15 mm→ⓒ左拱腰外侧裂缝宽度达到设计允许值0.15 mm→ⓓ左拱腰78.75°环向螺栓屈服→ⓔ拱顶11.25°环向螺栓屈服→ⓕ结构变形过大并失去承载性能.管片卸载损伤过程承载临界状态对应断面变形如表7所示,承载状态与竖向断面变形关系如图15所示.
表7 双侧卸载工况下达到不同承载临界状态的管片变形
图15 承载状态发展与竖向断面变形关系
为了研究卸载工况下淤泥地层盾构隧道结构变形与受力的对应关系,结合数值计算与模型试验的结果,用中环管片竖向收敛与管片隧道直径的比值(ΔD/D)来表示隧道关键节点的中环竖向断面收敛率.卸载工况下隧道断面收敛率随侧向荷载变化如图16所示(损伤过程编号与第1.4.1节和第2.2节相同),图中无下划线文字描述代表模型试验,有下划线文字描述代表数值计算.
由图16对比分析可知,2种研究手段存在一定差异,主要表现在:模型试验前期变形较大,这是由于存在拼装误差,因此在前期加卸载阶段接缝逐渐闭合、管片压紧,而且由于对接头处的传力橡胶进行简化,因此管片整体性较差,并且块与块之间约束力较弱,使得前中期结构变形偏大;模型试验中管片开裂较早,这是由于采用的混凝土相似材料为砂浆,因此抗拉性能较弱,易在较低应力状态下开裂.
图16 竖向收敛率随侧向荷载变化
上述差异主要限于试验条件,但仍有较为一致的现象,主要为:管片混凝土损伤的主要部位为拱顶内侧和左拱腰外侧;管片破坏的最终形式为结构变形过大,并且两者变形结果接近.这些相同点在一定程度上验证了数值计算结果的合理性,考虑数值计算更为严谨,数据分析精细化程度更高,以计算结果为主要依据分析结构破损过程,提炼变形控制标准,更加有利于运营安全.
综上所述,以数值计算结果为依据,得出不同侧向荷载与顶、底荷载比值(q2/q1)下结构的变形和破损过程,主要为以下4个阶段:
(1)弹性阶段(q2/q1=0.72),即初始状态至实际埋深阶段.管片呈“横鸭蛋”变形模式,变形呈线性增加,结构无裂损发生.在实际埋深条件下,接缝张开最大位置为左拱腰,拱底内侧为混凝土和钢筋受拉最不利位置.
(2)裂缝产生及发展阶段(q2/q1=0.72~0.50),即实际埋深至混凝土裂缝宽度达设计允许值0.15 mm阶段.中环管片在拱底内侧和左拱腰外侧混凝土抗拉强度先后达2.74 MPa并产生微裂缝.随着侧向卸载的增加,拱底内侧和左拱腰外侧钢筋拉应力相继达136 MPa,裂缝宽度增加至设计允许值0.15 mm.
(3)错台发展和接缝张开阶段(q2/q1=0.50~0.17),即裂缝宽度达设计允许最大值至螺栓屈服阶段.裂缝的发展直接削弱了管片整体刚度,变形加剧,接缝张开量增大.当左拱腰78.75°接缝张开5.34 mm时,该位置环向螺栓首先屈服,随后拱顶11.25°环向螺栓在接缝张开3.50 mm时屈服.
(4)破坏阶段(q2/q1=0.17~0),即螺栓屈服至变形过大并失去承载性能阶段.螺栓屈服后,管片承载力和整体刚度显著降低,当侧向卸载100%时,结构变形过大并失去承载性能.
从盾构隧道运维的角度出发,根据卸载工况下管片损伤演变过程,提出基坑开挖施工时的变形控制指标建议值,如表8所示.
表8 管片结构损伤节点与断面收敛关系
(1)淤泥地层地铁盾构隧道在1.5D埋深条件下,侧向压力与顶部荷载比值为0.72~0.50时,管片结构处于裂缝产生及发展阶段,断面收敛增长速率及构件承载性能降低程度对双侧卸载的敏感度低;侧向压力与顶部荷载比值从0.50减小至0的过程中,管片承载状态从错台发展和接缝张开阶段发展至破坏阶段,断面收敛增长速率及构件承载性能降低程度对双侧卸载的敏感度高.因此,软土地层条件下,在邻近地铁双侧基坑开挖工程中,从管片变形控制及构件承载性能角度考虑,建议将施工中的隧道双侧荷载控制在0.5倍顶、底荷载以上,以保证隧道的安全运营.
(2)淤泥地层错缝拼装盾构隧道在双侧卸载工况下,管片拱底内侧混凝土裂缝宽度达到设计允许值0.15 mm的收敛变形值为3.29‰D;管片左拱腰外侧混凝土裂缝宽度达到设计允许值0.15 mm的收敛变形值为7.67‰D;管片左拱腰78.75°环向螺栓屈服时收敛变形值为10.67‰D;管片拱顶11.25°环向螺栓屈服时收敛变形值为11.21‰D.结构破坏形式表现为:拱底内侧裂缝宽度达到设计允许值0.15 mm→左拱腰外侧裂缝宽度达到设计允许值0.15 mm→左拱腰78.75°环向螺栓屈服→拱顶11.25°环向螺栓屈服→变形过大并垮塌.