B080302应用FLAC软件对煤柱稳定性影响因素的模拟.docVIP

应用FLAC软件对煤柱稳定性影响因素的模拟 张念超 杜祥海 张延水 【中国矿业大学，江苏 徐州 221008】 摘 要 应用二维程序FLAC3.3模拟了不同煤柱宽度下巷道的位移场与应力场，以及破坏范围，得出了煤柱尺寸与巷道变形的关系曲线。采用三维程序FLAC3D建立了大型数值计算模型，分析了煤柱宽高比对煤柱稳定性的影响及在采动系数为2时，不同煤柱宽度对煤柱稳定性的影响。 关键词 FLAC3.3 FLAC3D 煤柱尺寸 煤柱稳定性 ----------------------------------------------------------------------- 1 概 述 FLAC是国外近十几年兴起并发展起来的一种用于工程力学计算的显式有限差分程序。该程序可模拟土、岩石等材料的力学行为，适应于多种材料模式与边界条件的非规则区域的连续介质求解问题。在求解过程中，FLAC采用动力松弛法，不需要形成刚度矩阵，避免了直接求解大型联立方程组，便于在微机上求解较大规模的工程问题，与有限元法相比，FLAC在处理几何非线性和大变形问题上具有明显优势。FLAC程序中采用多种类型的力学本构模型，如摩尔一库仑弹塑性模型、应变硬化／软化模型、节理模型等，来模拟非线性、大变形及不可逆等地质材料的变形特征。FLAC还采用四种结构单元，包括梁、锚杆、桩及支柱单元，可以模拟支护结构。 FLAC程序采用的Lagrangian法是一种分析连续体非线性大变形问题的数值方法。这种方法遵循连续介质的假设，利用差分格式，按时步积分求解，采用拖带坐标系，随着结构形状的变化不断更新坐标，允许连续介质有大的变形。 FLAC的求解过程是：首先离散结构为若干单元(Zone)和若干网格点(Gridpoint)相互连接的连续网格体系。在此体系中，对单元与网格点应满足下列平衡方程： （1-1） 为获得公式(1-1)的真实解，在尽可能少的虚拟时间里，划分多个时步，按照图1步骤对网格中单元与网格点循环往复计算，直至不平衡力消失，各点的位移与应力分别趋向常数。 图1 FLAC 计算过程 2 模拟实验条件 为了全面、系统地了解不同尺寸煤柱的强度、变形和应力分布等影响煤柱稳定性的因素，可采用FLAC3.3（二维）、FLAC3D（三维）进行了模拟，以便合理优化设计煤柱尺寸。 本模拟实验选取具有代表性的山西潞安矿区常村矿S2-3工作面为具体实验对象，该工作面的技术条件如下： S2-3轨道顺槽沿3#煤层底板掘进，巷道埋深为345～365m。3#煤层平均厚度6.05m，抗压强度为13.3MPa。煤层直接顶为黑灰色泥岩，厚度4.4m，抗压强度为51.3 MPa；老顶为细砂岩和中砂岩，厚度4.5m。3#煤层直接底板为中细砂岩，厚度平均为3.03m。 原岩应力测量结果表明：最大水平主应力为14.01MPa，方向为北偏西35.7°；最小水平主应力为6.96MPa；垂直主应力为9.00MPa。巷道一侧为S2-2工作面采空区，煤柱宽度为5m。 巷道断面为矩形，宽4m，高3.2m，掘进断面积为12.8m2。巷道支护形式为锚、网、梁+锚索联合支护。顶板锚杆为直径22mm、长2400mm的左旋无纵筋螺纹钢筋，树脂加长锚固，锚杆排距为800mm，每排6根锚杆，间距700mm；锚索直径15.24mm，长8300mm，树脂加长锚固。每排一根锚索，在距小煤柱侧1.3m处垂直顶板，排距1.6m。巷帮锚杆为直径22mm长2000mm的左旋无纵筋螺纹钢锚杆，树脂加长锚固，靠工作面侧，锚杆间距800mm，每排4根锚杆；小煤柱侧锚杆间距为600mm，每排5根锚杆，锚杆排距为800mm。采用钢筋托梁与金属网护帮护顶。 3 FLAC3.3二维数值模拟 3.1 模拟方案设计 模拟该工作面巷道的顶板下沉量和两帮移近量分别在不受采动影响、煤柱宽度为30m、煤柱宽度为15m、煤柱宽度为10m、煤柱宽度为5m及煤柱宽度为3m时的变化情况。 3.2 模拟结果分析 煤柱宽度与巷道围岩变形的关系模拟结果如下表所列，绘制折线图如图2所示。 煤柱宽度与巷道围岩变形的关系 煤柱宽度（m） 30 15 10 5 3 顶板下沉量（mm） 284 341 349 305 277 两帮移近量（mm） 536 677 742 778 753 图2 煤柱宽度与巷道围岩变形的关系 由表和图2分析可得，当煤柱宽度为30m时，巷道受到明显的采动影响，导致顶板下沉量比不受采动影响增加67mm，两帮移近量增加129mm。随着煤柱尺寸减小，巷道变形逐渐增加。当煤柱宽度为15m时，顶板下沉又增加了57mm，两帮移近量增加了141mm。当回采工作面距巷道10m时，与煤柱宽度15m时相比，顶板下沉量增加了8mm，变化不大；两帮移近量增加64mm，其中煤柱侧帮位移增加