第1章 实验研究背景

  本实验聚焦于煤矿井下煤层卸压与冲击地压防治这一核心工程领域。在众多常用的冲击地压防治技术中,钻孔卸压技术凭借施工工艺简便、设备要求较低、成本可控且卸压针对性强等优势,被广泛应用于各类煤层条件的卸压工程,成为保障煤矿安全生产的核心主流技术之一。

  现有学术研究与工程实践充分表明,钻孔布置参数的合理性直接决定卸压效果的优劣。国内外相关学者围绕钻孔卸压技术开展了大量研究,重点聚焦钻孔直径、钻孔深度、钻孔间距等单一参数的影响规律,通过室内物理试验、现场工业性试验等多种手段,形成了一系列可供参考的参数设计标准。但在实际工程应用中,钻孔排布形式(如线性排布、对称排布、交错排布等)对煤体应力传递路径、力学性能劣化规律的作用机制更为复杂,当前针对该方向的系统性研究仍存在显著不足。

  基于上述研究现状与工程实际痛点,本设计以“钻孔排布形式与煤体强度的关联规律”为核心研究方向,依托数值模拟技术开展系统性实验研究。通过精准构建不同钻孔排布组合的煤体试件模型,重点探究不同排布形式对煤体应力应变特性、应变场分布形态及塑性区演化规律的影响机制,量化不同排布形式下煤体峰值强度、弹性模量等关键力学参数的演化特征,最终明确兼顾卸压效果与煤体稳定性的最优钻孔排布方案。

第2章 实验方法介绍

  本实验核心采用Abaqus数值模拟方法开展研究,该方法的理论基石为材料力学理论与有限元分析原理,二者共同支撑起虚拟场景下煤体力学响应的精准模拟。从技术适配性来看,煤体在宏观受力分析尺度下可视为连续介质,其内部应力、应变的传递与分布符合材料力学的基本假设,这为数值模拟的适用性提供了理论前提。具体而言,该方法借助Abaqus专业工程模拟软件,先构建与实际煤体试件物理力学特性(如密度、强度、变形特性等)吻合的模型,再通过有限元离散化技术,将完整的煤体模型剖分为大量相互关联的微小单元体,实现对整体煤体在外部荷载作用下全过程力学响应的动态模拟。

  该数值模拟方法的核心实施流程与基本原理可细化为三个关键步骤:第一步是建立实体模型。需先通过查阅文献及同类煤体基础试验数据,获取目标煤体的关键物理力学参数,包括弹性模量、泊松比等,随后基于这些参数建立合适的实体模型。第二步是添加边界条件与载荷。煤柱在轴向压缩的荷载下,一方面,在煤柱上表面收到均布的轴向压力,另一方面,在煤柱的下表面需要限制竖直方向的位移与水平面两个方向的转角。第三步是数据求解与结果可视化,借助模拟软件的核心计算模块,基于算法求解各单元体的应力、应变等关键力学参数,再通过后处理模块将抽象的数值结果转化为直观的应力、应变、位移场云图等可视化成果,为后续分析提供清晰的数据支撑。

第3章 模拟实验方案

  模型尺寸:本实验采用模拟软件构建标准长方体煤体试件模型,尺寸精准设定为50mm×50mm×100mm(长×宽×高)。中间布置至少2个钻孔,孔径8mm。

  本构模型:煤体作为典型的脆性岩石类材料,在受力过程中会经历弹性变形。基于此,本实验选用弹性本构模型,具体设定如下:弹性模量2.5GPa、泊松比0.3等,上述参数组合可确保模型精准复刻目标煤体的力学特性。

  边界条件:对模型底部施加全约束,模拟煤体下部岩层的支撑作用;模型顶部不施加约束,设为自由受力面,为后续竖向加载提供受力空间。该边界条件的设计核心是还原煤体在地下开采过程中的单向受力状态(主要承受上覆岩层产生的竖向地应力,侧向应力相对较小),有效排除侧向额外应力对实验结果的干扰,确保模拟场景与工程实际受力环境的一致性,为后续加载试验的准确性提供保障。

  加载方式:结合煤矿井下地应力缓慢施加的实际工况,本实验采用单轴压缩加载方式。加载方向沿模型竖向,大小为1MPa的压强载荷。与煤体实际承受的竖向地应力方向一致。通过该加载方案,可获取煤体在受压过程中的完整荷载—位移曲线,为后续力学参数分析提供基础数据。

  试验变量:为精准探究钻孔排布形式对煤体强度的影响规律,本实验采用控制变量法,将钻孔排布形式设为唯一试验变量,其余参数(如钻孔直径、模型尺寸等)均保持一致。结合工程实际中常用的钻孔布置形式,共设计5种典型钻孔排布组合,具体方案如下:

  本文中使用的Abaqus版本是2025。

