DOI： 10.11817/j.issn.1672-7207.2019.08.020

不同方向冲击加载作用下含裂纹巷道的破坏行为

周磊^{1, 2},朱哲明^{1, 2},范勇^{1, 2},刘瑞峰^{1, 2},邓帅^{1, 2},董玉清^{1, 2}

（1. 四川大学 深地科学与工程教育部重点实验室，四川 成都， 610065；

2. 四川大学 建筑与环境学院，四川 成都， 610065）

摘要:为了研究来自不同方向的冲击加载作用对含裂纹巷道围岩破坏方式的影响，采用青砂岩作为试验材料，进行相应的模型试验研究。首先，在巷道围岩的圆拱顶处设置单裂纹以模拟围岩内天然的裂纹缺陷；随后，采用冲击试验机进行动态加载试验，通过应变计确定巷道模型内预制裂纹起裂时间及扩展速度，并进行相应的数值模拟分析，得到裂纹的扩展路径规律与巷道围岩的破坏方式；最后，将数值模拟结果与试验结果进行对比分析。研究结果表明：当冲击载荷来源于拱顶时，巷道围岩破坏主要为裂纹起裂、扩展，拱底发生拉伸破坏；当冲击载荷来源于侧帮时，巷道围岩仅在两侧拱肩处发生破坏。巷道围岩内预制裂纹的动态起裂韧度与冲击加载方向偏向角有关。在相同动态载荷作用下，巷道断面拱肩与拱脚压应力较大，容易产生拉伸破坏，动态稳定性较小。

关键词:巷道；复合型裂纹；冲击加载方向；起裂时间；破坏方式；动态起裂韧度

中图分类号:TU45    文献标志码:A     文章编号:1672-7207（2019）08-1921-09

Failure behaviors of cracked tunnel under impact from different loading orientations

ZHOU　Lei^{1, 2}, ZHU　Zheming^{1, 2}, FAN　Yong^{1, 2}, LIU　Ruifeng^{1, 2}, DENG　Shuai^{1, 2}, DONG　Yuqing^{1, 2}

(1. Key Laboratory of Deep Underground Science and Engineering of Ministry of Education, Sichuan University, Chengdu 610065, China；

2. College of Architecture and Environment, Sichuan University, Chengdu 610065, China)

Abstract: In order to investigate the effect of different impact loading orientations on the failure patterns of surrounding rock mass of cracked tunnel, green sandstone was chosen as the experimental material to simulate the tunnel, and impact experiments were conducted. A crack emanating from tunnel roof was set as the model to simulate the natural crack around a tunnel, and then the drop weight testing machine was used as the loading device. In the process of dynamic loading test, strain gauges were applied to determine crack initiation time and crack propagation velocity, and then corresponding numerical simulation were carried out. The properties of crack propagation and the failure mode of tunnel model were obtained. Finally, the numerical results were compared with the test results. The results show that when the impact loads come from the tunnel roof, the failure of the tunnel mainly includes crack initiation and propagation, and tensile failure of the tunnel bottom. When the impact loads come from the sidewall, the failure of the tunnel occurs at the arch shoulder of two sides. The dynamic initiation toughness of the pre-crack in the surrounding rock mass is greatly related to the impact loading direction. Under the same impact loading, the stress of the arch shoulder and arch foot of tunnel are larger than that of other regions, the arch shoulder and arch foot are liable to tensile failure, and the dynamic stability is relatively low.

Key words: tunnel; mixed-mode crack; impact loading direction; crack initiation time; failure mode; dynamic initiation toughness

