DOI： 10.11817/j.issn.1672-7207.2019.08.015

高频冲击中心孔岩石损伤模型及试验研究

吴万荣,姚壵,刘智

（中南大学 机电学院，湖南 长沙，410083）

摘要:为了提高矿山开采钻孔效率，基于冲击诱导复合破岩原理，建立高频冲击中心孔岩石损伤模型，利用霍普金森试验装置对冲击中心孔后花岗岩试样进行研究，分析高频冲击中心孔后围岩损伤特性。研究结果表明：在同一冲击载荷和一定中心孔半径范围内，冲击产生的围岩损伤半径比(围岩损伤区域半径与中心孔半径之比)随中心孔半径增大而减小；对于某一确定的冲击载荷，可得到相应的损伤半径比，进而由最终钻孔直径确定中心孔半径；在高频冲击中心孔条件下，本文建立的损伤模型结果与试验结果有较好的一致性，模型能够反映高频冲击中心孔岩石的损伤特性。

关键词:高频冲击；中心孔；分离式霍普金森压杆；损伤模型；岩石力学

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

Experimental study on a damage model of rock with center hole under high frequency impact

WU　Wanrong, YAO　Zhuang, LIU　Zhi

(School of Mechanical and Electrical Engineering, Central South University, Changsha 410083, China)

Abstract: In order to improve the drilling efficiency, a damage model of rock with center hole under high-frequency impact was established based on the principle of composite rock breaking induced by shock disturbance. The granite specimen behind the center hole was studied using split Hopkinson pressure bar (SHPB), and the damage characteristics of surrounding rock of the center hole under high-frequency impact were analyzed. The results show that under the same impact load and within a certain radius range, the ratio of the radius of the surrounding rock damage area to the central hole radius (damage radius ratio) decreases with the increase of the central aperture. For a certain impact load, the corresponding damage radius ratio can be obtained, and the radius of the center hole is determined by the final bore diameter. Under high-frequency impact conditions, results from the established damage model is in good consistency with the experimental results from the rock with center hole. The proposed model can reasonably reflect the damage characteristics of the rock with center hole under high-frequency impact.

Key words: high-frequency impact; center hole; split Hopkinson pressure bar (SHPB); damage model; rock mechanics

在大型硬岩矿山开采中，常采用冲击钻孔或牙轮碾压切削钻孔。然而这些钻孔工艺方法在进行硬岩大直径深孔作业时，存在效率较低、钻齿磨损严重等问题。针对这些问题，娄磊等^[1]提出了一种冲击诱导切削复合破岩的钻孔工艺。采用冲击钻头高频冲击钻凿中心孔，围岩会在冲击扰动载荷的影响下产生裂纹损伤，削弱围岩强度；采用牙轮钻头对损伤区域进行切削，从而提高钻孔的效率，降低钻孔能耗，减少钻头磨损。复合破岩最终钻孔孔径、钻头切削能耗与效率受到冲击中心孔周围岩石损伤区域及损伤程度的影响，损伤区域范围及损伤程度又受到冲击中心孔半径的影响。为了确定复合破岩钻孔工艺参数，需要研究冲击中心孔半径对围岩损伤的影响；而提高破岩效率、降低破岩能耗也涉及到隧道开挖、河道水坝开挖、基础工程等方面^[2]。MEGLIS等^[3]应用超声层析成像对隧道开挖的围岩损伤进行现场测试，得到围压损伤程度和损伤分布规律。LU等^[4]在考虑轴向原岩应力影响的条件下，对圆形隧道开挖进行了弹塑性分析。朱晶晶等^[5-6]通过改进的霍普金森压杆(SHPB)装置，对岩石在不同冲击载荷下的力学及损伤特性进行了研究。林大能等^[7]通过循环冲击试验分析了损伤与围压、冲击载荷以及冲击次数之间的关系^[7]。陈乐求等^[8]研究了初始损伤对岩石抗压力学性质的影响。周辉等^[9-10]基于岩石统一能量屈服准则和最小耗能原理建立了损伤本构模型。马念杰等^[11-12]分别对圆形巷道围岩塑性区分布规律以及一般形态进行了研究。由于复合破岩适用于硬质矿山开采，因此，本文作者以质地坚硬、分布广泛、常用于试验研究的花岗岩为研究对象，采用电锤模拟冲击钻头对试样进行高频冲击，再采用SHPB装置对高频冲击中心孔后的试样进行试验研究。在此基础上，基于Weibull分布建立适用于高频冲击中心孔试样的损伤本构模型，对不同冲击中心孔半径的损伤程度及损伤区域分布规律进行研究，以期为复合破岩钻具的结构设计以及高效、节能地进行硬岩大直径深孔采矿工程作业提供参考。

