CN114722573A - 一种多种类多枚杀爆弹对面目标毁伤评估方法 - Google Patents
一种多种类多枚杀爆弹对面目标毁伤评估方法 Download PDFInfo
- Publication number
- CN114722573A CN114722573A CN202210246280.9A CN202210246280A CN114722573A CN 114722573 A CN114722573 A CN 114722573A CN 202210246280 A CN202210246280 A CN 202210246280A CN 114722573 A CN114722573 A CN 114722573A
- Authority
- CN
- China
- Prior art keywords
- damage
- target
- explosive
- grid
- fragment
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000012634 fragment Substances 0.000 claims abstract description 74
- 239000002360 explosive Substances 0.000 claims abstract description 54
- 238000004364 calculation method Methods 0.000 claims abstract description 42
- 238000011156 evaluation Methods 0.000 claims abstract description 39
- 230000000694 effects Effects 0.000 claims abstract description 21
- 238000004088 simulation Methods 0.000 claims abstract description 16
- 230000035939 shock Effects 0.000 claims description 33
- 238000004880 explosion Methods 0.000 claims description 30
- 238000005070 sampling Methods 0.000 claims description 9
- 230000003068 static effect Effects 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 7
- 102100037651 AP-2 complex subunit sigma Human genes 0.000 claims description 6
- 101000806914 Homo sapiens AP-2 complex subunit sigma Proteins 0.000 claims description 6
- 230000009471 action Effects 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 3
- 238000005422 blasting Methods 0.000 claims description 2
- 238000004445 quantitative analysis Methods 0.000 claims description 2
- 238000011158 quantitative evaluation Methods 0.000 claims description 2
- 229910000831 Steel Inorganic materials 0.000 claims 1
- 239000000463 material Substances 0.000 claims 1
- 239000010959 steel Substances 0.000 claims 1
- 238000012512 characterization method Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 206010054949 Metaplasia Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 108010069898 fibrinogen fragment X Proteins 0.000 description 1
- -1 hydrogen Chemical class 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开的一种多种类多枚杀爆弹对面目标毁伤评估方法，属于杀爆战斗部对面目标毁伤效果评估领域。本发明将面目标划分成用于毁伤评估仿真的结构化网格，并对编号后的网格面赋予物理属性；在计算破片场与面目标网格微元交会情况时通过撞击点坐标计算得到各枚破片所在的网格编号，无需对每个网格与所有破片交会情况进行计算，提高评估效率。本发明能够充分考虑多枚杀爆弹与面目标弹目交会条件与不同杀爆弹的真实威力场参数，同时将面目标易损性赋予网格属性中，支持多杀爆弹耦合毁伤效果计算，提高毁伤幅员评估的合理性、精确性与实用性；此外通过引入结构化网格使多枚杀爆弹打击面目标精确毁伤幅员计算耗时难题得到解决，显著提高毁伤评估效率。
Description
技术领域
本发明属于杀爆战斗部对面目标毁伤效果评估领域，涉及一种基于结构化网格的多种类多枚杀爆战斗部对面目标毁伤幅员快速评估方法。
背景技术
毁伤评估主要用于解决弹药毁伤效能和目标毁伤效果的评价与估量的问题，尤其是近些年来，面向实战的演习、训练需求强烈以及信息化作战对毁伤方案精准定制提出了新的要求，必然牵引毁伤评估技术成为未来重要的、前沿的和学科交叉的发展领域。可以预见，随着毁伤效能评估技术的不断发展，必然可以产生“设计战争”的效果，提升作战效能，促进未来战争转型。毁伤评估主要涉及战斗部威力评估、弹药毁伤效能评估和目标毁伤效果评估三个方面。本发明提出的杀爆弹毁伤幅员计算方法可为杀爆弹动爆威力评估、多种类多枚杀爆弹联合打击毁伤效能评估、杀爆弹对面目标毁伤效果评估提供有力支撑，填补多杀爆弹耦合毁伤效果的评价方法空白。
目前，杀爆弹威力多通过破片质量、数量、速度及其分布和冲击波超压及其分布，即威力场特征参数进行描述，体现为威力数据集合的形式。这样的威力表征方法能够在一定程度上反映杀爆弹的毁伤性能，并间接地反映其毁伤目标的能力，但无法实现以归一化的度量指标定量表征与评定其综合威力，也就无法从定量的角度对比分析不同杀爆弹对同一目标以及同一杀爆弹对不同目标的毁伤能力差别。
对于已有的综合性威力度量指标——杀伤半径(密集杀伤半径和有效杀伤半径)和杀伤面积：前者只针对人员目标，没有考虑冲击波毁伤效应，另外事实上也可能存在两种弹杀伤半径相同而杀伤半径内外的毁伤威力存在着较大差别的现象，因此具有很大的局限性；后者是从目标分布和杀伤目标数量的角度定义，用于表征与度量战斗部本征功能和综合威力的含义不明确，特别是空中爆炸和毁伤空中目标需要考虑三维威力场结构时，则无法给出答案。
本发明提出的采用毁伤幅员定量表征杀爆弹综合威力的原理和方法，给出了基于战斗部威力场模型和目标易损性模型求解毁伤幅员的方法与模型，可对杀爆弹综合威力评估提供可靠指标依据。
目前已有的基于微元法的杀爆战斗部对面目标毁伤面积计算方法是先建立面目标的等效模型，划分网格后，使用遍历的方法，分别求每个网格面与破片迹线的交点、距离炸点位置等，以此来求面目标各网格微元上毁伤元数据(破片枚数、动能、冲击波超压峰值等)。这种方法在网格数量较多的时候，每个网格都要与破片轨迹线求交点，且每个网格都要判断是否达到毁伤阈值，因此毁伤评估过程会耗费大量的时间及计算资源。
以往的多枚杀爆弹毁伤面积快速计算方法为了节省时间以及计算资源，都采用毁伤半径的方式简化等效杀爆弹的威力场，此简化方法对于冲击波场的等效具备一定可靠性，但无法反应真实的破片场散布(月牙形、元宝形等)，故毁伤幅员计算结果缺乏精确性以及可靠性，无法满足精确毁伤评估需求。
发明内容
本发明主要目的是提供一种多种类多枚杀爆弹对面目标毁伤评估方法，能够充分考虑多枚杀爆弹与面目标弹目交会条件与不同杀爆弹的真实威力场参数，同时将面目标易损性赋予网格属性中，提高评估结果的合理性、精确性与实用性；此外通过引入结构化网格使多枚杀爆弹打击面目标精确毁伤幅员计算耗时难题得到解决，进而显著提高毁伤评估效率。
本发明的目的是通过下述技术方案实现的。
本发明公开的一种多种类多枚杀爆弹对面目标毁伤评估方法，包括如下步骤：
步骤1：建立面目标几何模型，并对建立的面目标几何模型进行目标区域微元划分。
建立面目标的几何模型。建立目标坐标系，将坐标系原点定在等效矩形的中心点，Y轴垂直于地面，X轴指向等效矩形较长边的方向，若两边长度相等则任取一边作为长边，Z轴通过右手定则获得。面目标为长宽分别为L、W的等效矩形阵地。
将面目标划分成用于毁伤评估仿真的结构化网格，对网格按照XZ平面位置顺序进行编号，赋予网格“毁伤状态”属性以及“目标类型”属性。目标区域微元化划分相关参量如下：
Δaim＝f(La,Wa,Na1,Na2,Δdamagecriterion) (1)
其中：La为目标区域长度；Wa为目标区域宽度；Na1表示在长度方向划分网格大小；Na2表示在宽度方向划分网格大；Δdamagecriterion为不同种类目标毁伤判据，所述不同种类目标毁伤判据包括破片比动能、冲击波超压峰值。
作为优选，目标微元的毁伤判据数据通过查阅文献、数值仿真或试验验证获取，此种赋予微元属性的方法扩展性强，能够简化目标易损性，且能够精细化区别目标各部位易损性，即改变微元对应目标类型即可改变对应毁伤判据。
步骤2：根据弹目交会条件与引信参数计算杀爆弹炸点坐标。
根据命中精度和末端弹道参数抽样得到末端弹道线方程，由于杀爆战斗部的末端弹道近似为一条直线，因此根据弹道线方向向量DPene＝(DX,DY,DZ)与线上一点PL＝(XL,YL,ZL)确定目标坐标系下的弹道线方程。
其中弹道线方向向量DPene＝(DX,DY,DZ)通过杀爆弹落角α与方位角β计算得到。
其中，入射角α是指末端弹道线与XOZ平面的夹角，取值为[0,90°]；方位角β是指末端弹道线在XOZ平面上的投影线与X轴的夹角，取值为[-180°,180°]，逆时针为正。
线上一点PL＝(XL,YL,ZL)通过瞄准点PAim＝(XAim,YAim,ZAim)与命中精度参数CEP进行满足均值为(XAim,ZAim)，方差为(σ1，σ2)的二维正态分布抽样计算得到。考虑弹药CEP，使毁伤幅员计算更贴合实际作战。其中σ1，σ2相等，计算公式为：
σ1＝σ2＝σ＝CEP/1.1774 (3)
根据已有参数能够构建两枚杀爆弹的末端弹道线方程，用于后续的基于引信参数的炸点位置计算。
根据输入杀爆战斗部的实际参数设定模型的引信参数和杀爆弹末端弹道线方程，计算得到杀爆弹炸点坐标PQ＝(XQ,YQ,ZQ)。
步骤3：计算目标区域微元上的毁伤元载荷数据。
通过文件读入各枚杀爆弹的静爆威力场数据，其中静爆威力场数据文件包含破片编号、弹体坐标系下每枚破片XYZ方向上的初始速度与初始位置、破片质量、破片速度衰减系数、炸点不同距离处的冲击波超压峰值与比冲量。
威力场数据文件中的参数通过经验公式、数值仿真或试验收集获得，通过标准格式威力场文件使得能够同时计算多枚不同种类杀爆弹打击效果，并且数据来源广泛，拓展性强。
当杀爆弹为威力场相同的同一型杀爆弹，威力场数据来源于经验公式计算。若为不同类型杀爆弹则只需读取不同的相对应的威力场文件即能够实现多种类杀爆弹协同毁伤幅员评估。
引入标准威力数据文件与采用毁伤半径等效战斗部毁伤威力相比，能够使毁伤效果仿真更接近真实打击情况，毁伤幅员评估结果更具实际参考价值，为“察-打-评”一体化作战提供有力支撑。
将威力场数据里的破片在弹体坐标系下的初始位置以及飞散方向转换到目标坐标系下。
弹体速度矢量与破片静爆飞散速度矢量相叠加即能够得到每枚破片在战斗部动爆情况下速度大小与方向。
对于破片的运动轨迹能够近似为一射线，即已知初始点与飞散方向能够构造射线方程，计算该射线与面目标的交会情况，并根据交会点坐标(射线与四边形求交得到)结合结构化网格特性计算出交会网格微元编号。具体计算方法如下：
其中iE为网格长方向上编号，jE为网格宽方向上编号，L为面目标区域长，W为面目标区域宽，DL为长方向上网格划分数量，DW为宽方向上网格划分数量，XExp为交会点X坐标，ZExp为交会点Z坐标。
引入结构化网格后每枚破片只需计算一次轨迹线与面目标交点即可根据交点坐标求得交会网格微元编号，与以往传统的网格法计算毁伤幅员需遍历计算轨迹线与所有网格微元(四边形)交点相比节省大量计算资源(网格划分越细计算效率提升越大)。将精确毁伤幅员计算时间缩减到秒级(传统考虑真实动爆威力场的毁伤幅员计算方法为分钟/小时级)。
根据每枚破片的初始速度大小、飞行距离以及速度衰减系数能够通过速度衰减公式计算得到破片作用与目标时的速度大小。
Vx＝V0e-αx (5)
其中Vx为剩余速度，V0为初始速度，α为破片速度衰减系数，x为破片飞行距离。
根据冲击波威力场参数(炸点不同距离出冲击波超压峰值与比冲量)以及各杀爆弹炸点坐标通过线性插值计算得到各微元上冲击波威力参数，即计算得到目标区域微元上的毁伤元载荷数据。
步骤4：根据目标易损性数据(文件读入)以及微元的“目标类型”属性获得网格微元上的破片与冲击波毁伤阈值，并选取毁伤准则。
作为优选，采用破片动能与冲击波超压峰值作为毁伤准则，或只需更改易损性文件格式拓展采用破片质量-速度联合判据、破片动能-枚数联合判据、冲击波超压-比冲量联合判据作为毁伤准则。
步骤5：根据步骤3计算得到作用在微元上的毁伤元载荷数据与毁伤阈值作比较判断微元是否毁伤。
作为优选，微元“毁伤状态”属性取值为0或1，其中0代表未被毁伤，1代表被毁伤，初始状态全部设为0。
步骤6：求和所有被毁伤网格微元总面积，计算火箭杀爆弹对目标区域毁伤幅员与毁伤比。
列于式(2)、式(3)，式中SDi为每个达到毁伤判据的网格面积，Saim为目标区域总面积。
ΔDamageArea＝∑SDi (6)
其中SDi为每个达到毁伤判据的网格面积
其中Saim为目标区域总面积。
步骤7：将多次抽样计算获得的毁伤幅员求期望得到弹药对目标的毁伤幅员，并得出毁伤幅员方差，即实现多种类多枚杀爆弹对面目标毁伤评估。
还包括步骤8：基于弹药真实威力场，根据步骤8得到的杀爆弹对面目标毁伤评估结果，能够作为杀爆弹综合威力定量评估与分析指标、杀爆弹对不同类型面目标毁伤能力评价指标、面目标毁伤效果评价与估量指标，用于支撑杀爆弹的毁伤规划或毁伤方案制订和目标毁伤效果预测，不仅能够支撑精确打击，高效毁伤，还能够优化弹药威力性能，为打击目标选择合理弹药提供依据。
有益效果：
1、本发明公开的一种多种类多枚杀爆弹对面目标毁伤评估方法，将面目标划分成用于毁伤评估仿真的结构化网格，对网格按照二维空间位置顺序进行编号，赋予网格“毁伤状态”以及“易损性参数”属性，并对每个网格的面、线、点分别进行编号，并对编号后的网格面赋予物理属性。在计算破片场与面目标网格微元交会情况时通过撞击点坐标计算得到各枚破片所在的网格编号，无需对每个网格与所有破片交会情况进行计算，节省计算资源，进而提高评估效率。
2、本发明公开的一种多种类多枚杀爆弹对面目标毁伤评估方法，采用弹体坐标系下每枚破片X、Y、Z方向上的初始速度与初始位置、破片质量、破片速度衰减系数作为破片威力场参数输入，距离炸点不同距离下的冲击波超压峰值、比冲量作为冲击波威力场参数输入，对杀爆弹威力场参数进行综合全面考虑。与传统的采用战斗部毁伤半径表征威力场的毁伤面积/幅员计算方法相比对杀爆弹的威力评估更加全面可靠。
3、本发明公开的一种多种类多枚杀爆弹对面目标毁伤评估方法，目标易损性数据通过目标类型、破片比动能阈值、冲击波超压峰值阈值等参数输入算法，赋予面目标网格微元目标类型参数计算时即可获取其毁伤判据参数，对每一微元赋予易损性属性可拓展面目标更精细的目标易损性(不同部位易损程度不同)，使幅员评估结果更具实用性与合理性。
4、本发明公开的一种多种类多枚杀爆弹对面目标毁伤评估方法，能够为作战时的火力打击提供一个预期的毁伤结果，为后续的作战方案提供支持。另外，通过设置不同的初始参数，利用每次仿真得到的毁伤幅员结果进行对比，能够在作战前为多杀爆战斗部的火力规划提供支撑，优化打击方案。
5、本发明公开的一种多种类多枚杀爆弹对面目标毁伤评估方法，支持多杀爆弹耦合毁伤效果计算，计算结果相对于杀伤半径等传统战斗部指标更具实战参考价值，可为弹药威力性能论证提供支撑。
6、本发明公开的一种多种类多枚杀爆弹对面目标毁伤评估方法，引入结构化网格后每枚破片只需计算一次轨迹线与面目标交点即能够根据交点坐标求得交会网格微元编号，与以往传统的网格法计算毁伤面积需遍历计算轨迹线与所有网格微元(四边形)交点相比节省大量计算资源(网格划分约细计算效率提升约大)，显著提升毁伤评估效率，能够将精确毁伤幅员计算时间缩减到秒级。
附图说明
图1为本发明公开的一种多种类多枚杀爆弹对面目标毁伤评估方法流程图。
图2为弹目交会分析流程图。
图3为面目标结构化网格划分及编号示意图。
图4为杀爆弹威力场散布以及网格毁伤状态示意图。
具体实施方式
下面结合附图实例对本发明进行详细描述
如图1所示，本实施例公开的一种多种类多枚杀爆弹对面目标毁伤评估方法，计算实例参数如下表所示。
本实施例公开的一种多种类多枚杀爆弹对面目标毁伤评估方法，具体实现步骤如下：
步骤1建立面目标几何模型以及目标区域微元划分，如图3所示，具体如下：
建立面目标的几何模型。建立目标坐标系，将坐标系原点定在等效矩形的中心点，Y轴垂直于地面，X轴指向等效矩形较长边的方向，若两边长度相等则任取一边作为长边，Z轴通过右手定则获得。此实例中面目标为长宽分别为L、W的等效矩形阵地，其中L＝100m，W＝50m。
将面目标划分成用于毁伤评估仿真的结构化网格，对网格按照XZ平面位置顺序进行编号，赋予网格“毁伤状态”属性以及“目标类型”属性。目标区域微元化划分相关参量如下：
Δaim＝f(La,Wa,Na1,Na2,Δdamagecriterion) (1)
其中：La为目标区域长度；Wa为目标区域宽度；Na1表示在长度方向划分网格大小；Na2表示在宽度方向划分网格大；Δdamagecriterion为不同种类目标毁伤判据(破片比动能、冲击波超压峰值等)。
此实例中采用网格划分参数为：La＝100m，Wa＝50m，Na1＝1m，Na2＝1m，目标类型为“轻型装甲车辆集群”，对应破片动能阈值采用2160J，冲击波超压峰值阈值采用0.02Mpa。目标微元的毁伤判据数据可通过查阅文献、数值仿真、试验验证等方式获取，此种赋予微元属性的方法扩展性强，可简化目标易损性也可精细化区别目标各部位易损性(改变微元对应目标类型即可改变对应毁伤判据，例如一装甲车辆前部与顶部抗破片侵彻能力不同)。
步骤2：如图2所示根据弹目交会条件与引信参数计算杀爆弹炸点坐标，具体如下：
根据命中精度和末端弹道参数抽样得到末端弹道线方程，由于杀爆战斗部的末端弹道近似为一条直线，因此根据弹道线方向向量DPene＝(DX,DY,DZ)与线上一点PL＝(XL,YL,ZL)确定目标坐标系下的弹道线方程。
其中弹道线方向向量DPene＝(DX,DY,DZ)通过杀爆弹落角α与方位角β计算得到。
其中，入射角α是指末端弹道线与XOZ平面的夹角，取值为[0,90°]；方位角β是指末端弹道线在XOZ平面上的投影线与X轴的夹角，取值为[-180°,180°]，逆时针为正。
此实例中杀爆弹1落角α1＝60°，方位角β1＝90°，杀爆弹2落角α1＝75°，方位角β1＝0°。
线上一点PL＝(XL,YL,ZL)通过瞄准点PAim＝(XAim,YAim,ZAim)与命中精度参数CEP进行满足均值为(XAim,ZAim)，方差为(σ1，σ2)的二维正态分布抽样计算得到。考虑弹药CEP，可以使毁伤幅员计算更贴合实际作战，结果更加具备实际参考价值。其中σ1，σ2相等，计算公式为：
σ1＝σ2＝σ＝CEP/1.1774 (3)
此实例中两枚杀爆弹瞄准点分别为(0，-5，0)，(10，5，0)(均为目标坐标系下坐标，坐标原点为面目标中心)。根据已有参数可以构建两枚杀爆弹的末端弹道线方程，用于后续的基于引信参数的炸点位置计算。
根据输入杀爆战斗部的实际参数设定模型的引信参数(炸高)和杀爆弹末端弹道线方程计算杀爆弹炸点坐标PQ＝(XQ,YQ,ZQ)。
此实例中杀爆弹1炸高为4m，杀爆弹2炸高为8m，根据两枚杀爆弹引信参数(炸高)与已有末端弹道线方程，即可得到两杀爆弹实际炸点坐标。
步骤3：计算目标区域微元上的毁伤元载荷数据，具体计算步骤如下：
通过文件读入各枚杀爆弹的静爆威力场数据，其中静爆威力场数据文件包含破片编号、弹体坐标系下每枚破片XYZ方向上的初始速度与初始位置、破片质量、破片速度衰减系数、炸点不同距离处的冲击波超压峰值与比冲量。
威力场数据文件中的参数可通过经验公式、数值仿真、试验收集等途径获得，通过此种标准格式威力场文件使得能够同时计算多枚不同种类杀爆弹打击效果，并且数据来源广泛，可拓展性强。
此实例中两枚杀爆弹为同一型杀爆弹(威力场相同)，威力场数据来源于经验公式计算。若为不同类型杀爆弹则只需读取不同的相对应的威力场文件即可实现多种类杀爆弹协同毁伤幅员计算。
引入标准威力数据文件与采用毁伤半径等效战斗部毁伤威力相比可使毁伤效果仿真更接近真实打击情况，毁伤幅员计算结果更具实际参考价值，可为“察-打-评”一体化作战提供有力支撑。
将威力场数据里的破片在弹体坐标系下的初始位置以及飞散方向转换到目标坐标系下。
弹体速度矢量与破片静爆飞散速度矢量相叠加即可得到每枚破片在战斗部动爆情况下速度大小与方向。
对于破片的运动轨迹可以近似为一射线(已知初始点与飞散方向可构造射线方程)，计算该射线与面目标的交会情况，如图4所示并可根据交会点坐标(射线与四边形求交得到)结合结构化网格特性计算出交会网格微元编号。具体计算方法如下：
其中iE为网格长方向上编号，jE为网格宽方向上编号，L为面目标区域长，W为面目标区域宽，DL为长方向上网格划分数量，DW为宽方向上网格划分数量，XExp为交会点X坐标，ZExp为交会点Z坐标。
引入结构化网格后每枚破片只需计算一次轨迹线与面目标交点即可根据交点坐标求得交会网格微元编号，与以往传统的网格法计算毁伤幅员需遍历计算轨迹线与所有网格微元(四边形)交点相比节省大量计算资源(网格划分约细计算效率提升约大)。将精确毁伤幅员计算时间缩减到秒级(传统考虑真实动爆威力场的毁伤幅员计算方法为分钟/小时级)。
根据每枚破片的初始速度大小、飞行距离以及速度衰减系数可以通过速度衰减公式计算得到破片作用与目标时的速度大小。
Vx＝V0e-αx (5)
其中Vx为剩余速度，V0为初始速度，α为破片速度衰减系数，x为破片飞行距离。
根据冲击波威力场参数(炸点不同距离出冲击波超压峰值与比冲量)以及各杀爆弹炸点坐标通过线性插值计算得到各微元上冲击波威力参数。
此实例中杀爆弹毁伤元(破片、冲击波)与面目标(轻型装甲车辆集群)交会效果如图4所示。
步骤4：根据目标易损性数据(文件读入)以及微元的“目标类型”属性获得网格微元上的破片与冲击波毁伤阈值，本发明采用破片动能与冲击波超压峰值作为毁伤准则，可拓展采用破片质量-速度联合判据、破片动能-枚数联合判据、冲击波超压-比冲量联合判据等只需更改易损性文件格式即可。
此实例中轻型装甲车辆集群破片动能阈值采用2160J，冲击波超压峰值阈值采用0.02Mpa(均来自文献调研)。
步骤5：根据步骤3计算得到作用在微元上的毁伤元载荷数据与毁伤阈值作比较判断微元是否毁伤，微元“毁伤状态”属性取值为0或1，其中0代表未被毁伤，1代表被毁伤，初始状态全部设为0。
步骤6：求和所有被毁伤网格微元总面积，计算火箭杀爆弹对目标区域毁伤幅员与毁伤比，列于式(2)、式(3)，式中SDi为每个达到毁伤判据的网格面积，Saim为目标区域总面积。
ΔDamageArea＝∑SDi (6)
其中SDi为每个达到毁伤判据的网格面积
其中Saim为目标区域总面积。
步骤7将多次抽样计算获得的毁伤幅员求期望得到弹药对目标的毁伤幅员，并给出方差。
此实例中两枚杀爆弹联合打击下对该面目标造成毁伤幅员为1094m2，耗时0.5秒。计算时耗满足现代作战指挥需求。
步骤8：基于弹药真实威力场，根据步骤8得到的杀爆弹对面目标毁伤评估结果，通过此实例可以说明本发明在实际作战中当侦察已知敌方目标大概区域范围后通过快速仿真计算得出毁伤幅员与毁伤范围，基于弹药真实威力场可为指挥系统提供合理可靠的毁伤效果预测，为下一步打击决策提供合理支撑，符合现代战争中精确打击，高效毁伤的作战理念。另一方面可为弹药威力性能优化提供参考，为打击目标合理弹药选择提供依据。
以上所述的具体描述，对发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (5)
1.一种多种类多枚杀爆弹对面目标毁伤评估方法，其特征在于：包括如下步骤，
步骤1：建立面目标几何模型，并对建立的面目标几何模型进行目标区域微元划分；
建立面目标的几何模型；建立目标坐标系，将坐标系原点定在等效矩形的中心点，Y轴垂直于地面，X轴指向等效矩形较长边的方向，若两边长度相等则任取一边作为长边，Z轴通过右手定则获得；面目标为长宽分别为L、W的等效矩形阵地；
将面目标划分成用于毁伤评估仿真的结构化网格，对网格按照XZ平面位置顺序进行编号，赋予网格“毁伤状态”属性以及“目标类型”属性；目标区域微元化划分相关参量如下：
Δaim＝f(La,Wa,Na1,Na2,Δdamagecriterion) (1)
其中：La为目标区域长度；Wa为目标区域宽度；Na1表示在长度方向划分网格大小；Na2表示在宽度方向划分网格大；Δdamagecriterion为不同种类目标毁伤判据，所述不同种类目标毁伤判据包括破片比动能、冲击波超压峰值；
步骤2：根据弹目交会条件与引信参数计算杀爆弹炸点坐标；
根据命中精度和末端弹道参数抽样得到末端弹道线方程，由于杀爆战斗部的末端弹道近似为一条直线，因此根据弹道线方向向量DPene＝(DX,DY,DZ)与线上一点PL＝(XL,YL,ZL)确定目标坐标系下的弹道线方程；
其中弹道线方向向量DPene＝(DX,DY,DZ)通过杀爆弹落角α与方位角β计算得到；
其中，入射角α是指末端弹道线与XOZ平面的夹角，取值为[0,90°]；方位角β是指末端弹道线在XOZ平面上的投影线与X轴的夹角，取值为[-180°,180°]，逆时针为正；
线上一点PL＝(XL,YL,ZL)通过瞄准点PAim＝(XAim,YAim,ZAim)与命中精度参数CEP进行满足均值为(XAim,ZAim)，方差为(σ1，σ2)的二维正态分布抽样计算得到；考虑弹药CEP，使毁伤幅员计算更贴合实际作战；其中σ1，σ2相等，计算公式为：
σ1＝σ2＝σ＝CEP/1.1774 (3)
根据已有参数能够构建杀爆弹的末端弹道线方程，用于后续的基于引信参数的炸点位置计算；
根据输入杀爆战斗部的实际参数设定模型的引信参数和杀爆弹末端弹道线方程，计算得到杀爆弹炸点坐标PQ＝(XQ,YQ,ZQ)；
步骤3：计算目标区域微元上的毁伤元载荷数据；
通过文件读入各枚杀爆弹的静爆威力场数据，其中静爆威力场数据文件包含破片编号、弹体坐标系下每枚破片XYZ方向上的初始速度与初始位置、破片质量、破片速度衰减系数、炸点不同距离处的冲击波超压峰值与比冲量；
威力场数据文件中的参数通过经验公式、数值仿真或试验收集获得，通过标准格式威力场文件使得能够同时计算多枚不同种类杀爆弹打击效果，并且数据来源广泛，拓展性强；
当杀爆弹为威力场相同的同一型杀爆弹，威力场数据来源于经验公式计算；若为不同类型杀爆弹则只需读取不同的相对应的威力场文件即能够实现多种类杀爆弹协同毁伤幅员评估；
引入标准威力数据文件与采用毁伤半径等效战斗部毁伤威力相比，能够使毁伤效果仿真更接近真实打击情况，毁伤幅员评估结果更具实际参考价值，为“察-打-评”一体化作战提供有力支撑；
将威力场数据里的破片在弹体坐标系下的初始位置以及飞散方向转换到目标坐标系下；
弹体速度矢量与破片静爆飞散速度矢量相叠加即能够得到每枚破片在战斗部动爆情况下速度大小与方向；
对于破片的运动轨迹能够近似为一射线，即已知初始点与飞散方向能够构造射线方程，计算该射线与面目标的交会情况，并根据交会点坐标结合结构化网格特性计算出交会网格微元编号；具体计算方法如下：
其中iE为网格长方向上编号，jE为网格宽方向上编号，L为面目标区域长，W为面目标区域宽，DL为长方向上网格划分数量，DW为宽方向上网格划分数量，XExp为交会点X坐标，ZExp为交会点Z坐标；
引入结构化网格后每枚破片只需计算一次轨迹线与面目标交点即可根据交点坐标求得交会网格微元编号，与以往传统的网格法计算毁伤幅员需遍历计算轨迹线与所有网格微元交点相比节省大量计算资源；将精确毁伤幅员计算时间缩减到秒级；
根据每枚破片的初始速度大小、飞行距离以及速度衰减系数能够通过速度衰减公式计算得到破片作用与目标时的速度大小；
Vx＝V0e-αx (5)
其中Vx为剩余速度，V0为初始速度，α为破片速度衰减系数，x为破片飞行距离；
根据冲击波威力场参数(炸点不同距离出冲击波超压峰值与比冲量)以及各杀爆弹炸点坐标通过线性插值计算得到各微元上冲击波威力参数，即计算得到目标区域微元上的毁伤元载荷数据；
步骤4：根据目标易损性数据以及微元的“目标类型”属性获得网格微元上的破片与冲击波毁伤阈值，并选取毁伤准则；
步骤5：根据步骤3计算得到作用在微元上的毁伤元载荷数据与毁伤阈值作比较判断微元是否毁伤；
步骤6：求和所有被毁伤网格微元总面积，计算火箭杀爆弹对目标区域毁伤幅员与毁伤比；
列于式(2)、式(3)，式中SDi为每个达到毁伤判据的网格面积，Saim为目标区域总面积；
ΔDamageArea＝∑SDi (6)
其中SDi为每个达到毁伤判据的网格面积
其中Saim为目标区域总面积；
步骤7：将多次抽样计算获得的毁伤幅员求期望得到弹药对目标的毁伤幅员，并得出毁伤幅员方差，即实现多种类多枚杀爆弹对面目标毁伤评估。
2.如权利要求1所述的一种多种类多枚杀爆弹对面目标毁伤评估方法，其特征在于：还包括步骤8：基于弹药真实威力场，根据步骤8得到的杀爆弹对面目标毁伤评估结果，能够作为杀爆弹综合威力定量评估与分析指标、杀爆弹对不同类型面目标毁伤能力评价指标、面目标毁伤效果评价与估量指标，用于支撑杀爆弹的毁伤规划或毁伤方案制订和目标毁伤效果预测，不仅能够支撑精确打击，高效毁伤，还能够优化弹药威力性能，为打击目标选择合理弹药提供依据。
3.如权利要求1或2所述的一种多种类多枚杀爆弹对面目标毁伤评估方法，其特征在于：目标微元的毁伤判据数据通过查阅文献、数值仿真或试验验证获取，此种赋予微元属性的方法扩展性强，能够简化目标易损性，且能够精细化区别目标各部位易损性，即改变微元对应目标类型即可改变对应毁伤判据。
4.如权利要求1或2所述的一种多种类多枚杀爆弹对面目标毁伤评估方法，其特征在于：采用破片动能与冲击波超压峰值作为毁伤准则，或只需更改易损性文件格式拓展采用破片质量-速度联合判据、破片动能-枚数联合判据、冲击波超压-比冲量联合判据作为毁伤准则。
5.如权利要求1或2所述的一种多种类多枚杀爆弹对面目标毁伤评估方法，其特征在于：微元被赋予多种属性，包含“毁伤状态”属性，取值为0或1，其中0代表未被毁伤，1代表被毁伤，初始状态全部设为0，“目标类型”属性，对应易损性文件中目标类型编号，并可拓展赋予“结构材料”属性，如“装甲钢”、“混凝土”。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210246280.9A CN114722573A (zh) | 2022-03-14 | 2022-03-14 | 一种多种类多枚杀爆弹对面目标毁伤评估方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210246280.9A CN114722573A (zh) | 2022-03-14 | 2022-03-14 | 一种多种类多枚杀爆弹对面目标毁伤评估方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114722573A true CN114722573A (zh) | 2022-07-08 |
Family
ID=82236786
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210246280.9A Pending CN114722573A (zh) | 2022-03-14 | 2022-03-14 | 一种多种类多枚杀爆弹对面目标毁伤评估方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114722573A (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115270528A (zh) * | 2022-09-28 | 2022-11-01 | 中国兵器科学研究院 | 一种分布弹药的方法、装置、设备和介质 |
| CN115661209A (zh) * | 2022-12-27 | 2023-01-31 | 中国人民解放军国防科技大学 | 基于事件相机的爆炸威力场毁伤元测量方法及系统 |
| CN116364272A (zh) * | 2023-03-06 | 2023-06-30 | 南京理工大学 | 一种面向爆炸的人员损伤测试评估方法 |
| CN117113787A (zh) * | 2023-10-25 | 2023-11-24 | 北京理工大学 | 一种含能破片冲击-释能耦合效应的确定方法 |
| CN117150825A (zh) * | 2023-10-31 | 2023-12-01 | 北京理工大学 | 一种装甲类目标最大易损方向的获取方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110955953A (zh) * | 2019-07-25 | 2020-04-03 | 北京理工大学 | 基于结构化网格的多种类侵爆弹对建筑目标毁伤评估方法 |
| CN112685944A (zh) * | 2021-01-07 | 2021-04-20 | 北京理工大学 | 一种同类多枚杀爆弹对面目标瞄准点规划方法 |
| CN113609574A (zh) * | 2021-07-08 | 2021-11-05 | 中国人民解放军火箭军工程大学 | 基于等毁线的多轴特种车辆轮胎系统抗毁伤能力表征方法 |
-
2022
- 2022-03-14 CN CN202210246280.9A patent/CN114722573A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110955953A (zh) * | 2019-07-25 | 2020-04-03 | 北京理工大学 | 基于结构化网格的多种类侵爆弹对建筑目标毁伤评估方法 |
| CN112685944A (zh) * | 2021-01-07 | 2021-04-20 | 北京理工大学 | 一种同类多枚杀爆弹对面目标瞄准点规划方法 |
| CN113609574A (zh) * | 2021-07-08 | 2021-11-05 | 中国人民解放军火箭军工程大学 | 基于等毁线的多轴特种车辆轮胎系统抗毁伤能力表征方法 |
Non-Patent Citations (3)
| Title |
|---|
| 李向荣;王国辉;徐峰;: "杀爆弹对主战坦克毁伤评估仿真", 装甲兵工程学院学报, no. 02, 15 April 2009 (2009-04-15) * |
| 梁斌;冯高鹏;周婕群;李俊承;: "离散杆和EFP组合战斗部对飞机毁伤效应研究", 弹箭与制导学报, no. 01, 15 February 2017 (2017-02-15) * |
| 高鹏;徐豫新;李可;杨晓红;王树山;: "火箭杀爆弹动爆毁伤幅员计算方法与可视化仿真实现", 兵工学报, no. 2, 15 December 2016 (2016-12-15) * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115270528A (zh) * | 2022-09-28 | 2022-11-01 | 中国兵器科学研究院 | 一种分布弹药的方法、装置、设备和介质 |
| CN115661209A (zh) * | 2022-12-27 | 2023-01-31 | 中国人民解放军国防科技大学 | 基于事件相机的爆炸威力场毁伤元测量方法及系统 |
| CN116364272A (zh) * | 2023-03-06 | 2023-06-30 | 南京理工大学 | 一种面向爆炸的人员损伤测试评估方法 |
| CN116364272B (zh) * | 2023-03-06 | 2023-11-28 | 南京理工大学 | 一种面向爆炸的人员损伤测试评估方法 |
| CN117113787A (zh) * | 2023-10-25 | 2023-11-24 | 北京理工大学 | 一种含能破片冲击-释能耦合效应的确定方法 |
| CN117113787B (zh) * | 2023-10-25 | 2024-01-05 | 北京理工大学 | 一种含能破片冲击-释能耦合效应的确定方法 |
| CN117150825A (zh) * | 2023-10-31 | 2023-12-01 | 北京理工大学 | 一种装甲类目标最大易损方向的获取方法 |
| CN117150825B (zh) * | 2023-10-31 | 2024-02-13 | 北京理工大学 | 一种装甲类目标最大易损方向的获取方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114722573A (zh) | 一种多种类多枚杀爆弹对面目标毁伤评估方法 | |
| CN110955953B (zh) | 基于结构化网格的多种类侵爆弹对建筑目标毁伤评估方法 | |
| CN104143028B (zh) | 一种基于仿真实验的装甲装备毁伤规律分析方法 | |
| CN102314537B (zh) | 一种掠飞击顶灵巧弹药毁伤概率计算方法 | |
| CN114492170A (zh) | 一种典型目标易损性毁伤评估方法及系统 | |
| CN109857145B (zh) | 一种基于迭代预测命中点的增程型拦截弹预测制导方法 | |
| CN109063331B (zh) | 一种小脱靶量破片战斗部对再入弹头毁伤效应的评估方法 | |
| CN104850747A (zh) | 一种弹药威力杀伤体积评价方法 | |
| KR101997387B1 (ko) | 음향 센서를 이용한 탄착점 추정 방법 및 장치 | |
| CN112464451B (zh) | 基于作战仿真系统的防空导弹武器命中概率修正方法 | |
| CN116305988A (zh) | 一种基于引信建模的战斗部威力仿真系统及方法 | |
| CN114692318A (zh) | 飞机冲击动力学测试用格栅式燃油箱结构抗毁伤评估方法 | |
| CN113609574B (zh) | 基于等毁线的多轴特种车辆轮胎系统抗毁伤能力表征方法 | |
| CN108920827A (zh) | 连发射击命中率试验方法 | |
| CN112836298A (zh) | 一种防空导弹对直升机目标毁伤效能评估模型的构建方法 | |
| Peng et al. | Terminal attitude selection method of missile attack aircraft | |
| CN115345035B (zh) | 一种体目标多向打击空间瞄准点优化方法 | |
| CN110929372A (zh) | 一种子母弹对体目标作用效能的评估方法 | |
| Chen et al. | Simulation analysis of equipment battle damage test layout design | |
| CN117150868B (zh) | 一种用于毁伤评估的有限元数值仿真模型生成方法及系统 | |
| Peng et al. | A new damage probability calculation of submunition to target based on intersection volume algorithm | |
| Qiankun et al. | Research on the measurement of muzzle shock wave pressure field for a naval gun | |
| CN112182898B (zh) | 一种动能拦截子母式再入弹头末端遭遇判别方法 | |
| Chen et al. | Epipolar multitarget velocity probability data association algorithm based on the movement characteristics of blasting fragments | |
| Zhang et al. | Mathematical calculation method for damage probability of projectile hitting ground target |
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 |