基于收敛－约束原理的岩质隧道初衬安全系数计算

《基于收敛－约束原理的岩质隧道初衬安全系数计算》由会员分享，可在线阅读，更多相关《基于收敛－约束原理的岩质隧道初衬安全系数计算（8页珍藏版）》请在装配图网上搜索。

1、基于收敛约束原理的岩质隧道初衬平安系数计算摘要：以越岭公路岩质隧道为研究对象，构建隧道围岩与支护结构的平安系数。首先，基于围岩开挖后二次应力场的空间效应，确定隧道模型的边界条件、尺寸要求及应力条件。结合实摘要：以越岭公路岩质隧道为研究对象，构建隧道围岩与支护结构的平安系数。首先，基于围岩开挖后二次应力场的空间效应，确定隧道模型的边界条件、尺寸要求及应力条件。结合实际的隧道工程施工工艺，建立具备开挖支护功能的隧道三维数值模拟研究模型，以该模型作为研究隧道开挖与支护结构耦合分析的根底。基于数值模拟计算，分别构建围岩收敛曲线及纵向变形规律曲线。其次，建立隧道围岩与支护系统的收敛-约束图，并基于收敛-

2、约束原理建立岩质隧道初衬平安系数求解方法。最后，通过案例分析，详细展示上述过程。结果说明该平安系数法可以量化表征围岩开挖与支护过程中的隧道平安性问题。同时，假设以平安系数为量化分析指标，那么能够对隧道施工工艺进行优化。本文源自常州大学学报(自然科学版) 2022-01-28?常州大学学报?杂志办刊宗旨：立足科学开展前沿，秉承科学精神，以贮存和传播科技信息、繁荣科学文化、促进学术交流、发现和培养人才、推动科技进步为使命。关键词：隧道工程;围岩;平安系数;收敛-约束原理;收敛曲线隧道工程所处环境十分复杂，使得隧道围岩的开挖工艺、支护方法及开挖与支护间的动态耦合分析较为繁琐，相关过程具有一定的不确定

3、性1-2。加之中国相关标准对支护结构的设计建议仍以半经验法为主3-4，没有成熟的平安系数概念，使得在隧道工程中难以量化评估隧道围岩的稳定性或支护结构的平安性，给广阔设计人员带来了一定的困惑。现有尝试中，一些研究通过理论解析法或模型试验法来解决上述问题。如在考虑支护反力作用的情况下，结合Schwarz交替法获得隧道原位扩挖围岩力学特性的解析算法5，该研究丰富了围岩力学响应机制。再者，在围岩大变形控制中，可以通过现场试验或模型试验6-7，研究隧道开挖中的围岩渐进性破坏过程及位移和应力变化规律，从而优化施工方案，对保障围岩平安具有较强的实际意义。此外，在数值模拟中，通过分析隧道围岩破坏模式的差异性，

4、获得了与试验方法吻合较好的数值模拟结果8，其建模方法及模拟过程具有较好的借鉴性。同样，在围岩变形的平安分析中，有研究针对软弱质岩石环境下的围岩变形机制建立极限状态方程9，可以为岩质隧道工程设计及平安分析提供参考依据。此外，掌子面稳定岩体的最小平安厚度计算方法10的提出，可以为隧道穿越断层破碎带的稳定性及平安防护提供清晰的求解方案，进一步解答围岩稳定相关机理，具有显著的实际工程价值。同时，结合强度折减法11，评估围岩参数与隧道工程稳定性间的演化关系，对控制施工进度、优化支护方案、提高开挖效率有一定的指导意义。然而，在理论研究中，隧道模型大多以圆形开挖断面为主，而现实中隧道为非圆形断面居多，与实际