1 高频冲击中心孔损伤区域分析

冲击钻头在轴向力的作用下，不断向前推进，经过一段时间后，岩石在冲击钻头作用下，便逐渐卸荷为0。钻头冲击岩石的应力状态可近似看作布西涅斯克问题^[13]，在冲击载荷P的作用下，冲击中心孔壁(r=r₀，其中r₀为中心孔半径，r为径向半径)处的径向应力为

(1)

式中：μ为泊松比；为中心孔壁处的径向应力；r₀为中心孔半径；r为径向半径。

1.1　弹塑性区应力分析

塑性区(r₀1)满足岩石塑性条件Mohr- Coulomb塑性准则^[14-15]：

(2)

式中：c和φ分别为内聚力和内摩擦角；σ_rp与σ_θp分别为塑性区径向与切向应力；r₁为塑性区半径。在围岩塑性区取一微分体abcd进行受力分析，如图1所示。图中，σ_r为微分体径向应力，σ_θ为微分体切向应力，τ_θr和τ_rθ分别为沿径向和切向的剪应力。

图1　围岩塑性区微元应力分析

Fig. 1　Stress analysis of micro-element in plastic zone of surrounding rock

根据单元体径向和切向受力平衡可得

(3)

式(3)为轴对称条件下的应力平衡方程。联立式(2)和(3)，根据边界条件，可得塑性区应力为：

(4)

(5)

对于弹性区，假设原岩应力为σ₀，弹塑性区交界处的径向应力为σ_r1，则围岩弹性区(r₁≤r)的应力^[16]为

(6)

式中：σ_re与σ_θe分别为弹性区径向与切向应力。

1.2　冲击损伤半径及变形分析

弹性区应力和与塑性区应力之和分别为：

(7)

(8)

弹塑性交界处径向应力和与切向应力之和相等^[17]，则有

(9)

由式(9)即可求出塑性区半径r₁。

总应变包含可恢复的弹性应变与不可恢复的应变，根据弹塑性理论，围岩塑性区总应变为

(10)

式中：E和G分别为弹性模量和剪切模量；u_r和u_θ分别为塑性区径向与切向变形；f为塑性函数，弹性变形阶段f=0。由式(10)可知体应变ε_v为

(11)

对式(11)积分，可得塑性区总变形Δ为

(12)

则由损伤面积定义的钻头冲击产生的围岩损伤D₀为

(13)

式中：R为试样半径。

2 高频冲击中心孔损伤本构模型

2.1　损伤本构模型的建立

损伤体服从Weibull损伤分布^[18-19]规律，概率密度函数p(F)的表达式为

(14)

式中：F为微元体破坏强度的分布变量；m和F₀为 Weibull分布参数，反映岩石材料的力学性质。统计损伤变量D为

(15)

式中：为总微元数。联立式(14)和(15)得

(16)

由式(16)可知：损伤变量与岩石微元体强度有关。考虑到Drucker-Prager破坏准则具有参数形式简单、在岩石材料中应用广泛等优点，基于Drucker-Prager破坏准则的岩石微元强度为

(17)

式中：；为应力张量的第一不变量；为应力偏量的第二不变量；k为与材料内聚力、内摩擦角相关的常数。结合胡克定律可得：

(18)

(19)

式中：和分别为最大和最小主应力，为主应变。

在没有围压的情况即一维应力状态下，有，(为试样沿轴向加压方向的应变)。此时，联立式(17) ~ (19)得

(20)

则高频冲击中心孔岩石的损伤本构方程为

(21)

2.2　损伤本构模型参数计算

对应力-应变曲线有

(22)

由式(21)和式(22)可得：

(23)

式中：；和分别为试验应力-应变曲线上的峰值应力和峰值应变。由式(23)求出损伤模型的参数m和F₀，进而确定损伤本构模型。

3 试验结果分析

3.1　花岗岩试样制备及力学参数

从开采自同一矿场的花岗岩中选用完整性和均质性较好的作为试样。根据岩石力学试验性能测试要求，试样直径×高为50 mm×50 mm。对每个试样两端进行打磨，保证平行度与垂直度小于0.02 mm。然后，用实测冲击力为3.49 kN、冲击频率为50 Hz、转速为1 000 r/min的电锤对花岗岩旋转冲击，形成不同直径的中心孔。花岗岩力学参数见表1。

