一种自然破片战斗部空爆冲击波和高速破片先后作用的临界爆距计算方法
【专利摘要】本发明涉及一种自然破片战斗部空中爆炸下冲击波与高速破片先后作用的临界爆距计算方法，它有七大步骤：步骤1，计算破片的初始速度；步骤2，计算等效裸装药大小；步骤3，计算战斗部的极限膨胀半径；步骤4，推导破片在空气中飞行距离与时间的关系；步骤5，采用步骤2的计算结果得到空气中冲击波波阵面的传播距离与时间的关系；步骤6，在时间相同，冲击波波阵面传播距离与破片飞行距离相等的条件下，计算得到冲击波和高速破片相遇时运动的距离；步骤7，将步骤3和步骤6的结果求和，得到自然破片战斗部空中爆炸下冲击波和高速破片先后作用的临界爆距值。本发明能较准确地计算出自然破片战斗部空爆冲击波和高速破片先后作用的临界爆距值，在毁伤及防护等领域均具有较高的实用价值和广阔的应用前景。
【专利说明】
一种自然破片战斗部空爆冲击波和高速破片先后作用的临界 爆距计算方法
技术领域
[0001] 本发明涉及一种自然破片战斗部空中爆炸下毁伤元素形成特性及传播(运动）规 律的计算方法，尤其涉及自然破片战斗部空中爆炸下冲击波和高速破片先后作用临界爆距 的动力学计算方法，可应用于毁伤和防护等领域，以合理评估战斗部毁伤效果或为防护结 构的设计提供参考依据等。
【背景技术】
[0002] 破片战斗部空中爆炸下会同时产生冲击波和高速破片两种毁伤载荷。另外，由于 炸药本身化学反应，还会形成一定的爆轰产物。相比较而言，接触爆炸对结构造成的毁伤要 更为严重，但是在实战中非接触爆炸出现的概率要远大于接触爆炸。因此，对于空中非接触 爆炸的情形，冲击波和高速破片则是两种主要的毁伤元素。
[0003] 破片战斗部空中爆炸下初始冲击波的传播速度和高速破片初始飞行速度相差很 大，又由于两者在空气中的速度衰减效应差异明显，因而导致冲击波和破片在传播(运动） 过程中存在相遇的问题，此相遇位置即为冲击波和高速破片先后作用的临界爆距值。爆炸 初期，形成的初始冲击波传播速度很高，要远高于形成的破片的初始速度;但冲击波在空气 中的速度衰减非常快，而破片则慢得多。因而，对于战斗部空爆而言，在临界爆距之前，冲击 波运动在破片之前;而超过临界爆距值后，破片始终运动在冲击波之前。
[0004] 当目标离破片战斗部较近时，其爆炸产生的冲击波会先作用在目标上，使目标产 生变形或破损，随后破片对目标形成穿甲破坏；当目标离破片战斗部较远时，爆炸产生的破 片会先作用于目标，形成穿甲破坏，随后冲击波对目标形成进一步的毁伤;而当目标处于临 界爆距位置处时，冲击波和破片会同时作用于目标，两者的能量会形成联合增强效应，可能 会使目标产生更为严重的破坏。对于战斗部设计者来说，如果能准确地掌握该战斗部的临 界爆距值，则对进一步提高战斗部毁伤威力的设计具有指导意义;而在防护领域，如果明确 了防御目标战斗部的临界爆距值，则对于尽可能减小战斗部爆炸毁伤效果的防护结构的设 计具有参考价值。因此，如何能够较为准确且方便快捷地得到战斗部爆炸下冲击波与破片 先后作用的临界爆距值，成为毁伤及防护等领域深入研究的前提和关键。
[0005] 然而，由于问题的复杂性，破片战斗部空中爆炸下的冲击波和高速破片先后作用 的临界爆距计算问题是毁伤和防护等领域计算的难点，一直以来没有一种明确的、可用于 工程实际的分析计算方法。目前对于破片战斗部空爆下冲击波和高速破片先后作用的计算 主要采用数值仿真的方法。但由于涉及壳体材料、炸药和空气等多种介质的流固耦合问题， 使得有限元数值仿真计算需耗费大量的计算资源和时间，且由于有限元仿真计算过程受网 格大小的影响较大，其精确度和可靠性需要验证;而无网格数值仿真计算同样存在耗费大 量计算资源和时间等问题。因此，急需提出一种相应的计算方法。
【发明内容】
[0006] 本发明所要解决的技术问题是提出一种自然破片战斗部空爆冲击波和高速破片 先后作用的临界爆距计算方法，能较好地符合自然破片战斗部空中爆炸下冲击波和高速破 片这两种毁伤元素的形成特点和传播(运动)规律。
[0007] 为了解决上述问题，本发明提出了一种自然破片战斗部空爆冲击波和高速破片先 后作用的临界爆距计算方法，包括以下步骤：
[0008] 步骤1，根据自然破片战斗部爆炸形成高速破片的理论，计算出高速破片的初始速 度；
[0009] 步骤2,根据战斗部装药的总内能等于爆轰产物的内能和动能及高速破片动能等 能量之和的总能量守恒关系，计算得到战斗部的等效裸装药大小；
[0010] 步骤3:忽略爆轰产物的绕流效应引起的能量损失，结合能量守恒原理，计算得到 战斗部的极限膨胀半径；
[0011] 步骤4,根据冲量守恒原理和牛顿第二定律，由破片在空气中飞行所受到的阻力公 式，推导得到破片在空气中飞行距离与所需时间的关系式；
[0012] 步骤5,根据爆炸动力学理论，采用步骤2的计算结果得到空气中冲击波波阵面的 传播距离与时间的关系式；
[0013] 步骤6,在时间相同，冲击波波阵面传播距离与破片飞行距离相等的条件下，计算 得到自然破片战斗部空中爆炸下冲击波和高速破片在空中相遇的位置；
[0014] 步骤7,在步骤3和步骤6的基础上，得到自然破片战斗部空中爆炸下冲击波和高速 破片先后作用的临界爆距值。
[0015] 进一步地，该步骤1中所述的高速破片初始速度vo计算公式为：
[0017] 式中：co为战斗部的装药量，ms为战斗部壳体质量，E为战斗部装药的特征能量。
[0018] 进一步地，若战斗部的形状为球形，则步骤1中所述的高速破片的初始速度VQ计算 公式为：
[0020]进一步地，若已知战斗部装药的爆速D，则战斗部装药的特征能量E可根据下式计 算得到：
[0022]进一步地，该步骤2中所述的战斗部等效裸装药大小为：
[0024]式中:Efk为破片的总动能，Et〇tai为战斗部装药总内能。
[0025]进一步地，破片的总动能Efk为：
[0026] = 0.5msv〇
[0027] 进一步地，战斗部装药的总内能Etotal为：
[0028] Etotai = CO ? Qv
[0029] 式中：Qv为炸药的爆热。
[0030] 进一步地，该步骤3所述的战斗部极限膨胀半径Ro需通过求解下式得到：
[0032]式中：pf和hf分别为战斗部壳体厚度和材料密度，为装药密度，ro为战斗部壳体 的初始半径。
[0033]进一步地，该步骤4中所述的破片在空气中飞行距离Rf与所需时间tf的关系为：
[0035]式中:Kf为阻力系数，mf为破片的设计质量。
[0036] 进一步地，阻力系数Kf为：
[0037] K{ = -cfpaA
[0038] 式中：cf为与破片形状相关的迎风阻力系数，pa为空气的质量密度，j为破片的平 均迎风面积。
[0039]进一步地，与破片形状相关的迎风阻力系数《为： '0.97， 球形
[0040] cf =<1.24, 立方体或规则矩形 1.50， 不规则矩形
[0041] 进一步地，破片的平均迎风面积]由下式计算：
[0042] A = a- mf
[0043]式中:a为破片的形状系数。
[0044] 进一步地，破片形状系数a的取值为(单位:l(T3m2/kg2/ 3): '3.07, 球形 3.09, 立方体
[0045] a = < 3.55, 规则矩形 4.42, 不规则矩形
[0046] 进一步地，破片的设计质量mf通过以下公式计算：
[0047] mf = 2BX ? hshdl t/inner(l + /?sheli /^inner)
[0048]式中：Bx为炸药的爆炸常数(单位kgV2/mV6)，h sheii和dinner分别为战斗部壳体的平 均厚度和平均内径。
[0049]进一步地，该步骤5中所述采用步骤2的计算结果得到空气中冲击波波阵面的传播 距呙Rs与时间ts的关系式为：
[0050] Rs = (〇.22 Xl04 -CO°eq16/3 -tj1116
[0051]式中即为步骤2中战斗部等效裸装药大小的计算结果。
[0052]进一步地，该步骤6中所述的，在时间ts和tf相同，冲击波波阵面传播距离R s与破片 飞行距离Rf相等的条件下，通过求解以下方程得到自然破片战斗部空中爆炸下破片追赶上 冲击波波阵面的时间tm:
[0054]进一步地，可求得自然破片战斗部空中爆炸下高速破片追赶上冲击波波阵面的距 离1，可采用下面两式中的任意一式计算：
[0057] 进一步地，该步骤7中所述自然破片战斗部空中爆炸下冲击波和高速破片先后作 用的临界爆距值Rct为：
[0058] Rcr = R〇+Rm
[0059] 式中：Ro为战斗部的极限膨胀半径，心为冲击波和高速破片空中相遇时运动(传播） 的距离。
[0060] 与现有技术相比，本发明具有以下优点：
[0061] 本发明提出的一种自然破片战斗部空爆冲击波和高速破片先后作用的临界爆距 计算方法通过动力学计算能够简便且较为准确地得到自然破片战斗部空中爆炸下冲击波 和高速破片先后作用的临界爆距值大小，从而能够为相关武器战斗部或防护结构等设计及 优化提供方便快捷且可靠的参考依据。
【附图说明】
[0062] 图1是本发明的一种自然破片战斗部空中爆炸下冲击波和高速破片先后作用临界 爆距的动力学计算方法的流程图。
[0063] 图2是本发明实施例中冲击波波阵面传播速度和高速破片飞行速度的衰减曲线。 [0064]图3是本发明实施例中冲击波波阵面传播距离和高速破片飞行距离的变化曲线。
[0065] 图4是本发明实施例中冲击波波阵面传播速度及其传播距离随时间的变化曲线。
[0066] 图5是本发明实施例中高速破片飞行速度及其飞行距离随时间的变化曲线。
【具体实施方式】
[0067] 为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实 施方式对本发明作进一步详细的说明。
[0068] 本发明目的是要提出可应用于工程实际的一种自然破片战斗部空中爆炸下冲击 波与高速破片先后作用的临界爆距计算方法。
[0069] 请参照图1，为解决上述问题，本发明基于自然破片战斗部空中爆炸所产生的冲击 波和高速破片的形成特点及传播(运动)规律，提出了一种自然破片战斗部空中爆炸下冲击 波与破片先后作用的临界爆距动力学计算方法，包括以下步骤：
[0070] 步骤1，根据自然破片战斗部爆炸破片形成理论，计算爆炸形成的高速破片的初始 速度；
[0071] 步骤2,根据战斗部装药的总内能等于爆轰产物的内能和动能及高速破片动能等 能量之和的总能量守恒关系，计算得到战斗部的等效裸装药大小；
[0072] 步骤3:忽略爆轰产物的绕流效应引起的能量损失，结合能量守恒原理，计算得到 战斗部的极限膨胀半径；
[0073] 步骤4,根据冲量守恒原理和牛顿第二定律，由破片在空气中飞行所受到的阻力公 式，推导得到破片在空气中飞行距离与所需时间的关系式；
[0074] 步骤5,根据爆炸动力学理论，采用步骤2的计算结果得到空气中冲击波波阵面的 传播距离与时间的关系式；
[0075]步骤6,在时间相同，冲击波波阵面传播距离与破片飞行距离相等的条件下，计算 得到自然破片战斗部空中爆炸下冲击波和高速破片在空气中的相遇距离；
[0076] 步骤7,在步骤3和步骤6的基础上，得到自然破片战斗部空中爆炸下冲击波和高速 破片先后作用的临界爆距值。
[0077] 下面结合实例对本发明的一种自然破片战斗部空中爆炸下冲击波和高速破片先 后作用的临界爆距计算方法进行详细说明。实施例为国外某一自然破片战斗部空中爆炸的 情形。该破片战斗部的炸药当量为34.93kg，爆炸所产生的破片设计质量及其飞行参数为：
[0078] 表1爆炸所产生的破片设计质量及其飞行参数
[0080] 下面是具体步骤：
[0081] 步骤1，根据自然破片战斗部爆炸破片形成理论，计算出爆炸形成的高速破片的初 始速度；
[0082] 根据前述步骤1中的计算公式可得出破片的初始速度为887m/s。
[0083] 步骤2,根据战斗部装药的总内能等于爆轰产物的内能和动能及高速破片动能等 能量之和的总能量守恒关系，计算得到战斗部的等效裸装药大小；
[0084] 将战斗部装药的总内能减去战斗部产生的破片的总动能，然后再将剩下的总能量 转化为裸装药质量，根据前述步骤2中相应的转化公式得到战斗部的等效裸装药为 23.75kg〇
[0085] 步骤3:忽略爆轰产物的绕流效应引起的能量损失，结合能量守恒原理，计算得到 战斗部的极限膨胀半径；
[0086] 根据前述步骤3中的方程，代入各参数求解得到战斗部的极限膨胀半径值约为 0.28m〇
[0087]步骤4,根据冲量守恒原理和牛顿第二定律，由破片在空气中飞行所受到的阻力公 式，推导得到破片在空气中飞行距离与所需时间的关系式；
[0088]根据前述步骤3,计算得到破片在空气中飞行时的阻力系数Kf = -0.3 XI. 25 X 19 ? 72 X 10-4 = -7 ? 4 X 10-4(kg/m) 〇
[0089] 进一步地，根据前述步骤4中的关系式，得到破片在空气中飞行距离与飞行时间的 关系式为=887.-393.4.4+232.6.? f3 (m)。
[0090] 步骤5,根据爆炸动力学理论，采用步骤2的计算结果得到空气中冲击波波阵面的 传播距离与时间的关系式；
[0091] 根据前述步骤4的关系式，得到空气中冲击波波阵面传播距离与所需时间的关系 式为：i?s= 125.1 .，s° 568 (m)〇
[0092]步骤6,在时间相同，冲击波波阵面传播距离与破片飞行距离相等的条件下，计算 得到破片战斗部空中爆炸下冲击波与破片先后作用的临界爆距值大小。
[0093]根据前述步骤6,通过联立步骤4和步骤5中飞行(传播)距离与时间的关系式，在时 间相同距离相等的条件下，求得破片追赶上冲击波波阵面的时间为10.8ms。
[0094] 进一步地，根据步骤4或步骤5中的关系式，求得该战斗部空中爆炸下冲击波与高 速破片相遇时运动(传播)的距离大小为9.56m。
[0095] 步骤7,在步骤3和步骤6的基础上，得到自然破片战斗部空中爆炸下冲击波和高速 破片先后作用的临界爆距值。
[0096] 进一步地，根据步骤7中关系式，得到该自然破片战斗部空中爆炸下冲击波和高速 破片先后作用的临界爆距值为9.84m。文献中实验值为10.2m，计算结果与文献实测值间的 相对误差约为-3.5%。考虑到实验测试的精度误差，计算结果与文献中的实验值吻合较好， 能够满足工程应用的需要。
【主权项】
1. 一种自然破片战斗部空爆冲击波和高速破片先后作用的临界爆距计算方法，其特征 在于该方法由以下步骤实现:步骤1，计算破片的初始速度;步骤2，计算等效裸装药大小;步 骤3,计算战斗部的极限膨胀半径;步骤4,推导破片在空气中飞行距离与时间的关系；步骤 5,采用步骤2的计算结果得到空气中冲击波波阵面的传播距离与时间的关系；步骤6,在时 间相同，冲击波波阵面传播距离与破片飞行距离相等的条件下，计算得到冲击波和高速破 片相遇时运动的距离;步骤7,将步骤3和步骤6的结果求和，得到自然破片战斗部空中爆炸 下冲击波和高速破片先后作用的临界爆距值。2. 如权利要求1所述的计算方法，其特征是步骤2、3的顺序可以任意更换。3. 如权利要求1所述的计算方法，其特征是步骤4、5的顺序可以任意更换。4. 如权利要求1所述的计算方法，其特征是步骤6中所述的先后作用临界爆距值大小Rct 等于办与1之和，Ro表示战斗部的极限膨胀半径，Rm表示冲击波和高速破片相遇时所运动的 距离。
【文档编号】G06F17/50GK105912744SQ201610145789
【公开日】2016年8月31日
【申请日】2016年3月15日
【发明人】陈长海, 朱锡, 侯海量, 胡年明
【申请人】中国人民解放军海军工程大学