受复杂采动影响沿空掘巷煤柱变形破坏机理研究.pdfVIP

168 ^铡jlI舛枝 2013年第5+期 受复杂采动影响沿空掘巷煤柱变形破坏机理研究 柳泉 (充州煤业股份有限公司鲍店煤矿，山东鄄城273500) 摘要借助大型数值模拟软件，对沿空掘巷煤柱经受邻近工作面回采和自身工作面回采两次采动过程中的受力、变形时空演变规律进行数 值仿真研究，着重分析了煤柱应力、位移变化、超前应力影响范围等规律。说明煤柱在采动中是从燥柱边缘向煤柱内部逐步过渡破坏。研究结 果为工程应用提供了理论基础和技术指导。 关键词 沿空掘巷 复杂采动 破坏机理 数值模拟 中图分类号TD322文献标识码A doi：10．3969／j．issn．1005—2801．2013．05．99 本文基于ANSYS大型数值模拟软件，按照一次采 a一材料参数； 动、二次采动的顺序对煤柱的应力、变形时空演变规律 G一煤岩体弹性剪切模量； 进行数值仿真研究，分析了沿空掘巷煤柱应力、位移变 K一弹性体积变形模量。 化、超前应力影响范围以及煤柱变形破坏规律，尤其是 该模型服从广义的VonMises屈服条件，即，=甜， 搞清沿空掘巷煤柱的破坏演变机理，为今后合理留设 +缸一C=0，因此D—P准则比较适合岩土类材料。 煤柱，对工程应用提供了理论基础和技术指导。 2模拟结果分析 1数值分析模型 2．1 一次采动煤柱受力、位移动态演变特征分析 本文假设回采煤柱宽度为15m，高3m，工作面埋2．1．1 上顺槽掘进煤柱受力与位移特征分析 深约400m，上赋等效压应力约为6MPa。采用通用有 掘进过程中，煤体的应力、位移变化显著；上顺槽 限元软件ANSYS建立三维数值分析模型，模型长、宽、拐角发生应力集中现象，最大压应力达10．7MPa，且位 高分别为100m、120m和45rn。模型分上下临近工作 移变化较大，呈现以上顺槽起始位置为中心的半圆形 面，边界侧面采用法向位移约束，底面竖直方向和水平 轮廓。靠近已开挖上顺槽部分的煤柱应力、位移变化 方向位移均约束；顶面为自由面。模型煤(岩)体的物 显著，且有一段超前影响。为进一步探讨采动过程中 理力学参数如表1所示。按照一次采动、二次采动的 的沿空掘巷煤柱稳定性，在模拟中对上顺槽掘进24m、 顺序对煤柱的应力、变形动态演变过程进行了模拟。 45m、63nl段为研究对象，对煤柱内部应力、位移变化情 表1煤(岩)体物理力学参数 况进行了监测。监测线沿着顺槽掘进方向，取距开切 层厚 弹性模量E 泊松比 密度 内聚力 摩擦角 眼45m处的煤柱上平面，该煤柱侧向支承压力和垂直 岩层名称 (qPa) LL (MPa) (o) (m) (kz／m3) 位移曲线见图1。 砂质泥岩 4．75 19．7l 0．257 2535 4 16 粗砂岩 5．25 19．78 0．18 2400 15 45 (1)随上顺槽不断掘进，煤柱支承压力变化主要 中砂岩 7．52 35．26 0．16 2550 17．2 23 集中在距上顺槽侧向距离为0—1．5m范围内，煤柱在 细砂岩 8 24．38 O．17 2440 3．1 42 煤 5 0．94 0．303 1447 1．5 26 距上顺槽侧向距离1．5m处压应力最大，最大值约为 粉砂岩 14．48 30．38 o．15 舢 2．4 45 7．8MPa，塑性区宽度小于I．5m。随着距上顺槽侧向距 离的增加，煤柱内应力逐渐趋向于初始地应力 煤(岩)体材料采用Dracker—Prier弹塑性本构 模型，本构方程如式(1)所示。