岩样单轴拉伸应变局部化及全程应力-应变曲线

基于梯度塑性理论，考虑了应变局部化，提出单轴拉伸条件下岩样全程应力一应变曲线解析解。沿试样长度方向的总应变的微分被认为由两部分构成：可恢复的弹性应变的微分和不可恢复的塑性应变的微分。根据虎克定律，弹性应变的微分依赖于应力的微分及弹性模量。塑性应变的微分与应力的微分、试样高度、软化模量及局部化带的厚度有关，局部化带的厚度由梯度塑性理论确定。根据总应变的微分等于弹性及塑性应变的微分之和这一假设，在峰后线性应变软化本构关系的情形下，得到全程应力一应变曲线的解析解。通过与DeBorst及Muhlhaus基于梯度塑性理论得到的数值解对比，分别验证局部化带内部塑性应变分布解析解及内部长度对全程应力-应变曲线的影响。研究有关的本构参数（弹性模量及软化模量）及试样高度对全程应力-应变曲线的影响。     

加载速度对随机缺陷岩样破坏过程的影响

在单轴平面应变压缩条件下,采用FLAC模拟加载速度对具有随机材料缺陷岩样破坏过程的影响。采用编写的FISH函数规定随机缺陷及统计发生破坏的单元数。密实岩石服从莫尔–库仑剪破坏与拉破坏复合的破坏准则,破坏之后呈现应变软化–理想塑性行为。缺陷破坏之后经历理想塑性行为。随着加载速度的提高,强度及其对应的应变提高,峰后应力–变形曲线变得平缓。在加载过程中,每10个时步内破坏单元数–轴向应变曲线中存在破坏单元数有显著增加的3个区段。随着加载速度的提高,该曲线的区段2及3变得开阔,区段2的峰值降低,区段3的峰值提高。在初始加载阶段和缺陷全部破坏之前,加载速度较高时的破坏单元总数–轴向应变曲线比较平缓,这是由于加载速度较高时试样内部的裂隙传播和应力转移不充分,当应变小于一定值时,在应变相同时,破坏单元数较少。破坏单元总数–应变曲线表明,随着加载速度的提高,试样最终遭受到更严重的破坏。     

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

采动引起的断层错动是诱发矿震和断层冲击矿压的主要原因之一。在现有的数值模拟研究中，断层上测点的数量较少，以致于对断层各处力学状态全貌的了解较为困难。而且，针对煤层应力分布的研究通常仅限于垂直应力（支承压力），而对其它应力尚缺乏了解。为研究采动条件下断层附近应力时空分布，采用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时，断层上盘煤层垂直、水平和切应力受断层影响较显著的区域通常分别在距断...     

岩层运动并行计算系统StrataKing及应用

为了克服国外串行商业软件的功能有限、计算规模小和计算效率低的问题，历时10余年，自主研发了拉格朗日元与离散元耦合的岩层运动GPU并行计算系统StrataKing。介绍了Strata-King的计算流程、基本原理和功能，各计算模块均在GPU上执行。通过模拟三点弯梁的开裂过程，检验了StrataKing的正确性。在静水压力条件下，当单元数目较多时，巷道围岩的裂纹以更加弯曲、柔美而非突兀、生硬状向围岩深部发展，这与塑性力学中圆环形的各处均破坏的塑性区有所不同。在采动条件下，工作面煤壁超前支承压力和水平应力均具有峰值，后者与岩层之间的黏结作用有关，这与基于连续介质假定的厚壁筒径向应力特征有所不同。由此，StrataKing的大规模计算能力、丰富的功能和行业鲜明特色得到了进一步呈现。