UHMWPE的应变率效应及其对超高速碰撞特性的影响

石景富; 于东; 徐铧东; 刘蕾; 苗常青

doi:10.11858/gywlxb.20220666

UHMWPE的应变率效应及其对超高速碰撞特性的影响

doi: 10.11858/gywlxb.20220666

石景富^1,,
于东²,
徐铧东¹,
刘蕾¹,
苗常青^1,*, ,

1.
哈尔滨工业大学特种环境复合材料技术国家级重点实验室, 黑龙江哈尔滨　150001
2.
哈尔滨工业大学机电学院机械设计系, 黑龙江哈尔滨　150001

详细信息

作者简介:
石景富（1994－），男，博士研究生，主要从事聚合物材料力学性能研究. E-mail：shijf@stu.hit.edu.cn

通讯作者:
苗常青（1972－），男，博士，教授，主要从事柔性复合材料研究. E-mail：miaocq@hit.edu.cn

中图分类号: O347.3; O414.19
计量
- 文章访问数: 281
- HTML全文浏览量: 126
- PDF下载量: 41
出版历程
- 收稿日期: 2022-09-28
- 修回日期: 2022-11-04
- 网络出版日期: 2023-04-17
- 刊出日期: 2023-06-05

Strain Rate Effect of UHMWPE and Its Influence on Hypervelocity Impact Performance

SHI Jingfu^1
,,
YU Dong²,
XU Huadong¹,
LIU Lei¹,
MIAO Changqing^{1,*
, ,}

1.
National Key Laboratory of Science and Technology on Advanced Composites in Special Environments, Harbin Institute of Technology, Harbin 150001, Heilongjiang, China
2.
Department of Mechanical Design, School of Mechanical and Electrical Engineering, Harbin Institute of Technology, Harbin 150001, Heilongjiang, China

摘要

摘要: 为分析超高分子量聚乙烯（ultra-high molecular weight polyethylene, UHMWPE）的应变率效应及其对超高速碰撞特性的影响规律，采用万能材料试验机和分离式霍普金森拉杆对UHMWPE纤维束进行静、动态拉伸实验，获得了不同应变率下材料的应力-应变关系，并进一步开展了UHMWPE纤维织物的超高速碰撞数值模拟。结果表明，UHMWPE的拉伸模量和强度均随应变率的升高而逐渐增大。随着材料应变率敏感系数的增大，防护结构对弹丸动能的吸收率呈现先减小后增大的趋势。
- 超高分子量聚乙烯 /
- 应变率 /
- 动态拉伸 /
- 超高速碰撞
Abstract: To analyze the strain rate effect and its influence on the dynamic response of ultra-high molecular weight polyethylene (UHMWPE) subjected to hypervelocity impact, the stress-strain relation of UHMWPE at different strain rates loading were obtained through universal material testing machine and Hopkinson bar experimental system. Furthermore, strain rate effect on the hypervelocity impact response was analyzed by numerical simulation. The results show that the tensile modulus and strength increase with the rise of strain rate. With the increase of the strain rate sensitivity index, the energy absorption ratio of the protective structure against the projectile decreases at the beginning followed by an increasing trend.
- ultra-high molecular weight polyethylene (UHMWPE) /
- strain rate /
- dynamic tensile /
- hypervelocity impact

HTML全文

弹药在储存、运输以及使用过程中，一旦遭受意外热刺激，则易发生燃烧、爆炸甚至爆轰，因此有必要对炸药的热安全性进行合理评估。在确保弹药安全的前提下提高炸药威力是目前亟待解决的问题。在混合炸药配方中添加氧化剂可改善其氧平衡，较大幅度地提升爆热。高氯酸铵（NH₄ClO₄，ammonium perchlorate，AP）具有密度较高、含氧量高、爆轰产物为气体、机械感度低等优点，是混合炸药配方的重要候选组分^[1]。研究含有AP的浇注炸药的热响应特性对其推广应用以及新型钝感高聚物黏结炸药（plastic bonded explosives，PBX）研发具有重要意义。

国内外对含AP的含能材料热分解特性以及热安全性开展了较多研究。刘子如等^[2]研究了AP-奥克托今（HMX）混合物的热分解特性，发现AP的存在使HMX的微商热重（derivative thermogravimetry，DTG）峰温大幅下降，证明AP与HMX在热分解中存在强烈的相互作用。Li等^[3]研究了铝粉的加入对AP-HMX混合物热分解特性的影响，发现足量的铝粉可以消除AP与HMX之间的相互催化作用。陈中娥等^[4]认为，AP热分解形成的孔隙是AP-HTPB推进剂慢速烤燃响应剧烈的原因。李苗苗等^[5]对不同AP含量的HTPE推进剂进行了差示扫描量热法测试及慢速烤燃实验，发现随着AP含量的降低，推进剂的热稳定性和热安全性略有提升。Kim等^[6]对含有AP的PBX炸药进行了快速和慢速烤燃实验，结果表明，在慢速烤燃中AP不参与点火，而在快速烤燃中点火由AP与硝胺类炸药共同参与。由此可见，AP在热刺激下的自热反应受到添加剂和温升速率等因素的影响，但是目前仍缺乏含AP的PBX炸药的热刺激响应特性相关研究^[7]。

在炸药的热刺激响应特性实验中，温升速率对炸药的点火时间、点火位置以及装药的温度梯度等产生很大的影响。王沛等^[8]对GHL炸药进行了慢速烤燃实验以及4种不同温升速率下的数值模拟，发现点火位置随着温升速率的增大从炸药内部移向炸药边缘。代晓淦等^[9]进行了不同温升速率下PBX-2炸药的烤燃实验和数值模拟，发现温升速率越大，炸药中心温度在烤燃过程中的变化越小。刘静等^[10]对模块装药进行了慢速烤燃实验和数值模拟，发现随着温升速率的提高，烤燃响应区域由1个中心环形区域变成关于中心对称的2个环形区域。由此可见，温升速率对炸药烤燃中的点火位置以及点火区域形状的影响较大，对获取热刺激下炸药的响应特性具有重要意义。