表1　花岗岩力学参数

Table 1　Mechanical parameters of granite

3.2　试验设备及方案

试验采用大直径SHPB装置，入射杆与试样等截面积，以冲击的方式进行恒应变率半正弦波加载^[20-21]，采用DL-850E示波器及CS-1D超动态应变仪采集试验数据。试样冲击时无轴压和围压，对试样进行预冲击，调整至合适的冲击气压，用同种冲击气压冲击试样1~2次，测得每次的应力-应变曲线，用声波检测仪测得每次冲击前、后的纵波波速。试样具体力学参数及试验结果如表2所示。

表2　花岗岩试样力学参数与试验结果

Table 2　Mechanical parameters and experimental results of granite specimens

3.3　试验结果及分析

图2所示为无冲击作用时花岗岩试样的应力-应变曲线。由图2可以看出：铣削加工不同中心孔半径试样的应力-应变曲线与完整标准试样曲线的峰值强度与形状基本一致，而由表2可知：二者的纵波波速也基本相同。由此可见，在一定范围内，中心孔仅由旋转切削加工，无冲击作用时，不会对围岩产生损伤，且不受孔半径的影响。

图3所示为冲击作用下试样的应力-应变曲线，图4所示为试样峰值应力与中心孔半径的关系。由图3可知：在高频冲击作用下试样峰值应力明显下降，且其纵波波速也减小(见表2)，其中试样2-1相较于试样3-1和试样4-1峰值应力下降更多，且试样3-1和试样2-1之间的峰值应力下降幅度比试样4-1和试样3-1之间的峰值应力下降幅度大，说明试样2-1产生的围岩损伤更大。当冲击载荷不变时，在一定中心孔半径范围内，随着中心孔半径减小，其产生的径向应力越大，造成的围岩损伤也越大，围岩损伤程度增大的幅度也在上升。

图2　无冲击作用时花岗岩试样应力-应变曲线

Fig. 2　Stress-strain curves of granite specimens without impact loads

图5所示为不同中心孔半径试样冲击试验结果与损伤本构模型计算结果对比。由图5可知：本构模型计算结果与试验结果有较好的一致性，说明建立的高频冲击中心孔损伤模型能够较好地反映带有围岩损伤的岩样强度与应变的关系，验证了模型的正确性。由图5还可以看出：3种不同中心孔半径试样在第2次冲击试验后，峰值应力均减小，且中心孔半径小的试样峰后应力-应变曲线更平缓。这意味着中心孔半径小的试样再次受到相同应力作用时，会产生更大的应变，更容易破坏。

图3　冲击作用下试样应力-应变曲线

Fig. 3　Stress-strain curves of granite specimens under impact loads

图4　试样峰值应力与中心孔半径关系

Fig. 4　Relationship between peak stress and center hole radius of specimens

图5　冲击试样试验结果与模型计算结果对比

Fig. 5　Comparison of experimental and model calculation results of specimens under impact loads

图6　r₁/r₀及冲击围岩损伤D₀与中心孔半径r₀的关系

Fig. 6　Relationship between r₁/r₀ , impact damage D₀ and central hole radius r₀

在同一冲击载荷下，高频冲击中心孔造成的围岩损伤D₀，围岩损伤区域半径与中心孔半径的比值r₁/r₀与中心孔半径的关系如图6所示。其中，图6 (b)中的试验结果根据试验实测的纵波波速与密度通过波阻抗^[22]定义的损伤公式计算得出。由图6可见：试验结果与模型计算结果相符，模型能够反映高频冲击中心孔的围岩损伤程度。从图6还可以看出：在同一冲击载荷作用下，D₀和r₁/r₀都随中心孔半径增大而减小，且减小趋势在不断减缓；当中心孔半径增大到一定程度时，几乎不会对围岩产生损伤。对于半径较小的中心孔，同样的冲击载荷在孔壁处产生的径向应力大，可能会造成中心孔内壁围岩破坏脱落的现象，此时，减小中心孔径，损伤区域与中心孔半径比也不再明显增长。对于某一确定的冲击载荷，可得到相应的r₁/r₀，进而根据最终的钻孔直径确定所需的中心孔半径r₀。

4 结论

