深埋隧道拱部围岩潜在塌落破坏上限有限元分析
隧道开挖后围岩受到扰动而应力重新分布[1]，形成局部落石[2]、剪切破坏[3]等危害隧道施工安全的破坏形式，且多以拱部岩体塌落形式发生[4−5]。隧道拱部塌落和围岩压力研究常用理论分析、数值模拟和现场监测等方法。如普罗托奇雅阔诺夫推出平衡拱理论；太沙基建立松散介质理论以确定衬砌垂直压力公式和滑动面[6]。目前，岩土极限分析成为解决隧道围岩稳定性与破坏模式问题的有力手段。如王永甫等[7]建立了圆形隧洞塌落破坏模式滑移线方程。武军等[8]以颗粒流椭球体理论改进太沙基松动压力公式。SAHOO 等[9]采用上限有限元分析多参数条件下圆形隧道围岩压力。LI等[10]运用Hoek-Brown 破坏准则及上限定理，明确了多圆弧截面深埋隧道坍塌机理。FRALDI 等[11]用变分法极限分析建立了矩形隧道塌落力学模型。HUANG 等[12]获得了孔隙水压力下隧道塌落数值解。刘俊等[13]采用强度折减法分析了浅埋隧道塌落破坏。白维仕等[14]以滑移线方式评价了黄土隧道坍塌拱及承载性能。孙雁军等[15]运用运动单元上限有限元揭示了隧道工作面滑移线网塌落形态。耿永旺[16]采用离散元分析了黄土隧道塌落拱形态特征。现有文献多倾向于研究浅埋隧道，如YANG等[17]对不同跨度影响下的浅埋隧道破坏进行了系列分析。通常隧道跨度增加时，围岩扰动范围更广，当隧道埋深较大时，拱部易形成潜在的塌落和松动范围，这种现象对于深埋土质隧道更为显著。本文应用运动单元上限有限元法，研究极限状态下不同跨度深埋隧道拱部的潜在塌落与松动范围，分析隧道稳定性系数和破坏规律，为施工技术方案和风险防范措施制定提供一定的参考和借鉴。
1 问题描述与假设条件
图1为本文建立的隧道拱部潜在塌落破坏分析模型图。假设模型满足如下条件：
1) 隧道围岩为均质各向同性材料，破坏服从摩尔-库伦屈服准则，容重为γ，内摩擦角为ϕ，剪胀角为ѱ，黏聚力为c；
2) 隧道处于未支护毛洞状态，且只考虑高度为h，宽度为B(半宽为b)的上台阶，隧道轮廓简化为椭圆形；
3) 隧道埋深足够大，荷载仅考虑围岩自重，模型依对称性取右侧一半。
陪伴英的是自己丈夫。今年82岁的丈夫走路歪歪扭扭，每走一步都给人重心失衡的感觉，一不小心就要翻跟头似的。老人咳嗽得十分厉害，每咳一次都像即将散架。英住在肿瘤中心四楼病房，丈夫不敢乘坐电梯，从来没有出过远门的老人害怕进了电梯会出不来，每次下楼老人都是扶着栏杆一步步往下挪。一日三餐，丈夫小心翼翼地给英喂饭，老人老眼昏花，双手发颤，几次都把米饭散落在地上。英吃上几口就不想吃了，实在是没有胃口，她不停地摇头。丈夫性子马上就来了，他呵斥了英几句，她又乖乖地一口一口吃，强迫自己将米饭咽下去。
如图1 所示，以b/h 表示无量纲跨度参数，其中h 不 变，b/h 依 次 取1.00，1.25，1.50，1.75，2.00。围岩内摩擦角ϕ 取35°～50°(按1°递增)。以ѱ/ϕ 表示剪胀角的影响，取0.0，0.2，0.4，0.6，0.8，1.0共6组值。
图1 隧道拱顶潜在塌落破坏分析模型
Fig.1 Analysis model diagram of potential collapse failure of tunnel vault
图1中l2为隧道埋深，设置l2/h大于20。l1表示模型水平方向范围，可取l1=l2。模型左侧边界Q1Q2水平速度为0(u=0)，下部边界Q3Q4 和右侧边界Q4Q5水平、竖向速度均为0(u=0，v=0)，隧道轮廓Q2Q3为自由边界。
考虑隧道拱部潜在塌落破坏由围岩自重引起，故选用无量纲量容重系数γh/c 作为评价指标。于是，极限状态下对应于临界容重γcr，容重系数γcrh/c 与隧道跨高比b/h，内摩擦角ϕ 和剪胀参数ѱ/ϕ 关联，即：
2 上限有限元模型建立
上限有限元模型划分为刚性单元，单元间设置速度间断线[18]。每个单元具有水平和竖向2 个自由度，极限状态下仅速度间断线上发生能量耗散。每个节点位置为变量，计算过程可实现间断线的优化调整，进一步通过多次网格更新形成由有效间断线网组成的近似滑移线网破坏模式。
4.维持平衡计分卡的长期效益。为了适应企业宏观环境、行业环境和自身条件的变化，平衡计分卡的评价指数也要与时俱进。因此，虎豹集团需要根据企业自身内部和外部环境的变化，对平衡计分卡的评估指数进行企业战略目标校验，确保平衡计分卡设计框架的长久运行。在企业信息化建设的基础上，随时报告并总结平衡计分卡的评价效果，收集和分析平衡计分卡的战略指标数据，建立平衡计分卡报告评价体系。同时，要建立平衡计分卡的定期维护团队，对考核结果及运行现状进行自我检查，从而形成长期有效的平衡计分卡创新绩效管理体系。
对于隧道拱部潜在塌落破坏分析，可将围岩临界容重γcr作为目标函数，由上限定理其非线性规划目标函数表达式为：
式中：Pd,i为第i条间断线上的耗散能，nd为速度间断线的总数。
非线性规划的约束条件如式3(a)～3(h)所示：
式(3(a))为速度间断线非线性约束，而ξ′i，ξ″i为间断线上的过程变量。式(3(b))～(3(f))分别为Q1Q5和Q2Q3上的几何边界条件，式中ng1～ng5对应于各个边界上的节点数目。式(3(g))约束单元面积为正值，其中ne为单元总数；式(3(h))表示求解临界容重γcr 时需施加的体力约束。Q1Q2，Q3Q4 和Q4Q5上的速度边界条件见上节。
应用能量等效屈服面将最初的土体强度参数ϕ和c 按剪胀角ѱ 的不同转化为等效参数ϕ*和c*[19]，如式(4)和式(5)所示：
在导师名上，有一部分论文授予名单上导师名为中文名，但在纸本上导师名为英文名；但反过来授予名单上导师名为英文名，纸本上导师名为中文名字的比较罕见。
3 临界容重系数γcrh/c计算结果
3.1 计算结果对比
当隧道断面为圆形时，临界容重γcr已有相关研究，对应于本文b/h=1的情况，以下均以γcrh/c进行对比分析。取ϕ=30°，ѱ=30°，对应的数据如表1。
表1 临界容重系数γcrh/c计算结果对比
Table 1 Comparison of calculation results of critical bulk
density coefficient γcrh/c
SAHOO等[9]3.15 ZHANG等[20]2.39 YANG等[17]2.35 BANERGEE等[21]1.95本文解2.50
可看出，本文γcrh/c 上限解处于SAHOO 和ZHANG Jian 计算结果之间，与SAHOO, ZHANG Jian ,YANG Feng 和BANERGEE 的计算结果差距分别为26.0%，4.6%，6.0%和22.0%。比较而言，本文数值偏大，这源于隧道按上台阶半断面的情况考虑，即半圆形状对应的潜在塌落破坏，故得出γcrh/c 计算数值稍大，论证运动单元上限有限元解答可靠。
3.2 临界容重系数γcrh/c规律分析
利用运动单元上限有限元对隧道塌落破坏时的临界地层容重进行分析，按上述不同b/h，ϕ 和ѱ/ϕ对应获得临界容重系数γcrh/c关系曲线见图2。
图2(a)为ѱ/ϕ=0 时，γcrh/c 上限解与内摩擦角ϕ以及b/h 的关系曲线。由图可知，γcrh/c 上限解随内摩擦角ϕ的变化较小，呈现出先增长后缓慢下降的趋势，而γcrh/c上限解随隧道跨度的增加而减小。
图2(f)为ψ/ϕ = 1时，γcrh/c上限解与内摩擦角ϕ以及b/h 的关系曲线。分析可知，隧道跨度增大将显著降低破坏时的临界容重γcr，如内摩擦角45°时，b/h=2 对应的临界容重系数γcrh/c 仅为b/h=1 时的39.5%，说明隧道跨度增加地层稳定性显著下降。此外，不同于图2(a)ψ/ϕ = 0 情况，ψ/ϕ = 1 时随内摩擦角ϕ 增长，γcrh/c 增长近似于线性增长关系。当隧道跨度大时，γcrh/c随内摩擦角ϕ增长的斜率变小。
通常，ψ/ϕ 介于0～1 之间，其值越接近于1 则剪胀效应越显著。对比图2(a)～2(f)可看出，γcrh/c上限解随剪胀参数ψ/ϕ 增大而增大；在这个过程中，γcrh/c 上限解曲线逐渐上扬，其数值与内摩擦角的关联度变大；反之，剪胀参数ψ/ϕ 较小时，γcrh/c与ϕ关系曲线转变为一略微上凸的抛物线。
图2 临界容重系数γcrh/c与影响因素关系
Fig.2 Relationship between critical bulk density coefficient γcrh/c and influencing factors
4 围岩潜在塌落破坏模式
4.1 围岩潜在塌落破坏模式对比
运动单元上限有限元程序可得到隧道拱部围岩潜在塌落破坏模式。选取ϕ=35°，ψ/ϕ = 0 和1 对应的破坏模式绘制如图3。
居民通过“环保大作战”回答环保小问题换取游戏券，在游戏中了解环保知识；“宝圈套环&垃圾分类”将垃圾模型丢入相应的干垃圾桶、湿垃圾桶、有害垃圾桶、可回收垃圾桶中，让居民有直观的印象。“环保易物”，居民将废报纸、废旧电池、二手书籍兑换长寿花。“文明养宠宣传”，居民签下文明养宠承诺书，将文明养宠文化在生活中实践，当天共签署承诺书50份。
图3展示了速度值分布云图，该值经过归一化处理，取值范围处于0～1 之间。图3 中还叠加有效速度间断线组成的滑移线网，并以虚线表示按太沙基围岩压力理论绘制的破坏面。
2.3 预防交叉感染 患儿因化疗致骨髓抑制，免疫力低下，另一方面，由于患儿年幼，家人强烈要求陪护，因此，需严格做好保护性隔离护理工作。谢绝来人探访，告知家属，有呼吸系统感染的家人避免照顾患儿;紫外线消毒房间，每日2次;定时通风，保持室内空气新鲜;遵医嘱及时准确使用升白细胞药物，每日3次观察体温变化。
图3 围岩潜在塌落破坏模式的相对速度矢量值云图(ϕ=35°)
Fig.3 Relative velocity vector nephogram of potential collapse failure of surrounding rock in tunnel arch(ϕ=35°)
由滑移线网看出，本文得到的破坏模式与太沙基理论的破坏面大致吻合。但本文破坏范围限制于一定高度范围内，这符合深埋隧道潜在塌落破坏限于局部的共识。由速度云图可知，可见耗散能占比前80%的块体单元集中于隧道拱部较小区域，该处可认为围岩潜在塌落范围，这与含平衡拱在内诸多理论相符[12,22−23]。在滑移线范围内还存在速度很小块体，此时间断线上耗散能很小，用上限原理仅能论证处于极限状态，与发生潜在塌落破坏还存在一些差异。
人民币汇率形成机制改革首先要基于我国国情和社会经济发展基本情况，同时改革方法要符合主动性、可控性以及渐进性原则。要充分考虑市场供求关系以及汇率水平。通过合理的控制汇率浮动空间增加人民币弹性。从而提升人民币汇率在维持市场内外部平衡中的作用。
很多人认为推迟退休年龄会增加就业压力。从劳动力市场的角度出发，推迟退休年龄意味着更多的岗位对于新增劳动力的需求量不大，这样就会导致劳动力市场供过于求的矛盾的出现，增加就业压力，失业率会明显提高，根据国家统计局2017年《中国统计年鉴》的数据显示，我国2016年城镇就业人口数41428万人，城镇登记失业人数近1000万人，城镇登记失业率达到4.02%。我国城镇失业人口总数还是很大的，就业形势十分严峻，就业压力逐年增长。
4.2 围岩潜在塌落破坏模式的参数影响分析
为揭示不同参数取值条件下隧道拱部围岩潜在塌落破坏几何形态的演变规律，绘制系列滑移线网破坏模式如图4所示。
1) 内摩擦角ϕ的影响
图4(a)～4(d)为圆形隧道ψ/ϕ = 0 和1，ϕ=35～50°时计算得到滑移线网破坏模式。可看出，当ϕ减小时，临界容重系数γcrh/c 随之减小，隧道上方破坏范围快速增加。例如，ψ/ϕ = 1 时，当内摩擦角ϕ 从40°降低到35°时，潜在塌落高度增长了59%，临界容重系数减少了19%。通常，内摩擦角降低围岩更易发生剪切破坏，引起更大范围塌落和松动趋势。
图4(e)展示了ϕ=35°～50°时，隧道拱部破坏范围外轮廓图，可看出塌落范围增长关联了内摩擦角减小，破坏模式顶端滑动面的倾角等于围岩内摩擦角。
2) 隧道跨度和高度比b/h的影响
图4(f)～4(i)列出了ψ/ϕ = 0 和1，ϕ=45°，b/h=1.25～2.0 对应滑移线网破坏模式。由图可知，当b/h 增大，即隧道跨度增大时，上方潜在塌落土体范围不仅在横向上增长，在竖向范围也有较大增长，这也是临界容重系数随隧道跨度增长而快速下降的原因。
图4(j)展示了b/h=1.0～2.0 对应的隧道拱部破坏范围外轮廓图。对比图4(e)可知，潜在塌落破坏范围受内摩擦角和隧道跨度两方面的影响，内摩擦角的改变引起塌落体形状，而不同跨度隧道对应的潜在塌落范围虽有改变，但外侧轮廓差异并不大。
“部门联动”“一呼百应”“执法统一”……吴建生直言这三个词是“平安西江”建设中令他感到最欣喜的。“虽然广东海事局有统一的执法标准和规定，但是以前辖区各分支局在实际执法时的标准、执法力度各不相同，而且各分支局之间合作机制不完善，一旦发生事故，分支局与分支局的沟通交流也存在一定的障碍，通过共建'平安西江’行动，这些问题统统都被解决了。”
3) 剪胀参数ѱ/ϕ的影响
图4(g)和图4(k)～4(n)为ϕ=45°，b/h=1.5，ψ/ϕ=0～1时的滑移线网破坏模式，可以看出随着剪胀参数增大，剪胀效应使得外侧围岩受到剪胀变形土体所产生的压力，使围岩塌落破坏范围减小，上限解增大。如ψ/ϕ 从0 增至1 时，塌落范围减小56%，上限解增长83%。
“可是，怎么整呢？写信吧，不行，现在的稿纸上都有号码，能查出来。打电话吧，又怕听出来，还是不行。后来，我想起几年前从单位拿回来的十来张有点宣纸性质的16开纸，这种纸现在市面上没有了，用它不易被查到，安全。我用左手写的字，最后连我自己都认不出来。
图4(o)展示了ϕ=45°，b/h=1.5 时不同ψ/ϕ 对应的外轮廓，可知塌落范围随ψ/ϕ 增大而减小趋势，随ψ/ϕ 增大而放缓，如ψ/ϕ 从0 到0.2，塌落面积下降32.1%，而ψ/ϕ 从0.8 到1，塌 落 面 积 仅 下 降2.9%。
建章立制，将资金的分配、管理、使用等全过程的纳入制度的笼子，是做好乡镇财政资金监管的最根本保障，也是最有效的途径。乡镇财政部门应按照上级财政部门的要求，特别是财政部《关于切实加强乡镇财政资金监管工作的指导意见》（财预〔2012〕28号）等文件，进一步做好制度的梳理和完善工作，有计划、有步骤地开展乡镇财政资金监管的制度建设，逐步构架好“乡镇资金监管的制度牢笼”。通过制度的完善，把乡镇财政资金监管的责任层层分解到岗位、到人员，努力实现乡镇财政资金的全方位、无缝隙监管。
图4 隧道拱部围岩潜在破坏模式演变规律
Fig.4 Evolution law of potential failure mode of surrounding rock at tunnel arch
5 结论
1)当ψ/ϕ = 1 时，临界容重系数γcrh/c 随ϕ 的增长而线性增长，随b/h的增长而快速降低；当ψ/ϕ =0 时，γcrh/c 随ϕ 增长的变化不显著，而与b/h 呈反比例关系。
2)内摩擦角ϕ增长时，隧道拱顶围岩潜在塌落范围增长迅速，且塌落形状有所变化；b/h 增长时塌落范围增长放缓，且塌落拱外轮廓保持不变。
3)当ϕ增长时，围岩破坏范围离拱顶稍远区域相对速度变得很小，b/h 增长使得拱顶上方速度集中区域范围增长，同时γcrh/c迅速减小。
参考文献：
[1] 张顶立,台启民,房倩.复杂隧道围岩安全性及其评价方法[J].岩石力学与工程学报,2017,36(2):270−296.
ZHANG Dingli, TAI Qimin, FANG Qian. Safety of complex surrounding rock of tunnels and related evaluation method[J]. Chinese Journal of Rock Mechanics and Engineering,2017,36(2):270−296.
[2] 曲海锋,杨重存,朱合华,等.公路隧道围岩压力研究与发展[J].地下空间与工程学报,2007(3):536−543.
QU Haifeng, YANG Zhongcun, ZHU Hehua, et al.Research and development on surrounding rock pressure of road tunnel[J]. Chinese Journal of Underground Space and Engineering,2007(3):536−543.
[3] 刘高,聂德新,韩文峰.高应力软岩巷道围岩变形破坏研究[J].岩石力学与工程学报,2000(6):726−730.
LIU Gao, NIE Dexin, HAN Wenfeng. Deformation and failure of surrounding rocks of roadway in high stressed soft rocks[J]. Chinese Journal of Rock Mechanics and Engineering,2000(6):726−730.
[4] 《中国公路学报》编辑部.中国隧道工程学术研究综述·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.
[5] 孙振宇,张顶立,房倩,等.基于超前加固的深埋隧道围岩力学特性研究[J].工程力学,2018,35(2):92−104.
SUN Zhenyu, ZHANG Dingli, FANG Qian, et al.Research on the mechanical property of the surrounding rock of deep-buried tunnel based on the advanced reinforcement[J].Engineering Mechanics,2018,35(2):92−104.
[6] TERZAGHI K.Rock defects and loads on tunnel supports[R]. Proctor RV,White TL, editors. Rock tunnelling with steel supports, vol. 1. Youngstown, OH: Commercial Shearing and Stamping Co.,1946:17−99.
[7] 王永甫,王成,唐晓松.基于非关联流动法则的圆形隧洞滑移线解答[J].应用数学和力学,2013,34(12):1285−1290.
WANG Yongfu, WANG Cheng, TANG Xiaosong. Slip line solution for circular tunnels based on the nonassociated flow rule[J]. Applied Mathematics and Mechanics,2013,34(12):1285−1290.
[8] 武军,廖少明,张迪.基于颗粒流椭球体理论的隧道极限松动区与松动土压力[J].岩土工程学报,2013,35(4):714−721.
WU Jun, LIAO Shaoming, ZHANG Di. Loosening zone and earth pressure around tunnels in sandy soils based on ellipsoid theory of particle flows[J]. Chinese Journal of Geotechnical Engineering,2013,35(4):714−721.
[9] SAHOO J P, KUMAR J. Stability of long unsupported twin circular tunnels in soils[J]. Tunnelling and Underground Space Technology,2013,38:326−335.
[10] LI Yongxin, YAO Cong. Collapse mechanism of deep tunnels with three-centered arch cross section[J]. Journal of Central South University,2016,23(12):3293−3301.
[11] FRALDI M, GUARRACINO F. Limit analysis of collapse mechanisms in cavities and tunnels according to the Hoek-Brown failure criterion[J]. International Journal of Rock Mechanics and Mining Sciences, 2009, 46(4):665−673.
[12] HUANG Fu, YANG Xiaoli. Upper bound limit analysis of collapse shape for circular tunnel subjected to pore pressure based on the Hoek–Brown failure criterion[J].Tunnelling and Underground Space Technology, 2011,26(5):614−618.
[13] 刘俊,刘新荣,赖勇,等.不同覆跨比下浅埋软弱隧道的破坏模式[J]. 中南大学学报(自然科学版), 2016, 47(5):1744−1751.
LIU Jun, LIU Xinrong, LAI Yong, et al. Failure mode of shallow-buried weak tunnel under different thicknessspan ratios[J]. Journal of Central South University(Science and Technology),2016,47(5):1744−1751.
[14] 白维仕,李荣建,赵学勐,等.基于滑移线网络法的黄土隧道坍塌拱分析及其承载评价[J].防灾减灾工程学报,2020,40(1):132−138.
BAI Weishi, LI Rongjian, ZHAO Xuemeng, et al. Based on loess arch sliding line field visualization program design and parameter sensitivity analysis[J]. Journal of Disaster Prevention and Mitigation Engineering, 2020,40(1):132−138.
[15] 孙雁军,阳军生,罗静静,等.隧道工作面稳定性与滑移线网破坏模式研究[J].岩土工程学报,2019,41(7):1374−1380.
SUN Yanjun, YANG Junsheng, LUO Jingjing, et al.Stability and mesh-like collapse mechanism of tunnel face[J]. Chinese Journal of Geotechnical Engineering, 2019,41(7):1374−1380.
[16] 耿永旺.黄土隧道围岩破坏模式离散元分析[J].科协论坛(下半月),2013(2):9−10.
GENG Yongwang. Discrete element analysis of surrounding rock failure mode of loess tunnel[J]. Forum of the Association for Science and Technology (Second Half of the Month),2013(2):9−10.
[17] YANG Feng, SUN Xinlei, ZHENG Xiangcou, et al.Stability analysis of a deep buried elliptical tunnel in cohesive-frictional (c-ϕ) soils with a nonassociated flow rule[ ]. Canadian Geotechnical Journal, 2017, 54(5): 736−741.
[18] 杨峰,赵炼恒,张箭,等.基于刚体平动运动单元的上限有限元研究[J]. 岩土力学, 2014, 35(6): 1782 −1786,1808.
YANG Feng, ZHAO Lianheng, ZHANG Jian, et al.Investigation on finite element upper bound solution based on rigid translatory moving element[J]. Rock and Soil Mechanics,2014,35(6):1782−1786.
[19] DRESCHER A, DETOURNAY E. Limit load in translational failure mechanisms for associative and nonassociative materials[J]. Géotechnique, 1993, 43(3): 443−456.
[20] ZHANG Jian, YANG Feng, YANG Junsheng, et al.Upper-bound stability analysis of dual unlined elliptical tunnels in cohesive-frictional soils[J]. Computers and Geotechnics,2016,80:283−289.
[21] BANERJEE S K, CHAKRABORTY D. Seismic stability of a long unlined circular tunnel in sloping ground[J].Canadian Geotechnical Journal,2016,53(8):1346−1352.
[22] 殷洪波,刘正初,郭永发,等.饱和粉细砂地层隧道围岩塌落机理及施工控制研究[J]. 建筑机械, 2018(6): 47−52,5.
YIN Hongbo, LIU Zhengchu, GUO Yongfa, et al. Study on the collapse mechanism and construction control of saturated fine sand tunnel[J]. Construction Machinery,2018(6):47−52,5.
[23] 曾中林,左学贤,王志斌,等.二维圆形截面浅埋硐室双层围岩塌落形态上限变分分析[J]. 自然灾害学报,2019,28(5):96−104.
ZENG Zhonglin, ZUO Xuexian, WANG Zhibin, et al.Upper bound variation analysis of two-layer rock collapse overlying a two dimensional circular shallow buried chamber[J]. Journal of Natural Disasters, 2019,28(5):96−104.