本研究将开展新型浇注炸药GOLA-1的慢速烤燃实验和数值模拟，分析GOLA-1炸药在不同条件下的热刺激响应特性及其在慢速烤燃中点火的主导因素，探讨温升速率对小型烤燃弹点火位置的影响。

1. 实验与数值模拟方法

1.1 小型烤燃弹实验

基于小型烤燃弹实验装置对GOLA-1炸药进行慢速烤燃实验。图1为小型烤燃弹实验测试系统示意图。该装置主要包括壳体、端盖、炸药、环形加热器、测温热电偶以及控温系统^[11]。烤燃弹的外径为58 mm，长度为106 mm，壳体厚度为4 mm，药柱的直径为50 mm，长度为88 mm。药柱与端盖之间有高12 mm的空气域，用以缓解炸药受热膨胀对壳体的压力。端盖中心以及炸药上表面至炸药中心均开有直径为1 mm的小孔，用于放置K型热电偶，监测温升过程中炸药中心点的温度变化。壳体和端盖材质均为45钢，两者采用螺纹连接。装药壳体外侧包裹环形加热器，在慢速烤燃实验中温升速率为1.0 和1.5 K/min。加热器外部使用岩棉包裹进行保温，防止热量散失影响实验结果。在壳体与加热器之间放置两根K型热电偶，用来控制加热速率以及测量壳体外壁面的温度变化。

图 1 小型烤燃弹实验装置示意图

Figure 1. Schematic diagram of small-scale cook-off bomb experiment device

下载: 全尺寸图片幻灯片

1.2 烤燃数值模拟

烤燃实验只能获得温度监测点的温度-时间关系，无法获取炸药内部的温度分布、点火位置、点火温度以及各组分的质量分数等信息。数值模拟可以较方便地获取以上信息，从而对炸药的热刺激响应特性进行较全面的评估。基于小型烤燃弹实验装置，建立如图2所示的数值模拟模型，其中端盖及壳体、空气域以及炸药分别用粉色、蓝色和红色表示。模型用六面体网格划分，网格尺寸为0.5 mm。

图 2 小型烤燃弹的计算网格模型

Figure 2. Computational grid model of small-scale cook-off bomb

下载: 全尺寸图片幻灯片

在直角坐标系下，热传导方程的形式为^[12]

$\rho {c_p}\frac{{\partial T}}{{\partial t}} = \lambda {\nabla ^2}T+S$

(1)

式中：ρ为密度，kg·m⁻³；c_p为比定压热容，J·kg⁻¹·K⁻¹；T为温度，K；λ为热导率，W·m⁻¹·K⁻¹；S为炸药的自热反应源项，W·m⁻³。

空气域与炸药之间的辐射热流q_r采用P1辐射模型描述^[13]

${q_{\rm r}} = \frac{1}{{3\left( {a+{\sigma _{\rm s}}} \right) - C{\sigma _{\rm s}}}}\nabla G$

(2)

式中：a为吸收系数，m⁻¹；σ_s为散射系数，m⁻¹；C为线性各向异性相位函数系数；G 为入射辐射。

GOLA-1为HMX基浇注炸药，主要成分为HMX，其他成分包括Al、 AP、黏结剂（Kel-F）。组分HMX、Al、AP、Kel-F的质量比为44∶25∶20∶11，密度为1.837 g/cm³。对于炸药中的主要成分HMX，在慢速烤燃过程中，随着温度上升，HMX首先由β相转变为δ相，而后热分解生成最终的气态产物。基于过渡态理论，HMX在烤燃过程中的4步化学反应动力学模型^[14]表示为：(1) β-HMX ↔ δ-HMX（反应1，单分子吸热反应）；(2) β-HMX+δ-HMX ↔ δ-HMX（反应2，双分子吸热反应）；(3) δ-HMX→产物（反应3，单分子吸热反应）；(4) δ-HMX+产物→产物（反应4，双分子放热反应）。反应1和反应2描述β相HMX转变为δ相HMX的晶型转变过程，反应3和反应4描述最终产物生成的点火过程。各步反应速率方程为

${r_1} = {w_\beta }\frac{{kT}}{h}\exp \left( {\frac{{T\Delta {S_{1{\text{f}}}} - {E_{1{\text{f}}}}}}{{RT}}} \right) - {w_\delta }\frac{{kT}}{h}\exp \left( {\frac{{T\Delta {S_{1{\text{r}}}} - {E_{1{\text{r}}}}}}{{RT}}} \right)$

(3)

${r_2} = {w_\beta }{w_\delta }\frac{{kT}}{h}\left[ {\exp \left( {\frac{{T\Delta {S_{2{\text{f}}}} - {E_{2{\text{f}}}}}}{{RT}}} \right) - \exp \left( {\frac{{T\Delta {S_{2{\text{r}}}} - {E_{2{\text{r}}}}}}{{RT}}} \right)} \right]$

(4)

${r_3} = {w_\delta }{Z_3}\exp \left( { - \frac{{{E_3}}}{{RT}}} \right)$

(5)

${r_4} = {w_\delta }{w_{\rm p} }{Z_4}\exp \left( { - \frac{{{E_4}}}{{RT}}} \right)$

(6)

式中：r₁～r₄为反应1～反应4的反应速率；w_β和w_δ分别为β相和δ相HMX的质量分数，w_p为产物的质量分数；k为玻尔兹曼常数；h为普朗克常数；R为普适气体常数；ΔS为活化熵，J∙mol⁻¹∙K⁻¹；E为活化能，J∙mol⁻¹；下标f代表正反应，r代表逆反应；Z₁～Z₄为反应1～反应4的指前因子，s⁻¹。

