三维连续-非连续并行计算方法及其在岩爆过程模拟中的应用

随着深部岩石工程的发展,岩爆变得越发严重。在岩爆的数值模拟方面,连续方法和非连续方法均具有一定的局限性。兼具二者优势的连续-非连续方法更具优势,且正在快速发展。基于CUDA对自主开发的三维拉格朗日元与离散元耦合连续-非连续方法进行了GPU并行加速。为了探索岩爆的机理和过程,模拟了不同静水压力、侧压系数和单元数目（最多达100×10⁴）条件下圆形洞室围岩V形坑的演化规律和单元弹射现象。考察了洞室围岩中裂纹的定量演化规律。研究表明：当静水压力较大时,基于芬纳公式的支护设计偏于不安全。由于V形坑的位置发生改变,V形坑的平均最大深度随着静水压力的增加先缓慢增加后快速增加。关于洞室围岩V形坑的模拟结果能与有关的实验结果、数值结果和现场观测结果吻合。上述研究很好地体现了岩爆并行计算较串行计算和商业软件计算的优势。     

正断层上盘开采断层附近应力时空分布数值模拟

采动引起的断层错动是诱发矿震和断层冲击矿压的主要原因之一。在现有的数值模拟研究中，断层上测点的数量较少，以致于对断层各处力学状态全貌的了解较为困难。而且，针对煤层应力分布的研究通常仅限于垂直应力（支承压力），而对其它应力尚缺乏了解。为研究采动条件下断层附近应力时空分布，采用FLAC3D建立了正断层上盘开采的数值模型。模型中共包括13个岩层和1个45°倾角的正断层，通过计算获得了工作面从断层上盘推进过程中断层各处及附近煤层应力的时空分布规律。监测了断层上盘断层上的22个节点以获取断层各处力学状态全貌。研究发现，在工作面推进过程中，断层上存在1个压紧区及1个或多个松动区；断层上的压紧区较为安全，压紧区上、下方松动区的正应力越来越低，当工作面邻近断层时该区的断层区段多趋于危险，该区范围最大为78 m；断层上盘煤层紧邻工作面煤壁的垂直、水平和切应力峰值分别位于煤壁前方7.5～10.5、10.5～15.5和11.5～12.5 m；断层下盘煤层3种应力最大值均上升或有上升趋势。当工作面煤壁距断层的水平距离由40 m减小至20 m时，断层上盘煤层垂直、水平和切应力受断层影响较显著的区域通常分别在距断...