Experimental study on influence of granite grain size on rockburst
-
摘要: 为探究晶粒尺寸对硬脆性岩石岩爆的影响,利用真三轴岩爆试验系统,对细中、中粗两种不同晶粒尺寸的含预制圆孔花岗岩开展了岩爆模拟试验。试验结果表明:在相同的加载过程中,细中晶粒花岗岩出现板裂静态脆性破坏,而中粗晶粒花岗岩出现岩爆动力破坏;细中晶粒花岗岩早期声发射活动较弱,大破裂、低主频事件在时空分布上较集中,特征应力较大,而中粗晶粒花岗岩早期声发射活动较活跃,大破裂、低主频事件在时空分布上较离散,特征应力较小,碎屑破碎程度更高。晶粒尺寸对花岗岩的岩爆倾向性具有重要影响,晶粒尺寸较大的硬脆性岩石的岩爆倾向性更强。深部地下岩体工程的岩爆倾向性评价中,除强度和脆性外,晶粒尺寸也是需要考虑的重要因素。Abstract: In order to explore the influence of granite grain size on rockburst, cubic granite specimens with an opening and different grain sizes (fine to medium and medium to coarse) were used to conduct the rockburst tests using the true triaxial rockburst testing system. The experimental results show that the failure process of the fine to medium-grained granite is mainly composed of brittle failure. However, rockburst failure (dynamic failure) dominates the failure for the medium to coarse-grained granite. The acoustic emission (AE) activity in the early loading stage is weak for the fine to medium-grained granite, and the low-frequency large-rupture events are concentrated in time and space, and the characteristic stress is higher. However, the AE activity in the early loading stage is stronger for the medium to coarse-grained granite, and the low-frequency, large-rupture events are more discrete in time and space, and the characteristic stress is lower, and the fragments are broken more. The grain size has an important influence on the rockburst proneness of granites. The hard-brittle rock with coarser grain size has a stronger rockburst proneness. In addition to strength and brittleness, grain size is an important factor to be considered in rockburst proneness evaluation of deep underground rock mass engineering.
-
Key words:
- rockburst /
- spalling /
- granite /
- grain size
-
铝蜂窝材料是一种二维多孔材料, 具有密度小、成本低、变形能力强的特点, 是一种优良的吸能材料, 广泛应用于各种轻质结构和缓冲吸能结构中[1-3]。
近年来, 很多学者对铝蜂窝的静态变形规律做了深入研究, 但对其动态力学性能仍然了解较少。已有研究表明, 当基体材料本身率不敏感时, 蜂窝材料在面外动态加载下仍存在明显的强化效应。W.Goldsmith[4]对蜂窝结构进行了动态压缩实验, 测得蜂窝动态强化率为20%~50%。H.Zhao等[5]对铝蜂窝进行了SHPB实验, 测得其动态强化率约40%。对于多孔材料的动态强化, 目前学界主要有3种解释:封闭气体压缩理论[2], 冲击波理论[6-7]和横向惯性效应理论[8-10]。其中, 横向惯性效应理论可以较好地解释中等速度冲击下, 蜂窝材料动态强化这一现象。该理论认为, 结构在动态加载下, 由于惯性的存在, 其横向运动需要一定的时间才会与纵向加载的速度相协调, 这种运动的不协调性导致结构的纵向压缩应变较静态时有所增加。具有应变硬化特性的材料则会出现强化现象。
为了进一步研究基体材料对铝蜂窝动态力学性能的影响, 验证横向惯性效应理论, 本文中从铝蜂窝结构中取出蜂窝壁, 对基体材料进行力学性能测试。对几何参数相同而基体材料不同的2种铝蜂窝分别进行单轴面外动态和静态压缩实验。最后结合横向惯性效应理论分析基体材料对铝蜂窝动态强化行为的影响。
1. 实验方法
选择3003H18和5052H18铝蜂窝进行基体材料实验、蜂窝的静态压缩和动态压缩实验。2种铝蜂窝的几何尺寸完全相同, 高度30 mm, 蜂窝芯格平均边长a=3 mm, 铝箔厚度t=60 μm。
1.1 基体材料实验
构成蜂窝的铝箔与块体铝合金的力学性能有很大不同, 用块体铝合金的本构关系描述蜂窝基体材料力学行为存在较大误差。因此有必要获得蜂窝壁的本构关系, 为描述蜂窝结构的力学行为奠定基础。
从2种铝蜂窝上分别截取蜂窝壁进行静态拉伸实验。考虑到线切割加工可能会影响蜂窝壁的力学性能, 从蜂窝结构上沿壁板交界的折痕取下矩形铝箔试样用于实验, 结果表明, 其精度完全满足实验要求。截取的蜂窝壁为矩形, 长30 mm、宽3 mm、厚60 μm。由于蜂窝壁太薄, 制成哑铃状试样难以保证精度, 直接拉伸矩形铝箔又会因夹具夹持造成应力集中, 所以实验难度较大。此外, 夹持式引伸仪因试样太薄无法使用, 视频引伸仪也因试样太短而无法直接使用。考虑到以上问题, 将截取的铝蜂窝壁两端分别粘接2个与其宽度相同的铝箔, 制成拉伸试样, 用粘接的方式代替夹具直接对试样的夹持, 如图 1所示。粘接而成的试样结构紧凑, 不受夹具的遮挡, 因此可以利用视频引伸仪来测定试样的应变。实验机采用INSTRON5848型微拉伸实验机, 加载速度为0.1 mm/min(向上运动)。
1.2 蜂窝压缩实验
将蜂窝静态压缩试样加工成六棱柱结构, 该试样共有19个封闭的胞元, 形状如图 2所示。试验机采用CSS88010万能电子试验机。加载速度设为0.1 mm/s(向下运动)。
铝蜂窝动态压缩试样与静态试样完全相同。采用SHPB系统进行测试, 平均加载速度分为10 m/s和28 m/s两种。SHPB系统由发射系统、撞击杆、入射杆、透射杆和信号采集系统组成。试样处于入射杆和透射杆之间。撞击杆以一定的速度撞击入射杆。应力波经入射杆传播后在入射杆端发生反射和透射, 贴在杆子上的应变片可以记录下入射、透射和反射的信号。
铝蜂窝材料分散性较大, 为使试样截面内包含多个胞元, 试样截面要足够大, 因此实验采用大直径压杆。另外, 铝蜂窝材料的波阻抗比常规压杆材料的波阻抗低很多, 因此采用低波阻抗的PMMA压杆, 以提高试样和压杆波阻抗的匹配性。杆子直径30 mm, 入射杆和透射杆的长度均为2 m。为了在有限的发射气压下完成较高速度的加载, 选取290 mm长的撞击杆进行实验, 结果表明采用该长度的撞击杆可以使试样发生足够大的变形, 满足后续研究的需要。通过入射、反射和透射信号得到试样两个端面上的力和速度的信息, 其计算公式如下[5]:
Fin(t)=SbE[εi(t)+εr(t)] (1) Fout (t)=SbEεt(t) (2) vin(t)=c0[εi(t)−εr(t)] (3) vout (t)=c0εt(t) (4) 式中:Fin、Fout、vin和vout分别为试样表面的力和速度。Sb、E和c0分别为压杆的横截面积、杨氏模量和弹性波速。εi(t)、εr(t)和εt(t)分别为入射、反射和透射脉冲。
2. 实验结果分析
2.1 基体材料实验
蜂窝壁拉伸的实验结果如图 3所示。实验结果表明, 5052H18铝箔的强度和破坏应变均高于3003H18铝箔。3003H18的破坏应变在1.3%左右; 5052H18的破坏应变明显高于3003H18, 但表现出一定的分散性, 这可能与5052H18试样边界的粗糙度不同有关。
2.2 蜂窝压缩实验
定义平均应力σ=F/S来反映蜂窝在稳态压溃过程中的载荷大小, 其中F为试验机载荷, S为铝蜂窝试样的横截面积。蜂窝试样横截面积的计算采用等效面积法, 即每一个“Y”形结构对应一个三角形区域的面积, 如图 4所示。等效面积的计算公式为
S=3√3na2/2 (5) 式中:S为等效截面积, n为截面中“Y”型单元的个数, a为六边形边长。19孔蜂窝试样有54个“Y”形单元, 试样等效的截面积为631.8 mm2。后文动态实验也用此方法得到试样的平均应力。由于铝蜂窝的压溃模式为非均匀的渐进式屈曲, 所以本文中以平均应力-位移曲线来描述蜂窝结构的变形行为。对每种材料的每种加载状态都进行了多次实验, 实验结果的一致性较好, 图 5为2种铝蜂窝的平均应力位移曲线。
蜂窝试样胞元较少, 分散性可能比较大, 所以对不同胞元数目的铝蜂窝试样进行了准静态压缩实验, 以验证采用19孔蜂窝作为试样的合理性。我们定义平均应力σ在不同位移x处的平均值为平台应力σp。实验结果(图 6)表明, 胞元数目对铝蜂窝的平台应力σp影响较小。所以, 19孔蜂窝试样在力学性能上具有一般性, 基本可以代表此类蜂窝结构。
实验结果表明, 铝蜂窝结构静态压缩过程分为3个阶段:第1阶段从初始加载点到初始应力峰值。在该阶段内, 蜂窝最初保持挺直的状态, 随着载荷的增大, 出现了整体的弹性或弹塑性屈曲, 此时沿着整个蜂窝高度方向上存在均布的波纹, 当达到初始应力峰值时, 蜂窝局部发生压溃。其本质是结构的失稳。第2阶段为平台段, 蜂窝以渐进的方式发生屈曲(图 7), 压溃区域逐渐扩展。该阶段内, 载荷明显小于初始峰值, 并稳定在某一应力水平附近。需要指出的是, 渐进屈曲和整体的波纹屈曲是2种不同的屈曲模式。第3阶段为压实段, 蜂窝进入密实化阶段, 载荷明显增大。由图 7可以看到, 铝蜂窝在动态和准静态加载下的屈曲模式完全相同(由于动态实验撞击杆的长度有限, 所以动态实验的铝蜂窝没有完全压实)。图 8为铝蜂窝动态压缩实验的波形图(5052H18, 28 m/s)。图 9为2种基体材料蜂窝的动态和准静态实验结果。
从图 9中可以看到, 5052H18铝蜂窝的平台应力高于3003H18铝蜂窝, 这与基体材料实验的结果(图 3)一致。3003H18和5052H18铝蜂窝渐进屈曲产生的褶皱均为12个, 基体材料力学性能的差异并未使褶皱波长出现显著差异。
定义蜂窝的动态强化量Δσ为蜂窝动态加载下的平台应力与准静态加载下的平台应力之差, 动态强化率γ为动态强化量Δσ与静态平台应力之比。综合准静态实验和动态实验的结果, 发现2种铝蜂窝均存在明显的动态强化现象。3003H18铝蜂窝在动态加载下的强化量和强化率都比5052H18铝蜂窝高(表 1)。从实验结果来看, 随着加载速度的增加, 2种铝蜂窝的动态强化率和动态强化量都有所增加。
表 1 不同基体材料蜂窝的平台应力Table 1. Plateau stress of honeycombs with different materials材料/加载速度 σp/MPa Δσ/MPa γ/% 3003H18 (0.1 mm/s) 1.51 3003H18 (10 m/s) 2.21 0.70 46 3003H18 (28 m/s) 2.40 0.89 59 5052H18 (0.1 mm/s) 1.91 5052H18 (10 m/s) 2.44 0.53 28 5052H18 (28 m/s) 2.66 0.75 39 研究表明, 铝合金的率敏感性与其静态流变应力有关[11], 流变应力较大的铝合金几乎没有应变率敏感性。这样, 对于流变应力较大的3003H18和5052H18铝合金, 认为其本身率不敏感。所以我们认为铝蜂窝的动态强化与蜂窝材料的几何结构有关。该现象可以用横向惯性效应理论给予定性的解释, H.Zhao等[9]用刚性杆-塑性铰模型对具有动态强化特性的结构进行了描述, 如图 10所示。其中Mp为塑性铰力矩, θ为初始缺陷角, v为横向位移, σs为材料屈服强度, L为杆长, b、t为杆子矩形截面边长。由力的平衡原理得到杆子塑性失稳的临界外载为
F=Mp|2δθ2δv|=σsbt24Lsinθ (6) 集中质量以加速度a沿x轴正方向运动, 它与杆子的轴向力满足以下关系:
2Nsinθ=ma (7) 在静态压缩下横向运动与纵向运动几何协调。然而在动态加载下, 整个结构沿壁板纵向的压缩载荷要受杆子塑性流动应力大小的限制
ma⩽2σsbtsinθ (8) 这意味着受到横向惯性的影响, 蜂窝壁横向运动对纵向压缩的协调能力有限。所以, 动态压缩下铝蜂窝壁塑性失稳时的纵向应变较静态时有所增大, 即铝蜂窝壁在纵向发生更大的塑性变形。
假设5052H18在塑性阶段的应变硬化率小于3003H18, 而2种铝蜂窝的动态应变增量相同(均为Δε), 那么3003H18铝蜂窝的动态强化量就会高于5052H18铝蜂窝。但由于实验中3003H18铝箔过脆, 所以未能获得足够长的强化段曲线。我们猜想, 在压缩情况下, 3003H18和5052H18符合以下定性的关系, 如图 11所示。另外, 5052H18铝合金的静态强度较3003H18铝合金高, 这使得在应变强化量相同的情况下, 5052H18铝蜂窝的动态强化率要比3003H18铝蜂窝低。
3. 结论
对几何参数完全相同而基体材料分别为3003H18和5052H18铝合金的2种铝蜂窝进行了动态和静态的力学性能测试, 获得了2种蜂窝在面外方向上动态和静态的平台应力, 并得到了各自的动态强化量和动态强化率。同时对以上2种铝蜂窝的基体材料进行了力学性能测试, 得到了其应力应变曲线。利用横向惯性效应模型对2种铝蜂窝的动态强化率的差异进行了分析, 得到如下结论:
(1) 铝蜂窝在面外单轴动态加载下的变形模式与静态加载时相同。
(2) 3003H18铝蜂窝和5052H18铝蜂窝存在明显的动态强化效应。3003H18铝蜂窝的动态强化率(59%)要高于5052H18铝蜂窝的动态强化率(39%)(平均加载速度为28 m/s)。这一现象可能与基体材料的应变硬化率有关。
(3) 基体材料对蜂窝整体力学性能有着重要的影响。5052H18铝蜂窝的平台应力在准静态时比3003H18铝蜂窝高26%, 在28 m/s动态加载下比3003H18高11%。
-
表 1 花岗岩基本物理、力学参数与矿物成分
Table 1. Basic physical and mechanical parameters and mineral composition
属地 晶粒 密度/
(kg·m−3)单轴强度/
MPa弹性模量/
GPa纵波速度/
(km·s−1)矿物成分 晶粒直径/
mm非均匀性 广东肇庆 细中 2 687 135 32.2 4.5 39%斜长石
28%石英
18%钾长石
7%黑云母
3%角闪石
5%其他0.6~5.0 低 广西梧州 中粗 2 680 110 30.6 5.3 54%钾长石
30%石英
10%斜长石
4%黑云母
2%其他2~22 高 表 2 加载过程中细中晶粒花岗岩的特征应力
Table 2. Characteristic stresses in the fine to medium-grained granite specimen during loading
特征应力 时间/s σZ/MPa σθ/MPa σθ/σθmax 起裂应力σci 1 547.64 74.00 212.00 0.55 明显颗粒弹射 2 301.12 113.72 331.16 0.86 损伤应力σcd 2 563.56 127.50 372.50 0.96 明显板裂化 2 650.32 132.08 386.24 1.00 表 3 加载过程中中粗晶粒花岗岩特征应力
Table 3. Characteristic stresses of the medium to coarse-grained granite specimen during loading
特征应力 时间/s 竖向应力σZ/MPa σθ/MPa σθ/σθmax 起裂应力σci 1 145.52 58.00 164.00 0.46 明显弹射 1 838.16 72.76 208.28 0.59 损伤应力σcd 2 241.36 113.00 329.00 0.93 明显岩爆 2 385.72 121.01 353.03 1.00 -
[1] 何满潮, 赵菲, 杜帅, 等. 不同卸载速率下岩爆破坏特征试验分析 [J]. 岩土力学, 2014, 35(10): 2737–2747; 2793. DOI: 10.16285/j.rsm.2014.10.001.HE Manchao, ZHAO Fei, DU Shuai, et al. Rockburst characteristics based on experimental tests under different unloading rates [J]. Rock and Soil Mechanics, 2014, 35(10): 2737–2747; 2793. DOI: 10.16285/j.rsm.2014.10.001. [2] ZHAO X G, WANG J, CAI M, et al. Influence of unloading rate on the strainburst characteristics of Beishan granite under true-triaxial unloading conditions [J]. Rock Mechanics and Rock Engineering, 2014, 47(2): 467–483. DOI: 10.1007/s00603-013-0443-2. [3] 苏国韶, 胡李华, 冯夏庭, 等. 低频周期扰动荷载与静载联合作用下岩爆过程的真三轴试验研究 [J]. 岩石力学与工程学报, 2016, 35(7): 1309–1322. DOI: 10.13722/j.cnki.jrme.2015.1249.SU Guoshao, HU Lihua, FENG Xiating, et al. True triaxial experimental study of rockburst process under low frequency cyclic disturbance load combined with static load [J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(7): 1309–1322. DOI: 10.13722/j.cnki.jrme.2015.1249. [4] 朱万成, 左宇军, 尚世明, 等. 动态扰动触发深部巷道发生失稳破裂的数值模拟 [J]. 岩石力学与工程学报, 2007, 26(5): 915–921. DOI: 10.3321/j.issn:1000-6915.2007.05.007.ZHU Wancheng, ZUO Yujun, SHANG Shiming, et al. Numerical simulation of instable failure of deep rock tunnel triggered by dynamic disturbance [J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(5): 915–921. DOI: 10.3321/j.issn:1000-6915.2007.05.007. [5] 宫凤强, 李夕兵, 刘希灵. 三维动静组合加载下岩石力学特性试验初探 [J]. 岩石力学与工程学报, 2011, 30(6): 1179–1190.GONG Fengqiang, LI Xibing, LIU Xiling. Preliminary experimental study of characteristics of rock subjected to 3D coupled static and dynamic loads [J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(6): 1179–1190. [6] 苏国韶, 陈智勇, 蒋剑青, 等. 不同加载速率下岩爆碎块耗能特征试验研究 [J]. 岩土工程学报, 2016, 38(8): 1481–1489. DOI: 10.11779/CJGE201608016.SU Guoshao, CHEN Zhiyong, JIANG Jianqing, et al. Experimental study on energy dissipating characteristics of rockburst fragments under different loading rates [J]. Chinese Journal of Geotechnical Engineering, 2016, 38(8): 1481–1489. DOI: 10.11779/CJGE201608016. [7] SU Guoshao, JIANG Jianqing, ZHAI Shaobin, et al. Influence of tunnel axis stress on strainburst: an experimental study [J]. Rock Mechanics and Rock Engineering, 2017, 50(6): 1551–1567. DOI: 10.1007/s00603-017-1181-7. [8] SU Guoshao, ZHAI Shaobin, JIANG Jianqing, et al. Influence of radial stress gradient on strainbursts: an experimental study [J]. Rock Mechanics and Rock Engineering, 2017, 50(10): 2659–2676. DOI: 10.1007/s00603-017-1266-3. [9] 张艳博, 徐东强. 岩爆在不同岩石中的模拟实验 [J]. 河北理工学院学报, 2002, 24(4): 8–11. DOI: 10.3969/j.issn.1674-0262.2002.04.002.ZHANG Yanbo, XU Dongqiang. The analysis of the experiment of rock burst on different rocks [J]. Journal of Hebei Institute of Technology, 2002, 24(4): 8–11. DOI: 10.3969/j.issn.1674-0262.2002.04.002. [10] 赵康, 赵红宇, 贾群燕. 岩爆岩石断裂的微观结构形貌分析及岩爆机理 [J]. 爆炸与冲击, 2015, 35(6): 913–918. DOI: 10.11883/1001-1455(2015)06-0913-06.ZHAO Kang, ZHAO Hongyu, JIA Qunyan. An analysis of rockburst fracture micromorphology and study of its mechanism [J]. Explosion and Shock Waves, 2015, 35(6): 913–918. DOI: 10.11883/1001-1455(2015)06-0913-06. [11] 楼沩涛. 干燥和水饱和花岗岩的动态断裂特性 [J]. 爆炸与冲击, 1994, 14(3): 249–254.LOU Weitao. Dynamic fracture behaviour of dry and waterlogged granites [J]. Explosion and Shock Waves, 1994, 14(3): 249–254. [12] 冯夏庭, 陈炳瑞, 张传庆. 岩爆孕育过程的机制、预警与动态调控[M]. 北京: 科学出版社, 2013: 159–167. [13] 赵红亮, 周又和. 深埋地下洞室断裂型岩爆机理的数值模拟 [J]. 爆炸与冲击, 2015, 35(3): 343–349. DOI: 10.11883/1001-1455-(2015)03-0343-07.ZHAO Hongliang, ZHOU Youhe. Numerical simulation on mechanism of fractured rock burst in deep underground tunnels [J]. Explosion and Shock Waves, 2015, 35(3): 343–349. DOI: 10.11883/1001-1455-(2015)03-0343-07. [14] 周辉, 孟凡震, 张传庆, 等. 深埋硬岩隧洞岩爆的结构面作用机制分析 [J]. 岩石力学与工程学报, 2015, 34(4): 720–727. DOI: 10.13722/j.cnki.jrme.2015.04.008.ZHOU Hui, MENG Fanzhen, ZHANG Chuanqing, et al. Effect of structural plane on rockburst in deep hard rock tunnels [J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(4): 720–727. DOI: 10.13722/j.cnki.jrme.2015.04.008. [15] LAN H X, MARTIN C D, HU B. Effect of heterogeneity of brittle rock on micromechanical extensile behavior during compression loading [J]. Journal of Geophysical Research Atmospheres, 2010, 115(B1): B01202. DOI: 10.1029/2009jb006496. [16] WAWERSIK W R, FAIRHURST C. A study of brittle rock fracture in laboratory compression experiments [J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1970, 7(5): 561–575. DOI: 10.1016/0148-9062(70)90007-0. [17] TANG C A, Kaiser P K. Numerical simulation of cumulative damage and seismic energy release during brittle rock failure: Part I: fundamentals [J]. International Journal of Rock Mechanics and Mining Sciences, 1998, 35(2): 113–121. DOI: 10.1016/S0148-9062(97)00009-0. [18] 张翀, 舒赣平. 颗粒形状对颗粒流模拟双轴压缩试验的影响研究 [J]. 岩土工程学报, 2009, 31(8): 1281–1286. DOI: 10.3321/j.issn:1000-4548.2009.08.020.ZHANG Chong, SHU Ganping. Effect of particle shape on biaxial tests simulated by particle flow code [J]. Chinese Journal of Geotechnical Engineering, 2009, 31(8): 1281–1286. DOI: 10.3321/j.issn:1000-4548.2009.08.020. [19] 黄润秋, 王贤能. 岩石结构特征对岩爆的影响研究 [J]. 地质灾害与环境保护, 1997, 7(2): 15–20.HUANG Runqiu, WANG Xianneng. Influence of the structure character of rock on rockburst intensity [J]. Journal of Geological Hazards and Environment Preservation, 1997, 7(2): 15–20. [20] 高秉璋, 洪大卫, 方宗斌. 花岗岩类区1∶5万区域地质填图方法指南[M]. 武汉: 中国地质大学出版社, 1991: 16–17. [21] 陈炳瑞, 冯夏庭, 明华军, 等. 深埋隧洞岩爆孕育规律与机制: 时滞型岩爆 [J]. 岩石力学与工程学报, 2012, 31(3): 561–569. DOI: 10.3969/j.issn.1000-6915.2012.03.014.CHEN Bingrui, FENG Xiating, MING Huajun, et al. Evolution law and mechanism of rockburst in deep tunnel: time delayed rockburst [J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(3): 561–569. DOI: 10.3969/j.issn.1000-6915.2012.03.014. [22] HU Xiaochuan, SU Guoshao, CHEN Guanyan, et al. Experiment on rockburst process of borehole and its acoustic emission characteristics [J]. Rock Mechanics and Rock Engineering, 2019, 52(3): 783–802. DOI: 10.1007/s00603-018-1613-z. [23] CAI M, KAISER P K, MORIOKA H, et al. FLAC/PFC coupled numerical simulation of AE in large-scale underground excavations [J]. International Journal of Rock Mechanics and Mining Sciences, 2007, 44(4): 550–564. DOI: 10.1016/j.ijrmms.2006.09.013. [24] CHANG S H, LEE C I. Estimation of cracking and damage mechanisms in rock under triaxial compression by moment tensor analysis of acoustic emission [J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(7): 1069–1086. DOI: 10.1016/j.ijrmms.2004.04.006. [25] 许金余, 刘石. 大理岩冲击加载试验碎块的分形特征分析 [J]. 岩土力学, 2012, 33(11): 3225–3229. DOI: 10.16285/j.rsm.2012.11.005.XU Jinyu, LIU Shi. Research on fractal characteristics of marble fragments subjected to impact loading [J]. Rock and Soil Mechanics, 2012, 33(11): 3225–3229. DOI: 10.16285/j.rsm.2012.11.005. [26] PENG J, WONG L N Y, TEH C I. Influence of grain size heterogeneity on strength and microcracking behavior of crystalline rocks [J]. Journal of Geophysical Research: Solid Earth, 2017, 122(2): 1054–1073. DOI: 10.1002/2016JB013469. [27] DIEDERICHS M S, KAISER P K, EBERHARDT E. Damage initiation and propagation in hard rock during tunnelling and the influence of near-face stress rotation [J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(5): 785–812. DOI: 10.1016/j.ijrmms.2004.02.003. -