可变形定向破片战斗部在不同展开模式下的性能数值模拟

杜伊杨; 黄炫宁; 杨鸿铭; 周晴; 李伟兵

doi:10.11858/gywlxb.20210887

可变形定向破片战斗部在不同展开模式下的性能数值模拟

doi: 10.11858/gywlxb.20210887

南京理工大学智能弹药技术国防重点学科实验室, 江苏南京　210094

基金项目: 国家级大学生创新创业训练计划（202010288007Z）

详细信息

作者简介:
杜伊杨（1999-），男，本科，主要从事战斗部设计研究. E-mail：duyiyang@njust.edu.cn

通讯作者:
李伟兵（1982-），男，博士，教授，博士生导师，主要从事高效毁伤技术研究.E-mail：njustlwb@163.com

中图分类号: O383; TJ410.33
计量
- 文章访问数: 1764
- HTML全文浏览量: 483
- PDF下载量: 43
出版历程
- 收稿日期: 2021-09-30
- 修回日期: 2021-10-15
- 录用日期: 2021-10-20
- 刊出日期: 2022-05-30

Numerical Simulation of Deformable Directional Fragment Warhead Performance in Different Expansion Modes

Ministerial Key Laboratory of ZNDY, Nanjing University of Science and Technology, Nanjing 210094, Jiangsu, China

摘要

摘要: 为使战斗部具有多种定向毁伤模式并实现一定程度上的可控毁伤，提出了一种扇形装药的可变形定向破片战斗部，该战斗部可实现轴向展开和侧向展开2种模式。采用AUTODYN软件进行破片场的数值模拟。首先，基于战斗部单元体分析获得了距离轴心25 mm处的最佳起爆点位置；其次，对整个战斗部进行分析，在轴向展开模式下分析了轴向展开角度对破片飞散速度、破片数目和破片空间分布的影响，发现轴向展开角在60°～75°范围内毁伤效果较佳；最后，在侧向展开模式下分析了整个战斗部的破片速度和破片空间分布情况，结果表明破片具有明显的定向飞散特性。
- 定向战斗部 /
- 轴向展开 /
- 侧向展开 /
- 破片空间分布
Abstract: In order to enable warhead to have multiple directional damage modes and to achieve a certain degree of controllable damage, a deformable directional fragment warhead with fan-shaped charge is proposed. It can realize an axial-expanding mode and an lateral-expanding mode. The AUTODYN software is used to carry out the numerical simulation of the fragment field. Firstly, the position of the initiation point at 25 mm away from the central axis is determined based on the analysis of the unit constituting the warhead. Secondly, for the axial-expanding mode, the influence of the axial-expansion angle on the fragment flying speed, the quantity and the spatial distribution of fragments are obtained. The results show that the damage effect of the warhead is better when the axial-expansion angle is in the range of 60°–75°. Thirdly, for the lateral-expansion mode, the velocity and the spatial distribution of fregments are analyzed. It showes that the fragments under the lateral-expansion mode have directional flying quality.
- directional warhead /
- axial-expansion /
- lateral-expansion /
- fragment spatial distribution

HTML全文

传感器技术、通信技术、计算机技术构成了现代信息技术的3大支柱。压电式传感器以频响高、体积小、质量轻、频带宽等特点^[1-3]，成为高技术领域中制备小型化及多功能化器件的重要元件之一。薄膜型压电传感器的敏感单元为铁电聚合物薄膜，在电子、超声、水声、红外、导航、生物等多个领域应用广泛。早在19世纪40年代人们便开始了高分子聚合物压电性能研究^[4-6]，但直到70年代才将其制成压电薄膜传感器。压电聚合物通常为非导电性高分子材料，从原理上看没有可移动的电子电荷，但在某些特定的条件下（经延展拉伸、极化等），带负电荷的引力中心可以被改变，从而成为具有压电性的高分子压电薄膜。常见的高分子压电薄膜有聚氟乙烯（PVF）、聚偏二氟乙烯（PVDF）、聚氯乙烯（PVC）、尼龙11等。

动高压聚合物薄膜测量技术的核心内容之一是敏感单元压电性能研究^[7-9]。PVDF传感器的敏感单元是PVDF压电薄膜。PVDF薄膜属于半结晶的高分子聚合物，晶区至少存在4种晶型结构（α相、β相、γ相和δ相），其压电性能直接与极性β相相关。在常温常压下，通过PVDF溶液挥发或熔融结晶得到的PVDF初始膜主要以非极性α相为主；对PVDF初始膜进行单轴拉伸或高电压极化处理，可产生具有更多压电β相的PVDF压电薄膜。

