HfZrTiTaNb系高熵合金的冲击反应释能定量确定

郭孜涵; 陈闯; 涂益良; 唐恩凌

doi:10.11858/gywlxb.20230817

HfZrTiTaNb系高熵合金的冲击反应释能定量确定

doi: 10.11858/gywlxb.20230817

沈阳理工大学辽宁省瞬态物理力学与能量转换材料重点实验室, 辽宁沈阳　110159

基金项目: 国家自然科学基金（12102272）

详细信息

作者简介:
郭孜涵（2001－），女，硕士研究生，主要从事冲击动力学研究. E-mail：guozihan1301@163.com

通讯作者:
陈　闯（1987－），男，博士，教授，主要从事爆炸与冲击动力学研究. E-mail：chenchuang517@126.com

中图分类号: O382
计量
- 文章访问数: 172
- HTML全文浏览量: 96
- PDF下载量: 45
出版历程
- 收稿日期: 2023-12-18
- 修回日期: 2024-01-02
- 网络出版日期: 2024-01-29
- 刊出日期: 2024-02-05

Quantitative Determination of Impact Reaction Energy Release for HfZrTiTaNb Based High-Entropy Alloys

Key Laboratory of Transient Physical Mechanics and Energy Conversion Materials of Liaoning Province, Shenyang Ligong University, Shenyang 110159, Liaoning, China

摘要

摘要: 作为一种新型含能材料，高熵合金在高速冲击过程中会释放大量能量，具有重要的应用价值。采用二级轻气炮系统加载真空环境下的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.
- high-entropy alloys /
- impact reaction energy release /
- energy distribution /
- unit mass release energy

HTML全文

高熵合金（high-entropy alloys，HEAs）具有优异的力学性能（高强度、高应变率强化效应以及强剪切自锐特性）和反应释能特性，被广泛应用于反应材料领域^[1–2]。高熵合金在高速撞击下破碎，诱发侵彻穿孔、高温气体、碎片云、火焰的传播与热传导等多种力学、热学、光学和化学现象并释放能量，该过程涉及复杂的物理化学变化，具有多阶段、多场耦合、多影响因素和高瞬态响应等特征。Zhang等^[3]采用电弧真空熔炼法制备了一种新型HfZrTiTa_0.53高熵合金，在高速撞击下，HfZrTiTa_0.53高能弹丸与空气发生反应释放出大量的能量。Geantă等^[4]采用口径为7.62 mm的燃烧穿甲子弹对4种不同化学成分的AlFeCrCoNi军用高熵合金试样进行了穿甲试验，样品表现出良好的韧性和力学性能，在高速撞击过程中产生了显著的侵爆破坏效应。王睿鑫^[5]通过弹道试验对NbZrTiTa高熵合金进行了反应释能特性评估，结果表明，在冲击过程中，靶室内的温度和压力急剧升高产生超压，随着侵彻速度的增加，弹丸在靶室内产生的超压增大，最高达0.18 MPa。Ren等^[6]研究了TiZrNbV高熵合金的压缩行为和冲击能释放特性，弹道试验后反应产物的化学反应程度和反应速率表明，高温剖面引起了严重的氧化反应，导致了冲击载荷下的化学能释放。Ma等^[7]研究了Ce含量对Al_0.5NbZrTi_1.5Ta_0.8动态压缩力学特性的影响，结果表明，含Ce的Al_0.5NbZrTi_1.5Ta_0.8高熵合金在弹道试验中可以穿透6 mm厚的钢板，并以719 m/s的速度点燃靶后方的棉花，其释能效果随着Ce含量的增加而增强。

高熵合金作为一种新型含能材料，相较于传统含能材料具有更好的释能效果。例如，在真空环境下，TiZrHfCu_0.3高熵合金弹丸冲击轴承钢靶板的释能能力为6.68 kJ/g^[8]，而Al/Teflon弹丸撞击钢板的释能能力为4.663 kJ/g^[9]。分析弹丸在冲击过程中的能量流向，定量评价其释能量对于高熵合金的工程应用具有重要意义。本研究按照固定的摩尔比制备HfZrTiTaNb系高熵合金，采用二级轻气炮系统加载真空环境中的HfZrTiTaNb系高熵合金弹丸，对GCr15轴承钢靶板进行冲击实验，利用红外热像仪、超压传感器以及瞬态光纤高温计测量闪光辐射温度、气体超压、火焰传播速度和容器壁温升等响应参数的演化过程，分析高熵合金弹丸撞击靶板过程中的能量流向，定量计算容器内混合气体的焓、闪光辐射能、准密闭容器壁吸收的能量、从入弹孔喷出气体的焓以及靶板的变形能，揭示不同元素及其含量对高熵合金释能量的影响。

1. 能量分布

高熵合金在高速冲击作用下会释放大量能量，其包含复杂的能量传递和转换过程。基于能量守恒定律，系统总能量由试件的动能E_v和试件与靶板撞击产生的化学能E_r构成。在高熵合金释能结束后，系统总能量主要转化为5种形式，分别为容器内混合气体的焓E_e（容器内气体的内能与压力功之和）、闪光辐射能E_f（物体在冲击条件下释放的可见光能量）、准密闭容器壁吸收的能量E_t、从入弹孔喷出气体的焓E_l（反应过程中从弹孔中喷出气体所携带的能量）以及靶板的变形能E_d（靶板受弹丸撞击后通过变形吸收的能量）。图1为准密闭容器实验系统中的能量分布。