3.1 方案1(钻孔1×2)

  2个钻孔沿模型竖直中轴线分布,钻孔中心间距为30mm,距离上下边界为35mm;

  方案1 Abaqus源文件点此处下载。

图3.1 方案1(钻孔1×2)试件截面草图

图3.1 方案1(钻孔1×2)试件截面草图

  用Abaqus进行仿真,绘制的应力、应变、位移云图如下。

图3.2 方案1(钻孔1×2)试件应力云图(等轴视图)

图3.2 方案1(钻孔1×2)试件应力云图(等轴视图)

图3.3 方案1(钻孔1×2)试件应力云图(前视图)

图3.3 方案1(钻孔1×2)试件应力云图(前视图)

图3.4 方案1(钻孔1×2)试件应变云图(等轴视图)

图3.4 方案1(钻孔1×2)试件应变云图(等轴视图)

图3.5 方案1(钻孔1×2)试件应变云图(前视图)

图3.5 方案1(钻孔1×2)试件应变云图(前视图)

图3.6 方案1(钻孔1×2)试件位移云图(等轴视图)

图3.6 方案1(钻孔1×2)试件位移云图(等轴视图)

图3.7 方案1(钻孔1×2)试件位移云图(前视图)

图3.7 方案1(钻孔1×2)试件位移云图(前视图)

3.2 方案2(钻孔1×3)

  3个钻孔沿模型竖直中轴线分布,钻孔中心间距为30mm,距离上下边界为20mm;

  方案2 Abaqus源文件点此处下载。

图3.8 方案2(钻孔1×3)试件截面草图

图3.8 方案2(钻孔1×3)试件截面草图

  用Abaqus进行仿真,绘制的应力、应变、位移云图如下。

图3.9 方案2(钻孔1×3)试件应力云图(等轴视图)

图3.9 方案2(钻孔1×3)试件应力云图(等轴视图)

图3.10 方案2(钻孔1×3)试件应力云图(前视图)

图3.10 方案2(钻孔1×3)试件应力云图(前视图)

图3.11 方案2(钻孔1×3)试件应变云图(等轴视图)

图3.11 方案2(钻孔1×3)试件应变云图(等轴视图)

图3.12 方案2(钻孔1×3)试件应变云图(前视图)

图3.12 方案2(钻孔1×3)试件应变云图(前视图)

图3.13 方案2(钻孔1×3)试件位移云图(等轴视图)

图3.13 方案2(钻孔1×3)试件位移云图(等轴视图)

图3.14 方案2(钻孔1×3)试件位移云图(前视图)

图3.14 方案2(钻孔1×3)试件位移云图(前视图)

3.3 方案3(钻孔2×1)

  2个钻孔沿模型水平中轴线分布,钻孔中心间距为30mm,距离左右边界为10mm;

  方案3 Abaqus源文件点此处下载。

图3.15 方案3(钻孔2×1)试件截面草图

图3.15 方案3(钻孔2×1)试件截面草图

  用Abaqus进行仿真,绘制的应力、应变、位移云图如下。

图3.16 方案3(钻孔2×1)试件应力云图(等轴视图)

图3.16 方案3(钻孔2×1)试件应力云图(等轴视图)

图3.17 方案3(钻孔2×1)试件应力云图(前视图)