5、状况有一定差异性。同时，隧道围岩开挖与支护系统间的动态耦合机理尚未形成统一意见，为促进施工工艺的优化，相关研究亟待加强。整体上，隧道工程的复杂性使其一直未有诸如房屋建筑结构所用的理论清晰、方法简易的平安系数求解方法，以致在隧道的设计、施工与后期评价等环节中，难以量化表述、分析或评估支护结构的平安性或围岩的稳定性，给隧道工程的设计与施工带来一定的困扰。本文以表征隧道围岩与支护结构间动态耦合关系的收敛-约束方法为根底，通过构建隧道平安系数，对支护结构的平安性及围岩的稳定性进行量化分析。该方法在评估支护设计与优化施工工艺方面具有一定的实际意义。主要工作有：结合围岩开挖扰动时力学响应机制，建立岩质隧道

6、三维数值模拟研究模型;根据数值模拟分析，确定隧道围岩与支护结构间收敛-约束关系;根据收敛曲线与支护特征曲线的力学响应机理，建立隧道围岩稳定及支护结构平安性的量化分析方法;最后，通过一个开挖断面为曲墙拱形的岩质隧道工程案例，对上述方法进行完整展示与说明。1数值模型及支护结构特征函数1.1数值模型隧道围岩开挖与支护结构的安装显然是一个三维动态过程。一个大尺寸的隧道数值模型可以获得较为精确的结果。基于隧道洞室开挖后岩体二次应力场的空间效应(如图1所示)，发现沿隧道径向大约46倍半径范围外，由于应力重分布作用而造成的影响逐渐趋向稳定，切向应力根本恢复到初始状态。因此，为提高数值模拟的计算效率，文中拟建

7、的大尺寸隧道模型的开挖断面边界取5倍半径。隧道图1岩体二次应力的影响范围横断面的网格划分采用放射式方法，即靠近隧道断面中心，网格越密集，远离隧道区域那么较为疏松。沿隧道纵向的网格划分尽量保持较好的长宽比，使各个单元尽量均衡化。模型中岩体采用实体单元，衬砌采用壳单元。在隧道横断面(x-z向)，两边及底部设置固定边界条件，顶部为自由边界条件，在各边界施加相应的应力场v，hx，如图2中所示。在纵断面(y-z向)，各范围的边界条件与横断面保持一致。在模型建立时，依据实际的隧道施工工艺，通过设置开挖步进行围岩开挖的模拟过程。在模拟开挖过程中，根据设置的监测点，自动记录围岩的收敛位移。按照文献12，当隧道

8、掘进是在一个非均匀的初始应力场时，初始应力场可由地层的两向的水平应力h1，h2确定，可参照文献13-14所建议的求解方法。1.2支护结构特征函数RABCEWICZ15、朱浮声16等提出了巷道按组合拱原理进行加固设计的方法，分别可以获得喷射混凝土、锚杆、金属网及岩石承压拱的支护抗力。文献17-18建议了复合支护结构的极限抗力。基于上述成果，结合各支护材料的刚度、支护单元的抗力，按照组合支护理论建立支护结构特征函数：式中：psho，pbol，proc，pnet为喷射混凝土、锚杆、岩石承压拱、金属网的支护抗力;ksho，ulim，Esho，sho，tsho为喷射混凝土的刚度、变形量、弹性模量、泊松比

9、及厚度;r为隧道半径;p为复合支护系统的极限抗力。式(2)那么为基于支护结构单元的材料属性建立的支护结构特征曲线。在确定组合刚度k及极限变形值ulim后，那么式(2)可以获得求解。2隧道结构平安系数建立方法2.1收敛-约束原理通常，将反映围岩变形特征的曲线称为收敛曲线。在支护结构安装后，作用在支护结构上的荷载与支护变形的关系曲线称为支护结构特征曲线。收敛-约束原理是以圆形洞室为根底的弹塑性分析方法，主要解释围岩收敛变形与支护抗力间的动态耦合关系，从而为工程设计与实践提供理论依据。根据围图3收敛-约束原理示意图Fig.3Schematicdiagramofconvergence-confinem