在煤矿开采过程中，巷道围岩随时可能会受到各种冲击载荷作用，例如落石冲击、掘进爆破、顶板剥落等，这些动力载荷将诱发巷道围岩周围的天然裂纹与节理等缺陷的孕育、扩展与贯通，最终导致巷道发生破坏与失稳坍塌。当来自不同方向的冲击载荷加载于巷道断面时，将会导致巷道断面发生不同形式的破坏，因此，研究不同方向冲击载荷作用下含预制缺陷巷道的破坏方式具有重要意义。人们对巷道在爆炸与冲击载荷等作用下的破坏行为进行了试验与数值分析研究。ZHU等^[1]采用岩石破坏过程分析(RFPA)软件对巷道围岩周围的动态扰动载荷作用进行了分析，并对侧压力系数和动态扰动波形对地下洞室周围破坏区发育的影响进行了数值模拟。郭东明等^[2-3]采用爆炸加载动态焦散线测试系统研究了爆炸载荷作用下倾斜裂纹对巷道围岩破坏机理的影响。王建国等^[4]采用自制砂浆材料研究了不同节理倾角对冲击应力波的动态响应，并对在节理交界面上冲击应力波的波动特性、能量传递和耗散进行了分析。朱晶晶等^[5]利用大直径霍普金森压杆(SHPB)试验装置对砂岩在单次冲击与重复多次冲击载荷作用下的动力学特性进行了分析，从岩石的细观裂纹扩展和能量吸收的角度分析了岩石的破坏过程。巷道围岩岩体内部含有各种天然节理等缺陷，为了更好地分析巷道围岩内裂纹缺陷对围岩强度评价的影响，采用试验-数值法对含裂纹巷道模型中裂纹尖端的动态应力强度因子进行计算，随后根据试验测定的起裂时间确定动态起裂韧度。WANG等^[6-7]采用裂纹扩展计、应变片并结合试验-数值法计算了平台巴西圆盘在冲击载荷作用下的动态起裂韧度与扩展韧度，并分析了其随裂纹扩展速度的变化规律。FAYE等^[8]研究了三点弯曲有机玻璃板(PMMA)试件在冲击载荷作用下动态起裂韧度的变化趋势，分析了动态加载率对动态起裂韧度的影响。ZHANG等^[9-10]采用高速摄影系统与裂纹扩展计测试系统分析了含缺口半巴西圆盘试件(NSCB)在动载荷作用下裂纹扩展速度的变化趋势，并根据准静态法与数字图像相关法则计算了相应的动态断裂韧度。ZHAO等^[11-13]利用高速摄影采集系统研究了爆炸载荷作用下裂纹尖端应力场的演化规律，并对裂纹尖端的动态断裂韧度进行计算。关于在冲击载荷作用下缺陷岩石内裂纹扩展规律的研究，采用数值模拟手段实现冲击加载过程的重现是一个较好的研究方式，因此，本文作者采用有限差分法AUTODYN软件进行动态断裂行为的数值模拟，该软件能够很好地解决几何非线性和材料非线性动态相关问题^[14-16]。ZHU等^[17-20]采用AUTODYN软件研究了爆炸载荷作用下巷道断面对缓慢卸载P波的动态响应，并验证了该软件计算结果的有效性。WONG等^[21]采用累积损伤模型研究了冲击载荷作用下含共线双裂纹岩石试件的破坏行为，并对裂尖的破坏应力场及破坏规律进行了分析。目前针对巷道围岩在动态载荷作用下的相关研究尚处于初步阶段，对于不同方向加载对巷道围岩破坏的影响，人们大多采用准静态法进行研究^[22-24]。本文基于霍普金森压杆试验原理，采用落锤冲击试验机进行巷道模型试件的冲击加载试验，研究不同冲击加载方向对含预制裂纹缺陷巷道模型破坏方式的影响，并采用应变计测定起裂时间及裂纹扩展速度；采用AUTODYN与ABAQUS软件并结合试验-数值法，验证巷道模型破坏行为并计算动态起裂韧度，分析不同方向冲击载荷作用对巷道断面稳定性的影响。

1 试验方案设计

1.1　试验设备

采用自行研制的中低速落锤冲击试验机进行动态加载，试验装置见图1。冲击试验机主要包含5个部分即冲击锤、入射板、透射板、数据采集系统和阻尼器。落锤冲击试验机依据SHPB原理进行安装，入射板与透射板均采用LY12CZ铝合金。该材料的弹性模量为71.7 GPa，泊松比为0.3，密度为2 850 kg/m³，基于一维应力波传播理论，理论纵波波速为5 015.76 m/s，实测纵波波速为5 006 m/s。入射板长度l_i=3 m，透射板长度l_t=2 m，冲击锤的高度h_s=300 mm，三者宽度×厚度均为300 mm×30 mm，入射板上粘贴有3张应变片(BX120-20AA)，距试件上端面的距离分别为s_i1=2500 mm，s_i2=1500 mm，s_i3=500 mm。经过多次试验验证，最后选择s_i3作为入射端信号采集端。在透射板中间位置s_t=1 m处粘贴应变片作为透射端信号的采集。

图1　落锤冲击试验机

Fig. 1　Drop weight testing apparatus

冲击锤下落高度可以根据试验设计方案进行设置，通过绳索将冲击锤向上提升，提升高度范围为0~10 m，高度采用红外线测距仪进行测试。随后根据下式计算:

