CN111221350B - 吸气式高超声速飞行器巡航导弹的弹道设计方法及系统 - Google Patents

吸气式高超声速飞行器巡航导弹的弹道设计方法及系统 Download PDF

Info

Publication number
CN111221350B
CN111221350B CN201911402958.2A CN201911402958A CN111221350B CN 111221350 B CN111221350 B CN 111221350B CN 201911402958 A CN201911402958 A CN 201911402958A CN 111221350 B CN111221350 B CN 111221350B
Authority
CN
China
Prior art keywords
aircraft
mach number
cruise
cruising
attack
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
Application number
CN201911402958.2A
Other languages
English (en)
Other versions
CN111221350A (zh
Inventor
罗钦钦
涂正光
曾庆伟
潘霏
毛金娣
万方
郭春霞
许琦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
General Designing Institute of Hubei Space Technology Academy
Original Assignee
General Designing Institute of Hubei Space Technology Academy
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by General Designing Institute of Hubei Space Technology Academy filed Critical General Designing Institute of Hubei Space Technology Academy
Priority to CN201911402958.2A priority Critical patent/CN111221350B/zh
Publication of CN111221350A publication Critical patent/CN111221350A/zh
Application granted granted Critical
Publication of CN111221350B publication Critical patent/CN111221350B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Navigation (AREA)

Abstract

本发明公开了一种吸气式高超声速飞行器巡航导弹的弹道设计方法及系统,涉及导弹制导技术领域,该方法包括基于飞行器发动机的性能,在发动机工作窗口内,对飞行器的马赫数、高度和质量,按照预设密度进行等间距采样,得到采样数据;并根据采样数据中,马赫数和质量固定时,不同高度对应的燃料每秒消耗量,将最小燃料每秒消耗量对应的高度记为最优巡航高度,得到最优巡航高度集合;设定巡航马赫数随时间的变化规律、巡航马赫数随飞行器质量的变化规律、飞行器初始马赫数、飞行器初始高度和初始弹道倾角,作为仿真条件;基于设定的仿真条件进行弹道仿真计算,得到巡航导弹射程的最优巡航弹道。本发明得到的弹道射程较优。

Description

吸气式高超声速飞行器巡航导弹的弹道设计方法及系统
技术领域
本发明涉及导弹制导技术领域,具体涉及一种吸气式高超声速飞行器巡航导弹的弹道设计方法及系统。
背景技术
吸气式高超声速飞行器是采用冲压发动机为动力,在特定工作窗口内巡航飞行的一种新动力形式的飞行器,其在军事上的应用也是未来导弹发展方向之一,其在高超声速飞行器上的应用也是未来航空技术发展的主要方向之一。当前,世界上的主要航天军事强国都在全力开展相关研究。
由于现有技术手段限制,导致飞行器存在尺寸较小、空间有限、携带燃料不足等问题,使得飞行器的巡航导弹射程不足以支持军事应用的战略需求,因此需要采取特定的巡航弹道设计方法,使导弹射程最大化。现有巡航弹道设计方法中,通常采取等高平飞策略,满足升重平衡和推阻平衡,但是由于巡航导弹设计空间小,其射程提高有限,很难进行较大的提升。
发明内容
针对现有技术中存在的缺陷,本发明的目的在于提供一种吸气式高超声速飞行器巡航导弹的弹道设计方法及系统,得到的弹道射程较优。
为达到以上目的,本发明提供的一种吸气式高超声速飞行器巡航导弹的弹道设计方法,包括以下步骤:
基于飞行器发动机的性能,在发动机工作窗口内,对飞行器的马赫数、高度和质量,按照预设密度进行等间距采样,得到采样数据;
基于采样数据计算飞行器巡航状态下的燃料每秒消耗量,并根据采样数据中,马赫数和质量固定时,不同高度对应的燃料每秒消耗量,将最小燃料每秒消耗量对应的高度记为最优巡航高度,得到最优巡航高度集合;
设定巡航马赫数随时间的变化规律、巡航马赫数随飞行器质量的变化规律、飞行器初始马赫数、飞行器初始高度和初始弹道倾角,作为仿真条件;
基于设定的仿真条件进行弹道仿真计算,得到巡航导弹射程的最优巡航弹道。
在上述技术方案的基础上,
所述发动机的性能包括工作窗口、推力和比冲;
所述工作窗口满足以下公式:
{(H,Ma)|Ma∈[Mamin,Mamax],H∈[Hmin(Ma),Hmax(Ma)]}
其中,H为高度,Ma为马赫数,min表示最小值,max表示最大值;
所述推力满足以下公式:
其中,P表示推力,Isp表示比冲,φ表示发动机油门大小,表示燃料每秒消耗量,α表示飞行器攻角。
在上述技术方案的基础上,所述基于采样数据计算飞行器巡航状态下的燃料每秒消耗量,具体步骤包括:
在发动机的工作窗口内,确定飞行器等高平飞巡航状态下满足的升重平衡条件和推阻平衡条件,所述升重平衡条件和推阻平衡条件满足以下公式:
然后求解得到:
其中,m表示飞行器初始质量,g表示重力加速度,ρ表示大气密度,v表示飞行器的空速,Cl表示飞行器升力系数,αequ表示配平攻角,φequ表示配平油门,fα表示配平攻角为与马赫数、高度和飞行器初始质量相关的函数,fφ表示配平油门为与马赫数、高度和飞行器初始质量相关的函数;
计算巡航状态下的燃料每秒消耗量,计算公式为:
其中,表示巡航状态下的燃料每秒消耗量。
在上述技术方案的基础上,所述基于设定的仿真条件进行弹道仿真计算,得到巡航导弹射程的最优巡航弹道,具体步骤包括:
S401:在弹道仿真计算中,根据当前马赫数,质量和高度,从最优巡航高度集合中插值计算对应的最优巡航高度;
S402:计算飞行器当前的配平攻角和配平油门,并根据飞行器当前高度和最优巡航高度的偏差、飞行器当前马赫数与标称巡航马赫数的偏差、以及弹道倾角,对配平攻角和油门进行修正,得到攻角和油门,所述标称巡航马赫数为设定巡航马赫数随时间的变化规律,或巡航马赫数随飞行器质量的变化规律;
S403:根据攻角和油门计算发动机推力、气动升力、气动阻力和燃料每秒消耗量,得到飞行器的下一步状态;
S404:重复步骤S401~S403,直至燃料耗尽,从而得到巡航导弹射程的最优巡航弹道。
在上述技术方案的基础上,所述攻角和油门的计算公式为:
其中,αt表示攻角,φt表示油门,k1、k2和k3为正反馈系数,Hopt_t表示最优巡航高度,Maref_t表示标称巡航马赫数,αequ_t表示当前时刻t的配平攻角,Mat表示当前时刻t的飞行器马赫数,Ht表示当前时刻t的飞行器高度,φequ_t表示当前时刻t的配平油门。
本发明提供的一种吸气式高超声速飞行器巡航导弹的弹道设计系统,包括:
采样模块,其用于基于飞行器发动机的性能,在发动机工作窗口内,对飞行器的马赫数、高度和质量,按照预设密度进行等间距采样,得到采样数据;
记录模块,其用于基于采样数据计算飞行器巡航状态下的燃料每秒消耗量,并根据采样数据中,马赫数和质量固定时,不同高度对应的燃料每秒消耗量,将最小燃料每秒消耗量对应的高度记为最优巡航高度,得到最优巡航高度集合;
设定模块,其用于设定巡航马赫数随时间的变化规律、巡航马赫数随飞行器质量的变化规律、飞行器初始马赫数、飞行器初始高度和初始弹道倾角,作为仿真条件;
仿真模块,其用于基于设定的仿真条件进行弹道仿真计算,得到巡航导弹射程的最优巡航弹道。
在上述技术方案的基础上,
所述发动机的性能包括工作窗口、推力和比冲;
所述工作窗口满足以下公式:
{(H,Ma)|Ma∈[Mamin,Mamax],H∈[Hmin(Ma),Hmax(Ma)]}
其中,H为高度,Ma为马赫数,min表示最小值,max表示最大值;
所述推力满足以下公式:
其中,P表示推力,Isp表示比冲,φ表示发动机油门大小,表示燃料每秒消耗量,α表示飞行器攻角。
在上述技术方案的基础上,所述基于采样数据计算飞行器巡航状态下的燃料每秒消耗量,具体过程包括:
在发动机的工作窗口内,确定飞行器等高平飞巡航状态下满足的升重平衡条件和推阻平衡条件,所述升重平衡条件和推阻平衡条件满足以下公式:
然后求解得到:
其中,m表示飞行器初始质量,g表示重力加速度,ρ表示大气密度,v表示飞行器的空速,Cl表示飞行器升力系数,αequ表示配平攻角,φequ表示配平油门,fα表示配平攻角为与马赫数、高度和飞行器初始质量相关的函数,fφ表示配平油门为与马赫数、高度和飞行器初始质量相关的函数;
计算巡航状态下的燃料每秒消耗量,计算公式为:
其中,表示巡航状态下的燃料每秒消耗量。
在上述技术方案的基础上,所述基于设定的仿真条件进行弹道仿真计算,得到巡航导弹射程的最优巡航弹道,具体过程包括:
A:在弹道仿真计算中,根据当前马赫数,质量和高度,从最优巡航高度集合中插值计算对应的最优巡航高度;
B:计算飞行器当前的配平攻角和配平油门,并根据飞行器当前高度和最优巡航高度的偏差、飞行器当前马赫数与标称巡航马赫数的偏差、以及弹道倾角,对配平攻角和油门进行修正,得到攻角和油门,所述标称巡航马赫数为设定巡航马赫数随时间的变化规律,或巡航马赫数随飞行器质量的变化规律;
C:根据攻角和油门计算发动机推力、气动升力、气动阻力和燃料每秒消耗量,得到飞行器的下一步状态;
D:重复过程A~过程C,直至燃料耗尽,从而得到巡航导弹射程的最优巡航弹道。
在上述技术方案的基础上,所述攻角和油门的计算公式为:
其中,αt表示攻角,φt表示油门,k1、k2和k3为正反馈系数,Hopt_t表示最优巡航高度,Maref_t表示标称巡航马赫数,αequ_t表示当前时刻t的配平攻角,Mat表示当前时刻t的飞行器马赫数,Ht表示当前时刻t的飞行器高度,φequ_t表示当前时刻t的配平油门。
与现有技术相比,本发明的优点在于:改变了传统巡航飞行等高模式,根据当前马赫数选择燃料秒耗量最小所对应的最优巡航高度,使得每一时刻燃料消耗最小,弹道射程最优,设计变量缩减到马赫数一维,设计难度明显降低,设计自由度显著增加,且本发明方法得到的巡航弹道状态变化缓慢,飞行过程近似满足升重平衡和推阻平衡,与传统巡航状态类似,利于工程实现。
附图说明
图1为本发明实施例中一种吸气式高超声速飞行器巡航导弹的弹道设计方法的流程图。
具体实施方式
以下结合附图及实施例对本发明作进一步详细说明。
参见图1所示,本发明实施例提供的一种吸气式高超声速飞行器巡航导弹的弹道设计方法,包括以下步骤:
S1:基于飞行器发动机的性能,在发动机工作窗口内,对飞行器的马赫数、高度和质量,按照预设密度进行等间距采样,得到采样数据。发动机的性能包括工作窗口、推力和比冲
本发明实施例中,
工作窗口满足以下公式:
{(H,Ma)|Ma∈[Mamin,Mamax],H∈[Hmin(Ma),Hmax(Ma)]}
其中,H为高度,Ma为马赫数,min表示最小值,max表示最大值;
推力满足以下公式:
其中,P表示推力,Isp表示比冲,φ表示发动机油门大小,表示燃料每秒消耗量,α表示飞行器攻角。
按照预设密度进行等间距采样,具体为:在发动机工作窗口马赫数和高度范围内,选择等间距马赫数采样点{Mai|Ma1=Mamin,Man=Mamax,Mai<Mai+1,i=1,…,n},等间距高度采样点{Hj|H1=Hmin,Hm=Hmax,Hj<Hj+1,j=1,…,m},在空载质量me和满载质量mf之间等间距采样{mk|m1=me,ml=mf,mk<mk+1,k=1,…,l}。
S2:基于采样数据计算飞行器巡航状态下的燃料每秒消耗量,并根据采样数据中,马赫数和质量固定时,不同高度对应的燃料每秒消耗量,将最小燃料每秒消耗量对应的高度记为最优巡航高度,得到最优巡航高度集合;
本发明实施例中,基于采样数据计算飞行器巡航状态下的燃料每秒消耗量,具体步骤包括:
在发动机的工作窗口内,确定飞行器等高平飞巡航状态下满足的升重平衡条件和推阻平衡条件,所述升重平衡条件和推阻平衡条件满足以下公式:
然后求解得到:
其中,m表示飞行器初始质量,g表示重力加速度,ρ表示大气密度,v表示飞行器的空速,Cl表示飞行器升力系数,αequ表示配平攻角,φequ表示配平油门,fα表示配平攻角为与马赫数、高度和飞行器初始质量相关的函数,fφ表示配平油门为与马赫数、高度和飞行器初始质量相关的函数;
计算巡航状态下的燃料每秒消耗量,计算公式为:
其中,表示巡航状态下的燃料每秒消耗量。
S3:设定巡航马赫数随时间的变化规律、巡航马赫数随飞行器质量的变化规律、飞行器初始马赫数、飞行器初始高度和初始弹道倾角,作为仿真条件;
S4:基于设定的仿真条件进行弹道仿真计算,得到巡航导弹射程的最优巡航弹道。该步骤具体为:
S401:在弹道仿真计算中,根据当前马赫数,质量和高度,从最优巡航高度集合中插值计算对应的最优巡航高度;
S402:计算飞行器当前的配平攻角和配平油门,并根据飞行器当前高度和最优巡航高度的偏差、飞行器当前马赫数与标称巡航马赫数的偏差、以及弹道倾角,对配平攻角和油门进行修正,得到攻角和油门,所述标称巡航马赫数为设定巡航马赫数随时间的变化规律,或巡航马赫数随飞行器质量的变化规律;
S403:根据攻角和油门计算发动机推力、气动升力、气动阻力和燃料每秒消耗量,得到飞行器的下一步状态;
S404:重复步骤S401~S403,直至燃料耗尽,从而得到巡航导弹射程的最优巡航弹道。
攻角和油门的计算公式为:
其中,αt表示攻角,φt表示油门,k1、k2和k3为正反馈系数,Hopt_t表示最优巡航高度,Maref_t表示标称巡航马赫数,αequ_t表示当前时刻t的配平攻角,Mat表示当前时刻t的飞行器马赫数,Ht表示当前时刻t的飞行器高度,φequ_t表示当前时刻t的配平油门。
本发明实施例中,在标准条件下使用三自由度弹道模型计算弹道。标准条件包括:a)地球模型为旋转椭球体;b)地球表面的重力加速度为9.80665m/s2;c)大气条件为国家标准大气;d)所有高度上无风;e)符合瞬时平衡原理。
本发明实施例的吸气式高超声速飞行器巡航导弹的弹道设计方法,改变了传统巡航飞行等高模式,根据当前马赫数选择燃料秒耗量最小所对应的最优巡航高度,使得每一时刻燃料消耗最小,弹道射程最优,设计变量缩减到马赫数一维,设计难度明显降低,设计自由度显著增加,且本发明方法得到的巡航弹道状态变化缓慢,飞行过程近似满足升重平衡和推阻平衡,与传统巡航状态类似,利于工程实现。
本发明实施例提供的一种吸气式高超声速飞行器巡航导弹的弹道设计系统,包括采样模块、记录模块、设定模块和仿真模块。采样模块用于基于飞行器发动机的性能,在发动机工作窗口内,对飞行器的马赫数、高度和质量,按照预设密度进行等间距采样,得到采样数据;记录模块用于基于采样数据计算飞行器巡航状态下的燃料每秒消耗量,并根据采样数据中,马赫数和质量固定时,不同高度对应的燃料每秒消耗量,将最小燃料每秒消耗量对应的高度记为最优巡航高度,得到最优巡航高度集合;设定模块用于设定巡航马赫数随时间的变化规律、巡航马赫数随飞行器质量的变化规律、飞行器初始马赫数、飞行器初始高度和初始弹道倾角,作为仿真条件;仿真模块用于基于设定的仿真条件进行弹道仿真计算,得到巡航导弹射程的最优巡航弹道。
本发明实施例中,发动机的性能包括工作窗口、推力和比冲;
所述工作窗口满足以下公式:
{(H,Ma)|Ma∈[Mamin,Mamax],H∈[Hmin(Ma),Hmax(Ma)]}
其中,H为高度,Ma为马赫数,min表示最小值,max表示最大值;
所述推力满足以下公式:
其中,P表示推力,Isp表示比冲,φ表示发动机油门大小,表示燃料每秒消耗量,α表示飞行器攻角。
本发明实施例中,基于采样数据计算飞行器巡航状态下的燃料每秒消耗量,具体过程包括:
在发动机的工作窗口内,确定飞行器等高平飞巡航状态下满足的升重平衡条件和推阻平衡条件,所述升重平衡条件和推阻平衡条件满足以下公式:
然后求解得到:
其中,m表示飞行器初始质量,g表示重力加速度,ρ表示大气密度,v表示飞行器的空速,Cl表示飞行器升力系数,αequ表示配平攻角,φequ表示配平油门,fα表示配平攻角为与马赫数、高度和飞行器初始质量相关的函数,fφ表示配平油门为与马赫数、高度和飞行器初始质量相关的函数;
计算巡航状态下的燃料每秒消耗量,计算公式为:
其中,表示巡航状态下的燃料每秒消耗量。
本发明实施例中,基于设定的仿真条件进行弹道仿真计算,得到巡航导弹射程的最优巡航弹道,具体过程包括:
A:在弹道仿真计算中,根据当前马赫数,质量和高度,从最优巡航高度集合中插值计算对应的最优巡航高度;
B:计算飞行器当前的配平攻角和配平油门,并根据飞行器当前高度和最优巡航高度的偏差、飞行器当前马赫数与标称巡航马赫数的偏差、以及弹道倾角,对配平攻角和油门进行修正,得到攻角和油门,所述标称巡航马赫数为设定巡航马赫数随时间的变化规律,或巡航马赫数随飞行器质量的变化规律;
C:根据攻角和油门计算发动机推力、气动升力、气动阻力和燃料每秒消耗量,得到飞行器的下一步状态;
D:重复过程A~过程C,直至燃料耗尽,从而得到巡航导弹射程的最优巡航弹道。
本发明实施例中,攻角和油门的计算公式为:
其中,αt表示攻角,φt表示油门,k1、k2和k3为正反馈系数,Hopt_t表示最优巡航高度,Maref_t表示标称巡航马赫数,αequ_t表示当前时刻t的配平攻角,Mat表示当前时刻t的飞行器马赫数,Ht表示当前时刻t的飞行器高度,φequ_t表示当前时刻t的配平油门。
本发明不局限于上述实施方式,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围之内。本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (4)

