杀爆类战斗部破片流密度计算方法技术

本发明专利技术公开了一种杀爆类战斗部破片流密度计算方法。该方法包括以下步骤：根据战斗部结构参数，计算战斗部静止爆炸时的破片飞散角范围，确定轴平面内破片的飞散角

Calculation method of flow density for fragmentation fragments of warhead

The invention discloses a method for calculating the flow density of fragmentation fragments of warhead fragments. The method includes the following steps: according to the structure parameters of warhead, the range of fragment divergence angle is calculated when the warhead explodes at rest, and the divergence angle of fragment in the axis plane is determined.

【技术实现步骤摘要】
杀爆类战斗部破片流密度计算方法
本专利技术武器系统毁伤效能评估
，特别是一种杀爆类战斗部破片流密度计算方法。
技术介绍
在进行杀爆类战斗部毁伤效能评估时，通常采用射击线方法，可到单个毁伤元对目标的毁伤效能，通过一定的方法将战斗部产生的所有毁伤元的毁伤效能合成为整个战斗部对目标的毁伤效能，计算简单、快捷。这对分析单枚破片对目标或部件的毁伤效能具有很高的精确性，但对于多枚破片同时毁伤目标时忽略了毁伤元的累积效应。比如，就单枚破片而言，无法毁伤目标时毁伤概率为0，利用射击线方法计算得出战斗部产生的同样破片无论击中目标的数量多少，对目标的毁伤概率均为0，这往往是不符合实际情况的。我国现有的毁伤评估系统多数为对毁伤元计算方法进行简化的射击线方法，很难准确描述多毁伤元的累积效应。
技术实现思路
本专利技术的目的在于提供一种计算方便、精确度高的杀爆类战斗部毁伤效能评估方法。实现本专利技术目的的技术解决方案为：一种杀爆类战斗部破片流密度计算方法，包括以下步骤：步骤1、根据战斗部结构参数，计算战斗部静止爆炸时的破片飞散角范围；步骤2、根据破片弹道信息和物理参数，计算战斗部爆炸后破片的空间分布及矢量信息；步骤3、根据弹目交会关系，计算战斗部破片流相对于目标的运动关系；根据破片的弹道信息，得到破片平均相对速度的方程，进而得到破片在相对运动中的坐标xm,ym,zm；步骤4、计算杀爆类战斗产生的破片流密度，战斗部静爆时所形成的破片流密度和破片场密度分别为进而求出相对运动破片流密度为为在经过战斗部纵轴平面内破片的平均飞散方向角，ηn为破片初始位置向量与导弹坐标系Z轴夹角，Nw为破片总数，D为破片落点与炸点距离，Sn为战斗部静爆时破片在半径为D的球面上分布场面积。本专利技术与现有技术相比，其显著优点：本专利技术采用破片流方法用于杀爆类战斗部毁伤效能评估，充分考虑了毁伤元的累积效应，能够提高杀爆类战斗部毁伤效能评估的准确性。附图说明图1为本专利技术杀爆类战斗部破片流密度计算方法的流程图。图2为破片的飞散范围。图3为破片矢量。图4为相对运动的破片流。图5为破片相对速度矢量。图6为相对运动飞散锥。图7为变形的破片流和目标的关系。图8为破片场的立体角。图9为战斗部破片分布场。图10为静态和动态破片流之间的关系。图11为破片场和目标舱段的交会示意图。具体实施方式下面结合附图和实施例对本专利技术作进一步说明。本专利技术一种杀爆类战斗部破片流密度计算方法，包括以下步骤：步骤1、根据战斗部结构参数，计算战斗部静止爆炸时的破片飞散角范围。如图2所示，静止状态下战斗部爆炸形成的破片流可用下列参数描述：Nw破片总数；qw单个破片的质量；αH破片的弹道系数；V0破片初速；在经过战斗部纵轴平面内破片的平均飞散方向角；在轴平面内破片的飞散角；破片流密度(单位面积上的破片数)。其中参数和可确定轴平面内破片的飞散角范围：步骤2、根据破片弹道信息和物理参数，计算战斗部爆炸后破片的空间分布及矢量信息。破片流的初始条件由向量确定，其在C·K·Om坐标中的位置用角ηn表示(图3)。对于结构确定的战斗部，破片参数Nw、qw、V0、基本上是确定，破片在范围内的分布具有一定的随机性，但不同类型的战斗部服从不同的分布规律(如均匀分布、正态分布)，在仿真时将此区间分为多个小的子区间，给出子区间的破片比例，来模拟破片不同的分布规律。步骤3、根据弹目交会关系，计算战斗部破片流相对于目标的运动关系。相对运动的破片流是指对目标的相对运动，也即处在目标上看到的破片流运动。如向量描述绝对运动中的破片弹道，则向量描述其相对弹道(如图4所示)。将方程两边除以破片飞行时间t，得到破片平均相对速度的方程：式中向量的值和其在C·K·Om中的位置(该位置用角和ηm表示，见图4)，很容易通过其在该坐标系轴上的投影求得：当攻角αm＝0时考虑到将上述关系代入上式得：当攻角αm＝0时如果方程的左部和右部都乘以破片的飞行时间t，则得到破片在相对运动中的坐标xm,ym,zm。当攻角αm＝0时xm＝xmy+Vttcosχ＝xn+Vmεt+Vttcosχym＝ymy＝ynzm＝zmy-Vttsinχ＝zn-Vttsinχ和相对运动中的飞行距离D＝Vt。xmy，ymy，zmy为在绝对运动中破片的坐标，而xn，yn，zn为在静态下，爆炸后经过时间t的坐标。这样，如果在静止中起爆战斗部时，破片具有某些固定值V0、ηn，则在给定的Vm、Vt、χ时，将由半径向量的末端划出其相对弹道，该向量的位置在每一瞬间由角和ηm表示。如果大气阻力很小，则系数ε≈1，从而，破片的相对速度V的量值和方向均为常值，这就意味着破片的相对弹道是直线和向量V方向相重合。当大气阻力大时，破片的相对弹道是非直线的。随着破片的绝对飞行距离Dy的增大，或者等价飞行时间t的增大，值ε减小。从而，随着破片飞行时间的增大其平均相对速度V将减小，并且向量相对目标方向的角和ηp将发生变化。在分析相对运动的破片流问题时，用导弹的相对速度向量Vr描述比较方便，即在坐标系OXrYrZr中，而不在导弹坐标系中。半径向量以及破片平均相对速度向量在C·K·Or中的位置将用角ηr表示(图5)：根据弹体坐标系和弹联相对坐标系的转换关系可得：显然为了求解方便，在研究破片流时，经常在xrzr平面内画出破片流的切面，在该平面内包含有速度向量三角形这里不考虑大气阻力。当在迎头和尾追航向上导弹攻击目标时，当sinχ＝0，导弹相对速度向量和导弹纵轴方向一致(图5)。这时，在相对运动中的破片流相对导弹纵轴和向量仍是对称的，并限制在两圆锥内，圆锥的母线倾角和如图6所示，假设导弹在C.K.O中的xoz平面内通过。如果相对目标战斗部在A点起爆，则破片流脱离目标。为了使破片落在目标上，必须沿导弹的相对弹道略晚起爆战斗部，例如在B点，。当sinχ≠0时，导弹以交叉航向攻击目标，破片流相对导弹纵轴是不对称的。在和方向不重合的向量作用下产生“变形”。在图6上列出了在xryr平面内，破片流的切面图。可以看出该切面的严重不对称性。破片流的变形将使导弹接近目标的所有可能情况下，引战配合设计更加困难。破片流的变形程度将取决于V0和Vr之比；当V0＞＞Vr时，破片流相对导弹纵轴在任意和间的夹角εp下非对称性都不大。在图6上，导弹在和目标固连的C.K.O的xoz平面内从目标的左边通过，当在A点起爆战斗部时，破片流将脱离目标，为了保证破片流覆盖目标，应略晚起爆，例如在B点。步骤4、计算杀爆类战斗部产生的破片流密度。定义飞散在单位立体角内的破片相对数量称为破片流密度，其值为：式中，为立体角内飞散的破片数；Nw为破片总数；η表示立体角相对炸点的方向角。破片流密度的概念适用于任何一种破片流(静止、绝对和相对运动中起爆战斗部时)，因此，在上式中略去参数和η的下标。下面对具体的破片流应对参数和η给出相应的下标。定义落在离炸点距离D垂直于破片飞行方向的单位面积上的破片数为破片场密度下面分析和之间的关系。设破片场密度已知，则落在垂直于半径向量的单位面积ΔS上的破片平均数为(图8)：对应面积ΔS的立体角为(见图8)：可得：可写为：可见，破片场的密度反比于离炸点距离的平方。战斗部静爆时，将破片流密度和破片场密度记为和常在研制和设计战斗部时，力求使静止下起爆战斗部本文档来自技高网...