各组分质量分数的微分方程为

$\frac{{{\text{d}}{w_\beta }}}{{{\text{d}}t}} = - {r_1} - {r_2}$

(7)

$\frac{{{\text{d}}{w_\delta }}}{{{\text{d}}t}} = {r_1}+{r_2} - {r_3} - {r_4}$

(8)

$\frac{{{\text{d}}{w_{\text{p}}}}}{{{\text{d}}t}} = {r_3}+{r_4}$

(9)

单位时间内HMX热分解反应的反应源项S_HMX为

${S_{\rm {HMX}}} = {\rho _{\rm {HMX}}}{Q_1}{r_1}+{\rho _{\rm {HMX}}}{Q_2}{r_2}+{\rho _{\rm {HMX}}}{Q_3}{r_3}+{\rho _{\rm {HMX}}}{Q_4}{r_4}$

(10)

式中：ρ_HMX为HMX的密度，kg∙m⁻³；Q₁～Q₄为反应1～反应4的反应热，J∙kg⁻¹。

Kel-F为含氟含氯黏结剂，分解时吸收热量，使用单步化学反应动力学模型描述其热分解过程。其热分解反应速率方程为^[15]

${r_5} = {Z_5}\exp \left( { - \frac{{{E_{{\rm a}5}}}}{{RT}}} \right)$

(11)

式中：r₅、Z₅、E_a5分别为Kel-F热分解反应的反应速率、指前因子和活化能。单位时间内Kel-F热分解反应的反应源项S_Kel-F为

${S_{\text {Kel-F}}} = {\rho _{\text {Kel-F}}}{Q_5}{r_5}$

(12)

式中：ρ_Kel-F为Kel-F的密度，Q₅为Kel-F热分解反应的反应热。

AP为氧化剂，其分解反应为放热反应，采用一步反应描述其热分解过程^[16]

${\text{N}}{{\text{H}}_{\text{4}}}{\text{Cl}}{{\text{O}}_{\text{4}}} \to {{\text{O}}_{\text{2}}}{\text+}{{\text{H}}_{\text{2}}}{\text{O+HCl+}}{{\text{N}}_{\text{2}}}$

AP在240 ℃左右由斜方晶型（α-AP）转变为立方晶型（β-AP）^[17]。在慢速烤燃的温升过程中，AP的温度几乎均在240 ℃以下，因此采用α-AP的热分解反应速率方程^[18]，即

${r_6} = \frac{{{\text{d}}\alpha }}{{{\text{d}}t}} = {Z_6}\exp \left( { - \frac{{{E_{{\rm a}6}}}}{{RT}}} \right){\alpha ^n}{\left( {1 - \alpha } \right)^m}$

(13)

式中：r₆、Z₆、E_a6分别为α-AP热分解反应速率、指前因子和活化能；α为AP的热分解转化率，初始阶段为零，完全分解为1；n和m为反应动力学参数，对于α-AP，n=1，m=1.38。根据文献[19]，AP的热分解开始于150 ℃左右，因此在UDF用户子程序中设α在150 ℃以下为零，150 ℃以上为10⁻⁴，以启动对AP分解反应计算。设AP的密度为ρ_AP，则单位时间内AP热分解反应的反应源项S_AP为

${S_{\rm {AP}}} = {\rho _{\rm {AP}}}{Q_6}{r_6}$

(14)

式中：Q₆为AP热分解反应的反应热。

因为Al不参与点火前的热分解反应，所以不考虑铝粉的热分解反应源项，则炸药的总自热反应源项S为

$S = 0.44{S_{{\rm{HMX}}}}+0.11{S_{{\text{Kel-F}}}}+0.20{S_{{\rm{AP}}}}$

(15)

根据叠加原理，炸药的导热系数λ和比热容C分别表示为

$\lambda = 0.44{\lambda _{{\text{HMX}}}}+0.11{\lambda _{{\text{Kel-F}}}}+0.20{\lambda _{\rm {AP}}}$

(16)

$C = 0.44{C_{{\text{HMX}}}}+0.11{C_{{\text{Kel-F}}}}+0.20{C_{\rm {AP}}}$

(17)

表1和表2列出了炸药组分HMX、Kel-F和AP的反应动力学参数，表3列出了材料的物性参数。采用Fluent软件进行计算，设置温升边界条件为1.0和1.5 K/min，通过编写用户自定义函数UDF，将热导率、边界条件以及反应动力学参数加载到主程序中进行计算。

表 1 HMX的反应动力学参数

Table 1. Reaction kinetic parameters of HMX

Reaction No.	Z/s⁻¹	ΔS_f/ (J∙mol⁻¹∙K⁻¹)	E_f/(kJ∙mol⁻¹)	ΔS_r/ (J∙mol⁻¹∙K⁻¹)	E_r/(kJ∙mol⁻¹)	Q/(J∙g⁻¹)
1		123.0	204	89.0	189	−25
2		−41.7	102	−75.2	865	−25
3	3.16×10¹⁶		200			−120
4	1.50×10¹⁴		173			3200

下载: 导出CSV

| 显示表格

表 2 Kel-F和AP的反应动力学参数

Table 2. Reaction kinetic parameters of Kel-F and AP

Component	Reaction No.	Z/s⁻¹	E/(kJ∙mol⁻¹)	Q/(J∙g⁻¹)
Kel-F	5	9.93×10¹⁷	272.0	−5870
AP	6	6.85×10¹²	146.3	297

下载: 导出CSV

| 显示表格

表 3 材料物性参数

Table 3. Parameters of material properties

