结合重庆同茂大道金竹湾市政交通隧道与电力隧道共建工程,在理论分析的基础上,针对提出的两种施工方案,建立二维有限元模型,对其交通–电力隧道共建方案下围岩的塑性区半径、弹塑性交界处的应力以及弹塑性位移进行对比分析,提出较佳的共建方案。 Based on the theoretical analysis of the Jinzhuwan municipal traffic tunnel and power tunnel co-construction project of Tongmao Avenue in Chongqing, a two-dimensional finite element model is established for the proposed two construction schemes, and the traffic-electricity tunnel co-construction scheme is selected. The plastic zone radius of the rock, the stress at the elastoplastic interface and the elastoplastic displacement are compared and analyzed, and a better co-construction scheme is proposed.
史宁强1,葛婷婷2,靳晓光3*
1林同棪国际工程咨询(中国)有限公司,重庆
2华南理工大学广州学院,广东 广州
3重庆大学土木工程学院,重庆
收稿日期:2019年4月29日;录用日期:2019年5月14日;发布日期:2019年5月21日
结合重庆同茂大道金竹湾市政交通隧道与电力隧道共建工程,在理论分析的基础上,针对提出的两种施工方案,建立二维有限元模型,对其交通–电力隧道共建方案下围岩的塑性区半径、弹塑性交界处的应力以及弹塑性位移进行对比分析,提出较佳的共建方案。
关键词 :共建隧道,方案优化,理论分析,数值模拟
Copyright © 2019 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
随着山地城市交通隧道和其他公用设施隧道的快速发展,传统的架空电力线路影响城市的空中景观及地块使用。此前国内的高压电缆隧道一般采用独立建设的模式,不仅要新开辟通道,而且建设代价高昂。将电力隧道和交通隧道合二为一,不仅节约土地和成本,而且可以共用一个监控系统。因此市政隧道和电力隧道共建的设计模式近年来在我国开始得到有效的推广、应用。目前关于共建隧道的研究方法和实际的施工经验还相对匮乏,关于共建隧道的研究主要借鉴单孔大断面隧道及双孔近接隧道的研究成果。
张晓军 [
金竹湾1#隧道左线与电力隧道共建,经综合分析提出两种方案,一是电力隧道建在城市交通隧道仰拱的正下方(简称方案一);二是电力隧道建在隧道仰拱的左侧(简称方案二)。如图1所示。市政隧道与电力隧道共建后,共建段隧道路面以上限界及净空尺寸不变,仅需对仰拱进行加深,非共建段隧道仰拱半径为21 m。方案一仰拱开挖较非共建段隧道加深1.8 m (开挖深度达到3.5 m),仰拱半径为8.89 m;方案二仰拱开挖较非共建段隧道加深2.8 m,仰拱半径为1.1 m。
在隧道围岩变形与破坏的简化理论分析中,常把边界形状复杂的洞室,如直墙拱形、曲墙拱形等接近圆形断面的洞室变换为形状简单的洞室形状,即圆形洞室,这种方法称为等代圆法。将非圆形的大跨度隧道等代为圆形隧道,一般取等代圆半径为隧道高度h和跨度b之和的四分之一。两方案隧道的几何尺寸分别如图2所示。
图1. 交通–电力共建隧道示意图
图2. 交通–电力共建隧道几何尺寸(单位:cm)
共建方案只在隧道左线,右线不变。右线隧道的h1 = 10.62 m,b1 = 15.29 m。方案一和方案二中共建隧道的高度和跨度分别为h2 = 12.42 m,b2 = 15.33 m和h3 = 13.39 m,b3 = 15.29 m。则右线隧道的等代圆半径为r1 = 6.4775 m;方案一中共建隧道的等代圆半径为r2 = 6.9375 m;方案二中共建隧道的等代圆半径为r3 = 7.17 m。如图3所示。
图3. 两方案等代圆几何尺寸
圆形洞室围岩的塑性区半径为:
R 0 = r 0 [ ( 1 − sin φ ) c cot φ + σ z c cot φ + p α ] 1 − sin φ 2 sin φ (1)
圆形洞室围岩弹性区与塑性区交界处的应力公式为:
{ σ r p = ( p α + c cot φ ) ( 1 − sin φ ) c cot φ + σ z c cot φ − c cot φ σ θ p = ( p α + c cot φ ) ( 1 + sin φ ) c cot φ + σ z c cot φ − c cot φ (2)
圆形洞室围岩中的弹塑性位移为:
u R 0 e = R 0 2 ( 1 + μ ) E r ( σ z − σ R 0 ) , r ≥ R 0 (3)
u p = R 0 2 ( 1 + μ ) E r ( σ z − σ R 0 ) , r 0 ≤ r ≤ R 0 (4)
围岩的塑性区半径和开挖半径成正比,开挖半径越大塑性区半径R0越大。以上三个隧道等代圆半径的关系为r3 > r2 > r1,则可得其塑性区半径的大小关系为R3 > R2 > R1。弹性区与塑性区交界处的应力只与围岩本身的物理力学性质c、φ及支护抗力有关,三个隧道的弹塑性交界处应力相同。围岩中任意一点的位移和洞室的塑性区半径成正比,塑性区半径越大,其径向位移越大,则有u3 > u2 > u1。
对隧道塑性区半径、弹塑性交界处的应力及围岩的弹塑性位移进行对比分析可以看出,方案一明显优于方案二。
针对以上两种方案,建立施工过程的二维有限元模型进一步进行分析,以确定较佳共建方案。模型采用理想弹塑性材料,屈服准则采用摩尔-库伦准则。方案一的数值模型有3701个单元,3556个节点;方案二的数值模型有3729个单元,3729个节点。所有单元类型采用四边形单元,整个模型范围为168 m × 99 m,两隧道的净矩为19 m。喷射混凝土采用梁单元。隧道的开挖和支护通过单元的激活和钝化功能来实现。
在有限元模型中,围岩和支护结构参数的选取依据隧道地勘报告,并结合《公路隧道设计规范》(JTGD70-2004)选取。具体的参数值见表1和表2。
岩土名称 | 重度(KN/m3) | 弹性模量(MPa) | 泊松比μ | 内聚力c (kPa) | 内摩擦角φ (˚) | 抗拉强度(kPa) |
---|---|---|---|---|---|---|
中等风化砂质泥岩 | 25.6 | 1080 | 0.34 | 370 | 32.17 | 150 |
中等风化砂岩 | 24.9 | 3030 | 0.30 | 1000 | 40.74 | 382 |
表1. 岩体力学参数表
名称 | 弹性模量(MPa) | 泊松比μ | 重度(kN/m3) | 备注 |
---|---|---|---|---|
喷混支护 | 30000 | 0.3 | 25 | 初期支护:C30厚0.22 m |
表2. 支护结构物理力学参数表
隧道施工采用新奥法,台阶法开挖,主要施工步序如图4所示。
图4. 共建隧道施工步序示意图
1) 围岩应力特征
隧道在开挖前,围岩处在应力平衡的稳定状态中。伴随着隧道的开挖,隧道围岩位于开挖边界的部分约束被解除,围岩失去了初始的应力平衡状态,发生应力状态的改变,并逐渐形成新的应力平衡状态。隧道开挖结束后,两方案围岩应力云图如图5~图6所示,可以看出,两方案围岩应力分布均在合理范围内。
图5. 隧道最大主应力云图
图6. 隧道最小主应力云图
由图5~图6可以看出:
①方案一和方案二的最大主应力均为压应力,其最大值差别不大。方案一中应力集中主要发生在右隧道的拱脚和左隧道的拱腰处,其值分别达到为2.74 MPa和2.47 MPa;方案二中应力集中主要发生在左隧道和右隧道的拱脚处,其值均达到2.79 MPa。
②方案一和方案二的最小主应力中均出现了拉应力,两方案中的拉应力均分布在围岩的上部和隧道的拱底处。方案一中拉应力的最大值发生在右隧道的拱底处,其值约为0.06 MPa,方案二中其最大拉应力出现在电力隧道和仰拱连接处,约为0.13 MPa,此处易发生拉裂破坏。
③隧道的最大和最小主应力,总体上表现出拱脚最大,其次是拱腰,再次是拱顶和拱底。
2) 支护结构特征
隧道支护结构轴力图如图7所示。
图7. 隧道支护结构轴力图
由图7可以得出:两方案中右隧道的支护结构顶部和两侧受压,底部受拉;左隧道支护结构底部有一部分受压。方案一中初期支护的最大压力出现在拱顶,约为898.9 kN,最大拉力出现在右隧道的仰拱处,约为1005.7 kN;方案二中最大压力出现在共建处,约为1189.3 kN,最大拉力出现在左隧道仰拱处,约为1388.8 kN。
3) 围岩位移特征
隧道的开挖引起对围岩的扰动,图8~图9分别为隧道的竖向位移和水平位移云图。
图8. 竖向位移变化云图
图9. 水平位移变化云图
从图8~图9可以看出,1) 两方案中左右隧道的拱顶均发生沉降,拱底均出现向上隆起。方案一及方案二中拱顶下沉值基本相同;但方案一中拱底隆起最大值约为1.22 cm,方案二中约为1.30 cm。2) 两方案中最大水平位移均出现在左隧道拱腰处,方案一中最大水平位移值为6.44 mm,方案二中其值为8.06 mm。
4) 围岩塑性区特征
两方案隧道施工完成后围岩塑性区的特征如图10所示。
图10. 围岩塑性区分布图
可以看出,方案一的围岩塑性区呈对称分布,右隧道塑性区主要分布在拱脚和拱腰处,左隧道塑性区主要分布在拱腰处,塑性应变最大值为6.07 × 10−5,围岩的塑性区发展较均匀;方案二的围岩塑性区呈不对称分布,右隧道与左隧道的塑性区均主要分布在拱腰和拱脚处,但共建隧道处没有塑性区分布,塑性应变最大值约为5.72 × 10−5。方案一和方案二的塑性应变整体较小,且其值相差不大。
对比隧道围岩的最大和最小主应力,方案一中最大主应力值小于方案二,且分布较均匀;方案二中在电力隧道和仰拱连接处出现了拉应力,此处易发生拉裂破坏;从两方案锚杆来看,方案一和方案二相差不大;从位移来看,两方案中拱顶下沉位移基本相同,但方案二拱底隆起位移和水平位移的最大值均比方案一大;两方案中塑性区整体应变均较小,且其值相差不大,两种方案没有明显区别。
综合以上对方案一和方案二支护结构的应力,围岩的应力、位移及塑性区的分布特征的对比分析,可以得到共建方案一优于方案二,即电力隧道建于左隧道仰拱的正下方的方案较优。
本文以重庆同茂大道金竹湾交通隧道与电力隧道共建工程为依托,通过理论计算和数值模拟,对电力隧道建在城市交通隧道左隧道仰拱的正下方及仰拱的左侧两种方案进行了系统地对比分析,得出较佳的设计施工方案,即电力隧道建在城市交通隧道仰拱的正下方较好。
本文研究获林同棪国际工程咨询(中国)有限公司科技计划项目资助。
史宁强,葛婷婷,靳晓光. 交通–电力暗挖隧道共建方案优化研究 Optimization on Co-Construction Scheme of the Traffic and Power Tunnel[J]. 土木工程, 2019, 08(03): 710-719. https://doi.org/10.12677/HJCE.2019.83083