(1)

式中：v为冲击锤下落的冲击速度；g为重力加速度 (g= 9.8 m/s²)；h_s为冲击锤冲击高度。

为了获得较为理想的试验结果，在冲击锤与入射板之间粘贴铜合金片作为波形整形器。经过前期的试验对比，该措施能够延长加载波形的上升沿，减小波传播过程中的弥散效应^[25]。为了减小巷道模型试件与入射板、透射板之间的摩擦效应对动态压缩试验的影响，在试件与入射板、透射板之间涂抹少量的凡士林作为润滑剂。

1.2　试件的制作

为了研究不同冲击加载方向对巷道模型断裂行为的影响，采用颗粒均匀、结构稳定的青砂岩作为模型材料来制作试件，将巷道模型理想化为连续介质力学模型，该材料力学参数如下：动态泊松比μ_d=0.165，动态弹性模量 E_d=13.58 GPa，单轴抗压强度σ_U=22.08 MPa，单轴抗拉强度σ_t=2.08 MPa，密度为ρ=2 265 kg/m³，膨胀波波速 c_d＝2 563 m/s，畸变波波速 c_s＝1 639 m/s，瑞雷波波速 c_R＝1 457.6 m/s，模型长×宽×高为350 mm×300 mm×30 mm。为了模拟巷道围岩内含有天然缺陷等情况，在巷道断面对称轴线上设置预制单裂纹。巷道断面与加载方向的夹角为θ，巷道高度h = 60 mm，宽度w = 50 mm，圆弧半径r = 25 mm，裂纹长度2a = 50 mm。巷道模型示意图如图2所示。试件首先采用高压水刀进行冲切，形成含巷道断面形状的矩形岩石板材，随后采用厚度为1 mm的钢锯条进行裂纹粗加工，加工裂纹长度至49 mm，最后采用厚度为0.2 mm的美工刀片进行1 mm裂纹尖端的锐化，使试件能够较好地满足国际岩石力学学会对裂纹宽度的制作要求。巷道模型试件以15°为梯度增量，依次设置θ为0°，15°，30°，45°，60°，75°和90°的7组试验模型，每组试验制作3个试件，共计21个巷道模型试件。

图2　巷道模型示意图

Fig. 2　Diagram of tunnel model

IRWIN^[26]按力学性质特征将裂纹分为张开型裂纹、滑开型裂纹和撕开型裂纹这3种基本类型。为了区分巷道模型内预制裂纹的破坏模式，基于IWRIN断裂力学理论(大部分结构主要是张开型裂纹与滑移型裂纹破坏)，引入拉伸应力强度因子K_I与剪切应力强度因子K_II对预制裂纹的断裂模式进行区分，随后基于图2中的巷道模型试件建立相应的数值模型，结合ABAQUS软件对裂纹尖端的准静态载荷作用下应力强度因子进行计算，设置加载应力为22.08 MPa，设置模型材料参数与砂岩参数一致，最终得到裂纹尖端K_I和K_II随静态载荷加载方向(即倾角θ)的变化，如图3所示(其中负值表示裂纹相对错动方向)。从图3可以看出：当θ = 0°时，K_II和K_I 均为0 MPa·m^1/2 ，表明此时为纯I型裂纹；当巷道断面倾斜角θ增大时，K_II先逐渐增大，随后减小，而K_I绝对值也逐渐增大，说明此时裂纹断裂模式为I/II复合型裂纹，需要考虑I/II复合型断裂模式对巷道动态稳定性的影响。

图3　准静态应力强度因子随倾角θ的变化

Fig. 3　Change of quasi-static stress intensity factor with θ

1.3　试验数据分析

首先对冲击试验机采集的应变数据进行分析，利用入射板与透射板上粘贴的应变片，采集得到巷道断面倾角θ = 45°，冲击速度v = 8.271 m/s时入射端与透射端的应变波形，如图4所示。从图4可以得知入射波与反射波能够较好地分离开，且透射波也能够很好地进行数据采集，波形试验数据与SHPB试验装置基本一致^[27-29]，说明冲击试验机能够实现岩石材料在中低速冲击载荷作用下的动态加载试验。