10、entprinciple岩收敛曲线及支护结构特征曲线，那么可以建立隧道围岩开挖与支护结构间的收敛-约束图，如图3所示。文中，收敛曲线由数值模拟计算而得。在收敛-约束关系中，关于支护特征曲线起始位置确实定，如图3中的B点，那么是支护时机确实定问题，对施工工艺的优化有重要意义。确定方法为依据施工工艺，在岩体开挖后到支护结构安装时，所对应的围岩收敛值。图3中，plim为支护结构的极限抗力。u0为支护结构安装时围岩已产生的径向变形，由围岩爆破开挖后所产生的二次应力产生。同时，收敛曲线与支护特征曲线的交汇点A(uequ，pequ)，即为围岩变形与支护抗力之间动态耦合的平衡点。2.2支护结构平安系数计算图

11、3中，plim为支护结构极限承载力，点A所对应的纵坐标值pequ为围岩处在平衡状态时所施加在支护结构上的变形压力。那么根据平安系数的定义(通常为结构抗力与荷载之比)，平安储藏p及支护结构平安系数fs分别为：式(6)为根据收敛-约束原理建立的隧道支护结构平安系数，假设对该式进行求解，那么还需获取隧道的纵向变形规律曲线，用该曲线或函数确定图3中B点的位置，即支护时机确实定问题。2.3纵向变形规律曲线建立方法围岩开挖断面的纵向变形规律函数(LongitudinalDeformationProfile，文中简称LDP曲线)能直观地反映围岩径向变形随掌子面推进而动态变化的关系，对施工有很好的指导作用，如

12、图4所示。LDP曲线可由模拟的围岩开挖掘进过程获得。此外，建议设置较多数量的模拟开挖循环过程，才能获得较为完整的LDP曲线，所以模型轴向长度建议取2030倍洞径。每个循环中，围岩掘进进尺依据实际取13m。每个开挖循环后，布置的监测点自动记录断面的径向位移。将离掌子的间距设置为横坐标，径向位移为纵坐标，那么依据监测点数据可以建立隧道开挖的LDP曲线。2.4支护时机与施工优化分析LDP曲线能够指导支护结构的安装时机。将其参加到收敛-约束图中，那么能够确定不同施工工艺下，支护结构特征曲线(如图5中Sr，1，Sr，2，Sr，3)在收敛-约束图中的起始位置u1，u2，u3。在图5中，假定l1，l2，l3

13、分别代表3种施工工艺，即分别为1，2，3个施工循环后再进行支护结构的安装。那么各掘进进尺在LDP曲线上的对应点分别为Qr，1，Qr，2，Qr，3，这3点的横坐标即支护结构安装时隧道径向位移分别为图5支护结构安装时机示意图Fig.5Schematicdiagramofinstallationtimeofsupportingstructure(u1，0)，(u2，0)，(u3，0)。支护特征曲线与收敛曲线的相交点所对应的纵坐标值即为求解平安系数时的平衡点支护压力pequ，1，pequ，2，pequ，3。对于施工优化分析，根据极限状态下的支护抗力pequ，1，pequ，2，pequ，3，结合图5，不

14、难发现当l1，l2，l3不同，支护效果也不尽一致。总体上，3个位置所对应的平衡点值由大到小为pequ，1，pequ，2，pequ，3。该关系说明经施工工艺的优化，即经屡次开挖循环后再进行支护结构的安装，所需的最小支护抗力逐渐减少，那么根据计算式(5)、式(6)，相应的平安充裕量及平安系数显著提升。具体过程在本文案例分析中进行详细介绍。2.5隧道围岩变形控制方案为利用岩体自稳能力，允许爆破后岩体出现一定的变形，释放局部变形压力。但是，为保障隧道施工平安，支护结构安装前的变形量应有严格的约束。量化表征时，图5中u1，u2，u3作为3个不同的支护系统安装起始位置，代表了支护结构安装时围岩已发生的径向

15、位移量，对其值应严格控制。即隧道平安定量描述方法除满足式(6)要求外，为约束围岩变形量、保障岩体稳定性，还应满足式(7)要求。式中：equ，l为围岩的径向应变;uequ，l为采用不同施工工艺时围岩径向位移量;R为隧道半径。式(7)中，根据研究成果，径向应变equ，l的值不宜大于1%19-20。所以，基于不同的施工工艺，那么可以获得不同的平安系数。结合式(7)，以平安系数为量化分析指标，那么可以对施工工艺进行优化。3工程案例分析3.1工程背景某岩质隧道的开挖高度、跨度等信息如图6所示。该隧道最大覆盖层厚度约350m。岩体的弹性模量为1.7GPa、泊松比为0.33、黏聚力为0.45MPa、内摩擦角

