基于赋权灰色关联分析的药型罩的结构优化

吴波; 崔耀中; 蒙国往; 徐世祥; WU Bo; CUI Yao-zhong; MENG Guo-wang; XU Shi-xiang

网刊加载中。。。

使用Chrome浏览器效果最佳，继续浏览，你可能不会看到最佳的展示效果，

确定继续浏览么?

复制成功，请在其他浏览器进行阅读

基于赋权灰色关联分析的药型罩的结构优化 PDF

- ORCID：
吴波 ^1,2
✉
- ORCID：
崔耀中 ¹
✉
- ORCID：
蒙国往 ¹
- ORCID：
徐世祥 ¹

1. 广西大学土木建筑工程学院，广西南宁 530004； 2. 东华理工大学土木与建筑工程学院，江西南昌 330013

中图分类号： TJ51； TD235； O212.6

最近更新：2022-09-08

DOI：10.11943/CJEM2022126

摘要

为研究圆锥形药型罩的锥角、锥径、壁厚3个参数对聚能射流效果的影响程度并优化药型罩结构，设计三因素四水平的正交试验，采用ABAQUS软件对聚能射流破甲进行仿真模拟，利用赋权灰色关联法分析射流破甲正交试验结果，得到药型罩的锥角、锥径、壁厚与射流破甲效果的关联度。结果表明：锥角对聚能射流破甲效果影响最大，锥径、壁厚的影响次之。最佳聚能射流破甲效果的药型罩参数组合为：圆锥角61.2°、圆锥半径18 mm、壁厚1.05 mm。最佳参数组合下，聚能金属射流最大速度为5855 m∙s^-1，侵彻靶板深度为59.43 mm，侵彻靶板开口直径为8.24 mm。研究结果可为聚能射流的侵彻应用提供相关理论指导。

图文摘要

Effect of liner structural parameters on shaped charge jet penetration are demonstrated by numerical and weighted gray correlation analysis. A well‑performance shaped charge numerical model for the formation of jet and penetration. It is simple and quick to determine the optimal structural parameters of the liner， which can greatly improve the jet formation and penetration performance， and can be quickly applied to practical projects.

关键词

聚能射流; 仿真模拟; 正交试验; 参数优化; 灰色关联度

0 引言

聚能药包是一种具有非常高穿透能力的爆炸装置，在生产建设中有着广泛的应用，如兵器工程破甲与销毁废弃哑弹、岩土工程开挖及工民建筑切割拆除^［

1-6］。聚能药包利用了空穴效应原理^{［参考文献 20

百度学术}20］，使得爆轰产物向空穴轴线聚集，进而压垮药型罩，形成一股超高动能的金属射流。国内外研究者从理论、实验、数值对聚能射流展开了研究，聚能射流形成机制与侵彻性能是相关研究中最受关注两个方向。Tamer等^{［参考文献 7

百度学术}7］对Allison & Vitalli经验公式^{［参考文献 8

百度学术}8］进行修正，认为靶板屈服强度的提高可以显著降低聚能射流穿透深度，修正的经验公式和数值模拟结果和实验结果有着很好的一致性。徐文龙等^{［参考文献 9

百度学术}9］对超聚能射流形成过程机理做了理论研究，得到了超聚能射流速度及质量分布公式，并进行了数值模拟，对比了超聚能射流速度与传统聚能射流速度分布情况，结果与理论公式基本一致。王成等^{［参考文献 10

百度学术}10］运用数值模拟对一种Al药型罩环形聚能装药结构的长径比和起爆方式进行了分析，得到了此类型聚能装药结构较好的起爆直径。Dehestani等^{［参考文献 11

百度学术}11］采用欧拉法模拟了聚能射流形成与侵彻，并探究了炸高和药型罩壁厚对聚能射流穿透深度效率的影响。Zaki等^{［参考文献 12

百度学术}12］研究了在给定聚能射流结构几何形状条件下，药型罩的材料和炸药类型对聚能射流侵彻能力的影响。黄炳瑜等^{［参考文献 13

百度学术}13］开展了双层含能药型罩K装药射流形成及侵彻性能试验，研究发现Al/Ni‑Cu双层药型罩的侵彻性能相比传统Cu‑Cu药型罩大幅提高。

以上研究表明，药型罩参数设计是聚能装药结构设计中至关重要的一环。在固定的爆发器中，药型罩的锥径、锥角对聚能穴高度与药柱底部直径起决定性作用，直接影响有效炸药层的范围，对爆轰波传播与爆炸能量利用率有重要影响。药型罩压垮与射流形成过程也会受药型罩的结构影响，如果药型罩的锥角过小，则射流没有充足的发育时间，会抑制射流形成并降低射流的侵彻性能。药型罩的壁厚直接决定了形成金属射流的质量，但在固定的装药结构下，有效层炸药爆炸后能量的利用效率有一个极值，药型罩过薄，炸药能量不能完全利用，药型罩过厚，炸药能量不足以支持高速射流的形成。因此，要设计出最佳射流形成与侵彻性能的聚能装药结构，有必要优化药型罩的结构参数，确定其最佳参数组合。Dehestani^［