为了准确地测定巷道内预制裂纹的起裂时间及裂纹扩展速度等参数，并为后期的动态起裂韧度计算提供前提条件，在裂纹尖端处及裂纹扩展路径中粘贴2张应变片即BX120-50AA SG₁和SG₂，2张应变片之间的距离s=20 mm。该方法已经被广泛用于测定预制裂纹起裂时间^[30-31]。当预制裂纹开始起裂时，应变片SG₁将被拉断，应变片两端的电压将突然增大，形成台阶电压，进而可以精确地确定裂纹的起裂时间t_i；当应变片SG₂被拉断时，同理可以得到裂纹扩展20 mm的扩展时间t_p，从而可以初步确定裂纹的扩展速度v_p=20/(t_p-t_i)。应变片电压随时间的变化如图5所示。从图5可以看出：巷道内预制裂纹的起裂时间t_i=262.52 μs，裂纹的扩展时间t_p=315.82 μs，裂纹的平均扩展速度v_p= 375.23 m/s。

图4　入射板与透射板上应变随时间的变化

Fig. 4　Change of strain of the incident and transmission plates with time

图5　应变片电压随时间的变化

Fig. 5　Change of strain gauge voltage with time

2 动态断裂韧度计算

由于巷道围岩内通常含有较多的天然裂纹与节理等缺陷，这将极大地影响巷道结构的动态稳定性，而岩石材料的动态断裂韧度是评价岩石在冲击载荷作用下抵抗断裂性能的重要依据^[32-33]，因此，巷道模型结构围岩内部裂纹尖端动态起裂韧度的计算就变得尤为重要。

2.1　动态载荷计算

采用试验-数值法计算巷道模型试件内裂纹的动态起裂韧度时，为了达到与试验荷载相同的加载条件，首先对模型试件上、下两端的应力分别进行计算，然后根据SHPB试验数据处理原则，对图4中的应变进行应力转换^[34-36]：

(2)

(3)

式中：σ₁(t)为试件上端所受载荷，可由入射波应变时程曲线ε_i(t)与反射波ε_r(t) 应变时程曲线叠加得到；σ₂(t)为试件下端所受载荷，可由透射波应变时程曲线ε_t(t)得到；E_i和E_t分别为入射板与透射板的弹性模量；A_i，A_t和A_s分别为入射板、透射板、巷道模型试件的加载横截面积。

模型试件上、下两端的应力荷载随时间的变化如图6所示。将图6中结果代入数值分析软件中进行分析，并结合试验所测的裂纹起裂时间t_i确定动态起裂韧度。

图6　巷道模型试件上、下端的冲击载荷

Fig. 6　Impact loads of tunnel model sample on both ends of the upper and lower sides

图7　2D数值模型网格单元

Fig. 7　Typical mesh of 2D numerical model

图8　动态应力强度因子曲线(θ =0°)

Fig. 8　Dynamic stress intensity factor curves (θ =0°)

2.2　动态应力强度因子计算

本文采用ABAQUS软件结合试验-数值法进行动态应力强度因子的计算。为了验证该软件计算的精准性，对有限板含中心预制单裂纹在Heaviside阶跃函数载荷形式作用下的动态应力强度因子进行计算。选取均匀的时间步长Δt= 0.2 μs，采用隐式的动态有限元数值方法模块(ABAQUS/Standard)进行求解，随后与文献[37]中的有限差分法计算结果及文献[38]中的边界元法计算结果进行对比，发现这3种方法计算结果较一致^[39]，说明ABAQUS软件对于冲击载荷作用下动态应力强度因子的计算具有较高的准确性。

图9 动态应力强度因子曲线(θ =45°)

Fig. 9　Dynamic stress intensity factor curves (θ =45°)

基于图2中的巷道模型尺寸，建立2D有限元数值模型，该模型裂纹尖端采用6节点的CPS6单元，其余区域采用8节点的CPS8单元进行划分，网格单元总数为8 036个。将裂纹理想化为尖锐裂纹，采用1/4节点奇异单元描述裂纹尖端应力场和应变场的奇异性，并在裂纹尖端进行网格细化。2D数值模型网格单元如图7所示。加载时间设为600 μs，每1 μs设置为1个加载步，随后结合试验测试的起裂时间t_i确定巷道模型试件的动态起裂韧度，并对其变化规律进行分析。