图 1 准密闭容器实验系统中的能量分布

Figure 1. Energy distribution in a quasi-closed container experiment system

下载: 全尺寸图片幻灯片

弹丸的动能为

$E_{\text{v}}=\frac{1}{2}mv_{\mathrm{s}}^2$

(1)

式中：m为弹丸的质量，v_s为弹丸的撞击速度。

真空条件下HfZrTiTaNb系高熵合金冲击释能系统的能量分配如图2所示，其关系式可表示为

图 2 能量分配关系

Figure 2. Energy allocation relationship

下载: 全尺寸图片幻灯片

${E_{\text{r}}}+{E_{\text{v}}} = {E_{\text{e}}}+{E_{\text{f}}}+{E_{\text{t}}}+{E_{\text{l}}}+{E_{\text{d}}}$

(2)

2. 冲击反应释能实验

利用二级轻气炮加载HfZrTiTaNb系高熵合金弹丸，实验前通过真空泵将靶室空气抽出，使靶室环境接近真空状态。将准密闭容器置于靶室内部，容器进弹孔对准靶室入弹口，使得弹丸能够进入并撞击容器底部的GCr15轴承钢靶板。图3为高熵合金弹丸冲击反应释能实验加载与测试系统。

图 3 高熵合金弹丸冲击反应释能实验加载与测试系统

Figure 3. Experimental loading and testing system of high-entropy alloys projectile impact reaction energy release

下载: 全尺寸图片幻灯片

准密闭容器的内径为100 mm，长度为500 mm，壁厚为4 mm，容器两端用法兰盘密封，容器侧面加工一长方形视窗，尺寸为400 mm×150 mm，用于观测反应释能过程中的火焰传播。冲击反应释能实验系统由加载部分和测试设备组成，加载部分包括一级气室、二级气室、高压锥段、发射管和膨胀室等，测试设备包括瞬态光纤高温计、超压传感器、红外热像仪、高速摄像机和电磁测速系统等，通过脉冲触发器连接测试设备以实现同步触发。在冲击反应释能实验中，瞬态光纤高温计采集容器内的闪光信号，超压传感器记录容器内的气体超压信号，红外热像仪记录容器外壁温度的演化过程，高速摄像机记录火焰的流向，电磁测速系统测量高熵合金弹丸的撞击速度，示波器采集记录数据。图4为主要实验设备。

图 4 实验设备

Figure 4. Experimental equipment

下载: 全尺寸图片幻灯片

为了比较不同元素及其含量对高熵合金弹丸冲击反应释能的影响，开展了5组实验，实验基本参数列于表1。

表 1 实验基本参数

Table 1. Basic parameters of the experiment

No.	Specimen	Size/(mm×mm)	Mass/g	Impact velocity/(km·s⁻¹)	Density/(g·cm⁻³)
1	HfZrTiTaNbCu_0.2	$\varnothing$ 4.63×4.68	0.786	1.35	9.98
2	HfZrTiTaNbCu_0.8	$\varnothing$ 4.96×4.69	0.918	1.30	10.13
3	HfZrTiTaNbAl_0.2	$\varnothing$ 4.95×4.83	0.899	1.31	9.68
4	HfZrTiTaNbAl_0.8	$\varnothing$ 4.81×5.15	0.853	1.57	9.12
5	HfZrTiTaNbAl_0.8	$\varnothing$ 4.76×5.11	0.832	1.59	9.15

下载: 导出CSV

| 显示表格

3. 能量计算

3.1 容器内的混合气体焓

在测试闪光辐射温度之前，必须对瞬态光纤高温计进行标定。将光纤探头放置于太阳光模拟器的标准光源下，记录各通道的标定电压和其他标定结果。本研究中，使用400、500、600和700 nm的波长测量高熵合金弹丸释能过程中的闪光信号。瞬态光纤高温计各通道的标定数据如表2所示，其中N_r为光谱辐照度。

表 2 瞬态光纤高温计各通道的标定数据

Table 2. Calibration data for each channel of the transient fiber pyrometer

Channel	Wave length/nm	Impedance/kΩ	Voltage/V	N_r/(μW·cm⁻²·nm⁻¹)
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

下载: 导出CSV

| 显示表格

准密闭容器内气体热焓^[10]可以表示为

${E_{\text{e}}} = pV+{n_2}{c_V}\Delta {T_{\text{e}}}$

(3)

式中：p为容器内气体的超压，Pa；V为准密闭容器的体积，V=3.274×10⁻³ m³；c_V为气体平均比定容热容，其中c_V=3R，R=8.314 J/(mol·K)为理想气体常数；n₂为气体的物质的量，mol；ΔT_e为准密闭容器内气体的温度变化，K。

沿靶板法线方向距容器底部10 mm处放置瞬态高温计光纤探头，确保光纤探头能完整记录准密闭容器内的闪光辐射变化过程，采集到的光信号经过传感器转化为电信号，存储在示波器中。闪光辐射强度I_exp可表示为

