扩展功能
文章信息
- 卢晓颖
- LU Xiao-ying
- 考虑时间效应的深埋小净距隧道围岩压力分析
- Analysis on Surrounding Rock Pressure on Deep Buried Neighborhood Tunnel Considering Time-dependency Effect
- 公路交通科技, 2021, 38(10): 100-106, 143
- Journal of Highway and Transportation Research and Denelopment, 2021, 38(10): 100-106, 143
- 10.3969/j.issn.1002-0268.2021.10.013
-
文章历史
- 收稿日期: 2020-12-08
囿于地形地貌、地质条件、选线要求等条件的限制, 公路工程中常采用小净距形式的隧道布置方案。而小净距隧道的围岩压力分布、开挖工法、支护技术等问题, 也是学界与业界关注的热点[1-2]。
李鹏飞等[3]基于普氏平衡拱理论, 提出了深埋非对称小净距隧道围岩压力分析模型, 并推导了深埋情况的围岩压力计算公式。李然等[4]考虑中岩柱承载作用, 提出了3孔小净距隧道围岩压力计算方法。孙振宇等[5]通过对小净距隧道围岩压力进行统计分析, 并基于隧道实际滑裂面提出了小净距隧道围岩压力的计算方法。李磊等[6]采用模型试验与数值模拟相结合, 对小净距隧道围岩变形规律进行研究, 结果表明后行洞对先行洞的卸载扰动使先行洞承受偏压荷载。谭忠盛等[7]通过现场对比试验分析, 认为两隧道间距较小时, 后行洞采用CD法施工对先行洞的影响较小。张恒等[8]基于室内模型试验, 分析双洞小净距隧道施工中的隧道变形规律, 认为起拱线以上的围岩变形随双洞净距增大而减小。
另一方面, 大量工程实践表明围岩具有显著的时间效应, 许多学者通过数值模拟方法对隧道围岩压力及其时间效应展开了研究。彭超等[9]和卢向勇等[10]通过数值模拟方法, 分析隧道围岩压力随时间的变化规律, 并对比了经典弹塑性和考虑蠕变效应的拱顶沉降计算结果。钱文喜等[11]同时考虑围岩埋深及蠕变效应的影响, 认为二衬所承担的围岩蠕变压力随着埋深增加而减小。耿大新等[12]研究围岩蠕变作用对公路隧道衬砌裂缝发展影响, 发现围岩蠕变效应是导致衬砌裂缝发育的主要原因之一。徐国文等[13]建立锚杆-围岩复合体流变模型, 分析了围岩蠕变作用对裂损衬砌长期安全性的影响。任松等[14]考虑围岩蠕变特性, 对隧道施工全过程进行动态模拟分析, 发现拱腰位置受围岩蠕变效应的影响最大。Xu等[15]采用颗粒离散元方法分析了蠕变作用下隧道二衬的破坏过程, 认为破坏过程可分为弹性阶段、初始损伤阶段、极限阶段和失稳阶段。Zheng等[16]针对隧道泥质砂岩展开蠕变试验, 建立了非线性蠕变模型, 并将其嵌入FLAC3D平台中模拟围岩的蠕变效应。
上述关于围岩压力时间效应的研究大部分仅针对初支, 二衬作为安全储备, 其围岩压力分布的相关研究仍较为稀缺。本研究以龙兴岭隧道为背景, 在FLAC3D数值平台上建立精细化数值模型, 分析小净距隧道围岩压力分布模式。进一步考虑围岩蠕变的时间效应, 探讨小净距隧道复合式衬砌支护结构围岩压力及初支/二衬荷载分担比随时间的变化规律, 并与现场监测结果相互验证。
1 工程背景龙兴岭隧道位于福建省平潭综合试验区, 为4洞小净距隧道, 全长500 m, 分别为2个机动车主洞、2个非机动车人行辅洞。主洞采用曲墙三心圆内轮廓, 其净空尺寸为14.2 m×9.46 m(宽×高), 左右主洞内轮廓间的最小距离为7.31 m。辅洞采用曲墙单心圆内轮廓, 其净空尺寸为7.69 m×6.4 m(宽×高), 左侧(右侧)主辅洞内轮廓间的最小距离为9.66 m。
根据相关勘察设计文件[17], 隧址地貌属构造、剥蚀形成的残坡积台地, 覆盖层厚度整体较薄, 隧道进口、出口处风化层厚度较大。隧道洞身围岩以Ⅳ级为主, 以深埋段K0+900断面为例, 其地质横断面如图 1所示, 由上到下依次为坡积粉质黏土、全风化凝灰岩、砂土状强风化凝灰岩、碎块状强风化凝灰岩、微风化凝灰岩。主洞埋深为45.1 m(大于2.5倍等效荷载高度, 属于深埋隧道), 采用上下台阶法开挖, 主洞初支采用C25喷射混凝土+I18钢拱架, 厚0.26 m, 并辅以系统锚杆, 二衬为C35模筑钢筋混凝土, 厚0.55 m。辅洞埋深为46.2 m(亦属深埋隧道), 采用全断面开挖。辅洞初支采用C25喷射混凝土+钢格栅, 厚0.26 m, 并辅以系统锚杆, 二衬为C35模筑钢筋混凝土, 厚0.45 m。
![]() |
图 1 龙兴岭隧道地质横断面(K0+900)(unit: m) Fig. 1 Geological cross-section of Longxingling Tunnel (at K0+900) (单位: m) |
|
为探讨围岩压力随时间的变化规律, 在初支外侧和二衬外侧埋设土压力盒, 对其所承担的围岩压力进行持续观测。以深埋段K0+900断面为例, 在左主洞(右主洞)拱顶、拱肩和边墙处的型钢拱架外侧埋设5只土压力盒, 如图 2 (a)所示; 在相同位置处的二衬钢筋网外侧埋设5只土压力盒, 并浇注于二衬内, 如图 2 (b)所示。各土压力盒的测线集束汇总后, 从横向排水盲管引出至边沟, 做好线头保护, 在此后的1 a间持续观测初支与二衬所承担的围岩压力。
![]() |
图 2 土压力盒埋设 Fig. 2 Embedment of earth pressure cells |
|
2 深埋段小净距隧道的数值模拟 2.1 数值模型
以龙兴岭小净距隧道K0+900断面为原型, 在FLAC3D平台上建立数值模型, 其整体尺寸为292 m×126 m×50 m(宽度×高度×进深), 如图 3所示。地层和二衬采用4节点或6节点实体单元模拟, 共计419 919个实体单元, 142 273节点; 初支采用3节点Liner单元模拟, 共计24 410个单元, 13 143个节点。模型边界条件为: 模型顶面(地表面)取自由边界, 其侧面为法向位移约束边界, 底面为全约束边界。岩土体采用Mohr-Coulomb本构模型, 衬砌和仰拱回填采用线弹性本构模型, 其物性参数分别如表 1和表 2所示[17]。
![]() |
图 3 深埋大断面小净距隧道的数值模型 Fig. 3 Numerical model of deep buried large sectional neighborhood tunnel |
|
参数 | 坡积粉质黏土 | 全风化凝灰岩 | 强风化凝灰岩 | 强风化凝灰岩 | 微风化凝灰岩 |
厚度/m | 4 | 3.2 | 10.8 | 6 | 102 |
密度/(kg·m-3) | 1 800 | 1 800 | 1 900 | 2 000 | 2 200 |
弹性模量/Pa | 1.35e7 | 5.29e7 | 9e7 | 1.5e9 | 1.6e8 |
泊松比 | 0.35 | 0.324 | 0.3 | 0.3 | 0.258 |
黏聚力/Pa | 2.5e4 | 2e4 | 2e4 | 8e4 | 2.5e5 |
内摩擦角/(°) | 11.7 | 25 | 30 | 35 | 40 |
参数 | 初支(辅/主洞) | 二衬(辅洞) | 二衬(主洞) | 仰拱回填 |
厚度/m | 0.26 | 0.45 | 0.55 | 2.3 |
密度/(kg·m-3) | 2 500 | 2 500 | 2 500 | 2 500 |
弹性模量/Pa | 2.8e10 | 4.82e10 | 4.82e10 | 2.57e10 |
泊松比 | 0.2 | 0.34 | 0.34 | 0.28 |
2.2 双向接触的Liner单元
初支与围岩、初支与二衬的相互作用如图 4所示, 通过liner单元side1的link连接实现初支与围岩的相互作用, 通过liner单元side2的link连接实现初支与二衬的相互作用。Liner单元的接触参数取值如表 3所示。
![]() |
图 4 围岩、初支、二衬之间的相互作用示意图 Fig. 4 Schematic diagram of interaction among surrounding rock, first support and second lining |
|
参数 | Side1 | Side2 |
法向连接刚度/(N·m-1) | 1e11 | 1e11 |
切向连接刚度/(N·m-1) | 1e11 | 1e5 |
切向连接黏聚力/Pa | 1.5e4 | 1.5e4 |
切向连接残余黏聚力/Pa | 0 | 0 |
切向连接摩擦角/(°) | 15 | 15 |
2.3 隧道开挖支护过程与测点布置
隧道开挖支护过程的模拟大致分为以下3步: (1)隧道开挖前的初始地应力计算; (2)辅洞采用全断面法开挖, 每步开挖进尺2 m, 每步开挖后及时进行辅洞的初支施工, 辅洞贯通后进行仰拱开挖回填与二衬施工; (3)主洞采用上、下台阶法进行开挖, 其施工步序与辅洞相同。
数值模型中的测点布置如图 5所示。辅洞布设14处监测点, 主洞布设18处监测点, 以监测初支(liner单元side1)与二衬(liner单元side2)上的法向围岩压力和切向围岩压力。
![]() |
图 5 开挖工法及测点布置图 Fig. 5 Excavation method and layout of measuring points |
|
2.4 考虑时间效应的Burgers-Mohr本构模型
开挖支护全部施作完成后, 将微风化凝灰岩地层改为可考虑时间效应的Burgers-Mohr本构模型, 以考察围岩压力随时间的变化趋势。Burgers-Mohr本构模型内置于FLAC3D平台的CREEP模块中, 由Kelvin模型和Maxwell模型串联组成, 如图 6所示。
![]() |
图 6 Burgers-Mohr本构模型示意图 Fig. 6 Schematic diagram of Burgers-Mohr constitutive model |
|
根据前人关于凝灰岩流变特性的研究[18], Burgers-Mohr本构模型中与时间相关的参数取值如下: 开尔文剪切模量Ek和黏滞系数ηk分别为4.68 GPa和1.10×105 GPa·s, 麦克斯韦剪切模量Em和黏滞系数ηm分别为9.23×10-2 GPa和1.10×107 GPa。
3 深埋小净距隧道围岩压力分布 3.1 辅洞开挖后围岩压力分布左、右辅洞开挖及支护完成后, 读取初支与二衬法向压力, 绘制其围岩压力分布图(图 7)。需要说明的是, 切向围岩压力值较小, 对围岩压力分布影响较小, 故本研究不予考虑。
![]() |
图 7 初支与二衬的围岩压力分布(左右辅洞)(单位: kPa) Fig. 7 Distribution of surrounding rock pressure on first support and second lining (left and right auxiliary tunnels)(unit: kPa) |
|
由图 7可知, 辅洞开挖完成后, 左、右辅洞初支围岩压力关于设计中线近似呈对称分布。从拱顶至拱肩至边墙, 初支围岩压力呈现先增大后减小的趋势, 最大值出现在拱肩处, 约为250~270 kPa; 边墙以下拱脚至拱底呈逐渐减小趋势。左、右辅洞二衬围岩压力基本为0, 仅仰拱附近存在较小的由仰拱自重引起的围岩压力。
3.2 主洞开挖后围岩压力分布左、右主洞开挖及支护完成后, 同样绘制其围岩压力分布图(图 8)。初支围岩压力的分布形态与辅洞基本相同, 最大值出现在外侧拱肩附近, 约为430~470 kPa; 二衬(除仰拱附近外)围岩压力也基本为0。
![]() |
图 8 初支与二衬的围岩压力分布(辅洞与主洞)(单位: kPa) Fig. 8 Distribution of surrounding rock pressure on first support and second lining (auxiliary tunnels and main tunnels) (unit: kPa) |
|
另一方面, 受主洞(后行洞)开挖扰动的影响, 辅洞(先行洞)初支所承担的围岩压力减小, 二衬所承受的围岩压力增大(内侧尤为明显), 但二者之和基本不变。
4 考虑时间效应的围岩压力分析进一步地, 考虑围岩蠕变的时间效应, 分析小净距隧道支护结构围岩压力及其初支/二衬荷载分担比随时间的变化规律。
4.1 蠕变1 a后的围岩压力分布开挖支护全部施作完成后, 将微风化凝灰岩地层改为可考虑时间效应的Burgers-Mohr本构模型, 计算蠕变1 a后初支与二衬上的围岩压力分布, 如图 9 (a)~(d)所示。与图 8对比可知, 围岩蠕变1 a后, 围岩压力由初支向二衬转移, 初支承担的围岩压力减小, 二衬承担的围岩压力增大, 但二者之和基本保持不变。
![]() |
图 9 初支与二衬的围岩压力分布及荷载分担比(蠕变计算1 a) Fig. 9 Distribution and sharing ratio of surrounding rock pressure on first support and second lining (creep calculation of 1 year after) |
|
进一步地, 将二衬(初支)所承受的围岩压力占总围岩压力比值, 定义为二衬荷载分担比(初支荷载分担比), 绘制初支/二衬荷载分担比, 如9 (e)~ (h)所示。可知主洞上台阶部分(拱顶至边墙)的二衬荷载分担比约为30%~45%, 下台阶部分(边墙至拱脚)的二衬荷载分担比约为20%~35%, 仰拱部分的二衬荷载分担比在25%~40%之间。辅洞的初支/二衬荷载分担比, 其发展规律与主洞类似, 不再赘述。
对初支与二衬所承担的围岩压力展开为期1 a的持续观测, 绘制其围岩压力分布及其荷载分担比, 如图 10所示。可知现场监测所得围岩压力的分布形态与数值模拟结果相似, 围岩压力最大值同样出现在拱肩处, 二衬荷载承担比在30%~50%之间, 与数值模拟结果相符。另外需要说明的是, 土压力盒在隧道开挖支护之后进行埋设, 因此前期释放的围岩压力无法测得, 故现场实测围岩压力普遍小于数值模拟结果。
![]() |
图 10 初支与二衬的围岩压力分布及其荷载分担比(现场实测) Fig. 10 Distribution and sharing ratio of surrounding rock pressure on first support and second lining (field monitoring) |
|
4.2 围岩压力随时间的变化规律
按隧道设计使用年限, 计算蠕变100 a内初支和二衬承担围岩压力的经时变化。以左主洞为例, 绘制监测点1, 3, 6, 9, 11上的围岩压力变化规律如图 11所示。可知拱顶、拱肩处的初支围岩压力在蠕变1 a内随时间大幅度减小, 而1 a后则继续缓慢减小; 二衬围岩压力在蠕变1 a年内大幅度增长, 而1 a后则继续缓慢增加。两侧边墙处的围岩压力则基本不随时间发生变化。
![]() |
图 11 初支与二衬围岩压力的经时变化(左主洞) Fig. 11 Surrounding rock pressure on first support and second lining varying with time (left main tunnel) |
|
综上所述, 随计算蠕变时间的增长, 围岩压力逐渐由初支向二衬转移, 且二者之和基本保持不变。蠕变稳定后, 主洞上台阶部分(拱顶至拱腰)的二衬荷载分担比约为25%~45%, 主洞下台阶部分(仰供至边墙)的二衬荷载分担比约为20%~35%。
5 结论以平潭综合实验区龙兴岭隧道为工程背景, 在FLAC3D平台上对4洞小净距隧道的开挖支护过程展开细致的数值模拟, 重点关注作用在初支与二衬上的围岩压力分布及其随时间变化规律, 得到以下主要结论。
(1) 复合式衬砌支护结构围岩压力总体上关于设计中线近似呈对称分布, 但小净距隧道后行洞开挖扰动引起先行洞二衬上围岩压力增大。
(2) 考虑围岩、初支、二衬之间的相互作用及围岩压力时间效应, 初支围岩压力随时间呈下降趋势, 二衬围岩压力随时间呈上升趋势, 但二者之和基本保持不变。
(3) 蠕变稳定后, 主洞上台阶部分(拱顶至拱腰)的二衬荷载分担比约为25%~45%, 主洞下台阶部分(仰供至边墙)的二衬荷载分担比约为20%~35%。
[1] |
何川, 李玉文, 姚勇, 等. 公路小净距隧道[M]. 北京: 人民交通出版社, 2015. HE Chuan, LI Yu-wen, YAO Yong, et al. Road Neighborhood Tunnel[M]. Beijing: China Communications Press, 2015. |
[2] |
《中国公路学报》编辑部. 中国隧道工程学术研究综述·2015[J]. 中国公路学报, 2015, 28(5): 1-65. Editorial Department of China Journal of Highway and Transport. Review on China's Tunnel Engineering Research: 2015[J]. China Journal of Highway and Transport, 2015, 28(5): 1-65. |
[3] |
李鹏飞, 王帆, 聂雄, 等. 深埋非对称连拱隧道围岩压力计算方法研究[J]. 岩土工程学报, 2016, 38(9): 1625-1629. LI Peng-fei, WANG Fan, NIE Xiong, et al. Methods for Calculating Rock Pressure of Symmetrical Multi-arch Deep Tunnels[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(9): 1625-1629. |
[4] |
李然, 张顶立, 房倩, 等. 深埋3孔小净距隧道围岩压力计算方法及其工程应用[J]. 中国铁道科学, 2020, 41(2): 81-90. LI Ran, ZHANG Ding-li, FANG Qian, et al. Calculation Method for Surrounding Rock Pressure of Deeply Buried Closely Spaced Triple Tunnels and Its Engineering Application[J]. China Railway Science, 2020, 41(2): 81-90. |
[5] |
孙振宇, 张顶立, 房倩, 等. 浅埋小净距公路隧道围岩压力分布规律[J]. 中国公路学报, 2018, 31(9): 84-94. SUN Zhen-yu, ZHANG Ding-li, FANG Qian, et al. Distribution of Surrounding Rock Pressure of Shallow Highway Tunnels with Small Spacing[J]. China Journal of Highway and Transport, 2018, 31(9): 84-94. |
[6] |
李磊, 谭忠盛, 郭小龙, 等. 挤压陡倾千枚岩地层小净距隧道大变形研究[J]. 岩石力学与工程学报, 2019, 38(2): 276-286. LI Lei, TAN Zhong-sheng, GUO Xiao-long, et al. Research on Large Deformation of Tunnels with Small Intervals in Squeezing Steeply Dipping Phyllite Strata[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(2): 276-286. |
[7] |
谭忠盛, 孟德鑫, 石新栋, 等. 大跨小间距黄土隧道支护体系及施工方法研究[J]. 中国公路学报, 2015, 28(11): 82-89. TAN Zhong-sheng, MENG De-xin, SHI Xin-dong, et al. Research on Support System and Construction Method for Large-section Small Spacing Tunnel in Loess[J]. China Journal of Highway and Transport, 2015, 28(11): 82-89. |
[8] |
王飞飞, 江学良, 杨慧, 等. 浅埋偏压小净距隧道加速度响应规律试验与数值模拟研究[J]. 振动与冲击, 2017, 36(17): 238-247. WAGN Fei-fei, JIANG Xue-liang, YANG Hui, et al. Tests and Numerical Simulation for Acceleration Response Laws of a Shallow Buried Small Spacing Tunnel with Asymmetrical Pressure[J]. Journal of Vibration and Shock, 2017, 36(17): 238-247. |
[9] |
彭超, 童小东. 大断面公路隧道围岩压力时间效应分析[J]. 青岛理工大学学报, 2019, 40(5): 12-17. PENG Chao, TONG Xiao-dong. Analysis of Time Effect of Surrounding Rock Pressure in Large-section Highway Tunnel[J]. Journal of Qiaodao Technological University, 2019, 40(5): 12-17. |
[10] |
卢向勇, 陈伟庚, 邓皇适, 等. 考虑围岩蠕变特性的隧道仰拱开挖时序性研究[J]. 铁道科学与工程学报, 2020, 17(4): 900-907. LIU Xiang-yong, CHEN Wei-geng, DENG Huang-shi, et al. Study on Time Series of Tunnel Inverted Arch Excavation Considering Creep Characteristics of Surrounding Rock[J]. Journal of Railway Science and Engineering, 2020, 17(4): 900-907. |
[11] |
钱文喜, 耿大新, 梁国卿. 围岩蠕变对运营隧道衬砌安全性的影响[J]. 公路交通科技, 2020, 37(9): 90-96. QIAN Wen-xi, GENG Da-xin, LIANG Guo-qing. Influence of Surrounding Rock Creep on Safety of Lining of Operating Tunnel[J]. Journal of Highway and Transportation Research and Development, 2020, 37(9): 90-96. |
[12] |
耿大新, 钱文喜, 梁国卿. 围岩蠕变对运营隧道衬砌裂缝发育的影响研究[J]. 华东交通大学学报, 2018, 35(3): 8-14. GENG Da-xin, QIAN Wen-xi, LIANG Guo-qing. Influence of Wall Rock Creep on the Lining Crack Development in the Operating Tunnel[J]. Journal of East China Jiaotong University, 2018, 35(3): 8-14. |
[13] |
徐国文, 何川, 汪耀, 等. 流变荷载作用下隧道裂损二次衬砌结构安全性能研究[J]. 土木工程学报, 2016, 49(12): 114-123. XU Guo-wen, HE Chuan, WNAG 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. |
[14] |
任松, 王小书, 欧阳汛. 基于砂岩长期蠕变的隧道施工全过程动态模拟[J]. 中国矿业大学学报, 2019, 48(5): 1003-1011. REN Song, WANG Xiao-shu, OUYANG Xun. Dynamic Simulation Study of Construction Process of Simian Mountain Tunnel Based on Creep Characteristics of Sandstone[J]. Journal of China University of Mining & Technology, 2019, 48(5): 1003-1011. |
[15] |
XU G W, HE C, YANG Q H, et al. Progressive Failure Process of Secondary Lining of a Tunnel under Creep Effect of Surrounding Rock[J]. Tunnelling and Underground Space Technology, 2019, 90: 76-98. |
[16] |
ZHENG A C, HUANG F, TANG Z D, et al. Stability Analysis of Neighborhood Tunnels with Large Section Constructed in Steeply Jointed Rock Mass[J]. Mathematical Problems in Engineering, 2020(10): 1-14. |
[17] |
福建省交通规划设计院. 平潭综合实验区麒麟大道东段工程龙兴岭隧道设计文件[R]. 福州: 福建省交通规划设计院, 2017. Fujian Provincial Transportation Planning and Design Institute. Longxingling Tunnel Design Document for Eastern Section of Qilin Avenue in Pingtan Comprehensive Experimental Area[R]. Fuzhou: Fujian Provincial Transportation Planning and Design Institute, 2017. |
[18] |
李习平. 凝灰岩地层大断面隧道围岩流变特性及变形控制技术研究[D]. 长沙: 中南大学, 2013. LI Xi-ping. Research on Rheological Properties and Deformation Control Technology of Surrounding Rock in Tuff Formation Large-section Tunnel[D]. Changsha: Central South University, 2013. |