Component	Density/(kg∙m⁻³)	C/(J·kg⁻¹·K)	λ/(W·m⁻¹·K)
HMX	1850	1004.26	0.5358
Kel-F	2020	1000.43	0.0527
AP	1950	1255	0.2760
Al	2719	871	1.39
Steel	8030	502.48	16.27
Air	1.225	1006.43	0.0242

下载: 导出CSV

| 显示表格

2. 结果与讨论

2.1 烤燃实验

图3为不同温升速率下GOLA-1炸药点火后烤燃弹照片。图3(a)显示：点火后烤燃弹壳体膨胀变形，端盖被冲飞，烤燃弹底部断裂脱落，壳体螺纹完整；壳体表面有明显的烧灼痕迹，壳体内有少量未反应炸药。图3(b)显示了1.5 K/min温升速率下GOLA-1炸药点火前后烤燃弹的对比。结果显示，点火后烤燃弹壳体膨胀变形，端盖被冲飞，壳体底部与侧面有明显的烧灼痕迹，壳体内有少量未反应炸药。由实验现象判定，当温升速率为1.0和1.5 K/min时，GOLA-1炸药的热刺激响应等级均为燃烧，说明炸药在慢速烤燃条件下的热安全性较好。对比两种温升速率下炸药的响应程度，可以发现，温升速率为1.0 K/min时炸药的慢速烤燃响应更剧烈。

图 3 点火后烤燃弹实物

Figure 3. Cook-off bombs after ignition

下载: 全尺寸图片幻灯片

2.2 数值模拟

基于实验进行慢速烤燃数值计算，在1.0和1.5 K/min温升速率下实验与数值模拟得到的炸药中心测温点温度-时间曲线如图4和图5所示。可以看出，实验与数值模拟得到的温升曲线基本重合，点火温度和点火时间也基本吻合，说明数值模型及参数比较合理。图4显示，9900 s时实验测得的温升曲线的温升放缓，而后恢复正常温升速率。分析认为，此阶段的温升放缓是由于HMX在由β相转变为δ相的过程中吸收热量，温度约为451 K。在慢速烤燃实验中，炸药于11946 s时发生点火反应，点火时药柱中心点温度为488 K；计算得到的点火时间为11790 s，与实验结果的相对偏差为1.3%。图5显示：HMX相变造成的温升速率放缓出现在约6800 s，温度约为451 K；实验结果显示，炸药在约8207 s发生点火反应，点火时刻药柱中心温度约为480 K；计算得到的点火时间为8070 s，与实验结果的相对偏差为1.7%。由此可见，温升速率对点火时间的影响较大，随着温升速率的上升，点火时间缩短，然而温升速率对HMX的相变温度以及GOLA-1炸药的点火温度没有影响。

图 4 GOLA-1炸药中心点的温度-时间曲线（1.0 K/min）

Figure 4. Temperature-time curves at the center point of GOLA-1 explosive (1.0 K/min)

下载: 全尺寸图片幻灯片

图 5 GOLA-1炸药中心点的温度-时间曲线（1.5 K/min）

Figure 5. Temperature-time curves at the center point of GOLA-1 explosive (1.5 K/min)

下载: 全尺寸图片幻灯片

1.0和1.5 K/min温升速率下不同时刻炸药的温度分布如图6和图7所示。相比于炸药，壳体的比热容小而热导率大，因此壳体温度高于炸药温度，热量由壳体传递到装药内部。装药的部分热量散失至烤燃弹上部的空气域中，造成装药上部温度较低。图6(d)和图7(d)为点火时烤燃弹的温度分布，可见，点火区域位于装药底部的环形区域。这是由于药柱底部受到侧面壳体与底面的共同加热，使得装药底部温度相对较高。温度的升高使得炸药的自热分解反应速率上升，造成热量进一步累积，如此循环往复，致使点火最终发生在装药底部的环形区域。对比两者发现，在1.5 K/min的温升速率下装药的温度梯度较大，点火位置相对靠外。分析认为，这是由于炸药的导热系数远低于钢壳体，而温升速率相对较大，热量来不及传递到炸药内部造成的。

图 6 1.0 K/min温升速率下不同时刻GOLA-1炸药的温度分布

Figure 6. Temperature distribution of GOLA-1 explosive at different time with the heating rate of 1.0 K/min

下载: 全尺寸图片幻灯片

图 7 1.5 K/min温升速率下不同时刻GOLA-1炸药的温度分布

Figure 7. Temperature distribution of GOLA-1 explosive at different time with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

1.0 K/min温升速率下不同位置的温升曲线和质量分数曲线如图8、图9和图10所示。从图8可以看出，HMX在约9300 s发生相变，相变温度与中心点温度相同，均为451 K。对比中心测温点与点火点的温升曲线可以发现，点火点处HMX相变所造成的温度平台较小，且持续时间较短。点火点处的热量累积速率大于热量传递速率，温度上升较快，相同时刻点火点处的温度高于其他位置，从而使HMX相变时间早且速率快。临近点火时，由于HMX热分解反应释放大量的热量，造成点火点的温度急剧升高，点火温度约为517 K。图9和图10给出了中心测温点和点火点的HMX质量分数变化曲线。在点火时刻，β-HMX已经全部转变为δ-HMX，相变过程基本完成。点火时中心测温点仅有微量的HMX热分解产物，而点火位置则有质量分数为0.2的产物。说明在点火位置，第4步强放热反应已经开始进行，释放大量热量，导致温度急剧上升，主导了点火的发生。

图 8 1.0 K/min温升速率下点火点处的温度-时间曲线

Figure 8. Temperature-time curve of the ignition position with the heating rate of 1.0 K/min

下载: 全尺寸图片幻灯片

图 9 1.0 K/min温升速率下中心测温点处HMX的质量分数变化曲线