${I_{\exp }} = \frac{{{h_{\exp }}}}{{{h_{\text{c}}}}} \frac{{l_0^2}}{{{l^2}}} \frac{{{N_{\text{r}}}(\lambda )}}{{2\text{π} (1 - \cos \theta )}}$

(4)

式中：I_exp为实验测量的闪光辐射强度，μW/(nm·cm²)；h_exp为示波器测量的电压幅值，V；h_c为实验前标定的电压值，V；l₀为标定距离，即实验前标定时光纤头部距离太阳光模拟器的距离，l₀=10 mm；l为光纤探头与靶板之间的距离，l=10 mm；N_r(λ)为标准光源给定的光谱辐照度，μW/(nm·cm²)；θ为光纤探头的孔径角，θ=37°。

根据普朗克热平衡辐射理论^[11]，实测闪光辐射强度与闪光辐射温度之间满足

$I\left( {\lambda ,T} \right) = \frac{1}{2}\varepsilon \tan^2{\left( {\frac{\theta }{2}} \right)}\frac{{{C_1}}}{{{\lambda ^5}}}\exp \left( { - \frac{{{C_2}}}{{\lambda T}}} \right)$

(5)

式中：I为理论计算的闪光辐射强度，μW/(nm·cm²)；T为闪光辐射温度，K；λ为波长，nm；C₁为第一辐射常数，C₁=3.7418×10⁻¹⁶ W·m²；C₂为第二辐射常数，C₂=1.4388×10⁻² m·K；ε为灰体表面发射率。

根据普朗克黑体辐射定律，采集4种不同波长（400、500、600和700 nm）的实验数据，计算瞬态闪光温度

$K(T) = \frac{{{I_1}{I_2}}}{{{I_3}{I_4}}} - {\left( {\frac{{{\lambda _2}{\lambda _4}}}{{{\lambda _1}{\lambda _3}}}} \right)^5} \exp \left[ {\frac{T}{{{C_2}}}\left( {{\lambda _2}+{\lambda _4} - {\lambda _1} - {\lambda _3}} \right)} \right]$

(6)

式中：I₁、I₂、I₃和I₄为不同波长对应的闪光辐射强度，μW/(nm·cm²)；λ₁、λ₂、λ₃和λ₄为每个通道对应的波长，nm。

图5显示了实验1～实验5的原始闪光数据。根据式(4)计算闪光辐射强度，辐射强度时程曲线如图6所示。

图 5 实验中光纤高温计测试的原始电压信号

Figure 5. Original voltage signals tested by the fiber pyrometer in the experiments

下载: 全尺寸图片幻灯片

图 6 实验中光纤高温计测试的闪光辐射强度时程曲线

Figure 6. Time-history curves of flash radiation intensity measured by fiber pyrometer in the experiments

下载: 全尺寸图片幻灯片

根据式(6)计算闪光辐射温度，图7显示了实验1～实验5准密闭容器内的气体超压和闪光辐射温度时程曲线。表3列出了实验中超压和闪光辐射温度的峰值。

图 7 准密闭容器中的超压和闪光辐射温度时程曲线

Figure 7. Time-history curves of overpressure and flash radiation temperature in the quasi-closed container

下载: 全尺寸图片幻灯片

表 3 实验中的超压和闪光辐射温度峰值

Table 3. Peak of overpressure and flash radiation temperature in the experiment

No.	Specimen	Overpressure/MPa	Temperature/K
1	HfZrTiTaNbCu_0.2	0.22	1260
2	HfZrTiTaNbCu_0.8	0.24	1370
3	HfZrTiTaNbAl_0.2	0.21	1120
4	HfZrTiTaNbAl_0.8	0.20	1228
5	HfZrTiTaNbAl_0.8	0.20	1281

下载: 导出CSV

| 显示表格

将测得的超压p、容器内的气体温度变化ΔT_e、准密闭容器体积V、平均比定容热容c_V代入式(3)，计算实验1～实验5中准密闭容器内的混合气体焓E_e，分别为761、834、729、698和701 J。受准密闭容器约束，向预留孔运动的反应产物在容器内壁发生反射并相互作用，因此，实验1～实验5的超压在达到峰值后振荡下降。在实验1、实验3和实验4中，弹丸的碎片部分在撞击靶板时发生反应，其闪光辐射温度达到最大峰值后还出现了多次较小峰值。

3.2 闪光辐射能

闪光辐射能量E_f^[12]可表示为

$E_{\text{f}}=\int_{t_1}^{t_2}\mathrm{d}t\int_{\lambda_1}^{\lambda_2}\frac{2hc^2}{\lambda^5}\frac{1}{\exp[hc/(\lambda k_{\mathrm{B}}T)]-1}\mathrm{d}\lambda$

(7)

式中：t为时间，s；c为光速，c=3×10⁸ m/s；h为普朗克常数，h=6.62×10⁻³⁴ J·s；k_B为玻尔兹曼常数，k_B=1.38×10⁻²³ J/K。

