冲击地压扰动响应失稳理论并行计算.pdfVIP

第50卷第1期煤炭学报Vol.50No.1

2025年1月JOURNALOFCHINACOALSOCIETYJan.2025

冲击地压扰动响应失稳理论并行计算

1233

潘一山，王学滨，郑一方，陈双印

(1.辽宁大学物理学院,辽宁沈阳110036；2.辽宁工程技术大学计算力学研究所,辽宁阜新123000；

3.辽宁工程技术大学力学与工程学院,辽宁阜新123000)

摘要：目前，冲击地压理论研究已经完成了从定性分析到定量分析的转变。巷道围岩临界应力计

算是巷道安全性评价的重要依据。鉴于冲击地压问题的极度复杂性，在理论上继续取得突破极为

困难。基于理论公式的巷道围岩临界应力计算，无法考虑更复杂的实际情况，例如非圆形巷道、

非静水压力和复杂岩层结构。冲击地压理论和数值计算相结合具有更加广阔的应用前景，能使冲

击地压理论进一步走向实际应用，这是极有价值的发展方向。这方面研究成果的成功取得依赖于

数值计算技术的快速发展。研究将当今较先进的岩层运动并行计算系统StrataKing(一种自主开发

的以拉格朗日元与离散元耦合方法为基础的非线性断裂力学GPU并行计算方法)与冲击地压扰动

响应失稳理论相结合，首次提出了圆形巷道扰动响应失稳理论的数值模拟方法。该方法的思想是

将非线性断裂力学数值分析方法中的Ⅱ型断裂能设定为中间变量，从而建立了静水压力条件下圆

形巷道围岩临界应力与冲击能指数之间的关系。为了获取冲击能指数的数值解，采用了仅出现一

个剪切面的理想岩样进行单轴压缩数值试验，以排除其他因素对应力−应变曲线峰后倾向于直线部

分斜率的影响。对于高角度剪切破裂，提出了将非标准岩样的计算结果转换成标准岩样的结果的

折算方法。折算后冲击能指数的范围为0.17～13.52，位于全国131个冲击地压矿井的调研数据之

内。巷道围岩临界应力的计算结果是理论结果的0.4～2.5倍，这与针对全国20个冲击地压矿井的

调研数据(临界应力的修正系数普遍大于1，甚至接近8)定性相符，从局部化破坏围岩比均匀破坏

围岩的承载力高的角度进行了解释。冲击地压与局部化的关系过去有讨论，扰动响应失稳理论与

局部化过去并无关系。通过局部化，扰动响应失稳理论与冲击地压之间在破坏机理上产生了密切

的关联。StrataKing可为冲击地压矿井巷道安全性评价提供强大的算力支撑。

关键词：冲击地压；定量分析；扰动响应失稳理论；冲击能指数；局部化；并行计算；临界应力

中图分类号：TD324文献标志码：A文章编号：0253−9993(2025)01−0081−11

Parallelcomputingofdisturbanceresponseinstabilitytheoryofrockbursts

1233

PANYishan,WANGXuebin,ZHENGYifang,CHENShuangyin

(1.SchoolofPhysics,LiaoningUniversity,Shenyang110036,