技术摘要：
本发明公开了一种隧道开挖上覆岩层开裂预报方法及围岩支护方法，属于隧道施工技术领域。本发明包括以下步骤：1)根据地质勘探结果确定正交各向异性弹性参数；2)确定隧道施工参数；3)建立有限元模型；4)初始地应力计算；5)根据开挖和支护步序开展扩展有限元仿真计算，获  全部
背景技术：
隧道是构筑在水下、地下或山体中的供车辆行驶的专门建筑，相比一般建筑，其施 工难度大、安全性要求高。现行常用的隧道开挖方法为新奥法，该法以控制爆破为开挖手 段，以锚杆、喷射混凝土为主要支护方法，通过现场监测以及数值计算获得围岩的变形量， 从而动态修正设计参数和变动施工方法。 在隧道工程施工过程中所进行的围岩变形监控量测，是新奥法隧道施工技术的重 要组成部分。隧道围岩变形监测成果，不仅可用于反馈修正支护参数，而且对于评价围岩稳 定性、保障现场人员与设备的安全具有重要意义。现有的用于计算围岩变形量的方法主要 是常规有限元方法，然而，控制爆破过程中易造成周围岩体的破裂，形成较大裂隙，而常规 有限元方法的模拟计算基于连续介质假设，没有考虑裂缝强间断面对围岩变形量的影响， 无法反映上覆岩层离断破坏对围岩变形量的影响，不能反映新奥法施工的真实状态，因此 对于设计参数的变更和施工方法的调整易造成不利影响，从而威胁工程安全。 经检索，中国专利申请号为2019102597490的申请案公开了一种隧道极限位移的 确定方法，该申请案的隧道极限位移确定方法包括建立有限元数值模型；确定围岩主要力 学参数；隧址区地应力采用自重地应力场；确定关键监测点；模拟开挖与支护过程；校核模 型合理性；将围岩强度除以折减系数，模拟开挖与支护过程；若计算未达到停止条件，则返 回上一步并将折减系数增大，直至达到停止条件；提取达到停止条件的前一个折减步的关 键点位移；将所得的关键点位移乘以可量测系数作为隧道极限位移，供监控量测预警使用。 该申请案通过围岩强度折减获得隧道围岩的极限状态和极限位移，能考虑隧道断面型式、 埋深、施工工序与支护结构等因素对隧道极限位移的影响，但其没有考虑岩体材料的各向 异性特征，无法反映上覆岩层离断破坏对围岩变形量的影响。
技术实现要素：
1.要解决的问题 本发明的目的在于克服现有隧道开挖围岩变形模拟预测方法的准确性相对较低， 不能反映新奥法施工真实状态的不足，提供了一种基于扩展有限元的隧道开挖上覆岩层开 裂预报方法及围岩支护方法。本发明综合考虑了开挖扰动、岩层破裂等因素的影响，从而能 够得到具有工程指导意义的围岩变形量和裂缝扩展路径的预报结果，因而能够为新奥法施 工过程中动态修正支护施工设计参数提供合理的数据支撑。 2.技术方案 为了解决上述问题，本发明所采用的技术方案如下： 本发明的一种基于扩展有限元的隧道开挖上覆岩层开裂预报方法，视岩体为正交 5 CN 111553101 A 说　明　书 2/5 页 各向异性材料，利用扩展有限元方法计算出隧道开挖过程中的应力场和位移场，进而计算 得到上覆岩层的裂缝扩展路径，具体包括以下步骤： 步骤1、根据地质勘探结果确定岩体力学参数； 步骤2、确定隧道施工参数； 步骤3、确定隧道开挖工程的模型尺寸，并创建有限元模型； 步骤4、计算初始地应力； 步骤5、根据开挖步序开展扩展有限元仿真计算，得到应力场和位移场，根据加权 平均最大主拉应力准则确定裂缝扩展方向，从而获得上覆岩层开裂路径和围岩变形量。 更进一步的，所述步骤5具体包括以下步骤： 步骤5.1、利用扩展有限元法描述裂缝，模型内任意点x的扩展有限元位移场u(x) 由下式表示： 其中Sall、Sfrac和Stip分别是常规节点集合、跳跃增强节点集合以及裂尖增强节点 集合；uI、aI和 分别是常规节点自由度、跳跃增强节点自由度和裂尖增强节点自由 度； 为x点的单元形函数；H(x)为跳跃增强函数；Fl(x)为正交各向异性材料裂尖增强 函数，其表达式为： 其中，r为θ裂尖局部极坐标系坐标，pk和tk表达式分别如下： 上式中，skx、sky为sk的实部和虚部，sk经求解下式获得： S 4 3 211s -2S16s (2S12 S66)s -2S26s S22＝0 其中，Sij为材料柔度矩阵系数； 步骤5.2、根据裂缝在模型中所处位置，确定增强节点和增强单元； 步骤5.3、获得各单元的刚度矩阵，并组集总体刚度矩阵； 步骤5.4、通过对单元刚度矩阵乘以系数10-6实现开挖过程的模拟，通过对单元刚 度矩阵乘以系数106实现衬砌支护过程的模拟，通过修改材料参数实现岩体加固过程的模 拟； 步骤5.5、解线性方程组，获得节点位移，并计算高斯点应力； 步骤5.6、根据加权平均最大主拉应力准则确定裂缝扩展方向。 更进一步的，所述步骤5.2具体是指：若某单元被裂缝完全贯穿，则将该单元及其 附属节点标记为跳跃增强单元和跳跃增强节点；若某单元被裂缝部分穿过，则将该单元及 其附属节点标记为裂尖增强单元和裂尖增强节点。 6 CN 111553101 A 说　明　书 3/5 页 更进一步的，所述步骤5.3中单元刚度矩阵计算式为： 其中，i、j表示单元节点编号，且 上式中的B为形函数导数矩阵，D为各向异性材料应力应变关系矩阵，Ωe表示单元 积分区域。 更进一步的，所述步骤5.6具体步骤为： 步骤5.6.1、裂尖应力状态由以裂尖为中心、半径为r的圆形区域内其他高斯点的 应力状态加权平均求得，采用下式计算裂尖加权应力： 其中，σi为第i个高斯点的应力，wi为第i个高斯点的权重； 步骤5.6.2、根据加权应力计算主应力σ1，若σ1大于岩层的抗拉强度，则裂缝发生扩 展；其中扩展步长Δl和扩展角θ分别按照下式进行计算： 其中， 表示全部增强单元的平均面积；γ为裂缝倾角， 为最大主应力方向与x 轴的夹角。 更进一步的，所述权重wi按照下式进行计算： 其中li为第i个高斯点到裂尖的距离。 更进一步的，所述跳跃增强函数H(x)在裂缝面两侧分别取值1和-1。 更进一步的，所述岩体力学参数包括各向异性弹性模量、剪切弹性模量、泊松比、 密度和抗拉强度，所述隧道施工参数包括隧道断面参数、开挖步序、支护结构。 本发明的一种隧道开挖围岩支护方法，采用本发明的方法对隧道开挖上覆岩层开 裂路径和围岩变形量进行预报计算，并根据计算所得岩层开裂路径和围岩变形量调整开挖 步序或制定支护方案，重复步骤4至5，直至上覆岩层稳定。 3.有益效果 相比于现有技术，本发明的有益效果为： (1)本发明的一种基于扩展有限元的隧道开挖上覆岩层开裂预报方法，通过将岩 体看作是正交各向异性材料，综合考虑了开挖扰动、岩层破裂等的影响，采用扩展有限元方 7 CN 111553101 A 说　明　书 4/5 页 法对隧道开挖围岩变形量和裂缝扩展路径进行模拟预报，其模拟结果能够反映新奥法施工 的真实状态，因此可用于指导实际隧道施工，为新奥法施工过程中动态修正支护施工设计 参数提供了数据支撑。 (2)本发明的一种基于扩展有限元的隧道开挖上覆岩层开裂预报方法，通过对有 限元模拟计算的具体计算公式进行优化，尤其是对裂尖增强函数以及裂尖加权应力计算公 式进行优化，从而保证能够精确捕捉裂缝的扩展路径，有利于保证覆岩层开裂预报结果的 准确性。 (3)本发明的一种隧道开挖围岩支护方法，通过对围岩变形量和裂缝扩展路径模 拟计算方法进行优化，有效提高其预报的准确性，然后依据预报结果不断对开挖步序或支 护方案进行优化调整，因此能够有效保证隧道开挖围岩支护的准确性和安全性。 附图说明 图1是本发明的基于扩展有限元的隧道开挖上覆岩层开裂预报方法及围岩支护方 法的流程图。 图2是本发明实施例所提供的有限元模型示意图。 图3是本发明实施例所提供的裂缝增强单元和增强节点示意图。 图4是本发明实施例所提供的未考虑上覆岩层开裂情况下的围岩竖直方向位移云 图。 图5是本发明实施例所提供的考虑了上覆岩层开裂情况下的围岩竖直方向位移云 图。 图6是本发明实施例所提供的围岩等效位移曲线。