1、第 卷 第 期 年 月自 然 灾 害 学 报 .收稿日期:修回日期:基金项目:国家自然科学基金项目()辽宁省“兴辽英才计划”项目():()“”()作者简介:李利萍()女教授博士主要从事矿山灾害力学方面的研究:.文章编号:()./.动载扰动下倾斜煤岩组合块体超低摩擦效应研究李利萍李明会潘一山(.辽宁工程技术大学 力学与工程学院辽宁 阜新 .辽宁大学 环境学院辽宁 沈阳)摘 要:随着我国煤矿相继进入深部开采三高一扰动复杂应力环境使深部煤岩体结构趋于破碎为块体煤岩块体界面多为含倾角结构面 冲击扰动作用下煤岩块体易沿煤岩界面产生失稳滑移诱发超低摩擦型冲击地压 利用有限差分软件建立砂岩煤砂岩组合块体以工
2、作块体水平位移、加速度幅值作为块体超低摩擦效应表征参数模拟不同倾角的煤岩组合块体受扰动冲击作用时工作块体超低摩擦效应特征参数的变化情况 考虑倾角大小、加载方式、冲击幅值和水平冲击延迟时间等因素影响得到含倾角煤岩组合块体超低摩擦效应影响规律 研究结果表明:动载扰动下含倾角煤岩组合块体煤岩块体倾角越大加速度峰值越大位移错动越明显超低摩擦效应越显著 随垂直扰动幅值的增加工作块体加速度振荡越剧烈趋于平衡收敛时间延长垂直扰动始终是诱发超低摩擦效应的主要因素 垂直扰动与水平冲击共同作用下模型底部、中部、顶部监测点超低摩擦效应强度依次增加随着倾角的增加位移差值不断扩大 工作块体位移振荡剧烈程度随水平冲击延迟
3、周期增加而增大在同一位移振荡周期内波峰点水平冲击比均值点与波谷点水平冲击更易发生超低摩擦效应关键词:超低摩擦效应煤岩组合含倾角块体动载扰动水平位移加速度幅值中图分类号:文献标识码:(.):.第 期李利萍等:动载扰动下倾斜煤岩组合块体超低摩擦效应研究 .:引言深部煤岩体结构在高地应力、高地温、高渗透压的复杂地质环境中逐渐破碎成块体 超低摩擦效应是指强冲击扰动作用于深部岩块介质系统时由于垂直扰动周期性的作用方式块体间压密程度随时间不断发生变化当岩块间处于疏松状态时其摩擦力减至最小此时在块系岩体结构接触面极易发生滑动失稳从而诱发超低摩擦型冲击地压对于超低摩擦效应的研究 等通过建立分析树脂玻璃块系岩体
4、一维动力模型首次提出了超低摩擦效应的存在并用实验以验证 等通过对地下断层岩性的分析发现地质断层间接触面软弱夹层对超低摩擦效应有重要影响 吴昊等在 等建立的块体模型基础上揭示了超低摩擦现象的出现是块系岩体间法向冲击力重新分布和动摩擦系数不断变化的结果 王明洋等引入冲击能量因子揭示了深部岩体超低摩擦现象、准共振及其动力事件相互关联的本质 何满潮等利用测量技术研究块系花岗岩的超低摩擦现象研究表明持续的周期扰动使层理裂隙间距累积变宽动力扰动荷载使花岗岩块体之间摩擦力降低 李利萍等通过自主研发的超低摩擦实验装置以工作块体水平位移作为反映超低摩擦效应的重要参数发现工作块体动能具有集聚效应、短时峰值特性和间
5、歇性具有由高频向低频转移的频移现象对于倾斜煤岩层发生冲击地压的研究潘一山等通过建立深部岩体的黏滑失稳模型分析了在开采过程中由于断层活化而导致冲击地压间歇性的一般规律 刘思妤等认为冲击地压是完整硬岩以失稳响应表达的自适应过程其中动力干扰并不是应力接近强度时的简单激发作用它可在二次应力远低于强度的情况下诱发冲击地压 穆贵清等以大安山矿为背景认为煤岩冲击是围岩应力演化的结果受煤岩物理力学性质、地质构造以及开采技术等多重因素的影响 文畅平综合考虑冲击地压发生的自然因素、地质因素基于 判别法对比分析了煤层厚度和煤层倾角对冲击倾向性的敏感性强弱 宋义敏等采用岩石直剪摩擦法对直剪滑动模型进行断层冲击试验研究
6、发现高应力环境下较小的侧向压力可使模型发生滑移 王襄禹等对红岭煤矿近断层采动巷道进行模拟发现巷道断层侧出现明显应力集中与剪切滑移现象 韩科明等模拟上下盘开采发现断层面关键层相对于界面滑移方式取决于侧向挤压力的大小而是否引起断层煤柱活化则取决于正应力与切应力的大小 王正义等基于动静组合加载相似模拟实验发现急倾斜煤岩冲击显现以顶板侧煤体为主底板侧煤体次之属于整体失稳型冲击地压在上覆岩层压力与顶板断裂导致的垂直扰动与水平冲击作用下煤岩块体接触面间摩擦力减弱从而诱发超低摩擦型冲击地压 采用数值模拟的方法建立不同倾角大小的砂岩煤砂岩组合块体模型以水平位移、加速度幅值等表征块体超低摩擦效应特征参数研究煤岩
7、层倾角与超低摩擦效应关系旨在对超低摩擦型冲击地压监测防治提供一定的理论依据含倾角块体模拟.研究背景煤矿冲击地压一般发生在深部高地应力下给巷道设施、人员安全带来严重威胁 沈阳焦煤红阳三矿发自 然 灾 害 学 报第 卷生超低摩擦型冲击地压如图 所示在上覆岩层压力与顶板断裂导致的冲击扰动作用下巨量煤体沿分离界面滑移最终完全堵塞巷道 与通常的冲击地压现象不同巷道顶板、底板整体保持完好仅煤体向巷道侧方产生滑移导致此次冲击地压事故发生图 沈阳焦煤红阳三矿超低摩擦示意图.国内外学者通过理论推导和实验手段对超低摩擦现象的研究开展了大量工作而前人对超低摩擦效应的研究所考虑的工况均为水平煤岩层未考虑倾斜煤岩层超低
8、摩擦效应的影响机制 复杂力学环境使深部煤岩体结构趋于破碎然而这些煤岩块体并不全是近水平排列的 实际上在对我国冲击地压矿井生产条件统计分析中倾斜煤层占 以上 因此对于超低摩擦效应机理研究考虑倾角的影响显得非常必要也更符合实际现场工况 结合图 超低摩擦效应发生实际背景以图 所示倾斜煤岩块体为研究对象考虑顶板断裂或采掘活动引起的垂直扰动和水平冲击作用于倾斜块体模型分析深部地下煤岩开采发生剪切滑移超低摩擦效应的影响机理 图 倾斜煤岩层示意图 图 数值模型图.模型构建根据图 所示倾斜煤岩层理划分从上至下将基本顶、直接顶、煤层、直接底和基本底假设成 个块体利用 软件建立整体倾斜的煤岩组合块体摩尔库伦塑性模
9、型 块体模型由 个尺寸大小相等的块体堆叠而成中间工作块体为煤块其余为砂岩块体 块体尺寸均为 各块体网格划分为 块体数值模型如图 所示在施加冲击载荷前首先固定模型底端限制底面移动并在块体模型上方施加轴压方向垂直于顶部块体表面向下用以维持块体模型稳定 块体所受重力与中心线夹角即等效为整个煤岩块体模型倾角参考实际开采工况与文献模拟倾角大小取、并设置 个监测点分别位于砂岩、煤块、砂岩 的形心位置如图 所示分析受外部扰动冲击后煤与砂岩块体的超低摩擦效应参数变化特征.模型参数与加载方式考虑模型的加载方式有 种仅垂直扰动作用与垂直和水平冲击共同作用试验过程中通过改变块体模型重力加速度方向模拟不同倾角的块体模
10、型 通过设置块体模型倾角替换为改变后重力方向与初始重力方向夹角减小模拟运算量使之更易收敛第 期李利萍等:动载扰动下倾斜煤岩组合块体超低摩擦效应研究图 单轴压缩实验.根据国家标准选取直径 、高 的砂岩和煤柱体各 组通过现有实验装置依次对试件进行单轴压缩实验如图 所示 根据监测得到的应力应变曲线求出弹性模量 与泊松比 组实验参数取平均值 由式()、式()得到煤与砂岩块体的体积模量 与剪切模量 结合已有实验数据块体模型试验参数整理为表 表 为设置的块体接触面参数()()()()表 砂岩与煤块力学参数 岩性参数砂岩煤密度/(/)体积模量/.剪切模量/.抗拉强度/.黏聚力/.内摩擦角/表 块体接触面参数
11、 界面参数砂岩与砂岩砂岩与煤法向刚度/剪切刚度/抗拉强度/.黏聚力/摩擦角/()各块体间设置硬接触面即是在荷载作用下物块可以产生滑移和分离的真实接触面界面参数法向刚度和剪切刚度取界面最硬相邻等效刚度的 倍换算关系如式():()()式中:为块体体积模量()为块体剪切模量()为法向刚度为剪切刚度D为接触面连接网格的最小尺寸.冲击与扰动幅值以准静态弹性接触的 理论为依据忽略模型的振动将冲击能 转化为幅值、作用时间 的半正弦形式作用力有关冲击载荷()的函数形式见式()其与力学转化见图:()()()式中:为动力载荷幅值()为扰动频率()冲击能量强度跟冲击扰动幅值密切相关而冲击扰动频率反映了模型所受冲击能
12、量的集中程度 根据课题组过去对水平工况块体的研究表明垂直扰动幅值越大频率越低超低摩擦效应越显著 因此在以前实验数据基础上取扰动频率 为 垂直扰动的幅值设定五组其幅值大小分别为、自 然 灾 害 学 报第 卷图 冲击能量力学转化形式.可行性验证对建立数值模型进行可行性验证在倾角 情形下对煤岩组合块体施加垂直扰动与水平冲击得到工作块体煤块 水平位移振动曲线 通过与已有实验结果对比由图()、()发现由于加载的扰动量级不同且模拟与实验本身存在差异得出数值模拟与前人实验结果的数据并不完全相同但变化趋势是一致的 组试验的工作块体在受水平冲击后出现了不同程度的滑移水平位移的变化随着扰动加载出现准正弦的周期性曲
13、线并逐渐趋于稳定达到平衡 因此从试验结果与理论分析的角度验证了所建立数值模型的可行性图 试验结果与已有研究对比.试验结果分析对于临近工作面的煤岩块体若遇冲击载荷强扰动时块体间产生周期性的挤压与分离使得块体沿倾斜煤岩层接触面产生滑移易诱发超低摩擦效应 因此考虑上覆岩层压力与顶板断裂等导致的垂直扰动与水平冲击对倾斜煤岩块体影响研究受垂直扰动和水平冲击作用下倾斜煤岩组合块体超低摩擦特征参数变化情况对进一步揭示煤岩超低摩擦型冲击地压具有重要意义.垂直扰动影响分析工作块体在轴压与垂直扰动的共同作用下其运动状态发生改变在垂直扰动施加瞬间块体在水平方向上发生瞬时滑动 由于扰动周期性的上下振动使得块体间也随之
14、产生周期性的挤压与分离导致块体在振动过程中滑移 随时间的推移块体水平位移幅度逐渐减小并逐渐稳定在一个新的平衡位置 针对顶板断裂对煤岩顶部造成强动载这一工况通过改变块体模型倾角与扰动幅值大小研究仅垂直扰动作用下煤岩层倾斜角度对超低摩擦效应的影响得到垂直扰动作用下工作块体水平位移与加速度响应规律如下:)工作块体水平位移随垂直扰动幅值与煤岩块体模型倾角增大而增大根据不同幅值的垂直扰动冲击作用下工作块体煤块 水平位移随倾角变化关系可知如图 所示当倾角一定时加载垂直扰动幅值越大工作块体水平位移振幅也越大此时若给予较小量级的水平冲击扰动煤块极易产生超低摩擦效应 而在同一扰动幅值下随着块体模型倾斜角度增大工
15、作块体水平位移也随之增大工作块体水平位移随倾角增加可近似成线性关系 这表明煤岩块体发生超低摩擦滑动失稳的可能性将随着扰动幅值与倾角的增加而急剧上升第 期李利萍等:动载扰动下倾斜煤岩组合块体超低摩擦效应研究图 垂直扰动作用下工作块体水平位移变化曲线.)工作块体水平方向与垂直方向加速度峰值随倾角增加而增加倾角 时产生摩擦自锁现象与块体接触面摩擦角等参数有关对超低摩擦效应有抑制作用在强扰动作用下块系煤岩体接触面间发生剧烈变化块体加速度在受到冲击作用瞬间产生剧烈波动随后逐步衰减为正弦扰动波 考虑含倾角煤岩块体受垂直扰动作用时水平加速度响应情况 设定在同一幅值、不同倾角块体模型探究其加速度历时响应曲线黑
16、色实线表示水平方向加速度响应曲线红色虚线表示垂直方向加速度响应曲线 通过对比分析工作块体的加速度历时响应曲线如图 所示倾角为、时水平加速度最大峰值分别为.、.、.、.、.、./垂直加速度最大峰值分别为.、.、.、.、.、./结果表明除 倾角外工作块体水平方向与垂直方向加速度峰值均随着倾角增加而增加 在不同倾角加速度响应曲线中倾角为 时水平方向加速度最大峰值明显小于 与 出现一个极小值点而垂直方向加速度与此相反其最大峰值明显大于与 出现一个极大值点 结合块体模型接触面设定的摩擦角参数为 分析原因极有可能是煤块底部产生了摩擦自锁现象图 不同倾角下工作块体加速度历时响应曲线.)倾角 恒定时随垂直扰动
17、幅值增加加速度受迫振动时间延长工作块体趋于稳定达到平衡状态时间延长根据加速度响应情况可将加速度历时曲线分为块体受扰动作用冲击振荡阶段随后衰减达到平衡收敛阶段 如图 所示由于正弦波动不稳定将加速度峰值点 /确定为判定 个阶段分界点并观察冲击自 然 灾 害 学 报第 卷振荡阶段区间变化 可以发现在倾角一定时随着垂直扰动幅值增加冲击能量衰减所需时间增加表现为块体水平与垂直加速度冲击振荡阶段时间延长块体趋于稳定达到平衡状态时间延长 这一结果与文献结论“作用于块体的法向力分别经历受迫振动和自由振动阶段变化形式为类正弦波其幅值不断衰减”相符验证了模拟可靠性图 不同扰动幅值下工作块体加速度历时响应曲线.综上
18、分析当煤岩块体所受垂直扰动载荷越大时超低摩擦效应强度越大且持续时间越长 在分析煤岩块体倾角与超低摩擦效应关系时不能忽略摩擦自锁因素影响这一因素与接触面摩擦角有关对块体间超低摩擦效应有抑制作用 研究发现煤岩块体倾角越大位移错动越明显产生滑动失稳的可能性也越大超低摩擦效应越显著.垂直扰动与水平冲击共同作用影响分析工程实践表明回采扰动、机械振动和顶板断裂等扰动载荷是导致超低摩擦效应的重要影响因素 工作面采动应力与顶板断裂引起垂直与水平方向块系煤岩体发生反复振动从而诱发导超低摩擦型冲击地压因此模拟垂直扰动与水平冲击共同作用时从不同延迟冲击时刻、垂直扰动与水平冲击对比、不同监测点变化趋势等角度研究含倾角
19、煤岩组合块体的超低摩擦效应特征参数的响应情况分析得到规律如下:)工作块体位移曲线重新振荡剧烈程度随水平冲击延迟周期增加而增大垂直扰动与延迟水平冲击共同作用结果除了位移振荡数值上的相互叠加仍保持周期性、收敛性特征 此外 种扰动共同作用效果还受到水平冲击时间点的影响同一周期不同时刻作用效果不同在位移振荡曲线同一周期内波峰点处施加水平冲击比均值点与波谷点处施加水平冲击位移振荡程度更剧烈对于倾斜煤岩组合块体模型工作块体的水平位移跟垂直扰动与水平冲击密切相关 但由于深部地下煤岩层结构复杂垂直扰动与水平冲击往往并不同时作用在倾斜煤岩层 前人研究表明垂直扰动与水平冲击的间隔时间对超低摩擦效应有很大影响 在垂
20、直扰动与水平冲击共同作用时垂直扰动会改变煤岩块体界面压紧或疏松状态使得块体界面产生周期性的挤压或分离对水平冲击能否引发超低摩擦效应以及引发超低摩擦效应强度大小至关重要 因此考虑垂直扰动与延迟水平冲击共同作用设置块体模型倾角为保持垂直扰动幅值与水平冲击幅值恒定研究水平冲击延迟时间对工作块体水平位移的影响规律测得不同延迟时间的水平冲击对块体水平位移的变化情况如图 所示至 为不同位移振荡周期至 时刻为位移振荡曲线不同周期内的波谷点可以看出 延迟水平冲击作用于不同周期波谷点处对煤岩块体水平位移影响效果不同 工作块体受垂直扰动与水第 期李利萍等:动载扰动下倾斜煤岩组合块体超低摩擦效应研究图 不同延迟时间
21、冲击下工作块体水平位移的变化情况.平冲击共同作用下其侧向水平位移产生了正弦形式波动 无水平冲击时块体水平位移波动幅值逐渐减小并收敛当垂直扰动作用一段时间后再对块体施加水平冲击会使得原本应该减小的波动幅值又开始产生剧烈波动在作用D.约一个周期后又逐渐恢复平稳状态 对比垂直扰动单独作用(无水平冲击)与垂直水平冲击扰动共同作用两者收敛的趋势是一致的但不同的是延迟的水平冲击会使得本应该趋于平衡的曲线又重新振荡起来但振荡的幅度受延迟时间的影响 可以看出至 时刻工作块体位移曲线重新振荡的剧烈程度随水平冲击延迟周期的增加呈线性递增趋势为了更清晰地描述延迟时间对于超低摩擦效应的影响将各周期内的波峰、均值、波谷
22、点的位移差值进行整理得到关系如图 所示可以发现上述现象不仅体现在波谷点处同一周期内不同延迟时间冲击点水平冲击加载作用均不相同 根据图 位移曲线周期划分建立如图 所示新振荡曲线与原振荡曲线最大差值随振荡周期变化关系分别在位移振荡曲线的波谷、均值、波峰处施加幅值恒定的水平冲击 通过对比可以发现位移振荡剧烈程度并不是简单随延迟时间增加而增加例如延迟水平冲击作用在 周期波峰点比 周期波谷与均值点位移振荡程度更剧烈 不同时刻施加水平冲击作用效果有显著差异波谷、均值、波峰处的水平冲击对位移振荡效果依次增大这种差异体现在 周期尤为明显 这表明作用在不同时间点上对超低摩擦效应影响效果不同而在 周期波峰点的水平
23、冲击诱导超低摩擦效应的概率会比同周期的均值点和波谷点更高 同时可以发现工作块体新位移振荡曲线与原曲线最大差值随振荡周期近似成线性增加但三者斜率不同有交汇趋势 图 不同冲击点位移最大差值与振荡周期变化关系 图 不同冲击幅值对工作块体水平位移的影响.)水平冲击与垂直扰动均对块体超低摩擦效应有不同程度影响在保持煤岩块体模型稳定下随冲击扰动幅值增加垂直扰动比水平冲击对超低摩擦效应影响更持续是超低摩擦现象发生的主导因素通过 组情况对比在倾角 下一组保持垂直扰动幅值不变改变其水平冲击幅值另一组保持水平冲击幅值不变改变其垂直扰动幅值考察垂直与水平冲击分别对超低摩擦效应的作用效果 根据水平位移随垂直扰动与水平
24、冲击变化曲线(图)可知当垂直扰动幅值从 依次增加到、块体水平位移近似成线性增加当水平冲击幅值从 依次增加到、块体水平位移近似成抛物线增加 显然对比分析 组不同试验加载方案可以看出水平冲击与垂直扰动均对块体的超低摩擦效应有不同程度的影响曲线后半段表明随着幅值的增加水平冲击扰动对超低摩擦效应的贡献逐渐减弱而垂直扰动始终是自 然 灾 害 学 报第 卷超低摩擦现象发生的主导因素)含倾角煤岩块体模型中在垂直扰动与水平冲击共同作用下越靠近施加垂直扰动的块体超低摩擦效应越明显图 块体模型顶部、中部、底部监测点水平位移变化情况.对比煤岩组合块体模型砂岩(顶部)、煤块(中部)与砂岩(底部)个监测点关于倾角变化的
25、位移参数可知如图 所示随着倾角的增加 个监测点水平位移随之增大 底部监测点水平位移从.增加到.中部监测点从.增加到.顶部监测点水平位移从.增加到.可以发现 个监测点间水平位移的增长幅度随着倾角变化差距越来越大在倾角 时 顶 部 与 底 部 监 测 点 位 移 相 差 达 到 最 大 为.这表明在垂直扰动与水平冲击共同作用下倾斜煤岩层的超低摩擦效应从下至上出现了差异性砂岩块体 超低摩擦效应响应情况比煤块 和砂岩块体 更敏感且随倾角的增加这种差异有不断放大的趋势 块体模型中越靠近垂直扰动施加点垂直扰动作用效果是越大的 块体界面间由于垂直扰动的作用产生挤压或分离的效果随扰动的增强越发明显这可以从结论
26、)和)得到 因此越靠近垂直扰动施加点块体超低摩擦效应越显著主要表现为剪切力所引起的块体顶部超低摩擦效应结合工程实际结果表明对于倾斜煤岩层除了受到采空区上覆岩层压力影响还受到工作面煤层倾斜因素影响 倾斜煤岩层使得采空区上覆岩层重心偏向下侧造成上端头支撑压力高于下端头支撑压力而在重力作用下工作面侧后方煤岩体向下端头方向滑移 上端头始终承受采空区段悬臂顶板压力使得工作面上端头围岩应力始终高于下端头围岩应力这也是倾斜或急倾斜煤岩层巷道冲击地压发生在上端头频次多于下端头频次的原因结论研究含倾角煤岩块体模型超低摩擦响应规律时考虑冲击幅值、倾角大小、延迟时间、加载方式等因素对超低摩擦效应的影响分析不同倾角煤
27、岩组合块体受冲击扰动时超低摩擦特征参数位移与加速度的响应规律得出以下结论:)含倾角煤岩组合块体仅受垂直扰动作用时随扰动幅值的增加工作块体水平位移随之增大块体由振动到平衡收敛的时间延长超低摩擦效应越显著 垂直扰动对超低摩擦效应的影响更持续是发生超低摩擦现象的主导因素)垂直扰动与水平冲击共同作用时随块体模型倾角的增加工作块体水平位移随之增大块体模型顶部、中部和底部监测点水平位移的增加幅度依次减弱越靠近垂直扰动施加的块体超低摩擦效应越显著)延迟的水平冲击作用使得趋于平衡收敛的位移曲线又重新开始振荡随水平冲击延迟周期增加位移振荡程度越剧烈垂直扰动与水平冲击共同作用结果除了数值上的叠加规律还具有周期性与
28、收敛性特征 共同作用还受到延迟水平冲击作用时刻的影响同一位移振荡周期内波峰点水平冲击比波谷与均值点作用效果更显著更易诱发超低摩擦效应参考文献:.():.():.():.():.第 期李利萍等:动载扰动下倾斜煤岩组合块体超低摩擦效应研究 .:.吴昊 方秦 王洪亮.深部块系岩体超低摩擦现象的机理分析.岩土工程学报 ():.():.()王明洋 李杰 李凯锐.深部岩体非线性力学能量作用原理与应用.岩石力学与工程学报 ():.():.()何满潮 王炀 刘冬桥 等.基于二维数字图像相关技术的块系花岗岩超低摩擦效应实验研究.煤炭学报 ():.():.()李利萍 潘一山.深部煤岩超低摩擦效应能量特征试验研究.
29、煤炭学报 (增刊):.():.()潘一山 王来贵 章梦涛 等.断层冲击地压发生的理论与试验研究.岩石力学与工程学报 ():.():.()刘思妤 徐则民.基于动静应力耦合的深埋隧道岩爆灾害控制.自然灾害学报 ():.():.()穆贵清 刘晓斐 胡少斌 等.大安山矿工作面开采煤岩冲击危险性动态评价研究.自然灾害学报 ():.():.()文畅平.基于 判别分析法的冲击地压预测与危险性分级.自然灾害学报 ():.():.()宋义敏 马少鹏 杨小彬 等.断层冲击地压失稳瞬态过程的试验研究.岩石力学与工程学报 ():.():.()王襄禹 柏建彪 李磊 等.近断层采动巷道变形破坏机制与控制技术研究.采矿与安
30、全工程学报 ():.():.()韩科明 于秋鸽 张华兴 等.上下盘开采影响下断层滑移失稳力学机制.煤炭学报 ():.():.()王正义 窦林名 何江.动静组合加载下急倾斜特厚煤层开采煤岩冲击破坏特征研究.采矿与安全工程学报 ():.():.()鞠文君 卢志国 高富强 等.煤岩冲击倾向性研究进展及综合定量评价指标探讨.岩石力学与工程学报 ():.():.()李利萍 吴金鹏 鞠翔宇 等.围压与冲击扰动作用下组合煤岩超低摩擦效应分析.地质力学学报 ():.():.()李利萍 唐垒 潘一山 等.应力波扰动对砂岩块体超低摩擦效应影响分析.岩石力学与工程学报 (增刊):.():.()陈育民 徐鼎平./基础与工程实例.版.北京:中国水利水电出版社:./.:.().():.王洪亮 葛涛 王德荣 等.块系岩体动力特性理论与实验对比分析.岩石力学与工程学报 ():.():.()