Figure 9. Mass fraction curves of HMX at the center point with the heating rate of 1.0 K/min

下载: 全尺寸图片幻灯片

图 10 1.0 K/min温升速率下点火点处HMX的质量分数变化曲线

Figure 10. Mass fraction curves of HMX at the ignition position with the heating rate of 1.0 K/min

下载: 全尺寸图片幻灯片

1.5 K/min温升速率下不同位置的温升曲线和质量分数曲线如图11、图12和图13所示。图11中的点火点温度-时间曲线特征与图8基本相同，点火温度约为518 K。图12和图13显示：相比于温升速率为1.0 K/min的情况，点火时刻中心测温点处β-HMX的质量分数为0.04，δ-HMX的质量分数为0.96，表明部分β-HMX尚未转变为δ-HMX，第4步强放热反应尚未发生；而点火位置的HMX已完全相变，并且HMX分解的第4步强放热反应已经开始进行，有质量分数为0.19的产物生成，主导了点火的发生。

图 11 1.5 K/min温升速率下点火点处的温度-时间曲线

Figure 11. Temperature-time curve of the ignition position with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图 12 1.5 K/min温升速率下中心测温点处HMX的质量分数变化曲线

Figure 12. Mass fraction curves of HMX at center point with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图 13 1.5 K/min温升速率下点火点处HMX的质量分数变化曲线

Figure 13. Mass fraction curves of HMX at ignition position with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图14给出了温升速率为1.5 K/min时中心测温点处AP的转化率-时间曲线。转化率α与分解速率dα/dt均随着温度的升高而增大。在点火时刻，AP的转化率为1.06×10⁻⁴，表明此时中心点的AP基本不发生分解。7000 s左右，AP的分解速率放缓，而后快速上升，这是HMX在晶型转变过程中吸收热量而导致温升放缓造成的。图15为点火点处AP的转化率-时间曲线。可以看出：转化率在临近点火时刻迅速上升；点火时刻AP的转化率为2.26×10⁻⁴，相较于中心点处的转化率高约1倍。不论是中心测温点还是点火点，发生分解的AP的量均很少，表明在慢速烤燃过程中AP基本不参与反应，点火反应由HMX主导。这与Kim等^[6]得到的含硝胺类及AP浇注炸药慢速烤燃实验结论一致。

图 14 1.5 K/min温升速率下中心测温点的AP转化率曲线

Figure 14. Conversion rate curve of AP at the center point with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

图 15 1.5 K/min温升速率下点火点的AP转化率曲线

Figure 15. Conversion rate curve of AP at ignition position with the heating rate of 1.5 K/min

下载: 全尺寸图片幻灯片

2.3 温升速率对点火位置的影响

将实验与数值模拟相结合，分析温升速率对炸药热刺激响应特性的影响。在1.0和1.5 K/min两种温升速率的基础上，设定不同的温升速率，对GOLA-1炸药进行计算。当温升速率为0.5和0.4 K/min时，炸药点火位置的变化如图16所示。0.5 K/min时，点火位置仍为装药底部的环形区域，相比1.0和1.5 K/min的情况，环形区域缩小，并向装药的上部及中心移动。这是由于当温升速率下降时，炸药整体受热比较均匀，装药内部的温度梯度减小，相比于温升速率为1.0 K/min的情况，相同加热温度下炸药内部高温区域的体积较大，炸药的自热反应速率上升，炸药内部积累热量，导致点火发生在装药底部的环形区域。当温升速率降至0.4 K/min时，点火位置的形状发生改变，由环形变为扁的椭球形，并且位置移动至药柱下部中心。

图 16 不同温升速率下点火位置的对比

Figure 16. Comparison of ignition positions with different heating rates

下载: 全尺寸图片幻灯片

为进一步分析温升速率对GOLA-1点火位置的影响，以0.05 K/min为间隔，进行0.1～0.4 K/min 7个不同温升速率的烤燃数值模拟。设温升速率为β，以壳体外侧底面的中心位置为坐标零点，沿烤燃弹中心轴线从壳体底面指向端盖为正方向，点火位置x为温度最高区域的中心，单位为m。温升速率对点火位置的影响如表4所示。随着温升速率的下降，点火位置向药柱上部移动。根据表4绘制点火位置随温升速率的散点图，如图17所示。观察图17发现，温升速率β与点火位置x之间存在线性关系，拟合关系式为 $x = - 0.052\,79\,\beta +0.041\,71$ ，相关系数 R²=0.966。

表 4 温升速率与点火位置的关系

Table 4. Heating rate vs. ignition position

No.	β/(K·min⁻¹)	x/m
1	0.40	0.0212
2	0.35	0.0220
3	0.30	0.0273
4	0.25	0.0287
5	0.20	0.0298
6	0.15	0.0334
7	0.10	0.0374

下载: 导出CSV

| 显示表格

图 17 温升速率对点火位置的影响

Figure 17. Effect of heating rate on ignition position

下载: 全尺寸图片幻灯片

3. 结　论

(1) 对新型HMX基含AP浇注炸药GOLA-1进行慢速烤燃实验，结果显示，当温升速率为1.0和1.5 K/min时炸药的烤燃响应等级均为燃烧，并且1.0 K/min下壳体的破坏更严重，响应更剧烈。从热刺激响应角度来看，温升速率越低，GOLA-1炸药的响应越剧烈。

(2) 开展了GOLA-1炸药烤燃数值模拟，得到的温升曲线与实验结果吻合较好。当温升速率为1.0 K/min时，GOLA-1炸药的点火时间与实验结果的相对偏差为1.3%，1.5 K/min时为1.7%，证明建立的烤燃模型及参数较合理。数值模拟结果显示，点火反应主要由HMX的第4步强放热反应主导，AP在慢速烤燃过程中基本不参与点火反应。