将式(7)的闪光辐射强度先对波长积分，再对闪光持续时间积分，得到闪光辐射能量。图8显示了实验1和实验5的闪光辐射能量曲线，可以看出，实验1和实验5的闪光辐射能量分别从1.08和1.10 ms开始上升，并在1.19和1.14 ms达到最大值，与瞬态光纤高温计的监测结果一致。实验1～实验5的闪光辐射能分别为0.52、0.57、0.46、0.51和0.50 J。

图 8 实验1和实验5闪光辐射能量时程曲线

Figure 8. Time history curves of flash radiation energy of experiment 1 and experiment 5

下载: 全尺寸图片幻灯片

3.3 准密闭容器壁吸收的能量

反应产物在热辐射、热对流和热传导的协同作用下加热容器壁，准密闭容器的内能增加，可以表示为^[12]

${E_{\text{t}}} = {c_3}{m_3}\Delta {T_{\text{t}}}$

(8)

式中：c₃为准密闭容器材料的比定压热容，c₃=0.46 kJ/(kg·K)；m₃为准密闭容器的质量，m₃=6.618 kg； $\Delta {T_{\text{t}}}$ 为准密闭容器的温升，K。

实验中，由红外摄像机记录准密闭容器的外壁温度变化情况，通过配套软件分析得到实际观测区域的温度变化。图9显示了准密闭容器的实际观测区域和红外视图。

图 9 准密闭容器的观测区域和红外视图

Figure 9. Observation area and infrared view of quasi-closed container

下载: 全尺寸图片幻灯片

图10(a)显示了准密闭容器的外壁温度时程曲线。弹丸冲击靶板后，观测区域吸收冲击反应产生的热量，导致准密闭容器的外壁温度急剧上升。实验1～实验5中，实际观测区域的温度分别上升0.35369、0.50296、0.25286、0.45865和0.46616 K。根据式(8)计算实验1～实验5中准密闭容器壁吸收的能量E_t，分别为1077、1531、769、1396和1419 J。图10(b)显示了t=2 s时实验1～实验5中准密闭容器壁的红外视图。

图 10 准密闭容器的外壁温度时程曲线和典型时刻的红外视图

Figure 10. Temperature time history curve of the outer wall of quasi-closed container and the infrared view at typical moments

下载: 全尺寸图片幻灯片

3.4 反应产物经进弹孔消耗的能量

图11～图15显示了准密闭容器内HfZrTiTaNb系高熵合金弹丸高速冲击轴承钢靶板的反应过程。

图 11 准密闭容器内的反应产物流动（实验1：HfZrTiTaNbCu_0.2弹丸）

Figure 11. Flow of reaction products in the quasi-closed container (Experiment 1: HfZrTiTaNbCu_0.2 projectile)

下载: 全尺寸图片幻灯片

图 12 准密闭容器内的反应产物流动（实验2：HfZrTiTaNbCu_0.8弹丸）

Figure 12. Flow of reaction products in the quasi-closed container (Experiment 2: HfZrTiTaNbCu_0.8 projectile)

下载: 全尺寸图片幻灯片

图 13 准密闭容器内的反应产物流动（实验3：HfZrTiTaNbAl_0.2弹丸）

Figure 13. Flow of reaction products in the quasi-closed container (Experiment 3: HfZrTiTaNbAl_0.2 projectile)

下载: 全尺寸图片幻灯片

图 14 准密闭容器内的反应产物流动（实验4：HfZrTiTaNbAl_0.8弹丸）

Figure 14. Flow of reaction products in the quasi-closed container (Experiment 4: HfZrTiTaNbAl_0.8 projectile)

下载: 全尺寸图片幻灯片

图 15 准密闭容器内的反应产物流动（实验5：HfZrTiTaNbAl_0.8弹丸）

Figure 15. Flow of reaction products in the quasi-closed container (Experiment 5: HfZrTiTaNbAl_0.8 projectile)

下载: 全尺寸图片幻灯片

在撞击靶板过程中弹丸反应形成的高温气体产物从形成弹坑位置以不同的速度向准密闭容器的预留孔流动，并从预留孔中逸出。忽略摩擦和热交换，流动参数在垂直截面上均匀分布，该流动可简化为一维定常等熵流动。

在一维定常等熵流动中，流场的控制方程组为^[13]