1.一种吸气式高超声速飞行器巡航导弹的弹道设计方法,其特征在于,包括以下步骤:
基于飞行器发动机的性能,在发动机工作窗口内,对飞行器的马赫数、高度和质量,按照预设密度进行等间距采样,得到采样数据;
基于采样数据计算飞行器巡航状态下的燃料每秒消耗量,并根据采样数据中,马赫数和质量固定时,不同高度对应的燃料每秒消耗量,将最小燃料每秒消耗量对应的高度记为最优巡航高度,得到最优巡航高度集合;
设定巡航马赫数随时间的变化规律、巡航马赫数随飞行器质量的变化规律、飞行器初始马赫数、飞行器初始高度和初始弹道倾角,作为仿真条件;
基于设定的仿真条件进行弹道仿真计算,得到巡航导弹射程的最优巡航弹道;
其中,
所述发动机的性能包括工作窗口、推力和比冲;
所述工作窗口满足以下公式:
{(H,Ma)|Ma∈[Mamin,Mamax],H∈[Hmin(Ma),Hmax(Ma)]}
其中,H为高度,Ma为马赫数,min表示最小值,max表示最大值;
所述推力满足以下公式:
Figure FDA0004105724090000011
其中,P表示推力,Isp表示比冲,φ表示发动机油门大小,
Figure FDA0004105724090000012
表示燃料每秒消耗量,α表示飞行器攻角;
其中,所述基于采样数据计算飞行器巡航状态下的燃料每秒消耗量,具体步骤包括:
在发动机的工作窗口内,确定飞行器等高平飞巡航状态下满足的升重平衡条件和推阻平衡条件,所述升重平衡条件和推阻平衡条件满足以下公式:
Figure FDA0004105724090000021
然后求解得到:
Figure FDA0004105724090000022
其中,m表示飞行器初始质量,g表示重力加速度,ρ表示大气密度,v表示飞行器的空速,Cl表示飞行器升力系数,αequ表示配平攻角,φequ表示配平油门,fα表示配平攻角为与马赫数、高度和飞行器初始质量相关的函数,fφ表示配平油门为与马赫数、高度和飞行器初始质量相关的函数;
计算巡航状态下的燃料每秒消耗量,计算公式为:
Figure FDA0004105724090000023
其中,
Figure FDA0004105724090000024
表示巡航状态下的燃料每秒消耗量;
其中,所述基于设定的仿真条件进行弹道仿真计算,得到巡航导弹射程的最优巡航弹道,具体步骤包括:
S401:在弹道仿真计算中,根据当前马赫数,质量和高度,从最优巡航高度集合中插值计算对应的最优巡航高度;
S402:计算飞行器当前的配平攻角和配平油门,并根据飞行器当前高度和最优巡航高度的偏差、飞行器当前马赫数与标称巡航马赫数的偏差、以及弹道倾角,对配平攻角和油门进行修正,得到攻角和油门,所述标称巡航马赫数为设定巡航马赫数随时间的变化规律,或巡航马赫数随飞行器质量的变化规律;
S403:根据攻角和油门计算发动机推力、气动升力、气动阻力和燃料每秒消耗量,得到飞行器的下一步状态;
S404:重复步骤S401~S403,直至燃料耗尽,从而得到巡航导弹射程的最优巡航弹道。
2.如权利要求1所述的一种吸气式高超声速飞行器巡航导弹的弹道设计方法,其特征在于,所述攻角和油门的计算公式为:
Figure FDA0004105724090000031
其中,αt表示攻角,φt表示油门,k1、k2和k3为正反馈系数,Hopt_t表示最优巡航高度,Maref_t表示标称巡航马赫数,αequ_t表示当前时刻t的配平攻角,Mat表示当前时刻t的飞行器马赫数,Ht表示当前时刻t的飞行器高度,φequ_t表示当前时刻t的配平油门。
3.一种吸气式高超声速飞行器巡航导弹的弹道设计系统,其特征在于,包括:
采样模块,其用于基于飞行器发动机的性能,在发动机工作窗口内,对飞行器的马赫数、高度和质量,按照预设密度进行等间距采样,得到采样数据;
记录模块,其用于基于采样数据计算飞行器巡航状态下的燃料每秒消耗量,并根据采样数据中,马赫数和质量固定时,不同高度对应的燃料每秒消耗量,将最小燃料每秒消耗量对应的高度记为最优巡航高度,得到最优巡航高度集合;
设定模块,其用于设定巡航马赫数随时间的变化规律、巡航马赫数随飞行器质量的变化规律、飞行器初始马赫数、飞行器初始高度和初始弹道倾角,作为仿真条件;
仿真模块,其用于基于设定的仿真条件进行弹道仿真计算,得到巡航导弹射程的最优巡航弹道;
其中,
所述发动机的性能包括工作窗口、推力和比冲;
所述工作窗口满足以下公式:
{(H,Ma)|Ma∈[Mamin,Mamax],H∈[Hmin(Ma),Hmax(Ma)]}
其中,H为高度,Ma为马赫数,min表示最小值,max表示最大值;
所述推力满足以下公式:
Figure FDA0004105724090000041
其中,P表示推力,Isp表示比冲,φ表示发动机油门大小,
Figure FDA0004105724090000042
表示燃料每秒消耗量,α表示飞行器攻角;
其中,所述基于采样数据计算飞行器巡航状态下的燃料每秒消耗量,具体过程包括:
在发动机的工作窗口内,确定飞行器等高平飞巡航状态下满足的升重平衡条件和推阻平衡条件,所述升重平衡条件和推阻平衡条件满足以下公式:
Figure FDA0004105724090000043
然后求解得到:
Figure FDA0004105724090000044
其中,m表示飞行器初始质量,g表示重力加速度,ρ表示大气密度,v表示飞行器的空速,Cl表示飞行器升力系数,αequ表示配平攻角,φequ表示配平油门,fα表示配平攻角为与马赫数、高度和飞行器初始质量相关的函数,fφ表示配平油门为与马赫数、高度和飞行器初始质量相关的函数;
计算巡航状态下的燃料每秒消耗量,计算公式为:
Figure FDA0004105724090000051
其中,
Figure FDA0004105724090000052
表示巡航状态下的燃料每秒消耗量;
其中,所述基于设定的仿真条件进行弹道仿真计算,得到巡航导弹射程的最优巡航弹道,具体过程包括:
A:在弹道仿真计算中,根据当前马赫数,质量和高度,从最优巡航高度集合中插值计算对应的最优巡航高度;
B:计算飞行器当前的配平攻角和配平油门,并根据飞行器当前高度和最优巡航高度的偏差、飞行器当前马赫数与标称巡航马赫数的偏差、以及弹道倾角,对配平攻角和油门进行修正,得到攻角和油门,所述标称巡航马赫数为设定巡航马赫数随时间的变化规律,或巡航马赫数随飞行器质量的变化规律;
C:根据攻角和油门计算发动机推力、气动升力、气动阻力和燃料每秒消耗量,得到飞行器的下一步状态;
D:重复过程A~过程C,直至燃料耗尽,从而得到巡航导弹射程的最优巡航弹道。
4.如权利要求3所述的一种吸气式高超声速飞行器巡航导弹的弹道设计系统,其特征在于,所述攻角和油门的计算公式为:
Figure FDA0004105724090000053
其中,αt表示攻角,φt表示油门,k1、k2和k3为正反馈系数,Hopt_t表示最优巡航高度,Maref_t表示标称巡航马赫数,αequ_t表示当前时刻t的配平攻角,Mat表示当前时刻t的飞行器马赫数,Ht表示当前时刻t的飞行器高度,φequ_t表示当前时刻t的配平油门。
CN201911402958.2A 2019-12-30 2019-12-30 吸气式高超声速飞行器巡航导弹的弹道设计方法及系统 Active CN111221350B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911402958.2A CN111221350B (zh) 2019-12-30 2019-12-30 吸气式高超声速飞行器巡航导弹的弹道设计方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911402958.2A CN111221350B (zh) 2019-12-30 2019-12-30 吸气式高超声速飞行器巡航导弹的弹道设计方法及系统