(3) 对不同温升速率下GOLA-1炸药的烤燃过程进行了数值模拟，结果表明：当温升速率为0.5 K/min时，点火仍发生在装药底部的环形区域，并向装药中心及上部移动；当温升速率下降至0.4 K/min时，点火区域的形状发生变化，由环形变为椭球形；随着温升速率继续下降，点火位置沿着药柱中心轴线向上部移动；当温升速率处在0.1～0.4 K/min之间时，温升速率与点火位置之间存在线性关系。

图 1 准静态拉伸测试：(a) 拉伸测试前， (b) 拉伸测试后

Figure 1. Quasi-static tensile test: (a) before the tensile test, (b) after the tensile test

下载: 全尺寸图片幻灯片

图 2 夹具示意图

Figure 2. Schematic diagram of the fixture

下载: 全尺寸图片幻灯片

图 3 UHMWPE纤维束动态拉伸试件

Figure 3. Specimens of UHMWPE fiber bundles for dynamic tensile

下载: 全尺寸图片幻灯片

图 4 动态拉伸波形

Figure 4. Waveform obtained by dynamic tensile

下载: 全尺寸图片幻灯片

图 5 准静态拉伸应力-应变曲线

Figure 5. Stress-strain curves obtained by quasi-static tensile

下载: 全尺寸图片幻灯片

图 6 动态拉伸应力-应变曲线

Figure 6. Stress-strain curves obtained by dynamic tensile

下载: 全尺寸图片幻灯片

图 7 模量和强度随应变率的变化曲线

Figure 7. Modulus and strength versus strain rate

下载: 全尺寸图片幻灯片

图 8 纤维织物纱线编织单胞模型

Figure 8. Unit cell model of fiber fabric yarn weaving

下载: 全尺寸图片幻灯片

图 9 单层纤维织物几何模型

Figure 9. Geometric model of single-layer fiber fabric

下载: 全尺寸图片幻灯片

图 10 纤维织物超高速碰撞数值模型

Figure 10. Numerical model of fiber fabric subjected to hypervelocity impact

下载: 全尺寸图片幻灯片

图 11 模拟结果与实验结果的对比：(a) 弹丸速度，(b) 动能吸收

Figure 11. Comparison between the numerical results and experiment: (a) projectile velocity, (b) energy absorption

下载: 全尺寸图片幻灯片

图 12 不同时刻的碎片云图

Figure 12. Debris cloud map at different moments

下载: 全尺寸图片幻灯片

图 13 纤维织物应力云图

Figure 13. Stress nephogram of fiber fabric

下载: 全尺寸图片幻灯片

图 14 弹丸动能吸收率随应变率敏感系数变化曲线

Figure 14. Energy absorption ratio versus strain-rate sensitivity coefficient

下载: 全尺寸图片幻灯片

参考文献(32)