【技术保护点】
1.一种杀爆类战斗部破片流密度计算方法，其特征在于包括以下步骤：步骤1、根据战斗部结构参数，计算战斗部静止爆炸时的破片飞散角范围；步骤2、根据破片弹道信息和物理参数，计算战斗部爆炸后破片的空间分布及矢量信息；步骤3、根据弹目交会关系，计算战斗部破片流相对于目标的运动关系；根据破片的弹道信息，得到破片平均相对速度的方程，进而得到破片在相对运动中的坐标xm,ym,zm；步骤4、计算杀爆类战斗产生的破片流密度，战斗部静爆时所形成的破片流密度和破片场密度分别为

【技术特征摘要】
1.一种杀爆类战斗部破片流密度计算方法，其特征在于包括以下步骤：步骤1、根据战斗部结构参数，计算战斗部静止爆炸时的破片飞散角范围；步骤2、根据破片弹道信息和物理参数，计算战斗部爆炸后破片的空间分布及矢量信息；步骤3、根据弹目交会关系，计算战斗部破片流相对于目标的运动关系；根据破片的弹道信息，得到破片平均相对速度的方程，进而得到破片在相对运动中的坐标xm,ym,zm；步骤4、计算杀爆类战斗产生的破片流密度，战斗部静爆时所形成的破片流密度和破片场密度分别为进而求出相对运动破片流密度为为在经过战斗部纵轴平面内破片的平均飞散方向角，ηn为破片初始位置向量与导弹坐标系Z轴夹角，Nw为破片总数，D为破片落点与炸点距离，Sn为战斗部静爆时破片在半径为D的球面上分布场面积。2.根据权利要求1所述的杀爆类战斗部破片流密度计算方法，其特征在于：所述步骤1中，根据轴平面内破片的飞散角公式计算战斗部静止爆炸时的破片飞散角...

【专利技术属性】
技术研发人员：周兰伟，宫小泽，李向东，姚志军，司凯，
申请(专利权)人：南京理工大学，中国人民解放军六三八五零部队，
类型：发明
国别省市：江苏,32