Three-dimensional crack propagation behaviors of transparent brittle materials under blasting load
-
摘要: 岩石等脆性材料在爆炸荷载作用下裂纹扩展行为通常难以捕捉。基于爆破损伤理论,利用有机玻璃(PMMA)脆性材料透明特性进行爆破模型试验,借助高速摄影技术和CT扫描系统深入探究爆破荷载作用下动态断裂行为和三维裂纹演化规律,结合三维扫描技术揭示裂纹三维形态和破裂面形貌特征。结果表明:在多段爆破能量持续作用下,脆性材料的裂纹存在多次激发扩展的情况;爆炸冲击波产生的初始裂纹数密度高,呈鱼鳞状形貌从“崎岖”向“微波”过渡,平整度提高。其中,裂纹面高程方差值随距离的增大从0.796降低至0.586;最大高度从3.2 m降低至2.8 mm,降低了12.5%。随着到爆炸中心距离的增大,介质由压剪破坏向张拉破坏转变,裂纹分布区分形维数和模型损伤度降低。Abstract: The crack propagation behavior of brittle materials, such as rock, is often challenging to capture under explosive loading conditions. To address this issue, model experiments were conducted based on the theory of explosive damage, utilizing transparent polymethyl methacrylate (PMMA) as a surrogate material to simulate the fracture response of brittle materials. High-speed photography and CT scanning were employed to investigate the dynamic fracture process and three-dimensional crack evolution under blast loading. In addition, 3D scanning technology was used to reconstruct the morphology of cracks and characterize the fracture surface features. The results indicate that under the sustained action of multi-stage explosive energy, cracks undergo repeated initiation and propagation. Initial cracks induced by shock waves exhibit high density and a “fish scale” pattern, primarily concentrated around the blast hole. In contrast, secondary cracks driven by detonation gases have a lower density and extend outward in “ear-shaped” or “dagger-shaped” forms. As the distance from the explosion center increases, the crack surface morphology transitions from rugged to microwave-like textures, with improved flatness. Specifically, the elevation variance of the fracture surface decreases from 0.796 to 0.586, while the maximum height reduces from 3.2 mm to 2.8 mm, representing a 12.5% reduction. Moreover, the failure mode of the material shifts from compressive-shear to tensile failure with increasing distance from the explosion center. This shift is accompanied by a decline in both the fractal dimension of the crack distribution and the overall damage degree of the model.
-
Key words:
- brittle materials /
- blasting load /
- dynamic fracture /
- three-dimensional crack /
- topographical feature
-
钻爆技术凭借其效率高、成本低及施工灵活等优势,已广泛应用于矿山开采、隧道掘进及城市地下空间开发等工程实践。通过合理布设炮孔,钻爆法能够有效释放炸药能量,实现对岩体的定向破碎。在爆破设计过程中,装药量、孔间距及起爆延期等关键参数的优化配置,是提高爆破作业效率、控制爆区形态及降低环境扰动的关键所在。因此,准确掌握爆炸载荷作用下裂纹扩展的空间分布特征和演化规律,对于科学制定爆破参数具有重要理论意义和工程价值。爆炸诱导裂纹的产生、扩展及贯通过程直接决定了岩体的破碎程度和爆破效果;而由爆破应力波与爆生气体共同作用引发的次生裂纹,则对炮孔周边保留岩体的完整性和稳定性产生重要影响,需在设计和施工过程中加以控制。
爆炸载荷作用下的动态断裂问题是爆破破岩、油气开采及工程防护等领域的核心科学问题。近年来,众多学者围绕脆性材料中的爆生裂纹扩展行为开展了大量研究,取得了一系列重要成果。基于PMMA材料的透明特性,杨仁树等[1-2]和岳中文等[3]利用数字激光动态焦散技术对二维爆破模型进行了系统试验,重点探讨了微差延时条件下裂纹动态强度因子及扩展方向的变化规律。李清等[4]和郭洋等[5]借助LS-DYNA、ABAQUS等数值模拟软件,深入分析了单孔爆破条件下孔底裂纹扩展行为及孔间应力场叠加效应对裂纹贯通的影响。
在复杂加载条件下,静载与爆炸荷载的叠加对裂纹行为产生显著影响。如杨立云等[6]、马泗洲等[7]、岳中文等[8]和葛进进等[9]通过施加围压模拟深部岩体环境,揭示了主应力对爆生裂纹扩展方向和速度的控制机制;刘健等[10]则从装药结构出发,探讨了煤岩介质中反射拉应力对裂纹扩展的主导作用。此外,马军等[11]构建了装药不耦合系数与主裂纹长度、扩展速度之间的量化关系;Daehnke等[12]和Forquin等[13]利用叠氮化铅作为爆破材料,研究了PMMA样本中动态断裂的传播,发现应力波的传播速度迅速高于裂纹的扩展速度,并且大多数断裂是由于爆生气体引起的加压而发生的。
围绕预制缺陷对裂纹行为的影响,已有研究表明裂纹扩展路径受应力波方向和几何缺陷控制。沈世伟等[14]通过人工预制缺陷,探究了不同缺陷角度情况下爆生裂纹的扩展规律,发现预制裂纹对爆生裂纹的扩展具有导向作用,同时,裂纹扩展前期的动力主要是压缩波和剪切波作用时间差而形成的应力差,扩展后期的动力主要是爆生气体的准静压作用,如图1所示。王雁冰等[15-16]重点研究了预制水平和竖直缺陷对爆生裂纹扩展的影响,发现预制裂纹的发育依赖爆炸应力波的作用,爆生主、次裂纹扩展的加速度和裂纹尖端的动强度因子呈现振荡变化趋势。与此同时,郭洋等[17]利用数值模拟技术和动态焦散实验技术深入探究了与柱状药包呈0°和90°角预制缺陷下爆生裂纹的行为特征,提出了预制裂纹近、远端翼裂纹的扩展规律。Jeong等[18]利用PMMA圆柱体试件,结合数值模拟技术探究了脆性材料在爆炸荷载作用下的动态断裂行为,发现爆生气体驱动裂纹发育比冲击波压裂更具优势,此外,裂纹扩展的程度在很大程度上取决于炸药装药量。
尽管学者们对二维平面裂纹的扩展规律已进行了较深入的探讨,但真实爆破过程中裂纹扩展具有显著三维特征。爆炸能量释放剧烈、过程极为短暂,孔壁附近瞬时高压难以直接捕捉,且裂纹扩展是典型的三维非稳态行为,传统方法难以观测其时空演化特征。因此,目前针对三维爆生裂纹扩展机制的研究仍相对薄弱,缺乏系统的试验数据和定量的描述手段。
众多学者在爆炸裂纹扩展理论与形态等方面已做了大量的研究,并取得了相应的成果。但大多从二维层面进行裂纹分析,忽略了爆炸作用的三维事件特征;同时,传统物理模型试验并不能将掏槽爆破裂纹的扩展过程可视化。因此,本文中将PMMA作为模型材料,基于其高透光性且在爆炸荷载作用下裂纹扩展形态与岩石类似的特性,探究掏槽爆破三维裂纹形态特征及扩展规律。使用重氮二硝基苯酚(DDNP)炸药为试验提供爆破动载能量,结合三维扫描点云设备、CT扫描系统将爆后不同区域裂纹面三维形貌展现出来,同时对爆后掏槽腔体形态进行三维重构,分析多段起爆下透明非岩石材料三维裂纹分布规律以及形态特征。探讨各区域裂纹形态与应力场耦合关系,揭示冲击波与爆生气体主导作用机制的空间演化规律。其中,高速摄影用于获取爆生裂纹萌生及扩展的时序演化图像,捕捉不同时间节点的裂纹分布特征;点云重建技术则用于后期对试件表面及内部裂纹形貌的三维建模,进一步提取各区域裂纹的几何特征参数,如扩展长度、数量密度和空间分布形态等。在此基础上,引入分形维数、偏度及峰度等统计参数,对不同区域裂纹面形貌特征进行定量分析,以期为爆炸载荷作用下脆性材料三维裂纹扩展的识别、分类和定量表征提供实验依据。
1. 试 验
模型尺寸为20 cm×20 cm×20 cm,如图2(a)所示。试验过程中利用高速摄影仪捕捉爆破过程中裂纹产生、扩展过程,相机设置在离模型2 m处,以10 000 s−1的采样频率和682×764的分辨率进行捕捉。同时,利用时序控制系统精确实现各炮孔间的延期时间,整个试验系统如图2(b)所示。炮孔布置于模型顶部中心位置,炮孔布置如图2(c)所示,所有炮孔分3段起爆,段间延期时间为10 ms。炮孔和空孔深8 cm,炮孔直径为4 mm,空孔直径为8 mm,DDNP作为试验用炸药,采用连续耦合装药,单孔装药量为70 mg,装药高度为7 mm,填塞高度为20 mm,采用正向起爆方式,装药结构如图2(c)所示,实际装药效果如图2(d)所示。
2. 材料的失效
2.1 裂纹演化过程
试验过程中利用固定装置控制各孔炸药在同一水平,以保证爆破能量集中释放。同时,将同段起爆炮孔串联至同一起爆通道,不同段炮孔并联至不同起爆通道。当爆炸荷载作用至炮孔孔壁时,炮孔周围径向裂纹萌生并向周围扩展,扩展过程如图3所示。随着首爆孔的起爆,爆炸能量迅速充斥整个炮孔,在0.2 ms时炮孔周围开始形成径向裂纹。并在爆生气体的进一步驱动下,部分裂纹扩展至空孔并形成贯穿裂纹面,初始槽腔形成。随后,在炮孔底部形成轴向扇面状裂纹,随着深度的增加,夹制作用增大,导致轴向裂纹面扩展范围明显小于径向裂纹。在20.8 ms时,裂纹扩展的能量迅速释放,裂纹扩展停止。
首爆孔爆破时,会在炮孔周围产生密集的径向裂纹。同时,炮孔底部会向下形成一个锥形的环状三维裂纹,如图3(c)所示。在10.0~10.8 ms时段内,第2段爆孔起爆释放能量继续作用在裂纹上,应力波和爆生气体的联合作用导致炮孔周围径向裂纹的迅速扩展,贯穿形成破裂块体。其中,径向裂纹扩展范围明显小于轴向裂纹,径向裂纹扩展速度也相对较小。由于首爆孔在模型上部形成贯穿裂纹,第2段炮孔起爆后,部分能量通过贯穿裂纹溢出,导致炮孔上部径向裂纹扩展范围、速度较小。而首爆孔起爆后仅在炮孔底部形成环形裂纹,并无出现贯穿裂纹,当第2段炮孔起爆后,爆炸作用在裂纹尖端,进而导致轴向裂纹迅速扩展并贯穿至试件边界。但炮孔底部并未呈现出炸药直接作用下的初始裂纹,仅出现爆生气体驱动产生的裂纹,说明在此装药条件下,爆炸产生的高压气体是炮孔底部区域损伤破坏的主要原因。
第3段炮孔起爆后,爆炸能量继续推进径向裂纹扩展,并形成倒三角状裂纹群,裂纹贯穿整个掏槽区域并形成最终的爆破槽腔。在20.4 ms时,由于轴向破裂面已经形成,爆炸能量迅速从破裂面处外泄,导致炮孔周围垂直于破裂面方向仅有少量裂纹产生。
2.2 裂纹的三维形态
多段起爆下爆炸能量持续作用在介质上,裂纹存在着多次被激发、扩展的情况[19-20]。如图4(a)所示,在多段起爆条件下,初始裂纹迅速在炮孔壁周围产生,爆生气体驱动裂纹继续扩展,形成圆弧状或耳状裂纹。由于模型边界效应等因素影响,轴向主裂纹贯穿整个模型并形成破裂面。对裂纹进行编号如图4(b)所示,在第2段炮孔爆破作用下,1号裂纹和2号裂纹连结并扩展至整个模型边界,形成最大的破裂面。同时在1号裂纹起始区域形成衍生裂纹3号和4号,衍生裂纹长度明显小于1号主裂纹长度。由于前序炮孔爆破在模型内部已经产生大量的贯穿式裂纹,第3段炮孔起爆后,部分爆炸能量从裂纹处逃逸,导致炮孔周围裂纹扩展较短,如5号、6号裂纹。
炸药爆炸引起的断裂过程可分为2个阶段:炸药爆炸瞬间,以冲击波形式向周围释放大量能量,作用于介质,在炮孔周围形成放射鳞片状的初始裂纹。这类裂纹主要集中于炮孔附近,裂纹长度较短,但单位区域内裂纹数密度较高,呈现密集分布的网状形态。在本次试验中,经过图像测量得知,炮孔周围初始裂纹长度主要集中在3~8 mm之间,如图4(c)~(d)所示。(2)爆生气体驱动阶段,在高压气体持续加载下,初始裂纹得以扩展,同时诱导远场产生新的裂纹。
此阶段主要形成2类典型裂纹形态:一类为整体形态呈圆弧状的裂纹,从炮孔位置沿径向向外扩展,形似耳朵,裂纹长度较长,边界相对光滑。此类裂纹多形成于爆炸中心周围应力波多向叠加区域,反映出强烈的局部张应力控制作用,常见于第2段和第3段爆炸耦合作用显著的区域。另一类为形似匕首的裂纹,通常起始于初始裂纹,在爆生气体定向驱动下向远场延伸,具有明显的裂纹尖端。该类裂纹路径较长,扩展方向性强,端部呈半圆状或楔形尖端,单位区域密度较低。匕首形裂纹是爆炸裂纹贯穿发展的主要路径,对介质的结构性破坏具有决定性作用,如图4(d)所示。
结合3段起爆的时序特征,第2段与第3段爆炸所激发的应力波在时间和空间上与第1段爆炸形成的应力场发生耦合与干涉,这种多源应力波叠加作用直接影响局部应力集中区域的形成,从而控制了耳朵形与匕首形裂纹的空间分布及其三维形貌演化。
2.3 CT扫描和三维裂纹重建
利用CT扫描和三维重建技术,可以将试件内部爆后三维裂隙网络进行可视化、数字化处理,从而可以直观分析爆后三维裂纹的破坏形态和破坏方式,并方便对裂纹进行定量化处理,试验采用X射线源CT扫描系统,如图5所示。扫描厚度为0.2 mm,图像尺寸为20 cm×20 cm,从xOy、xOz和yOz三面对模型进行扫描,如图6(a)所示,每个方向得到
1000 张原始CT截面切片,如图6(b)所示。同时,利用Avizo软件对爆炸后的PMMA进行二值化处理,并得到了内部裂隙的三维重建模型,如图6(c)所示。从图6(c)可以看到爆后试件内部爆炸腔体的三维裂纹形态,观察到2条裂纹从装药区域向外延伸至试件边界,将试件一分为二。同时,在爆区的另一侧裂纹向试件边界延伸,如图6(d)所示。根据高速影像记录得知,炮孔底部首先形成匕首状裂纹,随即相互贯穿形成断裂面向模型边界延伸。从试件内部裂纹结构透视图可以看出,试件炮孔周边区域爆炸应力波作用产生细小裂纹并相互贯穿。由于扫描精度的原因,并不能将试件内部细小裂纹完全展现,利用AVIZO数据处理软件提取不同切片孔隙度,得到不同切片位置对应体积分数,如图6(e)所示。体积分数为1说明在此切片位置裂隙贯穿至整个模型,此时切片孔隙度为1。从图6(e)可以看出,装药区域附近孔隙度达到1左右,非装药区平均孔隙度约为0.97,说明装药区域裂纹完全贯穿至整个试件,而非装药区裂纹扩展过程中出现部分未完全贯通现象。结合图6(d)可以看出,装药区域由于爆炸冲击波直接作用,形成高密度细小裂纹相互贯穿,同时空孔的存在,加剧孔间裂纹的相互贯穿并最终形成爆破腔体。而爆生气体作用驱动形成的裂纹扩展范围广,是模型产生块体的主要原因。
3. 断裂面的三维形貌
断裂面的三维形貌是指介质在发生断裂时,断裂面的形状、特征以及表面结构[18, 21-22]。通过对断裂面的三维形貌进行分析,能够揭示裂纹孕育机制,探究不同区域内介质的断裂模式及爆破荷载作用下裂纹演化特征。
3.1 激光扫描和点云数据处理
爆后模型使用FreeScan combo计量级手持三维扫描仪对破裂面进行扫描,收集不同区域的裂纹面点云数据并进行数字化处理,最终生成裂纹面的三维形貌色谱图,具体过程见图7。三维扫描系统由扫描主机、激光镜头、数据传输线、数据处理软件以及高性能笔记本组成。设置扫描点距精度为0.1 mm,扫描频率为2 250 000 s−1,扫描面幅为500 mm×500 mm。为了方便处理扫描文件,将扫描获取的点云文件导入Geomagic studio软件进行降噪和补洞填充处理,并进行封装。将封装文件导入Matlab软件,利用Height算法提取并计算曲面高度,使用Colormap语句将曲面不同高度赋予色谱,形成裂纹面三维形貌色谱。
3.2 损伤破裂面形貌变化
由于试件不同区域的破坏模式和破裂角度不同,裂纹面三维形貌差异较大,因此区域选取必须能够反映整个裂纹面形态特征,同时尽量保证各区域大小一致,如图8(a)所示,其中区域1和区域2为爆炸近区,区域3和区域4为爆炸中远区,区域5为爆炸远区。图8(b)~(f)中不同颜色反映爆生裂纹表面平整度,颜色反差越大,裂纹面高度差越大,裂纹表面起伏越明显。从图8(b)~(c)可以看出爆炸近区(区域1和区域2)的裂纹面起伏较明显。统计并分析区域内点云相对高度值得到,区域1中裂纹面相对高差为3.2 mm,区域内高度在2 mm以上的面积占比71.6%;区域2中裂纹面相对高度差为5.1 mm,区域内高度在3.5 mm以上的面积占比14.3% 。虽然点云相对高度占比差异较大,但在区域内有完整的起伏波动,裂纹面均呈现较深下凹和明显的尖锐凸起,整体表现为显著的崎岖状。在爆炸中远区(区域3和区域4),裂纹面起伏的波长明显增大,逐渐形成阶梯状,区域3中裂纹面相对高差为4.75 mm,区域内高度在3.5 mm以上的面积占比8.2%,绝大部分裂纹面处于较低位置,且裂纹面波动较小;区域4中裂纹面相对高差为3.5 mm,且区域内高度在2.0 mm以上的面积占比46.2%,裂纹面并无出现明显完整波动周期。随着三维裂纹尖端继续扩展,裂纹面逐渐平缓并无明显波浪状起伏,同时,裂纹面上呈现出一些细微波纹状径向裂纹,此时,区域5中裂纹面相对高差为2.8 mm,且区域内高度在1.5 mm以上的面积占比41.8%。
爆炸初始阶段,冲击波作用在周围介质上时,爆炸冲击波强度远高于介质的动态抗压强度,炮孔周围产生压缩破坏。随着距离的增大,冲击波衰减为强度较低的应力波,并持续形成压缩区,导致介质在切向上形成拉伸变形,介质沿切向破裂在炮孔附近形成细密的径向裂纹群。裂纹面形貌起伏较明显,表面较粗糙。随后,在应力波和爆生气体共同作用下初始裂纹持续扩展并贯通。
3.3 裂纹面的高差和平整度
通过提取三维断裂面上微凸体的高度和方差定量描述断面三维形貌变化规律,评价掏槽爆破不同距离处裂纹特征[23-25]。选取断裂表面最低点作为参考点,求出断面上任一点高度:
Hx=Nx−NminNmin=min(N1,N2,N3,⋯,Nn−1,Nn) (1) 式中:Hx为破裂面上任一点的相对高度,Nx为破裂面上任一点绝对高度。
为了表征裂纹面高度的变化规律,断裂面上所有点高度Hi的平均值
ˉH 和断裂面微凸体高程方差σ2分别通过下式计算:¯H=1nn∑i=1Hi (2) σ2=1nn∑i=1(Hi−¯H)2 (3) 通过式(1)~(3)可以得到掏槽爆破裂纹面不同位置表面凹凸的相对高度和方差,如图9所示。随着到炮孔距离的增大,断裂面微凸体的最大高度和方差均呈现先增大后减小的趋势,波动幅度较大。其中在炮孔周围区域1到区域2处,最大高度从3.2 mm增大至5.1 mm,增大了59.4%,方差从0.796 mm2增大至0.817 mm2,增大了2.6%,如图9所示。说明在爆炸近区裂纹面点云数据离散性较大,裂纹面起伏较明显,易出现较大断层。由于爆炸冲击波直接作用在炮孔周边,周围介质裂纹的产生方式主要以压拉破坏为主,断面平整度较差。随着距离的增大,区域2到区域3高度维持在4.75 mm。但方差降低至0.648 mm2,降幅达到20.7%,表明距离由爆炸近区变为中远区过程中,裂纹面的平整度提高。随着距离的增大和时间的增长,爆炸冲击波衰减为应力波,同时爆生气体充满裂纹内部,此时裂纹的扩展由应力波和爆生气体共同驱动,裂纹的产生逐渐由拉伸破坏向剪切破坏转变。区域4的最大高度降低为3.5 mm,与区域3相比降低了26.3%,对应的方差为0.728 mm2。随着距离的增大,应力波逐渐衰减,爆生气体压力降低,驱动裂纹扩展的能力减弱,区域5破裂面最大高度为2.8 mm,方差减少至0.586 mm2,减小了19.5%。随着到爆炸区域距离的增大,裂纹的产生方式逐渐发生变化,裂纹面运移空间受到影响,微凸体相对高度明显减小,裂纹面整理均匀性和平整度逐渐提高。
3.4 裂纹面高度频率直方图分布
为了更细致、全面地分析断口的整体形貌,引入了2个特征参数,即偏度和峰度,分析了断口上特征区域内所有微凸体的高度分布特征。偏度(Hsk)是数据分布偏斜方向和程度的度量,可以用来描述断裂面上点云高度频率分布的不对称性,偏度可以通过下式[26]求得:
Hsk=1nn∑i=1(Hi−¯Hσ)3 (4) 峰度(Hku)是用来表征概率密度分布曲线在平均值处的峰值高低,可以描述断裂面上点云高度分布的平坦状态,峰度可以用下式[26]表示:
Hku=1nn∑i=1(Hi−¯Hσ)4−3 (5) 式中:σ为点云高度标准差。
如图10(a)~(b)所示,在爆破近区域1、2时,破裂面上微凸体高度频率分布的偏度均大于1,表明区域内平均高度以上的面积大于平均高度以下的面积,裂纹面高度分布不均匀性明显。且在区域1中峰度Hku为4.35,区域2中峰度Hku为0.47,表明区域2中数据集中度较高,裂纹面起伏程度高于区域1。当距离增大时,如图10(c)~(d)所示,在爆炸中远区域3、区域4,断裂面的形成逐渐从以拉伸破坏为主向剪切破坏转变,高度频率分布偏度减小,峰度均小于3,裂纹面高度分布均匀性增强,破裂面平整度提高。当距离继续增大时,如图10(e)所示,区域5中偏度增大,峰度继续减小,破裂面上低于高度平均值的凹凸体数量增加,但凹凸体的高度差减小,裂纹面形态变得更加平整。
如图10(a)~(b)所示,在爆炸近区区域1和区域2内,裂纹面微凸体的高度频率分布偏度均大于1,说明高度分布存在明显的右偏性,即表面高于平均高度的区域占比显著,表征裂纹面起伏特征存在较强不均匀性。进一步对比可知,区域1中峰度Hku为4.35,而区域2中峰度Hku仅为0.47。结合图8可知,尽管两者均位于近爆区,但区域2的高度分布更集中,表面微观起伏程度更剧烈,突起构造更密集。
随着测点到炮孔距离的增大,进入区域3和区域4(图10(c)~(d)),偏度呈下降趋势,峰度均小于3,显示裂纹面高度分布趋于对称且波峰尖锐性降低。该变化反映断裂模式由炮孔周围以拉伸主导的高能快速破坏逐步向剪切主导的能量耗散性破坏过渡,裂纹面整体趋于平整。
当进一步远离爆源至区域5(图10(e))时,偏度略有回升而峰度继续降低,说明裂纹面局部下陷构造增多,尽管低于平均高度的凹陷体数量增加,但其整体高度差异变小,表面起伏被进一步抹平,表现为一种局部扰动增强而整体平整性进一步提升。
4. 裂纹分形维数及破坏模式
4.1 裂纹分形维数
计盒维数法提供了一种定量分析复杂几何形状和现象的有效手段,通过盒子覆盖和数量关系,能够揭示物体的内在结构和维度特性[27-29],通过公式(6)计算得到:
lgNδk=−Dflgδk+b (6) 式中:
Nδk 为含有裂纹区域的盒子数目,Df 为区域裂纹场的分形维数,δk 为裂纹区域网格边长,b为拟合参数。将裂纹图案用边长为δk 的立体盒子进行覆盖,并记录覆盖所需盒子的数量Nδk ,同时改变盒子边长,并观察盒子数量的变化。通过线性回归的方式得到裂纹场的斜率Df ,即为该形状的分形维数。试件损伤度
ω 与分形维数的关系表达式[30]为:ω=Df−D0Dmaxf−D0 (7) 式中:D0为材料爆前内部初始损伤面积的分形维数,
Dmaxf 为材料达到的最大损伤面积时的分形维数,对于平面问题,Dmaxf =2。图11为模型不同位置爆后CT扫描切片裂纹的分形维数拟合曲线,从图中可以看出,模型上部分形维数为1.535,根据式(7)得到模型损伤度为0.535;模型中部裂纹分布区分形维数为1.466,模型损伤度为0.466,其中损伤度下降了12.9%。模型底部裂纹分布区分形维数为1.354,对应的模型损伤度为0.354,比模型中部损伤度下降了24.1%。可以看出,z轴方向上,不同高度裂纹分布区分形维数和损伤度有显著差异,但整体呈现出降低趋势,爆生裂纹复杂程度下降。这是由于爆炸冲击波直接作用产生了大量的微裂隙场,并相互贯穿形成损伤,随着爆炸中心距离的增大,爆炸冲击波作用减弱,介质由压剪破坏为主逐渐向张拉破坏转变。
4.2 三维裂纹扩展模式
脆性材料在爆破载荷作用下的破碎过程通常可分为以下3个阶段:(1)在爆炸初期,炸药周围介质首先受到高强度冲击波的作用,产生瞬时高压,导致局部介质发生压缩破碎并形成初始裂隙;(2)随着应力波的传播与反射,介质内部在环向拉应力及反射波诱导的张应力作用下,部分初始裂隙进一步扩展,引发介质的二次破裂过程;(3)随后,爆生气体的高速膨胀进一步推动裂隙贯通,促使岩体解体形成一定粒径的破碎块体,表现出成块或成片剥离的破坏特征。按照介质中裂隙尖端与附加应力间的关系可以将裂隙分为3种主要状态,即张开型裂纹、滑开型裂纹以及撕开型裂纹[31]。
Schöllmann等[32]、Richard等[33]和Wang等[34]在传统裂纹扩展准则基础上进行改进,将裂纹扩展过程中裂纹偏转角φ0和ψ0对裂纹扩展的影响考虑在内,形成新的三维裂纹扩展准则。该准则基于以下2个前提:
(1)裂纹从裂纹尖端径向扩展至垂直于最大主应力σʹ1的方向。σʹ1为裂纹尖端周围单位圆柱面上的最大主应力,其由σφ,σz以及τφz决定,如图12所示,图中P点为特征点。
(2)当裂纹尖端附近最大主应力σʹ1达到介质破坏强度σʹ1C时,裂纹将发生不稳定扩展。
因此,当爆炸冲击波瞬时产生极高的压荷载作用在炮孔周围介质上,其荷载远大于材料临界损伤值,炮孔周围介质在高压环境下产生压剪损伤,主要发生模式Ⅱ破坏,裂纹向着径向方向传播。此模式一般发生在爆炸近区如炮孔壁附近,如图13所示,图中σ为裂纹尖端所受应力。随着后续炮孔的起爆和应力波叠加效应,裂纹尖端受到不同方向应力作用,产生模式Ⅲ和复合模式受力情况,介质发生压剪损伤的同时伴随着纵向撕开型损伤,因此裂纹呈不规则形状。随着时间推移,应力波作用减弱,爆生气体均匀充斥在裂纹中,并以气楔膨胀作用驱动裂纹继续以张开方式向外传播,此时裂纹尖端受力如模式Ⅰ,主要发生张拉破坏。
通过观察裂纹扩展形态得到模式Ⅰ、Ⅱ情况下,炮孔一侧沿裂纹扩展方向距离明显小于垂直于裂纹扩展方向上的距离,其裂纹横纵长度比约在0.5以下。而模式Ⅲ和复合模式情况下,沿裂纹扩展方向上距离明显大于垂直裂纹扩展方向上的距离。由于试验过程中装药轴向长度大于径向长度,整体药柱呈柱状。爆炸气体和能量释放主要向炮孔延伸方向传递,因此轴向上的能量释放率较高,应力水平高于径向方向,导致裂纹优先在此方向进行扩展。另外,在炮孔边缘和爆破近区应力集中同样会引起裂纹沿轴向扩展,如图13所示。试验结果说明装药结构以及药柱形态对介质断裂和裂纹扩展有一定影响。现场生产过程中,可通过合理的改变装药结构和药柱形态达到改变爆破效果的目的。
5. 结 论
采用DDNP炸药为试验提供爆破荷载,利用PMMA材料脆性透明特性进行了多段爆破动态断裂损伤试验,探究了多段爆破作用下脆性材料的裂纹扩展特性,借助高速摄影技术捕捉爆破过程中三维裂纹扩展行为,结合扫描点云技术得到不同区域裂纹三维形貌特征,得到的结论如下。
(1)多段起爆下,爆炸冲击波在首爆孔周围形成密集径向初始裂纹,后续爆破荷载驱动裂纹逐步贯通。爆生气体在孔底区域诱发轴向裂纹快速扩展,是非装药区域材料破坏的主导因素。
(2)爆后装药区形成高密度放射状裂纹,非装药区裂纹呈耳状或匕首状沿径向扩展,贯穿性强,主要由爆生气体驱动,是试件破裂的关键破坏形式。
(3)随着距爆心距离的增加,裂纹面三维形貌特征表现为最大高度和方差先升后降,峰度和偏度整体减小,表明远离爆心区域后裂纹面趋于平整,表面微凸结构减弱。
(4)炮孔周围裂纹多以模式Ⅱ启动,受爆生气体驱动后以模式Ⅰ扩展贯通;随着距离的增大,裂纹受应力耦合影响,以模式Ⅲ及复合型形式扩展,材料破坏形态复杂。
-
-
[1] 杨仁树, 丁晨曦, 王雁冰, 等. 爆炸应力波与爆生气体对被爆介质作用效应研究 [J]. 岩石力学与工程学报, 2016, 35(S2): 3501–3506. DOI: 10.13722/j.cnki.jrme.2016.0066.YANG R S, DING C X, WANG Y B, et al. Action-effect study of medium under loading of explosion stress wave and explosion gas [J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(S2): 3501–3506. DOI: 10.13722/j.cnki.jrme.2016.0066. [2] 杨仁树, 丁晨曦, 杨国梁, 等. 微差爆破的爆生裂纹扩展特性试验研究 [J]. 振动与冲击, 2017, 36(24): 97–102. DOI: 10.13465/j.cnki.jvs.2017.24.015.YANG R S, DING C X, YANG G L, et al. Tests for blasting induced crack propagation characteristics of short-delay blasting [J]. Journal of Vibration and Shock, 2017, 36(24): 97–102. DOI: 10.13465/j.cnki.jvs.2017.24.015. [3] 岳中文, 张士春, 邱鹏, 等. 切缝药包微差爆破爆生裂纹扩展机理 [J]. 煤炭学报, 2018, 43(3): 638–645. DOI: 10.13225/j.cnki.jccs.2017.0889.YUE Z W, ZHANG S C, QIU P, et al. Mechanism of explosive crack propagation with slotted cartridge millisecond blasting [J]. Journal of China Coal Society, 2018, 43(3): 638–645. DOI: 10.13225/j.cnki.jccs.2017.0889. [4] 李清, 徐文龙, 郭洋, 等. 柱状炮孔端部爆生裂纹动态扩展力学行为研究 [J]. 岩石力学与工程学报, 2019, 38(2): 267–275. DOI: 10.13722/j.cnki.jrme.2018.1096.LI Q, XU W L, GUO Y, et al. Study on mechanical behaviors of crack dynamic propagation at the end of cylinder blastholes [J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(2): 267–275. DOI: 10.13722/j.cnki.jrme.2018.1096. [5] 郭洋, 李清, 杨仁树, 等. 三维模型柱状药包爆生裂纹扩展规律研究 [J]. 振动与冲击, 2020, 39(10): 133–140,184. DOI: 10.13465/j.cnki.jvs.2020.10.018.GUO Y, LI Q, YANG R S, et al. Study on crack propagation law of cylindrical charges in three-dimensional models [J]. Journal of Vibration and Shock, 2020, 39(10): 133–140,184. DOI: 10.13465/j.cnki.jvs.2020.10.018. [6] 杨立云, 丁晨曦, 郑立双, 等. 初始静态压应力场中爆生裂纹的扩展行为 [J]. 岩土工程学报, 2018, 40(7): 1322–1327. DOI: 10.11779/CJGE201807020.YANG L Y, DING C X, ZHENG L S, et al. Evolution of blasting cracks in different static compression fields [J]. Chinese Journal of Geotechnical Engineering, 2018, 40(7): 1322–1327. DOI: 10.11779/CJGE201807020. [7] 马泗洲, 刘科伟, 杨家彩, 等. 初始应力下岩体爆破损伤特性及破裂机理 [J]. 爆炸与冲击, 2023, 43(10): 105201. DOI: 10.11883/bzycj-2023-0151.MA S Z, LIU K W, YANG J C, et al. Blast-induced damage characteristics and fracture mechanism of rock mass under initial stress [J]. Explosion and Shock Waves, 2023, 43(10): 105201. DOI: 10.11883/bzycj-2023-0151. [8] 岳中文, 田世颖, 张士春, 等. 单向围压作用下切缝药包爆破爆生裂纹扩展规律的研究 [J]. 振动与冲击, 2019, 38(23): 186–195. DOI: 10.13465/j.cnki.jvs.2019.23.027.YUE Z W, TIAN S Y, ZHANG S C, et al. Expanding law of cracks formed by slotted cartridge blast under unidirectional confining pressure [J]. Journal of Vibration and Shock, 2019, 38(23): 186–195. DOI: 10.13465/j.cnki.jvs.2019.23.027. [9] 葛进进, 徐颖, 程琳, 等. 深部岩石爆破主裂纹扩展方向与地应力的关系 [J]. 振动与冲击, 2023, 42(4): 54–64. DOI: 10.13465/j.cnki.jvs.2023.04.007.GE J J, XU Y, CHENG L, et al. Relationship between in-situ stress and propagation direction of main cracks induced by blasting [J]. Journal of Vibration and Shock, 2023, 42(4): 54–64. DOI: 10.13465/j.cnki.jvs.2023.04.007. [10] 刘健, 刘泽功, 高魁, 等. 不同装药模式爆破载荷作用下煤层裂隙扩展特征试验研究 [J]. 岩石力学与工程学报, 2016, 35(4): 735–742. DOI: 10.13722/j.cnki.jrme.2015.0865.LIU J, LIU Z G, GAO K, et al. Experimental study of extension characters of cracks in coal seam under blasting load with different charging modes [J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(4): 735–742. DOI: 10.13722/j.cnki.jrme.2015.0865. [11] 马军, 汪旭光, 李祥龙, 等. 不耦合装药刻痕爆破裂纹的动态力学特征及损伤分形规律实验 [J]. 兵工学报, 2023, 44(12): 3676–3686. DOI: 10.12382/bgxb.2022.1270.MA J, WANG X G, LI X L, et al. Experiment on dynamic mechanical characteristics and damage fractal law of crack in decoupled charge scratch blasting [J]. Acta Armamentarii, 2023, 44(12): 3676–3686. DOI: 10.12382/bgxb.2022.1270. [12] DAEHNKE A, ROSSMANITH H P, NAPIER J A L. Gas pressurisation of blast-induced conical cracks [J]. International Journal of Rock Mechanics and Mining Sciences, 1997, 34(3/4/): 263.e1–263.17. DOI: 10.1016/S1365-1609(97)00282-7. [13] FORQUIN P, NASRAOUI M, RUSINEK A, et al. Experimental study of the confined behaviour of PMMA under quasi-static and dynamic loadings [J]. International Journal of Impact Engineering, 2012, 40/41: 46–57. DOI: 10.1016/j.ijimpeng.2011.09.007. [14] 沈世伟, 李国良, 李冬, 等. 不同角度预制裂隙条件下双孔爆破裂纹扩展规律 [J]. 煤炭学报, 2019, 44(10): 3049–3057. DOI: 10.13225/j.cnki.jccs.2018.1362.SHEN S W, LI G L, LI D, et al. Crack propagation law of two-hole blasting under different angles of prefabricated fissure [J]. Journal of China Coal Society, 2019, 44(10): 3049–3057. DOI: 10.13225/j.cnki.jccs.2018.1362. [15] 王雁冰, 商禹智, 石震鑫, 等. 定向断裂双孔爆破含缺陷介质裂纹扩展的动焦散试验 [J]. 爆破, 2018, 35(1): 15–20,48. DOI: 10.3963/j.issn.1001-487X.2018.01.003.WANG Y B, SHANG Y Z, SHI Z X, et al. Dynamic caustics experiment on crack propagation in defective medium by directional breaking with double hole blasting [J]. Blasting, 2018, 35(1): 15–20,48. DOI: 10.3963/j.issn.1001-487X.2018.01.003. [16] 王雁冰, 杨仁树, 丁晨曦, 等. 双孔爆炸应力波作用下缺陷介质裂纹扩展的动焦散试验 [J]. 煤炭学报, 2016, 41(7): 1755–1761. DOI: 10.13225/j.cnki.jccs.2015.0370.WANG Y B, YANG R S, DING C X, et al. Dynamic caustics experiment on crack propagation of defective medium under the effect of explosive stress waves of double holes [J]. Journal of China Coal Society, 2016, 41(7): 1755–1761. DOI: 10.13225/j.cnki.jccs.2015.0370. [17] 郭洋, 李清, 徐文龙, 等. 条形药包爆破预制贯通裂纹动态断裂过程研究 [J]. 岩土力学, 2018, 39(10): 3882–3890. DOI: 10.16285/j.rsm.2018.0219.GUO Y, LI Q, XU W L, et al. Dynamic fracture process of a pre-crack under linear charge explosion [J]. Rock and Soil Mechanics, 2018, 39(10): 3882–3890. DOI: 10.16285/j.rsm.2018.0219. [18] JEONG H, JEON B, CHOI S, et al. Fracturing behavior around a blasthole in a brittle material under blasting loading [J]. International Journal of Impact Engineering, 2020, 140: 103562. DOI: 10.1016/j.ijimpeng.2020.103562. [19] XIE L X, LU W B, ZHANG Q B, et al. Damage evolution mechanisms of rock in deep tunnels induced by cut blasting [J]. Tunnelling and Underground Space Technology, 2016, 58: 257–270. DOI: 10.1016/j.tust.2016.06.004. [20] PETROV Y V. Structural-temporal approach to modeling of fracture dynamics in brittle media [M]//ZHAO J, LI J C. Rock Dynamics and Applications-State of the Art. London: CRC Press, 2013: 101-110. [21] TAO J, YANG X G, LI H T, et al. Numerical investigation of blast-induced rock fragmentation [J]. Computers and Geotechnics, 2020, 128: 103846. DOI: 10.1016/j.compgeo.2020.103846. [22] BANADAKI M M D, MOHANTY B. Numerical simulation of stress wave induced fractures in rock [J]. International Journal of Impact Engineering, 2012, 40/41: 16-25. DOI: 10.1016/j.ijimpeng.2011.08.010. [23] CHEN B, SHEN B T, ZHANG S C, et al. 3D morphology and formation mechanism of fractures developed by true triaxial stress [J]. International Journal of Mining Science and Technology, 2022, 32(6): 1273–1284. DOI: 10.1016/j.ijmst.2022.09.002. [24] FENG X T, KONG R, ZHANG X W, et al. Experimental study of failure differences in hard rock under true triaxial compression [J]. Rock Mechanics and Rock Engineering, 2019, 52(7): 2109–2122. DOI: 10.1007/s00603-018-1700-1. [25] ZHU S, ZHENG J H, ZHU Z D, et al. Experiments on three-dimensional flaw dynamic evolution of transparent rock-like material under osmotic pressure [J]. Tunnelling and Underground Space Technology, 2022, 128: 104624. DOI: 10.1016/j.tust.2022.104624. [26] LIU Y, HUANG D, CEN D F, et al. Tensile strength and fracture surface morphology of granite under confined direct tension test [J]. Rock Mechanics and Rock Engineering, 2021, 54(9): 4755–4769. DOI: 10.1007/s00603-021-02543-7. [27] 左进京, 杨仁树, 龚敏, 等. 分段装药爆炸应变场与裂隙场分布规律 [J]. 爆炸与冲击, 2023, 43(3): 035101. DOI: 10.11883/bzycj-2022-0333.ZUO J J, YANG R S, GONG M, et al. On the distribution of explosion strain field and fracture field in segment charge [J]. Explosion and Shock Waves, 2023, 43(3): 035101. DOI: 10.11883/bzycj-2022-0333. [28] 周星源, 岳中文, 金庆雨, 等. 反射爆炸应力波作用下动静裂纹的贯通机理[J/OL]. 爆炸与冲击, 1–16. DOI: 10.11883/bzycj-2024-0409.ZHOU X, YUE Z, JIN Q, et al. Experimental study on the influence of the reflected explosiveStress waves on the dynamic crack propagation characteristics[J/OL]. Explosion and Shock Waves, 1–16. DOI: 10.11883/bzycj-2024-0409. [29] 王雁冰, 付代睿, 李杨, 等. 不同耦合介质爆破裂纹动态扩展特性研究 [J]. 力学与实践, 2023, 45(5): 1021–1032. DOI: 10.6052/1000-0879-23-327.WANG Y B, FU D R, LI Y, et al. Study on dynamic propagation characteristics of blasting cracks in different coupling media [J]. Mechanics in Engineering, 2023, 45(5): 1021–1032. DOI: 10.6052/1000-0879-23-327. [30] 李成孝, 张渊通, 安晨. 单侧开半圆孔PMMA试件Ⅰ型和Ⅰ-Ⅱ混合型裂纹动态扩展及数值模拟研究 [J]. 矿业科学学报, 2020, 5(5): 490–501. DOI: 10.19606/j.cnki.jmst.2020.05.003.LI C X, ZHANG Y T, AN C. Study on the dynamic propagation and numerical simulation of mode I and mixed mode I-II cracks in PMMA specimens with unilateral semicircular holes [J]. Journal of Mining Science and Technology, 2020, 5(5): 490–501. DOI: 10.19606/j.cnki.jmst.2020.05.003. [31] 孙强, 李雪东, 姚腾飞, 等. 基于DIC的爆炸加载下脆性材料裂纹扩展规律的试验研究 [J]. 爆炸与冲击, 2019, 39(10): 103102. DOI: 10.11883/bzycj-2018-0308.SUN Q, LI X D, YAO T F, et al. Experimental study on crack propagation of brittle materials based on DIC under explosive loading based on DIC [J]. Explosion and Shock Waves, 2019, 39(10): 103102. DOI: 10.11883/bzycj-2018-0308. [32] SCHÖLLMANN M, RICHARD H A, KULLMER G, et al. A new criterion for the prediction of crack development in multiaxially loaded structures [J]. International Journal of Fracture, 2002, 117(2): 129–141. DOI: 10.1023/A:1020980311611. [33] RICHARD H A, SCHRAMM B, SCHIRMEISEN N H. Cracks on Mixed Mode Loading - Theories, experiments, simulations [J]. International Journal of Fatigue, 2014, 62: 93–103. DOI: 10.1016/j.ijfatigue.2013.06.019. [34] WANG J, REN L, XIE L Z, et al. Maximum mean principal stress criterion for three-dimensional brittle fracture [J]. International Journal of Solids and Structures, 2016, 102/103: 142–154. DOI: 10.1016/j.ijsolstr.2016.10.009. -