矿震冲击波作用下巷道围岩长期稳定性数值模拟方法

1.本发明涉及煤炭开采工程技术和数值模拟技术领域，特别涉及一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法。


背景技术：

2.冲击地压作为深部煤炭开采作业中常见的动力破坏现象，其发生往往伴随强矿震的产生。强矿震产生时，震动波形式的动载经过巷道上覆岩层，向巷道围岩空间传播变性能，到达巷道围岩时震动波会产生一定的能量衰减。在衰减后的震动波和井下采掘空间的静载应力耦合作用下，其中，静载应力主要为地压和支承压力组成，巷道围岩的材料性质不断劣化，耦合作用造成的时效损伤不断累积，巷道围岩最终达到变形临界值而发生失稳破坏。
3.矿震冲击波作用下巷道围岩长期稳定性分析包括两个方面：
①
震动波对巷道围岩的冲击作用，其作用时间短，时间在毫秒级；
②
巷道围岩在静载应力下的长期稳定性，作用时间长。其中，在冲击载荷和蠕变应力两个不同应变率的载荷下，巷道围岩的长期稳定性分析和数值模拟非常复杂，需要考虑岩石动力学过程、岩石静力学及流变学过程相互耦合作用引起巷道围岩发生失稳破裂的复杂力学效应。因此，建立巷道围岩长期稳定性分析和寿命评估方法，并应用于深部煤炭开采作业中，对深部煤炭资源的安全和高效开采意义重大。
4.目前，在深部煤炭地下开采过程围岩稳定性评估及数值模拟方法中，大多仅考虑了冲击载荷对巷道围岩稳定性的影响。但深部煤炭开采作业中，巷道围岩长期处于“高静载应力”和“冲击载荷”的共同作用下，其非稳态、非线性变形、损伤和能量演化规律十分复杂。因此，仅考虑冲击载荷对巷道围岩稳定性的影响，难以满足巷道围岩稳定性数值模拟求解的实际需求。为了能够综合考虑冲击载荷和岩石蠕变对巷道围岩稳定性的影响，以及实现在冲击载荷和蠕变应力两个不同应变率下的巷道围岩损伤参数的表征，本发明提出了一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法。


技术实现要素：