1) 高频冲击中心孔会对围岩造成损伤；在同一冲击载荷作用下，随着中心孔半径的减小，岩样峰值应力下降，损伤区域与中心孔半径比r₁/r₀增大，冲击产生的围岩损伤增大；半径比的增长速率随着中心孔半径的减小而增大，当中心孔半径减小到一定值后，r₁/r₀不再明显增长；当冲击载荷不变时，中心孔半径r₀增大到一定值后，半径比为1，不再对围岩产生损伤。

2) 根据冲击载荷与建立的损伤模型得到r₁/r₀，进而由最终的钻孔直径确定所需的中心孔半径r₀，可为复合凿岩钻具参数设计及硬岩大直径深孔采矿提供参考。

3) 建立的损伤模型适应于高频冲击中心孔条件，模型计算结果与试验结果有较好的一致性；模型能够反映冲击中心孔岩样的应力、冲击围岩损伤与应变的关系，验证了模型的合理性。

参考文献：

[1] 娄磊, 吴万荣, 王送来, 等. 冲击诱导损伤岩石的侵入特性[J]. 东北大学学报(自然科学版), 2017, 38(10): 1501-1506.

LOU Lei, WU Wanrong, WANG Songlai, et al. Intrusive characters of damage rock by shock disturbance[J]. Journal of Northeastern University(Natural Science), 2017, 38(10): 1501-1506.

[2] 黄理兴, 陈奕柏. 我国岩石动力学研究状况与发展[J]. 岩石力学与工程学报, 2003, 22(11): 1881-1886.

