Quantitative Determination of Impact Reaction Energy Release for HfZrTiTaNb Based High-Entropy Alloys
-
摘要: 作为一种新型含能材料,高熵合金在高速冲击过程中会释放大量能量,具有重要的应用价值。采用二级轻气炮系统加载真空环境下的HfZrTiTaNb系高熵合金弹丸,对轴承钢进行了冲击实验,测量了闪光辐射温度、气体超压、火焰传播速度和容器壁温升等响应参数的演化过程,分析了高熵合金弹丸撞击靶板过程中的能量流向,定量计算了高熵合金在密闭容器内冲击反应的混合气体焓、闪光辐射能、容器壁的吸收能、喷出气体焓以及撞击靶板产生的变形能,得到了不同元素及其含量对高熵合金释能量的影响。结果表明,高熵合金弹丸冲击反应释放的能量主要被准密闭容器壁吸收。随着Cu或Al含量的增加,HfZrTiTaNb系高熵合金的单位质量释能量增加。在相近的撞击速度下,含Cu高熵合金的单位质量释能量比含Al高熵合金大。Abstract: As a new type of energetic material, high-entropy alloys will release a large amount of energy during high-speed impact, which has important application value. A two-stage light gas gun system was used to load the HfZrTiTaNb high-entropy alloys projectile under vacuum environment, and the impact experiment of bearing steel was carried out. The evolution process of response parameters such as flash radiation temperature, gas overpressure, flame propagation velocity and temperature rise of container wall was measured. The energy flow direction during the impact reaction of the high-entropy alloys projectile impacting the target plate was analyzed. The enthalpy of the mixed gas, the flash radiation energy, the absorption energy of the container wall, the enthalpy of the ejected gas and the deformation energy generated by the impact on the target during the impact reaction of the high-entropy alloys in a closed container were quantitatively calculated. The effects of different elements and their contents on the energy release of high-entropy alloys were obtained. The results showed that the energy released by the impact reaction of high-entropy alloys projectiles was mainly absorbed by the quasi-closed container wall. With the increase of Cu or Al content, the unit mass release energy of HfZrTiTaNb based high-entropy alloys increased. At similar impact velocities, the high-entropy alloys containing Cu released more energy per unit mass than the high-entropy alloys containing Al.
-
针对现代战场环境毁伤目标多样化的需求,具备多毁伤元打击能力的复合战斗部(Composite warhead)应运而生。复合战斗部将多种不同的战斗部(如破片战斗部、爆炸成型弹丸和空心装药战斗部等)结合在一起,可以同时对多种目标进行有效打击,战斗部利用效率高、通用化程度高、适应性强,已成为目前战斗部技术发展的一个重要方向。
近年来,针对复合毁伤战斗部的研究逐渐增多。王利侠等[1]利用新型薄型波形控制器、半预制壳体和精密破甲战斗部技术,设计了破甲/杀伤多用途战斗部,使单兵战斗部在基本保持原有穿深威力的基础上,增加了杀伤功能。李兴隆等[2]则在此基础上,进一步研究了战斗部装填系数(装药质量与壳体质量之比)对破片杀伤半径和成型射流破甲深度的影响。然而,更多的复合战斗部还是通过在原有聚能战斗部上增设预制破片,实现聚能侵彻和破片杀伤相结合。张俊等[3]提出了一种轴向爆炸成型弹丸(Explosively formed projectile,EFP)、周向预制破片的新型破甲杀伤复合战斗部结构,并仿真分析了该结构下EFP的成型性能与破片场的杀伤威力。龚柏林等[4]提出了一种将预制破片周向分布在药型罩端面的复合战斗部结构,使预制破片与EFP均沿轴线飞行。基于此类结构,常江等[5]对药型罩设计参数进行了针对性的优化;董晓亮等[6]分析了预制破片类型、直径、材料、间距等对复合毁伤元战斗部发散角的影响;Ji等[7]分析了Ta-2.5W合金药型罩和钨铬合金预制破片的可行性及性能。此外,陈闯等[8]基于战斗部毁伤威力可控技术,进一步提出了一种外围采用刻槽壳体且可展开不同角度的EFP战斗部,并对EFP和预控破片群相结合的复合破片场成型过程进行了数值模拟研究,分析了中心与外围装药量、外围装药起爆位置、展开角度对复合破片场成型特性的影响。总体而言,当前研究多集中于战斗部结构本身的影响及优化,较少从起爆能量的角度分析其毁伤效能。
为了进一步提高复合战斗部的能量输出效率,本研究利用LS-DYNA数值仿真软件,研究起爆点位置、起爆直径和起爆点数量对复合战斗部所形成的聚能侵彻体、自然破片和预制破片3种毁伤元成型及性能的影响,讨论实现战斗部毁伤威力可调的技术路径。
1. 战斗部计算模型及研究方案
1.1 战斗部结构模型
本研究中的复合战斗部主要包括主装药、药型罩、预制破片、壳体等部分,可形成聚能侵彻体、预制破片和自然破片3种毁伤元。其装药口径Dk为100 mm,装药高度为90 mm。药型罩采用弧锥结合结构,锥角155°,曲率半径55 mm,罩厚4.3 mm;预制破片为方块,边长为8 mm,共9排324枚;壳体厚度为4 mm,具体结构如图1所示。
主装药为JH-2炸药,采用High_Explosive_Burn模型和JWL状态方程描述。JWL状态方程表示为
p=A(1−ωR1V)e−R1V+B(1−ωR2V)e−R2V+ωe0V 式中:A、B、R1、R2和
ω 为输入参数,V为相对体积,e0为炸药材料的初始比内能。炸药材料参数如表1[9]所示,其中: ρe 为炸药密度,D为爆速,pC-J为C-J爆轰压力。空气选用Null材料模型和Grüneisen状态方程描述。对于压缩材料,Grüneisen状态方程表示为
p=ρ0C2μ[1+(1−γ02)μ−a2μ2][1−(S1−1)μ−S2μ2μ+1−S3μ3(μ+1)2]2 对于膨胀材料,则式(2)进一步表征为
p=ρ0C2μ+(γ0+aμ)E 式中:C为材料的静态体积声速;S1、S2、S3为系数;
γ0 为Grüneisen常数;μ=ρa/ρ0−1 ,ρa 为空气密度,ρ0 为空气初始密度;a为γ0 的一阶体积修正,E为空气的初始比内能。空气材料参数如表2[9]所示。金属部件主要包括药型罩、预制破片和壳体,所用材料分别为紫铜、钨合金和45钢。毁伤元成型过程属于典型的高温、高应变过程,因此选择Johnson-Cook本构模型和Grüneisen状态方程描述。Johnson-Cook本构方程式表示为
σy=(σ+b¯εpn)(1+cln˙ε∗)(1−T∗m) 式中:
σy 为材料的流动屈服强度,σ 为准静态下的屈服应力,b为应变硬化系数,c为应变率相关系数,m为温度相关系数,n为应变硬化指数,¯εp 为等效塑性应变,˙ε∗ 为无量纲化后的等效塑性应变率,T∗ 为无量纲温度。金属材料参数见表3[10],其中:ρ 为材料密度,G为剪切模量。1.2 数值仿真方案
利用ANSYS/LS-DYNA软件对复合毁伤元战斗部的成型过程进行数值模拟。考虑到模型为轴对称结构,取1/4结构建立三维模型。数值模拟过程采用ALE算法,其中药型罩、炸药、空气选择Euler网格;预制破片及壳体采用Lagrange网格,置于Euler网格中。Euler网格与Lagrange网格之间采用流体/结构耦合方式[9]。具体的网格划分为:药型罩径向60等分,厚度划分为6份,共16 200个网格;装药径向同样60等分,轴线划分为50份,共135 000个网格;空气与装药、药型罩相接位置划分方式一致;其他位置网格大小沿EFP成型方向分别取1.5、2.0和3.0 mm,共1 035 600个网格;预制破片网格大小为0.2 mm,共5 184个网格;考虑自然破片的形成,壳体网格大小取1.0 mm,共62 928个网格。
预制破片和壳体之间采用自动面-面接触方式,钨球定义为自动单面接触[6]。为避免压力在边界面反射,在整个空气计算域的周围表面定义无反射边界。考虑到战斗部作用过程属于典型爆炸问题,时间步长取0.67 μs。
本研究主要分析起爆点位置、起爆直径和起爆点数量对复合战斗部毁伤元成型及性能的影响,涉及到的起爆形式主要包括:端面单点起爆、端面多点起爆、轴线多点同步起爆、轴线阵列多点同步起爆等,起爆方式如图2所示。
单点起爆指装药端面中心点处起爆,多点起爆除特别说明外默认为8点起爆。起爆位置如图2(a)所示,分别取装药顶点(P1)、装药高度60 mm处(P2)、装药高度30 mm处(P3)和药型罩顶点(P4)4个位置;对于起爆直径,如图2(b)所示,分别取装药端面中心点、Dk/4、 Dk/2、3Dk/4和装药直径Dk处;对于起爆数量,如图2(c)所示,分别取端面单点、4点、6点、8点和轴线1 × 3点、4 × 3点、6 × 3点、8 × 3点。
2. 起爆方式对复合毁伤元的影响
对于该复合战斗部结构,其毁伤元主要包括聚能侵彻体、自然破片和预制破片3种形式。对于聚能侵彻体,主要成型参数为头部速度、头尾速度差和长径比,此处长径比为侵彻体总长与最大直径之比;对于破片,则主要分析其成型状态及速度分布,以探究起爆方式对毁伤元成型的影响。
2.1 起爆点位置的影响
为了研究起爆点位置对毁伤元成型的影响,分别研究起爆点位置为P1、P2、P3、P4时中心单点起爆及端面多点起爆对复合战斗部毁伤元成型的影响。端面多点起爆采用8点同步起爆,且起爆直径为装药直径;轴线多点起爆为P1、P2、P3处3点同步起爆,采用P1-2-3来表示。通过数值模拟计算得到150 μs时刻各方案聚能侵彻体成型参数及形态,如图3所示,此时聚能侵彻体成型已基本稳定,形态不再发生明显变化。
由图3可知,随着起爆点位置逐渐向药型罩靠近,单点起爆形成的聚能侵彻体头部速度由2 434 m/s减小至2 165 m/s,而多点起爆形成的聚能侵彻体头部速度则由3 640 m/s减小至不足2 000 m/s。单点起爆时聚能侵彻体头尾速度差较小,始终保持在50 m/s以下,形成了典型的EFP;多点起爆时形成的聚能侵彻体头尾速度差更大,甚至形成了头尾速度差约为1 100 m/s的聚能杆式侵彻体(Jetting projectile charge,JPC)。起爆点位置离药型罩越远,形成的聚能侵彻体头部长度越大,其头尾速度差也越大,越有利于JPC的成型。轴线多点起爆时,其聚能侵彻体成型状态与P3处相近。这说明对于聚能侵彻体,轴线多点起爆时,距药型罩最近的起爆点受到的影响更大。在进行聚能毁伤元形态仿真研究的同时,对预制破片的飞散情况也进行了研究,150 μs时刻预制破片的速度参数如图4所示。
由图4可知,不论中心单点起爆还是端面多点起爆,起爆点位置由装药端面向药型罩逼近的过程中,预制破片的最高速度均呈现出先增大后减小的规律。整体来看,虽然多点起爆与单点起爆的最低速度相近,但多点起爆的速度分布带更大,均在200 m/s以上,对预制破片最高速度的提升作用显著,速度最高可提升100 m/s以上。对于复合战斗部,壳体形成的自然破片也至关重要,150 μs时刻自然破片成型参数及形态如图5所示。
由图5可知,与预制破片的规律相似,起爆点位置由装药端面向药型罩逼近的过程中,自然破片平均速度也呈现出先增大后减小的规律,且两种起爆模式下仍然是P2处起爆时速度最高,分别达到728和738 m/s。对于破片质量分布,相较于端面多点起爆,中心单点起爆时壳体破碎形成的自然破片分布更均匀;而端面起爆时破片质量更大,也更完整,有利于提高破片威力。
2.2 起爆直径的影响
为了研究起爆直径对毁伤元成型的影响,分别研究了不同起爆直径时端面8点起爆和轴线阵列多点起爆对复合战斗部毁伤元成型的影响,轴线阵列起爆时P1、P2、P3同步24点(3 × 8)起爆。起爆直径分别取0 mm(中心单点)、Dk/4、Dk/2、3Dk/4和Dk。通过数值模拟计算得到150 μs时刻各方案聚能侵彻体成型参数及形态,如图6所示。
由图6可知,端面起爆时形成的侵彻体头尾速度均大于轴向阵列起爆时所形成的侵彻体,相较于轴向阵列起爆形成的EFP,其更易形成JPC。随着起爆直径的增大,两种起爆模式下形成的聚能侵彻体头部速度和头尾速度差均增大,杆型头部逐渐增长,长径比大幅提高,毁伤元由EFP形态逐渐转变为JPC形态。而端面起爆增益效果更为明显,头部速度从2 434 m/s提高到3 640 m/s的同时,头尾速度差也从零增大至570 m/s;侵彻体长径比达到2.94,是中心点起爆时长径比的4.8倍。150 μs时刻预制破片速度参数如图7所示。
由图7可知,不论是装药端面起爆还是轴向阵列起爆,当起爆直径增大时,预制破片的最高速度均增大,分别提高63和86 m/s;但端面起爆时预制破片的最小速度逐渐减小至585 m/s,轴向阵列起爆时预制破片的最小速度却逐渐增大至621 m/s,且轴向阵列起爆的速度分布带更大,对预制破片最高速度的提升作用明显,最高提升速度达到100 m/s以上。150 μs时刻自然破片成型参数及形态如图8所示。
由图8可知,与预制破片不同,随着起爆直径的增大,端面起爆自然破片速度先减小后增大,3Dk/4处速度最低约690 m/s;轴向3处圆形阵列起爆自然破片的速度逐渐增大,速度提高约50 m/s,达到729 m/s。起爆直径为Dk/2时,两种起爆方式下自然破片的速度近乎一致,均为696 m/s。装药端面起爆直径存在一个最值,此时破片平均速度最低,且随着起爆直径的增大,大尺寸破片增加。
2.3 起爆点数量的影响
为了研究起爆点数量对毁伤元成型的影响,分别研究了不同起爆点数量时装药端面及轴向阵列起爆对复合战斗部毁伤元成型的影响。装药端面起爆时,起爆点数量分别取中心单点、端面4点、6点和8点;轴向阵列起爆时,取3点(1 × 3)、12点(4 × 3)、18点(6 × 3)和24点(8 × 3)。通过数值模拟计算得到150 μs时刻各方案聚能侵彻体成型参数及形态,如图9所示。
由图9可知,端面单点和轴向阵列起爆容易形成EFP,而端面多点起爆更易形成JPC。随着起爆点数量的增加,聚能侵彻体头部速度增大,端面起爆和轴向阵列起爆的速度增益分别达到1 000 m/s以上和200 m/s;长径比也随之增大,分别由0.62和0.51增大至2.58和0.96,但随着起爆点数量增加,长径比增大速度逐渐放缓,端面起爆点数量由4个增加到8个时,长径比仅提高约10%;轴向阵列起爆形成的侵彻体头尾速度增长较为缓和,头尾速度和长度均远小于端面起爆形成的侵彻体,最大头尾速度差不超过120 m/s,只能形成较典型的EFP。150 μs时刻预制破片的速度参数如图10所示。
由图10可知,随着起爆点数量的增加,端面起爆时预制破片的最大速度由743 m/s逐渐提高到803 m/s,而最小速度则由620 m/s减小到585 m/s;轴向阵列起爆时预制破片的最大速度及最小速度均先增大后减小,最大时分别为831和626 m/s。整体来看,端面起爆预制破片速度区间随着起爆点数量的增加而增大,速度区间由123 m/s增大至211 m/s,增长70%以上;但轴向阵列起爆的最大速度更大,达到831 m/s,对预制破片的最高速度提高了10%,作用更加显著。壳体自然破片在150 μs时刻的成型参数及形态如图11所示。
由图11可知,与预制破片的规律相似,随着起爆点数量增加,轴向阵列起爆时破片速度先增大后缓慢减小。但端面起爆时形成的自然破片平均速度变化不大,始终为710 m/s左右,上下浮动在10 m/s以内。壳体质量分布方面,随着起爆点数量增加,壳体破片质量增大,且破片更完整,更有利于破片威力的提高。总体来看,端面起爆时起爆点数量的作用不明显,轴线阵列起爆时增加起爆点数量有利于破片速度的提高。
3. 结 论
通过数值仿真研究了起爆点位置、直径和起爆点数量对复合战斗部所形成的聚能侵彻体、自然破片和预制破片3种毁伤元成型及性能的影响,得出以下主要结论。
(1)随着起爆点位置向药型罩靠近,聚能侵彻体长径比和头部速度均呈现减小趋势,端面8点起爆时尤为明显,且头尾速度差明显减小;单点起爆时位置的影响较小,P1-2-3同时起爆的效果并不明显,与P3处起爆形成的EFP相近。预制破片均在P2处起爆时达到最大破片速度,且端面8点起爆的速度增益更为明显。对于壳体形成的自然破片,其平均速度也在P2处达到最大;多点起爆位置对应壳体处形成的破片最大;P1-2-3同时起爆虽不能提高破片速度,但有助于使自然破片的断裂更均匀。
(2)随着起爆直径的增大,端面起爆和轴线阵列起爆对聚能侵彻体的影响一致,其长径比、头部速度和头尾速度差均随之增大,速度增益最高可达50%;端面起爆更有利于聚能侵彻体的拉伸,甚至可以形成JPC。预制破片的最大速度均呈上升趋势,但最小速度呈相反的变化趋势。对于壳体形成的自然破片,装药端面起爆直径存在一个最值,此时破片平均速度最低,且随着起爆直径的增大,大尺寸破片增加;而轴线阵列起爆直径的增大有利于破片速度的提高。
(3)随着起爆点数量的增加,聚能侵彻体长径比、头部速度和头尾速度差同样随之增大,速度增益与起爆直径的影响一致,也形成了JPC,但随着起爆点数量的增加,速度增大的趋势逐渐减缓。预制破片的速度变化与起爆直径的影响相近,但具有平台现象,整体来看,两种起爆方式对预制破片速度的增益相近。对于壳体形成的自然破片,端面起爆点数量的变化对破片平均速度的影响不大,仅对破片质量分布有一定影响,轴线阵列起爆点数量的增加有利于破片速度的提高。
-
表 1 实验基本参数
Table 1. Basic parameters of the experiment
No. Specimen Size/(mm×mm) Mass/g Impact velocity/(km·s−1) Density/(g·cm−3) 1 HfZrTiTaNbCu0.2 ∅4.63×4.68 0.786 1.35 9.98 2 HfZrTiTaNbCu0.8 ∅4.96×4.69 0.918 1.30 10.13 3 HfZrTiTaNbAl0.2 ∅4.95×4.83 0.899 1.31 9.68 4 HfZrTiTaNbAl0.8 ∅4.81×5.15 0.853 1.57 9.12 5 HfZrTiTaNbAl0.8 ∅4.76×5.11 0.832 1.59 9.15 表 2 瞬态光纤高温计各通道的标定数据
Table 2. Calibration data for each channel of the transient fiber pyrometer
Channel Wave length/nm Impedance/kΩ Voltage/V Nr/(μW·cm−2·nm−1) 1 400 10 2.934 281.7 2 500 10 4.400 235.5 3 600 1 0.880 347.1 4 700 1 6.046 280.7 表 3 实验中的超压和闪光辐射温度峰值
Table 3. Peak of overpressure and flash radiation temperature in the experiment
No. Specimen Overpressure/MPa Temperature/K 1 HfZrTiTaNbCu0.2 0.22 1260 2 HfZrTiTaNbCu0.8 0.24 1370 3 HfZrTiTaNbAl0.2 0.21 1120 4 HfZrTiTaNbAl0.8 0.20 1228 5 HfZrTiTaNbAl0.8 0.20 1281 表 4 冲击反应释能结果
Table 4. Impact reaction energy release results
No. Specimen Impact velocity/
(km·s−1)Mass/g Energy release/kJ Energy release per
unit mass/(kJ·g−1)1 HfZrTiTaNbCu0.2 1.35 0.786 1.32 1.68 2 HfZrTiTaNbCu0.8 1.30 0.918 1.76 1.92 3 HfZrTiTaNbAl0.2 1.31 0.899 0.88 0.98 4 HfZrTiTaNbAl0.8 1.57 0.853 1.19 1.40 5 HfZrTiTaNbAl0.8 1.59 0.832 1.21 1.45 表 5 HfZrTiTaNb系高熵合金的原子尺寸差
Table 5. Atomic size difference of HfZrTiTaNb based high-entropy alloys
x δ/% HfZrTiTaNbCux HfZrTiTaNbAlx 0.2 4.92 4.84 0.8 6.73 6.37 -
[1] 刘丘林, 刘允中, 王艳群. 高熵合金的研究现状和应用前景 [J]. 粉末冶金工业, 2017, 27(6): 64–69. doi: 10.13228/j.boyuan.issn1006-6543.20160115LIU Q L, LIU Y Z, WANG Y Q. Research status and application prospect of high-entropy alloy [J]. Powder Metallurgy Industry, 2017, 27(6): 64–69. doi: 10.13228/j.boyuan.issn1006-6543.20160115 [2] ZHANG W R, LIAW P K, ZHANG Y. Science and technology in high-entropy alloys [J]. Science China Materials, 2018, 61(1): 2–22. doi: 10.1007/s40843-017-9195-8 [3] ZHANG Z R, ZHANG H, TANG Y, et al. Microstructure, mechanical properties and energetic characteristics of a novel high-entropy alloy HfZrTiTa0.53 [J]. Materials & Design, 2017, 133: 435–443. doi: 10.1016/j.matdes.2017.08.022 [4] GEANTĂ V, VOICULESCU I, STEFĂNOIU R, et al. Dynamic impact behaviour of high entropy alloys used in the military domain [J]. IOP Conference Series: Materials Science and Engineering, 2018, 374: 012041. doi: 10.1088/1757-899X/374/1/012041 [5] 王睿鑫. NbZrTiTa高熵合金的组织结构演变及结构释能特性研究 [D]. 长沙: 国防科技大学, 2018.WANG R X. Microstructure evolution and energetic structural properties of NbZrTiTa high-entropy alloy [D]. Changsha: National University of Defense Technology, 2018. [6] REN K R, LIU H Y, CHEN R, et al. Compression properties and impact energy release characteristics of TiZrNbV high-entropy alloy [J]. Materials Science and Engineering: A, 2021, 827: 142074. doi: 10.1016/j.msea.2021.142074 [7] MA Y S, ZHOU L, ZHANG K C, et al. Effects of cerium doping on the mechanical properties and energy-releasing behavior of high-entropy alloys [J]. Materials, 2022, 15(20): 7332. doi: 10.3390/ma15207332 [8] CHEN C, GAO R K, GUO K, et al. Quantitative determination of impact release energy for TiZrHfX0.3 multicomponent materials in vacuum environment [J]. International Communications in Heat and Mass Transfer, 2022, 133: 105958. doi: 10.1016/j.icheatmasstransfer.2022.105958 [9] TANG E L, WANG Q C, CHEN C, et al. Quantitative reactive release energy models of Al/Teflon material during high velocity impact [J]. International Journal of Energy Research, 2022, 46(15): 23069–23082. doi: 10.1002/er.8606 [10] 张宝平, 张庆明, 黄风雷. 爆轰物理学 [M]. 北京: 兵器工业出版社, 2009.ZHANG B P, ZHANG Q M, HUANG F L. Detonation physics [M]. Beijing: Weapon Industry Press, 2009. [11] HAN Y F, TANG E L, HE L P, et al. Evolutionary characteristics of thermal radiation induced by 2A12 aluminum plate under hypervelocity impact loading [J]. International Journal of Impact Engineering, 2019, 125: 173–179. doi: 10.1016/j.ijimpeng.2018.11.013 [12] CHEN C, TANG E L, ZHU W J, et al. Modified model of Al/PTFE projectile impact reaction energy release considering energy loss [J]. Experimental Thermal and Fluid Science, 2020, 116: 110132. doi: 10.1016/j.expthermflusci.2020.110132 [13] AMES R. Vented chamber calorimetry for impact-initiated energetic materials [C]//The 43rd AIAA Aerospace Sciences Meeting and Exhibit. Reno: AIAA, 2005: 275–279. [14] 高学平. 高等流体力学 [M]. 天津: 天津大学出版社, 2005.GAO X P. Advanced fluid mechanics [M]. Tianjin: Tianjin University Press, 2005. [15] 唐恩凌, 徐名扬, 张庆明, 等. 超高速撞击厚靶过程的能量分配研究 [J]. 固体力学学报, 2016, 37(2): 152–160. doi: 10.19636/j.cnki.cjsm42-1250/o3.2016.02.005TANG E L, XU M Y, ZHANG Q M, et al. Study on partitioning of energy in hypervelocity impact on thick target [J]. Chinese Journal of Solid Mechanics, 2016, 37(2): 152–160. doi: 10.19636/j.cnki.cjsm42-1250/o3.2016.02.005 [16] 段春争, 秦泗伟. 淬硬GCr15钢在高温和高应变率下的动态再结晶动力学模型 [J]. 金属热处理, 2017, 42(2): 34–38. doi: 10.13251/j.issn.0254-6051.2017.02.008DUAN C Z, QIN S W. Dynamic recrystallization kinetics model of hardened GCr15 steelat high temperature and high strain rate [J]. Heat Treatment of Metals, 2017, 42(2): 34–38. doi: 10.13251/j.issn.0254-6051.2017.02.008 [17] 张路明, 马胜国, 李志强, 等. Al x CoCrFeNi高熵合金力学性能的分子动力学模拟 [J]. 高压物理学报, 2021, 35(5): 052201. doi: 10.11858/gywlxb.20210730ZHANG L M, MA S G, LI Z Q, et al. Mechanical properties of Al x CoCrFeNi high-entropy alloy: a molecular dynamics study [J]. Chinese Journal of High Pressure Physics, 2021, 35(5): 052201. doi: 10.11858/gywlxb.20210730 -