中墙顶部回填混凝土弹性模量对连拱隧道稳定性的影响

陈赟，马云锋，唐前松，周勇

(1.长沙理工大学 交通运输工程学院，湖南 长沙，410004；2.湖南省高速公路管理局，湖南 长沙，410001)

连拱隧道和分离式隧道是目前修建公路隧道的 2种主要形式。由于连拱隧道线路占地面积小，在公路隧道中的应用越来越多。但是，连拱隧道跨度较大，其开挖和支护相对于分离式隧道要复杂得多，因此，对施工技术提出了更高的要求[1]。中墙是连拱隧道衬砌体系中的关键部位，一方面，要支承中导洞上部的围岩压力，与围岩一起组成拱形受力体系，发挥围岩的自承作用；另一方面，它又是内衬的支点，要承受内衬通过耳墙传来的压力。因此，中墙的受力状态非常复杂，其稳定性已成为工程技术人员关注的焦点。夏才初等[2]利用有限元方法计算了连拱隧道的中墙应力，并与实测结果进行对比，认为在计算中墙应力时，必须考虑中墙顶部围岩由于扰动所产生的附加围岩压力。时亚昕等[3]认为围岩类别对中墙的力学行为影响显著。刘新荣等[4]对黄土地区连拱隧道中墙受力情况进行了分析，认为偏压有可能使中墙产生不均匀沉降和扭转等问题。曹云钦等[5−6]对不同的中墙结构形式进行了计算分析，提出复合曲中墙结构有利于改善中墙的受力状况。曾胜等[7−8]研究了施工工序对连拱隧道中墙的影响。由此可见，中墙的受力情况十分复杂，并对连拱隧道稳定性产生重要影响。为了使中墙有效发挥作用，施工中通常利用回填混凝土密实填充中墙顶部与围岩之间因超挖产生的空隙，有效传递来自中墙上部围岩的压力。然而，因受到实际施工过程中施工水平及其他多种因素的影响，无法保证中墙顶部与围岩之间的空隙得到完全密实填充，这必然对连拱隧道的结构受力体系产生影响。中墙顶部回填混凝土填充不密实有多种表现，其中最直接的表现是对围岩支撑刚度减弱，在宏观力学行为上体现为回填混凝土材料弹性模量降低。在此，本文作者针对连拱隧道中墙顶部回填混凝土不同的密实度情况，采用有限元方法模拟计算回填混凝土不同弹性模量下隧道围岩的变形与衬砌结构的受力情况，综合研究中墙顶部回填混凝土对连拱隧道稳定性的影响。

1 工程概况

某高速公路连拱隧道全长300 m，隧道进出口处有3～8 m残坡积层，其余地段为花岗岩及其风化层直接出露，围岩等级为Ⅳ级；洞身段以弱～微风化混合花岗岩为主，有少部分凝灰熔岩，局部辉长(绿)岩脉岩侵入，围岩等级为Ⅱ～Ⅳ级，致密而坚硬。整个隧道按新奥法进行设计和施工。连拱隧道的双洞不对称，分成大洞(右洞)和小洞(左洞)。采用复合衬砌，以锚杆湿喷混凝土、钢筋网等为初期支护，并辅以钢支撑、注浆小导管等支护措施，充分调动和发挥围岩的自承能力，在监控量测信息的指导下施筑二次衬砌。该连拱隧道衬砌结构形式如图1所示，结构参数见表1。

图1 某高速公路连拱隧道衬砌结构形式Fig.1 Lining structure of expressway tunnel

表1 隧道结构参数Table 1 Parameters of lining structure mm

2 计算模型

2.1 有限元模型建立

在有限元模型中，岩体和中墙采用4节点的平面四边形单元模拟，局部采用 3节点的三角形单元模拟；锚杆采用两维杆单元模拟；喷层、二衬(内衬)和仰拱采用两维梁单元进行模拟。同时，考虑到内外衬之间有可能错动，在两者之间设立了接触面单元，允许二衬(内衬)和初衬(外衬)之间有一定的滑动。有限元模型的建模范围取隧道跨度的6倍以上，以减少边界效应。模型顶面无约束，以模拟地面变形；左右两侧水平位移约束，竖向自由；底面同时约束水平和竖向位移。由于该隧道围岩主要为坚硬花岗岩，结构完整，呈弹性性质，故本研究中仅进行线弹性分析，并通过自重力计算得到初始应力场。有限元网格剖分如图2所示，材料计算参数见表2。

