煤柱稳定性的弹塑性数值分析

《煤柱稳定性的弹塑性数值分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、短壁开采煤柱稳定性的粘弹塑性分析作者：秦昊李青锋短壁开采技术是解决“三下”开采难题的主要途径和方法之一，而短壁开采中确定合理的采宽与留宽是控制岩层稳定性的关键。本文以五阳煤矿村庄下短壁开采为研究对象，针对保护煤柱的几何尺寸及其力学特性，采用ANSYS数值软件分析煤岩弹塑性本构关系下影响煤柱稳定性的各个因素，进而引入粘弹性本构关系分析时间效应下煤柱稳定性，并确定合理开采方案。分析所得结果可为开采实践提供指导性意见。2[摘要][关键词]短壁开采；弹塑性；粘弹性；数值模拟3[煤柱强度与载荷的估计方法]Wilson提出了“两区约束理论”。该理论基于以下4点假设：(1)煤柱由核区和屈服

2、区组成，已破坏的屈服区包围核区并对核区形成约束，核区处于三轴应力状态，大体符合弹性法则。(2)煤柱的屈服应力为侧压的倍，与煤体的内摩擦角有关：，一般可取。(3)煤柱边缘的无约束垂直应力，屈服区水平约束力由外往里渐增，至与核区交界面时为最大即等于原岩自重应力，屈服区宽度为(4)一旦核区内部达到峰值应力，核区弹性状态将逐渐消失，煤柱将失稳。故煤柱的稳定性极限为核区平均应力。式中M，H为煤柱的高度和埋深ffbsin1sin1tan-+=4[煤柱设计]考虑现场各影响因素，依据煤柱设计的一般方法进行修改评价，煤柱合理设计核区率必须满足如下公式：核区宽度必须满足如下公式：式中—煤柱塑性区

3、宽度；—煤柱宽度。5数值计算模型[煤柱稳定性的弹塑性数值分析]边界条件为：模型两侧为滑动铰支约束，下表面固支约束，上表面自由并承受与覆岩厚度相应的垂直应力，即q=γH，γ为上覆岩层平均容重，H为采场至地表深度。考虑到减少边界效应，取5个采掘面进行数值模拟。图1几何边界条件图2网格剖面在计算中采用plane82单元，在重点分析部位如煤柱处设置较密的网格划分，其它区域的网格稀疏些，模型共7~8万个单元，具体网格剖分如图2所示。6[煤柱稳定性的弹塑性数值分析]数值计算方案在弹塑性数值分析中采用如下3种计算方案：（1）保持所留煤柱宽度L=24m不变，改变开采长度a，分析弹塑性条件下煤

4、柱的变形与煤柱间距变化之间的关系。分别取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五种计算模型。7（1）开采长度的影响图3煤柱竖向侧面X方向位移变化图4煤柱横向截面等效应力变化图3从上至下依次对应开采宽度a=25m，30m，35m，4

5、0m，45m，50m六个计算结果。煤柱侧面的变形趋势是一致的，侧面向采空区一侧突出，但随开采宽度的增大，侧面突出的程度愈是显著。图4由下至上，对应开采宽度a=25m，30m，35m，40m，45m，50m六个计算结果。横向截面应力在不同开采方案下变化整体趋势一致，两侧的应力值大于内部，长度的增大导致煤柱整体等效应力的提高；同时随长度的增大，煤柱横截面中两侧的等效应力大幅度提高，对应着煤柱两侧的塑性区域扩张，弹性区缩小。[煤柱稳定性的弹塑性数值分析]8（2）煤柱留设宽度的影响图5不同煤柱宽度第二煤柱竖向侧面X方向位移变化图5从下至上依次对应煤柱宽度L=15m，18m，21m，2

6、4m，27m，30m六个计算结果。煤柱侧面的变形趋势一致，侧面都是向采空区一侧突出；但随煤柱宽度的增大，突出的程度愈是微弱，变形减小。在开采宽度不变的情况下，随煤柱长度增大，煤柱的应力集中现象得到很大的解决，整个煤柱的塑性区域不断收缩，煤柱整体Y方向的位移也受其影响。对于整个工作面而言，整个区域的应力值也是随其有降低的趋势，并直接导致了煤岩上覆各层的应力值变小，变形平缓，有利于采掘面的推进。[煤柱稳定性的弹塑性数值分析]9（3）煤岩性质的影响图6煤柱竖向侧面X方向位移变化图7煤柱横向截面等效应力变化图6从下往上依次对应弹模E=1GPa，2GPa，3GPa，4GPa，5GPa五

7、个计算结果。煤柱侧面的变形大体是一致的，侧面都是向采空区一侧突出。但随弹模的增大，突出的程度是减弱的，变形的程度是在降低的。图7从下往上，对应弹模E=1GPa，2GPa，3GPa，4GPa，5GPa五个计算结果。截面应力在不同开采方案下变化整体趋势相对一致。但弹性模量增大导致煤柱两侧等效应力提高愈大，相反其内部等效应力却愈小。这是由于煤岩弹模增大导致煤柱两侧的应力相对集中，但变形不大，进而间降低了内部应力值大小。同时发现仅管等效应力在数值上相对有较大的增大趋势，但弹模改变对塑性区和弹性区的分布改变不大。