HUANG Lixing, CHEN Yibai. Rock dynamics in China: past, present and future[J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(11): 1881-1886.

[3] MEGLIS I L, CHOW T, MARTIN C D, et al. Assessing in situ microcrack damage using ultrasonic velocity tomography[J]. International Journal of Rock Mechanics and Mining Sciences, 2005, 42(1): 25-34.

[4] LU Aizhong, XU Guisheng, SUN Feng, et al. Elasto-plastic analysis of a circular tunnel including the effect of the axial in situ stress[J]. International Journal of Rock Mechanics and Mining Sciences, 2010, 47(1): 50-59.

[5] 朱晶晶, 李夕兵, 宫凤强, 等. 单轴循环冲击下岩石的动力学特性及其损伤模型研究[J]. 岩土工程学报, 2013, 35(3): 531-539.

ZHU Jingjing, LI Xibing, GONG Fengqiang, et al. Dynamic characteristics and damage model for rock under uniaxial cyclic impact compressive loads[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(3): 531-539.

[6] 宫凤强, 李夕兵, 刘希灵. 一维动静组合加载下砂岩动力学特性的试验研究[J]. 岩石力学与工程学报, 2010, 29(10): 2076-2085.

GONG Fengqiang, LI Xibing, LIU Xiling. Experimental study of dynamic characteristics of sandstone under one-dimensional coupled static and dynamic loads[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(10): 2076-2085.

[7] 林大能, 陈寿如. 循环冲击荷载作用下岩石损伤规律的试验研究[J]. 岩石力学与工程学报, 2005, 24(22): 4094-4098.

LIN Daneng, CHEN Shouru. Experimental study on damage evolution law of rock under cyclical impact loadings[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(22): 4094-4098.

[8] 陈乐求, 张家生, 陈俊桦, 等. 初始损伤对脆性岩石抗压力学性质的影响[J]. 中南大学学报(自然科学版), 2017, 48(2): 484-490.

CHEN Leqiu, ZHANG Jiasheng, CHEN Junhua, et al. Influences of initial damage on mechanics of brittle rock under compressed stress[J]. Journal of Central South University(Science and Technology), 2017, 48(2): 484-490.

[9] 周辉, 李震, 朱国金, 等. 基于岩石统一能量屈服准则的硬岩损伤模型[J]. 岩土力学, 2016, 37(3): 609-615, 624.

ZHOU Hui, LI Zhen, ZHU Guojin, et al. A damage model for hard rock based on unified energy yield criterion of rock[J]. Rock and Soil Mechanics, 2016, 37(3): 609-615, 624.

[10] 孙梦成, 徐卫亚, 王苏生, 等. 基于最小耗能原理的岩石损伤本构模型研究[J]. 中南大学学报(自然科学版), 2018, 49(8): 2067-2075.

SUN Mengcheng, XU Weiya, WANG Susheng, et al. Study on damage constitutive model of rock based on principle of minimum dissipative energy[J]. Journal of Central South University(Science and Technology), 2018, 49(8): 2067-2075.

[11] 马念杰, 李季, 赵志强. 圆形巷道围岩偏应力场及塑性区分布规律研究[J]. 中国矿业大学学报, 2015, 44(2): 206-213.

MA Nianjie , LI Ji, ZHAO Zhiqiang. Distribution of the deviatoric stress field and plastic zone in circular roadway surrounding rock[J]. Journal of China University of Mining & Technology, 2015, 44(2): 206-213.

[12] 郭晓菲, 马念杰, 赵希栋, 等. 圆形巷道围岩塑性区的一般形态及其判定准则[J]. 煤炭学报, 2016, 41(8): 1871-1877.

GUO Xiaofei, MA Nianjie, ZHAO Xidong, et al. General shapes and criterion for surrounding rock mass plastic zone of round roadway[J]. Journal of China Coal Society, 2016, 41(8): 1871-1877.

[13] 刘宗平. 冲击凿岩工具及其理论基础[M].北京：地质出版社, 1987: 107-110.

LIU Zongping. Percussive drilling tool and its theoretical foundation [M]. Beijing：Geological Press，1987: 107-110.

[14] LISJAK A, GRASSELLI G. A review of discrete modeling techniques for fracturing processes in discontinuous rock masses[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2014, 6(4): 301-314.

[15] LISJAK A, GRASSELLI G, VIETOR T. Continuum-discontinuum analysis of failure mechanisms around unsupported circular excavations in anisotropic clay shales[J]. International Journal of Rock Mechanics and Mining Sciences, 2014, 65: 96-115.

[16] 李民庆. 岩体力学的力学基础[M]. 长沙: 湖南科学技术出版社, 1979: 145-148.

LI Minqing. Mechanical basis of rock mechanics[M]. Changsha: Hunan Science and Technology Press, 1979: 145-148.

[17] ALEJANO L R, RODRGUEZ-DONO A, VEIGA M. Plastic radii and longitudinal deformation profiles of tunnels excavated in strain-softening rock masses[J]. Tunnelling and Underground Space Technology, 2012, 30: 169-182.

[18] 曹文贵, 赵明华, 刘成学. 基于Weibull分布的岩石损伤软化模型及其修正方法研究[J]. 岩石力学与工程学报, 2004, 23(19): 3226-3231.

CAO Wengui, ZHAO Minghua, LIU Chengxue. Study on the model and its modifying method for rock softening and damage based on Weibull random distribution[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(19): 3226-3231.

[19] 赵怡晴, 刘红岩, 吕淑然, 等. 基于宏观和细观缺陷耦合的节理岩体损伤本构模型[J]. 中南大学学报(自然科学版), 2015, 46(4): 1489-1496.

ZHAO Yiqing, LIU Hongyan, LU Shuran, et al. Damage constitutive model of jointed rock mass based on coupling macroscopic and mesoscopic flaws[J]. Journal of Central South University(Science and Technology), 2015, 46(4): 1489-1496.

[20] LI X B, LOK T S, ZHAO J, et al. Oscillation elimination in the Hopkinson bar apparatus and resultant complete dynamic stress-strain curves for rocks[J]. International Journal of Rock Mechanics and Mining Sciences, 2000, 37(7): 1055-1060.

[21] 李夕兵, 宫凤强, ZHAO J, 等. 一维动静组合加载下岩石冲击破坏试验研究[J]. 岩石力学与工程学报, 2010, 29(2): 251-260.

LI Xibing, GONG Fengqiang, ZHAO J, et al. Test study of impact failure of rock subjected to one- dimensional coupled static and dynamic loads[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(2): 251-260

[22] 金解放, 李夕兵, 殷志强, 等. 循环冲击下波阻抗定义岩石损伤变量的研究[J]. 岩土力学, 2011, 32(5): 1385-1393, 1410.

JIN Jiefang, LI Xibing, YIN Zhiqiang, et al. A method for defining rock damage variable by wave impedance under cyclic impact loadings[J]. Rock and Soil Mechanics, 2011, 32(5): 1385-1393, 1410.

(编辑  伍锦花)

收稿日期： 2018 -08 -16; 修回日期： 2018 -10 -16

基金项目(Foundation item)：国家自然科学基金资助项目(51774340) (Project(51774340) supported by the National Natural Science Foundation of China)

通信作者：吴万荣，博士，教授，从事液压凿岩技术与装备研究；E-mail: wwr@csu.edu.cn