本研究首先探讨PVDF传感器的设计和制作工艺，包括初始膜制备、单轴拉伸取向诱导压电相、高压热极化试验、封装技术等；然后采取溶液等温结晶法制备初始PVDF薄膜，对初始膜进行单轴拉伸试验，分析不同拉伸倍率、拉伸速率和拉伸温度对PVDF薄膜中结晶相的影响；接着在高电场交流极化实验中，分析不同电场频率和幅值对PVDF薄膜位移电流的影响，获得薄膜本征位移电流曲线；最后对传感器的形状、尺寸、厚度进行设计，封装完成PVDF压电薄膜传感器。希望本研究结果有助于提高PVDF薄膜材料在动态冲击过程中的检测能力和可靠性，推进响应快、测量范围宽、灵敏度高的聚合物薄膜测量技术在动高压冲击实验中的应用。

1. PVDF敏感单元设计

PVDF压力计上、下电极的重合部分为其敏感单元。当外力使PVDF压力计的敏感单元发生变形时，薄膜敏感单元表面就会产生与应力大小成正比的电荷量。本研究设计的敏感单元尺寸为3 mm × 3 mm，薄膜厚度选用 $20 \;{\text{μ}}{\rm{m}}$ 或 $30 \;{\text{μ}}{\rm{m}}$ 。PVDF传感器总厚度小于0.1 mm，有利于提高响应时间。敏感单元电极结构如图1所示。

图 1 敏感单元形状设计

Figure 1. Sensitive element design

下载: 全尺寸图片幻灯片

1.1 电极制作和极化

敏感单元使用双层电极结构，下电极Pt的厚度为50 nm，上电极Au的厚度为250 nm。敏感单元的尺寸直接影响PVDF薄膜的铁电性能测试，因此要精确控制敏感单元的尺寸。为了获得特定形状的内电极金属结构，使用掩膜法在PVDF薄膜衬底上的特定位置沉积金属电极结构。为了提高敏感单元面积精度，采用双面掩膜板结构，上、下掩膜板的对位精度控制在 ± 0.05 mm以内。PVDF薄膜的居里温度约为170 ℃，对热敏感，在镀电极过程中温度的上升会引起电极与PVDF薄膜间的热应力增大，影响电极和PVDF薄膜的黏结度，且对薄膜性能造成影响，因此在沉积金属电极时，需对沉积温度进行严格控制。上电极Au的导电性好，但与PVDF薄膜的黏合度较差，而金属Pt与PVDF薄膜的黏合度较好，但是活性较低，所以采用Au/Pt双层金属电极结构，能够发挥两种金属的各自优点，使电极结构的附着力强、导电性好。

在溅射金属薄膜电极过程中，有两个溅射参数需要重点关注，即溅射功率和工作压强。当溅射功率较大时，电离出的Ar⁺的能量较大，轰击靶材溅射出更多的粒子，粒子的迁移速率变大，易形成连续电极结构；但溅射功率过高时，溅射出的高能粒子将会产生反溅射作用，对薄膜衬底造成损害。当工作压强增大时，溅射粒子向薄膜运动的过程中受到室内气体散射的几率增大，薄膜上沉积的粒子能量减少，且平均粒径减小。过高的工作压强会降低沉积速率，不利于形成连续电极；而过低的工作压强会使沉积的电极表面粗糙度增大。本实验采用多次间歇溅射，溅射温度为室温。所得的PVDF薄膜衬底沉积金属薄膜内电极结构如图2（a）所示，可见，金属电极与薄膜无褶皱，电极边界清晰。

图 2 PVDF薄膜电极结构和热极化装置

Figure 2. Schematic diagram of PVDF film electrode and thermal polarization

下载: 全尺寸图片幻灯片

对于PVDF薄膜，一般采用拉伸工艺获得含量较多的压电β相，但此时薄膜中的偶极子排列无序，宏观上不显电性，需要经过人工极化处理。通过施加高电压可使薄膜中的偶极子沿外电场方向有序排列，当去除电场后，一部分偶极取向瞬间消失，另一部分偶极取向保留下来，称为剩余极化强度，它直接反映PVDF薄膜的压电性。PVDF是偶极电荷型压电铁电聚合物，极化以偶极子的转向优先，常采用热极化方法和电晕放电方法极化，此时偶极子取向由沉积在电介质中的空间电荷场诱导产生。电晕极化是在薄膜表面注入新的导电载流子，使薄膜的电导增强，击穿场强降低。热极化方法具有设备简单、操作方便、极化彻底等优点。因为空气的击穿场强较低，所以本研究采用热极化方法对PVDF薄膜进行高压热极化处理，如图2（b）所示。在样品盒内充满绝缘油，本研究选用硅油。硅油的介电常数高，将硅油作为绝缘介质时，可避免高电压的边缘飞弧，可施加极高的极化电压。热极化过程中伴随着分子链、离子、空间电荷的运动，导致薄膜表面弯曲变形，影响敏感面积的稳定极化，因此极化过程需要对样品施加一定的压强。