[1]	常利军, 黄星源, 袁圣林, 等. 压缩载荷下UHMWPE纤维复合材料层合板的力学性能与失效分析 [J]. 高压物理学报, 2023, 37(1): 014102. doi: 10.11858/gywlxb.20220633 CHANG L J, HUANG X Y, YUAN S L, et al. Mechanical properties and failure analysis of UHMWPE fiber composite laminates under compressive load [J]. Chinese Journal of High Pressure Physics, 2023, 37(1): 014102. doi: 10.11858/gywlxb.20220633
[2]	董澎, 王柯, 李军方, 等. 超高分子量聚乙烯烧结制品的链缠结调控及其对性能影响 [J]. 高分子学报, 2020, 51(1): 117–124. doi: 10.11777/j.issn1000-3304.2020.19159 DONG P, WANG K, LI J F, et al. Chain entanglement regulation of sintered ultrahigh molecular weight polyethylene and its effect on properties [J]. Acta Polymerica Sinica, 2020, 51(1): 117–124. doi: 10.11777/j.issn1000-3304.2020.19159
[3]	付杰, 李伟萍, 黄献聪, 等. 新型超高分子量聚乙烯膜材料防弹性能及机理 [J]. 兵工学报, 2021, 42(11): 2453–2464. doi: 10.3969/j.issn.1000-1093.2021.11.019 FU J, LI W P, HUANG X C, et al. Bullet-proof performance and mechanism of new ultra-high molecular weight polyethylene film [J]. Acta Armamentarii, 2021, 42(11): 2453–2464. doi: 10.3969/j.issn.1000-1093.2021.11.019
[4]	莫根林, 刘静, 金永喜, 等. 超高分子量聚乙烯纤维防护机理研究综述 [J]. 兵器装备工程学报, 2021, 42(10): 23–28. doi: 10.11809/bqzbgcxb2021.10.004 MO G L, LIU J, JIN Y X, et al. Review on protective mechanism of UHMWPE fiber [J]. Journal of Ordnance Equipment Engineering, 2021, 42(10): 23–28. doi: 10.11809/bqzbgcxb2021.10.004
[5]	袁子舜, 陆振乾, 许玥, 等. 超高分子量聚乙烯纤维平纹织物-单向布混合堆叠板的防弹机制[J]. 复合材料学报, 2022, 39(6): 2707−2715. YUAN Z S, LU Z Q, XU Y, et al. Ballistic mechanism of the hybrid panels with UHMWPE woven fabrics and UD laminates [J].Acta Materiae Compositae Sinica, 2022, 39(6): 2707−2715.
[6]	苗常青, 徐铧东, 靳广焓, 等. 纤维编织材料超高速撞击特性实验研究 [J]. 高压物理学报, 2019, 33(2): 024203. doi: 10.11858/gywlxb.20180654 MIAO C Q, XU H D, JIN G H, et al. Experimental study of hypervelocity impact characteristics for fiber fabric materials [J]. Chinese Journal of High Pressure Physics, 2019, 33(2): 024203. doi: 10.11858/gywlxb.20180654
[7]	苗常青, 徐铧东, 杜明俊, 等. 芳纶/环氧纤维复合材料超高速撞击特性研究 [J]. 实验力学, 2019, 34(4): 609–615. doi: 10.7520/1001-4888-17-231 MIAO C Q, XU H D, DU M J, et al. On the hypervelocity impact characteristics of aramid/epoxy fiber composite [J]. Journal of Experimental Mechanics, 2019, 34(4): 609–615. doi: 10.7520/1001-4888-17-231
[8]	张鹏, 王志军, 马武伟, 等. 高速弹体侵彻钢/陶瓷/超高分子量聚乙烯纤维/钢实验 [J]. 兵器材料科学与工程, 2016, 39(5): 104–109. doi: 10.14024/j.cnki.1004-244x.20160826.007 ZHANG P, WANG Z J, MA W W, et al. Experiment on high velocity projectile penetrating composite target of steel/ceramic/ultra-high molecular weight polyethylene fiber/steel [J]. Ordnance Material Science and Engineering, 2016, 39(5): 104–109. doi: 10.14024/j.cnki.1004-244x.20160826.007
[9]	张宝玺, 哈跃, 邓云飞, 等. 超高速撞击Kevlar纤维布填充防护结构研究 [J]. 高压物理学报, 2013, 27(1): 105–112. doi: 10.11858/gywlxb.2013.01.015 ZHANG B X, HA Y, DENG Y F, et al. Optimal structural design of stuffed shields with Kevlar fiber clothes against hypervelocity impact [J]. Chinese Journal of High Pressure Physics, 2013, 27(1): 105–112. doi: 10.11858/gywlxb.2013.01.015
[10]	石景富. 聚合物超高速碰撞特性及应变率效应分析 [D]. 哈尔滨: 哈尔滨工业大学, 2021. SHI J F. Analysis of polymer hypervelocity impact characteristics and strain rate effect [D]. Harbin: Harbin Institute of Technology, 2021.
[11]	赵荣国, 陈朝中, 罗文波, 等. 聚合物材料SHPB实验关键问题 [J]. 固体力学学报, 2011, 32(Suppl 1): 134–144. ZHAO R G, CHEN C Z, LUO W B, et al. Key problems of SHPB experiments used for polymeric materials [J]. Acta Mechanica Solida Sinica, 2011, 32(Suppl 1): 134–144.
[12]	林玉亮, 卢芳云, 卢力. 高应变率下硅橡胶的本构行为研究 [J]. 高压物理学报, 2007, 21(3): 289–294. doi: 10.3969/j.issn.1000-5773.2007.03.012 LIN Y L, LU F Y, LU L. Constitutive behaviors of silicone rubber at high strain rates [J]. Chinese Journal of High Pressure Physics, 2007, 21(3): 289–294. doi: 10.3969/j.issn.1000-5773.2007.03.012
[13]	王庭辉, 宋顺成, 王明超, 等. 高强度纤维束的动态拉伸性能 [J]. 西南交通大学学报, 2008(5): 638–642. doi: 10.3969/j.issn.0258-2724.2008.05.016 WANG T H, SONG S C, WANG M C, et al. Dynamic tensile properties of high strength fiber bundles [J]. Journal of Southwest Jiaotong University, 2008(5): 638–642. doi: 10.3969/j.issn.0258-2724.2008.05.016
[14]	YANG B, XIONG T, XIONG J. Statistical tensile strength for high strain rate of aramid and UHMWPE fibers [J]. Journal of Materials Engineering, 2006(5): 46–50.
[15]	CHEN L, ZHENG K, FANG Q. Effect of strain rate on the dynamic tensile behaviour of UHMWPE fiber laminates [J]. Polymer Testing, 2017, 63: 54–64. doi: 10.1016/j.polymertesting.2017.07.031
[16]	栗建桥, 宋卫东, 宁建国. 超高速撞击产生的等离子体特性研究 [J]. 高压物理学报, 2013, 27(4): 542–548. doi: 10.11858/gywlxb.2013.04.012 LI J Q, SONG W D, NING J G. A study on characteristics of plasma generated by hypervelocity impact [J]. Chinese Journal of High Pressure Physics, 2013, 27(4): 542–548. doi: 10.11858/gywlxb.2013.04.012
[17]	林健宇, 罗斌强, 徐名扬, 等. 铝弹丸超高速撞击防护结构的研究进展 [J]. 高压物理学报, 2019, 33(3): 030112. doi: 10.11858/gywlxb.20190774 LIN J Y, LUO B Q, XU M Y, et al. Progress of aluminum projectile impacting on plate with hypervelocity [J]. Chinese Journal of High Pressure Physics, 2019, 33(3): 030112. doi: 10.11858/gywlxb.20190774
[18]	DHOTE D, VERMA P N. Investigation of hole formation by steel sphere impacting on thin plate at hypervelocity [J]. Thin-Walled Structures, 2018, 126: 38–47.
[19]	张祎, 王玉林, 石景富, 等. 纤维织物超高速碰撞热-力学模型与分析 [J]. 宇航学报, 2021, 42(11): 1475–1482. doi: 10.3873/j.issn.1000-1328.2021.11.014 ZHANG Y, WANG Y L, SHI J F, et al. Modeling and analysis of thermo-mechanical behavior for fiber fabric under hypervelocity impact [J]. Journal of Astronautics, 2021, 42(11): 1475–1482. doi: 10.3873/j.issn.1000-1328.2021.11.014
[20]	徐铧东, 于东, 王玉林, 等. 预张力纤维织物超高速碰撞热-力学特性分析 [J]. 爆炸与冲击, 2022, 42(5): 053301. doi: 10.11883/bzycj-2021-0307 XU H D, YU D, WANG Y L, et al. Thermo-mechanical characteristics of pre-tensioned fiber fabrics subjected to hypervelocity impact [J]. Explosion and Shock Waves, 2022, 42(5): 053301. doi: 10.11883/bzycj-2021-0307
[21]	SHINTATE K, SEKINE H. Numerical simulation of hypervelocity impacts of a projectile on laminated composite plate targets by means of improved SPH method [J]. Composites Part A: Applied Science & Manufacturing, 2004, 35(6): 683–692.
[22]	ZHAO S, SONG Z, ESPINOSA H D. Modelling and analyses of fiber fabric and fabric-reinforced polymers under hypervelocity impact using smooth particle hydrodynamics [J]. International Journal of Impact Engineering, 2020, 144: 103586. doi: 10.1016/j.ijimpeng.2020.103586
[23]	徐铧东, 王玉林, 刘蕾, 等. 纤维织物FEM-SPH耦合单胞模型及超高速碰撞特性 [J]. 复合材料学报, 2021, 38(9): 3123–3132. doi: 10.13801/j.cnki.fhclxb.20201231.001 XU H D, WANG Y L, LIU L, et al. A fiber fabric unit-cell model based on FEM-SPH coupling algorithm and application on analyses of hypervelovcity impact [J]. Acta Materiae Compositae Sinica, 2021, 38(9): 3123–3132. doi: 10.13801/j.cnki.fhclxb.20201231.001
[24]	KATZ S, GROSSMAN E, GOUZMAN I, et al. Response of composite materials to hypervelocity impact [J]. International Journal of Impact Engineering, 2008, 35(12): 1606–1611. doi: 10.1016/j.ijimpeng.2008.07.032
[25]	HEIMBS S, WAGNER T, VIANA J, et al. Comparison of impact behaviour of glass, carbon and dyneema composites [J]. Journal of Mechanical Engineering Science, 2019, 233(3): 951–966.
[26]	ROGERS J A, MOTE A, MEAD P T, et al. Hypervelocity impact response of monolithic UHMWPE and HDPE plates [J]. International Journal of Impact Engineering, 2022, 161: 104081.
[27]	WANG H X, WEERASINGHE D, MOHOTTI D, et al. On the impact response of UHMWPE woven fabrics: experiments and simulations [J]. International Journal of Mechanical Sciences, 2021, 204: 106574.
[28]	CHA J H, KIM Y H, SARATH K, et al. Ultra-high-molecular-weight polyethylene as a hypervelocity impact shielding material for space structures [J]. Acta Astronautica, 2020, 168: 182–190.
[29]	杨鹏飞, 汪洋, 夏源明. 基于Hopkinson杆的材料高应变率拉伸实验技术 [J]. 实验力学, 2011, 26(6): 674–679. YANG P F, WANG Y, XIA Y M. Experimental technique of high strain-rate tension based on Hopkinson bar [J]. Journal of Experimental Mechanics, 2011, 26(6): 674–679.
[30]	朱德举, 张晓彤, 张怀安. 动态拉伸载荷下应变率和温度对Kevlar 49芳纶纤维布增强环氧树脂复合材料力学性能的影响 [J]. 复合材料学报, 2016, 33(3): 459–468. ZHU D J, ZHANG X T, ZHANG H A. Effects of strain rate and temperature on mechanical properties of Kevlar 49 aramid fabric reinforced epoxy polymers under dynamic tensile loading [J]. Acta Materiae Compositae Sinica, 2016, 33(3): 459–468.
[31]	JOHNSON G R, COOK W H. A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures [J]. Engineering Fracture Mechanics, 1983, 21: 541–548.
[32]	彭建祥. Johnson-Cook本构模型和Steinberg本构模型的比较研究 [D]. 绵阳: 中国工程物理研究院, 2006. PENG J X. Comparative study of Johnson-Cook constitutive model and Steinberg constitutive model [D]. Mianyang: China Academy of Engineering Physics, 2006.

