CN108563893B - 杀爆类战斗部破片流密度计算方法 - Google Patents
杀爆类战斗部破片流密度计算方法 Download PDFInfo
- Publication number
- CN108563893B CN108563893B CN201810373342.6A CN201810373342A CN108563893B CN 108563893 B CN108563893 B CN 108563893B CN 201810373342 A CN201810373342 A CN 201810373342A CN 108563893 B CN108563893 B CN 108563893B
- Authority
- CN
- China
- Prior art keywords
- fragment
- warhead
- fragments
- calculating
- explosion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
Description
技术领域
本发明武器系统毁伤效能评估技术领域,特别是一种杀爆类战斗部破片流密 度计算方法。
背景技术
在进行杀爆类战斗部毁伤效能评估时,通常采用射击线方法,可到单个毁伤 元对目标的毁伤效能,通过一定的方法将战斗部产生的所有毁伤元的毁伤效能合 成为整个战斗部对目标的毁伤效能,计算简单、快捷。这对分析单枚破片对目标 或部件的毁伤效能具有很高的精确性,但对于多枚破片同时毁伤目标时忽略了毁 伤元的累积效应。比如,就单枚破片而言,无法毁伤目标时毁伤概率为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、根据破片弹道信息和物理参数,计算战斗部爆炸后破片的空间分布 及矢量信息。
对于结构确定的战斗部,破片参数Nw、qw、V0、基本上是确定,破片 在范围内的分布具有一定的随机性,但不同类型的战斗部服从不同的分布规 律(如均匀分布、正态分布),在仿真时将此区间分为多个小的子区间,给出子区 间的破片比例,来模拟破片不同的分布规律。
步骤3、根据弹目交会关系,计算战斗部破片流相对于目标的运动关系。
相对运动的破片流是指对目标的相对运动,也即处在目标上看到的破片流运 动。
当攻角αm=0时
当攻角αm=0时
如果方程的左部和右部都乘以破片的飞行时间t,则得到破片在相对运动中 的坐标xm,ym,zm。
当攻角αm=0时
xm=xmy+Vttcosχ=xn+Vmεt+Vttcosχ
ym=ymy=yn
zm=zmy-Vttsinχ=zn-Vttsinχ
和相对运动中的飞行距离D=Vt。
xmy,ymy,zmy为在绝对运动中破片的坐标,而xn,yn,zn为在静态下,爆 炸后经过时间t的坐标。
如果大气阻力很小,则系数ε≈1,从而,破片的相对速度V的量值和方向 均为常值,这就意味着破片的相对弹道是直线和向量V方向相重合。
当大气阻力大时,破片的相对弹道是非直线的。随着破片的绝对飞行距离Dy的增大,或者等价飞行时间t的增大,值ε减小。从而,随着破片飞行时间的增 大其平均相对速度V将减小,并且向量相对目标方向的角和ηp将发生变化。
在分析相对运动的破片流问题时,用导弹的相对速度向量Vr描述比较方便, 即在坐标系OXrYrZr中,而不在导弹坐标系中。
根据弹体坐标系和弹联相对坐标系的转换关系可得:
显然
如图6所示,假设导弹在C.K.O中的xoz平面内通过。如果相对目标战斗部 在A点起爆,则破片流脱离目标。为了使破片落在目标上,必须沿导弹的相对 弹道略晚起爆战斗部,例如在B点,。
当sinχ≠0时,导弹以交叉航向攻击目标,破片流相对导弹纵轴是不对称的。 在和方向不重合的向量作用下产生“变形”。在图6上列出了在xryr平面内, 破片流的切面图。可以看出该切面的严重不对称性。破片流的变形将使导弹接近 目标的所有可能情况下,引战配合设计更加困难。破片流的变形程度将取决于V0和Vr之比;当V0>>Vr时,破片流相对导弹纵轴在任意和间的夹角εp下非 对称性都不大。
在图6上,导弹在和目标固连的C.K.O的xoz平面内从目标的左边通过,当 在A点起爆战斗部时,破片流将脱离目标,为了保证破片流覆盖目标,应略晚 起爆,例如在B点。
步骤4、计算杀爆类战斗部产生的破片流密度。
对应面积ΔS的立体角为(见图8):
可得:
可写为:
可见,破片场的密度反比于离炸点距离的平方。
式中,Sn为战斗部静爆时破片在半径为D的球面上分布场面积。
对于径向破片战斗部,杀伤元素分布在球面上(图9所示),因此
得到:
显然
取垂直于半径向量V0εt的单位面积ΔS(见图10)。落在该面积上的破片平均 数:
由于ΔS面积很小,可以认为所有落在其上的破片飞行时间t是相等的。因 此,在相对运动中所有落在上面的破片还沿向量移动相等的距离,这 等价于面积ΔS在时间t内平行地从An点移动到A0点。这样在静态中飞散在立体 角Δηn中的破片在相对运动中将限制在立体角Δηm中,该立体角支 撑在同一个面积ΔS上,但该ΔS位于离炸点距离D=Vt的A0点上。据此
可得:
所以有
得到
cosv=lVlV0+mVmV0+nVnV0
最终得到
当sinχ=sinq=0(迎头或尾追攻击目标),ηn=ηm从而
如果不计空气阻力(ε=1),则:
这时
这样,如果已知导弹接近目标的条件Vm、Vt、αm、χ(或q),静止下起爆战 斗部时的破片流密度以及向量可以求出破片的相对速 度V和表示其在相对运动中飞行方向的角ηm;也可以求出在向量方向上的 破片流密度。
实施例
本实施例应用于杀爆类战斗部破片流密度的计算,具体如下:
Claims (3)
1.一种杀爆类战斗部破片流密度计算方法,其特征在于包括以下步骤:
步骤1、根据战斗部结构参数,计算战斗部静止爆炸时的破片飞散角范围;
步骤2、根据破片弹道信息和物理参数,计算战斗部爆炸后破片的空间分布及矢量信息;
步骤3、根据弹目交会关系,计算战斗部破片流相对于目标的运动关系;根据破片的弹道信息,得到破片平均相对速度的方程,进而得到破片在相对运动中的坐标xm,ym,zm;具体如下:
其中,χ为截获角;Vm为导弹运动速度;Vt为目标运动速度;αm为导弹攻角;xn、yn、zn为破片在战斗静态下爆炸后经过时间t的坐标;ε为距离D处的破片与爆炸初始时速度之比;
相对运动时流密度根据以下计算公式得到:
3.根据权利要求1所述的杀爆类战斗部破片流密度计算方法,其特征在于:所述步骤2中,根据战斗部类型选择不同的分布规律来模拟破片不同的分布规律。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810373342.6A CN108563893B (zh) | 2018-04-24 | 2018-04-24 | 杀爆类战斗部破片流密度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810373342.6A CN108563893B (zh) | 2018-04-24 | 2018-04-24 | 杀爆类战斗部破片流密度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108563893A CN108563893A (zh) | 2018-09-21 |
CN108563893B true CN108563893B (zh) | 2022-08-09 |
Family
ID=63536571
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810373342.6A Active CN108563893B (zh) | 2018-04-24 | 2018-04-24 | 杀爆类战斗部破片流密度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108563893B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111125623B (zh) * | 2019-11-28 | 2023-09-19 | 北京电子工程总体研究所 | 一种超低空作战条件下的战斗部破片飞散轨迹计算的方法 |
CN110991030B (zh) * | 2019-11-28 | 2023-10-13 | 北京电子工程总体研究所 | 一种通用复杂引战系统仿真的方法及系统 |
CN112035780B (zh) * | 2020-09-04 | 2022-05-31 | 清华大学 | 一种导弹末制导阶段杀伤效果计算方法 |
CN112214896B (zh) * | 2020-10-14 | 2022-09-16 | 贵州理工学院 | 一种活性破片战斗部毁伤威力场的获取方法 |
CN112414232B (zh) * | 2020-11-05 | 2023-01-03 | 江西洪都航空工业股份有限公司 | 一种杀爆战斗部飞散角的生成方法和杀爆战斗部 |
CN112733366B (zh) * | 2021-01-12 | 2022-09-16 | 中国人民解放军陆军军事交通学院军事交通运输研究所 | 一种实时高爆弹破片毁伤仿真方法 |
CN113361085A (zh) * | 2021-05-28 | 2021-09-07 | 西安工业大学 | 一种战斗部破片密度计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6421116B1 (en) * | 1999-08-23 | 2002-07-16 | Bodenseewerk Geratetechnik Gmbh | Method for determining the relative movement between missile and target |
CN104850747A (zh) * | 2015-05-26 | 2015-08-19 | 中国人民解放军63867部队 | 一种弹药威力杀伤体积评价方法 |
CN107798208A (zh) * | 2016-08-28 | 2018-03-13 | 南京理工大学 | 对空目标导弹破片飞散最大毁伤算法 |
-
2018
- 2018-04-24 CN CN201810373342.6A patent/CN108563893B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6421116B1 (en) * | 1999-08-23 | 2002-07-16 | Bodenseewerk Geratetechnik Gmbh | Method for determining the relative movement between missile and target |
CN104850747A (zh) * | 2015-05-26 | 2015-08-19 | 中国人民解放军63867部队 | 一种弹药威力杀伤体积评价方法 |
CN107798208A (zh) * | 2016-08-28 | 2018-03-13 | 南京理工大学 | 对空目标导弹破片飞散最大毁伤算法 |
Non-Patent Citations (4)
Title |
---|
偏心起爆式定向战斗部破片分布及毁伤效能研究;吕胜涛;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20120415;第C032-17页 * |
基于射击线技术的杀爆战斗部杀伤面积计算;应国淼等;《兵器装备工程学报》;20161125(第11期);全文 * |
瞄准式战斗部杀伤装置结构设计研究;石志彬 等;《兵工学报》;20130331;第34卷(第3期);全文 * |
破片式战斗部对飞机类目标毁伤评估方法研究;司凯等;《弹道学报》;20171215(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108563893A (zh) | 2018-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108563893B (zh) | 杀爆类战斗部破片流密度计算方法 | |
CN107798208B (zh) | 对空目标导弹破片飞散最大毁伤算法 | |
Liu et al. | Formation of explosively formed penetrator with fins and its flight characteristics | |
CN108694281B (zh) | 一种基于多维角度破片散布侵彻的目标毁伤概率计算方法 | |
CN109063331B (zh) | 一种小脱靶量破片战斗部对再入弹头毁伤效应的评估方法 | |
CN112035780B (zh) | 一种导弹末制导阶段杀伤效果计算方法 | |
Ji et al. | Dynamic behaviors of multi-layered steel targets with air gaps subjected to the impact of EFP simulants | |
Massey et al. | Testing the maneuvering performance of a mach 4 projectile | |
Massey et al. | Combining experimental data, computational fluid dynamics, and six-degree of freedom simulation to develop a guidance actuator for a supersonic projectile | |
RU192765U1 (ru) | Разделяющаяся боевая часть с планирующими боевыми блоками | |
Chen et al. | Motion characteristics after ricochet: an experimental investigation | |
Li et al. | Modeling and calculation method of target damage based on multi-attitude flying projectile in space intersection | |
KR20170106743A (ko) | 잠수함발사탄도미사일의 비행특성을 시뮬레이션하는 방법 | |
Moldoveanu et al. | Aspects Regarding the Unguided Rocket Effectiveness | |
Bhange et al. | Technique to improve precision of kinetic energy projectiles through motion study | |
Jiang et al. | Dynamic multiprojectile attack and killing effects of detonation warheads | |
Liu et al. | Monte Carlo-Based Analysis and Experimental Validation of the Interception-Damage Probability of the New Active Interception Net | |
Selivanov et al. | Explosive propelling device for testing a debris protection of missiles and spacecraft in ground-based conditions | |
Massey et al. | A System-of-System Design of a Guided Projectile Mortar Defense System | |
Massey et al. | Combining Experimental Data, CFD, and 6-DOF Simulation to Develop a Guidance Actuator for a Supersonic Projectile | |
Herrmann | Air-to-air missile engagement analysis using the USAF Trajectory Analysis Program (TRAP) | |
Vedernikov et al. | Evolution and comparative analysis of group armor and aeroballistics of ancient and modern ruled and poly-wedge arrows | |
Wang et al. | Damage Simulation Analysis of Discrete Rods to SAM based on FEA | |
Fang et al. | Analysis on the damage efficiency caused by image terminal guided rocket to fort targets | |
Xiaochen et al. | Damage evaluation model and simulation study on functional systems of ultra low altitude aircraft |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |