基于离散元法的巷道围岩非线性大变形分析与控制.docxVIP

基于离散元法的巷道围岩非线性大变形分析与控制 1 深部巷道围岩支护机制 深部岩石位于复杂地质环境中。高肥力、高地形、高渗透和强时间效应会导致深体和平面岩石的组织结构、基本行为特征和工程反应。具体体现为在浅部表现为坚硬稳定的围岩,在深部一般都表现出显著的软岩特征。对深部巷道而言,由于存在着高赋存地应力和围岩体低强度之间的突出矛盾,巷道开挖后造成的应力集中和应力重分布将导致很大范围内的围岩进入峰后破裂碎胀状态。因此,深部巷道围岩多为破坏后的裂隙岩体或破裂块体,采用的各种支护结构与支护方式也主要用来限制破裂岩体的大变形。然而,由于支护机制没有搞清楚,长期以来深部地下工程支护设计还主要依靠工程经验来进行,已有的理论分析方法也大多局限于圆形洞室在静水压力条件下的解析解,且很难全面考虑支护结构对围岩体的加固作用。近年来,一些学者借助于有限单元法、有限差分法等连续介质数值模拟方法,对其支护机制及相应支护对策进行了大量研究,并已得到一些有价值的成果。但由于深部破裂围岩多为被结构面切割而成的镶嵌块体,具有非线性、非连续变形特征,因此采用以小变形假设为前提的连续介质力学方法来研究,既与实际工程围岩普遍处于破裂状态的情况严重不符,也无法模拟破裂围岩的非线性大变形力学特性。而基于离散单元法发展而来的UDEC(universal distinct element code)软件目前已广泛应用于采矿工程、隧洞工程、边坡工程等方面。与传统有限单元和有限差分法等软件相比,其优越性在于能够较好地表征岩体几何非连续特征及模拟节理裂隙的张开、块体滑移、旋转、垮落等大位移、大变形分析。 本文结合淮南顾北煤矿深部岩巷支护试验研究,采用离散元数值模拟与现场监控量测相结合的方法分析研究该类巷道非线性大变形机制及支护设计方案,以期为今后开展的类似工程提供有益参考。 2 断裂围岩的非线性大变形机制分析 2.1 试验巷道选择及试验段选择 淮南顾北煤矿-648 m水平围岩以炭质泥岩、砂质泥岩及弱胶结的粉砂岩为主,属于大规模松软围岩。井底车场自开始建设时起,巷道就不断产生大变形而破裂失稳,衬砌破坏严重,底臌变形大,因而不得不多次重复返修巷道,从而大大延长了建设周期,增加了工程投资。为了从根本上解决该矿-648 m水平巷道围岩稳定控制和支护的技术难题,选择-648 m南翼回风大巷作为试验巷道进行现场支护试验。 整个南翼回风大巷全长约2 000 m,区域地质条件复杂,围岩裂隙、层理发育,且其中有一段需要穿过F92大断层,该断层附近又次生有4条小断层,断层带岩层被切割得非常破碎,自稳能力极差,支护难度很大。由于是新建巷道,在其中选择一定长度的试验段进行支护试验研究,为使成果具有针对性和代表性,试验段位于F92断层及其次生断层的影响带内。 通过岩体结构面现场地质调查,可确定巷道横断面上主要存在两组节理,其中一组节理倾角为120°,间距为0.8 m;另一组倾角约50°,间距为0.5 m。同时,地应力测试结果显示,该处垂直主应力为17.7~17.8 MPa,最大水平主应力为19.1~20.2 MPa,侧压力系数平均为1.10,整体属高地应力场区。 2.2 块体材料模型 为分析深部巷道破裂围岩变形破坏过程及非线性大变形机制,本文采用离散元数值分析软件UDEC4.0进行数值模拟计算。UDEC是Itasca公司推出的二维离散单元软件,在该软件中岩体被视为裂隙节理切割并处于镶嵌状态的碎裂块体,且块体和节理材料分别服从自身的应力-应变本构关系。其中块体代表岩石材料,它可以是刚体,也可以是变形体,而块体边界则代表节理结构面。若块体是可变形的,则块体进一步划分为三角形常应变有限差分网格,计算时通过运动微分方程的数值积分获得三角有限变形区上各个节点的位移增量,然后依据块体材料模型,得到块体内的应力值。 (1)可变形块体模型 由于岩土类材料的抗拉强度远不及其抗压强度,因此,本文可变形块体材料采用了考虑抗拉强度的Mohr-Coulomb弹塑性本构模型,即当块体承受的拉应力超过其抗拉强度时,块体发生拉破坏。其在应力空间的表示见图1。 图1中A到B由Mohr-Coulomb屈服函数fs定义为 式中:?为内摩擦角;c为黏聚力;N?=(1+sin?)/(1-sin?)。 BC段由拉应力屈服函数ft定义为 式中:σt为抗拉强度。 (2)节理材料模型 UDEC软件提供了几种节理材料模型,如点接触库仑滑移模型、面接触库仑滑移模型、连续屈服模型等,其中面接触库仑滑移模型最适于地下工程岩体的开挖模拟,通常的节理张开、剪切屈服及剪胀效应在此模型中都能实现。 节理面接触行为以弹簧-滑块来模拟,根据应力-位移关系,其作用力分为法向应力和切向应力。 法向应力作用在节理的法向方向,假设其应力-位移关系是线性,且由法向刚度系数Kn控制,则