11］、Tamer^{［参考文献 7

百度学术}7］、Zaki^{［参考文献 12

百度学术}12］的研究中采用了圆锥形药型罩，虽然有考虑到药型罩的结构参数对射流侵彻性能的影响，但仅考虑到药型罩的某一个结构参数，独立分析此参数对聚能射流侵彻性能的影响，而药型罩的各参数综合作用下对聚能射流侵彻性能的影响还有待研究。

鉴于此，本研究设计了三因素四水平正交试验，采用ABAQUS软件^［

11］对聚能射流侵彻靶板进行数值模拟，运用赋权灰色关联分析对试验结果进行统计计算，研究圆锥形药型罩的锥角、锥径、壁厚三个参数组合作用下对聚能射流侵彻性能的影响，对一种实验用聚能射流装药结构的药型罩参数进行优化设计。

1 分析方法

1.1 灰色关联度法

灰色关联分析由邓聚龙^［

14］提出，此方法对于多目标响应问题有明显优势，其原理是根据曲线几何形状的贴合度来衡量数据序列间的关联程度。药型罩的结构参数优化问题是典型的多目标响应问题，通过计算出药型罩的各参数在不同水平下关联度并比较大小，找出最优药型罩的设计参数。灰色关联分析的具体步骤如下。

（1）设灰色关联评价系统由n个试验方案和m个评价指标组成，那么评价指标矩阵 $A = (a_{i j})_{n \times m}$ 表示为：

A = [\begin{matrix} a_{11} & a_{12} & \dots & a_{1 m} \\ a_{21} & a_{22} & \dots & a_{2 m} \\ ⋮ & ⋮ & ⋮ & ⋮ \\ a_{n 1} & a_{n 2} & \dots & a_{n m} \end{matrix}]

（1）

（2）将评价指标矩阵A归一化处理，归一化矩阵 $B = (b_{i j})_{n \times m}$ 中的元素由式（2）确定：

