地震作用下巷道破坏模拟与实验对比分析.docVIP

地震作用下巷道破坏模拟与实验对比分析 　　摘要：基于工程结构波动及结构动力学理论，建立煤矿采空区岩层的动力学方程，将小型室内相似振动台实验结合ANSYS有限元数值计算，对比分析巷道的地震动力破坏特征。研究发现：在地震作用下，在拱形和圆形巷道的几何中心水平线正负45°区域以及拱形巷道底角处的剪应力、主应力最显著，而在拱形巷道的两底角及顶板处易产生破坏，圆形巷道的顶板产生了裂纹，圆形巷道的最大主应力及剪应力分布偏向右侧；拱形巷道最大主应力偏向右侧，剪应力呈对称分布，拱形巷道的两底角及顶板易产生破坏。在巷道的抗震设计中，可对巷道受力薄弱处加强加固和采取抗震措施，提高巷道的抗震性能。 　　关键词：巷道；地震波；ANSYS；相似材料实验 　　中图分类号：P315.9文献标识码：A文章编号：1000-0666（2017）04-0661-07 　　0引言 　　煤炭是我国主要能源之一，分别占能源生产和消费总量的76%和69%，是我国国民经济和社会发展不可缺少的物资基础。我国煤炭资源丰富，探明储量达10 202亿吨，其中可开采储量1 891亿吨；目前我国的煤炭产业仍然处在“高开采、低利用、高排放”的状况之中（薛冰，2012）。煤矿不断开采致使煤矿不断深入地下，煤矿巷道的安全问题显得尤为重要。一般情况下，岩体地下工程由于受到岩体介质强约束作用，其抗震性能明显优于地上结构（Lanzano et al，2008）。尽管如此，国内外的震害调查统计发现，由于地震的作用，地下巷道也出现了损伤、破坏。陈晓祥和韦四江（2008）基于柴里矿二水固南轨上山实际的地质情况，应用UDEC数值模拟分析了初始地应力对巷道围岩稳定性的影响。吴多华等（2016）采用FLAC3D数值模拟分析了地震荷载对巷道稳定性的影响，指出地震作用会使巷道发生低臌和冒顶现象。姜耀东等（2005）研究了放炮震动下巷道围岩失稳机理，指出放炮震动波导致围岩产生裂纹并诱发了煤层和顶底板间滑动。马行东等（2005，2006）、李海波等（2005）通过研究西部强震区地震波的不同参数对地下洞室的动力响应，得出地震波的入射方向、空间不均匀性对地下洞室的动力响应的影响主要与地下洞室的埋置深度、地应力等因素有关。 　　目前我国关于煤矿巷道的研究多是静力荷载作用下的巷道支护方面或静载作用下巷道的破坏方面的研究。即使是动荷载作用下，巷道的响应研究多是采用数值模拟，而相似材料振动台试验方面的研究较少。本文结合ANSYS数值模拟计算分析和相似材料试验结果对比分析，对地震作用下巷道的响应问题进行研究。 　　1地震作用下煤矿巷道-围岩结构动力响应的分析方法 　　在地震波作用下煤矿巷道围岩介质在水平方向（x方向）上会发生一定的位移，假设在地震波作用下煤矿采空区围岩所产生的应力为σ，应变为ε，位移为μ（魏晓刚等，2016）。 　　在围岩介质的弹性变形阶段，岩体内纵向地震波传播的波动方程为：2u′t2=c202u′x2（1）式中：u′=u-u0为地震波引起围岩介质发生的位移；c0为地震波在岩体中的传播速度。 　　由工程结构波动理论可知，在围岩介质的岩体中c20=λ+2μρ，λ、μ均为拉梅常数，可以判断地震波的传播速度与岩体的材料性质密切相关。 　　对（1）式进行求解可以得到：μ′（x，t）=f（x-c0t）+g（x+c0t）（2）式中：f（x-c0t）为地震波入射纵波的波动方程，g（x+c0t）为地震波反射纵波的波动方程。地震作用下在围岩介质的弹性变形阶段，分别求解地震波作用下煤矿巷道围岩在水平方向上所产生的应力σ、应变ε、位移u的数值解： 　　对上式进行积分求解，可以得到地震作用下煤矿采空区围岩结构体系的应变能为：U=∫0∫A∫0∫ετσdεdA（4）通过分析煤矿巷道围岩在地震作用下所产生的应力、应变、位移可以较为直观地判断煤矿巷道-围岩介质的动力稳定性。由于部分煤矿巷道的埋置深度较大，此时煤矿巷道及围岩介质所承受的矿山岩层压力加大。在外界动力荷载的扰动下，即使是较小的扰动荷载，围岩介质所承受的应力都有可能超过自身的屈服应力，导致围岩进入塑性损伤状态，此时围岩发生动力失稳破坏的可能性较大。 　　2相似材料实验 　　2.1实验设计材料配合比参照了康希并和张建义（1988）相似材料的配比，砂∶石灰∶石膏∶水=5∶0.5∶0.5∶0.6，如图1所示。拟定模型的几何相似比为1/20，容重比为0.7，（模型容重在1.5～1.8g/cm3间较为合适），弹性模量相似比为1/20。抗压强度为1/20。选择的实际岩体抗压强度为20 MPa，则模拟材料的强度约为1.0 MPa。 　　模具设计为40 mm×40 mm×100 mm的长方体相似材料试块，磨具中心为巷道，巷道分两种形式，分别为半径为4 cm的圆形巷道和底宽