Publications (2)

Publication Number Publication Date
CN111221350A CN111221350A (zh) 2020-06-02
CN111221350B true CN111221350B (zh) 2023-05-02

Family

ID=70808288

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911402958.2A Active CN111221350B (zh) 2019-12-30 2019-12-30 吸气式高超声速飞行器巡航导弹的弹道设计方法及系统

Country Status (1)

Country Link
CN (1) CN111221350B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114610057A (zh) * 2022-02-25 2022-06-10 北京航空航天大学 一种高马赫飞行器机动突防策略设计方法
CN115328191B (zh) * 2022-07-15 2023-09-12 北京星途探索科技有限公司 一种掠地巡航靶标巡航控制方法、系统、设备及存储介质
CN115993847B (zh) * 2023-03-23 2023-06-06 北京航天众信科技有限公司 一种掠地巡航弹的定高方法、装置、介质、设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4490793A (en) * 1981-12-21 1984-12-25 Sperry Corporation Cruise speed control for aircraft performance management system
RU2537193C1 (ru) * 2013-09-27 2014-12-27 Федеральное государственное унитарное предприятие "Центральный научно-исследовательский институт машиностроения" (ФГУП ЦНИИмаш) Способ управления спуском космического аппарата в атмосфере планеты
CN110334439A (zh) * 2019-07-04 2019-10-15 南京航空航天大学 一种利用气动力辅助的超轨道速度飞行技术方案

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB812748A (en) * 1956-01-23 1959-04-29 Bendix Aviat Corp Improvements in computers for determining the optimum cruising altitude of an aircraft
BG43502A1 (en) * 1985-08-19 1988-06-15 Jjordan D Zaprjanov Method and system for control of fuel consumption
US7774106B2 (en) * 2006-12-22 2010-08-10 Pratt - Whitney Canada Corp. Cruise control FADEC logic
US8818696B2 (en) * 2011-03-23 2014-08-26 Ge Aviation Systems Llc Method and system for aerial vehicle trajectory management
EP3065019B1 (en) * 2015-03-04 2017-09-13 The Boeing Company Method for optimum maximum range cruise speed in an aircraft
CN110413000B (zh) * 2019-05-28 2020-10-13 北京航空航天大学 一种基于深度学习的高超声速飞行器再入预测校正容错制导方法
CN110276479B (zh) * 2019-05-31 2023-01-03 南京航空航天大学 飞机质量变化的巡航阶段燃油消耗预测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4490793A (en) * 1981-12-21 1984-12-25 Sperry Corporation Cruise speed control for aircraft performance management system
RU2537193C1 (ru) * 2013-09-27 2014-12-27 Федеральное государственное унитарное предприятие "Центральный научно-исследовательский институт машиностроения" (ФГУП ЦНИИмаш) Способ управления спуском космического аппарата в атмосфере планеты
CN110334439A (zh) * 2019-07-04 2019-10-15 南京航空航天大学 一种利用气动力辅助的超轨道速度飞行技术方案

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Zhifu Shi, Daizhong Zhou, Shunhong Wang, Lefei Pan.Simulation of Skip-Glide Trajectory for Hypersonic Vehicle in Near Space.Proceedings of 2012 System Simulation and Scientific Computing.2012,论文第250-258页. *
Zhi-kai Zhang ; Qing Guo.Terminal guidance method in glide attack phase of aerodynamic missile.Proceedings of the 10th World Congress on Intelligent Control and Automation.2012,论文第1150-1155页. *
刘星.高超声速巡航飞行器轨迹优化与制导方法研究.中国优秀硕士学位论文全文数据库工程科技II辑.2018,(第7期),论文第76-94页. *

