Study of the splitting crack propagation morphology using high-speed 3D DIC
-
摘要: 为探寻非钻孔条件下露天爆破大块二次破碎形态的控制方法, 应用线性聚能射流对圆柱混凝土模型试件进行侵彻实验, 使用High-speed 3D DIC(高速三维数字图像相关方法)方法分析试件劈裂发展过程的全场三维形变特征。研究结果表明, 数据分析区内劈裂裂纹扩展速度在4个区间内呈阶梯式变化趋势, 峰值速度为235.52 m/s, 平均速度为140.89 m/s;线性聚能射流侵彻对劈裂裂纹扩展有明显导向作用, 应力集中作用使得劈裂裂纹围绕线性射流侵彻对称轴扩展, 扩展方向变化幅度较小;在劈裂裂纹扩展速度突变的3个时刻, 劈裂裂纹路径产生了3处明显拐点, 在拐点处伴随有支裂纹的产生, 支裂纹的扩展距离均未超过5 cm;主应变集中带形状及分布位置决定了裂纹扩展路径及趋势, 拉应变集中先于裂纹出现, 试件呈现准静态劈裂形态, 劈裂面平整度较高。Abstract: To explore a method for control over the boulders splitting morphology, cylindrical concrete specimens were penetrated by a linear-shaped charge, the splitting cracks' propagation and development were photographed and the splitting development process of their 3D deformation characteristics were analyzed using high-speed 3D DIC. The results show that the velocity of the crack propagation within the data analysis region exhibits a tendency for a step-by-step increase. The peak velocity and the average velocity are respectively 235.52 m/s and 140.89 m/s. The impact of the linear shaped-charge jet plays a significant role in determining the splitting of the target so that the cracks of the target propagated symmetrically downward in an s-shape along the axis, and on its path the crack showed three obvious inflection points, where branch cracks were produced whose propagation distance was below 5 cm. The shape and location of the main strain concentration determines the crack propagation trend and path, and the tensile strain concentration occurs before the crack appears. The specimen exhibits a quasi-static splitting and the cracks have a fairly even distribution.
-
Key words:
- mechanics of explosion /
- crack /
- High-speed DIC-3D /
- concrete /
- velocity /
- strain
-
多功能含能结构材料(multifunctional energetic structural material, MESM)是综合利用化学能和动能以提高战斗部毁伤效能的新型功能材料[1]。不同于传统惰性破片单一的动能杀伤作用,当MESM破片以一定速度撞击靶板时,会激发剧烈的爆炸/燃烧反应,破片贯穿靶板后,通过动能和化学能的联合作用,实现对靶后目标产生更大的杀伤破坏[2]。
MESM对特定目标的毁伤是现今较活跃的领域之一,研究主要集中在其释能特性和侵彻毁伤效应上。Zhang等[3-4]和Xiong等[5-6]系统研究了Al/Ni基MESM的冲击释能特性及反应机理,并推导了温度控制冲击诱发化学反应的热化学模型[4, 6],结果表明MESM的释能反应特性与破片冲击引发的材料温升有关,反应规律符合Avrami-Erofeev方程。Wang等[7]和Luo等[8]等分别研究了高聚物、W/ZrMESM的冲击释能特性。Xu等[2,9]研究了PTFE/Al/W破片对装甲铝板的侵彻毁伤效应,通过弹道试验结果拟合了PTFE/Al/W破片对装甲铝板的极限穿透速度计算公式,并研究了破片对双层铝板的结构破坏作用。
W/ZrNiAlCu亚稳态合金复合材料是新一代MESM,具有较强的侵彻能力及释能特性,但其对轧制均质装甲(RHA)靶板的侵彻释能特性研究仍属空白。基于此,本文通过弹道枪侵彻实验和高速摄影,测量破片对RHA靶板的撞击速度和靶板贯穿形成的冲塞体速度,记录破片的侵彻释能过程,研究W/ZrNiAlCu亚稳态合金复合材料破片的侵彻释能特性,为该亚稳态合金复合材料的进一步研究和应用提供理论和实验依据。
1. 侵彻实验
1.1 实验方案
实验试样为W/ZrNiAlCu亚稳态合金复合材料破片,材料通过液态真空压力浸渗工艺制备:将熔体ZrNiAlCu亚稳态合金通过高纯氩气吹入由W颗粒粉末烧结而成的基体中,保压0.5 h保证熔体完全浸渗W基体孔隙,随后放入饱和食盐水中淬火,制备出W/ZrNiAlCu亚稳态合金复合材料。通过机加工,得到边长8 mm的立方体破片,破片密度为12.9 g/cm3。
侵彻实验的实验布置如图1所示。破片及弹托由14.5 mm弹道枪发射,弹托刻有凹槽,保证出枪口后弹托、破片分离;10.5 mm厚RHA靶板(GY4)中点位于射击线上,距弹道枪口5 m,一组测速靶置于RHA靶板前方,距靶板1 m,用以测量破片初速,另一组测速靶置于RHA靶板后方,距靶板1 m,用以测量靶板剪切形成的冲塞体速度,每组测速靶由相距0.5 m的两块断通靶纸组成,断通靶纸与计时器相连接;高速摄影机置于靶板侧方,用以记录W/ZrNiAlCu亚稳态合金复合材料破片对RHA靶板的侵彻释能过程。
1.2 实验结果
图2为4个典型撞击速度vi下W/ZrNiAlCu亚稳态合金复合材料破片对RHA靶板的侵彻释能过程,可以看出,破片撞击速度对破片侵彻过程中的释能现象有较大影响。图2(a)中,由于撞击速度较低,破片未贯穿靶板。破片撞击靶板引发了材料的释能反应,靶板前方可看到2束对称的火光并逐渐向外扩散。图2(b)~2(d)中,随着撞击速度增加,靶板前方光束变长,火光范围增大、亮度变高。图2(c)和图2(d)中均可观察到4道光束。当立方体破片撞击靶板时,破片的棱角处应力集中导致材料破碎并向外飞散,因此侧向观察可以看到2~4束对称的火光,火光由材料碎片燃烧产生,并伴随破片飞散而逐渐扩散。当破片撞击速度增大时,撞击能增加,碎片数量增多、初始动能变大、材料温度更高、释能反应更充分,因此靶板前方火光范围增大、亮度变高。
当撞击速度vi=1 067或1 259 m/s时(图2(b)和2(c)),靶板后方呈现一道圆柱型光束,光束顶端可以看到靶板剪切破坏形成的冲塞体。在撞击并贯穿靶板过程中破片破碎,部分碎片穿过靶板并继续飞行,破片碎片由于初始速度、质量、形状的差异,形成碎片云并呈圆柱状分布。当vi=1 508 m/s时(图2(d)),靶板后方圆柱型光束的末端出现椭球状光亮,并在圆柱型光柱消散后继续燃烧(t=3.75 ms)。当破片撞击速度较大时,部分破片碎片温度较高,贯穿靶板并与空气混合后发生燃烧,形成椭球状碎片云,该部分破片温度高,反应充分,燃烧时间较长,在靶后方形成了持续一定时间的椭球状光亮。在一定撞击速度下,W/ZrNiAlCu亚稳态合金复合材料破片既可有效贯穿10.5 mm厚的RHA靶板,又可在靶后发生明显的燃烧反应,是一种具有较强侵彻能力的多功能含能结构材料。需要注意的是,由于侵彻过程中破片的破碎不可忽略,不能使用传统惰性破片侵彻公式解释W/ZrNiAlCu亚稳态合金复合材料破片的侵彻规律。
表1为不同撞击速度vi下W/ZrNiAlCu亚稳态合金复合材料破片对RHA靶板的侵彻结果,m为破片质量,vp为测得的冲塞体速度,p为根据连续介质力学公式计算得到的破片理论冲击压力[10]。可以看到,vi=973 m/s时:破片贯穿靶板;而vi=995 m/s时,破片未贯穿靶板。考虑到破片质量差别及弹道极限速度vbl附近数据量较小,可以认为W/ZrNiAlCu亚稳态合金复合材料破片对10.5 mm厚RHA靶板的弹道极限速度vbl处于810~1 067 m/s范围内。
表 1 破片侵彻实验结果Table 1. Experimental results of fragments penetration编号 m/g vi/(m·s−1) vp/(m·s−1) p/GPa 是否穿透 1# 6.65 726 − 16.46 否 2# 6.54 810 − 18.63 否 3# 6.39 995 − 23.62 否 4# 6.68 973 34 23.01 是 5# 6.80 1 067 137 25.63 是 6# 6.72 1 144 176 27.83 是 7# 6.37 1 159 131 28.26 是 8# 6.65 1 240 195 30.64 是 9# 6.51 1 295 317 32.28 是 10# 6.56 1 359 397 34.21 是 11# 6.68 1 508 497 38.85 是 2. 理论分析
2.1 弹道极限
破片侵彻靶板过程中加速度历程较为复杂,因此采用质量守恒和能量守恒为基础的整体法求解弹道极限速度,做以下三点假设:
(1)破片首先与自由冲塞体发生塑性碰撞,此过程不考虑靶板的绝热剪切效应[11];
(2)破片撞击靶板后部分破碎,靶前飞散碎片不参与侵彻过程;
(3)冲塞体速度vp约为靶后碎片飞出靶板瞬间统计平均速度vr的1.35倍[11]。
破片贯穿靶板前后的能量守恒方程为:
12mv2i=12mf(mm+mpvi)2+12mbv2r+12mpv2p+Efn+Wp (1) 式中:mf为靶板前方飞溅的破片碎片质量之和,mb为破片贯穿靶板后残余质量,mp为靶板冲塞体质量,Efn为破片与靶板塑性碰撞损失能量,Wp为冲塞绝热剪切耗能。弹靶组合确定时,mf仅与撞击速度有关,采用二次多项式的形式拟合:
mf=m[a(viv∗bl)2+b(viv∗bl)+c] (2) 式中:a、b、c为待定系数,
v∗bl 为不考虑质量损失时的弹道极限速度,mf满足质量守恒m=mf+mb。式(1)中,Efn、Wp的表达式分别为[11]:Efn=mpm+mp⋅12mv2i (3) Wp=4dh2⋅Y√3 (4) 式中:d为破片边长;h为靶板厚度;Y=1.2 GPa,为靶板失效应力。由弹道极限速度的定义,当vi= vbl时,vp=vr=0,令a=b=c=0,将式(3)和(4)代入式(1)得vbl*=937.8 m/s;将表1实验数据代入式(1)拟合得到参数a= −2.02, b= 6.47, c=−4.41,解式(1)得vbl=987.1 m/s,在实验所确定的弹道极限速度范围内。
图3为W/ZrNiAlCu亚稳态合金复合材料破片侵彻10.5 mm厚RHA靶板的vi-vp曲线,其中蓝色实线为考虑破片质量损失所得曲线,红色实线为未考虑破片质量损失所得曲线,星型为实验数据,可以看到,vp与vi正相关,蓝色实线始终位于红色实线下方。破片侵彻能力与其初始动能有关,故破片撞击速度vi越大,冲塞体速度vp越大,侵彻能力增强。破片在侵彻靶板过程中的质量损失将导致动能降低,进而降低破片的侵彻能力,因此蓝色实线位于红色实线下方,理论弹道极限速度vbl>
v∗bl ,破片质量损失对侵彻能力有较大影响。图上星状点为实验数据,可以看到,理论计算结果与实验数据吻合较好,表明所做假设合理,式(1)可以较好地解释包含质量损失的破片对中厚靶板的侵彻规律。2.2 释能特性
由Grünesien物态方程和Rankin-Hugoniot能量方程得到固体p-V形式的物态方程[6]:
p(V)=Vγ(V)pc(V)−Ec(V)Vγ(V)−12(V0−V) (5) 式中:p为材料压力,V为比容,
γ 为Grünesien系数,pc、Ec分别为材料冷压和冷能。通过Meyer势,金属材料的冷压和冷能可以表示为材料常数Q-q的形式[6]:Ec(δ)=3Qρ0K{1qexp[q(1−δ−1/3)]−δ−1/3−1q+1} (6) pc(δ)=Qδ2/3{exp[q(1−δ−1/3)]−δ2/3} (7) 式中:
δ=ρ/ρ0K=V0K/V ,为材料的相对于0 K时密度的压缩度,ρ0K 、V0K 为0 K时的密度和比容。Dugdale -MacDonald关系式的Q-q形式为[6]:
γ(V)=16⋅q2δ−1/3exp[q(1−δ−1/3)]−6δqexp[q(1−δ−1/3)]−2δ (8) 采用胡金彪等给出的解析法可以求得材料常数Q、q[6]:
s[1+(s2−18γ20Ks−1)αVT0]=112⋅q2+6q−18q−2 (9) {C0[1+(2s−γ20K4−1)αVT0]}2=Q(q−2)3ρ0K (10) 式中:C0、s为Hugoniot参数,其中C0理论上等于材料的零压体积声速;
γ0K 为0 K时的Grünesien系数;αV 为材料的体积膨胀系数;T0为起始温度,即实验室室温,取298 K。Boslough根据固态材料绝热线和等熵线的关系,假定材料等容比热为常数,得到冲击波温度表达式[12]:
TH=T0exp(−∫VHV0VγdV)+VγpH−pSCV (11) 式中:TH为击波温度,pH为击波压力,V0、VH分别为初态、击波压缩状态下的比容,CV为材料的等容比热,pS为与冲击绝热线同一起始状态出发的等熵线压力,其解析式为[12]:
pS(η)=ρ0C20exp(γ0η)∫η0(1+(s−γ0))x(1−sx)3exp(−γ0x)dx (12) 式中:
η=1−V/V0 为压缩度,x为压缩度的积分变量,ρ0 为起始状态下材料密度,γ0 为初始Grünesien系数。将式(5)、式(11)联立,即可得到击波压力与击波温度的关系。假设含能结构材料的化学反应过程仅受击波温度控制,研究表明,含能结构材料在高升温速率下的固态反应动力学满足n维Avrami-Erofeev方程,经过整理,Avrami-Erofeev方程可以表示为击波温度对反应效率一次微分的形式[6]:
dTdy=RuT2Ea[12y−nln(1−y)+n−1n(1−y)[−ln(1−y)]] (13) 式中:y为反应效率,Ea为表观活化能,Ru为理想气体常数,n为与边界条件和反应机理有关的系数。
表2为计算所需W/ZrNiAlCu亚稳态合金复合材料参数,图4为破片冲击压力与反应效率关系的理论曲线,可以看到,在40 GPa冲击压力范围内,材料并未达到完全反应,当p=40 GPa,材料反应效率y=0.541;材料反应效率与冲击压力正相关,因此,随着冲击速度增加,破片冲击释能所产生的火光范围、火光亮度增大,与实验现象相符合。
表 2 W/ZrNiAlCu亚稳态合金复合材料参数[13]Table 2. Parameters of W/ZrNiAlCu metastable reactive alloy composite材料 C0/(m∙s−1) s γ0K αV/(10−5·K−1) CV/(J·g−1·K−1) Ea/(kJ·mol−1) n W/ZrNiAlCu 3 417 1.732 1.629 1.559 0.182 459.25 0.347 3. 结 论
(1)W/ZrNiAlCu亚稳态合金复合材料破片可以贯穿10.5 mm厚的RHA靶板,破片撞击、贯穿靶板过程中导致材料破碎并激发燃烧反应,靶前、靶后均可观察到大范围火光;随着撞击速度增加,火光范围增大,亮度提高,燃烧反应时间变长。
(2)推导了W/ZrNiAlCu亚稳态合金复合材料破片侵彻规律方程,得到了破片侵彻10.5 mm厚RHA靶板的理论弹道极限速度vbl=987.1 m/s,理论计算结果与实验数据吻合较好;破片侵彻过程中的质量损失对其侵彻能力有较大影响,在相同撞击速度下,W/ZrNiAlCu亚稳态合金复合材料破片侵彻形成冲塞体速度小于不考虑质量损失的理想情况下冲塞体速度。
(3)由材料的Grunesien物态方程和Avrami-Erofeev方程得到材料反应效率y与冲击压力p的理论关系,在15~40 GPa范围内,材料反应效率随着冲击压力的增加而增加,当p=40 GPa时,y=0.541,因此实验过程中材料反应并不完全,出现随着撞击速度增加,火光范围增大、亮度提高的现象。
-
[1] 钱伟长.穿甲力学[M].北京:国防工业出版社, 1984. [2] Alekseevskii V P. Penetration of a rod into a target at high velocity[J]. Combustion, Explosion, and Shock Waves, 1966, 2(2):63-66. doi: 10.1007-BF00749237/ [3] Sternberg J. Material properties determining the resistance of ceramics to high velocity penetration[J]. Journal of Applied Physics, 1989, 65(9):3417-3424. doi: 10.1063/1.342659 [4] 威廉·普·沃尔特斯, 乔纳斯·埃·朱卡斯.成型装药原理及其应用[M].王树魁, 贝静芬, 译.北京:兵器工业出版社, 1992. [5] 肖强强, 黄正祥, 祖旭东.双材质复合射流对混凝土的侵彻[J].爆炸与冲击, 2014, 34(4):457-463. doi: 10.11883/1001-1455(2014)04-0457-07Xiao Qiangqiang, Huang Zhengxiang, Zu Xudong. Penetration of jacketed jet into concret[J]. Explosion and Shock Waves, 2014, 34(4):457-463. doi: 10.11883/1001-1455(2014)04-0457-07 [6] Kirugulige M S, Tippur H V. Measurement of fracture parameters for a mixed-mode crack driven by stress waves using image correlation technique and high-speed digital hotography[J]. Strain, 2009, 45(2):108-122. doi: 10.1111/str.2009.45.issue-2 [7] Chevalier L, Calloch S, Hild F, et al. Digital image correlation used to analyze the multiaxial behavior of rubber-like materials[J]. European Journal of Mechanics-A: Solids, 2001, 20(2):169-187. doi: 10.1016/S0997-7538(00)01135-9 [8] Erdogan F, Sih G C. On the crack extension in plates under plane loading and transverse shear[J]. Journal of Basic Engineering, 1963, 85(4):519-525. doi: 10.1115/1.3656897 [9] 岳中文, 杨仁树, 郭东明, 等.爆炸应力波作用下缺陷介质裂纹扩展的动态分析[J].岩土力学, 2009, 30(4):949-954. http://d.old.wanfangdata.com.cn/Periodical/ytlx200904015Yue Zhongwen, Yang Renshu, Guo Dongming, et al. Dynamic analysis of crack propagation in media containing flaws under the explosive stress wave[J]. Rock and Soil Mechanics, 2014, 35(8):2048-2414. http://d.old.wanfangdata.com.cn/Periodical/ytlx200904015 [10] 肖正学, 张志呈, 郭学彬.断裂控制爆破裂纹发展规律的研究[J].岩石力学与工程学报, 2002, 21(4):546-549. doi: 10.3321/j.issn:1000-6915.2002.04.019Xiao Zhengxue, Zhang Zhicheng, Guo Xuebin. Research on crack developing law of rock fracture controlled blasting[J]. Chinese Journal of Rock Mechanics and Engineering, 2002, 21(4):546-549. doi: 10.3321/j.issn:1000-6915.2002.04.019 [11] 戴俊.岩石动力学特性与爆破理论[M].北京:冶金工业出版社, 2002. [12] 王汉军, 黄风雷, 张庆明.岩石定向断裂爆破的力学分析及参数研究[J].煤炭学报, 2003, 28(4):399-402. http://d.old.wanfangdata.com.cn/Periodical/mtxb200304014Wang Hanjun, Huang Fenglei, Zhang Qingming. Mechanics effect analysis and parameters study on borehole directional fracture blasting[J]. Journal of China Coal Society, 2003, 28(4):309-402. http://d.old.wanfangdata.com.cn/Periodical/mtxb200304014 [13] Chao Y J, Sutton M A, Peters W H, et al. Measurement of three-dimensional displacements in deformable bodies by digital image processing[C]//1989 SEM Spring Conference on Experimental Mechanics. Cambridge, 1989: 139-146. [14] Kolsky H, Shearman A C. Investigation of fractures produced by transient stress waves[J]. Research, 1949, 2(8):384-389. [15] 杨建华, 卢文波, 胡英国, 等.隧洞开挖重复爆炸荷载作用下围岩累积损伤特性[J].岩土力学, 2014, 35(2):512-518. http://d.old.wanfangdata.com.cn/Periodical/ytlx201402030Yang Jianhua, Lu Wenbo, Hu Yingguo, et al. Accumulated damage in surrounding rocks due to repeated blasting loads during blasting excavation of tunnels[J]. Rock and Soil Mechanics, 2014, 35(2):512-518. http://d.old.wanfangdata.com.cn/Periodical/ytlx201402030 [16] 严成增, 郑宏, 孙冠华, 等.基于数字图像技术的岩土材料有限元-离散元分析[J].岩土力学, 2014, 35(8):2048-2414. http://d.old.wanfangdata.com.cn/Periodical/ytlx201408038Yan Chengzeng, Zheng Hong, Sun Guanhua, et al. FDEM of geomaterials based on digital image technology[J]. Rock and Soil Mechanics, 2014, 35(8):2048-2414. http://d.old.wanfangdata.com.cn/Periodical/ytlx201408038 期刊类型引用(8)
1. 贾杰,智小琦,郝春杰,李劲,郭璐,柳星河. Zr基非晶破片对碳纤维复合靶及后效铝靶的侵彻试验研究. 高压物理学报. 2024(02): 130-139 . 百度学术
2. 杨林,于述贤,范群波. Zr_(77.1)Cu_(13)Ni_(9.9)非晶合金破片侵彻LY12铝合金及TC4钛合金靶板毁伤后效及机理对比研究. 北京理工大学学报. 2023(04): 417-428 . 百度学术
3. 熊玮,张先锋,李逸,谈梦婷,刘闯,侯先苇. 活性材料冲击压缩及反应行为模拟方法研究进展. 北京理工大学学报. 2023(10): 995-1015 . 百度学术
4. 芦永进,梁增友,邓德志,朱聪. 铜基非晶合金双层药型罩射流形成及侵彻性能. 火炮发射与控制学报. 2022(01): 14-20+28 . 百度学术
5. 郭志平,王飞,姜波,张杰,程波,王传婷,何勇. ZrCuAlNi合金的Taylor撞击断裂行为研究. 兵器装备工程学报. 2022(05): 185-190 . 百度学术
6. 尚春明,施冬梅,张云峰,徐雪涛. Zr基非晶合金毁伤研究进展. 兵器装备工程学报. 2020(07): 182-186 . 百度学术
7. 陈海华,张先锋,熊玮,刘闯,魏海洋,汪海英,戴兰宏. WFeNiMo高熵合金动态力学行为及侵彻性能研究. 力学学报. 2020(05): 1443-1453 . 百度学术
8. 何丽灵,张方举,颜怡霞,谢若泽,徐艾民,周燕良. Ti-6Al-4V弹体破坏模式对冲击反应的影响研究. 爆炸与冲击. 2020(12): 58-69 . 本站查看
其他类型引用(1)
-