CN108984808A - 引战配合仿真交会破片命中点参数计算方法 - Google Patents
引战配合仿真交会破片命中点参数计算方法 Download PDFInfo
- Publication number
- CN108984808A CN108984808A CN201710410853.6A CN201710410853A CN108984808A CN 108984808 A CN108984808 A CN 108984808A CN 201710410853 A CN201710410853 A CN 201710410853A CN 108984808 A CN108984808 A CN 108984808A
- Authority
- CN
- China
- Prior art keywords
- fragmentation
- point
- face element
- target
- impact
- 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.)
- Withdrawn
Links
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
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
技术领域
本发明涉及引战配合仿真技术,具体涉及一种引战配合仿真交会破片命中点参数计算方法。
背景技术
引战系统可视化仿真(Visual Simulation of Fuze-warhead Coordination)是指利用计算机技术的图形化、动态化和实时化模拟在导弹和目标的交会过程中引信利用环境信息、目标信息和制导信息,控制战斗部对目标的最佳毁伤,即引战配合过程。随着现代武器系统的发展,来自空中的威胁不断升级,目标种类增多、速度加快、机动性能不断提升,使得交会状态变得越来越复杂;同时,在实际运用之中,由于导弹的研制与试验成本高,研发周期长,每次试验会耗费大量的人力物力。仅依靠外场试验不仅价格昂贵、周期长,而且有些情况几乎是无法完成的,在引战配合过程中,充分利用计算机仿真技术,以图形化的方式逼真地显示导弹和目标的末端交会时,引信探测目标、战斗部起爆、破片动态飞散和破片命中目标各部位的过程,并通过计算引战配合效率和可视化显示命中结果来分析引战系统的性能,以优化引战系统的参数设计。
近年来引战配合可视化仿真技术朝着平台化、系统化方向快速发展。如何精确计算破片战斗部爆炸后的命中点参数(包括命中点坐标、命中部位、命中数、破片落速、落角、动能和命中时刻)是一个重大课题。为了提高命中点参数的计算精度和通用性,一些命中点参数算法被提出。最基本的两种方法是计算破片数目的投影球算法和按时间步长计算的命中点参数算法,它们分别满足了引战配合可视化仿真中的算法简单和实时计算的要求。
(1)投影球方法具有概念明确、思路清晰、算法相对简单之特点,具有较强的通用性。它是基于光学原理的物理概念,以点光源的光的扩散来模拟战斗部爆炸后破片的飞散过程,点光源在目标面上的投影区域即全部命中目标的破片数,但这样的方法无法计算命中点坐标等参数,致使引战配合效率计算的精度和可信度较差。
(2)按时间步长计算的命中点参数算法,通过对破片——目标面元顶点空间位置关系的实时计算,判断该枚破片是否命中目标,若命中,则计算命中点参数。但该方法没有考虑破片速度衰减对命中点参数的影响,且该方法在每个时间步长要对所有命中点参数进行遍历判断,时间复杂度高,且未能考虑速度衰减对命中点参数的影响,对命中点参数的计算存在着一定的偏差。
发明内容
本发明的目的在于提供一种引战配合仿真交会破片命中点参数计算方法,解决现有引战配合可视化命中点参数算法存在的精度不足、效率不高、以及因速度衰减而降低计算结果准确性等问题。
实现本发明目的的技术方案为:一种引战配合仿真交会破片命中点参数计算方法,包括以下步骤:
步骤1,对目标进行面元划分,得到目标三角面元模型;
步骤2,建立破片与目标面元交会的数学模型;
步骤3,计算目标面元所在平面与破片射线交点与命中点参数;
步骤4,判定命中点是否在三角形面元内,若位于三角形面元内,则破片命中该面元。
与现有技术相比,本发明的显著优点为:本发明提出了弹目交会破片命中点参数的改进型计算模型方法,通过在考虑速度衰减条件下,利用射击线法提出了一种通过一次遍历求解所有破片命中点参数的计算模型方法,在保证计算精度的基础上,提高了计算效率。
附图说明
图1为本发明引战配合仿真交会破片命中点参数计算方法流程示意图。
具体实施方式
结合图1,一种引战配合仿真交会破片命中点参数计算方法,包括以下步骤:
步骤1,对目标进行面元划分,得到目标三角面元模型;
提取出面元的3个顶点坐标,由A、B、C三点确定面元的法向量n′=sAB×εAC,并将n′单位化得到n;其中sAB、sAC分别为AB连线、AC连线的方向矢量。
步骤2,建立破片与目标面元交会的数学模型,具体为:
以战斗部爆炸时刻为0时刻,列出破片的轨迹方程:
式(1)中,X为t时刻破片空间坐标,s(t)为破片在t时刻所飞行的距离,为破片速度方向的单位矢量,X0为0时刻破片的空间坐标;
考虑破片在空气中飞行速度的衰减,其式为:
υf=υ0exp(-αs) (2)
υf为破片在t时刻的速率,υ0为破片初速速率,α为空气密度衰减系数,s为破片在t时刻所飞行的距离;
由(1)、(2)式得到破片飞行距离与时间的关系:
S(t)=[ln(1+αv0t)]/α (3)
接下来建立目标面元所在平面的运动方程;在0时刻,目标某面元所在的平面方程可表示为:
n·X+D0=0
D0=-n·Y0 (4)
式(4)中,n为目标面元法向单位向量;D0为一标量常数,表示坐标原点到目标面元的距离;Y0为三角形面元中任意一顶点的初始坐标。由该面元在0时刻的3个顶点坐标即可求得n和D0的值。
假设面元只有平动。面元的法向不变。面元所在平面的运动方程可以表示为:
n·(X-vtt)+D0=0 (5)
式(5)中,vt为目标在遭遇段的运动速度。
步骤3,计算目标面元所在平面与破片射线交点与命中点参数;步骤3具体为:
由式(1)和式(4)联立解得:
式(6)中,h(t)为t时刻面元与破片射线交点到破片初始位置X0的距离。
将h(t)化为:
式(7)中,皆可以由弹目交会参数求出。若存在某时刻t*使得s(t)=h(t),则说明破片可以命中目标面元所在的平面,即要求解方程:
将该方程的解代入式(5)或式(1)求得命中点X*的坐标:
通过破片速度方向向量、目标速度和面元方向向量求得落角:
步骤4,判定命中点是否在三角形面元内,若位于三角形面元内,则破片命中该面元;具体为:
用三角形面积法判断该命中点是否属于该目标面元,即判断某点与三角形任意两顶点面积之和是否等于该三角形的面积,如果等于则该点在三角形内,否则在三角形外;
假设三角形面元ΔA′B′C′,T′为破片与三角形面元平面的交点,判定下式:
S(ΔT′A′B′)+S(ΔT′B′C′)+S(ΔT′C′A′)≥S(ΔA′B′C′) (11)
若取等号则T′位于三角形面元内,破片命中该面元;并且以上获得的命中点参数即为该破片的命中点参数;否则,点T′位于面元外。
下面结合具体实施例对本发明作进一步说明。
实施例
如图1所示,基于引战配合可视化仿真软件的一种弹目交会破片命中点参数计算模型方法,具体步骤如下:
步骤1,将目标模型数据以xml格式导入引战配合可视化仿真软件,输入战斗部参数,交会参数,获得破片总数Ai,部位总数Aj,面元数Ajk。
步骤2,第i枚破片对j部位的所有面元进行一次遍历。
步骤3,用牛顿迭代法求解s(t)=h(t),若无解则直接判断为未命中,若有解则计算命中点参数,并将命中点参数保存于软件数据库的列表中。
步骤4,一次遍历完成后,j移向下一部位,直到完成第i枚破片对所有部位的遍历。
步骤5,重复步骤2、3、4,直到所有破片完成遍历,结束算法。
本发明通过设置不同衰减系数α,不同目标速度vt,同一交会参数的初始条件下仿真出了命中点参数的结算统计结果。结果表明,本算法对于尺寸较小或者高速运动的目标,更加符合实际情况,计算更加准确,为毁伤概率的计算提供了更加可靠地参数。
Claims (5)
1.一种引战配合仿真交会破片命中点参数计算方法,其特征在于,包括以下步骤:
步骤1,对目标进行面元划分,得到目标三角面元模型;
步骤2,建立破片与目标面元交会的数学模型;
步骤3,计算目标面元所在平面与破片射线交点与命中点参数;
步骤4,判定命中点是否在三角形面元内,若位于三角形面元内,则破片命中该面元。
2.根据权利要求1所述的引战配合仿真交会破片命中点参数计算方法,其特征在于,步骤1具体为:
对目标进行面元划分,得到目标三角面元模型;提取出面元的3个顶点坐标,由A、B、C三点确定面元的法向量n′=sAB×sAC,并将n′单位化得到n;其中sAB、sAC分别为AB连线、AC连线的方向矢量。
3.根据权利要求1所述的引战配合仿真交会破片命中点参数计算方法,其特征在于,步骤2具体为:
以战斗部爆炸时刻为0时刻,列出破片的轨迹方程:
式(1)中,X为t时刻破片空间坐标,s(t)为破片在t时刻所飞行的距离,为破片速度方向的单位矢量,X0为0时刻破片的空间坐标;
考虑破片在空气中飞行速度的衰减,其式为:
vf=v0exp(-αs) (2)
vf为破片在t时刻的速度,v0为破片初速速率,α为空气密度衰减系数,s为破片在t时刻所飞行的距离;
由(1)、(2)式得到破片飞行距离与时间的关系:
S(t)=[ln(1+αv0t)]/α (3)
在0时刻,目标某面元所在的平面方程表示为:
n·X+D0=0
D0=-n·Y0 (4)
式(4)中,n为目标面元法向单位向量,D0为一标量常数,表示坐标原点到目标面元的距离;Y0为三角形面元中任意一顶点的初始坐标;
由该面元在0时刻的3个顶点坐标即可求得n和D0的值;
假设面元只有平动,面元的法向不变,面元所在平面的运动方程表示为:
n·(X-vtt)+D0=0 (5)
式(5)中,vt为目标在遭遇段的运动速度。
4.根据权利要求1所述的引战配合仿真交会破片命中点参数计算方法,其特征在于,步骤3具体为:
由式(1)和式(4)联立解得:
式(6)中,h(t)为t时刻面元与破片射线交点到破片初始位置X0的距离;
将h(t)化为:
式(7)中,若存在某时刻t*使得s(t)=h(t),则说明破片可以命中目标面元所在的平面,即要求解方程:
将该方程的解代入式(5)或式(1)求得命中点X*的坐标:
通过破片速度方向向量、目标速度和面元方向向量求得落角:
5.根据权利要求1所述的引战配合仿真交会破片命中点参数计算方法,其特征在于,步骤4具体为:
用三角形面积法判断该命中点是否属于该目标面元,即判断某点与三角形任意两顶点面积之和是否等于该三角形的面积,如果等于则该点在三角形内,否则在三角形外;
假设三角形面元ΔA′B′C′,T′为破片与三角形面元平面的交点,判定下式:
S(ΔT′A′B′)+S(ΔT′B′C′)+S(ΔT′C′A′)≥S(ΔA′B′C′) (11)
若取等号则T′位于三角形面元内,破片命中该面元;并且以上获得的命中点参数即为该破片的命中点参数;否则,点T′位于面元外。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710410853.6A CN108984808A (zh) | 2017-06-04 | 2017-06-04 | 引战配合仿真交会破片命中点参数计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710410853.6A CN108984808A (zh) | 2017-06-04 | 2017-06-04 | 引战配合仿真交会破片命中点参数计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108984808A true CN108984808A (zh) | 2018-12-11 |
Family
ID=64501745
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710410853.6A Withdrawn CN108984808A (zh) | 2017-06-04 | 2017-06-04 | 引战配合仿真交会破片命中点参数计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108984808A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110991030A (zh) * | 2019-11-28 | 2020-04-10 | 北京电子工程总体研究所 | 一种通用复杂引战系统仿真的方法及系统 |
CN111125623A (zh) * | 2019-11-28 | 2020-05-08 | 北京电子工程总体研究所 | 一种超低空作战条件下的战斗部破片飞散轨迹计算的方法 |
CN112861330A (zh) * | 2021-01-25 | 2021-05-28 | 清华大学 | 一种基于matlab的导弹杀伤效果计算和可视化方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040068415A1 (en) * | 2002-04-22 | 2004-04-08 | Neal Solomon | System, methods and apparatus for coordination of and targeting for mobile robotic vehicles |
CN1779409A (zh) * | 2005-09-20 | 2006-05-31 | 湖南文理学院 | 近炸引信预制破片弹对反舰导弹毁伤效应的计算方法 |
-
2017
- 2017-06-04 CN CN201710410853.6A patent/CN108984808A/zh not_active Withdrawn
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040068415A1 (en) * | 2002-04-22 | 2004-04-08 | Neal Solomon | System, methods and apparatus for coordination of and targeting for mobile robotic vehicles |
CN1779409A (zh) * | 2005-09-20 | 2006-05-31 | 湖南文理学院 | 近炸引信预制破片弹对反舰导弹毁伤效应的计算方法 |
Non-Patent Citations (1)
Title |
---|
王马法 等: "弹目交会中破片命中点参数计算模型研究", 《弹道学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110991030A (zh) * | 2019-11-28 | 2020-04-10 | 北京电子工程总体研究所 | 一种通用复杂引战系统仿真的方法及系统 |
CN111125623A (zh) * | 2019-11-28 | 2020-05-08 | 北京电子工程总体研究所 | 一种超低空作战条件下的战斗部破片飞散轨迹计算的方法 |
CN111125623B (zh) * | 2019-11-28 | 2023-09-19 | 北京电子工程总体研究所 | 一种超低空作战条件下的战斗部破片飞散轨迹计算的方法 |
CN110991030B (zh) * | 2019-11-28 | 2023-10-13 | 北京电子工程总体研究所 | 一种通用复杂引战系统仿真的方法及系统 |
CN112861330A (zh) * | 2021-01-25 | 2021-05-28 | 清华大学 | 一种基于matlab的导弹杀伤效果计算和可视化方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104215935B (zh) | 一种基于决策加权融合的雷达炮弹目标识别方法 | |
CN108984808A (zh) | 引战配合仿真交会破片命中点参数计算方法 | |
CN104978693B (zh) | 一种基于国网gis电力应急处置态势模拟推演系统及方法 | |
CN110340890A (zh) | 一种空间机械臂全局无碰轨迹规划系统 | |
US20170278303A1 (en) | Method and system for virtual sensor data generation with depth ground truth annotation | |
CN103631153B (zh) | 可视化双模制导航空时敏炸弹半实物仿真方法 | |
KR102272411B1 (ko) | 합성 개구면 레이더 영상의 시뮬레이션 영상 데이터베이스의 식별률을 높이기 위한 인공신경망 학습 방법 및 그 장치 | |
CN105259913B (zh) | 确定飞行器自动着陆引导指令的方法及装置 | |
RU2571567C2 (ru) | Способ управления движением аэробаллистического летательного аппарата по заданной пространственной траектории | |
CN103699715A (zh) | 一种基于光滑粒子流体动力学和非线性有限元的流固耦合方法 | |
Yang et al. | Unity 3D production and environmental perception vehicle simulation platform | |
CN109325249A (zh) | 用于在虚拟测试环境中对机动车辆传感器建模的方法 | |
US20190381355A1 (en) | Sport range simulator | |
CN107152887A (zh) | 一种基于虚拟现实技术的模拟火箭筒操作训练方法和系统 | |
CN105787155A (zh) | 红外/毫米波制导的自主寻的导弹弹道分布式联合仿真框架 | |
US9646417B1 (en) | Augmented reality system for field training | |
CN105334753A (zh) | 一种近地告警视景仿真验证平台系统及其验证方法 | |
RU2432592C1 (ru) | Моделирующий комплекс для проверки системы управления беспилотного летательного аппарата | |
US10319132B2 (en) | Method and system for representing objects with velocity-dependent particles | |
CN110472185B (zh) | 圆形弹丸击打与地面呈倾斜角的窗户的弹道轨迹计算方法 | |
CN110472180B (zh) | 圆形弹丸击打与地面呈垂直的窗户的弹道轨迹计算方法 | |
CN105259789B (zh) | 一种探测仿真装置 | |
Schüβler et al. | Achieving efficient and realistic full-radar simulations and automatic data annotation by exploiting ray meta data from a radar ray tracing simulator | |
KR102578886B1 (ko) | 발사체 식별 방법, 발사체 식별 장치 및 상기 방법을 실행시키기 위하여 기록매체에 저장된 컴퓨터 프로그램 | |
KR101926107B1 (ko) | 항공 무장효과 분석 시뮬레이터 및 그 시뮬레이션 방법 |
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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20181211 |