分别选择巷道断面倾角θ = 0°与45°这2种工况作为典型研究对象进行分析，两者动态应力强度因子曲线分别如图8和图9所示。从图8可以看出：动态应力强度因子K_I变化较大，K_II变化较小且接近0 MPa·m^1/2 ，说明此时裂纹为纯I型裂纹破坏方式，这与第1.2节中裂纹断裂模式验证结果一致。结合试验测试的起裂时间t_i，可确定此时巷道模型试件的动态起裂韧度=1.63 MPa·m^1/2，该值小于动态应力强度因子曲线的最大值K_IMAX (2.29 MPa·m^1/2)，说明裂纹在并未达到拉伸应力最大值时就已经起裂，随后继续向前扩展；而在冲击载荷的动态惯性效应作用下，裂纹开始阶段出现在150 μs左右，此前K_I始终为负值，说明裂纹初始阶段处于压缩闭合状态，与准静态载荷作用下的应力场分布有较大区别。

图9所示为巷道断面倾斜角θ = 45°时的动态应力强度因子时程曲线。从图9可见：前期K_I与K_II均为负值，后变为正值；在应力波传至裂纹尖端处，裂纹开始处于压剪状态，随后处于拉剪状态，经过压缩与剪切共同作用促使裂纹尖端起裂，起裂时动态起裂韧度=-6.64 MPa·m^1/2，=-6.53 MPa·m^1/2，两者都小于各自时程曲线中动态应力强度因子绝对值的最大值，说明当冲击加载来源于不同方向时，巷道内预制裂纹缺陷断裂模式将发生改变，特别是剪切作用的影响将随之增大。

不同倾角θ下动态起裂韧度分别如表1所示。由于巷道断面倾斜θ = 75°与θ = 90°时，主要在巷道围岩两侧拱肩部位发生破坏，而预制裂纹没有发生起裂，因此，裂纹的起裂时间、裂纹平均扩展速度、动态起裂韧度无法进行相应地计算。

表1　不同倾角θ下动态起裂韧度

Table 1　Dynamic initiation toughness of different dip angles

3 巷道模型破坏行为分析

为了验证试验结果的准确性，对不同冲击载荷方向作用下的巷道模型进行数值模拟，采用有限差分法软件AUTODYN进行数值分析，内置显式非线性动力学分析系统。该软件能够解决动态载荷作用下几何非线性和材料非线性的动态断裂问题，内嵌多种材料数据库和动态本构模型，可满足爆炸与冲击载荷等作用下的研究要求。

基于图2中的巷道模型尺寸，建立7组相应的三维有限元数值模型，如图10所示。模型参数与本文第1.2节中的砂岩参数一致。模型中建立了冲击锤、试件、入射板、透射板、阻尼器等部件，以充分模拟试验条件下的载荷情况。整个数值模型含有153 310个六面体单元，网格单元边长为1 mm，冲击锤的下落冲击速度根据试验条件进行设置。

图10　3D数值模型网格单元示意图

Fig. 10　Diagram of 3D numerical model mesh unit

3.1　巷道围岩破坏行为分析

当θ分别为0°，45°和90°时，不同方向冲击载荷的数值模拟结果与试验结果对比如图11所示。从图11可知：当θ = 0°时，巷道围岩内的破坏方式主要是由预制裂纹的起裂及巷道底板处拉伸应力造成的拉伸破坏；当θ = 45°时，巷道围岩的破坏方式主要是预制裂纹的起裂、扩展及拱脚的破坏，说明此时拱脚区域稳定性较小；由于当θ = 45°时裂纹为I/II复合型裂纹，起裂以后在剪应力与拉应力共同作用下，裂纹走向逐渐偏向于最大主应力区，即巷道模型试件的对称轴区域。当θ = 90°时，巷道断面主要是两侧拱肩的破坏，而预制裂纹没有发生起裂及扩展，说明此时天然裂纹缺陷的影响较小，主要考虑巷道断面拱肩与拱脚区域对应力波的动态响应。由图11可知数值分析结果与试验结果较吻合，说明试验结果具有较高的准确性。总体来说，当冲击载荷来源于不同加载方向时，巷道断面的破坏形式也不同，需根据冲击载荷的方向来加强巷道结构的动态稳定性。

3.2　巷道围岩应力分析

为了分析冲击载荷作用对巷道围岩内部的影响，对巷道断面倾角θ = 0°时围岩周围的应力进行分析。在巷道断面周围设置高斯观测点，观测巷道围岩周围应力的变化，如图12所示。从图12可知：巷道的拱肩与拱脚压应力比其他区域的大，易造成拉伸破坏；在底板位置处，由于受到拉伸应力作用造成拉伸破坏，形成拉伸裂纹，应力突然降为零，与试验结果较为吻合；而在高斯观测点①和③处，压应力相对较小，说明观测点①和③处的巷道区域在竖向动态载荷作用下相对稳定，这与准静态加载作用下巷道围岩周围的应力分布情况不同^[40-41]。

