毕设矿压规律.docVIP

3.综放仰斜开采矿压显现规律研究 3.1 矿压显现规律数值计算 放顶煤综合机械化开采作为一种新的开采厚煤层的方法，目前还不能很好用一般的力学解析方法求解其力学变化与破坏过程。为了研究和掌握开采过程中顶板的力学变化过程及对工作面的力学作用等，运用合理的数值方法与软件进行计算和分析以及预测是目前行之有效的方法之一。本次研究采用计算软件FLAC3D 对工作面进行计算分析，其目的在于：通过上述研究获得在不同的推进度下支承压力的的破坏区域。 1.1 有限差分软件FLAC 3D简介 FLAC(Fast Lagrangian Analysis of Continue)是为地质工程和岩土工程应用而开发的连续介质显式有限差分计算机软件，主要适用于模拟计算岩土类工程和地质材料的力学行为，特别是岩土材料达到屈服极限后产生的塑性流动，材料通过单元和区域表示，根据研究对象的形状构成相应的网格结构。每个单元在外载和边界的约束条件下，按照约定的线性和非线性应力－应变关系产生力学响应。FLAC 软件建立在拉格朗日算法基础上，介于有限元和离散元之间，特别适用于模拟地质材料的大变形和扭曲转动。FLAC 程序设有多种本构模型，可模拟地质材料的高度非线性（包括应变软化和硬化）、不可逆剪切破坏和压密、粘弹性（蠕变）、孔隙介质的流固耦合、热力学耦合以及动力学行为等。另外，程序还设有边界单元，可以模拟断层、节理和摩擦边界的滑动、张开和闭合行为。支护结构，如衬砌、锚杆、可缩性支架或板壳等与围岩的相互作用也可以在FLAC 程序中进行模拟。同时用户可根据需要在FLAC 中创建自己的本构模型进行各种特殊的修正和补充。FLAC 程序具有强大的后处理功能，用户可以根据需要可以直接在屏幕上绘制或以文件的形式创建和输出打印多种形式的图形，使用者还可以根据需要，将若干个变量合并在同一个图形中进行研究和分析。由于FLAC 程序主要为地质工程应用而开发的岩石力学数值计算程序，程序包括反映地质材料力学行为效应的特殊数值。FLAC 设计有七种材料本构模型： 各向同性弹性材料模型； 横观各向同性弹性材料模型； 莫尔－库仑弹塑性材料模型； 应变软化、硬化塑性材料模型； 双屈服塑性材料模型； 节理材料模型； 空单元模型，可用来模拟地下开挖和煤层开采。 FLAC 采用显式算法来获得模型的全部运动方程（包括内变量）的时间步长解，从而可以追踪材料的渐进破坏直至整体垮落，这对研究采矿设计是非常重要的。此外，程序允许输入多种材料模型，也可以在计算过程中改变某个局部变量的材料参数，增强了程序使用的灵活性，极大地方便了模拟计算时的处理。基于上述计算功能和材料模型，FLAC 程序比较适合于采矿工程和地下工程的分析和设计。FLAC3D 程序的基本原理和算法与离散元相似，但它却像有限元那样适用于多种材料模式与边界条件的非规则区域的连续问题求解；在求解过程中，FLAC3D 采用了离散元的动态松弛法，不需要求解大型联立方程组（刚度矩阵）。另一方面，同以往的差分分析方法相比，FLAC3D 不但可以对连续介质进行大变形分析，而且能模拟岩体沿某一软弱面产生的滑动变形。FLAC3D 还能针对不同的材料特性，使用相应的本构方程来较真实的反映实际材料的动态行为。此外，该方法还可考虑锚杆、挡土墙等支护结构与围岩的相互作用。FLAC3D 在模型计算过程中，先是用运动方程由应力和力求出新的速度和位移。然后由速度推导出应变率，并由此求得新的应力。一个时步为沿回线的一个循环，如图-1 所示。 计算模型及参数确定 （1）以 工作面地质和开采技术条件为背景200m、100m和140m，共划分282000个六面体单元，生成网格节点295728个。X轴平行于工作面推进方向，采用Mohr-Coulomb本构模型。模型水平方向上左右两侧限制水平位移，模型底面限制水平移动和垂直移动，模型的上部施加上覆岩层的自重应力，建立模型图1。 （2）计算参数确定 根据现场取样和岩石力学试验结果，当载荷达到强度极限后，岩体产生破坏，在峰后塑性流动过程中，岩体残余强度随着变形发展逐步减小。因此，计算中采用莫尔—库仑（Mohr-Coulomb）屈服准则判断岩体的破坏；采用应变软化模型以反映煤体破坏后随变形发展残余强度逐步降低的性质 （3.1） 式中，σ1、σ3 分别是最大和最小主应力，c 和分别是粘结力和摩擦角。代表材料发生破坏的临界值，当0 时，材料将发生剪切破坏。在通常应力状态下，岩体的抗拉强度很低，因此可根据抗拉强度准则（σ3≥σT）判断岩体是否产生拉破坏。 根据现场地质调查和相关研究提供的岩石力学试验结果，考虑岩石的尺度效应，模拟计算采用的岩体力学参数表。 的垂直应力随深度线性变化，根据工作煤层埋藏深（H=- m）和