图3（a）显示了对 $20 \;{\text{μ}}{\rm{m}}$ 厚PVDF薄膜施加不同幅值极化电场时所得的位移电流曲线。可以看出，随着电场强度的增加，位移电流的电阻、电容效应越来越弱，说明350 MV/m的场强为 $20 \;{\text{μ}}{\rm{m}}$ 厚PVDF薄膜的最佳热极化条件，此时的位移电流和电滞回线如图3（b）所示。

图 3

$20 \;{\text{μ}}{\rm{m}}$ 厚PVDF薄膜的热极化

Figure 3. Effect of polarized electric field on

$20 \;{\text{μ}}{\rm{m}}$ thick PVDF film

下载: 全尺寸图片幻灯片

1.2 敏感单元封装

使用磁控溅射完成PVDF传感器内电极的制作后，内电极为Au/Pt双层金属电极结构，总厚度为300 nm。使用紫铜箔作为外电极，并用聚酰亚胺（PI）薄膜进行封装。因为PVDF本身为高阻元器件，外部保护膜的阻抗应大于传感器本身阻抗2个数量级以上，PI薄膜的电阻率为1 014～1 015 ${\text{Ω}}$ ·m，PVDF薄膜体积电阻率为1 011 ${\text{Ω}}$ ·m，所以PI薄膜厚度应不小于PVDF薄膜厚度的1/10即可满足使用要求。本研究使用的PVDF薄膜的厚度分别为 $20 \;{\text{μ}}{\rm{m}}$ 和 $30 \;{\text{μ}}{\rm{m}}$ ，PI封装薄膜厚度为 $25 \;{\text{μ}}{\rm{m}}$ 。图4为PVDF传感器所示3层结构封装示意图，上、下内电极的重叠部分为敏感单元，其尺寸为3 mm × 3 mm，内电极与外电极铜箔使用冷压技术连接，然后将低损耗同轴电缆与外电极铜箔焊接。

图 4 PVDF传感器结构

Figure 4. PVDF sensor structure

下载: 全尺寸图片幻灯片

2. 动态冲击设计和标定

2.1 测试系统组件选择

当外界压力使压电材料发生变形时，压电材料的表面会产生与应力大小成正比的电荷量，电介质内部的电偶极子排序改变上下表面的电势，与表面接触的两个电极上的电荷重新调整，以平衡PVDF表面电势，因此在电路中存在电荷流动。将PVDF压电薄膜等效为电容，测量所产生的电荷，可采用两种模式测量：电流模式和电荷模式。电荷法适合信号变化要求不高的数据采集模式，频响（上升时间）受电荷积分器响应特性的影响；电流法适合于高频响数据采集，尤其是压力高于20 GPa的数据采集，缺点是误差与测试电缆长度的关联性高。考虑到本研究在低压范围标定，电荷积分器满足频率响应，因此选取电荷模式。

电荷模式是将传感器并联一个匹配电容，将其产生的电荷Q(t)导入该电容，然后外接显示设备，输出与外界压力σ(t)直接对应的电压信号V(t)，该方法也称直接测量法，如图5(a)所示，其中R_a传感器内阻。电荷测试部分采用无源RC积分电路，如图5(b)所示。在无源RC积分电路中，若时间常数τ = RC足够大，则外加电压时，电容上的电压只能慢慢上升。所以在测量时间远小于时间常数τ的时间范围内，电路达到稳态的时间长，电容C两端的电压很小，输入电压主要降落在电阻R上，此时充电电流i_c和输出电压u₀(t)为

图 5 压电薄膜传感器电荷模式等效电路

Figure 5. The equivalent circuit of electric charge mode of PVDF film sensor

下载: 全尺寸图片幻灯片

${i_{\rm{c}}} = {u_{\rm{i}}}\left( t \right)/R$

(1)