b_{i j} = \{\begin{matrix} \frac{a_{i j} - m i n (a_{1 j}, a_{2 j}, \dots a_{n j})}{m a x (a_{1 j}, a_{2 j}, \dots a_{n j}) - m i n (a_{1 j}, a_{2 j}, \dots a_{n j})} 正 指 标 \\ \frac{m a x (a_{1 j}, a_{2 j}, \dots a_{n j}) - a_{i j}}{m a x (a_{1 j}, a_{2 j}, \dots a_{n j}) - m i n (a_{1 j}, a_{2 j}, \dots a_{n j})} 逆 指 标 \end{matrix}

（2）

式中， $i = 1,2, \dots, n$ ； $j = 1,2, \dots, m$ 。

评价指标矩阵A经过式（2）处理后得到归一化矩阵B：

B = [\begin{matrix} b_{11} & b_{12} & \dots & b_{1 m} \\ b_{21} & b_{22} & \dots & b_{2 m} \\ ⋮ & ⋮ & ⋮ & ⋮ \\ b_{n 1} & b_{n 2} & \dots & b_{n m} \end{matrix}]

（3）

（3）计算关联系数矩阵。选取指标最大值作为参考方案，参考方案 $Z = [z_{1}, z_{2}, \dots, z_{m}]$ 由 $z_{j}$ 组成，其中 $z_{j} = m a x (b_{1 j}, b_{2 j}, \dots b_{n j}), j = 1,2, \dots, m$ 。由参考方案矩阵求出序列差 $Δ_{i j} = |z_{j} - b_{i j}|$ ，进而得到两级最小差值 $Δ_{m i n}$ 和两级最大差值 $Δ_{m a x}$ 。关联系数矩阵ξ中的元素 $ξ_{i j}$ 由式（4）求得：

\{\begin{array}{l} Δ_{m i n} = \underset{i}{m i n} \underset{j}{m i n} |z_{j} - b_{i j}| \\ Δ_{m a x} = \underset{i}{m a x} \underset{j}{m a x} |z_{j} - b_{i j}| \\ ζ_{i j} = \frac{Δ_{m i n} + ρ Δ_{m a x}}{Δ_{i j} + ρ Δ_{m a x}} \end{array}

（4）

式中， $ρ$ 为分辨系数，取 $ρ = 0.5$ 。

关联系数矩阵ξ为：

ξ = [\begin{matrix} ξ_{11} & ξ_{12} & \dots & ξ_{1 m} \\ ξ_{21} & ξ_{22} & \dots & ξ_{2 m} \\ ⋮ & ⋮ & ⋮ & ⋮ \\ ξ_{n 1} & ξ_{n 2} & \dots & ξ_{n m} \end{matrix}]

（5）

1.2 熵值赋权法

熵值法通过各指标的差异程度来确定信息量的大小，从而确定各指标的权重^［

15］。传统灰色关联分析对各指标采取无差异化处理，不能反应各指标的权重差异，运用熵值赋权法能对传统灰色关联分析进行客观赋权，增强灰色关联分析的客观性和可信度。熵值赋权法的步骤如下：

（1）正规化处理评价指标矩阵 $A = (a_{i j})_{n \times m}$ ，计算矩阵A第j项指标下第i个方案的指标比重：

a_{i j}^{'} = \frac{a_{i j}}{\sum_{i = 1}^{n} a_{i j}}

（6）

式中， $i = 1,2, \dots, n$ ； $j = 1,2, \dots, m$ 。

（2）计算评价指标的信息熵值 $e_{j}$ ：

e_{j} = - \frac{1}{l n m} \sum_{i = 1}^{n} a_{i j}^{'} l n a_{i j}^{'}

（7）

（3）由各指标的信息熵值计算其客观权重 $α_{j} (j = 1,2, \dots, m)$ ：

α_{j} = \frac{1 - e_{j}}{\sum_{j = 1}^{m} (1 - e_{j})}

（8）

1.3 灰色关联赋权分析

基于正交试验、熵值赋权‑灰色关联分析的聚能射流药型罩的结构参数优化步骤如下：

（1）以药型罩的锥角、锥径、壁厚作为试验因素，以金属射流最大速度、金属射流侵彻靶板后最大完整长度、靶板侵彻开口半径、靶板侵彻长度作为评价指标，建立三因素四水平的正交试验，构建评价指标矩阵A。

（2）根据公式（2）~（5），求出关联系数矩阵 $ξ$ 。

（3）根据公式（6）~（8），计算各评价指标的客观权重值α_j，得到客观权重矩阵α。

（4）结合关联系数矩阵和客观权重矩阵，由式（9）计算灰色关联度矩阵ω：

ω = ξ α

（9）

2 仿真模拟

Dehestani等^［

11］运用ABAQUS对顶部单点起爆、圆锥形药型罩聚能射流装药结构侵彻靶板进行数值模拟，其实验室试验结果与数值模拟结果对比如图1所示。本研究采用的数值模拟方法主要参考Dehestani等^{［参考文献 11

百度学术}11］的研究。

a. measurement result of experiment

b. measurement result of simulation

图1 聚能射流侵彻靶板空腔对比^［

11］

Fig.1 Cavity created due to jet penetration into the target by using （a） experimental testing，（b） simulation process^［

11］

由图1可知，试验结果与数值模拟结果的空腔长度量测偏差为3.64%，小于工程一般可接受偏差5%^［

16］，认为实验与模拟的一致性较好，且二者空腔形状由外至里相似，故研究所采用的数值模拟方法正确。

2.1 聚能装置仿真模型

选取一种实验用的聚能射流装置，聚能装药结构1/2几何模型如图2所示。装药截面形状为轴对称型，药型罩对称布置于装药轴线。聚能装药结构具有几何对称性，采用1/18模型，降低运算时间。仿真模型如图3所示，x轴方向为聚能装药结构轴线方向。聚能装药结构起爆方式为单点起爆，起爆点位于装药顶部中央，实验中炸高一般为0.8倍~1.05倍药型罩口径，炸高为30 mm满足研究所有试验的有利炸高要求。假设模型中所有材料均为连续均匀介质，整个运算过程为绝热过程。仿真分析采用动力‑显示分析步，时长95 μs，环境温度300 K，所有材料均采用欧拉算法。

图2 1/2几何模型

Fig.2 1/2 model

图3 仿真模型

Fig.3 Simulation model

2.2 材料模型

炸药选择梯恩梯（TNT）炸药，由JWL状态方程^［

17］描述炸药爆炸后压力与体积之间的关系：

p = A (1 - \frac{ω}{R_{1} V}) e^{- R_{1} V} + B (1 - \frac{ω}{R_{2} V}) e^{- R_{2} V} + \frac{ω E}{V}

（10）

式中，p为爆轰压力，GPa；V为初始相对体积，A、B、ω、R₁、R₂为状态方程基本参数，E为初始单位体积内能，炸药参数如表1^［

18］所示。

表1 炸药材料参数^{［参考文献 18

百度学术}18］

Table 1 Explosive material parameters^{［参考文献 18

百度学术}18］

explosive

/kg∙m^-3

V_D

/m∙s^-1

/GPa

R₁

R₂

/J∙m^-3

TNT

1630

6930

374

3.75

4.15

0.9

0.35

6.0×10⁹

Note： 1） ρ is the density of explosive； 2） V_D is the detonation velocity； 3） p is the detonation pressure； 4） A，B，R₁，R₂，ω are the material constants measured for the explosive； 5） E is the detonation energy density.

由于药型罩的压垮过程、射流的形成过程、射流侵彻靶板的过程及爆发器外壳破碎的过程均在非常短的时间内完成，所有的金属材料在此过程中将发生很大的塑性变形，金属材料的力学行为将非常依赖于应变速率和温度。因此，选用Johnson‑cook材料模型^［

19］描述药型罩、装药壳体、靶板的力学行为。Johnson‑cook材料模型^{［参考文献 11

百度学术}11］是Mises各向同性塑性模型中的一种特殊类型，在此模型中，材料对应变速率和温度的行为依赖关系的本构方程为：

σ = [A + B {({\bar{ε}}^{p l})}^{n}] \cdot [1 + C l n (\frac{{\dot{\bar{ε}}}^{p l}}{{\dot{\bar{ε}}}_{0}})] \cdot [1 - {\bar{θ}}^{m}]

（11）

式中， ${\bar{ε}}^{p l}$ 为塑性应变， ${\dot{\bar{ε}}}^{p l}$ 为塑性应变率，A、B、C、m、n为测定的材料参数。参数 ${\bar{θ}}^{m}$ 的定义如下：

\bar{θ} = \{\begin{matrix} 0 & θ < θ_{t r a n s i t i o n} \\ \frac{θ - θ_{t r a n s i t i o n}}{θ_{m e l t} - θ_{t r a n s i t i o n}} & θ_{t r a n s i t o n} < θ < θ_{m e l t} \\ 1 & θ > θ_{m e l t} \end{matrix}

（12）

式中， $θ$ 是材料温度，K； $θ_{m e l t}$ 是材料熔点温度，K； $θ_{t r a n s i t i o n}$ 是材料转变温度，K。

药型罩的材料采用99.99%铜，装药壳体的材料采用AISI‑1045钢，靶板的材料采用含氮硬钢，三种材料的Johnson‑cook模型参数见表2^［

11］。

表2 药型罩、壳体、靶板的Johnson‑Cook材料模型参数^{［参考文献 11

百度学术}11］

Table 2 Johnson‑Cook material model parameters for the liner， case and target^{［参考文献 11

百度学术}11］

material	ρ / kg∙m^‑3	A / MPa	B / Mpa	n	m	c	E / GPa	ν	θ_transition / K	θ_melt / K
copper	8960	1425	168	0.12	1	0.032	125	0.34	293	1356
steel AISI‑1045	7850	553	600	0.234	1	0.013	210	0.3	293	1733
steel high‑nitrogen	7800	1500	2600	0.82	1.05	0.0009	200	0.3	293	1800

Note： 1） ρ is the density of explosive； 2） A is the quasi static yield stress； 3） B is the strain hardening coefficient； 4） n， m， c are material parameters； 5） E is the shear modulus； 6） ν is the Poisson ratio.

3 正交试验及结果

理想的聚能装药结构应具有以下特征：药型罩的质量大、有效炸药层能量充足、射流发育充分。圆锥形药型罩的结构参数对聚能装药结构的上述3个特征影响很大。汪旭光等^［

20］依据定常、理想、不可压缩的流体力学理论推导计算，认为药型罩的小锥角利于提高聚能射流速度，进而影响射流侵彻靶板穿孔直径与穿孔深度；药型罩的直径影响了有效炸药层与非有效炸药分布，前者参与形成聚能射流，后者保证有效炸药层稳定爆轰；药型罩的壁厚对射流性能和聚能威力会产生显著的影响，最佳壁厚随着药型罩的材料、锥角、直径、有无外壳而变化。聚能射流最大速度与最大完整长度是判断聚能射流形成优劣的2个指标；侵彻靶板穿孔直径与穿孔深度用以衡量聚能射流侵彻效果，穿孔直径越大、深度越长，射流杀伤效果越好。

为测定在固定爆发器、炸高30 mm条件下，圆锥角62.4°，圆锥半径17 mm，壁厚0.95 mm的铜制圆锥形药型罩在聚能装药爆轰作用下金属射流侵彻靶板杀伤效果并优化药型罩的结构提高射流形成与侵彻性能，设计正交试验，采取熵值赋权‑灰色关联度分析数值模拟结果，预期在药型罩参数的较小区间内确定最优参数组合，能较大程度提高射流形成与侵彻性能，结果与结论可快速应用于工程。

对上述实验用聚能射流装药结构的药型罩的结构参数设计正交试验，以药型罩的锥角、锥径、壁厚作为试验因素，以金属射流最大速度、金属射流侵彻靶板后最大完整长度、侵彻靶板开口半径、侵彻靶板长度作为评价指标，建立L16（4³）正交试验表，试验具体安排如表3所示，3因素4水平正交实验，一共16组试验。该实验用聚能射流装药结构的药型罩的结构参数为：圆锥角62.4°、圆锥直径34 mm、壁厚0.95 mm，正交试验表的水平取值已考虑上述3个特征并包含了原药型罩参数。由于所有的仿真模拟试验中，不存在射流因头、尾的速度差过大而断开的情况，且杵体在接触靶板后会对射流的完整形状造成直接影响，所以射流最大完整长度取杵体接触靶板前射流的长度。

表3 药型罩的结构正交试验表

Table 3 Orthogonal test table for liner

level	test factors
level	one half of the conic angle / （°）	cone radius / mm	thickness of liner / mm
1	30.6	16.5	0.90
2	31.2	17.0	0.95
3	31.8	17.5	1.00
4	32.4	18.0	1.05

试验结果表明，16组正交试验中，试验4与试验13射流形成及侵彻靶板的结果对比最明显，这两组试验的射流形成与侵彻靶板的仿真结果如图4所示。选取试验4的全过程速度场图用于展现聚能射流侵彻靶板的全过程，如图5所示。运用ABAQUS的后处理功能对四项评价指标进行测量，仿真模拟计算结果的测量结果如表4所示。

a. simulation results of test 4

b. simulation results of test 13

图4 射流形成与侵彻靶板的仿真结果

Fig.4 Simulation results of jet formation and penetration

图5 试验4聚能射流全过程速度场图

Fig.5 Velocity profile at different times of shaped charge jet formation in test 4

表4 正交试验方案及结果

Table 4 Numerical results of orthogonal tests

test number	test factor			evaluation index
test number	α/2/（°）	R/mm	T/mm	V_max/m∙s^-1	L_max/mm	A/mm	B/mm
1	30.6	16.5	0.90	5391	20.56	5.04	55.05
2	30.6	17.0	0.95	5536	20.68	5.04	55.95
3	30.6	17.5	1.00	5570	20.80	5.25	57.70
4	30.6	18.0	1.05	5855	20.75	4.12	59.43
5	31.2	16.5	0.95	5194	19.60	4.90	52.15
6	31.2	17.0	0.90	5338	19.98	4.76	53.85
7	31.2	17.5	1.05	5134	19.75	5.78	50.90
8	31.2	18.0	1.00	5287	20.05	5.66	54.10
9	31.8	16.5	1.00	5005	19.20	5.15	49.05
10	31.8	17.0	1.05	4910	19.25	5.30	48.40
11	31.8	17.5	0.90	5310	19.70	4.99	53.67
12	31.8	18.0	0.95	5291	19.55	5.17	53.10
13	32.4	16.5	1.05	4847	18.53	5.29	47.06
14	32.4	17.0	1.00	4947	19.05	5.16	49.18
15	32.4	17.5	0.95	5018	19.23	5.11	50.58
16	32.4	18.0	0.90	5274	19.58	5.06	53.13

Note： 1） α/2 is the half of the conic angle； 2） R is the conic radius； 3） T is the thickness of the liner； 4） V_max is the maximum velocity of the jet； 5） L_max is the maximum complete length of the jet after it has penetrated the target； 6） A is the opening radius of the target； 7） B is the length of the jet penetrating the target.

4 灰色关联分析

根据表4的正交试验结果，得到4项评价指标的16组试验值，进而得出评价指标矩阵A。由式（2）对评价指标矩阵A进行归一化处理，得到归一化矩阵B，再计算得出序列差矩阵，并通过式（4）计算得到关联系数矩阵ξ。

B = [\begin{matrix} 0.540 & 0.894 & 0.554 & 0.646 \\ 0.684 & 0.947 & 0.554 & 0.719 \\ 0.717 & 1.000 & 0.681 & 0.860 \\ 1.000 & 0.978 & 0.000 & 1.000 \\ 0.344 & 0.471 & 0.470 & 0.411 \\ 0.487 & 0.639 & 0.386 & 0.549 \\ 0.285 & 0.537 & 1.000 & 0.310 \\ 0.437 & 0.670 & 0.928 & 0.569 \\ 0.157 & 0.295 & 0.620 & 0.161 \\ 0.063 & 0.537 & 0.711 & 0.108 \\ 0.459 & 0.515 & 0.524 & 0.534 \\ 0.440 & 0.449 & 0.633 & 0.488 \\ 0.000 & 0.000 & 0.705 & 0.000 \\ 0.099 & 0.229 & 0.627 & 0.171 \\ 0.170 & 0.308 & 0.596 & 0.285 \\ 0.424 & 0.463 & 0.566 & 0.491 \end{matrix}]

ξ = [\begin{array}{r} 0.521 & 0.825 & 0.529 & 0.585 \\ 0.621 & 0.904 & 0.529 & 0.640 \\ 0.639 & 1.000 & 0.610 & 0.781 \\ 1.000 & 0.958 & 0.333 & 1.000 \\ 0.433 & 0.486 & 0.485 & 0.459 \\ 0.494 & 0.581 & 0.449 & 0.526 \\ 0.411 & 0.519 & 1.000 & 0.420 \\ 0.470 & 0.602 & 0.874 & 0.537 \\ 0.372 & 0.415 & 0.568 & 0.373 \\ 0.348 & 0.519 & 0.634 & 0.359 \\ 0.480 & 0.508 & 0.512 & 0.518 \\ 0.472 & 0.476 & 0.576 & 0.494 \\ 0.333 & 0.333 & 0.629 & 0.333 \\ 0.357 & 0.393 & 0.572 & 0.376 \\ 0.376 & 0.420 & 0.553 & 0.411 \\ 0.465 & 0.482 & 0.535 & 0.495 \end{array}]

通过式（6）对评价指标矩阵A进行正规化处理，再根据式（7）计算出各个评价指标的信息熵值 $e_{j}$ ，最后由各指标的信息熵计算得到对应的客观权重矩阵 $α = {[0.2523 0.2524 0.2489 0.2463]}^{T}$ ，根据关联系数矩阵和客观权重矩阵计算结果，通过式（9）式计算出16个试验的灰色关联度并进行排序，结果如表5所示。灰色关联度值越大，表明试验结果越好，试验因素取值组合越好。

表5 目标函数灰色关联度

Table 5 Grey correlation degree of objective function

test number	test factor			Grey correlation degree	sorting
test number	α/2/（°）	R/mm	T/mm	Grey correlation degree	sorting
1	30.6	16.5	0.90	0.616	5
2	30.6	17.0	0.95	0.672	3
3	30.6	17.5	1.00	0.758	2
4	30.6	18.0	1.05	0.823	1
5	31.2	16.5	0.95	0.466	11
6	31.2	17.0	0.90	0.512	7
7	31.2	17.5	1.05	0.587	6
8	31.2	18.0	1.00	0.620	4
9	31.8	16.5	1.00	0.432	14
10	31.8	17.0	1.05	0.465	12
11	31.8	17.5	0.90	0.505	8
12	31.8	18.0	0.95	0.504	9
13	32.4	16.5	1.05	0.407	16
14	32.4	17.0	1.00	0.425	15
15	32.4	17.5	0.95	0.440	13
16	32.4	18.0	0.90	0.494	10

Note： 1） α/2 is the one half of the conic angle； 2） R is the conic radius； 3） T is the thickness of the liner.

由目标函数的灰色关联度可知，试验1和试验8的灰色关联度分别为0.616，0.620、试验5和试验10的灰色关联度分别为0.466，0.465、试验11和试验12的灰色关联度分别为0.505，0.504，其灰色关联度相近，说明药型罩的结构参数不同，对射流形成与侵彻性能的影响也可能会相近。试验4的灰色关联度值最大，试验13的灰色关联度值最小，二者相差0.416，说明不同的药型罩的结构参数对射流形成与侵彻性能的影响可能会很大，研究最佳的药型罩的结构参数组合对提升射流形成与侵彻性能具有重要意义。

根据表5可得到三种试验因素在4个水平下的4个关联度值，平均后得到其平均关联度，进而确定药型罩的最优结构参数组合。三种试验因素不同水平的平均关联度见表6，平均关联度柱状图见图6。

表6 药型罩的结构参数4水平平均关联度

Table 6 4 level average correlation degree of structural parameters of the liner

test level	α/2 /（°）	correlation degree	R /mm	correlation degree	T /mm	correlation degree
level 1	30.6	0.717	16.5	0.480	0.90	0.532
level 2	31.2	0.546	17.0	0.519	0.95	0.521
level 3	31.8	0.477	17.5	0.573	1.00	0.557
level 4	32.4	0.442	18.0	0.610	1.05	0.571

Note： 1） α/2 is the half of the conic angle； 2） R is the conic radius； 3） T is the thickness of the liner.

图6 关联度柱状图

Fig.6 Histogram of correlation degree

由表6和图6可知，药型罩的锥角灰色关联度最大值为0.717，对应锥角61.2°；药型罩圆锥半径的灰色关联度最大值为0.610，对应圆锥半径18 mm；药型罩壁厚的灰色关联度最大值为0.571，对应壁厚为1.05 mm。因此，锥角对射流形成与侵彻性能的影响最大，锥径、壁厚的影响次之。聚能射流结构药型罩的锥角、锥径、壁厚的最佳参数组合为：锥角61.2°，圆锥半径18 mm，壁厚1.05 mm，试验4中的药型罩参数为：二分之一锥角30.6°，圆锥半径18 mm，壁厚1.05 mm，试验4的药型罩参数组合恰巧是最优参数组合。

由图6关联度柱状图与表6水平平均关联度可知：药型罩的锥角在61.2°~64.8°范围内，随着锥角增大，关联度呈下降趋势，与射流形成和侵彻性能呈负相关。药型罩的直径在33~36 mm范围内，随直径增大，关联度呈上升趋势，与射流形成和侵彻性能呈正相关。药型罩壁厚在0.9~1.05 mm范围内，随壁厚增大，关联度先下降再上升，壁厚为0.95 mm的灰色关联度值最小，壁厚为1.05 mm的灰色关联度值最大，高于壁厚为0.90 mm的灰色关联度值。因此在一定范围内，适当减小药型罩的锥角，增大药型罩的直径与壁厚，可提高聚能装药结构的射流形成与侵彻性能。

然而正交试验‑赋权灰色关联分析仍存在一定的局限性。由于试验数量的限制，试验因素的水平取值只能控制在一定的范围且个数不宜过多，因此得到的结果不具有全范围的代表性。若要得到更大范围内的药型罩的结构参数取值对聚能射流形成与侵彻性能的影响，可结合仿真模拟与大数据、机器学习技术。仿真模拟提供丰富的试验数据，大数据、机器学习技术对试验数据进行深度学习处理，二者互辅，能较好预测出药型罩的结构参数改变对聚能射流效果的影响。

5 结论

设计三因素四水平正交试验，运用ABAQUS软件建立聚能金属射流破甲模型，采用欧拉法计算分析16个试验的金属射流形成与侵彻靶板的过程，结合灰色关联分析和赋权法对聚能金属射流结构的药型罩参数进行优化，得出如下结论：

（1）在三因素四水平的装药结构顶部单点起爆试验条件下，圆锥形药型罩的最佳结构参数组合为：药型罩圆锥角61.2°，圆锥半径18 mm，药型罩壁厚1.05 mm。最佳参数组合下，聚能金属射流最大速度为5855 m∙s^-1，较试验平均值高11.65%；侵彻靶板深度为59.43 mm，较试验平均值高12.75%；射流最大完整长度为20.75 mm，较试验平均值高4.8%。侵彻靶板直径为8.24 mm，较试验平均值低19.37%。

（2）聚能金属射流药型罩的结构参数中，锥角对聚能金属射流的形成与侵彻性能影响最大，锥径、壁厚的影响次之，聚能金属射流装药结构的药型罩参数设计应重点考虑圆锥形药型罩的锥角。在一定范围内，适当减小药型罩的锥角，增大药型罩的直径和壁厚，可显著提高射流形成与侵彻性能。

（3）采用正交试验‑赋权灰色关联分析方法用于聚能金属射流药型罩的结构参数优化，能直观清晰地表现出药型罩参数对聚能射流形成与侵彻性能的影响程度大小。其计算简便，所需试验量较少，有较好的适用性与应用性。

参考文献

郗文博，同剑，赵云涛，等. 便携式非接触聚能射流引爆器销毁废旧弹药［J］. 工程爆破， 2020， 26（3）： 75-78. [百度学术]

XI Wen‑bo， TONG Jian， ZHAO Yun‑tao， et al. Destruction of abandoned ammunition by portable non‑contact shaped charge jet detonator［J］. Engineering Blasting， 2020， 26（3）： 75-78. [百度学术]

李晓杰，江德安，闫鸿浩，等. 大型聚能切割器对浅海废弃油井的爆破拆除［J］. 辽宁工程技术大学学报（自然科学版）， 2008， 27（S1）： 136-138. [百度学术]

LI Xiao‑jie， JIANG Dean， YAN Hong‑hao， et al. Blasting demolition of oil well discarded by larger linear cumulative cutter in shallow sea［J］. Journal of Liaoning Technical University（Natural Science）， 2008， 27（S1）： 136-138. [百度学术]

付艳恕，俞建魁. 公交车起火逃生通道的聚能射流开辟与扩充［J］. 高压物理学报， 2017， 31（2）： 162-168. [百度学术]

FU Yan‑shu， YU Jian‑kui. Exploring and enlarging escape way by shaped charge for ignited bus［J］. Chinese Journal of High Pressure Physics， 2017， 31（2）： 162-168. [百度学术]

熊祖钊，王成，钟冬望. 线性聚能装药切割岩石的数值模拟［J］. 火炸药学报， 2004， 27（4）： 29-32. [百度学术]

XIONG Zu‑zhao， WANG Cheng， ZHONG Dong‑wang. Numerical simulation of linear shaped charge on rock cutting［J］. Chinese Journal of Explosive & Propellants，2004，27（4）：29-32. [百度学术]

任新见，李林，汪剑辉. 一种销毁哑弹的聚能射流设计［J］. 爆破， 2008， 25（4）： 82-84. [百度学术]

REN Xin‑jian， LI Lin， WANG Jian‑hui. Design of cumulative jet to eliminate unexploded bomb［J］. BLASTING， 2008， 25（4）： 82-84. [百度学术]

安二峰，沈兆武，周听清，等. 一种新型聚能破甲战斗部及其发展趋势探讨［J］. 中国工程科学， 2004， 6（6）： 85-90. [百度学术]

AN Er‑feng， SHEN Zhao‑wu， ZHOU Ting‑qing， et al. A new model shaped charge warhead and the discussion of it's evolution［J］. Engineering Science， 2004， 6（6）： 85-90. [百度学术]

Tamer Elshenawy， Ahmed Elbeih， Q M Li. Influence of target strength on the penetration depth of shaped charge jets into RHA target［J］. International Journal of Mechanical Sciences， 2018， 136： 234-242. [百度学术]

ALLISON F E， VITALI R. A new method of computing penetration variables for shaped charge jet［R］. Ballistic Research Laboratory Report， No.1184， 1963. [百度学术]

徐文龙，王成，徐斌. 超聚能射流形成过程机理研究［J］. 兵工学报， 2018， 39（2）： 261-268. [百度学术]

XU Wen‑long， WANG Cheng， XU Bin. Investigation of hyper shaped charge jet formation theory［J］. Acta Armamentarii， 2018， 39（2）： 261-268. [百度学术]

王成，邓涛，徐文龙. 铝药型罩环形聚能射流的数值模拟［J］. 北京理工大学学报， 2019， 39（12）： 1211-1218. [百度学术]

WANG Cheng， DENG Tao， XU Wen‑long. Numerical simulation of annular shaped charge with Al liner［J］. Transactions of Beijing Institute of Technology， 2019， 39（12）： 1211-1218. [百度学术]

Dehestani P， Fathi AR， Mohammadi Daniali HR. Numerical study of the stand‑off distance and liner thickness effect on the penetration depth efficiency of shaped charge process［J］. Proceedings of the Institution of Mechanical Engineers， Part C： Journal of Mechanical Engineering Science， 2019， 233（3）： 1-10. [百度学术]

S Zaki， Emad Uddin， B Rashid， et al. Effect of liner material and explosive type on penetration effectiveness of shaped charge［J］. Proceedings of the Institution of Mechanical Engineers， Part L： Journal of Materials： Design and Applications， 2019， 233（7）： 1375-1383. [百度学术]

黄炳瑜，熊玮，张先锋，等. 双层含能药型罩K装药射流成型及侵彻性能试验［J］. 含能材料， 2021， 29（2）： 149-156. [百度学术]

HUANG Bing‑yu， XIONG Wei， ZHANG Xian‑feng， et al. Experimental study on jet formation and penetration performance of double⁃layered reactive liners with k‑charge［J］. Chinese Journal of Energetic Materials（Hanneng Cailiao）， 2021， 29（2）： 149-156. [百度学术]

邓聚龙.灰色控制系统［J］.华中工学院学报，1982，10（3）：9-18. [百度学术]

DENG Ju‑long. The grey control system［J］. Journal of Huazhong University of Science and Technology， 1982， 10（3）： 9-18. [百度学术]

吴波，路明，雷领，等. 基于正交试验、组合赋权‑灰色关联度法研究隧道施工最优方案［J］. 中国安全生产科学技术， 2019， 15（8）： 124-130. [百度学术]

WU Bo， LU Ming， LEI Ling， et al. Research on optimal scheme of tunnel construction based on orthogonal test and combined weighting‑grey correlation method［J］. Journal of Safety Science and Technology， 2019， 15（8）： 124-130. [百度学术]

WU Hao， HU Feng， FANG Qin. A comparative study for the impact performance of shaped charge jet on UHPC targets ［J］. Defence Technology， 2019， 15（4）： 506-518. [百度学术]

崔宇，马宏昊，沈兆武，等. 爆炸二极管的数值模拟设计［J］. 含能材料， 2017， 25（4）： 266-272. [百度学术]

CUI Yu， MA Hong‑hao， SHEN Zhao‑wu， et al. Numerical simulation and design of explosive diode［J］. Chinese Journal of Energetic Materials（Hanneng Cailiao）， 2017， 25（4）： 266-272. [百度学术]

程兵，汪海波，宗琦. 基于SPH‑FEM耦合法切缝药包爆破机理数值模拟［J］. 含能材料， 2020， 28（4）： 300-307. [百度学术]

CHENG Bing， WANG Hai Bo， ZONG Qi. Numerical simulation on blasting mechanism of slotted cartridge based on coupled SPH⁃FEM algorithm［J］. Chinese Journal of Energetic Materials（Hanneng Cailiao）， 2020， 28（4）： 300-307. [百度学术]

钱石川，甘强，任志伟，等. 爆炸箔起爆器发火阈值影响因素的数值模拟［J］. 含能材料， 2018， 26（3）： 248-254. [百度学术]

QIAN Shi‑chuan， GAN Qiang， REN Zhi‑we， et al. Numerical simulation of the factors affecting the ignition threshold of an exploding foil initiator［J］. Chinese Journal of Energetic Materials（Hanneng Cailiao）， 2018， 26（3）： 248-254. [百度学术]

汪旭光，郑炳旭，张正忠，等. 爆破手册［M］. 北京：冶金工业出版社， 2010： 868-870. [百度学术]

WANG Xu‑guang， ZHENG Bing‑xu， ZHANG Zheng‑zhong， et al. Handbook of blasting［M］. Beijing： Metallurgical Industry Press， 2010： 868-870. [百度学术]

首页

期刊简介

编委会

作者专区

审稿专区

道德声明

联系我们

进入EMF

基于赋权灰色关联分析的药型罩的结构优化 PDF

摘要

图文摘要

关键词

0 引言

1 分析方法

1.1 灰色关联度法

1.2 熵值赋权法

1.3 灰色关联赋权分析

2 仿真模拟

2.1 聚能装置仿真模型

2.2 材料模型

3 正交试验及结果

4 灰色关联分析

5 结论

参考文献

基于赋权灰色关联分析的药型罩的结构优化 PDF

摘要

图文摘要

关键词

0 引 言

1 分析方法

1.1 灰色关联度法

1.2 熵值赋权法

1.3 灰色关联赋权分析

2 仿真模拟

2.1 聚能装置仿真模型

2.2 材料模型

3 正交试验及结果

4 灰色关联分析

5 结 论

参考文献

0 引言

5 结论