图11　不同倾角θ下数值模拟结果分析

Fig. 11　Numerical simulation results of different dip angles

图12　巷道围岩测点分布及周围应力随时间的变化

Fig. 12　Distribution of monitoring points around the tunnel and change of surrounding stress of rock mass with time

4 结论

1) 当加载不同方向的冲击载荷时，巷道围岩的破坏形式有很大差别。当冲击载荷来源偏向于侧帮方向时，巷道围岩仅仅两侧拱肩破坏；当冲击载荷来源偏向于拱顶时，巷道断面主要是预制裂纹的起裂、扩展及拱脚或底板的破坏。

2) 当冲击载荷来源于不同方位时，围岩内天然预制裂纹缺陷的断裂模式也发生较大变化，从纯I型裂纹逐渐变为I/II复合型裂纹的破坏，并且动态断裂韧度与冲击载荷的加载方向有很大的关系。

3) 当冲击载荷来源于竖直方向时，巷道断面的拱肩与拱脚位置处的压应力较大，说明动态稳定性较小，并且在拱底容易受到拉伸应力作用而产生拉伸破坏，这与静态载荷作用下的破坏行为有较大的区别。

参考文献：

[1] ZHU W C, LI Z H, ZHU L, et al. Numerical simulation on rockburst of underground opening triggered by dynamic disturbance[J]. Tunnelling and Underground Space Technology, 2010, 25(5): 587-599.

[2] 郭东明, 闫鹏洋, 杨仁树, 等. 动静荷载下巷道围岩倾斜裂纹的动焦散试验[J]. 采矿与安全工程学报, 2016, 33(4): 668-675.

GUO Dongming, YAN Pengyang, YANG Renshu, et al. Dynamic caustics test on the inclined cracks in roadway surrounding rock subjected to static-dynamic load[J]. Journal of Mining & Safety Engineering, 2016, 33(4): 668-675.

[3] 郭东明, 刘康, 胡久羡, 等. 爆生气体对邻近硐室背爆侧预制裂纹影响机理[J]. 煤炭学报, 2016, 41(1): 265-270.

GUO Dongming, LIU Kang, HU Jiuxian, et al. Experimental study on the effect mechanism of the explosion gas on the precrack in the back-blasting of adjacent tunnel[J]. Journal of China Coal Society, 2016, 41(1): 265-270.

[4] 王建国, 梁书锋, 高全臣, 等. 节理倾角对类岩石冲击能量传递影响的试验研究[J]. 中南大学学报(自然科学版), 2018, 49(5): 1237-1243.

WANG Jianguo, LIANG Shufeng, GAO Quanchen, et al. Experimental study of jointed angles impact on energy transfer characteristics of simulated rock material[J]. Journal of Central South University(Science and Technology), 2018, 49(5): 1237-1243.

[5] 朱晶晶, 李夕兵, 宫凤强, 等. 冲击载荷作用下砂岩的动力学特性及损伤规律[J]. 中南大学学报(自然科学版), 2012, 43(7): 2701-2707.

ZHU Jingjing, LI Xibing, GONG Fengqiang, et al. Experimental test and damage characteristics of sandstone under uniaxial impact compressive loads[J]. Journal of Central South University(Science and Technology), 2012, 43(7): 2701-2707.

[6] WANG Q Z, YANG J R, ZHANG C G, et al. Sequential determination of dynamic initiation and propagation toughness of rock using an experimental–numerical–analytical method[J]. Engineering Fracture Mechanics, 2015, 141: 78-94.

[7] WANG Q Z, ZHANG S, XIE H P. Rock dynamic fracture toughness tested with holed-cracked flattened Brazilian discs diametrically impacted by SHPB and its size effect[J]. Experimental Mechanics, 2010, 50(7): 877-885.

[8] FAYE A, PARAMESWARAN V, BASU S. Dynamic fracture initiation toughness of PMMA: a critical evaluation[J]. Mechanics of Materials, 2016, 94: 156-169.

[9] ZHANG Q B, ZHAO J. Effect of loading rate on fracture toughness and failure micromechanisms in marble[J]. Engineering Fracture Mechanics, 2013, 102: 288-309.

[10] ZHANG Q B, ZHAO J. Determination of mechanical properties and full-field strain measurements of rock material under dynamic loads[J]. International Journal of Rock Mechanics and Mining Sciences, 2013, 60: 423-439.