Also Published As

Publication number Publication date
CN111221350A (zh) 2020-06-02

Similar Documents

Publication Publication Date Title
CN111221350B (zh) 吸气式高超声速飞行器巡航导弹的弹道设计方法及系统
CN107844128B (zh) 一种基于复合比例导引的高超声速飞行器巡航段制导方法
CN110334439B (zh) 一种利用气动力辅助的超轨道速度飞行技术方法
CN109190248B (zh) 一种用于滑翔飞行器的滑翔射程解析方法及解析系统
CN113589847B (zh) 一种柔性飞行器飞行半径确定方法
CN107315419B (zh) 一种实现飞艇远距离定点悬停的三段式控制策略
CN113569429B (zh) 短距垂直起降飞行器三维飞行包线计算方法及装置
CN114004021B (zh) 用于飞行管理系统性能管理的巡航燃油流量计算方法
CN116301028B (zh) 基于吸气式高超声速平台的多约束在线飞行轨迹规划中段导引方法
CN115017721B (zh) 一种飞机巡航特性辨识方法、装置及飞行控制系统
CN114384935B (zh) 一种无人航空载运飞行器多约束气动减速控制方法
CN116186904A (zh) 一种具有升力面在流体中运动的机械总体气动布局方法
CN114329976A (zh) 一种螺旋桨飞机巡航状态气动数据分析处理方法
CN111914344B (zh) 一种飞机极曲线试飞数据处理方法
CN114935277A (zh) 一种滑翔增程制导炮弹理想弹道的在线规划方法
CN109117544B (zh) 一种天地往返飞行器全轨迹的优化方法
CN117807915B (zh) 一种喷气式飞机给定重量最大航程巡航设计方法
CN112364433A (zh) 一种高效的固定翼飞行器飞行动力学模型配平方法
Relf Recent aerodynamic developments
CN112109921B (zh) 一种月面基地发射低轨道环月飞行器燃耗最省控制方法
Papoči et al. Modelling piston engine aircraft performance using pilot's operating handbook
CN113608547B (zh) 一种滑翔制导炸弹基于风补偿的机载火控方法及系统
CN115933431A (zh) 一种数字对抗仿真环境下的亚音速飞行器快速建模方法
CN115617057B (zh) 四倾转旋翼飞行器着舰末端纵向俯仰角同步控制方法
CN117761341A (zh) 一种无人机的远航速度计算方法

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