$\left\{\begin{gathered}\mathrm{d}m=A\mathrm{d}\rho\mathrm{d}u=\mathrm{const} \\ \mathrm{d}p=\rho u\mathrm{d}u \\ \mathrm{d}(h_{\text{j}}+0.5u^2)=0 \\ \end{gathered}\right.$

(9)

式中：A为弹孔的横截面积，m²；ρ为通过的气体密度，g/cm³；u为通过的气体流速，m/s；h_j为气体产物的焓，J。

反应产物经进弹孔实时喷射到准密闭容器外的质量流出速率为^[14]

$\frac{\mathrm{d}m_1}{\mathrm{d}t}=A\sqrt{\frac{\gamma m_1p_{\text{j}}}{V}}\left(\frac{2}{\gamma+1}\right)^{\frac{\gamma+1}{2(\gamma-1)}}$

(10)

式中：γ为绝热指数，γ=1.4；m₁为预留孔喷射气体的质量，g；p_j为压力，MPa；t为时间，s。

反应气体产物经进弹孔丢失的能量可表示为

${E_{\text{l}}} = n{c_p}\Delta T$

(11)

式中：n为喷出气体的物质的量，mol；c_p=4R为比定压热容。

根据式(11)计算实验1～实验5中反应气体产物经进弹孔丢失的能量E_l，分别为58.67、63.79、55.87、57.18和59.64 J。

采用Fluent软件模拟准密闭容器内实验1的反应气体流动过程^[8]，为提高计算效率，选择二维轴对称模型，加密入口和通道处的网格，网格模型如图16所示。将实验的气体运动速度和温度作为数值模拟的边界条件，求解弹孔溢出气体的质量时程曲线和准密闭容器内气体的温度时程曲线。图17显示了不同时刻准密闭容器内的数值模拟温度云图。

图 16 数值模拟网格模型

Figure 16. Grid model of numerical simulation

下载: 全尺寸图片幻灯片

图 17 不同时刻准密闭容器内的数值模拟温度云图

Figure 17. Numerical simulation of temperature cloud in the quasi-closed container at different moments

下载: 全尺寸图片幻灯片

图18(a)显示了实验中距容器底部10 mm处的数值模拟温度时程曲线，模拟的气体温度峰值（1252 K）比实验结果小8 K。模拟的从准密闭容器弹孔喷出的气体产物能量为50.88 J（见图18(b)），比实验结果低13.2%。

图 18 容器内气体温度及喷出气体质量的时程曲线

Figure 18. Time history curve of the gas temperature and the ejected gas mass in the container

下载: 全尺寸图片幻灯片

3.5 靶板变形能

冲击作用下，准密闭容器底部的靶板形成外围呈“翻唇”形的半球形弹坑，弹坑周围的靶板受应力波作用形成应变区。弹丸撞击靶板形成的弹坑及其剖面如图19所示。

图 19 弹丸撞击靶板形成的弹坑及其剖面示意图

Figure 19. Craters formed by the impact of projectile on the target plate and their cross-sections

下载: 全尺寸图片幻灯片

靶板材料的变形能包含体积应变能E_V和剪切应变能E_τ^[15]

${E}_{\text{d}}={E}_{\text{V}}+{E}_{\tau }=2\text{π} {{{\displaystyle \int }}^{\text{}}}_{ 0}^{+\infty }\sigma {\varepsilon }_{\text{V}}{r}^{2}\text{d}r+2\text{π} {{{\displaystyle \int }}^{\text{}}}_{ 0}^{+\infty }\tau \gamma {r}^{2}\text{d}r$

(12)

式中：σ为压缩应力，τ为剪切应力，ε_V为体积应变，γ为剪切应变。E_d由弹坑内部的能量E₁和孔洞外部的能量E₂组成

${E}_{1}=2\text{π} {{{\displaystyle \int }}^{\text{}}}_{ 0}^{{r}_{0}}\sigma {\varepsilon }_{\text{V}}{r}^{2}\text{d}r+2\text{π} {{{\displaystyle \int }}^{\text{}}}_{ 0}^{{r}_{0}}\tau \gamma {r}^{2}\text{d}r$

(13)

${E}_{2}=2\text{π} {\displaystyle {\int }_{ {r}_{0}}^{+\infty }\sigma {\varepsilon }_{\text{V}}{r}^{2}\text{d}r}+2\text{π} {{{\displaystyle \int }}^{\text{}}}_{ {r}_{0}}^{+\infty }\tau \gamma {r}^{2}\text{d}r$

(14)

弹丸侵彻靶板时，撞击点处的部分靶板熔化或者汽化，该过程产生的能量主要以闪光辐射的形式发散出去。弹坑内部吸收的能量由相应冲击条件下的破坏应力和温度决定，因此，弹坑内部能量E₁为

${E}_{1}=2\text{π} {{{\displaystyle \int }}^{\text{}}}_{ 0}^{{r}_{0}}{\sigma }_{\text{m}}{\varepsilon }_{\text{m}}{r}^{2}\text{d}r+2\text{π} {{{\displaystyle \int }}^{\text{}}}_{ 0}^{{r}_{0}}{\tau }_{\text{m}}{\gamma }_{\text{m}}{r}^{2}\text{d}r=2\text{π} {{{\displaystyle \int }}^{\text{}}}_{ 0}^{{r}_{0}}\left({\sigma }_{\text{m}}{\varepsilon }_{\text{m}}+{\tau }_{\text{m}}{\gamma }_{\text{m}}\right){r}^{2}\text{d}r$

(15)

式中：σ_m为破坏压缩应力，τ_m为剪切破坏应力，ε_m为破坏压缩体积应变，γ_m为剪切破坏应变。

将 ${\sigma _{\text{m}}}{\varepsilon _{\text{m}}}+{\tau _{\text{m}}}{\gamma _{\text{m}}} = {\sigma _{\text{u}}}{\varepsilon _{\text{u}}}$ 代入式(15)，可得

${E}_{1}=2\text{π} {{{\displaystyle \int }}^{\text{}}}_{ 0}^{{r}_{0}}{\sigma }_{\text{u}}{\varepsilon }_{\text{u}}{r}^{2}\text{d}r=\frac{2}{3}\text{π} {r}_{0}^{3}{\sigma }_{\text{u}}{\varepsilon }_{\text{u}}$

(16)

式中：σ_u为破坏压缩应力，ε_u为破坏应变。

弹坑外部能量E₂为

${E_2} = 2\text{π} \int_{{r_0}}^{+\infty} {\left( {\sigma {\varepsilon _{\text{V}}}+\tau \gamma } \right){r^2}{\text{d}}r} = 2\text{π} \int_{{r_0}}^{+\infty} {\sigma \varepsilon {r^2}{\text{d}}r}$

(17)

根据靶材硬度与应变之间的关系，可得应变区到孔洞中心的距离（r）与形变（ε）之间的关系

$\varepsilon {\text{ = }}{{\text{e}}^{ - {{B}}\left( {r - {{C}}} \right)}}$

(18)

式中：B和C为常数。当r=r₀（r₀为弹坑内部半径）时，ε=ε_u。所以

$\varepsilon_{\text{u}}\text{ = }\text{e}^{\mathit{\mathrm{-}{B}}\left(r_0-\sigma\right)}$

(19)

将 $\sigma {\text{ = }}{\sigma _{\text{u}}}r_0^2/{r^2}$ 和式(19)代入式(17)，得到弹坑外部的能量E₂

${E}_{2}=2\text{π} {{{\displaystyle \int }}^{\text{}}}_{ 0}^{+\infty }\sigma \varepsilon {r}^{2}\text{d}r=2\text{π} {{{\displaystyle \int }}^{\text{}}}_{ 0}^{+\infty }{\sigma }_{\text{u}}\frac{{r}_{0}^{2}}{{r}^{2}}{\text{e}}^{-{B}(r-{C})}{r}^{2}\text{d}r=2\text{π} {\sigma }_{\text{u}}{\varepsilon }_{\text{u}}\frac{{r}_{0}^{2}}{{\rho }_{\text{t}}}$

(20)

式中：ρ_t为靶板材料的密度，ρ_t=7.8×10³ kg/m³；σ_u为轴承钢靶板的破坏应力，σ_u=2.2 GPa；ε_u=0.1^[16]。实验1～实验5中，弹丸侵彻靶板形成的弹坑直径分别为13.16、12.05、11.80、11.24和10.92 mm。根据式(20)计算实验1～实验5中的靶板变形能，分别为138、107、101、86和80 J。

3.6 冲击反应释能总量

将E_v、E_e、E_f、E_t、E_l、E_d代入式(2)计算高熵合金的冲击反应释能E_r，再将E_r除以试件质量得到单位质量HfZrTiTaNb系高熵合金的反应释能量。从表4可以看出，随着Cu或Al含量的增加，HfZrTiTaNb系高熵合金的单位质量释能量增大。对比实验2、实验4和实验5，在相近的撞击速度下，含Cu高熵合金的单位质量释能量比含Al高熵合金大。

表 4 冲击反应释能结果

Table 4. Impact reaction energy release results

No.	Specimen	Impact velocity/ (km·s⁻¹)	Mass/g	Energy release/kJ	Energy release per unit mass/(kJ·g⁻¹)
1	HfZrTiTaNbCu_0.2	1.35	0.786	1.32	1.68
2	HfZrTiTaNbCu_0.8	1.30	0.918	1.76	1.92
3	HfZrTiTaNbAl_0.2	1.31	0.899	0.88	0.98
4	HfZrTiTaNbAl_0.8	1.57	0.853	1.19	1.40
5	HfZrTiTaNbAl_0.8	1.59	0.832	1.21	1.45

下载: 导出CSV

| 显示表格

3.7 释能分析

Cu和Al的原子尺寸不同，如图20所示，原子尺寸会影响合金晶格畸变的程度。

图 20 不同原子的半径

Figure 20. Radius of different atoms

下载: 全尺寸图片幻灯片

由于Cu的原子半径小于Al，当Cu被引入合金中时，会引起更大的合金族原子尺寸差异。原子尺寸差（δ）可表示为^[17]

$\delta = \sqrt {\sum\limits_{i = 1}^n {{c_i}\left(1 - {r_i}\Big/\sum\limits_{i = 1}^n {c_i}{r_i}\right)^2} }$

(21)

式中：c_i为第i类原子在合金内的原子比，r_i为第i类原子的原子半径。

HfZrTiTaNb系高熵合金的原子尺寸差列于表5。在合金元素摩尔比（x）相同的情况下，含Cu高熵合金的原子尺寸差比含Al的大。Cu元素的添加使合金更易产生晶格畸变和相变，在冲击载荷下合金更易破裂，从而释放更多的能量。

表 5 HfZrTiTaNb系高熵合金的原子尺寸差

Table 5. Atomic size difference of HfZrTiTaNb based high-entropy alloys

x	δ/%
x	HfZrTiTaNbCu_x	HfZrTiTaNbAl_x
0.2	4.92	4.84
0.8	6.73	6.37

下载: 导出CSV

| 显示表格

4. 结　论

基于二级轻气炮实验加载系统，通过超压传感器、瞬态光纤高温计、红外热像仪以及高速摄像机等实验设备，研究了真空环境下HfZrTiTaNb系高熵合金弹丸在撞击靶板过程中引发的力学、热学、光学和化学效应，计算了单位质量高熵合金冲击反应的释能量，得到以下主要结论。

(1) 在高熵合金弹丸撞击靶板过程中，系统释放的能量表现为靶板变形能、闪光辐射能、容器内的混合气体焓、准密闭容器壁吸收的能量和反应产物经进弹孔丢失的能量，其中准密闭容器壁吸收的能量占比最高。

(2) 随着Cu或Al含量的增加，HfZrTiTaNb系高熵合金的单位质量释能量增加。在相近的撞击速度下，含Cu高熵合金的单位质量释能量比含Al高熵合金高。

图 1 准密闭容器实验系统中的能量分布

Figure 1. Energy distribution in a quasi-closed container experiment system

下载: 全尺寸图片幻灯片

图 2 能量分配关系

Figure 2. Energy allocation relationship

下载: 全尺寸图片幻灯片

图 3 高熵合金弹丸冲击反应释能实验加载与测试系统

Figure 3. Experimental loading and testing system of high-entropy alloys projectile impact reaction energy release

下载: 全尺寸图片幻灯片

图 4 实验设备

Figure 4. Experimental equipment

下载: 全尺寸图片幻灯片

图 5 实验中光纤高温计测试的原始电压信号

Figure 5. Original voltage signals tested by the fiber pyrometer in the experiments

下载: 全尺寸图片幻灯片

图 6 实验中光纤高温计测试的闪光辐射强度时程曲线

Figure 6. Time-history curves of flash radiation intensity measured by fiber pyrometer in the experiments

下载: 全尺寸图片幻灯片

图 7 准密闭容器中的超压和闪光辐射温度时程曲线

Figure 7. Time-history curves of overpressure and flash radiation temperature in the quasi-closed container

下载: 全尺寸图片幻灯片

图 8 实验1和实验5闪光辐射能量时程曲线

Figure 8. Time history curves of flash radiation energy of experiment 1 and experiment 5

下载: 全尺寸图片幻灯片

图 9 准密闭容器的观测区域和红外视图

Figure 9. Observation area and infrared view of quasi-closed container

下载: 全尺寸图片幻灯片

图 10 准密闭容器的外壁温度时程曲线和典型时刻的红外视图

Figure 10. Temperature time history curve of the outer wall of quasi-closed container and the infrared view at typical moments

下载: 全尺寸图片幻灯片

图 11 准密闭容器内的反应产物流动（实验1：HfZrTiTaNbCu_0.2弹丸）

Figure 11. Flow of reaction products in the quasi-closed container (Experiment 1: HfZrTiTaNbCu_0.2 projectile)

下载: 全尺寸图片幻灯片

图 12 准密闭容器内的反应产物流动（实验2：HfZrTiTaNbCu_0.8弹丸）

Figure 12. Flow of reaction products in the quasi-closed container (Experiment 2: HfZrTiTaNbCu_0.8 projectile)

下载: 全尺寸图片幻灯片

图 13 准密闭容器内的反应产物流动（实验3：HfZrTiTaNbAl_0.2弹丸）

Figure 13. Flow of reaction products in the quasi-closed container (Experiment 3: HfZrTiTaNbAl_0.2 projectile)

下载: 全尺寸图片幻灯片

图 14 准密闭容器内的反应产物流动（实验4：HfZrTiTaNbAl_0.8弹丸）

Figure 14. Flow of reaction products in the quasi-closed container (Experiment 4: HfZrTiTaNbAl_0.8 projectile)

下载: 全尺寸图片幻灯片

图 15 准密闭容器内的反应产物流动（实验5：HfZrTiTaNbAl_0.8弹丸）

Figure 15. Flow of reaction products in the quasi-closed container (Experiment 5: HfZrTiTaNbAl_0.8 projectile)

下载: 全尺寸图片幻灯片

图 16 数值模拟网格模型

Figure 16. Grid model of numerical simulation

下载: 全尺寸图片幻灯片

图 17 不同时刻准密闭容器内的数值模拟温度云图

Figure 17. Numerical simulation of temperature cloud in the quasi-closed container at different moments

下载: 全尺寸图片幻灯片

图 18 容器内气体温度及喷出气体质量的时程曲线

Figure 18. Time history curve of the gas temperature and the ejected gas mass in the container

下载: 全尺寸图片幻灯片

图 19 弹丸撞击靶板形成的弹坑及其剖面示意图

Figure 19. Craters formed by the impact of projectile on the target plate and their cross-sections

下载: 全尺寸图片幻灯片

图 20 不同原子的半径

Figure 20. Radius of different atoms

下载: 全尺寸图片幻灯片

表 1 实验基本参数

Table 1. Basic parameters of the experiment

No.	Specimen	Size/(mm×mm)	Mass/g	Impact velocity/(km·s⁻¹)	Density/(g·cm⁻³)
1	HfZrTiTaNbCu_0.2	$\varnothing$ 4.63×4.68	0.786	1.35	9.98
2	HfZrTiTaNbCu_0.8	$\varnothing$ 4.96×4.69	0.918	1.30	10.13
3	HfZrTiTaNbAl_0.2	$\varnothing$ 4.95×4.83	0.899	1.31	9.68
4	HfZrTiTaNbAl_0.8	$\varnothing$ 4.81×5.15	0.853	1.57	9.12
5	HfZrTiTaNbAl_0.8	$\varnothing$ 4.76×5.11	0.832	1.59	9.15

下载: 导出CSV

表 2 瞬态光纤高温计各通道的标定数据

Table 2. Calibration data for each channel of the transient fiber pyrometer

Channel	Wave length/nm	Impedance/kΩ	Voltage/V	N_r/(μW·cm⁻²·nm⁻¹)
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

下载: 导出CSV

表 3 实验中的超压和闪光辐射温度峰值

Table 3. Peak of overpressure and flash radiation temperature in the experiment

No.	Specimen	Overpressure/MPa	Temperature/K
1	HfZrTiTaNbCu_0.2	0.22	1260
2	HfZrTiTaNbCu_0.8	0.24	1370
3	HfZrTiTaNbAl_0.2	0.21	1120
4	HfZrTiTaNbAl_0.8	0.20	1228
5	HfZrTiTaNbAl_0.8	0.20	1281

下载: 导出CSV

表 4 冲击反应释能结果

Table 4. Impact reaction energy release results

No.	Specimen	Impact velocity/ (km·s⁻¹)	Mass/g	Energy release/kJ	Energy release per unit mass/(kJ·g⁻¹)
1	HfZrTiTaNbCu_0.2	1.35	0.786	1.32	1.68
2	HfZrTiTaNbCu_0.8	1.30	0.918	1.76	1.92
3	HfZrTiTaNbAl_0.2	1.31	0.899	0.88	0.98
4	HfZrTiTaNbAl_0.8	1.57	0.853	1.19	1.40
5	HfZrTiTaNbAl_0.8	1.59	0.832	1.21	1.45

下载: 导出CSV

表 5 HfZrTiTaNb系高熵合金的原子尺寸差

Table 5. Atomic size difference of HfZrTiTaNb based high-entropy alloys

x	δ/%
x	HfZrTiTaNbCu_x	HfZrTiTaNbAl_x
0.2	4.92	4.84
0.8	6.73	6.37

下载: 导出CSV

参考文献(17)

[1]	刘丘林, 刘允中, 王艳群. 高熵合金的研究现状和应用前景 [J]. 粉末冶金工业, 2017, 27(6): 64–69. doi: 10.13228/j.boyuan.issn1006-6543.20160115 LIU 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 HfZrTiTa_0.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 TiZrHfX_0.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.005 TANG 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.008 DUAN 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_xCoCrFeNi高熵合金力学性能的分子动力学模拟 [J]. 高压物理学报, 2021, 35(5): 052201. doi: 10.11858/gywlxb.20210730 ZHANG L M, MA S G, LI Z Q, et al. Mechanical properties of Al_xCoCrFeNi high-entropy alloy: a molecular dynamics study [J]. Chinese Journal of High Pressure Physics, 2021, 35(5): 052201. doi: 10.11858/gywlxb.20210730

施引文献

资源附件(0)

访问统计

点击查看大图

图(20) / 表(5)

计量

文章访问数: 172
HTML全文浏览量: 96
PDF下载量: 45

1. 能量分布
2. 冲击反应释能实验
3. 能量计算
3.1 容器内的混合气体焓
3.2 闪光辐射能
3.3 准密闭容器壁吸收的能量
3.4 反应产物经进弹孔消耗的能量
3.5 靶板变形能
3.6 冲击反应释能总量
3.7 释能分析
4. 结　论

HfZrTiTaNb系高熵合金的冲击反应释能定量确定

doi: 10.11858/gywlxb.20230817

作者简介: 郭孜涵（2001－），女，硕士研究生，主要从事冲击动力学研究. E-mail：guozihan1301@163.com

通讯作者: 陈 闯（1987－），男，博士，教授，主要从事爆炸与冲击动力学研究. E-mail：chenchuang517@126.com

计量

出版历程

Quantitative Determination of Impact Reaction Energy Release for HfZrTiTaNb Based High-Entropy Alloys

1. 能量分布

2. 冲击反应释能实验

3. 能量计算

3.1 容器内的混合气体焓

3.2 闪光辐射能

3.3 准密闭容器壁吸收的能量

3.4 反应产物经进弹孔消耗的能量

3.5 靶板变形能

3.6 冲击反应释能总量

3.7 释能分析

4. 结 论

计量

出版历程

目录

1. 能量分布

2. 冲击反应释能实验

3. 能量计算

3.1 容器内的混合气体焓

3.2 闪光辐射能

3.3 准密闭容器壁吸收的能量

3.4 反应产物经进弹孔消耗的能量

3.5 靶板变形能

3.6 冲击反应释能总量

3.7 释能分析

4. 结 论

作者简介:
郭孜涵（2001－），女，硕士研究生，主要从事冲击动力学研究. E-mail：guozihan1301@163.com

通讯作者:
陈　闯（1987－），男，博士，教授，主要从事爆炸与冲击动力学研究. E-mail：chenchuang517@126.com

4. 结　论

4. 结　论