[11] ZHAO Yixin, GONG Shuang, HAO Xianjie, et al. Effects of loading rate and bedding on the dynamic fracture toughness of coal: Laboratory experiments[J]. Engineering Fracture Mechanics, 2017, 178: 375-391.

[12] YUE Zhongwen, SONG Yao, YANG Renshu, et al. Comparison of caustics and the strain gage method for measuring mode I stress intensity factor of PMMA material[J]. Polymer Testing, 2017, 59: 10-19.

[13] WANG Yanbing, YANG Renshu, ZHAO Gaofeng. Influence of empty hole on crack running in PMMA plate under dynamic loading[J]. Polymer Testing, 2017, 58: 70-85.

[14] JIANG Hongxiang, DU Changlong, LIU Zenghui. Theoretical and numerical investigation on rock fragmentation under high-pressure water-jet impact[J]. Iranian Journal of Science and Technology, Transactions of Civil Engineering, 2017, 41(3): 305-315.

[15] DEHGHAN BANADAKI M M, MOHANTY B. Numerical simulation of stress wave induced fractures in rock[J]. International Journal of Impact Engineering, 2012, 40/41: 16-25.

[16] JIANG Hongxiang, DU Changlong, LIU Songyong, et al. Numerical simulation of rock fragmentation under the impact load of water jet[J]. Shock and Vibration, 2014, 2014: 1-11.

[17] ZHU Zheming. Numerical prediction of crater blasting and bench blasting[J]. International Journal of Rock Mechanics and Mining Sciences, 2009, 46(6): 1088-1096.

[18] ZHU Zheming, MOHANTY B, XIE Heping. Numerical investigation of blasting-induced crack initiation and propagation in rocks[J]. International Journal of Rock Mechanics and Mining Sciences, 2007, 44(3): 412-424.

[19] ZHU Zheming, WANG Chao, KANG Jiming, et al. Study on the mechanism of zonal disintegration around an excavation[J]. International Journal of Rock Mechanics and Mining Sciences, 2014, 67: 88-95.

[20] ZHU Zheming, XIE Heping, MOHANTY B. Numerical investigation of blasting-induced damage in cylindrical rocks[J]. International Journal of Rock Mechanics and Mining Sciences, 2008, 45(2): 111-121.

[21] WONG L N Y, LI Huanqiang. Numerical study on coalescence of two pre-existing coplanar flaws in rock[J]. International Journal of Solids and Structures, 2013, 50(22/23): 3685-3706.

[22] ZHU Zheming, LI Yuanxin, XIE Jun, et al. The effect of principal stress orientation on tunnel stability[J]. Tunnelling and Underground Space Technology, 2015, 49: 279-286.

[23] 李永乾, 朱哲明, 胡荣. 主应力方向对隧道稳定性影响规律的研究[J]. 四川大学学报(工程科学版), 2012, 44(S1): 93-98.

LI Yongqian, ZHU Zheming, HU Rong. Study on the effect of principal stress orientation on tunnel stability[J]. Journal of Sichuan University(Engineering Science Edition), 2012, 44(S1): 93-98.

[24] 李元鑫, 朱哲明, 范君黎. 主应力方向对围岩稳定性的影响[J]. 岩土工程学报, 2014, 36(10): 1908-1914.