图3.17 方案3(钻孔2×1)试件应力云图(前视图)

图3.18 方案3(钻孔2×1)试件应变云图(等轴视图)

图3.18 方案3(钻孔2×1)试件应变云图(等轴视图)

图3.19 方案3(钻孔2×1)试件应变云图(前视图)

图3.19 方案3(钻孔2×1)试件应变云图(前视图)

图3.20 方案3(钻孔2×1)试件位移云图(等轴视图)

图3.20 方案3(钻孔2×1)试件位移云图(等轴视图)

图3.21 方案3(钻孔2×1)试件位移云图(前视图)

图3.21 方案3(钻孔2×1)试件位移云图(前视图)

3.4 方案4(钻孔2×2)

  4个钻孔于模型中心2×2分布,钻孔中心间距为30mm,距离左右边界为10mm,距离上下边界为35mm;

  方案4 Abaqus源文件点此处下载。

图3.22 方案4(钻孔2×2)试件截面草图

图3.22 方案4(钻孔2×2)试件截面草图

  用Abaqus进行仿真,绘制的应力、应变、位移云图如下。

图3.23 方案4(钻孔2×2)试件应力云图(等轴视图)

图3.23 方案4(钻孔2×2)试件应力云图(等轴视图)

图3.24 方案4(钻孔2×2)试件应力云图(前视图)

图3.24 方案4(钻孔2×2)试件应力云图(前视图)

图3.25 方案4(钻孔2×2)试件应变云图(等轴视图)

图3.25 方案4(钻孔2×2)试件应变云图(等轴视图)

图3.26 方案4(钻孔2×2)试件应变云图(前视图)

图3.26 方案4(钻孔2×2)试件应变云图(前视图)

图3.27 方案4(钻孔2×2)试件位移云图(等轴视图)

图3.27 方案4(钻孔2×2)试件位移云图(等轴视图)

图3.28 方案4(钻孔2×2)试件位移云图(前视图)

图3.28 方案4(钻孔2×2)试件位移云图(前视图)

3.5 方案5(钻孔2×3)

  6个钻孔于模型中心2×3分布,钻孔中心间距为30mm,距离左右边界为10mm,距离上下边界为20mm。

  方案5 Abaqus源文件点此处下载。

图3.29 方案5(钻孔2×3)试件截面草图

图3.29 方案5(钻孔2×3)试件截面草图

  用Abaqus进行仿真,绘制的应力、应变、位移云图如下。

图3.30 方案5(钻孔2×3)试件应力云图(等轴视图)

图3.30 方案5(钻孔2×3)试件应力云图(等轴视图)

图3.31 方案5(钻孔2×3)试件应力云图(前视图)

图3.31 方案5(钻孔2×3)试件应力云图(前视图)

图3.32 方案5(钻孔2×3)试件应变云图(等轴视图)

图3.32 方案5(钻孔2×3)试件应变云图(等轴视图)

图3.33 方案5(钻孔2×3)试件应变云图(前视图)

图3.33 方案5(钻孔2×3)试件应变云图(前视图)

图3.34 方案5(钻孔2×3)试件位移云图(等轴视图)

图3.34 方案5(钻孔2×3)试件位移云图(等轴视图)

图3.35 方案5(钻孔2×3)试件位移云图(前视图)

图3.35 方案5(钻孔2×3)试件位移云图(前视图)

第4章 总结

  提取各方案上表面竖直方向的位移,并与无孔时的解析解进行比较,见表3.1。

表3.1 各方案上表面竖直方向的位移

方案 轴向位移
无孔 -4.00E5
1x2 -4.13E5
1x3 -4.17E5
2x1 -4.19E5
2x2 -4.45E5
2x3 -4.69E5

  由表3.1可以看出,打孔可以有效减小轴向压缩时的位移。当孔的列数相同时,孔的行数越多,减小的轴向压缩时的位移越多,进而煤块的强度越强。当孔的行数相同时,孔的列数越多,减小的轴向压缩时的位移越多,进而煤块的强度越强。