图2 有限元网格Fig.2 FEM mesh

表2 围岩计算参数Table 2 Parameters of surrounding rock

2.2 施工步骤模拟

为考察实际施工过程中连拱隧道的稳定性，根据实际施工工序对连拱隧道的全施工过程进行数值模拟。连拱隧道的施工工序为：(0)初始状态→(1)中导洞开挖→(2)中导洞锚喷支护→(3)施作中墙→(4)两侧导洞开挖→(5)两侧导洞锚喷支护→(6)左洞开挖→(7)左洞锚喷支护→(8)左洞施作二衬→(9)左洞开挖仰拱→(10)左洞施作仰拱→(11)右洞开挖→(12)右洞锚喷支护→(13)右洞施作二衬→(14)右洞开挖仰拱→(15)右洞施作仰拱。由此得到模拟施工的关键步骤见图3。

图3 关键施工步骤Fig.3 Key procedures of construction

2.3 不同密实度的模拟

图4 中墙顶部的回填混凝土单元Fig.4 Elements of back-fill concrete above mid-wall

中墙顶部回填混凝土密实度不同，在宏观力学表现上主要表现为弹性模量不同。因此，有限元计算中可通过改变中墙顶部回填混凝土单元材料的弹性模量来模拟不同的密实程度，如图4所示。在计算分析中，假设6种不同的密实度，其对应的中墙顶部的回填混凝土的弹性模量分别为 C25标号混凝土弹性模量(20 GPa)的0%(不回填，完全脱离)，10%，25%，50%，75%，100%(完全密实)。

2.4 计算指标的确定

对连拱隧道稳定性的影响因素很多，如围岩变形、围岩应力分布、衬砌结构的内力等[9−16]，要寻求单一反映连拱隧道稳定性的指标非常困难。通过综合分析，在连拱隧道安全性评价中，隧道围岩的拱顶沉降，拱顶围岩拉应力和内衬弯矩最为关键；此外，偏压对中墙稳定性也有重要影响[4−6]。因此，本研究中，选取反映隧道围岩稳定性的指标即拱顶沉降和拱顶处最大拉应力、反映衬砌安全的指标即内衬拱顶最大正弯矩(假定隧道内侧受拉的弯矩为正弯矩)、反映偏压情况的指标即中墙左侧和右侧的最大压应力进行计算和比较分析。

3 结果与讨论

3.1 计算结果

在不同弹性模量下，分别对拱顶沉降和拱顶处最大拉应力、内衬拱顶最大正弯矩、中墙左侧和右侧的最大压应力进行计算，结果见表3。

表3 回填混凝土弹性模量不同时的计算结果Table 3 Calculated results under different elastic modulus of back-fill concrete

3.2 围岩稳定性

在回填混凝土不同弹性模量条件下，左洞和右洞的拱顶沉降及拱顶受拉区的围岩最大拉应力分别见图5和图6。从图5和图6可以看出：左右洞的拱顶沉降和围岩拉应力均随着回填混凝土弹性模量的增加而逐渐降低，但降低幅度逐渐趋缓；弹性模量比(即回填混凝土弹性模量与混凝土弹性模量之比)为100%时左、右洞的拱顶沉降约为弹性模量比为0时的75%左右；拱顶围岩拉应力的降幅更为明显，100%弹性模量时左洞的拉应力仅为弹性模量比为0时的20%左右。这说明回填混凝土弹性模量对连拱隧道围岩变形和受力的影响非常显著。当弹性模量比为0即中墙顶部与围岩完全脱离时，中墙上部的围岩失去了中墙的支撑，导致整个围岩系统无法形成一个良好的拱形受力体系以发挥围岩的自承作用，这对于围岩的稳定非常不利，在施工中应当尽量避免。

此外，从图5和图6还可以看到：右洞的拱顶沉降和拱顶拉应力均比左洞的大。这是因为所模拟的非对称连拱隧道，右洞跨度大于左洞跨度，因此，右洞开挖时释放的围岩压力荷载大于左洞压力荷载，产生的围岩变形和应力也相应地比左洞的高。

3.3 衬砌安全性

