急斜煤层开采石门巷道变形的数值模拟分析.pdfVIP

急斜煤层开采石门巷道变形的数值模拟分析 李 强，杜石文，陈小波，田献国，卢荣利 (中国矿业大学力学与建筑工程学院，江苏 徐州 221116) 直位移均约束，下部边界取U=V=0为全约束边界。 上部边界不约束 ，为 自由边界。上部边界以上的岩层 作为外载荷施加在模型的上边界上。 1。1．3计算模型网格图 在模拟计算中，将采区煤层回采分成采块2一l0依 次进行开采。石门巷道上方布置4个采块，下方布置4个 采块，每个采块垂直高度为 10m。中间采块为穿越煤层的 石门巷道所在的采块，石门巷道周围留有 2m厚煤层，该 采块垂直高度为7m。 0 引 言 计算模型网格图如图 1所示。 拉格朗日方法是岩土体力学行为分析的有效方法 ，可以准确地 模拟材料的屈服、塑性流动 、软化直至大变形 ，尤其在大变形分析以 及模拟施工过程等领域有其独到的优点。三维快速拉格朗日分析程 序FLAC3D具有强大的后处理功能。非常适用采矿工程中对地质体 的锚固分析等数值模拟。 快速拉格朗13分析程序有很强的分析功能，其基本思想是：通过 对三维介质的离散 ，使所有外力与内力集中于三维网络节点上，进而 将连续介质运动定律转化为离散节点上的牛顿定律；时间与空间的 导数采用沿有限空间与时间间隔线性变化的有限差分来近似计算； 将静力问题当作动力问题来求解，运动方程中惯性项用来作为达到 所求静力平衡的一种手段。 1 石门巷道数值模拟模型建立 1．1计算几何模型的确定 1．1．1计算模型范围的确定 石门巷道开挖后 ，包括顶底板在内的围岩原始应力要发生变化 ， 图1 模型网格图 产生变形、移动乃至破坏，其波及范围往往很大。按照圣维南原理 ，其 模拟过程中对每个采块开采结束，变形稳定后的位 影响范围仍有一定的限度。同时，考虑到便于数值模拟计算方便，按 移图进行分析，在加固的30m石门巷道里选取典型断 照实践经验 ，需要选取一个合适的影响范围。急斜煤层开采时，石门 面，在断面上设置数据提取点，分析其位移的变化 ，以了 巷道加固长度一般为30m，模拟煤层开采倾向长度 87m，考虑到空间 解煤层的开采对石门巷道加固段的影响，由此对石门巷 半无限体的对称性，计算模型模拟的空间长度 ×宽度 ×高度=70mX 道的支护结构进行完善。分析断面以及数据点的位置选 20m ×107m。 取见图2。 模型建立时首先模拟石门巷道围岩在 自重应力下的静力稳定 性 ，然后在此条件下进行煤层开采。为了提高计算速度，将石门巷道 从中心剖分，取其一侧进行计算。考虑计算机模拟速度及模拟结果的 C D 精确性 ，将石门巷道附近以及开采煤层顶底板的网格划分适当加密， 其中开采煤层部分采用空单元进行模拟。 1．1．2边界约束条件的确定