16、分比为30.85。该隧道为曲墙拱型开挖断面，钻爆法施工，同时岩体质量适中，具有较强的代表性。所以，选取该工况作为文中方法例如对象。依据文献17，可以求得该工程实例支护系统极限变形ulim=3.083mm，极限抗力plim=0.693MPa。3.2计算结果分析根据工程资料，建立该案例的数值模型，如图7(a)所示。经逐级减小支护抗力，结合监测点记录数据，那么可以绘制该工程实例的收敛曲线，如图7(b)所示。LDP曲线需要考虑隧道围岩开挖与支护过程的动态变化关系，因此，采用三维模型建立方法十分必要。根据前述二次应力影响范围，该工程实例模型及边界尺寸分别如图8所示。根据不同施工工艺所建立的LDP曲线如图

17、9所示。图9那么可以展示不同施工工艺下隧道断面的最大径向位移与掌子面间的动态变化关系。基于该关系，那么可以确定收敛-约束图中支护结构特征曲线的起始位置，即支护时机确定问题。例如，依托该工程案例，为比照分析，本文共计列举了3种施工工艺，分别为开挖循环13m。根据图9，选取7组数据，建立不同施工工艺下的隧道围岩与支护结构收敛-约束图，如图10所示。在图10中，可清晰的看到各施工工艺下的支护结构特征曲线，如(2，3;6m)，那么表示开挖循环为2m，在3个开挖循环后进行支护结构的安装(即未支护的岩体距离掌子面间距为6m)。所以，根据平安系数求解方法，结合收敛-约束图，那么可以求解各个施工工艺下的隧道结

18、构平安系数，见表1。3.3基于平安系数的施工工艺优化方法该案例中，有几组施工工艺的平安系数小于1，理论上该隧道支护结构处在危险状态，如(1，3)，(2，3)，(3，2)组。解决方式可以采用加大支护抗力或采用施工工艺优化方法。对于施工工艺的优化，假设在开挖进尺不变的情况下，待多个爆破开挖循环后再进行支护结构安装，即允许释放局部围岩压力。此时，正如表1计算的结果，可以发现在支护参数不变的情况下，隧道支护结构的平安系数获得逐步提高，如fs(1，6)fs(2，3)fs(3，2)。当然，在施工工艺的优化过程中，循环次数不能过多，即围岩变形不能过大。依据文中式(7)所列控制条件，即在采用的施工工艺下，未支

19、护的岩体径向应变不能大于1%。假设经支护优化设计后平安系数仍小于1，同时径向应变已超1%，那么建议重新设计支护单元及支护参数。所以，在平安系数及围岩允许变形量双重指标的控制下，可以有效利用围岩自稳能力，提高支护效力，促进施工工艺的进一步优化。4结论1)根据地下工程开挖扰动时围岩应力响应原理，确定数值模型的建模方法，梳理收敛-约束图构建的主要过程，并建议围岩收敛曲线、支护结构特征曲线、LDP曲线的数值模拟建立方法。2)基于数值模拟研究，确定隧道围岩与支护结构收敛-约束关系。根据收敛-约束图中支护结构与隧道收敛变形的平衡状态，并结合平安系数定义，提出隧道工程支护结构平安系数求解方法。3)通过工程实例分析，完整展示了上述方法的计算流程。研究说明，平安系数方法能够量化分析隧道支护结构的平安性及围岩的稳定性。隧道开挖掘进时，围岩爆破循环进尺、循环次数、支护安装时机等施工工艺对平安系数值影响较大。假设以平安系数为量化分析指标，那么能够对施工工艺进行合理优化。此外，根据围岩变形量、施工工艺的优化应结合围岩自身稳定能力而综合考虑。