图5 中墙顶部回填混凝土在不同弹性模量时的拱顶沉降Fig.5 Settlement of arch crown under different elastic modulus

图6 中墙顶部回填混凝土在不同弹性模量时的拱顶最大拉应力Fig.6 Maximum tensile stress of surrounding rock at arch crown under different elastic modulus

图7 中墙顶部回填混凝土在不同弹性模量时内衬的最大正弯矩Fig.7 Maximum bending moment of inner lining at arch crown under different elastic modulus

回填混凝土弹性模量变化对内衬最大弯矩的影响见图7。从图7可以看出：左、右洞的内衬最大弯矩受回填混凝土弹性模量的影响也非常显著，内衬弯矩随着弹性模量的增加而急剧降低；弹性模量比为100%时，左、右洞的内衬弯矩仅为弹性模量比为0时的17%左右。根据前面的分析，当弹性模量比为0时，中墙顶部与围岩完全脱离，围岩的自承能力较弱，围岩荷载集中作用在衬砌结构上，导致内衬弯矩增大，这对于衬砌结构的安全非常不利。另外，当弹性模量比达到60%以上时，内衬弯矩已基本保持不变。这表明：从隧道安全的角度考虑，应当使中墙顶部回填混凝土的弹性模量至少达到混凝土弹性模量标准值的60%以上才进行下一步施工。

3.4 偏压分析

回填混凝土弹性模量对中墙左右侧最大压应力的影响见图8。从图8可以看出：中墙最大压应力随回填混凝土弹性模量的增加而增加，且增加的幅度也随着弹性模量的增加而趋缓。这说明随着弹性模量的增加，中墙上方围岩通过与回填混凝土接触传递到中墙的压力也将增大；当密实度为0即中墙顶部与围岩脱离时，中墙没有受到上方围岩传来的压力，只承受耳墙传来内衬的压力，此时，中墙两侧最大压应力仅为回填混凝土完全密实时的33%左右。

图8 中墙顶部回填混凝土在不同弹性模量时中墙两侧最大压应力Fig.8 Maximum pressure in two sides of mid-wall under different elastic modulus

4 结论

(1)围岩的变形和应力随着中墙顶部回填混凝土弹性模量的增加而降低，当回填混凝土弹性模量比由 0(完全脱离)增加到 100%(完全密实)时，拱顶沉降和围岩拉应力分别降至75%和20%左右。

(2)衬砌最大弯矩随着中墙顶部回填混凝土弹性模量的增加而大幅降低；当回填混凝土弹性模量比由0(完全脱离)增加到100%(完全密实)时，内衬弯矩降至17%左右，当弹性模量比增加到60%以上时内衬弯矩基本保持不变。

(3)中墙应力随着回填混凝土弹性模量的增加而增加，弹性模量越低，中墙应力变化幅度越为显著；回填混凝土完全脱离时，中墙两侧最大压应力仅为回填混凝土完全密实时的33%左右。

(4)在实际施工中，应使中墙回填混凝土的弹性模量至少达到设计标号混凝土标准值的60%以上时方可进行隧道正洞开挖。

[1]陈少华,李勇.联拱隧道的结构分析[J].中国公路学报,2000,13(1): 48−51.CHEN Shao-hua,LI Yong.Structural analysis for a joined roof tunnel[J].China Journal of Highway and Transport,2000,13(1):48−51.

[2]夏才初,刘金磊.相思岭连拱隧道中墙应力研究[J].岩石力学与工程学报,2000,19(S): 1115−1119.XIA Chai-chu,LIU Jin-lei.Study on the middle wall stress of Xiangsilin doubled arch tunnel[J].Chinese Journal of Rock Mechanics and Engineering,2000,19(S): 1115−1119.

[3]时亚昕,王明年,李强.浅埋双连拱隧道中墙的合理选用[J].地下空间与工程学报,2005,1(6): 948−951.SHI Ya-xin,WANG Ming-nian,LI Qiang.Proper selection of the middle wall of a shallow double-linked-arch tunnel[J].Chinese Journal of Underground Space and Engineering,2005,1(6):948−951.

[4]刘新荣,孙辉,陈晓江,等.黄土连拱隧道中墙结构的数值模拟研究[J].地下空间与工程学报,2005,1(6): 837−840.LIU Xin-rong,SUN Hui,CHEN Xiao-jiang,et al.Numerical analysis on mid-partition structure of the loess double arch tunnel[J].Chinese Journal of Underground Space and Engineering,2005,1(6): 837−840.