5.本发明的目的在于提供一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法，以解决上述背景技术中提出的问题。
6.本发明提供的技术方案如下：
7.一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法，包括以下步骤：
8.s1，建立煤矿采掘区域的巷道围岩几何模型，网格化所述巷道围岩几何模型为若干个岩石单元；
9.s2，根据巷道围岩内发生的实际矿震事件参数，赋值所述岩石单元波动载荷；
10.s3，赋予所述巷道围岩几何模型边界条件，所述边界条件包括位移边界条件或载荷边界条件；
11.s4，赋予所述岩石单元静载应力；
12.s5，建立所述岩石单元在矿震波冲击作用下的蠕变-冲击耦合模型方程组；
13.s6，赋予所述巷道围岩几何模型材料性质；
14.s7，判断所述岩石单元是否发生损伤破坏，所述损伤破坏包括拉伸损伤和/或剪切损伤；
15.s8，根据步骤s5中建立的所述蠕变-冲击耦合模型方程组，通过下述步骤对所述蠕变-冲击耦合模型方程组的长期稳定性进行模拟：
16.a.先对所述蠕变-冲击耦合模型方程组进行稳态求解，使所述蠕变-冲击耦合模型方程组达到应力平衡；
17.b.再对所述蠕变-冲击耦合模型方程组进行瞬态求解，得出所述蠕变-冲击耦合模型方程组的应力场及损伤场随时间的变化规律；
18.c.若达到蠕变求解时间，所述蠕变-冲击耦合模型方程组的求解结束；若未达到蠕变求解时间，则进行应力传递和求解时间转化，求解冲击载荷作用下的所述岩石单元的应力场和损伤场，重复步骤a-步骤b至所述蠕变-冲击耦合模型方程组达到蠕变求解时间；
19.d.判定所述岩石单元是否发生损伤破坏,并计算所述岩石单元的损伤值大小；根据所述岩石单元损伤值的大小，折减所述岩石单元的拉伸损伤参数为极小值，所述岩石单元的剪切损伤保留残余强度；
20.e.重复步骤a-步骤d并进行骤迭代计算，计算完成后，判断所述巷道围岩几何模型是否破坏或是否达到计算时间；若未达到计算时间则继续计算，直至计算到设定的时步数或者所述巷道围岩几何模型完全破坏。
21.进一步的，步骤s2中，所述矿震事件参数包括矿震能量、矿震频率和矿震位置。
22.进一步的，步骤s5中，所述蠕变-冲击耦合模型方程组包括岩石单元的蠕变应变εc、岩石单元的动态屈服强度σ
dyn
随应变率的变化规律。
23.进一步的，所述岩石单元的蠕变应变εc方程为：式中，a为材料相关系数，m为应力指数，n为时间硬化指数，e为蠕变激活能，t为蠕变时间。
24.进一步的，所述岩石单元的动态屈服强度σ
dyn
随应变率的变化规律方程为：式中，c为材料常数，为应变率，σ
sta
为材料准静态屈服强度。
25.进一步的，步骤s7中，所述岩石单元发生拉伸损伤通过最大拉应力准则判断；所述岩石单元发生剪切损伤通过摩尔库伦准则判断。
26.进一步的，所述岩石单元发生拉伸损伤时，其应力状态函数f1满足：f1＝-σ
3-σ
t0
＝0；
27.所述岩石单元发生剪切损伤时，其应力状态函数f2满足：式中，σ1、σ3为第一、第三主应力，σ
t0
为单元单轴抗拉强度，为内摩擦角，σ
c0
为单元单轴抗压强度σ3。
28.进一步的，所述岩石单元发生损伤时，其损伤变量d按以下关系给出：
[0029][0030]
式中，ε
t0
和ε
c0
分别为当岩石单元发生拉伸损伤和剪切损伤时对应的最大拉伸主应变和最大压缩主应变。
[0031]
进一步的，所述模拟方法还包括以下步骤：建立所述岩石单元蠕变寿命预测方程，所述岩石单元的蠕变寿命预测方程为logwc＝logβ'-plogt,式中，wc为单位体积外载荷，即冲击和蠕变载荷，所作总的机械功；b'和p为在一定条件下材料的常数。
[0032]
本发明的技术效果：
[0033]
1、本发明描述的数值模拟方法，根据深部地下煤炭开采的工程实际情况，考虑地下煤炭开采过程中矿震冲击及蠕变的影响，建立蠕变-冲击耦合作用下岩石时效失稳模型。通过在工程现场获得基本的物理力学参数，以及矿震导致的巷道围岩冲击载荷和地应力水平，进而模拟深部巷道围岩在蠕变-冲击耦合作用下的应力场、变形场及损伤场的时空状态演化。因此，本发明提供的数值模拟方法，能真实的模拟矿震冲击波导致巷道围岩损伤的破坏过程及破坏范围，对深部地下巷道失稳破坏过程的模拟更加严谨和合理。
[0034]
2、本发明描述的数值模拟方法以及巷道围岩寿命预测模型，可以更真实地模拟地下巷道围岩在矿震冲击波和蠕变作用下的失稳破坏位置、范围及寿命。其中，根据弹性波理论，计算矿震波引起的纵波和横向动应力计算公式，得到震动波最终传递至工作面的能量和波动载荷；根据动量守恒定律，熵守恒定律以及热力学第一定律，得到恒温条件下巷道围岩的蠕变寿命预测的能量模型，因此，本发明描述的数值模拟方法以及巷道围岩寿命预测模型理论严密且更切合工程实际情况。
附图说明
[0035]
图1为不同均匀性系数岩石介质细观基元强度的weibull分布图；
[0036]
图2为单轴拉伸和压缩下岩石单元损伤本构模型图；
[0037]
图3为本发明所述方法的计算流程示意图；
[0038]
图4为岩石单元蠕变和冲击载荷加载示意图；
[0039]
图5为煤柱在冲击和蠕变作用下损伤示意图；
[0040]
图6为矿震波冲击作用下煤柱能量演化示意图。
具体实施方式
[0041]
下面结合附图对本发明的优选实施例进行详细阐述。
[0042]
结合附图1至附图6所示，本发明提供一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法，包括以下步骤。
[0043]
s1，对煤矿采掘区域进行地质参数探测，建立巷道围岩几何模型，网格化所述巷道
围岩几何模型为若干个岩石单元。
[0044]
s2，采用微震监测系统，在煤矿巷道处实地采集巷道围岩内部发生的矿震数据信息，所述矿震数据信息包括矿震能量ej、矿震频率和矿震位置。
[0045]
本步骤中所使用的微震监测系统，是矿山开采作业中普遍使用的微震监测系统，可在加拿大esg微震监测系统、南非iss综合微震系统、波兰sos微震监测系统中选用。本实施例使用微震监测系统的普遍功能，即对微震事件的能量、频率、震源半径、震级和空间位置进行识别和计算的功能，因此本实施例中不再赘述微震监测系统的作业原理。
[0046]
s2-1，计算巷道围岩内部发生的矿震的震级m
l
。
[0047]
震级m
l
通过矿震能量ej和下述波兰经验公式计算
[0048]
log ej＝1.8+1.9m
l
[0049]
对上述公式进行举例说明：
[0050]
104j能量转换为震级m
l
＝1.15；
[0051]
106j能量转换为震级m
l
＝2.2；
[0052]
108j能量转换为震级m
l
＝3.2。
[0053]
s2-2，计算巷道围岩的质点移动峰值速度。
[0054]
微震事件发生后，矿震能量以震动波的形式释放，震动波在巷道围岩传递的过程中会产生能量的耗散，剩余能量的数值即为巷道围岩的冲击动能。通过下述公式计算巷道围岩的质点移动峰值速度：
[0055]
log rv＝log 2.6+0.5log m
l
[0056]
式中，
[0057]
r为冲击点与震源点的距离，可根据地质条件和生产经验给出r值；
[0058]ml
为矿震的震级。
[0059]
s2-3，赋值巷道围岩几何模型岩石单元处的波动载荷。
[0060]
根据弹性波理论，通过下述公式计算矿震冲击波在岩石单元处的波动载荷
[0061]
σd＝ρc
psv[0062]
式中，
[0063]
ρ为巷道围岩介质密度；
[0064]cps
为矿震波波速，通过现场测定得到；
[0065]
v为巷道围岩的质点移动峰值速度。
[0066]
本步骤中，矿震频率决定了波动载荷的作用时间。
[0067]
s2-4，赋予巷道围岩几何模型边界条件，所述边界条件包括位移边界条件或载荷边界条件。其中，
[0068]
位移边界条件由质点峰值震动速度乘以求解时间计算得到；
[0069]
载荷边界条件由波动载荷计算公式得到。
[0070]
s3，赋予岩石单元静载应力。
[0071]
使用直接法和间接法，测量岩石单元在x、y、z三个轴向上的地应力水平，进而得出岩石单元的静载应力水平。
[0072]
其中，
[0073]
直接法包括但不限于水压致裂法、扁千斤顶法、刚性包体应力计法和声发射法；
[0074]
间接法包括但不限于套孔应立解除法、局部应力解除法、松弛应变测量法、地球物理探测法。
[0075]
s4，基于流变学、弹性力学、岩石力学、岩石动力学基础理论，建立岩石单元在矿震波冲击作用下的蠕变-冲击耦合模型方程组。
[0076]
所述蠕变-冲击耦合模型方程组包括：“岩石单元的蠕变应变ε
c”和“岩石单元的动态屈服强度σ
dyn
随应变率的变化规律”。
[0077]
s4-1，通过下述norton-bailey蠕变定律公式，计算岩石单元的蠕变应变εc：
[0078][0079]
s4-2，巷道围岩的应变包括弹性应变和蠕变应变。弹性应变中，应力与应变之间具有一一对应的路径依赖性，而蠕变应变不具备此规律，因此，需采用率律理论进行描述，通过下述公式计算岩石单元的蠕变应变率：
[0080][0081]
由于巷道围岩在x、y、z三个轴向上均受到应力作用，即巷道围岩处于多轴应力状态，因此，对于多轴应力状态，通过下述公式表述岩石单元的蠕变应变率：
[0082][0083]
式中，
[0084]
a为材料相关系数，m为应力指数，n为时间硬化指数，e为蠕变激活能，a、m、n是依赖于温度和应力的系数，通过对实验数据拟合得到；
[0085]
t为蠕变时间，s
ij
为偏应力张量，σe为有效应力张量。
[0086]
综上所述，通过下述公式表述岩石单元的动态屈服强度σ
dyn
随应变率的变化规律：
[0087][0088]
式中，
[0089]
c为材料常数，通过实验数据拟合得到；
[0090]
为应变率；
[0091]
σ
sta
为材料准静态屈服强度，通过实验测得。
[0092]
s5，赋予巷道围岩几何模型材料性质。
[0093]
如附图1所示，由于巷道围岩的煤岩体是一种天然非均质材料，本实施例中采用weibull分布函数来描述材料的非均质性，其公式如下：
[0094][0095]
式中：
[0096]
u为岩石介质基元体力学性质参数(强度、弹性模量等)；
[0097]
u0为基元体力学性质的平均值；
[0098]
φ为均匀性系数，反映了岩石介质的均质性；
[0099]
是岩石基元体力学性质u的统计分布密度。
[0100]
s6-1，判断岩石单元是否发生损伤破坏。
[0101]
岩石单元发生损伤破坏的方式为拉伸损伤和剪切损伤，其判断依据为：
[0102]
①
根据最大拉应力准则，岩石单元发生拉伸损伤时，岩石单元的应力状态函数f1满足：
[0103]
f1＝-σ
3-σ
t0
＝0
[0104]
②
根据摩尔库伦准则，岩石单元发生剪切损伤时，岩石单元的应力状态函数f2满足：
[0105][0106]
式中，
[0107]
σ1、σ3为第一、第三主应力；
[0108]
σ
t0
为单元单轴抗拉强度；
[0109]
为内摩擦角；
[0110]
σ
c0
为单元单轴抗压强度。
[0111]
s6-2，判定岩石单元发生损伤时，岩石单元的损伤变量d按以下关系给出：
[0112][0113]
其中，
[0114]
ε
t0
和ε
c0
分别为当单元发生拉伸损伤和剪切损伤时对应的最大拉伸主应变和最大压缩主应变。
[0115]
附图2示出了单轴拉伸和压缩下岩石单元损伤本构模型。
[0116]
s7，根据步骤s4中建立的蠕变-冲击耦合模型方程组，通过下述步骤对巷道围岩几何模型的长期稳定性进行数值模拟计算：
[0117]
a.先对所述蠕变-冲击耦合模型方程组进行稳态求解，使所述蠕变-冲击耦合模型方程组达到应力平衡；
[0118]
b.再对所述蠕变-冲击耦合模型方程组进行瞬态求解，得出所述蠕变-冲击耦合模型方程组的应力场及损伤场随时间的变化规律；
[0119]
c.若达到蠕变求解时间，所述蠕变-冲击耦合模型方程组的求解结束；若未达到蠕变求解时间，则进行应力传递和求解时间转化，求解冲击载荷作用下的所述岩石单元的应力场和损伤场，重复步骤a-步骤b至达到蠕变求解时间；
[0120]
d.判定所述岩石单元是否发生损伤破坏,并根据弹脆性的损伤本构关系，计算所述岩石单元的损伤值大小；根据所述岩石单元损伤值的大小，折减拉伸损伤岩石单元的参数为极小值，所述剪切损伤岩石单元保留残余强度；
[0121]
e.重复步骤a-步骤d进行骤迭代计算，计算完成后，判断所述巷道围岩几何模型是
否破坏或是否达到计算时间；若未达到计算时间则继续计算，直至计算到设定的时步数或者所述巷道围岩几何模型完全破坏。
[0122]
s8，建立基于能量理论的蠕变和冲击作用下岩石单元寿命预测方程。
[0123]
根据不可逆热力学原理，遵循动量守恒定律，熵守恒定律以及热力学第一定律，依据
[0124]
①
动量守恒定律公式
[0125][0126]
②
熵守恒定律公式
[0127][0128]
③
热力学第一定律公式
[0129][0130]
其中，
[0131]
①
动量守恒定律，即euler第一运动定律，该定律认为连续物体内任一部分物质的动量变化率等于作用于该部分物质上所有外力的主向量，式中σ
ij,j
为cauchy应力张量分量，ρ为物体密度，fi为体积力矢量分量，为连续体质点速率的牛顿时间导数；
[0132]
②
熵守恒定律，该定律表明了熵产生的来源，认为连续物体的熵增率与热通量的和等于其熵产生与热供给的和，式中σ
*
为单位质量内部熵产生，为单位质量熵的牛顿时间导数，γ为单位质量热供给，h为热通量，t为绝对温度(k)；
[0133]
③
热力学第一定律，即能量守恒定律，作用于系统上的功的增量δw加上系统接受的热的增量δq等于系统内能的增量de加上动能的增量，当不考虑体力矩和面力矩时，有
[0134][0135][0136][0137][0138]
则该公式转化为：
[0139][0140]
式中，为单位质量含有内能的牛顿时间导数，为应变分量的牛顿时间导数。
[0141]
根据上述公式，将岩石单元作为处于一定温度环境中的连续致密材料，其满足以上守恒定律，在承受蠕变载荷情况下，机械功w
t
＝σε，根据步骤s4-1中的norton-bailey蠕变定律公式，在等温条件下：
[0142]
ε＝aσntm[0143]
式中，
[0144]
a，n，m为材料常数，t为时间。
[0145]
根据步骤s4-2中的蠕变应变率计算公式，分两种情况，一种是恒应力情况，一种是变应力情况，考虑扰动应力作用，此处仅考虑变应力情况，即应力是时间的函数σ＝f(t)。
[0146]
通过公式推导可得在恒温条件下，岩石单元的蠕变寿命预测能量模型：
[0147]
log wc＝log b'-p log t
[0148]
式中，
[0149]
wc为单位体积外载荷(冲击和蠕变载荷)所作总的机械功，b'和p为在一定条件下材料的常数。
[0150]
本实施例给出的步骤，在现有的仿真模拟软件，例如ansys、matdem以及comsol等软件中均可以进行数值模拟。本实施例通过给出所描述的数值模拟方法步骤，以及通过附图3给出的程序逻辑步骤，使用上述仿真模拟软件进行设定程序的计算以及生成图文或展示数据，并未超出上述仿真模拟软件的基础功能。
[0151]
在工程实践中，例如对标准试件通过数值模拟方法进行后处理分析，如附图4和附图5所示，能够得出深部煤炭开采过程中煤柱在蠕变-冲击作用下的变形和损伤场分布规律；如附图6所示，通过分析蠕变和冲击载荷作用下损伤区域的位置和范围，能够明确该区域进行煤炭地下开采过程中，巷道岩石在矿震冲击波作用下的损伤破坏情况。通过标准试件数值模拟模型结果分析和工程实例数值模拟结果分析，可以明确煤柱在受矿震冲击波和蠕变作用下的破坏范围，确定了损伤破坏位置和区域。
[0152]
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何不经过创造性劳动想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书所限定的保护范围为准。