施引文献

资源附件(0)

访问统计

点击查看大图

图(14)

计量

文章访问数: 281
HTML全文浏览量: 126
PDF下载量: 41

1. 实验与数值模拟方法
1.1 小型烤燃弹实验
1.2 烤燃数值模拟
2. 结果与讨论
2.1 烤燃实验
2.2 数值模拟
2.3 温升速率对点火位置的影响
3. 结　论

UHMWPE的应变率效应及其对超高速碰撞特性的影响

doi: 10.11858/gywlxb.20220666

作者简介: 石景富（1994－），男，博士研究生，主要从事聚合物材料力学性能研究. E-mail：shijf@stu.hit.edu.cn

通讯作者: 苗常青（1972－），男，博士，教授，主要从事柔性复合材料研究. E-mail：miaocq@hit.edu.cn

计量

出版历程

Strain Rate Effect of UHMWPE and Its Influence on Hypervelocity Impact Performance

1. 实验与数值模拟方法

1.1 小型烤燃弹实验

1.2 烤燃数值模拟

2. 结果与讨论

2.1 烤燃实验

2.2 数值模拟

2.3 温升速率对点火位置的影响

3. 结 论

计量

出版历程

目录

1. 实验与数值模拟方法

1.1 小型烤燃弹实验

1.2 烤燃数值模拟

2. 结果与讨论

2.1 烤燃实验

2.2 数值模拟

2.3 温升速率对点火位置的影响

3. 结 论

作者简介:
石景富（1994－），男，博士研究生，主要从事聚合物材料力学性能研究. E-mail：shijf@stu.hit.edu.cn

通讯作者:
苗常青（1972－），男，博士，教授，主要从事柔性复合材料研究. E-mail：miaocq@hit.edu.cn

3. 结　论

3. 结　论