${u_0}\left( t \right) = \frac{1}{C}\int_0^t {{i_{\rm{c}}}} {\rm{d}}t = \frac{1}{C}\int_0^t {\frac{{{u_{\rm{i}}}\left( t \right)}}{R}{\rm{d}}t = \frac{1}{{RC}}\int_0^t {{u_{\rm{i}}}\left( t \right)} {\kern 1pt} {\rm{d}}t}$

(2)

式中：u_i(t)为输入电压。可以看出，输出电压u₀(t)近似与输入电压u_i(t)的时间积分值成比例。 $\int_0^t {{{[{u_{\rm{i}}}\left( t \right)} / {R]}}} {\rm{d}}t$ 为流过电阻R的电流积分，即电路中所产生的电荷Q(t)，所以该种测量方式为电流模式，此时电路中的电荷为

${u_0}\left( t \right) \cdot C{\rm{ = }}Q\left( t \right)$

(3)

测试元件中电容C的选择依据：RC积分电路中时间常数τ远远大于输入脉冲宽度（一般超过10倍）。计算时间常数τ ，若电容C的单位为F，R的单位为 ${\text{Ω}}$ ，则时间常数τ的单位为s。当C = $0.1 \;{\text{μ}}{\rm{F}}$ 时

$\tau = RC = 50 \times 0.1 \times {10^{ - 6}}\;{\rm{s}} = 5{\kern 1pt} 000{\kern 1pt} {\kern 1pt} {\kern 1pt} {\rm{ns}}$

(4)

对于所测系统而言，频响（上升时间）为50～150 ns，所以选用C = $0.1 \;{\text{μ}}{\rm{F}}$ 的电容，能够保证 $\tau \gg t$ （10倍以上）。

本研究选用50 ${\text{Ω}}$ 特征阻抗的同轴电缆，阻抗匹配的测试系统电路如图6所示，其中C₁为电缆电容，R₁为传感器泄漏电阻和电缆电阻，Z₀为同轴电缆特征阻抗（50 ${\text{Ω}}$ ），R_c是与同轴电缆阻抗匹配的电阻阻值（50 ${\text{Ω}}$ ），C为0.1 ${\text{μ}}{\rm{F}}$ 。

图 6 测试系统测量电路

Figure 6. Measurement circuit of PVDF film sensor test system

下载: 全尺寸图片幻灯片

2.2 PVDF传感器动态冲击压缩实验标定方法

采用对称碰撞方法（飞片和靶为同种材料）标定PVDF传感器的动态压缩曲线。所选材料为PMMA、LY12铝和45钢，将PVDF传感器粘在样品和飞片的撞击面上，如图7所示。飞片击靶速度u₀采用激光遮断法测量，利用飞片遮断测速环中不同路光纤的激光信号得到相应的时间差，测速环由3路光纤组成，每路光纤间距8.5 mm。

图 7 对称碰撞实验示意图

Figure 7. Experimental scheme of symmetrical impact

下载: 全尺寸图片幻灯片

按照应力波理论中的阻抗匹配法，对称碰撞过程中靶中冲击波阵面后粒子速度u_p是飞片击靶速度u₀的1/2，结合样品材料的冲击Hugoniot关系式，通过测量撞击速度u₀可确定标定压力值

$\sigma = {\rho _0}D{u_{\rm{p}}} = {\rho _0}\left( {{c_0} + \lambda {u_{\rm{p}}}} \right){u_{\rm{p}}} = \frac{1}{2}{\rho _0}\left( {{c_0} + \frac{1}{2}\lambda {u_0}} \right){u_0}$

(3)

式中：D为冲击波速度，km/s；ρ₀为样品初始密度，g/cm³；c₀为线性Hugoniot关系式中常系数，km/s； $\lambda$ 为线性Hugoniot关系式中常系数；u₀为飞片撞击速度，km/s；u_p为样品中的粒子速度，km/s。式（5）给出了撞击压力σ与飞片击靶速度u₀的对应关系。实验中所使用撞击材料的Hugoniot参数列于表1。

表 1 Al和45钢的Hugoniot参数

Table 1. Hugoniot parameters of Al and 45 steel

Material	${{\;\rho _0}}$ /(g·cm^–3)	c₀/(km·s^–1)	λ
Al	2.785	5.328	1.338
45 steel	7.850	4.483	1.332

下载: 导出CSV

| 显示表格

PVDF压电传感器的输出电荷由电荷积分器采集，然后由示波器输出电压得到聚合物薄膜敏感压力-电荷密度标定曲线。

2.3 30 μm厚PVDF传感器动态冲击压缩实验标定