技术特征：
1.一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法，其特征在于，包括以下步骤：s1，建立煤矿采掘区域的巷道围岩几何模型，网格化所述巷道围岩几何模型为若干个岩石单元；s2，根据巷道围岩内发生的实际矿震事件参数，赋值所述岩石单元波动载荷；s3，赋予所述巷道围岩几何模型边界条件，所述边界条件包括位移边界条件或载荷边界条件；s4，赋予所述岩石单元静载应力；s5，建立所述岩石单元在矿震波冲击作用下的蠕变-冲击耦合模型方程组；s6，赋予所述巷道围岩几何模型材料性质；s7，判断所述岩石单元是否发生损伤破坏，所述损伤破坏包括拉伸损伤和/或剪切损伤；s8，根据步骤s5中建立的所述蠕变-冲击耦合模型方程组，通过下述步骤对所述蠕变-冲击耦合模型方程组的长期稳定性进行模拟：a.先对所述蠕变-冲击耦合模型方程组进行稳态求解，使所述蠕变-冲击耦合模型方程组达到应力平衡；b.再对所述蠕变-冲击耦合模型方程组进行瞬态求解，得出所述蠕变-冲击耦合模型方程组的应力场及损伤场随时间的变化规律；c.若达到蠕变求解时间，所述蠕变-冲击耦合模型方程组的求解结束；若未达到蠕变求解时间，则进行应力传递和求解时间转化，求解冲击载荷作用下的所述岩石单元的应力场和损伤场，重复步骤a-步骤b至所述蠕变-冲击耦合模型方程组达到蠕变求解时间；d.判定所述岩石单元是否发生损伤破坏,并计算所述岩石单元的损伤值大小；根据所述岩石单元损伤值的大小，折减所述岩石单元的拉伸损伤参数为极小值，所述岩石单元的剪切损伤保留残余强度；e.重复步骤a-步骤d并进行骤迭代计算，计算完成后，判断所述巷道围岩几何模型是否破坏或是否达到计算时间；若未达到计算时间则继续计算，直至计算到设定的时步数或者所述巷道围岩几何模型完全破坏。2.根据权利要求1所述的一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法，其特征在于，步骤s2中，所述矿震事件参数包括矿震能量、矿震频率和矿震位置。3.根据权利要求1所述的一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法，其特征在于，步骤s5中，所述蠕变-冲击耦合模型方程组包括岩石单元的蠕变应变ε
c
、岩石单元的动态屈服强度σ
dyn
随应变率的变化规律。4.根据权利要求3所述的一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法，其特征在于，所述岩石单元的蠕变应变ε
c
方程为：式中，a为材料相关系数，m为应力指数，n为时间硬化指数，e为蠕变激活能，t为蠕变时间。5.根据权利要求3所述的一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法，
其特征在于，所述岩石单元的动态屈服强度σ
dyn
随应变率的变化规律方程为：式中，c为材料常数，为应变率，σ
sta
为材料准静态屈服强度。6.根据权利要求1所述的一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法，其特征在于，步骤s7中，所述岩石单元发生拉伸损伤通过最大拉应力准则判断；所述岩石单元发生剪切损伤通过摩尔库伦准则判断。7.根据权利要求6所述的一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法，其特征在于，所述岩石单元发生拉伸损伤时，其应力状态函数f1满足：f1＝-σ
3-σ
t0
＝0；所述岩石单元发生剪切损伤时，其应力状态函数f2满足：式中，σ1、σ3为第一、第三主应力，σ
t0
为单元单轴抗拉强度，为内摩擦角，σ
c0
为单元单轴抗压强度σ3。8.根据权利要求7所述的一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法，其特征在于，所述岩石单元发生损伤时，其损伤变量d按以下关系给出：式中，ε
t0
和ε
c0
分别为当岩石单元发生拉伸损伤和剪切损伤时对应的最大拉伸主应变和最大压缩主应变。9.根据权利要求1至8任一项权利要求所述的一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法，其特征在于，所述模拟方法还包括以下步骤：建立所述岩石单元蠕变寿命预测方程，所述岩石单元的蠕变寿命预测方程为：logw
c
＝logβ'-plogt,式中，w
c
为单位体积外载荷，即冲击和蠕变载荷，所作总的机械功；b'和p为在一定条件下材料的常数。

技术总结
本发明提供了一种矿震冲击波作用下巷道围岩长期稳定性数值模拟方法，通过建立巷道围岩几何模型，并将其网格化为岩石单元，根据实际测量的基本物理力学参数，赋予巷道围岩几何模型物理力学参数，并构建巷道围岩几何模型在矿震波冲击作用下的蠕变-冲击耦合模型方程组。通过对蠕变-冲击耦合模型方程组按照设定程序进行求解，判断巷道围岩几何模型是否破坏或是否达到计算时间来完成巷道围岩数值模拟过程。通过模拟深部巷道围岩在蠕变-冲击耦合作用下的应力场、变形场及损伤场的时空状态演化，本发明提供的数值模拟方法，能真实的模拟矿震冲击波导致巷道围岩损伤的破坏过程及破坏范围，对深部地下巷道失稳破坏过程的模拟更加严谨和合理。加严谨和合理。加严谨和合理。


技术研发人员：周广磊 孟祥军 陈连军 张广超 马俊鹏 王保齐 李志勇 周泽森
受保护的技术使用者：山东科技大学
技术研发日：2023.03.30
技术公布日：2023/10/6