战斗部破片对目标打击迹线的计算方法.pdf

第3 0 卷第4 期2 0 1 0 年7 月爆炸与冲击E X P I，O S I()NA N DS H O C KW A V E SV 0 1 3 0，N o 4J u l，2 0 1 0文章编号：1 0 0 1 1 4 5 5(2 0 1 0)0 4 0 4 1 3 0 6战斗部破片对目标打击迹线的计算方法黄广炎h 2，冯顺山2，刘沛清1(1 北京航空航天大学航空科学与工程学院，北京1 0 0 0 8 3；2 北京理工大学爆炸科学与技术国家重点实验室，北京1 0 0 0 8 1)摘要：提出了破片对目标打击迹线概念，通过建立预制或半预制破片战斗部对目标杀伤作用场参数的计算模型，以及在空间坐标系下破片各威力特性参数的相互转换得到了一种破片打击迹线的计算方法，并基于V C 和M A T L A B 联合编程实现了战斗部破片打击迹线作用场的分布计算。通过破片对目标打击迹线可定量分析战斗部所有破片的整个空间分布，并直观描述每个破片的飞散方向、能量衰减及对目标打击能力变化情况，从而可精确地评价破片对目标不同要害位置处的打击效果。该杀伤作用场参数计算模型和破片对目标打击迹线计算方法可用于战斗部起爆姿态选择和对不同目标杀伤效果的威力评价，为战斗部方案设计和战斗部终点毁伤效应分析计算提供了新的参考依据和技术途径。关键词：爆炸力学；打击迹线；M A T L A B；破片；战斗部；V C中图分类号：0 3 8 3 1国标学科代码：1 3 0 3 5 2 0文献标志码：Al引言破片是杀伤或杀爆战斗部毁伤的重要元素，战斗部爆炸后的破片初速、飞行轨迹、能量衰减等杀伤特性参数非常复杂。以往基本上根据实验判断毁伤效果，但攻击目标的差异、战斗部装药结构、起爆条件等相关因素的变化都必须要用不同的实验方案，耗费人力物力巨大。计算机仿真技术的发展，带来了战斗部毁伤领域新的发展契机。在系统直观描述破片飞行弹道及杀伤能力方面，钱立新等采用破片场仿真方法对大面积杀爆战斗部作用过程进行了威力评定，钱立新等 2 基于防空战斗部提出了破片打击迹线的概念，研究了可描述杀伤、聚焦、多P 和连续杆战斗部杀伤威力的打击迹线计算模型，蒋建伟等o 为模拟战斗部破片的宏观可控性和微观随机性，建立了描述广义杀伤战斗部起爆及破片飞散的射击迹线仿真模型，冯顺山等 43 对战斗部破片初速的轴线分布规律进行了实验研究，蒋建伟等口3 建立了破片杀伤威力参数的计算模型。本文中，基于V C 和M A T I A B，进一步提出破片对目标打击迹线的概念，探讨杀伤战斗部破片对典型目标打击迹线的计算模型，形成新的战斗部对特定目标的终点毁伤效应评估概念。2 破片对目标打击迹线概念破片对目标打击迹线与描述单个弹丸的运动轨迹及侵彻行为的射击线类似，破片对目标打击迹线是一条带有变化质量、速度和飞散方向的线段，该线段起点为战斗部起爆瞬间破片原始所处壳体微元的位置，终点为该破片飞行一段距离后能量衰减至临界有效杀伤特定目标的位置处，用以表述破片的运动轨迹和对欲打击目标的有效毁伤范围。因此，战斗部破片对目标打击迹线随打击目标的不同而不同。战斗部起爆后对目标的杀伤作用场可用多条破片对目标打击迹线描述，每条破片对目标打击迹线分别对应不同战斗部壳体位置处破片的运动轨迹和对特定目标打击能力变化情况。即使是同一战斗部，如所打击的目标对象发生变化，战斗部的破片对目标打击迹线分布情况也将不同。*收稿日期：2 0 0 9 0 3 1 7；修回日期：2 0 0 9 0 9 1 2基金项目：中国博士后科学基金项目(2 0 0 9 0 4 6 0 1 8 5)作者简介黄广炎(1 9 8 2)。男。博士，讲师。万方数据4 1 4爆炸与冲击第3 0 卷在已知战斗部装药结构、起爆姿态条件下，可以得到战斗部起爆驱动破片的离散初始威力特性参数，包括每个破片的质量、几何外形、三维空间坐标值、速度和飞散方向。以这些初始条件为基础，根据空气中外弹道方法计算所有破片的整个空气运行轨迹以及能量衰减情况，选择战斗部欲打击的目标并确定破片对该特定目标的毁伤准则，即可得到战斗部破片对目标打击迹线作用场的分布情况，从而精确评价破片对目标不同位置的打击效果。3 破片杀伤作用场威力参数计算模型破片对目标打击迹线的计算，需建立战斗部破片杀伤作用场各威力特性参数的计算模型。破片对目标的杀伤威力性能参数主要包括质量、初速、飞散方向、速度(能量)衰减等。取战斗部结构如图1 所示(预制破片壳体、轴对称装药、一端起爆)，分别讨论破片各威力参数。3 1 破片初速采用对战斗部破片初速格尼公式修正的计算方法H ，该方法适合各种装药结构、不同起爆方式下的图1 战斗部装药结构图F i g 1C h a r g es t r u c t u r eo fw a r h e a d战斗部破片初速计算，是一种广义破片初速计算方法。认为一般装药结构战斗部破片初速轴向分布函数主要与破片位置，和对应的质量比p 有关，即可o，=H(f d，8)V o破片速度值发生较大变化的位置在战斗部的两端，=H(f d，8)H z(孚，卢)u。式中：H，(f d)为起爆端修正系数，H z(L 一厂)d)为非起爆端修正系数。尼公式H 1(f d)=1 一A e-B：埘H z(生亍)=1 一c e _ D L 川。式中：A、B、C、D 为由实验数据确定的常数。假设沿弹轴各微元环形成的破片仍遵循格尼公式厄焉(2)通过X 射线测试结果修正格(3)(4)(5)于是一端起爆条件下，适用各种结构形状的杀伤战斗部的破片初速为础。，=-1-A 1e-B f l d(D-E l A z e-c(L-力似力 属、f l(f)i 而(6)3 2 破片飞散方向使用兰德一皮尔森对泰勒公式的修正计算方法，考虑了沿着金属表面紧挨着的两点，每一点都按照下面的方程做加速运动 o饥；u(1 一e 牮)(7)移、T、r 都随初始位置变化，得到破片飞散计算的修正公式为毋=a r c s i 吖爰一半一半1式中：艿，为破片飞散方向与壳体法线方向的夹角，白为破片初速轴向梯度，“为爆轰波扫过炸药壳体交界面的速度，“=D。c o s O，O s 为爆轰波传播方向与破片微元处壳体装药接触界面的夹角，r 为破片加速常数，r=O 4 d(f)(2 V o r)，d(，)为距起爆端厂处的炸药直径。万方数据第4 期黄广炎等：战斗部破片对目标打击迹线的计算方法4 1 5战斗部各破片飞散方向角为叩f 一1 f 一黾f式中：y，为炸药、破片接触界面的法线与装药轴线的夹角，妒，为破片飞散方向与战斗部轴线的夹角。3 3 子弹破片速度衰减破片在运动过程中因受空气阻力速度不断衰减，当速度减小到某一值时将不再满足杀伤作用的要求，丧失杀伤能力。假设破片飞行距离为s 时，破片速度衰减为哺3一r D 胛巩一l 3u 一v o r e l。，e 广一(9)式中：C o 为常数，对球形破片C-，一o 9 7，|D 为当地空气密度，m 为破片质量。3 4 杀伤作用场坐标系建立以及参数转换假设战斗部起爆作战几何空间如图2 所示，战斗部距地面高度为H，速度为，着角为口。O x y z 为地面坐标系，0。F 为战斗部轴线，0。为起爆点，P 为战斗部任一破片，破片P 相对战斗部的初速为7 3。，相对地面坐标系的破片初速为 U d。，飞散方向与战斗部轴夹角为9，。P 7 为在地面的着点。假如到地面时，破片初速仍满足对某目标的杀伤要求，则P P 7 为破片P 的有效射击迹线。为便于将各破片参数转换成地面坐标系下的数值，引入柱坐标系0-F R e 和体坐标系0。F 德，如图3 所示。图2 战斗部起爆几何空间F i g 2E x p l o s i v eg e o m e t r ys p a c eo fw a r h e a d破片P 的相对地面坐标系的初始坐标一s i n a0 1 f V p f lc。即oI 口唧I+o1 口辟J(a)C r l i n d r i c a lc o o r d i n a t ef广刁一丛L 拙1。(b)B o d yc o o r d i n a t efW1 I广刁?丛L 斗一C O S S l n o 一C O S C o s o S I m图3 战斗都坐标系F i g 3C o o r d i n a t e so fw a r h e a dS l n dC o s o 0(x e=f pC O S O 一R Pc o s s i n a_ y P 一一f P s i n a R P c o s s i n o+HIz P=R P s i n er 可艮=V o f s l n g f c o s a u o，c o s P，c o s s I n 盯十v d C O S 6 r_ u P y。V o r s i n g js i n a+V o f c o s c p f c o s c o s a+v ds i n a【口如一。，c o s 9，s i n A 私(1 0)(1 2)(1 3)龟l 奄善oc一|l：、I卜J度却如却速f始初的系标坐面地对相的P片破兽o砌W，|f 沁咖咿”仉善moCSJ J罾o砌善mo导一聊一口丹毋R 解r 仁榔万方数据爆炸与冲击第3 0 卷4 破片对目标打击迹线计算将破片杀伤作用场各威力特性参数的计算模型用V C 和M A T I。A B 联合编程，可针对不同装药结构和起爆方式的战斗部对目标的有效打击迹线作用场分布进行数值计算，计算步骤为输入战?l 部结构参数(炸药性质、破片材料、皱药绵构、起爆力式、端部情况等)l输入战斗部起爆繁什(士1 岛，姿态，着速)I依次，俑环计铬钾破片初盎f 威力笏数(形状、选择|1 标赝龌、舌u 速行分鼍、b 散J、初始啦杯J(人员，E 机、车辆等)循环汁箅再破”的终点毁伤符参数(速度、能垃寝减、举杯变化、毁伤能力判断)li f 算所仃破片的打击迹线(针对j 体H 标)4 1破片打击迹线作用场分布计算实例图4 杀伤战斗部结构示意图F i g 4S t r u c t u r eo ff r a g m e n tw a r h e a d计算图4 所示杀伤战斗部在6 种不同起爆姿态下对飞机目标打击迹线分布，计算中破片对飞机目标的杀伤标准使用8 0 0J c m 2 的比动能杀伤准则。杀伤战斗部壳体材料为3 5 C r M n S i A 4 2 C r，炸药装药为8 7 0 1，装药密度p=1 7 2g c m 3，战斗部长度为1 4 0m m，装药总质量为5 0 0g，起爆药柱为A 5 炸药，战斗部两端内置钨球破片，中间段为内刻槽壳体破片，一端起爆，战斗部两端端盖堵塞。沿战斗部轴向排布1 8 个典型破片(钨球破片和刻槽破片)，在战斗部表面标记母线，在该标记母线上，各典型破片质心位置呈直线排布口。图5 不同起爆姿态下战斗部对飞机目标打击迹线分布图F i g 5S h o t l i n e st op l a n eu n d e rd i f f e r e n te x p l o s i v ec o n d i t i o no fw a r h e a d万方数据第4 期黄广炎等：战斗部破片对目标打击迹线的计算方法4 1 7计算结果如图5 所示，图中每根细线分别代表相应战斗部壳体位置处破片在比动能满足对飞机的毁伤要求下的飞行弹道。计算结果比较清晰地表述了战斗部在不同起爆姿态下破片对飞机目标的杀伤作用场分布，直观地体现了战斗部所有破片的1 l 行轨迹和能量衰减情况。在作战空间，通过判断飞机各要害舱段所处化置区域破片打击迹线穿过的情况。就可以精确评估破片对飞机目标的毁伤效果，从而可对应调整战斗部最佳起爆姿态，达到对目标的最佳毁伤效果。4 2 与实验数据的对比杀伤战斗部静爆实验口 见图6。在距离杀伤战斗部心6m 处共布置1 0 块长1 5m、高21 T I、厚4m i D _ 的A 3 钢板(等效飞机目标)，形成1 3 0。左右的扇形拦截靶。战斗部非起爆端指向逆时针方向第2块靶下边缘处，战斗部表面标记母线位置对准靶板水平中心位置处。战斗部距离地面高度0 7 5I T I，以逆时针方向第2 块钢板下边缘为零经度线，相对战斗部中心每隔2。逆时针方向划经度线，在与战斗部标记母线同一水平高度的靶板位置画直线作为零纬度线，相对战斗部中心每隔1。划纬度线。破片对靶板打击迹线分布如图7 所示。图中，手为模拟目标靶板经度，为模拟目标靶板纬度。由图可知，计算的杀伤战斗部破片打击迹线与靶板位置交汇点坐标数据与实验中对靶板上的破片穿孔位置的统计数据基本吻合，较真实地体现了杀伤战斗部对日标的杀伤效果，验证了本文中所建立的破片杀伤作用场各威力参数计算模型和破片对目标打击迹线计算方法的正确性。图6 靶场实验布置F i g 6D i s p o s a lo fs e c t o rt a r g e te x p e r i m e n t图7 破片对靶板打击迹线分布F i g 7C a l c u l a t i n ga n de x p e r i m e n t a ls h o t l i n e st ot a r g e t5结论为了完善破片对目标打击迹线评估方法这种新的战斗部终点毁伤效应评估概念，在提出破片对目标打击迹线概念的基础上，建立系统的破片杀伤作用场各威力参数计算模型和破片对目标打击迹线计算方法，基于V C 和M A T I。A B 联合编程，对某战斗部在不同起爆姿态下对飞机目标打击迹线散布情况进行了计算，计算结果直观地反映了战斗部破片对目标的毁伤能力，为战斗部方案设计和战斗部终点毁伤效应分析提供了新的参考依据和技术途径。参考文献：1 Q I A NL i x i n，L I UT o n g，Z H A N GS h o u q i，e ta 1 F r a g m e n ts h o t l i n em o d e lf o ra i r d e f e n c ew a r h e a d J P r o p e l l a n t s，E x p l o s i v e s P y r o t e c h n i c s，2 0 0 0，2 5(2)：9 2 9 8 2 钱芷新。屈明。刘彤再入式大面积杀爆战斗部威力评定方法研究 J-I 弹道学报2 0 0 1，13(3)：4 2 4 6 Q I A NL i x i n，Q UM i n g，L I UT o n g L e t h a l i t ya s s e s s m e n ts t u d yo fr e e n t r yl a r g e-s c a l eb l a s t f r a g m e n t a t i o nw a r h e a d J J o u r n a lo fB a l l i s t i c s，2 0 0 1，1 3(3)：4 2 4 6 3 蒋建伟，卢永刚，钱立新射击迹线技术在战斗部破片场仿真中的应用 J 弹箭与制导学报，2 0 0 1。2 1(1)：2 9 3 4 J I A N G-J i a n-w e i，I UY o n g g a n g，Q I A NI i x i n A p p l i c a t i o no fs h o tl i n em o d e li ns i m u l a t i o no ff r a g m e n tw a r h e a d万方数据4 1 8爆炸与冲击第3 0 卷 J J o u r n a lo fP r o j e c t i l e s，R o c k e t s。M i s s i l e sa n dG u i d a n c e，2 0 0 1，2 1(1)：2 9 3 4 4 冯顺山，崔秉贵战斗部破片初速轴向分布规律的实验研究 J 兵工学报，1 9 8 7。1 1(4)：6 1 6 3 5 6 3 7 3F E N GS h u n s h a h，C U IB i n g g u i A ne x p e r i m e n t a li n v e s t i g a t i o nf o rt h ea x i a ld i s t r i b u t i o no fi n i t i a lv e l o c i t yo fs h e l l s J A c t aA r m a m e n t a r i i，1 9 8 7，1 1(4)：6 卜6 3 蒋建伟。冯顺山，何顺录一种用于评价和计算杀伤威力的数学模型 J 北京理工大学学报，1 9 9 2，1 2(2)：9 0 9 4 J I A N GJ i a n w e i，F E N GS h u n-s h a n，H ES h u n l u Am o d e lf o re v a l u a t i n ga n dc a l c u l a t i n gt h el e t h a lp o w e r J T r a n s a c t i o n so fB e i j i n gI n s t i t u t eo fT e c h n o l o g y，1 9 9 2，1 2(2)：9 0 9 4 隋树元，王树山终点效应学 M 北京：国防工业出版社，2 0 0 0 黄广炎，冯顺山基于战斗部微圆柱分析的破片飞散特性研究 C 第九届全国冲击动力学学术会议论文集2 0 0 9：1 6 5 1 6 9 Av i s u a lC(+)a n dM a t l a b-b a s e dc o m p u t a t i o n a lm e t h o df o rs h o t l i n e so fw a r h e a df r a g m e n t st oat a r g e t H U A N GG u a n g y a n l”，F E N GS h u n s h a n 2，L I UP e i q i n 9 1(1 S c h o o lo fA e r o n a t i cS c i e n c ea n dE n g i n e e r i n g，B e i j i n gU n i v e r s i t yo fA e r o n a u t i c&A s t r o n a u t i c，B e i ji n g1 0 0 0 8 3，C h i n a；2 S t a t eK e yL a b o r a t o r yo fE x p l o s i o nS c i e n c ea n dT e c h n o l o g y，B e i j i n gI n s t i t u t eo fT e c h n o l o g y，B e i j i n g1 0 0 0 8 1，C h i n a)A b s t r a c t：T h ee x i s t e n ts h o t l i n ec o n c e p tw a se x t e n d e db a s e do nt h ev i s u a lC(+)a n dM a t l a bc o d e s Ac o m p u t a t i o n a lm e t h o dw a sd e v e l o p e db ys e t t i n gu pt h ec o m p u t a t i o n a lm o d e lo ft h ed a m a g ef i e l dp a r a m e t e r s，w h i c hc o n s i s t i n go ft h ei n i t i a lv e l o c i t y，d i s p e r s i o nd i r e c t i o n，v e l o c i t ya t t e n u a t i o n，a n dS Oo n，f o rt h ep r e f o r m e do rh a l f p r e f o r m e df r a g m e n tw a r h e a da t t a c k i n gt h et a r g e t，a n dt h et r a n s f o r m a t i o no ft h ef r a g m e n tp o w e rp a r a m e t e r si nt h es p a c ec o o r d i n a t e A n dan e wc o m p u t a t i o n a lm e t h o dw a sa c h i e v e dt op r e d i c tt h es h o t l i n ef i e l dd i s t r i b u t i o no ft h ef r a g m e n tw a r h e a db yu s i n gt h ec o m b i n e dp r o g r a mw i t hv i s u a lC(+)a n dM A T L A B Ac o m p u t a t i o n a le x a m p l ew a sg i v e nb ya p p l y i n gt h ea c h i e v e dc o m p u t a t i o n a lm e t h o dt od e t e r m i n et h es h o t-l i n ef i e l dd i s t r i b u t i o no ft h ef r a g m e n tw a r h e a du n d e rd i f f e r e n te x p l o s i o nc o n d i t i o n sa t t a c k i n gt h ep l a n et a r g e t T h ec o m p u t a t i o n a lr e s u l t sd e s c r i b e dd i s t i n c t l yt h ed a m a g ef i e l dd i s t r i b u t i o no ft h ew a r h e a df r a g m e n t st ot h ep l a n ea n dr e f l e c t e dv i s u a l l yt h et r a j e c t o r i e sa n de n e r g ya t t e n u a t i o n so fa l lt h ew a r h e a df r a g m e n t s T h ea c h i e v e dc o m p u t a t i o n a lm e t h o dw a sv a l i d a t e db yt h ec o m p a r i s o nb e t w e e nt h ec o m p u t a t i o n a ld a t aw i t ht h ek n o w ne x p e r i m e n t a lr e s u i t s T h i sr e s e a r c hc a np r o v i d ean e wr e f e r e n c ef o rt h ew a r h e a dd e s i g na n dt e r m i n a Id a m a g ee f f e c ta n a l y s i s K e yw o r d s：m e c h a n i c so fe x p l o s i o n；s h o t l i n e；M A T I。A B；f r a g m e n t；w a r h e a d；V CR e c e i v e d17M a r c h2 0 0 9；R e v i s e d1 2S e p t e m b e r2 0 0 9C o r r e s p o n d i n ga u t h o r：H U A N GG u a n g y a h，h u a n g g y 0 0 6 7 b i t e d u c n(责任编辑丁峰)万方数据战斗部破片对目标打击迹线的计算方法战斗部破片对目标打击迹线的计算方法作者：黄广炎，冯顺山，刘沛清，HUANG Guang-yan，FENG Shun-shan，LIU Pei-qing作者单位：黄广炎,HUANG Guang-yan(北京航空航天大学航空科学与工程学院,北京,100083;北京理工大学爆炸科学与技术国家重点实验室,北京,100081)，冯顺山,FENG Shun-shan(北京理工大学爆炸科学与技术国家重点实验室,北京,100081)，刘沛清,LIU Pei-qing(北京航空航天大学航空科学与工程学院,北京,100083)刊名：爆炸与冲击英文刊名：EXPLOSION AND SHOCK WAVES年，卷(期)：2010,30(4)参考文献(7条)参考文献(7条)1.黄广炎;冯顺山 基于战斗部微圆柱分析的破片飞散特性研究 20092.隋树元;王树山 终点效应学 20003.蒋建伟;冯顺山;何顺录 一种用于评价和计算杀伤威力的数学模型期刊论文-北京理工大学学报 1992(02)4.冯顺山;崔秉贵 战斗部破片初速轴向分布规律的实验研究 1987(04)5.蒋建伟;卢永刚;钱立新 射击迹线技术在战斗部破片场仿真中的应用期刊论文-弹箭与制导学报 2001(01)6.钱立新;屈明;刘彤 再入式大面积杀爆战斗部威力评定方法研究期刊论文-弹道学报 2001(03)7.QIAN Li-xin;LIU Tong;ZHANG Shou-qi Fragment shot-line model for air-defence warhead 2000(02)本文链接：