采动影响下覆岩破坏过程分析.docxVIP

采动影响下覆岩破坏过程分析 1 总结 2.1 破坏单元的检查和测量 sppa2系统是一个值模拟工具，可以模拟岩石梁的变形、破坏和坠落过程。与其他通常的破坏模型一样，smpa2包括两个方面：电压分析和破坏分析。SFPA2D的应力分析采用有限元法,破坏分析则根据修正的库仑(Coulomb)准则(包含拉伸截断Tension cut-off)来检查材料中是否有破坏单元。对破坏单元则采用刚度特性退化(处理分离)和刚度重建(处理接触)的办法进行处理。采用空单元法模拟采动工作面的变化,令挖去部分中单元的自重和弹性参数为一个无意义的数。每一步开挖即相当于一次卸载,首先进行应力计算,然后根据破坏准则检查模型中是否有破坏单元。如果没有,继续下一步开挖,并作应力计算。如果有破坏单元,则根据单元的拉或剪破坏状态进行刚度退化处理,然后重新进行当前步的应力计算。重复上述过程,直至岩梁形成宏观裂纹,并最终冒落而形成整个覆岩破坏的动态发展过程。 2.2 数值模拟方法 地质学家把岩层定义为由同一岩性组成的、有两个平行或近于平行的界面所限制的层状岩石,即每一岩层具有相似的物理力学性质。若不考虑断层等大型地质构造的影响,则不连续面中对地表沉陷贡献最大的是层面,而且层面延展性强,分布范围广,因此可仅考虑分层效应。虽然岩体的每一个岩层具有复杂的力学性质,但许多文献表明,仅考虑弹性或者弹脆性,而不考虑弱面,二者结果相近,但与现场实测资料相差较大;若考虑弹性和弱面与考虑弹脆性和弱面,二者结果相近,且与现场实测资料吻合。开采引起的地表下沉与工程岩体尺寸相比可视为小变形,据此,在过去的研究中将每一岩层视为均质、各向同性线弹性体。然而,岩石具有明显的非线性本构,将整个岩层视为均质、各向同性的线弹性体有失准确。建立在统计方法理论基础上的数值模拟方法被看成是解决这一问题的有效途径,其特点在于引入简单的基元本构模型,同时运用统计的方法来反映由简单基元构成的整体的复杂力学行为。在这方面已有不少有益的成果,如DLNM模型、Rule-based模型,Micromechanics 模型,PFC模型等。他们均认识到材料的非均匀性对宏观力学性质的影响,并通过物理统计的方法在各自的模型中引入了非均匀性。 在本文的数值模拟中,设定岩层中每一个基元为均质、各向同性线弹脆性体,每一岩层为非均匀的,其岩石的力学性质(单元的弹性模量、强度等)服从韦伯分布φc(m,u),其中参数m为形状参数,反映岩石材料力学性质的均质度,m越大,表明岩石的性质越均质;u为反映岩石材料平均性质的参数,从而使每一岩层具有相似的力学性质。将开挖工程岩体介质概化为具有开挖孔洞的分层有相似力学性质的非均质各向同性的非线性空间体,如图1所示,其范围为水平方向无限大,垂直方向上至地表,向下延伸到煤层底部某一深度的厚坚硬岩层或对下位岩层的影响可以忽略的某一岩层处。 本文采用平面应变模型来模拟某矿覆岩的垮落过程。模型沿水平方向取300 m,沿垂直方向取100 m。煤层厚6 m,赋存深度84 m。模型共计12层,其主要岩体力学参数见表1所示。共划分为30 000单元,采用长壁式开采,开采强度5 m/d,垮落法顶板管理。具体的数值模拟模型和参数分布见图2所示。图中每一层的灰度代表岩层力学参数(弹性模量、抗压强度等)的大小,灰度越高,值越大。 3 数值模拟结果 图3是采动影响下覆岩垮落的动态发展过程的数值模拟结果,因篇幅有限,仅将其中具有代表性的7步弹性模量分布图和工作面推进过程中应力场的动态演化图的对比结果示于图中。 3.1 工作面顶次破断采空区 (1)直接顶冒落 开切眼形成以后,上覆岩层悬露,直接顶受开挖而引起的卸荷膨胀变形的影响,首先在中部产生离层,随着采场工作面的推进,在重力的作用下弯曲。岩梁悬露的跨度达到一定的跨度后,弯曲沉降发展到一定的限度,岩梁的端部首先开裂,当工作面推进到20 m左右,在岩梁的中部开裂形成“假塑性岩梁”,如图3(a)所示。当岩梁的沉降值超过“假塑性岩梁”允许沉降值时,悬露岩层即自行冒落,冒落形态存在不对称性,冒落高度为2 m,以后随着工作面的继续推进,直接顶随采随冒。如图3(b)所示。 (2)基本顶初次破断 随着采场工作面的进一步推进,基本顶在重力的作用下弯曲下沉。当工作面推进到45 m时,如图3(b)所示,基本顶出现大范围移动,采空区前后煤壁上方基岩的上部、中部出现拉裂隙。当工作面推进到50 m时,如图3(c)所示,基岩基本顶初次垮落,垮落形态表现出不对称性,基本顶破断岩块沿工作面煤壁切落,基本顶失稳运动有回转切落的特点。 (3)基本顶的周期性破断 基本顶初次破断后,当工作面推进到60 m时,如图3(d)所示,上部顶板出现离层、端部中部断裂现象。当工作面推进到70 m时,上位顶板在靠近工作面侧的1/4顶板跨度处又