药型罩锥角对破片飞散的影响.pdf
doi10. 3969/ j. issn. 1001- 8352. 2020. 01. 004 药型罩锥角对破片飞散的影响 * 韩文斌 张国伟 邵 彬 曹 鹏 中北大学机电工程学院(山西太原,030051) [摘 要] 为了研究药型罩锥角对某小长径比多功能战斗部破片飞散的影响,利用 ANSYS/ LS- DYNA 模拟了不同 的药型罩锥角条件下破片飞散速度和破片飞散角的变化。 仿真结果表明破片的整体飞散速度随着药型罩锥角的 增大而增加,锥角每增大 10 ,破片整体飞散速度增加 2%,因此破片杀伤威力得到提高;破片飞散角度随着药型罩 锥角的增加会略微增大,提升破片对目标的杀伤面积,且破片飞散角度会整体向上偏移。 后续通过静爆试验验证 了数值模拟结果的准确性。 可为小长径比多功能战斗部在调整破片杀伤威力方面的研究提供参考。 [关键词] 小长径比;药型罩锥角;破片飞散;数值模拟 [分类号] TJ410 Effect of Cone Angle of Shaped Liner on Fragment Dispersion HAN Wenbin, ZHANG Guowei, SHAO Bin, CAO Peng College of Mechanical and Electrical Engineering, North University of China (Shanxi Taiyuan, 030051) [ABSTRACT] In order to study the influence of the cone angle of shaped liner on the fragment dispersion of a small aspect ratio multi- functional warhead, ANSYS/ LS- DYNA was used to simulate the variation of the fragment dispersion velocity and the fragment dispersion angle under different cone angles of the liner.The results show that the overall disper- sion velocity of fragments increases with the increase of the cone angle of the liner.When the cone angle increases by 10 the overall dispersal velocity of fragments increases by 2%, and the lethality of fragments increases.The fragment disper- sion angle increases slightly with the increase of the cone angle of the liner, which enlarges the damage area of the fragment to the target, and the fragment dispersion angle will shift upward as a whole.The accuracy of numerical simulation results was verified by static explosion test. The research results can provide a reference for adjusting the fragment killing power of the multi- functional warhead with small aspect ratio. [KEYWORDS] small aspect ratio; cone angle of liner; fragment dispersion; numerical simulation 引言 现代战场上,面对防护能力日益增加的装甲目 标,兼具破甲、杀伤等功能于一体的多功能战斗部的 破甲威力和破片杀伤威力亟待提高。 战斗部爆炸 后,破片在空间的飞散分布是确定破片杀伤作用场 必须研究的一个重要问题 [1] 。 破片初速、飞散角是 衡量战斗部性能的重要参数 [2] 。 国内,随着仿真软 件的日益成熟,在破片飞散领域展开了较为广泛的 数值模拟研究,并取得了一定的研究成果。 李松楠 等 [3] 对不同起爆点位置对破片飞散方向的影响进 行了数值模拟研究;刘琛等 [4] 研究了偏心起爆方式 对棱柱型定向战斗部破片飞散影响的数值模拟仿 真。 国外对于破片飞散速度和角度也有广泛的研 究。 David 等 [5] 对布列赫的二维格尼(Gurney)方程 进行改进,提高了数值模拟过程中初始破片速度的 精度;Pankaj 等 [6] 对单约束的破片发射装置进行参 数化研究,得出装药量以及装药长径比对破片飞散 速度和角度的影响。 对于小长径比多功能战斗部而言,锥角对破片 飞散情况的影响并没有人做过具体的分析。 本文 中,在小长径比多功能战斗部的其他参量都不变的 情况下,改变药型罩锥角的大小,进行数值模拟分 析,找出药型罩锥角的变化对破片飞散的影响规律。 第49 卷 第 1 期 爆 破 器 材 Vol. 49 No. 1 2020 年 2 月 Explosive Materials Feb. 2020 * 收稿日期2019- 07- 04 第一作者韩文斌(1995 - ),男,硕士,主要从事弹药毁伤评估的研究。 E- mail1115083842@qq. com 通信作者张国伟(1968 - ),男,硕导,教授,主要从事弹药毁伤评估的研究。 E- mailzhangguowei1728@nuc. edu. cn 万方数据 1 破片飞散初速和飞散角的理论计算 多功能战斗部兼具破甲、杀伤等功能于一体,因 此战斗部威力包含破甲威力和破片杀伤威力两部 分。 其中,破片杀伤威力从破片飞散速度、破片飞散 角两方面来衡量。 1. 1 破片飞散初速 破片的飞散初速 v0是破片的主要性能指标之 一,是衡量杀伤战斗部毁伤威力的重要参数。 理论 分析计算时,假设破片与爆轰产物只在径向飞散, 忽略轴向运动;装药的能量全部消耗在爆轰产物与 破片上,忽略弹体破裂吸收的能量;破片初速相等, 忽略周围空气等介质吸收的能量。 根据能量平衡方 程,推导得出 Gurney 公式 v0= 2E m M 1 +m 2M 。(1) 式中m 为炸药质量;M 为壳体质量; 2E为 Gurney 常数或 Gurney 比能。 该方法由于使用简便且准确 度较高而得到了广泛应用。 每种炸药都有一定的 Gurney 比能,通过炸药理 论与实验得出其与爆速 D 的关系为 2E =D 8 。(2) 因此,已知炸药的爆速就可以求出 Gurney 比 能。 公式的前提条件是瞬时爆轰且不考虑轴向稀疏 波的影响。 式(1)适用于圆柱形弹,且只能给出最 大速度的近似值,不能给出速度的分布。 对于预制 破片的弹体,式(1)修正后为 v0=D 2 2m 5M 。(3) 1. 2 破片飞散角 破片飞散角是指战斗部爆炸后,在战斗部轴线 平面内,以质心为顶点所做的包含有效破片 90%的 锥角,记为 θ,也就是破片飞散图中(图 1)包含有效 破片 90%的两线之间的夹角。 θ = 40Qv D( 1 α- 1 2 ) 1 2 (cosθ2- cosθ1) 。(4) 式中Qv为炸药爆热;D 为炸药爆速;α为弹丸装填 系数;θ1和 θ2分别为爆轰波到达弹丸前、后端时爆 轰波的法向与壳体表面的夹角。 1 - 战斗部中心轴线;2 - 破片飞散方向;3 - 水平线。 图1 破片飞散角度简图 Fig. 1 Scattering angle sketch of fragment 2 战斗部的仿真建模 2. 1 模型的建立 通常而言,小长径比战斗部是指长径比小于0. 8 的战斗部结构。 本文中所研究的小口径多功能战斗 部的长径比为0. 7,具有破片杀伤和破甲侵彻功能, 轴对称,采用单层轴向预制钨球破片交错排列摆放。 整个模型由主装药、药型罩(采用球锥结合药型 罩)、壳体、预制钨球破片和空气组成。 如图 2 所 示。 其中,装药直径为 40. 7 mm,装药高度为 42. 0 mm,炸高为85. 0 mm,壳体厚度为 1. 5 mm;预制破 片为直径为 3. 0 mm 的钨球,共 11 层,每层 14 个。 如图3 所示。 起爆方式为装药顶端面中心起爆。 1 - 空气;2 - 起爆点;3 - 壳体; 4 - 破片;5 - 炸药;6 - 药型罩。 图 2 战斗部仿真结构模型 Fig. 2 Simulation structure model of warhead 图 3 破片摆放示意图 Fig. 3 Fragment display sketch 912020 年2 月 药型罩锥角对破片飞散的影响 韩文斌,等 万方数据 由于该战斗部为轴对称结构,故建立1/ 4 模型, 单位制为 cm- g- μ s。 炸药、药型罩和空气采用多物 质 ALE 算法;壳体、破片采用 Lagrange 算法,二者采 用 Constrained_Lagrange_In_Solid 进行流固耦合,这 样不容易出现网格畸变。 所有单元均采用8 节点实 体单元 Solid164 [7- 8] 。 2. 2 材料模型与参数 聚能装药采用 8701 炸药,密度 ρ = 1. 7 g/ cm 3, 爆速 De= 8 100 m/ s,爆压 pCJ= 30 GPa,采用 High_ Explosive_Burn 材料模型和 EOS_JWL 状态方程进行 描述。 此方程是 Jones- Wilkins- Lee 研究得到的,并 假定爆轰前沿以常速率传播。 JWL 状态方程定义 压力为相对体积 V 和单位体积的初始能量 E 的函 数 p = A(1 - ω R1V )e -R1V + B(1 - ω R2V)e -R2V + ω E V 。 (5) 式中参数 ω 、A、B、R1和 R2为表征炸药特性的常 数,该状态方程能很好地描述高能炸药,因为它在涉 及结构金属加速度的应用中可以确定炸药的爆轰压 力 [9- 12] 。 材料具体参数如表1 所示。 表1 8701 炸药的参数 Tab. 1 Material parameters of 8701 explosive A/ GPa B/ GPa ωVR1膊R2 56F. 406∫. 800]. 3610梃. 004\. 001珑. 30 药型罩的材料为紫铜,选用 Johnson_Cook 模型 和 EOS_JWL 状态方程进行描述;壳体采用钢材质, 用 Plastic_Kinematic 模型描述;空气采用空物质材 料 Null 描述,对应的状态方程为多线性状态方程; 预制破片采用钨材质,用 Plastic_Kinematic 模型描 述。 各个材料的主要参数如表2 所示 [13- 16] 。 表2 各个材料参数 Tab. 2 Material parameters 材料ρ/ (g cm - 3) 弹性模量泊松比 紫铜8u. 96846贩 空气1祆. 293 10 - 3 钢7u. 830210屯0妸. 28 钨18构. 03\. 450 . 3 3 数值模拟分析 3. 1 初始模型仿真结果 初始战斗部中,药型罩采用球锥结合结构,锥角 为 100 ,为便于统计破片飞散速度和破片飞散角, 破片层远离药型罩的层定义为第 1 层破片,靠近药 型罩的最下层破片定义为第 n 层破片。 选取战斗部 起爆点所在方向与战斗部中心水平线夹角为正,战 斗部药型罩所在方向与战斗部水平线夹角为负。 经过仿真计算后,得到初始仿真结果,见表 3。 表3 初始仿真结果杀伤威力统计 Tab. 3 Statistics of killing power of initial simulation results 破片层数破片飞散速度/ (m s - 1) 破片飞散角/ ( ) 1侣526434拻. 9 2侣639416拻. 6 3侣736410拻. 4 4侣82846{. 7 5侣88443{. 1 6侣92240{. 6 7侣9444- 1*. 7 8侣9484- 3*. 0 9侣9304- 5*. 2 10ж8814- 9*. 7 11ж7824- 23览. 1 3. 2 杀伤威力仿真结果 由于本次仿真实验主要研究该多功能战斗部在 金属射流成型较好的情况下药型罩锥角对破片飞散 的影响规律,因此,根据参考文献选取了药型罩锥角 为60 、70 、80 、90 、110 、120 6 个药型罩锥角进 行建模,数值模拟计算后,不同药型罩锥角杀伤威力 统计结果如图4、图5 所示。 预制破片战斗部在炸药起爆后,金属预制破片 在爆轰波和爆轰产物的作用下发生膨胀变形,并开 始逐渐加速,直到爆轰产物膨胀速度相对于预制破 片的速度可以忽略为止。 由于破片速度衰减程度远 图 4 破片飞散速度与破片层数的关系 Fig. 4 Relationship between fragment dispersion velocity and fragment layer number 02 爆 破 器 材 第 49 卷第 1 期 万方数据 图5 破片飞散角度与药型罩锥角的关系 Fig. 5 Relationship between fragment dispersion angle and cone angle of liner 低于爆轰产物和冲击波的衰减速度,因此,破片是远 距离目标的主要杀伤元素。 炸药的比冲量是指炸药质量与弹体质量的比 值。 一般来说,炸药质量增加,炸药威力大,炸药的 爆速和比冲量也大。 图4 为不同药型罩锥角对破片飞散速度和破片 层数的分布曲线。 从图 4 可以看出,随着药型罩锥 角的增加,破片整体飞散速度增加。 前 7 层破片的 飞散速度逐渐增加,第 8 层到最后一层破片飞散速 度逐渐减小;锥角每增加 10 ,破片整体飞散速度增 加2%。 其原因是当药型罩锥角增大且战斗部口径 不变时,药型罩的罩高降低,战斗部装药高度提高, 从而增加了炸药的总体质量 m。 从理论计算所得到 破片初速公式可以看出,在壳体质量 M 保持不变 时,药型罩锥角增大,炸药质量 m 和炸药爆速 D 得 到提高,破片整体初速增加。 从炸药的作用机理来 看,随着炸药质量和爆速的提高,产生的爆轰产物增 加,爆轰波爆速提高,作用在单个破片上的球形爆轰 波爆速和压强增加,促使破片的整体速度提高。 图5 为不同药型罩锥角对破片飞散角影响的情 况。 从图5 可以看出,随着药型罩锥角增大,破片飞 散角中线整体向上偏移,破片飞散角度增大,破片散 布更加分散。 其原因是药型罩锥角增大,需要增加 炸药质量;一般来说,炸药质量增加,炸药威力大,炸 药的爆速和比冲量也大。 由上述的破片飞散角度公 式可以得出,随着炸药爆速 D 的增加,破片整体飞 散角度增加。 由于受到边界效应的影响,第 1 层和最后一层 破片的飞散速度和飞散方向角变化过大,分布相对 比较离散。 因此,破片在飞散边缘的分布密度偏低。 4 试验验证 4. 1 试验方法 图6 所示为试验布置现场。 此次试验为静态打 靶,试验状态为在战斗部下方放置 100 mm 装甲 钢,用以测试破甲威力是否满足战绩指标,其中,炸 高为85 mm;在距离战斗部轴线 3 m 处有一个面积 为1m 2m、厚度为3 mm 的 Q235 钢靶,用以测试破 片飞散的杀伤威力和飞散分布,3 m 处的 Q235 钢靶 主要是模拟该小长径比多功能战斗部产生的破片对 轻型装甲车辆顶装甲的杀伤威力。 4. 2 试验结果 本次试验为锥角为 100 的小长径比多功能战 斗部打靶试验。 前面的试验显示,当锥角选取 100 时满足破甲威力的战绩指标。 试验结果为破片上 靶数量每层2 ~ 3 枚,共 9 层;上靶23 枚破片,其中 卡4 枚,透19 枚。 破片飞散角上线为16. 9 ;破片飞 散角中线为0. 8 ;破片飞散角下线为 - 9. 5 ;其中, 第1 层和最后一层破片由于破片飞散角度过大,均 未着靶。 剩余层数,破片飞散角度的计算与仿真结 果基本吻合,见表4。 因为战斗部中线高度在1. 5 m 左右,钢靶高度上线为2. 0 m,破片飞散时战斗部中 线上半部分的破片只有0. 5 m高度可以上靶,所以 (a)战斗部实物图 (b)试验现场 (c)破片的着靶情况 图6 试验布置现场及破片着靶 Fig. 6 Testing layout and fragments in target 122020 年2 月 药型罩锥角对破片飞散的影响 韩文斌,等 万方数据 表 4 试验数据与仿真结果对比 Tab. 4 Comparison of test and simulation results 破片层数试验破片飞散角/ ( )仿真破片飞散角/ ( ) 1媼34�. 9 2媼16D. 916�. 6 3媼10D. 710�. 4 4媼7.. 16珑. 7 5媼3.. 43珑. 1 6媼0.. 80珑. 6 7媼- 1\. 6- 1 . 7 8媼- 3\. 2- 3 . 0 9媼- 5\. 1- 5 . 2 10 - 9\. 5- 9 . 7 11 - 23,. 1 为能更好地统计破片飞散,上靶数量只统计钢靶上 方1. 0 m 处(即以1. 5 m 为对称中线,上下 0. 5 m, 共1. 0 m 高度)破片。 经对比发现,该初始仿真结果与上述试验结果 基本吻合,证明此仿真具有参考价值。 5 结论 1)对于小长径比多功能战斗部来说,在其他结 构不变、只改变药型罩锥角的前提下,药型罩锥角的 增大会提高破片的整体飞散速度。 2)对于小长径比多功能战斗部来说,在满足破 甲威力战绩指标的前提下,可通过适当增加药型罩 的锥角来增加破片的飞散方向角,增大破片杀伤范 围,进一步提高战斗部杀伤威力。 参 考 文 献 [1] 黄广炎,刘沛清,冯顺山.基于战斗部微圆柱分析的破 片飞散特性研究 [J].兵工学报,2010,31(增 1)215- 218. HUANG G Y, LIU P Q, FENG S S. Research on disper- sion characteristic of fragment based on micro- column analyses for warhead [J].Acta Armamentarii, 2010,31 (Suppl. 1)215- 218. [2] 孙韬,张国伟,齐蕾.预制破片飞散规律的数值模拟研 究 [J].兵器装备工程学报,2017,38(12)45- 47. SUN T,ZHANG G W,QI L. Numerical simulation study on the rule of prefabricated fragments flying[J]. Journal of Ordnance Equipment Engineering,2017,38(12)45- 47. [3] 李松楠,张国伟,崔小杰,等.起爆点位置对破片飞散 方向的影响研究 [J].兵器装备工程学报,2018,39 (11)49- 53. LI S N, ZHANG G W, CUI X J, et al.Influence of the location of the detonating point on the dispersion direction of fragments[J].Journal of Ordnance Equipment Engi- neering,2018,39(11)49- 53. [4] 刘琛,李元,李燕华,等. 偏心起爆方式对棱柱形定向 战斗部破片飞散规律的影响[J]. 含能材料,2017,25 (1)63- 68. LIU C, LI Y, LI Y H,et al. Influence of eccentric initia- tion ways on fragment dispersion rule of prismatic aimable warhead[ J].Chinese Journal of Energetic Materials, 2017,25(1)63- 68. [5] FELIX D, COLWILL I, STIPIDIS E.Real- time calcula- tion of fragment velocity for cylindrical warheads [J]. Defence Technology,2019,15(3)264- 271. [6] CHOUDHA P K, KUMARASWAMY A, DHOTE K D. Parametric study of single confined fragment launch ex- plosive device [J].Defence Technology,2019,15(2) 179- 185. [7] 张世林,黄德雨,樊雪锋,等.轴向驱动反应破片的作 用过程 [J].爆破器材,2016,45(4)49- 53. ZHANG S L, HUANG D Y, FAN X F, et al.Research on the mechanism of axial driving response fragments[J]. Explosive Materials, 2016,45(4)49- 53. [8] 赵飞扬,王志军,尹建平,等.预制破片/ MEFP 组合式 战斗部数值仿真[J].兵器装备工程学报,2017,38 (12)52- 56. ZHAO F Y, WANG Z J, YIN J P, et al.Numerical simulation of prefabricated fragment/ MEFP combat war- head[J].Journal of Ordnance Equipment Engineering, 2017,38(12)52- 56. [9] 印立魁,候秀成,赵太勇,等.对破片飞散参量计算模 型的分析和评价 [J].弹箭与制导学报,2017,37(4) 55- 59. YIN L K, HOU X C, ZHAO T Y, et al. Analysis and as- sessment on computational model of fragment emission pa- rameters [ J].Journal of Projectiles, Rockets, Missiles and Guidance,2017,37(4)55- 59. [10] 周亚萍,王凤英,刘天生,等.基于 LS- DYNA 对高速 运动战斗部中预制破片飞散规律的数值模拟 [J]. 科技通报,2017,33(5)94- 97,101. ZHOU Y P,WANG F Y,LIU T S,et al.Numerical simulation of preformed fragment distributing regularity with the warheads in high- speed movement based on LS- DYNA[J].Bulletin of Science and Technology, 2017, 33(5)94- 97,101. [11] 凌琦,何勇,何源. 定向战斗部破片飞散的数值模拟 与试验研究 [J]. 振动与冲击,2017,36(3)234- 241. 22 爆 破 器 材 第 49 卷第 1 期 万方数据 LING Q,HE Y,HE Y.Numerical simulation and tests for fragments dispersion of an aimed warhead[J].Jour- nal of Vibration and Shock, 2017,36(3)234- 241. [12] 卢世杰. 预制破片战斗部破片飞散的数值研究 [D]. 太原中北大学,2015. LU S J. Numerical study on damage effect of prefabrica- ted fragment warhead[D].TaiyuanNorth University of China,2015. [13] 印立魁,蒋建伟. 多层球形预制破片战斗部破片初速 场的计算模型 [J]. 含能材料,2014,22(3)300- 305. YIN L K, JIANG J W.Calculation model of initial ve- locity field on multilayered spherical fragments warhead [J].Chinese Journal of Energetic Materials, 2014,22 (3)300- 305. [14] XU Y, WANG X S, KONG C Z, et al.Effect of initia- tion manners on the scattering characteristics of semi- preformed fragment warhead [J].Defence Technology, 2018,14(5)578- 584. [15] AMES R G.A standardized evaluation technique for re- active warhead fragments[C] / / 23rd International Sym- posium on Ballistics.Tarragona International Ballistics Committee,2007 16- 20. [16] JEOUNG H A, KHANH H N, YONG B P, et al. A nu- merical study of the high- velocity impact response of a composite laminate using LS- DYNA [J].International Journal of Aeronautical and Space Sciences,2010,11 (3) 221- 226. (上接第17 页) [8] 陈玉琼,叶红齐,刘辉. 表面活性剂预处理铝粉包覆聚 合物后的分散和耐蚀性能[J]. 材料保护,2010,43 (6)16- 17,34. CHEN Y Q,YE H Q,LIU H. Dispersion and anticorrosion behavior of flaky aluminum powders pretreated with sur- factants and encapsulated with polymer[J]. Materials Pro- tection,2010,43(6)16- 17,34. [9]王敏. TMPTMA/ AA/ St 共聚包覆片状铝粉的研究[D]. 长沙中南大学,2008. [10] 陶俊,王晓峰,王彩玲,等. 含氟聚合物包覆铝粉的分 子动力学计算[J]. 火炸药学报,2014,37(6)24- 30. TAO J,WANG X F,WANG C L,et al. Molecular dy- namic calculation of aluminum power coated by fluoro- polymer[J]. Chinese Journal of Explosives & Propel- lants,2014,37(6)24- 30. [11] 王慧心,任慧,闫涛,等. 聚乙烯醇原位包覆铝粉结构 表征及活性铝含量测定[J]. 兵工学报,2019,40(7) 1373- 1380. WANG H X,REN H,YAN T,et al. Micro- structure and active aluminum content of aluminum powder in situ coa- ted by polyvinyl alcohol[J]. Acta Armamentarii,2019, 40(7)1373- 1380. [12] 郝洁,丁昂,黄敬晖,等. 高活性纳米铝粉 HTPB- TDI 包覆组装与表征[J]. 稀有金属,2018,42(2)168- 174. HAO J,DING A,HUANG J H,et al. Cladding assembly and characterization of HTPB- TDI coated high activity nano- aluminum powder [J]. Chinese Journal of Rare Metals,2018,42(2)168- 174. [13] ZHANG K,FAN J H,WU J Y,et al. Preparation and characterization of nano- aluminium/ polystyrene encap- sulation particales[J]. High Power Laser and Particle Beams,2015,27(2)024145. 张凯,范敬辉,吴菊英,等. 纳米铝粉/ 聚苯乙烯包覆 粒子的制备及表征[J]. 强激光与粒子束,2015,27 (2)024145. [14] 王敏,叶红齐,唐新德. 三元共聚法包覆片状铝粉 [J].粉末冶金材料科学与工程,2007,12(5)301- 304. WANG M,YE H Q,TANG X D. Study on flaky alumi- num powder coated by ternary copolymerization[J]. Ma- terials Science and Engineering of Powder Metallurgy, 2007,12(5)301- 304. [15] 肖春,祝青,谢虓,等. PDA 包覆铝粉及其在 HTPB 中 的分散稳定性[J]. 火炸药学报,2017,40(3)60- 63, 76. XIAO C,ZHU Q,XIE H,et al. Polydopamine coated on aluminum powers and its disperse stability in HTPB[J]. Chinese Journal of Explosives & Propellants,2017,40 (3)60- 63,76. 322020 年2 月 药型罩锥角对破片飞散的影响 韩文斌,等 万方数据