短壁开采煤柱稳定性的粘弹塑性分析解析 .ppt

短壁开采煤柱稳定性的粘弹塑性分析 作者：秦昊 李青锋 短壁开采技术是解决“三下”开采难题的主要途径和方法之一，而短壁开采中确定合理的采宽与留宽是控制岩层稳定性的关键。本文以五阳煤矿村庄下短壁开采为研究对象，针对保护煤柱的几何尺寸及其力学特性，采用ANSYS数值软件分析煤岩弹塑性本构关系下影响煤柱稳定性的各个因素，进而引入粘弹性本构关系分析时间效应下煤柱稳定性，并确定合理开采方案。分析所得结果可为开采实践提供指导性意见。 * [摘 要] [关键词] 短壁开采；弹塑性；粘弹性；数值模拟 * [煤柱强度与载荷的估计方法 ] Wilson提出了“两区约束理论”。该理论基于以下4点假设： (1)煤柱由核区和屈服区组成，已破坏的屈服区包围核区并对核区形成约束，核区处于三轴应力状态，大体符合弹性法则。 (2)煤柱的屈服应力为侧压 的 倍， 与煤体的内摩擦角有关： ，一般可取 。 (3)煤柱边缘的无约束垂直应力，屈服区水平约束力 由外往里渐增，至与核区交界面时为最大即等于原岩自重应力 ，屈服区宽度为 (4)一旦核区内部达到峰值应力，核区弹性状态将逐渐消失，煤柱将失稳。故煤柱的稳定性极限为核区平均应力 。 式中M，H为煤柱的高度和埋深 f f b sin 1 sin 1 tan - + = * [煤柱设计] 考虑现场各影响因素，依据煤柱设计的一般方法进行修改评价，煤柱合理设计核区率必须满足如下公式： 核区宽度必须满足如下公式： 式中 —煤柱塑性区宽度； —煤柱宽度。 * 数值计算模型 [煤柱稳定性的弹塑性数值分析] 边界条件为：模型两侧为滑动铰支约束，下表面固支约束，上表面自由并承受与覆岩厚度相应的垂直应力，即q=γH，γ为上覆岩层平均容重，H为采场至地表深度。考虑到减少边界效应，取5个采掘面进行数值模拟。 图1几何边界条件 图 2网格剖面 在计算中采用plane82单元，在重点分析部位如煤柱处设置较密的网格划分，其它区域的网格稀疏些，模型共7~8万个单元，具体网格剖分如图2所示。 * [煤柱稳定性的弹塑性数值分析] 数值计算方案 在弹塑性数值分析中采用如下3种计算方案： （1）保持所留煤柱宽度L=24m不变，改变开采长度a，分析弹塑性条件下煤柱的变形与煤柱间距变化之间的关系。分别取a=25m，30m，35m，40m，45m，50m六种计算模型； （2）开采长度a=25m不变，改变保护煤柱宽度L，分析弹塑性条件下煤柱的变形与煤柱宽度变化之间的关系。分别取L=15m，18m，21m，24m，27m，30m六种计算模型； （3）取a=25m，L=24m，改变煤岩的力学性质，分析煤柱的变形与煤岩性质之间的关系。分别取煤的弹模E=1GPa，2GPa，3GPa，4GPa，5GPa五种计算模型。 * （1）开采长度的影响 图3煤柱竖向侧面X方向位移变化 图4煤柱横向截面等效应力变化 图3从上至下依次对应开采宽度a=25m，30m，35m，40m，45m，50m六个计算结果。煤柱侧面的变形趋势是一致的，侧面向采空区一侧突出，但随开采宽度的增大，侧面突出的程度愈是显著。 图4由下至上，对应开采宽度a=25m，30m，35m，40m，45m，50m六个计算结果。横向截面应力在不同开采方案下变化整体趋势一致，两侧的应力值大于内部，长度的增大导致煤柱整体等效应力的提高；同时随长度的增大，煤柱横截面中两侧的等效应力大幅度提高，对应着煤柱两侧的塑性区域扩张，弹性区缩小。 [煤柱稳定性的弹塑性数值分析] * （2）煤柱留设宽度的影响 图5 不同煤柱宽度第二煤柱竖向侧面X方向位移变化 图5从下至上依次对应煤柱宽度L=15m，18m，21m，24m，27m，30m六个计算结果。煤柱侧面的变形趋势一致，侧面都是向采空区一侧突出；但随煤柱宽度的增大，突出的程度愈是微弱，变形减小。 在开采宽度不变的情况下，随煤柱长度增大，煤柱的应力集中现象得到很大的解决，整个煤柱的塑性区域不断收缩，煤柱整体Y方向的位移也受其影响。对于整个工作面而言，整个区域的应力值也是随其有降低的趋势，并直接导致了煤岩上覆各层的应力值变小，变形平缓，有利于采掘面的推进。 [煤柱稳定性的弹塑性数值分析] * （3）煤岩性质的影响 图6煤柱竖向侧面X方向位移变化