第39卷第4期 2 0 1 8年4月 兵 工 学 报 Vo1.39 No.4 Apr. 2018 ACTA ARMAMENTARII 单兵火箭弹炸高优化的数值模拟与试验研究 孙韬,张国伟,王一鸣,郭帅 (中北大学机电工程学院,山西太原030051) 摘要:为研究杀伤战斗部的空炸威力,建立战斗部有限元模型和破片着地坐标计算模型,利用 有限元分析软件ANSYS/LS—DYNA对战斗部动态起爆过程进行了数值模拟。根据模拟得到的破片 初速和飞散方向,结合破片在空气中的速度衰减公式以及对无防护有生目标的78 J动能杀伤标 准,计算出有效杀伤破片在地面的坐标,并对不同炸高条件下有效杀伤破片在地面目标靶上的分布 情况进行分析。研究结果表明:炸高为8 m左右时,战斗部满足杀伤威力指标要求;仿真结果与试 验结果基本吻合。 关键词:杀伤战斗部;空炸;数值模拟;试验研究 中图分类号:TJ415. 1 文献标志码:A 文章编号:1000-1093(2018)04-0681-07 DoI:10.3969/i.issn.1000—1093.2018.04.007 Numerical Simulation and Experimental Study on Optimization of Individual Rocket Burst Height SUN Tao,ZHANG Guo—wei,WANG Yi—ming,GUO Shuai (College of Mechatronic Engineering,North University of China,Taiyuan 030051,Shanxi,China) Abstract:A finite element model of warhead and a calculation model of ground coordinate of a fragment are established to give full play to the air burst power of antipersonnel warhead,and the dynamic detona- tion of warhead is numerically simulated by using ANSYS/LS—DYNA software.According to the initial velocity and flying direction of fragments in simulation,the velocity attenuation formula of fragment in air and the 78 J killing criterion of personnel,the coordinates of the effective lethal fragment on the ground are calculated,and the distribution of the effective lethal fragments on the ground targets at different burst heights is analyzed.The results show that the warhead meets the lethality requirements when burst height is about 8 m.The experimental results are in good agreement with the simulated results. Key words:antipersonnel warhead;air burst;numerical simulation;experimental study 0 引言 杀伤战斗部主要是在高能炸药爆炸作用下形成 大量高速破片,利用破片的高速撞击、引燃和引爆作 用毁伤目标 。对于地面目标,破片在地面的覆盖 收稿13期:2017-06-28 面积和分布密度是衡量战斗部威力的重要指标,而 炸高是影响破片在地面的覆盖面积和分布密度的重 要因素。随着炸高的增大,破片在地面的覆盖面积 增大、分布密度减小,因此如何获得最佳炸高、使战 斗部的杀伤威力得到充分发挥,具有一定的研究意 作者简介:孙韬(1992一),男,硕士研究生。E-mail:suntao1437@qq.com 通信作者:张国伟(1968一),男,教授,硕士生导师。E-mail:zha“ggu0wei1728@163.con 兵 工 学 报 义。对于炸高的研究,相关的理论 已经比较成 第39卷 速度为l30 m/s时破片的飞散情况进行数值模拟。 1.1模型结构及材料参数 熟,但在数值模拟方面还有些不足,目前仍停留在以 破片初速和破片飞散角 。衡量战斗部杀伤威力的 阶段,未考虑炸高的影响。 本文采用数值计算与半解析相结合的方法,对 如图2所示,战斗部仿真模型由空气、炸药、内 衬、预制破片及壳体组成。预制破片为直径5 mm 的钢球破片,质量0.5l2 g,共1 360个。图3所示为 预制破片的排列方式,共80层,每层17个,与实际 装填数相比,误差为1.5%.内衬壁厚2 mm,壳体壁 厚1.2 mm,战斗径比为2.5.起爆方式为头部 顶端中心起爆。 轴线l,j向 杀伤战斗部动态起爆 进行数值模拟,建立相应的 破片着地坐标的计算模型;通过破片飞行过程中的 速度衰减公式和无防护有生目标的78 J动能杀伤 判据,计算得到有效杀伤破片在地面的坐标,并对不 同炸高条件下有效杀伤破片在地面目标靶上的分布 情况进行分析;对比试验结果可知,二者数据吻合较 好,可为单兵火箭弹引信炸高的选择提供参考。 1 数值模拟 内衬 通过文献[8]获得的单兵火箭弹射击概况:初 速一般为120 m/s,密集杀伤半径大于15 m,射角过 大,火箭筒后面喷出的火焰会危及射手安全,射角过 小满足不了射程要求。采用质心运动方程组 对 火 ,●-●’-l-I-,,●-I- l_l-一,,--● l-I一, 破片 潍 腱勇苎 火箭弹不同射角、一定初速的外弹道进行仿真计算, 结合试验数据,选取与实际弹道轨迹拟合较好的外 弹道曲线,如图l所示。 图2 战斗部仿真结构模型 壳体 Fig.2 Simulation structure model of warhead 图l 外弹道曲线 Fig.1 Exterior ballistic curve 在火箭弹800~1 000 II1的考核射程范围以及合 理炸高范围(取2~10 Il1)内,其落姿(落角、落速)变 化情况如表1所示。表l中数据显示,落角和落速 的变化不大,因此,不考虑其变化对战斗部杀伤威力 图3 预制破片的排列方式 Fig.3 Arrangement of fragments 的影响,设定落速为130 m/s,落角为l0。,对战斗部 表1 火箭弹落姿 Fab.1 Rocket falling attitude 炸药使用TNT炸药,密度P=1.58 g/cm ,爆速 D =6 930 m/s,爆压P¨=20 GPa,采用HIGH—Ex— PLOSIVE—BURN模型及EOS—JWL状态方程,基本 形式 为 第4 l9】 火箭弹炸高优化的数值模拟与试验研究 p=以( 一 )e l+B( 一 )e I+ toE, (1) 式巾:, 为 力;E为爆轰产物的内能; 为爆轰产物 的相埘体积;A、B、R 、R:干¨60为待定常数,如表2 所永 、 表2 TN~I炸药的主要参数 I、aI1 2 Main palanloterS of TN…I 『人J{1ff¨做片采用钢质材料,壳体采Hj铝合金材 料,模J 为PI AsTIC—KINEMATIC,空气采胴MAT— NUI I {:炎, 、各材料主要参数如表3所示。 表3 各材料主要参数 I、Il}1.3 Main pal’LIIII ̄'tc'zS of mat( rials 1.2数值模拟结果分析 4所爪为不同[J1f刻破片的飞散情况。由图4 可知,破片飞敞较均匀,战 一部中部破片初速较大, 头部、尾邪破片初速较小 。 2 有效杀伤破片着地坐标计算 2.1破片着地坐标计算 盘l1 5所示,建j 坐标系Oxyz和坐标系Qdef, 图L}1,Ozx为射面,0¨为地面,0为原点,Q为炸点, h为炸I ,0 为落角,尺为破片到地面的飞行距离, 为破片 敞方位角, 为破片飞散方向 轴的 火角. ,为做片着地坐 点,,v为弹丸轴线与 轴的 交点,P轴,l Ozx内 、 … ,f】J髟△ /、_卡¨△OMN町得 MN!=QX + 一2 QNRcos , (2) |、 =OM:+D,v 一2 OM·ONcos 0, (3) 式中: = ; (4) R: I_: (5) 条纹等级/(irmI· J 1·IXs 1 【翊4破片飞敞情况 Fig.4 Dispersion of fragm ̄ nts 丽= ; (6) OM=htan‘D. (7) 由于同一层中的破片受力大小都相同,敞㈤一 【f】每个破片的 散方位角郜卡H¨,取每层破片的 其中一个破片进行研究,图6所示为沿战斗邸轴线 方向破片的初述变化和初速的轴向分速度变化曲 线。由图6可 ,忽略两端受稀疏波影响的破片,做 片初速基本分布在l 000~1 200 nl/s之间,具仃较 大的初始动能,初速的轴向分速度先减小、后增大 684 兵 工 学 报 第39卷 图5计算模型 Fig.5 Computational model 图6破片初速及初速轴向分速度沿轴向变化曲线 Fig.6 Variation curves of initial fragment velocity and its axial velocity along axial direction 如图7所示,根据破片初速 i及初速轴向分速 度 的关系: .COS|B= , (8) V 求得破片飞散方向与战斗部轴线方向的夹角,即破 片飞散方位角,结果如图8所示。 图8中舭段、LK段和KJ段分别对应战斗部 头部、中部和尾部,本文所建立的战斗部模型头部 和尾部为圆锥部、中部为圆柱部。由图8可知:圆 柱部破片飞散方位角差异小,破片集中飞散,杀伤 范围小,但分布密度大;圆锥部破片飞散方位角差 异大,破片飞散较分散,分布密度小,但杀伤范 围大。 由破片初速 ;与初速在e轴上的分速度 和轴 线方向分速度v 的关系可知: YiCOS COS +vdsin 0 , (9) 式中: 大于90。的破片飞向天空, 小于90。的破片 飞向地面,计算得到飞向地面的破片有679个。 轴线方向 破片飞散方向 弹丸 图7破片飞散方位角 Fig.7 Fragment dispersion azimuth angle 图8破片 散方位角沿轴向变化曲线 Fig.8 Variation curve of fragment dispersion azimuth angle along axial direction 将(4)式~(9)式代入(2)式和(3)式求出0,最 终得到破片的着地点坐标: =OMcos 0, (10) Y=±OMsin 0. (1 1) (11)式中的正负号由破片初速在/轴的分速度 ,的方向确定, ,为正则Y为正,vl为负则Y为负,由 于数量较多, 和 ,的模拟结果在此省略。 2.2有效杀伤破片判据 通过破片在空气中的速度衰减公式 ,计算破 片着地时的速度 : CDSpaR =Vie一 , (12) 式中:C。为破片阻力系数,对于球形破片C。=0.97; p 为空气密度;Js为破片迎风面积;m 为破片质量。 对于25 mm厚的松木靶,根据无防护有生目标 的78 J动能杀伤标准¨ ,由动能定理可得击穿松木 靶的破片最小着地速度: E =÷m , (13) 第4期 单兵火箭弹炸高优化的数值模拟‘j试验研究 式中:E 为破片动能。 指标要求。 由(I 2)式和(13)式计算得到满足击穿25 nlm 松木靶最小破片初速的破片,即有效杀伤破片。 2.3计算结果及分析 ’ 图9为不同炸高条件下有效杀伤破片在地面的 分布情况。以原点O为中心,将12块1 1"I1×2 m的 25 mm厚松木靶对称排布,图10为不同炸高条件下 目标靶上有效杀伤破片的分布情况。结合图10与 试验要求的杀伤威力指标(每块松木板上的有效穿 ’ : : 童 , d 童 曩., ∥ 毒. ● ●· ● . , · } 1. (口)炸高2111 (}lJ斜 I 4Il _孔数不小于3个)可知:炸高为2~6 m时,不能满足 每块松木板都有穿孔;炸高为10 nl时,松木板上的 穿孔数 能满足要求;炸高为8 Ill时,满足杀伤威力 。 + [ 鼍 ·: : H . i 。‘ ;=二 ● i :● : ,.1 ’ 善 t. 港 罐j . ● .二爱 曼 ● ● ● ‘。. 霹 · · 0 ∥ ‘ · - .h● - , )) , ( l1炸高2|I 1(1ll炸岛4I11 (flJ炸高6in ai Ihlrst heighl t)f 2 ril【bl Burst ight of4 n1 Ic)Burst height of6 I1 ^ - ●●k: · - . 。 .毒 , ‘ · 。 垂 直 : ‘ : :‘ :耋 . 一 鬻 墓 擘: + h:; ● 鼍 ? ' -‘ . 。 . 、- t - 嚣 ● i , 。I 繁 . t+● 、工 j ●● ¨● · : - ● ●0 ● ● r● - + t● I喜 .yy , (d}炸高8 nl ㈦炸高l0Ill (I”Bur ̄t hrighl of 8[11 (e1 Burst height of 10 131 图9 有效杀伤破片在地面的分布 Fig.9 D istribution of fragments on the ground 1 But st llPight t,f 2【1 J 1 lhll sl lll ighl f,f 4 I1 ’l:一 口 。≯ .●● __l ● ·- . I。 ㈧炸高6III 【III,k'r ̄'2i8Ill l。1 BUl‘st height f】f6 in 【(I1 Bursl ight nf 8 I1 - l . 一. ’ _J ‘I’‘ ●● Q .卜 一 』- ‘ L :. {¨炸高10Il 】J Bu rst hrighl ot 10 ill 图10 仃效杀伤破片住 标靶I 的分布 Fig.1 0 1)isi ̄’ibution of fragments f1lI the targets 3试验结果及分析 图l1所示为试验布置现场, 预计空炸点『F下 方放置了l2块l nl×2 m的25 illI11厚松小板 与静态试验 的 靶方式不同,动态试验根据 对各参量的分析,最终选取1 2块l nl×2 ill的 25 Inii1厚松木靶按图1 1所示排列,并设定杀伤威力 指标:每块松木板上有效穿孔数 小于3个。通过 设置不同的装填参数,对不同炸高条件下战斗部对 松木板的毁伤效果进行分析,最终确定有效炸高 范围。 本次共进行了24发弹的试验,炸高分布在 1.8~15.6 113之问。由试验结果町知:炸高为1.8~ 686 兵 工 学 报 第39卷 图1 1试验现场 Fig.1I Test site 6.2 m时,不能满足每个松木板都有穿孔;炸高为 8.9 15.6 m时,每个松木靶上都有穿孔,但穿孑L数 较少,不能满足杀伤威力指标要求;当炸高为6.8~ 8.5 m时,满足杀伤威力指标要求。部分数据如表4 所示。 表4 不同炸高下的毁伤效果 Tab.4 Damage effects at different heights of burst 炸高/m 毁伤效果 2.3 有穿孔的松木板共9个,穿孔数满足杀伤威力指标要求 4.3 有穿 L的松木板共10个,穿孔数满足杀伤威力指标要求 6.2 有穿孔的松木板共11个,穿孔数满足杀伤威力指标要求 7.9 每块松木板都有穿孔且穿孔数满足杀伤威力指标要求 每块松木板都有穿孔,其中4块松木板上的 穿孔数不满足杀伤威力指标要求 图12所示为对应表4中不同炸高条件下其中 一块松木板的穿孑L情况。图12中红色三角对应炸 高条件为7.9 m,共l6个;蓝色圆圈对应炸高条件 为10.4 m,共2个;炸高为2.3 m、4.3 m、6.2 m时, 此松木板上没有穿孔。 通过试验结果分析可知,满足杀伤威力指标的 有效炸高范围为6.8~8.5 m,试验结果与仿真结果 基本吻合。 4 结论 1)采用数值计算与半解析相结合的方法,对不 同炸高条件下有效杀伤破片在地面的分布情况进行 研究,得到最佳炸高在8 m左右,与试验结果吻合 较好。 2)针对不同的地面目标以及杀伤威力指标,有 效炸高范围不同。 3)研究结果可为单兵火箭弹引信炸高的选择 图l2 单块松木板穿孑L情况 Fig.1 2 Perforation of a piece of’pine WOO(] 提供参考。 参考文献(References) 『1] 李向东,钱建平,曹兵,等.弹药概论[M].北京:国防工业 出版社,2004:260—266. LI Xiang-dong,Q1AN Jian—ping,CAn Bing,et a1.Ammunition generality[M].Beijing:National Defense lndustq Press,2004: 260—266.(in Chinese) r11 。‘J 应国淼,董桂旭,肖师云.基于射击线技术的杀爆战斗部杀 伤面积计算[J].兵器装备工程学报,2016,37(I 1):20—23 YING Gun—miao.DONG Gut-XU,XIAO Shi-yun.Lethal area ral· culation of blast—fragmentation warhead based on firing line tech· nology[J].Journal of Ordnance Equipment Engineering,2016, 37(11):20—23.(in Chinese) [3] 金丽.赵捍东,曹红松,等.预制破片对地面人员的杀伤威力 分析计算[J].弹箭与制导学报,2006,26(4):157—159 JIN Li,ZHAO Han—dong,CAn Hong-song,et al,Analysis and calculation on the lethality of premade fragments to the personnel targets tm the ground[J].Journal of Projectiles,Rockets,Missiles and Guidance,2006,26(4):157—159.(in Chinese) [4] 杨云川,厉相宝,万仁毅.预制破片初速和飞散角的数值模 拟[j].弹箭与制导学报,2009.29(4):96—98. YANG Yun—chuan,LI Xiang—ban,WAN Ren—yi.Numerieal sinlU lation of the initial velocity and scattering angle of performed frag ments[J].Journal of Projectiles,Rockets,M issiles and Guid- ance,2009,29(4):96—98.(in Chinese) [5] 卢世杰.预制破片战斗部破片飞散的数值研究[D].太原:中 北大学,2015:I8—68. LU Shi-jie.Numerical study on damage effect of prefabricated fragment warhead[D].Taiyuan:Noah University of China, 2015:18—68.(in Chinese) 【6] 魏继锋,焦清介,吴成.预制破片战斗部试验与数值模拟研 究[J].弹箭与制导学报,2004,24(3):39—45. wEl Ji-feng。JIAO Qing-jie,wu Cheng.Experiment anti simula。 第4期 单兵火箭弹炸高优化的数值模拟与试验研究 687 tion on premade fragment warhead[J].Journal of Projectiles, Rockets.Missiles and Guidance,2004,24(3):39—45.(in Chinese) [7] Padraic E O,William W P,Charles E A.Dynamic launch process of performed fragments[J].Journal of Applied Physics,1988, 63(2):337—347. [8] 于庆金,钱伯成.轻武器射击实用手册[M].北京:出 版社,1999:305—590. YU Qing-jin,QIAN Bo-cheng.Practical manual on the shooting of light weapons[M].Beijing:The People’S Liberation Army Press, 1999:305—590.(in Chinese) [9] 徐明友.火箭外弹道学[M].哈尔滨:哈尔滨工业大学出版 社.2004:4—57. XU Ming—you.Rocket exterior ballistics[M].Harbin:Harbin In stitute of Technology Press,2004:4—57.(in Chinese) [1O] 高元浩.榴弹战斗部动爆条件下毁伤评估方法研究[D].太 原:中北大学,2016:31—46. GAO Yuan·hao.Damage assessment method of shrapnel warhead under dynamic explosion conditions[D].Taiyuan:Noah Univer- sity of China。2016:31—46.(in Chinese) [11] 时党勇,李裕春,张胜民.基于ANSYS/LS—DYNA8.1进行显 示动力分析[M].北京:清华大学出版社,2005. SHI Dang—yong。LI Yu-chun,ZHANG Sheng—min.Explicit dy— namic analysis using ANSYS/LS-DYNA8.1[M].Beijing: Tsinghua University Press,2005.(in Chinese) [12] Charon Y J.Estimation of velocity distribution of fragmenting warheads using a modified gurney method,AFIT/GAE/AA/79S 1[R].OH,US:Air Force Institute ofTechnology,1979:40— 42. [13] 魏惠之,朱鹤松,汪东晖,等.弹丸设计理论[M].北京:国 防工业出版社,1985:261—273. WEI Hui.zhi,ZHU He—song,WANG Dong-hui,et a1.Projeetile design theory[M].Beijing:National Defense Indust ̄Press, 1985:261—273.(in Chinese) [14] 张国伟.终点效应及其应用技术[M].北京:国防工业出版 社,2005:148—150. ZHANG Gun-wei.Terminal effect and its applied technology [M].Beijing:National Defense Industu Press,2005:148— 150.(in Chinese) [15] 高秀娟,金建祥,王林,等.弹丸杀伤半径计算方法研究 [J].弹箭与制导学报,2012,32(4):106—108. GAO Xiu-juan,JIN Jian-xiang,WANG Lin,et a1.Calculation method of killing radius[J].Journal of Projectiles,Rockets, Missiles and Guidance,2012,32(4):106—108.(in Chinese)