Damage Evolution and Dynamic Response of Anchorage Caverns with a Crack under Top Explosion
-
摘要: 基于相似模型试验,利用数值分析方法研究了含裂隙锚固洞室在顶爆下的损伤演化和动态响应规律,并探讨了裂隙倾角和长度对锚固洞室损伤分布和拱顶位移的影响。随着爆炸应力波的传播,在裂隙表面形成张拉损伤,并在裂隙尖端扩展出翼生裂纹,随后在锚固洞室附近的拱顶自由面、锚固区边界处以及底板中间依次形成较大拉伸损伤区域。随着裂隙倾角从0°加大到90°,拱顶锚固区及锚固区边界最大损伤区域面积先减小后增加,裂隙倾角为45°时面积最小。当裂隙的长度相同时,拱顶峰值位移随裂隙倾角的增大先增加后减小;当裂隙的倾角相同时,裂隙长度越长对应力波的阻碍作用越明显,裂隙长度为30 cm、倾角为45°时洞室破坏最小。Abstract: Based on the similar model test, the damage evolution and dynamic response law of anchored caverns with a crack under top explosion was studied by numerical analysis, and the influence of crack angle and length on the damage distribution and crown displacement of the anchor caverns was discussed.With the propagation of the explosion stress waves, tensile damage is formed on the surface of the fracture, and a winged crack extends at the tip of the fracture.Subsequently, a large tensile damage zone was formed successively at the vault free surface near the anchoring cavern, at the boundary of the anchoring zone, and in the middle of the floor.As the crack angle increases from 0° to 90°, the area of the largest damage zone at the boundary of the anchorage zone and the anchorage zone first decreases and then increases, and the area at the crack inclination angle of 45° is the smallest.When the cracks have the same length, the peak displacement of the crown increases first and then decreases with the increase of the crack dip angle.When the fracture angle is the same, the longer the crack length is, the more obvious the obstruction of the stress wave is, and the damage of the cavern with the facture 30 cm in length and 45° in inclination angle under top explosion is minimal.
-
Key words:
- underground engineering /
- blast loads /
- anchorage caverns /
- opening crack /
- damage evolution /
- numerical analysis
-
夹芯结构由于具有较好的吸能缓冲性能以及较高的比强度和比刚度,被广泛应用于航空航天、交通运输、航海等领域。近年来随着对夹芯板研究的不断深入,提出了各种桁架夹芯和褶皱型夹芯等形式。其中,层级结构思想引起了众多学者的关注。Lakes[1]通过研究发现,层级结构与非层级结构相比具有更加显著的刚度和强度,同时也具有足够的稳定性及可靠性。Zhang等[2-3]通过梯形波纹板夹芯结构的静态压缩和低速冲击实验,研究了不同芯层材料、不同载荷下结构的能量吸收和变形过程,结果表明不同芯层材料和冲击能量下结构的能量吸收性能不同。Hou等[4]对不同层数的梯形夹芯板结构进行了实验研究,发现不同层数结构的能量吸收性能不同。Meza等[5]对多层级桁架结构的演变模式进行了分析,提出了几种不同结构的破坏形式和变形模式。Kooistra等[6]研究了层级波纹板结构的6种失效模式,并给出了结构在不同相对密度及角度下的破坏机理图。Velea等[7]提出了一种制备二级夹芯结构的方法,将一级波纹板结构芯层替换为多孔填充芯层,通过实验和理论分析发现,二级结构有着良好的刚度和强度。Wu等[8]将金字塔型夹芯结构的芯层替换为层级夹芯结构,理论推导出了结构的等效刚度,通过实验研究了芯层几何参数对结构变形模式的影响,并总结出结构的6种不同变形模式。Wu等[9]设计了二级金字塔形夹芯结构,进行了理论分析和数值模拟,结果表明二级芯层结构具有较好的强度和稳定性。Liu等[10]将传统格栅夹芯结构芯层板替换为蜂窝夹芯结构,通过实验研究和数值模拟研究发现,在压缩载荷下二级芯层结构有着较好的能量吸收性能和稳定的变形性能,在轻量化和能量吸收方面具有更好的工程应用前景。
多孔结构及蜂窝结构有着良好的力学性能和能量吸收性能[11]。Jing等[12]通过实验研究了3种不同多孔芯层夹芯板结构在脉冲载荷下夹芯板的变形模式,结果表明不同芯层夹芯结构的变形模式不同,芯层和面板尺寸不同的夹芯结构的力学性能不同。Chen等[13-14]根据层级蜂窝结构的设计理论设计了不同芯层形式的层级蜂窝结构,并由理论分析得到结构变形模式,利用有限元对结构变形模式进行了验证,与传统蜂窝结构相比,层级蜂窝结构具有更好的能量吸收性能。Yin[15]、Qiao[16]等利用有限元软件对层级蜂窝结构进行了数值模拟,研究结果表明层级蜂窝结构的力学性能远优于传统蜂窝。
已有的研究主要集中于传统波纹板夹芯结构和层级蜂窝结构,关于多层级夹芯结构的力学性能与能量吸收的研究较少。本研究基于波纹板夹芯结构,设计了3种二级波纹板夹芯结构,通过理论分析和数值模拟相结合的方法研究一级和二级夹芯结构的变形模式与能量吸收性能,理论预测结构的临界失效载荷;建立不同芯层层数的二级波纹板夹芯结构的有限元模型,研究不同芯层厚度对多层级波纹板夹芯结构变形模式和能量吸收性能的影响,并与传统的一级夹芯波纹板结构进行对比,分析其比吸能与结构效率。
1. 层级夹芯结构的芯层设计和理论分析
1.1 芯层结构设计
根据层级蜂窝结构的设计理论,将一级梯形波纹板夹芯结构的芯层替换为三角形多孔夹芯结构,共设计了3种芯层结构,如图1所示。图1(a)为一级波纹板夹芯结构的芯层,斜板长度为la,顶板长度为l,斜板角度为
ω ,芯层板厚度为ta;二级结构的芯层如图1(b)、图1(c)和图1(d)所示,其中图1(b)为二级单层结构,图1(c)和图1(d)分别为二级双层和三层结构,二级结构芯层由大支撑面板和小支撑构成,水平方向小支撑为横向小支撑。二级结构中不同层数结构的尺寸参数相同,斜向大支撑面板长度为la,厚度为ta,小支撑长度为lb,厚度为tb,小支撑与外板角度均为θ ,外板与水平方向夹角为ω 。结构垂直于纸面方向的厚度为30 mm。1.2 理论分析
小支撑厚度与大支撑面板厚度相同的结构模型如图2所示,在外部载荷Fx和Fy作用下考虑横向小支撑的变形,基于小变形假设,由受力分析可以得到结构变形和力之间的关系,通过几何分析可以得到合力与变形之间的关系
{Fx2=Fncosω+FtsinωFy2=Fnsinω−Ftcosω (1) {δx=δncosω+δtsinωδy=δnsinω−δtcosω (2) 式中:Fx和Fy分别为沿x轴和y轴的外部载荷,
δx 和δy 为沿x轴和y轴的变形,Fn和Ft分别为夹芯夹层在外部载荷下的轴向压缩分量和剪切分量,δn 和δt 分别为夹芯夹层的轴向变形和剪切变形。当只有Fx作用时,
δy =0,δx ≠0,此时有Fn=EAδnl=Fx2cosω (3) δn=FnlEA (4) 式中:EA为夹芯夹层截面的抗拉刚度。由结构几何关系可知
δx=δncosω (5) 当只有Fy作用时,
δy ≠0,δx =0,此时有Fn=EAδnl=Fy2sinω (6) 同样有
δy=δnsinω (7) 当结构芯层为二级结构时,如图2(d)所示,(1)式~(7)式仍然成立。根据文献[17]给出的不同尺寸比下结构的失效机制,当
lb/la>√εy/πsinθ (其中εy 为材料屈服应变)时层级结构出现大支撑面板塑性屈服的失效模式。本研究中结构模型的尺寸满足上述条件,因此失效模式为大支撑面板发生塑性屈服,即大支撑面板达到材料屈服应力时结构失效。由(1)式和(2)式可得单胞结构大支撑面板在y方向上发生塑性屈服时结构的合力Fy′=(2n+2)σybtasinω (8) 式中:n为二级结构层数,b为结构垂直于纸面方向的厚度,σy为屈服应力。
根据结构变形模式,考虑横向小支撑的变形,如图3所示,Fb为横向小支撑在y方向的力,
δb 为δn 在y方向的分量,由此可知δb=Fbl3b3D (9) 弯曲刚度
D=Et3b12(1−ν2) (10) 式中:E为弹性模量,
ν 为泊松比。由此可以得到Fb=2nσyt3bsinω4(1−ν2)l2b (11) 式中:tb和lb分别为二级结构芯层小支撑厚度和长度。二级结构单胞在y方向的合力为
Fy=Fy′+Fb (12) 2. 有限元模型
有限元模型主要包括3部分:刚性压板、刚性支撑底板和波纹板夹芯结构。波纹板夹芯结构由芯层和上、下面板构成,上、下面板设置为刚体,厚度为2 mm。图4为二级双层结构有限元模型。波纹板夹芯结构和下部底板采用壳单元,刚性压板采用实体单元。刚性压板速度V=1 m/s,下部刚性底板固定,中部波纹板夹芯结构芯层与上、下两面板之间采用绑定约束。结构材料为铝合金6061-T6,材料参数如表1所示,其中
ρ 为密度。图5为材料应力-应变曲线[6]。结构单元类型采用四节点线性缩减积分壳单元(S4R),单元特征长度L=1 mm。芯层和上、下面板为自接触,刚性压板与上面板、底板与下面板为面面接触,静摩擦系数和动摩擦系数分别为0.20和0.15[18]。表 1 6061-T6铝合金材料性能Table 1. Material properties of aluminum 6061-T6Material ρ/(g·cm–3) σy/MPa E/GPa ν Al 6061-T6 2.700 251 69 0.33 2.1 准静态分析
利用LS-DYNA进行准静态分析时[19]:首先在模拟过程中结构动能要相对于内能较小,其次结构的力-位移曲线相对于速度的变化是稳定的。如图6(a)所示,结构在压板速度为0.5 m/s和1.0 m/s时动能相对于内能较小,结构整体内能受速度变化的影响较小;图6(b)为结构在压板速度为0.5 m/s和1.0 m/s时的力-位移曲线,可以看出两条曲线的重合度较高,速度变化对结构力-位移曲线的影响较小。通过对结构进行网格敏感性验证,并考虑到计算时间等因素,计算时取压板的速度为1.0 m/s,网格单元尺寸为1 mm。
2.2 能量吸收指标
为了评估结构的能量吸收性能[20],选取总能量吸收(EA)、平均载荷(Pm)、比吸能(
ξEA )、载荷效率(AE)4个评价指标。总能量吸收EA为结构在变形过程中吸收的能量
EA=∫δ0F(x)dx (13) 式中:
δ 为结构的变形位移(取δ =70 mm),F(x)为压缩过程中的瞬时载荷。平均载荷Pm定义为单位变形长度的能量吸收
Pm=EA/δ (14) 比吸能
ξEA 定义为结构单位质量的能量吸收ξEA=EA/m (15) 式中:m为结构总质量。比吸能越大,结构的能量吸收性能越好。
载荷效率定义为平均载荷与峰值载荷的比值
AE=Pm/Pk (16) 式中:Pk为峰值载荷。较好的吸能结构应该具有较大的载荷效率,并且力-位移曲线波动较小。
3. 模拟结果与分析
3.1 芯层板厚度的影响
建立一级和二级结构模型,其中一级结构芯层壁厚为2 mm,二级结构模型中ta=tb,即二级结构内部小支撑厚度与外部大支撑面板厚度相同。层级结构芯层板厚度分别为0.565、0.75、1.00、1.25 mm,结构几何参数如表2所示。图7给出了一级和二级层级结构在不同芯层板厚度下的力-位移曲线。由图7可知,二级结构平台力明显高于一级结构。一级结构在靠近芯层上部产生塑性铰,随着变形的增大未产生新的塑性铰,力-位移曲线变化平稳。二级结构随着变形增大发生折叠塑性变形,力-位移曲线会随之发生波动;在芯层厚度为1.25 mm时二级单层结构由于变形模式发生变化,力-位移曲线发生波动而随之上升;随着二级结构层数的增加,结构平台力增大;随着二级结构厚度的增大,其平台力增大。
表 2 模型尺寸Table 2. Size parameters of the modelType l/mm la/mm lb/mm ta/mm tb/mm θ/(°) ω/(°) First order 50 100 2 60 Second order-1 50 100 10 0.565/0.75/1.00/1.25 0.565/0.75/1.00/1.25 60 60 Second order-2 50 100 10 0.565/0.75/1.00/1.25 0.565/0.75/1.00/1.25 60 60 Second order-3 50 100 10 0.565/0.75/1.00/1.25 0.565/0.75/1.00/1.25 60 60 图8所示为不同芯层厚度下结构失效载荷理论值与模拟值的对比。对于数值模拟结果,当结构的承载力达到峰值时认为结构失效,即取力-位移曲线上的峰值载荷为结构失效载荷,可以看出理论值与模拟值吻合较好。当结构厚度为0.565 mm时,由于其厚度较小,在压缩载荷下的变形模式受屈曲影响较大,故理论值较模拟值偏大;当结构厚度为1.25 mm时,小支撑端部变形受边界条件的影响较大,理论值相对于模拟值偏小。
3.2 变形模式
图9所示为芯层板厚度t=2 mm时一级结构在不同应变下的变形模式,应变为结构在受载荷方向上的变形位移与结构上、下两面板之间原始距离之比。由图9可知:一级结构首先发生屈曲,在靠近芯层斜支撑上部产生塑性铰,随着变形的增大,结构发生失稳,从而产生非对称变形。图10所示为二级结构在应变为0.3时不同芯层厚度下的变形模式,可以看出:当厚度较小时,结构的变形发生在整体斜支撑与整体水平支撑的交界处,整体斜支撑内部的小支撑发生逐层折叠的塑性变形;当芯层厚度为1.25 mm时,结构斜支撑整体产生塑性铰而发生非对称变形;厚度变化对二级双层和三层结构的变形模式影响较小。在芯层厚度较小时,二级结构主要为小支撑整体逐层发生塑性变形,从而有更多能量转换为非弹性能,提高了结构整体的能量吸收能力。
3.3 能量吸收性能分析
二级结构变形模式主要为小支撑逐层折叠的塑性变形,相对于一级结构,其在变形过程中小支撑不断发生塑性变形产生塑性铰,从而有更多的能量转换为非弹性能,提高了结构的能量吸收性能。图11为不同芯层厚度结构的比吸能柱状图。由图11可知:层级波纹板夹芯结构的比吸能高于传统波纹板夹芯结构;随着结构芯层厚度增大,结构的比吸能不断提高;由于二级单层结构变形模式的变化,芯层厚度较小时二级单层结构的比吸能高于二级双层及三层结构,厚度较大时二级单层结构的比吸能低于二级双层及三层结构;二级双层结构的比吸能略高于二级三层结构。
图12为结构载荷效率随芯层厚度的变化关系。随着芯层厚度的增大,结构的载荷效率提高,二级单层结构由于变形模式的变化,载荷效率发生波动。由于二级单层结构在芯层厚度为1.25 mm时结构变形模式为斜支撑整体发生失稳,相对于小支撑发生折叠变形时其能量吸收性能降低,从力-位移曲线可以看出,结构平台力减小导致结构的载荷效率下降。二级双层和三层结构的载荷效率不断增加,二级三层结构的载荷效率高于二级双层结构。
4. 结 论
设计了不同层数的层级波纹板夹芯结构,利用数值模拟方法研究了在压缩载荷下层级夹芯结构的变形规律与能量吸收性能,理论推导了结构临界失效载荷公式;分析了结构参数对其变形模式和能量吸收性能的影响,并与一级结构进行了对比分析,得到如下研究结果。
(1)理论分析得到的结构失效载荷与数值模拟结果吻合较好。
(2)在准静态压缩载荷作用下,一级结构首先发生屈曲,随着变形增大,结构发生失稳而产生非对称变形;二级单层结构在厚度为1.25 mm时由于芯层整体产生塑性铰而失稳;二级双层和二级三层结构在整体斜支撑与整体水平支撑的交界处,整体斜支撑内部的小支撑发生逐层折叠的塑性变形。
(3)二级波纹板夹芯结构的比吸能显著大于一级结构;二级结构芯层小支撑发生逐层折叠的塑性变形时,结构能量吸收性能较好;随着芯层厚度的增大,二级结构的比吸能和载荷效率增加,芯层厚度较小时二级单层结构的比吸能高于二级双层和三层结构,二级双层结构的比吸能略大于二级三层结构。
-
表 1 CDP模型参数
Table 1. Parameters of CDP model
Density/(kg·m-3) E/GPa μ Dilation angle/(°) Eccentricity σb0/σc0 K Viscosity paramenter 1800 2.03 0.16 25 0.1 1.16 0.66667 0 -
[1] 俞缙, 钱七虎, 赵晓豹.岩体结构面对应力波传播规律影响的研究进展[J].兵工学报, 2009, 30(增刊2):308-316. http://d.old.wanfangdata.com.cn/Conference/7151837YU J, QIAN Q H, ZHAO X B.Research progress on effects of structural planes of rock mass on press wave propagation law[J]. Acta Armamentarii, 2009, 30(Suppl 2):308-316. http://d.old.wanfangdata.com.cn/Conference/7151837 [2] 顾金才, 陈安敏, 徐景茂, 等.在爆炸荷载条件下锚固洞室破坏形态对比试验研究[J].岩石力学与工程学报, 2008, 27(7):1315-1320. doi: 10.3321/j.issn:1000-6915.2008.07.003GU J C, CHEN A M, XU J M, et al.Model test study of failure patterns of anchored tunnel[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(7):1315-1320. doi: 10.3321/j.issn:1000-6915.2008.07.003 [3] 杨自友.锚固洞室的抗爆性能研究[D].合肥: 中国科学技术大学, 2008. http://cdmd.cnki.com.cn/Article/CDMD-10358-2008091576.htmYANG Z Y.Study on the anti-detonating characteristics of underground opening reinforced by rockbolts[D]. Hefei: University of Science & Technology of China, 2008. http://cdmd.cnki.com.cn/Article/CDMD-10358-2008091576.htm [4] 马海春.锚固洞室抗爆能力现场实验研究[D].合肥: 中国科学技术大学, 2011. http://cdmd.cnki.com.cn/Article/CDMD-10358-1011124968.htmMA H C.Research on explosion resistinhg capacity site tests of underground cavern reinforced by grouted rockbolts[D]. Hefei: University of Science and Technology of China, 2011. http://cdmd.cnki.com.cn/Article/CDMD-10358-1011124968.htm [5] 周纪军.爆炸动载对近区锚杆支护结构影响的试验研究[D].北京: 中国矿业大学, 2011. http://cdmd.cnki.com.cn/Article/CDMD-11413-1011184747.htmZHOU J J.Effects of blast loading on rockbolt anchorage structure near the working face in the model test[D]. Beijing: China University of Mining & Technology, 2011. http://cdmd.cnki.com.cn/Article/CDMD-11413-1011184747.htm [6] 王光勇, 顾金才, 陈安敏, 等.锚固洞室在顶爆作用下破坏形式及破坏过程研究[J].岩土工程学报, 2015, 37(8):1381-1389. http://d.old.wanfangdata.com.cn/Periodical/ytgcxb201508006WANG G Y, GU J C, CHEN A M, et al.Failure modes and process of tunnels reinforced by rockbolts under top explosion[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(8):1381-1389. http://d.old.wanfangdata.com.cn/Periodical/ytgcxb201508006 [7] CHANG X, WANG G Y, TANG C A, et al.Dynamic behavior of cement-mortar cavern reinforced by bars[J]. Engineering Failure Analysis, 2015, 55:343-354. doi: 10.1016/j.engfailanal.2015.07.020 [8] 徐干成, 袁伟泽, 顾金才, 等.地下洞库围岩外加固抗炸弹爆炸性能研究[J].岩石力学与工程学报, 2015, 34(9):1767-1776. http://www.cqvip.com/QK/96026X/201509/89837688504849534857484855.htmlXU G C, YUAN W Z, GU J C, et al.Explosive resistivity of anchored cavern surface rock[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(9):1767-1776. http://www.cqvip.com/QK/96026X/201509/89837688504849534857484855.html [9] 徐景茂, 顾金才, 陈安敏, 等.锚杆长度和间距对洞室抗爆性能影响研究[J].岩土力学, 2012, 33(11):3489-3496. http://d.old.wanfangdata.com.cn/Periodical/ytlx201211043XU J M, GU J C, CHEN A M, et al.Study of anti-explosion ability of reinforced tunnels with different anchor lengths and spacings[J]. Rock and Soil Mechanics, 2012, 33(11):3489-3496. http://d.old.wanfangdata.com.cn/Periodical/ytlx201211043 [10] 王卫华, 李夕兵, 周子龙, 等.不同应力波在张开节理处的能量传递规律[J].中南大学学报(自然科学版), 2006, 37(2):376-380. doi: 10.3969/j.issn.1672-7207.2006.02.032WANG W H, LI X B, ZHOU Z L, et al.Energy-transmitted rule of various stress waves across open joint[J]. Journal of Central South University (Science and Technology), 2006, 37(2):376-380. doi: 10.3969/j.issn.1672-7207.2006.02.032 [11] 柴少波, 李建春, 李海波.柱面波在节理岩体中的传播特性[J].岩石力学与工程学报, 2014, 33(3):523-530. http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb201403010CHAI S B, LI J C, LI H B.Propagation characteristics of cylindrical wave in joint rock masses[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(3):523-530. http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb201403010 [12] 郭东明, 闫鹏洋, 杨仁树, 等.爆破开挖中巷道围岩缺陷扩展的动焦散模型试验研究[J].采矿与安全工程学报, 2015, 32(5):728-734. http://d.old.wanfangdata.com.cn/Periodical/ksylydbgl201505005GUO D M, YAN P Y, YANG R S, et al.Dynamic caustics model experimental study on the defects extension of roadway surrounding rock when blasting excavation[J]. Journal of Mining & Safety Engineering, 2015, 32(5):728-734. http://d.old.wanfangdata.com.cn/Periodical/ksylydbgl201505005 [13] 郭东明, 刘康, 杨仁树, 等.爆炸冲击荷载下相邻巷道裂隙扩展机理模拟试验[J].振动与冲击, 2016, 35(2):178-183. http://d.old.wanfangdata.com.cn/Periodical/zdycj201602030GUO D M, LIU K, YANG R S, et al.Simulation experiments on the mechanism of adjacent tunnel crack extension induced by explosion load[J]. Journal of Vibration and Shock, 2016, 35(2):178-183. http://d.old.wanfangdata.com.cn/Periodical/zdycj201602030 [14] 刘红岩, 邓正定, 王新生, 等.节理岩体动态破坏的SHPB相似材料试验研究[J].岩土力学, 2014, 35(3):659-665. http://d.old.wanfangdata.com.cn/Periodical/ytlx201403009LIU H Y, DENG Z D, WANG X S, et al.Similar material test study of dynamic failure of jointed rock mass with SHPB[J]. Rock and Soil Mechanics, 2014, 35(3):659-665. http://d.old.wanfangdata.com.cn/Periodical/ytlx201403009 [15] YANG R S, WANG Y B, GUO D M, et al.Experimental research of crack propagation in polymethyl methacrylate material containing flaws under explosive stress waves[J]. Journal of Testing & Evaluation, 2016, 44(1):20140055. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=31c6583276e5680de2a1c0125b61e37c [16] 李超, 刘红岩, 阎锡东.动载下节理岩体破坏过程的数值试验研究[J].岩土力学, 2015, 36(增刊2):655-664. http://d.old.wanfangdata.com.cn/Conference/8776474LI C, LIU H Y, YAN X D.Numerical experiment of failure process of jointed rock mass under dynamic loading[J]. Rock and Soil Mechanics, 2015, 36(Suppl 2):655-664. http://d.old.wanfangdata.com.cn/Conference/8776474 [17] 杨仁树, 许鹏, 杨立云, 等.节理岩体中爆炸应力波传播规律的研究[J].金属矿山, 2016, 45(6):49-54. doi: 10.3969/j.issn.1001-1250.2016.06.009YANG R S, XU P, YANG L Y, et al.Study of regularity of explosive stress wave propagation in jointed rock mass[J]. Metal Mine, 2016, 45(6):49-54. doi: 10.3969/j.issn.1001-1250.2016.06.009 [18] 张秀丽, 焦玉勇, 刘泉声, 等.节理对爆炸波传播影响的数值研究[J].岩土力学, 2008, 29(3):717-721. doi: 10.3969/j.issn.1000-7598.2008.03.027ZHANG X L, JIAO Y Y, LIU Q S, et al.Numerical study on effect of joints on blasting wave propagation in rock mass[J]. Rock and Soil Mechanics, 2008, 29(3):717-721. doi: 10.3969/j.issn.1000-7598.2008.03.027 [19] DENG X F, ZHU J B, CHEN S G, et al.Numerical study on tunnel damage subject to blast-induced shock wave in jointed rock masses[J]. Tunnelling and Underground Space Technology, 2014, 43(6):88-100. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=c58ef1d2ad1e56ac09cf1d80e8ed3d35 [20] 王光勇, 顾金才, 陈安敏, 等.端部消波和加密锚杆支护洞室抗爆能力模型试验研究[J].岩石力学与工程学报, 2010, 29(1):51-58. http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb201001006WANG G Y, GU J C, CHEN A M, et al.Model test research on anti-explosion capacity of underground openings with end weakened by holes and anchor top reinforced by dense rock bolts[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(1):51-58. http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb201001006 [21] PRICE R H, BOYD P J, NOEL J S, et al.Relation between static and dynamic rock properties in welded and nonwelded tuff: SAND 94-0306C[R]. Albuquerque, NM: Sandia National Laboratories, 1994. [22] 王海兵, 张海波, 田宙, 等.岩石动力学计算中的网格效应及机理研究[J].兵工学报, 2016, 37(10):1828-1836. doi: 10.3969/j.issn.1000-1093.2016.10.009WANG H B, ZHANG H B, TIAN Z, et al.Mesh size effect and its mechanism research in numerical calculation of rock dynamics[J]. Acta Armamentarii, 2016, 37(10):1828-1836. doi: 10.3969/j.issn.1000-1093.2016.10.009 -