CN116579118A - 面向近距空战格斗的三维动态发射区在线快速解算方法 - Google Patents
面向近距空战格斗的三维动态发射区在线快速解算方法 Download PDFInfo
- Publication number
- CN116579118A CN116579118A CN202310067984.4A CN202310067984A CN116579118A CN 116579118 A CN116579118 A CN 116579118A CN 202310067984 A CN202310067984 A CN 202310067984A CN 116579118 A CN116579118 A CN 116579118A
- Authority
- CN
- China
- Prior art keywords
- missile
- target
- model
- boundary
- maneuver
- 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.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims description 15
- 238000000034 method Methods 0.000 claims abstract description 52
- 239000013598 vector Substances 0.000 claims description 47
- 230000008859 change Effects 0.000 claims description 23
- 230000009194 climbing Effects 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 12
- 238000005474 detonation Methods 0.000 claims description 8
- 239000000446 fuel Substances 0.000 claims description 8
- 238000010297 mechanical methods and process Methods 0.000 claims description 8
- 230000005226 mechanical processes and functions Effects 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 6
- 230000009189 diving Effects 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 4
- 230000006870 function Effects 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 230000003068 static effect Effects 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 238000010276 construction Methods 0.000 claims description 2
- 238000001816 cooling Methods 0.000 claims description 2
- 230000001419 dependent effect Effects 0.000 claims description 2
- 229920001971 elastomer Polymers 0.000 claims description 2
- 239000000806 elastomer Substances 0.000 claims description 2
- 230000008685 targeting Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 5
- 230000007423 decrease Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
本发明提供了一种面向近距空战格斗的三维动态发射区在线快速解算方法,通过构建导弹模型和目标机动模型,根据当前近距空战态势信息,结合导弹模型和目标机动模型,构建导弹可攻击边界解算模型;解算目标不同机动模式下的远边界、近边界和不可逃逸边界,确定近距空空导弹三维动态发射区。本发明直观地呈现了当前激烈近距空战对抗环境下导弹的三维形式的可攻击范围,与此同时依托武器装备参数建立模型并进行实时在线解算,保证可攻击边界的解算精度。采用三维动态发射区在线快速解算方法能够在强动态、高复杂的近距空战格斗环境下帮助己方人员准确、实时地掌握当前敌我态势并为后续近距空战决策提供数据支撑。
Description
技术领域
本发明涉及航空武器火力控制领域,具体为一种导弹三维动态发射区的解算方法。
背景技术
伴随着信息化战争的发展和变革,近距空战已发展为以近距空空导弹为主要攻击武器而进行的多回合、强对抗、高动态的近距格斗过程。在此过程中,近距空空导弹的使用成为飞行员掌控战场主动权,夺取战斗胜利的关键。导弹攻击区作为反映当前导弹性能和对敌攻击态势的重要指标,能够为导弹的相关决策提供有力的支撑和参考。
导弹攻击区的呈现形式一般可分为两种:(1)以目标机为中心;(2)以本机为中心(又称为动态发射区)。目前,导弹攻击区主要采用前者的形式予以呈现,并且多针对二维情况。该呈现形式下的攻击区样式难以为飞行员带来直观的对抗态势信息,同时存在应用局限性,当面向剧烈变化的近距格斗场景时难以满足快速决策的需求,同时影响决策质量。传统的导弹攻击区解算方法主要是多项式拟合法与查表插值法。其中,多项式拟合法通过构造多项式实现攻击区边界的拟合,具有解算速度快、易于实现的优点,但由于多采用固定模型函数回归分析的方法实现拟合,其拟合精度较差。查表插值法则是将影响攻击区的关键参数与对应边界距离制表,利用相关数据进行线性插值近似,该方法解算实时性较好,但面向三维空间的数据量激增,火控系统的内存占用问题十分突出。
发明内容
为了克服现有技术的不足,本发明提供一种面向近距空战格斗的三维动态发射区在线快速解算方法,以解决现有技术方法在近距空战格斗过程中难以准确掌握对敌实时攻击态势、导弹攻击区二维形式的实际应用存在较强局限性、导弹攻击区的解算精度不足等问题。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤一:针对近距空空导弹的制导打击过程进行建模,构建导弹模型;
(1)建立导弹发动机模型;
(2)建立导弹运动模型;
(3)建立导弹制导模型;
(4)建立导弹导引头模型;
(5)建立导弹引战配合模型;
步骤二:针对目标遭受导弹攻击情况下可能执行的机动逃逸策略,构建目标机动模型;
建立目标匀速飞行机动模型;
建立目标最大过载置尾机动模型;
建立目标最大过载左转/右转机动模型;
建立目标最大过载爬升/俯冲机动模型;
步骤三:根据当前近距空战态势信息,结合导弹模型和目标机动模型,构建导弹可攻击边界解算模型;
步骤3.1:建立可攻击边界搜索区域模型;
根据当前己方和目标的相对位置,以及己方的飞行速度方向,确定导弹动态发射区的解算搜索区域方程,具体为:
其中,表示动态发射区的解算搜索平面的法向量;/>表示己方飞机三轴速度所构成的速度矢量;/>表示己方和目标间的相对位置矢量;/>表示辅助向量,对应值为/>也表示辅助向量,对应值为/>
步骤3.2:建立可攻击边界搜索模型;
步骤四:解算目标不同机动模式下的远边界、近边界和不可逃逸边界,确定近距空空导弹三维动态发射区;
目标的6种机动动作为匀速飞行机动、最大过载置尾机动、最大过载左转机动、最大过载右转机动、最大过载爬升机动和最大过载爬升机动,分别记作MVR-1、MVR-2、MVR-3、MVR-4、MVR-5和MVR-6,按照步骤三中建立的目标机动模型,依次解算目标采用6种机动动作时,不同离轴发射角所对应的最远可攻击边界和最近可攻击边界,得到远边界矩阵DisL和近边界矩阵DisN:
其中,符号表示向下取整;/>表示导弹针对目标执行机动MVR-m时对应离轴发射角/>的可攻击远边界,m=1,2,…,6,表示采取不同机动动作所对应的动作编号,/>表征离轴发射角的遍历参数;/>表示导弹针对目标执行机动MVR-m时对应离轴发射角/>的可攻击近边界,确定导弹动态发射区的远边界、近边界和不可逃逸边界时,采用目标保持当前状态进行匀速飞行机动的可攻击远边界计算值、可攻击近边界计算值作为远边界值disFar和近边界值disNear;对于不可逃逸边界值,则根据公式(35)进行确定:
根据边界类型,按照离轴发射角遍历顺序依次连接对应边界距离的端点,形成导弹可攻击远边界、不可逃逸边界和近边界,最终得到当前态势情况下的近距空空导弹动态发射区。
所述步骤一中,构建导弹模型的具体步骤为:
(1)建立导弹发动机模型
导弹发动机模型包括导弹发动机推力模型和弹体重量变化模型,根据导弹发射飞行过程中导弹发动机推力与导弹飞行时间的变化关系,构建导弹发动机推力模型,发动机推力Pm具体为:
其中,P1和P2分别表示两个工作阶段发动机提供的动推力;Se表示发动机喷管出口处的横截面积;P0表示海平面大气静压强;PH表示导弹所处高度的大气静压强;tf表示导弹飞行时间;t1和t3分别为发动机两阶段的工作时间;t2表示两阶段切换过程中发动机的冷却时间;
与此同时,近距空空导弹飞行过程伴随着燃料损耗,整体重量随导弹飞行时间tf变化;因此,构建弹体重量变化模型,具体为:
其中,mweight表示导弹实时重量;m0表示导弹空重;mFuel0表示导弹发动机携带的燃料的初始重量;rFuel1和rFuel2分别表示第一阶段和第二阶段下的发动机燃料消耗速度;若导弹发动机采用单推力策略导弹发动机推力模型,则公式(2)中发动机燃料消耗速度rFuel1和rFuel2相等;
(2)建立导弹运动模型
导弹运动模型由反映弹体姿态变化的动力学模型和导弹运动的运动学模型组成,导弹动力学方程的构建考虑导弹重力Gm、空气阻力Qm、发动机推力Pm的影响;
导弹重力Gm的计算模型为:
Gm=mweightg (7)
其中,g表示重力加速度,取值9.8m/s2;
导弹遭受的空气阻力Qm的计算模型为:
其中,为导弹阻力系数,与导弹飞行马赫数Ma有关;Sm为导弹最大截面积,ρH为导弹所处高度H的大气密度;Vm为导弹的当前飞行速度;
结合导弹制导系统提供的导弹在俯仰方向和偏航方向上的需用过载nmy和nmz,以及导弹的飞行俯仰角θm和飞行方位角φm,构建近距空空导弹的动力学模型,具体为:
其中,符号·表示该变量的导数;
近距空空导弹的运动学方程为:
(3)建立导弹制导模型;
导弹制导模型包括导引律控制模型和导弹飞行姿态调整模型,根据导弹和目标的位置、速度信息以及相对运动参数,在俯仰方向和偏航方向上建立导引律控制模型,具体为:
其中,K1、K2分别表示俯仰方向和偏航方向上的导引比例系数,R表示导弹和目标间的相对距离,θL、φL分别表示目标线的倾角和偏角,(xm,ym,zm)为当前时刻导弹的位置坐标,(xt,yt,zt)为当前时刻目标的位置坐标;
依据导引律控制模型计算俯仰和偏航方向上的需用过载nmy和nmz,再结合公式(5),构建导弹飞行姿态调整模型为:
其中,θm_t、φm_t为当前时刻导弹的飞行俯仰角和飞行方位角;θm_t+△t、φm_t+△t为间隔时间△t后的导弹飞行俯仰角和飞行方位角更新值;
(4)建立导弹导引头模型
导弹导引头模型根据距离截获状态FlagD、角度截获状态FlagA、跟踪角速度约束情况Flagω和导弹工作时间约束情况FlagT,综合判断能否对目标保持截获跟踪状态;
导引头的距离截获条件表示为:
其中,Dmax表示导引头的最大可探测距离;数字1表示导引头对目标成功实现距离截获;数字0表示导引头对目标距离截获失败;最大可探测距离与当前目标进入角αentry大小相关,最大可探测距离随目标进入角αentry的递增先增大后减小,两者间的数值关系根据导引头具体性能数据,采用如下多项式形式进行拟合表征:
其中,Ci表示常数,i=1,2,…,n;表示目标进入角αentry的n次方;
导引头的角度截获条件表示为:
其中,AngleIR表示导引头的最大可探测角域;Angle表示目标方位角,根据弹目连线向量与导弹飞行速度向量/>进行计算,具体公式为/>数字1表示导引头对目标成功实现角度截获;数字0表示导引头对目标角度截获失败;
导引头的跟踪角速度约束条件表示为:
其中,ωmax表示导弹导引头对目标的最大跟踪角速度;ωR为目标线的旋转角速度,计算公式为其中/>为导弹速度向量,/>为目标速度向量,/>为弹目连线向量,Angle为目标方位角,数字1表示导引头跟踪角速度满足需求,能够成功跟踪目标;数字0表示由于跟踪角速度大小限制,导引头无法跟踪目标;
导引头的工作时间约束条件表示为:
其中,Tmax表示导弹可飞行的最长时间;数字1表示导引头处于正常工作时间内;数字0表示导引头超出正常工作时间,无法对目标进行截获;
综上,导弹对目标的截获跟踪状态Flagcapture表示为:
Flagcapture=FlagD∧FlagA∧Flagω∧FlagT (18)
其中,符号∧为逻辑与运算符,仅当FlagD、FlagA、Flagω和FlagT均为1时,导弹才能够保持对目标的截获跟踪,否则将会丢失目标;
(5)建立导弹引战配合模型
导弹引战配合模型根据弹目间的相对距离条件ConsD、相对速度条件ConsV和保险解除时间条件ConsT,判断导弹能够正常引爆并对目标造成杀伤;
弹目间的相对距离条件表示为:
其中,Reffective表示导弹的有效杀伤半径;数字1表示当前弹目距离在导弹有效杀伤半径内,能够对目标进行有效杀伤;数字0表示当前弹目距离不满足杀伤条件;
相对速度条件ConsV表示为:
其中,表示导弹当前的飞行速度向量;/>表示目标当前的飞行速度向量;Vfuse表示导弹引信起爆的最小遇靶相对速度条件阈值;数字1表示当前弹目相对速度满足引信起爆条件;数字0表示当前弹目相对速度不满足引信起爆条件;
保险解除时间条件ConsT表示为:
其中,TIns表示引信保险工作的解锁时间;数字1表示导弹飞行时长大于引信保险解锁时间,引信开始正常工作;数字0表示导弹飞行时长较短,引信还未解锁;
综上,成功触发导弹引信并引爆战斗部的条件ConsHit表示为:
ConsHit=ConsD∧ConsV∧ConsT (22)
即,仅当ConsD、ConsV和ConsT均为1时才能完成导弹引爆并杀伤目标。
所述步骤二中构建目标机动模型的具体步骤为:
(1)建立目标匀速飞行机动模型
根据目标当前飞行速度为Vt、飞行俯仰角为θt和飞行方位角为φt,其匀速飞行机动模型为:
其中,Vt+Δt、θt+Δt、φt+Δt分别表示时间间隔Δt后的飞行速度、飞行俯仰角和飞行方位角的更新值;
(2)建立目标最大过载置尾机动模型
目标以最大过载noverload,在目标线矢量与自身速度矢量/>所形成的平面内进行置尾机动,直至目标速度矢量与初始目标线矢量保持方向一致;置尾机动过程中,目标飞行速度的变化为:
其中,at为目标最大加速度;Vt_max为目标最大飞行速度;
置尾机动的旋转轴线向量为:
其中,cross()表示求解两向量叉积;||表示求解向量长度。
置尾机动过程中待旋转的角度aimAnglezw为:
目标在置尾机动过程中的飞行角速度ωzw为:
(3)建立目标最大过载左转/右转机动模型;
目标在当前状态下,以最大过载noverload进行左转或右转机动,直至自身速度矢量与初始目标线矢量/>形成垂直条件,随后保持该方向飞行;左转或右转机动过程中,目标飞行速度的变化按照公式(20)进行处理;
目标在左转机动过程中,转弯角速度ωL为:
目标在右转机动过程中,转弯角速度ωR为:
(4)建立目标最大过载爬升/俯冲机动模型;
目标在当前状态下,以最大过载noverload进行爬升机动,直至最终的飞行方向与初始目标线矢量和自身速度矢量/>所形成平面的法向量方向一致,或以最大过载noverload进行俯冲机动,直至最终的飞行方向与初始目标线矢量/>和自身速度矢量/>所形成平面的法向量方向相反,随后保持飞行状态;爬升/俯冲机动过程中,目标飞行速度的变化按照公式(20)进行处理;
目标在爬升机动过程中,最终的飞行速度的参考指向为:
目标爬升机动过程中的旋转轴线向量为:
爬升过程中的角度变化为:
目标在俯冲机动过程中,最终的飞行速度的参考指向为:
目标俯冲机动过程中的旋转轴线向量为:
俯冲过程中的角度变化为:
所述步骤3.2建立可攻击边界搜索模型的步骤如下:
步骤3.2.1:根据当前近距空战对抗态势,给定初始预估搜索边界的上界值和下界值为[a0,b0];
步骤3.2.2:在公式(32)确定的搜索区域内,根据导弹离轴发射角范围给定初始搜索角度值/>
步骤3.2.3:在给定离轴发射角上,基于当前敌我位置、速度信息,采用建立的导弹模型和目标机动模型,进行导弹制导打击仿真,依据导弹引战配合模型判定导弹最终能否命中杀伤目标,从而对搜索边界的上界值、下界值[ai,bi]进行更新,具体更新方式为:
其中,边界搜索点gp=ai+0.618(bi-ai);
步骤3.2.4:重复步骤3.2.3的边界搜索过程,直至可攻击边界的精度ε满足|bi-ai|≤ε,则以作为当前离轴发射角上的可攻击边界值;
步骤3.2.5:以Δτ为间隔进行离轴发射角更新,然后重复步骤3.2.3和3.2.4,直至完成离轴发射角遍历。
本发明的有益效果在于以己方飞机为中心,直观地呈现了当前激烈近距空战对抗环境下导弹的三维形式的可攻击范围,与此同时依托武器装备参数建立模型并进行实时在线解算,保证可攻击边界的解算精度。采用三维动态发射区在线快速解算方法能够在强动态、高复杂的近距空战格斗环境下帮助己方人员准确、实时地掌握当前敌我态势并为后续近距空战决策提供数据支撑。
附图说明
图1是本发明目标最大过载置尾机动示意图。
图2是本发明目标最大过载左转机动示意图。
图3是本发明目标最大过载右转机动示意图。
图4是本发明目标最大过载爬升机动示意图。
图5是本发明目标最大过载俯冲机动示意图。
图6是本发明导弹动态发射区解算结果。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
本发明的实例包括以下几个步骤:
(1)以己方飞机为原点,建立北天东坐标系求解动态发射区,设定当前近距空战对抗态势信息为:
(2)以某型近距空空导弹为参考,构建近距空空导弹模型;
导弹发动机推力的变化情况为:
飞行过程中,导弹重量变化遵循如下规律:
空气阻力Qm中的各项参量分别为:
Sm=0.02
导弹制导模型中的导引比例系数设置为K1=3、K2=4。
导弹导引头模型中的各项参数设置为:
其中,红外导引头最大探测距离Dmax采用关于目标进入角αentry的多项式拟合函数进行表示,具体为:
导弹引战配合模型中的各项参数设置为:
(3)基于典型战机参数,构建目标机动模型
目标各项机动参数设置为:
依照机动参数设置,仿真目标采用最大过载置尾机动、最大过载左转/右转机动和最大过载爬升/俯冲机动的飞行轨迹,分别如图1、图2、图3、图4和图5所示。
(4)解算动态发射区远边界、不可逃逸边界和近边界
根据当前近距空战对抗态势,设定远边界初始预估搜索边界为[2000m,20000m],远边界搜索精度为ε=50m,近边界初始预估搜索边界为[200m,2000m],近边界搜索精度为ε=20m,导弹离轴发射角搜索间隔Δτ=6°。依照动态发射区边界搜索流程,依次解算目标采用6种不同机动情况下的导弹可攻击边界。按照远边界、不可逃逸边界和近边界距离定义,得到导弹不同离轴发射角对应的各项边界值,如下所示。动态发射区边界搜索总共耗时165.38ms,满足近距空战格斗环境下的使用需求。
(5)绘制动态发射区进行直观呈现
根据当前态势信息,解算得到动态发射区所在区域平面的法向量为据此确定动态发射区样式,并同时结合解算得到的远边界、不可逃逸边界和近边界数据,得到最终三维动态发射区结果,如图6所示,可观察到当前目标处于己方飞机导弹武器的攻击范围之外,暂无法对其进行有效打击。/>
Claims (4)
1.一种面向近距空战格斗的三维动态发射区在线快速解算方法,其特征在于包括下述步骤:
步骤一:针对近距空空导弹的制导打击过程进行建模,构建导弹模型;
(1)建立导弹发动机模型;
(2)建立导弹运动模型;
(3)建立导弹制导模型;
(4)建立导弹导引头模型;
(5)建立导弹引战配合模型;
步骤二:针对目标遭受导弹攻击情况下可能执行的机动逃逸策略,构建目标机动模型;
建立目标匀速飞行机动模型;
建立目标最大过载置尾机动模型;
建立目标最大过载左转/右转机动模型;
建立目标最大过载爬升/俯冲机动模型;
步骤三:根据当前近距空战态势信息,结合导弹模型和目标机动模型,构建导弹可攻击边界解算模型;
步骤3.1:建立可攻击边界搜索区域模型;
根据当前己方和目标的相对位置,以及己方的飞行速度方向,确定导弹动态发射区的解算搜索区域方程,具体为:
其中,表示动态发射区的解算搜索平面的法向量;/>表示己方飞机三轴速度所构成的速度矢量;/>表示己方和目标间的相对位置矢量;/>表示辅助向量,对应值为 也表示辅助向量,对应值为/>
步骤3.2:建立可攻击边界搜索模型;
步骤四:解算目标不同机动模式下的远边界、近边界和不可逃逸边界,确定近距空空导弹三维动态发射区;
目标的6种机动动作为匀速飞行机动、最大过载置尾机动、最大过载左转机动、最大过载右转机动、最大过载爬升机动和最大过载爬升机动,分别记作MVR-1、MVR-2、MVR-3、MVR-4、MVR-5和MVR-6,按照步骤三中建立的目标机动模型,依次解算目标采用6种机动动作时,不同离轴发射角所对应的最远可攻击边界和最近可攻击边界,得到远边界矩阵DisL和近边界矩阵DisN:
其中,符号表示向下取整;/>表示导弹针对目标执行机动MVR-m时对应离轴发射角的可攻击远边界,m=1,2,…,6,表示采取不同机动动作所对应的动作编号,/>表征离轴发射角的遍历参数;/>表示导弹针对目标执行机动MVR-m时对应离轴发射角/>的可攻击近边界,确定导弹动态发射区的远边界、近边界和不可逃逸边界时,采用目标保持当前状态进行匀速飞行机动的可攻击远边界计算值、可攻击近边界计算值作为远边界值disFar和近边界值disNear;对于不可逃逸边界值,则根据公式(35)进行确定:
根据边界类型,按照离轴发射角遍历顺序依次连接对应边界距离的端点,形成导弹可攻击远边界、不可逃逸边界和近边界,最终得到当前态势情况下的近距空空导弹动态发射区。
2.根据权利要求1所述的面向近距空战格斗的三维动态发射区在线快速解算方法,其特征在于:
所述步骤一中,构建导弹模型的具体步骤为:
(1)建立导弹发动机模型
导弹发动机模型包括导弹发动机推力模型和弹体重量变化模型,根据导弹发射飞行过程中导弹发动机推力与导弹飞行时间的变化关系,构建导弹发动机推力模型,发动机推力Pm具体为:
其中,P1和P2分别表示两个工作阶段发动机提供的动推力;Se表示发动机喷管出口处的横截面积;P0表示海平面大气静压强;PH表示导弹所处高度的大气静压强;tf表示导弹飞行时间;t1和t3分别为发动机两阶段的工作时间;t2表示两阶段切换过程中发动机的冷却时间;
与此同时,近距空空导弹飞行过程伴随着燃料损耗,整体重量随导弹飞行时间tf变化;因此,构建弹体重量变化模型,具体为:
其中,mweight表示导弹实时重量;m0表示导弹空重;mFuel0表示导弹发动机携带的燃料的初始重量;rFuel1和rFuel2分别表示第一阶段和第二阶段下的发动机燃料消耗速度;若导弹发动机采用单推力策略导弹发动机推力模型,则公式(2)中发动机燃料消耗速度rFuel1和rFuel2相等;
(2)建立导弹运动模型
导弹运动模型由反映弹体姿态变化的动力学模型和导弹运动的运动学模型组成,导弹动力学方程的构建考虑导弹重力Gm、空气阻力Qm、发动机推力Pm的影响;
导弹重力Gm的计算模型为:
Gm=mweightg (7)
其中,g表示重力加速度,取值9.8m/s2;
导弹遭受的空气阻力Qm的计算模型为:
其中,为导弹阻力系数,与导弹飞行马赫数Ma有关;Sm为导弹最大截面积,ρH为导弹所处高度H的大气密度;Vm为导弹的当前飞行速度;
结合导弹制导系统提供的导弹在俯仰方向和偏航方向上的需用过载nmy和nmz,以及导弹的飞行俯仰角θm和飞行方位角φm,构建近距空空导弹的动力学模型,具体为:
其中,符号·表示该变量的导数;
近距空空导弹的运动学方程为:
(3)建立导弹制导模型;
导弹制导模型包括导引律控制模型和导弹飞行姿态调整模型,根据导弹和目标的位置、速度信息以及相对运动参数,在俯仰方向和偏航方向上建立导引律控制模型,具体为:
其中,K1、K2分别表示俯仰方向和偏航方向上的导引比例系数,R表示导弹和目标间的相对距离,θL、φL分别表示目标线的倾角和偏角,(xm,ym,zm)为当前时刻导弹的位置坐标,(xt,yt,zt)为当前时刻目标的位置坐标;
依据导引律控制模型计算俯仰和偏航方向上的需用过载nmy和nmz,再结合公式(5),构建导弹飞行姿态调整模型为:
其中,θm_t、φm_t为当前时刻导弹的飞行俯仰角和飞行方位角;θm_t+△t、φm_t+△t为间隔时间△t后的导弹飞行俯仰角和飞行方位角更新值;
(4)建立导弹导引头模型
导弹导引头模型根据距离截获状态FlagD、角度截获状态FlagA、跟踪角速度约束情况Flagω和导弹工作时间约束情况FlagT,综合判断能否对目标保持截获跟踪状态;
导引头的距离截获条件表示为:
其中,Dmax表示导引头的最大可探测距离;数字1表示导引头对目标成功实现距离截获;数字0表示导引头对目标距离截获失败;最大可探测距离与当前目标进入角αentry大小相关,最大可探测距离随目标进入角αentry的递增先增大后减小,两者间的数值关系根据导引头具体性能数据,采用如下多项式形式进行拟合表征:
其中,Ci表示常数,i=1,2,…,n;表示目标进入角αentry的n次方;
导引头的角度截获条件表示为:
其中,AngleIR表示导引头的最大可探测角域;Angle表示目标方位角,根据弹目连线向量与导弹飞行速度向量/>进行计算,具体公式为/>数字1表示导引头对目标成功实现角度截获;数字0表示导引头对目标角度截获失败;
导引头的跟踪角速度约束条件表示为:
其中,ωmax表示导弹导引头对目标的最大跟踪角速度;ωR为目标线的旋转角速度,计算公式为其中/>为导弹速度向量,/>为目标速度向量,/>为弹目连线向量,Angle为目标方位角,数字1表示导引头跟踪角速度满足需求,能够成功跟踪目标;数字0表示由于跟踪角速度大小限制,导引头无法跟踪目标;
导引头的工作时间约束条件表示为:
其中,Tmax表示导弹可飞行的最长时间;数字1表示导引头处于正常工作时间内;数字0表示导引头超出正常工作时间,无法对目标进行截获;
综上,导弹对目标的截获跟踪状态Flagcapture表示为:
Flagcapture=FlagD∧FlagA∧Flagω∧FlagT (18)
其中,符号∧为逻辑与运算符,仅当FlagD、FlagA、Flagω和FlagT均为1时,导弹才能够保持对目标的截获跟踪,否则将会丢失目标;
(5)建立导弹引战配合模型
导弹引战配合模型根据弹目间的相对距离条件ConsD、相对速度条件ConsV和保险解除时间条件ConsT,判断导弹能够正常引爆并对目标造成杀伤;
弹目间的相对距离条件表示为:
其中,Reffective表示导弹的有效杀伤半径;数字1表示当前弹目距离在导弹有效杀伤半径内,能够对目标进行有效杀伤;数字0表示当前弹目距离不满足杀伤条件;
相对速度条件ConsV表示为:
其中,表示导弹当前的飞行速度向量;/>表示目标当前的飞行速度向量;Vfuse表示导弹引信起爆的最小遇靶相对速度条件阈值;数字1表示当前弹目相对速度满足引信起爆条件;数字0表示当前弹目相对速度不满足引信起爆条件;
保险解除时间条件ConsT表示为:
其中,TIns表示引信保险工作的解锁时间;数字1表示导弹飞行时长大于引信保险解锁时间,引信开始正常工作;数字0表示导弹飞行时长较短,引信还未解锁;
综上,成功触发导弹引信并引爆战斗部的条件ConsHit表示为:
ConsHit=ConsD∧ConsV∧ConsT (22)
即,仅当ConsD、ConsV和ConsT均为1时才能完成导弹引爆并杀伤目标。
3.根据权利要求1所述的面向近距空战格斗的三维动态发射区在线快速解算方法,其特征在于:
所述步骤二中构建目标机动模型的具体步骤为:
(1)建立目标匀速飞行机动模型
根据目标当前飞行速度为Vt、飞行俯仰角为θt和飞行方位角为φt,其匀速飞行机动模型为:
其中,Vt+Δt、θt+Δt、φt+Δt分别表示时间间隔Δt后的飞行速度、飞行俯仰角和飞行方位角的更新值;
(2)建立目标最大过载置尾机动模型
目标以最大过载noverload,在目标线矢量与自身速度矢量/>所形成的平面内进行置尾机动,直至目标速度矢量与初始目标线矢量保持方向一致;置尾机动过程中,目标飞行速度的变化为:
其中,at为目标最大加速度;Vt_max为目标最大飞行速度;
置尾机动的旋转轴线向量为:
其中,cross()表示求解两向量叉积;||表示求解向量长度;
置尾机动过程中待旋转的角度aimAnglezw为:
目标在置尾机动过程中的飞行角速度ωzw为:
(3)建立目标最大过载左转/右转机动模型;
目标在当前状态下,以最大过载noverload进行左转或右转机动,直至自身速度矢量与初始目标线矢量/>形成垂直条件,随后保持该方向飞行;左转或右转机动过程中,目标飞行速度的变化按照公式(20)进行处理;
目标在左转机动过程中,转弯角速度ωL为:
目标在右转机动过程中,转弯角速度ωR为:
(4)建立目标最大过载爬升/俯冲机动模型;
目标在当前状态下,以最大过载noverload进行爬升机动,直至最终的飞行方向与初始目标线矢量和自身速度矢量/>所形成平面的法向量方向一致,或以最大过载noverload进行俯冲机动,直至最终的飞行方向与初始目标线矢量/>和自身速度矢量/>所形成平面的法向量方向相反,随后保持飞行状态;爬升/俯冲机动过程中,目标飞行速度的变化按照公式(20)进行处理;
目标在爬升机动过程中,最终的飞行速度的参考指向为:
目标爬升机动过程中的旋转轴线向量为:
爬升过程中的角度变化为:
目标在俯冲机动过程中,最终的飞行速度的参考指向为:
目标俯冲机动过程中的旋转轴线向量为:
以及:
公式(31)为俯冲过程中的角度变化。
4.根据权利要求1所述的面向近距空战格斗的三维动态发射区在线快速解算方法,其特征在于:
所述步骤3.2建立可攻击边界搜索模型的步骤如下:
步骤3.2.1:根据当前近距空战对抗态势,给定初始预估搜索边界的上界值和下界值为[a0,b0];
步骤3.2.2:在公式(32)确定的搜索区域内,根据导弹离轴发射角范围给定初始搜索角度值/>
步骤3.2.3:在给定离轴发射角上,基于当前敌我位置、速度信息,采用建立的导弹模型和目标机动模型,进行导弹制导打击仿真,依据导弹引战配合模型判定导弹最终能否命中杀伤目标,从而对搜索边界的上界值、下界值[ai,bi]进行更新,具体更新方式为:
其中,边界搜索点gp=ai+0.618(bi-ai);
步骤3.2.4:重复步骤3.2.3的边界搜索过程,直至可攻击边界的精度ε满足|bi-ai|≤ε,则以作为当前离轴发射角上的可攻击边界值;
步骤3.2.5:以Δτ为间隔进行离轴发射角更新,然后重复步骤3.2.3和3.2.4,直至完成离轴发射角遍历。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310067984.4A CN116579118B (zh) | 2023-02-06 | 2023-02-06 | 面向近距空战格斗的三维动态发射区在线快速解算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310067984.4A CN116579118B (zh) | 2023-02-06 | 2023-02-06 | 面向近距空战格斗的三维动态发射区在线快速解算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116579118A true CN116579118A (zh) | 2023-08-11 |
CN116579118B CN116579118B (zh) | 2024-02-13 |
Family
ID=87534698
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310067984.4A Active CN116579118B (zh) | 2023-02-06 | 2023-02-06 | 面向近距空战格斗的三维动态发射区在线快速解算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116579118B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2743479C1 (ru) * | 2020-03-25 | 2021-02-18 | Акционерное общество "Концерн радиостроения "Вега" | Способ и система определения наиболее благоприятных для атаки воздушных целей в режиме многоцелевого сопровождения |
CN113139331A (zh) * | 2021-03-29 | 2021-07-20 | 西北工业大学 | 一种基于贝叶斯网络的空空导弹态势感知与决策方法 |
RU2759057C1 (ru) * | 2020-09-08 | 2021-11-09 | Василий Васильевич Ефанов | Способ управления вооружением многофункциональных самолетов тактического назначения и система для его осуществления |
CN115186378A (zh) * | 2022-06-24 | 2022-10-14 | 西北工业大学 | 一种空战模拟环境中的战术控制距离实时解算方法 |
CN115238450A (zh) * | 2022-06-07 | 2022-10-25 | 西北工业大学 | 一种制导火箭弹行进间发射攻击区仿真计算方法 |
-
2023
- 2023-02-06 CN CN202310067984.4A patent/CN116579118B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2743479C1 (ru) * | 2020-03-25 | 2021-02-18 | Акционерное общество "Концерн радиостроения "Вега" | Способ и система определения наиболее благоприятных для атаки воздушных целей в режиме многоцелевого сопровождения |
RU2759057C1 (ru) * | 2020-09-08 | 2021-11-09 | Василий Васильевич Ефанов | Способ управления вооружением многофункциональных самолетов тактического назначения и система для его осуществления |
CN113139331A (zh) * | 2021-03-29 | 2021-07-20 | 西北工业大学 | 一种基于贝叶斯网络的空空导弹态势感知与决策方法 |
CN115238450A (zh) * | 2022-06-07 | 2022-10-25 | 西北工业大学 | 一种制导火箭弹行进间发射攻击区仿真计算方法 |
CN115186378A (zh) * | 2022-06-24 | 2022-10-14 | 西北工业大学 | 一种空战模拟环境中的战术控制距离实时解算方法 |
Non-Patent Citations (3)
Title |
---|
孔祥骏;商重阳;高正红;: "歼击机空战机动的仿真建模研究", 西北工业大学学报, no. 03, pages 18 - 23 * |
毕文豪 等: "新型远程空空导弹体系贡献度评估研究", 电光与控制, vol. 27, no. 8, pages 1 - 4 * |
王杰;丁达理;许明;韩博;雷磊;: "基于目标逃逸机动预估的空空导弹可发射区", 北京航空航天大学学报, no. 04, pages 85 - 97 * |
Also Published As
Publication number | Publication date |
---|---|
CN116579118B (zh) | 2024-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105910495B (zh) | 基于性能指标的面向效能的导弹武器系统设计方法 | |
CN111580556B (zh) | 一种时空约束下多无人机协同路径规划与制导方法 | |
US20080206718A1 (en) | Apparatus, method and computer program product for weapon flyout modeling and target damage assessment | |
CN107609307B (zh) | 一种考虑气弹和地球影响的远程火箭弹弹道分析方法 | |
CN111649624B (zh) | 一种空间微型精确制导武器控制方法 | |
CN111336871B (zh) | 一种基于迂回式飞行的垂直攻击制导方法 | |
CN112526872B (zh) | 带大落角约束中制导和末制导交接与制导信息的处理方法 | |
CN113625740A (zh) | 一种基于迁移学习鸽群优化的无人机空战博弈方法 | |
CN111176315B (zh) | 一种基于l1自适应控制的变体巡飞弹控制方法 | |
US4086841A (en) | Helical path munitions delivery | |
CN114035616B (zh) | 一种飞行器对移动目标打击控制方法及系统 | |
CN112731965B (zh) | 一种基于目标机动辨识的制导方法 | |
CN116579118B (zh) | 面向近距空战格斗的三维动态发射区在线快速解算方法 | |
CN116294811A (zh) | 基于多目标狼群算法的武装直升机火飞协同机制分析方法 | |
CN112729015B (zh) | 一种基于拦截弹导引律辨识的机动突防方法 | |
CN114153226A (zh) | 动态视线信息辅助的无人飞行器视场保持导引方法及系统 | |
Palumbo | Guest editor’s introduction: homing missile guidance and control | |
RU2602162C2 (ru) | Способ стрельбы реактивными снарядами реактивной системы залпового огня в условиях контрбатарейной борьбы | |
Zuoe et al. | Study on vertical attack penetration probability of anti-ship missile | |
Baba et al. | Air combat guidance law for an UCAV | |
RU2408832C1 (ru) | Способ стрельбы управляемым артиллерийским снарядом с лазерной полуактивной головкой самонаведения | |
Huang et al. | Modeling and simulation of the modified proportional guidance law of ship-to-air missiles | |
Zhao et al. | Study on fuze-guidance integration technology for improving air target striking capability of fortification storming/heat missiles | |
Li et al. | Anti-jamming Trajectory Planning of Infrared Imaging Air-to-air Missile | |
Nocoń et al. | Optimal compensator for anti-ship missile with vectorization of engine thrust |
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 |