LI Yuanxin, ZHU Zheming, FAN Junli. Effect of principal stress orientation on stability of surrounding rock of tunnels[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(10): 1908-1914.

[25] FREW D J, FORRESTAL M J, CHEN W. Pulse shaping techniques for testing elastic-plastic materials with a split Hopkinson pressure bar[J]. Experimental Mechanics, 2005, 45(2): 186-195.

[26] IRWIN G R. Analysis of stresses and strains near the end of a crack traversing a plate[J]. Journal of Applied Mechanics. 1957, 24: 109-114.

[27] 李地元, 成腾蛟, 周韬, 等. 冲击载荷作用下含孔洞大理岩动态力学破坏特性试验研究[J]. 岩石力学与工程学报, 2015, 34(2): 249-260.

LI Diyuan, CHENG Tengjiao, ZHOU Tao, et al. Experimental study of the dynamic strength and fracturing characteristics of marble specimens with a single hole under impact loading[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(2): 249-260.

[28] WANG Xiaomeng, ZHU Zheming, WANG Meng, et al. Study of rock dynamic fracture toughness by using VB-SCSC specimens under medium-low speed impacts[J]. Engineering Fracture Mechanics, 2017, 181: 52-64.

[29] ZHOU Lei, ZHU Zheming, WANG Meng, et al. Dynamic propagation behavior of cracks emanating from tunnel edges under impact loads[J]. Soil Dynamics and Earthquake Engineering, 2018, 105: 119-126.

[30] ZHANG Q B, ZHAO J. A review of dynamic experimental techniques and mechanical behaviour of rock materials[J]. Rock Mechanics and Rock Engineering, 2014, 47(4): 1411-1478.

[31] 杨井瑞, 张财贵, 周妍, 等. 用SCDC试样测试岩石动态断裂韧度的新方法[J]. 岩石力学与工程学报, 2015, 34(2): 279-292.

YANG Jingrui, ZHANG Caigui, ZHOU Yan, et al. A new method for determining dynamic fracture toughness of rock using SCDC specimens[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(2): 279-292.

[32] GUI Yilin, BUI H H, KODIKARA J, et al. Modelling the dynamic failure of brittle rocks using a hybrid continuum-discrete element method with a mixed-mode cohesive fracture model[J]. International Journal of Impact Engineering, 2016, 87: 146-155.

[33] ZHU Z M, XU W T, FENG R Q. A new method for measuring mode-I dynamic fracture toughness of rock under blasting loads[J]. Experimental Techniques, 2016, 40(3): 899-905.

[34] YANG Renshu, CHEN Jun, YANG Liyun, et al. An experimental study of high strain-rate properties of clay under high consolidation stress[J]. Soil Dynamics and Earthquake Engineering, 2017, 92: 46-51.

[35] YAO Wei, LIU H W, XU Ying, et al. Thermal degradation of dynamic compressive strength for two mortars[J]. Construction and Building Materials, 2017, 136: 139-152.

[36] 谢晓锋, 陶明, 吴秋红, 等. 不同形状砂岩的动态力学特性[J]. 中南大学学报(自然科学版), 2017, 48(9): 2441-2448.

XIE Xiaofeng, TAO Ming, WU Qiuhong, et al. Dynamic properties of sandstones with different shapes[J]. Journal of Central South University(Science and Technology), 2017, 48(9): 2441-2448.

[37] CHEN Y M. Numerical computation of dynamic stress intensity factors by a Lagrangian finite-difference method (the HEMP code)[J]. Engineering Fracture Mechanics, 1975, 7(4): 653-660.

[38] FEDELINSKI P, ALIABADI M H, ROOKE D P. Boundary element formulations for the dynamic analysis of cracked structures[J]. Engineering Analysis with Boundary Elements, 1996, 17(1): 45-56.

[39] 解德, 钱勒, 李长安. 断裂力学中的数值计算方法及工程应用.北京, 科学出版社; 2009:134-145.

XIE De, QIAN Le, LI Changan. Numerical calculation method and engineering application in fracture mechanics[M]. Beijing: Science Press, 2009: 134-135.

[40] 周磊, 朱哲明, 刘邦. 隧道周边不同位置径向裂纹对隧道围岩稳定性影响规律的研究[J]. 岩土工程学报, 2016, 38(7): 1230-1237.

ZHOU Lei, ZHU Zheming, LIU Bang. Influence of radial cracks on stability of surrounding rocks at different locations around tunnel[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(7): 1230-1237.

[41] 周磊, 朱哲明, 刘邦. 裂纹对直墙拱形隧道围岩损伤破坏模式的影响规律研究[J]. 岩土力学, 2017, 38(12): 3688-3697.

ZHOU Lei, ZHU Zheming, LIU Bang. Influence of cracks on surrounding rock damage-failure mode of straight wall arch tunnel[J]. Rock and Soil Mechanics, 2017, 38(12): 3688-3697.

(编辑  伍锦花)

收稿日期： 2018 -10 -18; 修回日期： 2018 -12 -18

基金项目(Foundation item)：国家自然科学基金资助项目(11672194)；四川省安全监管局安全生产科技项目(aj20170515161307)；四川省科技计划项目(2018JZ0036);(Project(11672194) supported by the National Natural Science Foundation of China; Project(aj20170515161307) supported by the Science and Technology Program of Sichuan Administration of Work Safety; Project (2018JZ0036) supported by the Science and Technology Plan of Sichuan Province)

通信作者：朱哲明，教授，博士生导师，从事岩石力学与采矿工程研究；E-mail:zhemingzhu@hotmail.com