对于敏感单元厚度为 $30 \;{\text{μ}}{\rm{m}}$ 的PVDF传感器，在0.3～10.0 GPa压力区间由电荷法得到的实验结果如表2所示，其中：低压段采用气炮加载方式，高压段采用火炮加载方式，U_max为电压峰值，(Q/A)_max为电荷密度峰值，d₃₃为压电常数，P_r为剩余极化强度。

表 2 电荷法实验结果

Table 2. Experimental results by charge method

Exp.No.	Material of flyer	u₀/(m·s⁻¹)	p/GPa	U_max/V	(Q/A)_max/( ${\text{μ}}$ C·cm⁻²）	d₃₃/(pC·N⁻¹)	P_r/（ ${\text{μ}}$ C·cm⁻²）
0830-1	LY12 Al	48.434	0.361 2	0.598 67	0.665 2	12.0	8.279 7
0830-2	45 steel	83.384	1.473 0	1.097 90	1.219 9	10.8	7.528 9
0830-3	45 steel	109.557	1.942 8	1.367 27	1.519 2	10.5	7.878 1
0830-4	45 steel	129.536	2.303 8	1.433 01	1.592 2	15.5	7.987 1
0831-1	LY12 Al	126.221	0.950 6	0.813 58	0.904 0	15.2	10.301 7
1128-1	LY12 Al	47.137	0.351 5	0.459 01	0.510 0	15.1	7.806 1
1030	LY12 Al	628.531	5.395 0	1.636 01	1.817 8	15.6	6.934 0
1026	PMMA	1 004.704	2.426 1	1.120 07	1.244 5	15.3	7.832 2
1107	LY12 Al	970.275	8.069 3	1.883 43	2.092 7	11.4	7.216 3
1108	LY12 Al	1 250.441	10.725 5	2.357 62	2.619 6	12.8	8.287 2

下载: 导出CSV

| 显示表格

PVDF传感器在一定压力范围内具有线性度，但是在较宽压力范围内不具备线性响应特性。本研究中，为了提高标定曲线的精度，选取幂函数、线性函数和多项式分别进行拟合。结果表明，幂函数的拟合精度最高，为此选择幂函数拟合标定曲线。采用幂函数，分别对表2中0.35～2.30 GPa和2.4～10.7 GPa压力区间的实验数据进行拟合。低压段拟合结果如图8所示，拟合曲线方程为

图 8

$30 \;{\text{μ}}{\rm{m}}$ -PVDF传感器低压段标定曲线

Figure 8. Calibration curve of

$30 \;{\text{μ}}{\rm{m}}$ -PVDF in low pressure zone

下载: 全尺寸图片幻灯片

$\frac{Q}{A} = \left( {1.001\;66 \pm 0.031\;18} \right) \cdot {\sigma ^{\left( {0.565\;11 \pm 0.045\;67} \right)}}$

(4)

其相关指数为0.989，剩余标准差为6.611%。高压段采用 $\dfrac{Q}{A} = a{\sigma ^b}$ 函数进行拟合，如图9所示，拟合曲线方程为

图 9

$30 \;{\text{μ}}{\rm{m}}$ -PVDF传感器高压段标定曲线

Figure 9. Calibration curve of

$30 \;{\text{μ}}{\rm{m}}$ -PVDF in high pressure zone

下载: 全尺寸图片幻灯片

$\frac{Q}{A} = \left( {0.779\;60 \pm 0.091\;40} \right) \cdot {\sigma ^{\left( {0.498\;43 \pm 0.056\;81} \right)}}$

(5)

其相关指数为0.979，剩余标准差为10.015%。

3. 结　论

介绍了PVDF传感器研制过程的关键技术和工艺；在此基础上采用气炮加载装置进行了低压段和高压段一维平面应力标定方法研究，使用电荷测量模式，得到了敏感单元厚度为 $30 \;{\text{μ}}{\rm{m}}$ 的PVDF传感器在0.35～2.30 GPa和2.4～10.7 GPa压力范围的冲击压力与电荷密度关系。为了得到最佳的标定曲线，使用幂函数分别对低压段和高压段进行拟合，相关指数接近1，且剩余标准差最小。在标定实验过程中，当冲击压力大于10 GPa时，敏感单元厚度为 $30 \;{\text{μ}}{\rm{m}}$ 的PVDF传感器的响应时间为37.2 ns，说明传感器具有较高的频率响应特性。

图 1 战斗部结构示意图

Figure 1. Schematic diagram of warhead structure

下载: 全尺寸图片幻灯片

图 2 展开模式示意图

Figure 2. Schematic diagram of expansion modes

下载: 全尺寸图片幻灯片

图 3 战斗部的简化模型示意图（单位：mm）

Figure 3. Schematic diagram of warhead model (Unit: mm)

下载: 全尺寸图片幻灯片

图 4 选取的起爆点位置示意图

Figure 4. Schematic diagram of initiation point locations

下载: 全尺寸图片幻灯片

图 5 轴向展开角示意图

Figure 5. Schematic diagram of axial-expansion angle

下载: 全尺寸图片幻灯片

图 6 侧向展开角示意图

Figure 6. Schematic diagram of lateral-expansion angle

下载: 全尺寸图片幻灯片

图 7 不同起爆点位置下破片的飞散速度

Figure 7. Fragments’ mean velocity for different initiation points

下载: 全尺寸图片幻灯片

图 8 10 μs时不同起爆点位置下壳体所受压强分布云图

Figure 8. Pressure distribution on warhead case for different initiaion points at a time delay of 10 μs

下载: 全尺寸图片幻灯片

图 9 不同起爆点位置下的破片数统计

Figure 9. Fragment quantity for different initiation points

下载: 全尺寸图片幻灯片

图 10 100 μs时单元体破片的飞散情况

Figure 10. Fragment dispersion of unit constituting warhead at different initiation points at a time delay of 100 μs

下载: 全尺寸图片幻灯片

图 11 不同展开角下破片数统计

Figure 11. Fragment quantity at different axial-expansion angles

下载: 全尺寸图片幻灯片

图 12 不同轴向展开角下破片位置的散点图

Figure 12. Scatter plots of fragments’ locations at different axial-expansion angles

下载: 全尺寸图片幻灯片

图 13 不同轴向展开角下破片的平均速度及其Z轴的分量

Figure 13. Mean velocities and velocities in Z axis of fragments at different expansion angles

下载: 全尺寸图片幻灯片

图 14 破片密集区

Figure 14. Dense area of fragments

下载: 全尺寸图片幻灯片

图 15 不同轴向展开角下破片的分布情况

Figure 15. Fragment dispersions at different axial-expansion angles

下载: 全尺寸图片幻灯片

图 16 破片盲区锥顶角示意图

Figure 16. Schematic diagram of the cone angle for the zone without fragment

下载: 全尺寸图片幻灯片

图 17 不同轴向展开角下破片盲区的锥顶角

Figure 17. Cone angles of zone without fragment at different axial-expansion angles

下载: 全尺寸图片幻灯片

图 18 侧平面破片飞散角示意图

Figure 18. Schematic diagram of fragments’ dispersion angle in lateral plane

下载: 全尺寸图片幻灯片

图 19 不同轴向展开角度下破片的侧平面飞散角

Figure 19. Fragments’ dispersion angles in lateral plane for different axial-expansion angles

下载: 全尺寸图片幻灯片

图 20 不同轴向展开角下破片速度随周向角的分布情况

Figure 20. Fragments’ velocity distribution in versus of $\phi$ for different axial-expension angles

下载: 全尺寸图片幻灯片

图 21 不同轴向展开角下破片速度随侧平面飞散角的分布情况

Figure 21. Fragments’ velocity distributions in versus of $\delta$ for different axial-expansion angles

下载: 全尺寸图片幻灯片

图 22 侧向展开模式下破片飞散情况

Figure 22. Fragment dispersion in lateral-expanding mode

下载: 全尺寸图片幻灯片

图 23 侧向展开模式下破片位置散点图

Figure 23. Scatter plot of fragments’ locations in lateral-expanding mode

下载: 全尺寸图片幻灯片

图 24 侧向展开模式下破片速度分布

Figure 24. Velocity distribution of fragments in lateral-expanding mode

下载: 全尺寸图片幻灯片

表 1 材料模型

Table 1. Material models

Component	Material	Equation of state	Strength model	Failure model
Explosive	Comp. B	JWL	Hydro
Casing	Steel 4340	Linear	Johnson-Cook	Principal strain

下载: 导出CSV

表 2 Comp. B炸药材料参数

Table 2. Material parameters of Comp. B

$\,\rho$ /(g·cm⁻³)	D/(m·s⁻¹)	p_CJ/GPa	A/GPa	B/GPa	R₁	R₂	ω
1.717	7980	29.5	524.2	7.678	4.2	1.1	0.34

下载: 导出CSV

表 3 4340钢材料参数

Table 3. Material parameters of steel 4340

$\,\rho$ /(g·cm⁻³)	I_YS/GPa	H_C/GPa	H_E	S_RC	R_SR	B_M/GPa	S_M/GPa
7.83	0.792	0.51	0.26	0.014	1	159	81.8

下载: 导出CSV

表 4 密集区破片的平均速度

Table 4. Mean velocity of fragments in dense area

$\,\beta$ /(°)	Mean velocity/(m∙s⁻¹)
$\,\beta$ /(°)	d=10 mm	d=20 mm	d=30 mm	d=40 mm	d=50 mm	d=60 mm
60	1180	1182	1225	1287	1311	1323
75	1302	1332	1343	1366	1376	1401
90	1413	1491	1503	1515	1524	1527

下载: 导出CSV

表 5 侧平面飞散角在不同取值区间的破片占比

Table 5. Distribution of fragments at different intervals of dispersion angle in the lateral plane

$\,\beta$ /(°)	$\delta$ _min/(°)	$\delta$ _max/(°)	$x{_\delta}$ /%
$\,\beta$ /(°)	$\delta$ _min/(°)	$\delta$ _max/(°)	$\delta$ <10°	10° −20°	20° −30°	30° −40°	40° −50°	50° −60°	60° −70°	70° −80°	80° −90°	$\delta$ >90°
15	−3.54	33.47	33.26	41.34	22.25	3.15	0	0	0	0	0	0
30	4.24	55.05	6.92	7.14	45.76	31.03	4.02	5.13	0	0	0	0
45	16.76	69.07	0	4.34	5.84	28.35	43.29	13.42	4.76	0	0	0
60	30.98	84.57	0	0	0	13.12	13.62	42.16	25.45	3.34	2.31	0
75	49.99	91.84	0	0	0	0	10.60	21.66	27.19	29.95	9.45	1.15
90	47.75	101.00	0	0	0	0	0.17	8.95	25.99	23.75	29.95	11.19

下载: 导出CSV

表 6 飞散角 $\gamma$ 和 $\eta$ 在不同取值区间内的碎片数占比

Table 6. Fragment distributions at different intervals of dispersion angles $\gamma$ and $\eta$

$x{_\gamma}$ /%									$\gamma$ _max/(°)
0°−5°	5°−10°	10°−15°	15°−20°	20°−25°	25°−30°	30°−35°	35°−40°	40°−45°	$\gamma$ _max/(°)
57.68	23.41	9.88	3.29	2.32	0.73	0.85	0.61	1.23	44.24

$x{_\eta}$ /%									$\eta$ _max/(°)
0°−5°	5°−10°	10°−15°	15°−20°	20°−25°	25°−30°	30°−35°	35°−40°	40°−45°	$\eta$ _max/(°)
54.63	30.00	5.24	2.32	5.98	1.59	0.24	0	0	31.34

下载: 导出CSV

参考文献(13)

[1]	崔瀚, 张国新. 定向战斗部研究现状及展望 [J]. 飞航导弹, 2019(3): 84–89.
[2]	李记刚, 余文力, 王涛. 定向战斗部的研究现状及发展趋势 [J]. 飞航导弹, 2005(5): 25–29. doi: 10.3969/j.issn.1009-1319.2005.05.008
[3]	张天光. 美英定向战斗部的研究与应用 [J]. 航空兵器, 2002, 9(3): 38–41. doi: 10.3969/j.issn.1673-5048.2002.03.013
[4]	SAVAGE J C, O'NEAL J K, BROWN R A, et al. Powered low cost autonomous attack system: cooperative, autonomous, wide-area-search munitions with capability to serve as non-traditional ISR assets in a network-centric environment [C]//Proceedings of SPIE 5791, Laser Radar Technology and Applications X. Orlando: SPIE, 2005.
[5]	谭多望, 张振宇, 王志兵, 等. 爆炸变形战斗部模型试验研究 [J]. 爆炸与冲击, 2002, 22(3): 247–251. doi: 10.3321/j.issn:1001-1455.2002.03.010 TAN D W, ZHANG Z Y, WANG Z B, et al. Investigation of a detonatively deformable warhead [J]. Explosion and Shock Waves, 2002, 22(3): 247–251. doi: 10.3321/j.issn:1001-1455.2002.03.010
[6]	MILLER C W. Warhead having selectable axial effects: United States Patent Application, 20150033971 [P]. 2015-02-05.
[7]	凌琦, 何勇, 何源, 等. 扇形复合装药驱动破片定向飞散的数值模拟 [J]. 高压物理学报, 2017, 31(05): 557–565. LING Q, HE Y, HE Y, et al. Numerical simulation of directed scattering of fragments driven by sector-shaped double-layer charge [J]. Chinese Journal of High Pressure Physics, 2017, 31(05): 557–565.
[8]	洪晓文, 李伟兵, 李文彬, 等. 展开角度及起爆位置对轴向展开式定向战斗部性能的影响 [J]. 含能材料, 2018, 26(5): 383–389. HONG X W, LI W B, LI W B, et al. Effect of expanding angle and initiation position on the performance of axial-expanding directional warhead [J]. Chinese Journal of Energetic Materials, 2018, 26(5): 383–389.
[9]	陈闯, 杨丽, 阮凌飞, 等. 展开式EFP战斗部复合破片场成型的数值模拟 [J]. 火炸药学报, 2019, 42(3): 295–302. CHEN C, YANG L, RUAN L F, et al. Numerical simulation on composite fragment field formation of expansion EFP warhead [J]. Chinese Journal of Explosives & Propellants, 2019, 42(3): 295–302.
[10]	赵宇哲, 李健, 马天宝. 展开式定向战斗部展开过程实验研究 [J]. 高压物理学报, 2016, 30(2): 116–122. ZHAO Y Z, LI J, MA T B. Experiment on spread processes of the spreadable aimed warhead [J]. Chinese Journal of High Pressure Physics, 2016, 30(2): 116–122.
[11]	耿荻, 马天宝, 宁建国. 定向战斗部爆炸驱动规律研究 [J]. 高压物理学报, 2013, 27(5): 685–690. doi: 10.11858/gywlxb.2013.05.005 GENG D, MA T B, NING J G. Study on laws of explosive driven behaviors of aimed warhead [J]. Chinese Journal of High Pressure Physics, 2013, 27(5): 685–690. doi: 10.11858/gywlxb.2013.05.005
[12]	马征, 宁建国, 马天宝. 展开型定向战斗部展开过程的理论计算与数值模拟 [J]. 弹箭与制导学报, 2007, 27(1): 111–114. doi: 10.3969/j.issn.1673-9728.2007.01.033 MA Z, NING J G, MA T B. Dynamic analysis and numerical simulation for spread process of a evolvable aimed warhead [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2007, 27(1): 111–114. doi: 10.3969/j.issn.1673-9728.2007.01.033
[13]	AUTODYN, Autodyn user’s manual [M]. Release 19.0, 2019.

施引文献

资源附件(0)

访问统计

点击查看大图

图(25) / 表(6)

计量

文章访问数: 1764
HTML全文浏览量: 483
PDF下载量: 43

1. PVDF敏感单元设计
1.1 电极制作和极化
1.2 敏感单元封装
2. 动态冲击设计和标定
2.1 测试系统组件选择
2.2 PVDF传感器动态冲击压缩实验标定方法
2.3 30 μm厚PVDF传感器动态冲击压缩实验标定
3. 结　论

可变形定向破片战斗部在不同展开模式下的性能数值模拟

doi: 10.11858/gywlxb.20210887

作者简介: 杜伊杨（1999-），男，本科，主要从事战斗部设计研究. E-mail：duyiyang@njust.edu.cn

通讯作者: 李伟兵（1982-），男，博士，教授，博士生导师，主要从事高效毁伤技术研究.E-mail：njustlwb@163.com

计量

出版历程

Numerical Simulation of Deformable Directional Fragment Warhead Performance in Different Expansion Modes

1. PVDF敏感单元设计

1.1 电极制作和极化

1.2 敏感单元封装

2. 动态冲击设计和标定

2.1 测试系统组件选择

2.2 PVDF传感器动态冲击压缩实验标定方法

2.3 30 μm厚PVDF传感器动态冲击压缩实验标定

3. 结 论

计量

出版历程

目录

1. PVDF敏感单元设计

1.1 电极制作和极化

1.2 敏感单元封装

2. 动态冲击设计和标定

2.1 测试系统组件选择

2.2 PVDF传感器动态冲击压缩实验标定方法

2.3 30 μm厚PVDF传感器动态冲击压缩实验标定

3. 结 论

作者简介:
杜伊杨（1999-），男，本科，主要从事战斗部设计研究. E-mail：duyiyang@njust.edu.cn

通讯作者:
李伟兵（1982-），男，博士，教授，博士生导师，主要从事高效毁伤技术研究.E-mail：njustlwb@163.com

3. 结　论

3. 结　论