[5]曹云钦,王小林.浅埋偏压连拱隧道中墙优化分析[J].岩土工程学报,2006,28(4): 537−540.CAO Yun-qin,WANG Xiao-lin.Optimal analysis of mid-board of shallow double-arch tunnels under unsymmetrical pressure[J].Chinese Journal of Geotechnical Engineering,2006,28(4):537−540.

[6]李强,王明年,李玉文.双跨连拱隧道两种中墙的空间力学效应分析[J].岩土力学,2006,27(4): 667−672.LI Qiang,WANG Ming-nian,LI Yu-wen.Three-dimensional mechanical effect analysis of two types of mid-board for double span-arch-connecting tunnel[J].Rock and Soil Mechanics,2006,27(4): 667−672.

[7]曾胜,覃庆通,阳军生.大跨度连拱隧道直中墙受力现场监测分析[J].岩土力学,2007,29(9): 2537−2541.ZENG Sheng,QIN Qing-tong,YANG Jun-sheng.In-situ monitoring behaviors of middle wall of large-span multi-arch tunnel[J].Rock and Soil Mechanics,2007,29(9): 2537−2541.

[8]王胜辉,杜小平,袁勇,等.不同开挖工序对连拱隧道中墙的影响[J].岩土力学,2004,25(7): 1102−1107.WANG Sheng-hui,DU Xiao-ping,YUAN Yong,et al.Effects of different excavating sequences on mid-partition of double-arched tunnel[J].Rock and Soil Mechanics,2004,25(7): 1102−1107.

[9]Whittaker B N,Reddish D J.Evaluation of tunnel stability in relation to rock stress and displacement aspects[J].International Journal of Rock Mechanics and Mining Sciences &Geomechanics Abstracts,1991,28(5): 177-185.

[10]Cai M.Influence of stress path on tunnel excavation response–numerical tool selection and modeling strategy[J].Tunnelling and Underground Space Technology,2008,23(6): 618−628.

[11]Nazimko V V,Peng S S,Lapteev A A,et al.Damage mechanics around a tunnel due to incremental ground pressure[J].International Journal of Rock Mechanics and Mining Sciences &Geomechanics Abstracts,1997,34(3/4): 1−14.

[12]杨小礼,李亮,刘宝琛.偏压隧道结构稳定性评价的信息优化分析[J].岩石力学与工程学报,2002,21(4): 484−488.YANG Xiao-li,LI Liang,LIU Bao-chen.Evaluation on structure stability of unsymmetrically loaded tunnels using the theory of information optimization analysis[J].Chinese Journal of Rock Mechanics and Engineering,2002,21(4): 484−488.

[13]靳晓光,李晓红,杨春和,等.深埋隧道围岩−支护结构稳定性研究[J].岩土力学,2005,26(9): 1473−1476.JIN Xiao-guang,LI Xiao-hong,YANG Chun-he,et al.Stability of surrounding rock-supports structure of deep buried tunnel[J].Rock and Soil Mechanics,2005,26(9): 1473−1476.

[14]许崇帮,夏才初,王华牢.特大断面连拱隧道中墙偏压机制及施工影响分析[J].岩石力学与工程学报,2010,29(S1):2819−2826.XU Chong-bang,XIA Cai-chu,WANG Hua-lao.Analysis of mechanism of mid-partition eccentric compression and construction effect on super-large section multi-arch tunnel[J].Chinese Journal of Rock Mechanics and Engineering,2010,29(S1): 2819−2826.

[15]杨小礼,眭志荣.浅埋小净距偏压隧道施工工序的数值分析[J].中南大学学报: 自然科学版,2007,38(4): 764−770.YANG Xiao-li,SUI Zhi-rong.Numerical simulation of construction sequence for shallow embedded bias tunnels with small clear distance[J].Journal of Central South University:Science and Technology,2007,38(4): 764−770.

[16]LI Xi-bin,ZHANG Wei,LI Di-yuan,et al.Influence of underground water seepage flow on surrounding rock deformation of multi-arch tunnel[J].Journal of Central South University of Technology,2008,15(1): 69−74.