CN109240323A - 一种实时解析构造的空天飞行器再入制导方法 - Google Patents
一种实时解析构造的空天飞行器再入制导方法 Download PDFInfo
- Publication number
- CN109240323A CN109240323A CN201811302375.8A CN201811302375A CN109240323A CN 109240323 A CN109240323 A CN 109240323A CN 201811302375 A CN201811302375 A CN 201811302375A CN 109240323 A CN109240323 A CN 109240323A
- Authority
- CN
- China
- Prior art keywords
- aircraft
- drag acceleration
- energy
- acceleration
- feature profile
- 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
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000010276 construction Methods 0.000 title claims abstract description 34
- 230000001133 acceleration Effects 0.000 claims abstract description 193
- 238000012937 correction Methods 0.000 claims abstract description 17
- 238000000926 separation method Methods 0.000 claims description 4
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 241000405147 Hermes Species 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
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)
- Traffic Control Systems (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种实时解析构造的空天飞行器再入制导方法,包括步骤如下:1)待飞行器弹道倾角为零后,确定飞行器能量‑阻力加速度走廊上界和走廊下界,2)生成飞行器能量阻力加速度特征剖面,3)修正飞行器能量阻力加速度特征剖面,获得飞行器能量阻力加速度特征剖面的修正量,4)解析构造再入制导模型,5)确定飞行器倾侧角指令,完成飞行器再入制导工作。本发明的制导方法自适应能力强,精度高,计算简单,易于工程实现。
Description
技术领域
本发明涉及一种实时解析构造的空天飞行器再入制导方法,属于制导控制技术领域,特别适用于空天飞行器、可重复使用运载器、高超声速飞行器的再入段的制导。
背景技术
空天飞行器(Aerospace Vehicle,ASV)是既能在航空领域又能在航天领域工作的新型飞行器,它结合了航空技术与航天技术。空天飞行器能像普通飞机一样水平起飞,以高超声速在大气层内飞行,并可直接加速进入地球轨道,成为航天飞行器,返回大气层后,像飞机一样在机场着陆。
再入段是空天飞行器返回着陆过程的一个重要阶段,是空天飞行器安全返回与顺利着陆的重要保障。再入段位于离轨段之后,末端能量管理段之前,其主要目的是逐步耗散飞行器的动能和势能量,控制飞行器以合适的高度、速度和航向到达末端能量管理窗口。典型的初期段飞行起始于速度约7511m/s、海拔高度约120000m,结束于末端能量管理窗口,终端标称速度约765m/s、终端标称高度约26700m。
空天飞行器再入段飞行高度和速度变化范围大,飞行状态变化剧烈。同时,再入动压、过载和热流密度等限制严重,决定了飞行走廊窄,还存高空大气环境和气动力系数的不确定性大的问题,过程约束条件强以及参数不确定性大的特点,制导控制难度大。
再入制导方法通常分为标准轨迹制导和预测校正制导两大类。标准轨迹制导方法通过设计导引律以实时跟踪事先确定的标准再入轨迹,对GNC控制器的性能要求比较低,但是较大的初始轨迹偏差以及再入过程环境扰动等因素的影响较大,可能会导致制导性能严重退化。预测校正制导方法则是通过预测再入飞行器的落点位置偏差来对制导指令进行实时修正,能够显著降低初始散布误差和飞行过程中各种偏差对制导性能的影响,提高制导精度。因此,预测校正制导成为空天飞行器再入制导方法的发展趋势。然而,目前大量的预测校正制导方法均建立在数值积分弹道方程的基础上,该类方法的计算量大,给工程实施造成了巨大的难度。
在工程实施中,航天飞机采用基于阻力加速度剖面的标准弹道再入制导方法计算量小,对于初始误差和过程干扰较小的再入任务可以实现较好的制导效果,并进一步应用在了Hermes飞行器、X-33飞行器。航天飞机采用的基于阻力加速度剖面的标准弹道再入制导方法难以摆脱鲁棒性和适应性差的固有缺陷。为满足新一代空天飞行器对自主性、精确性、安全性和可靠性的苛刻要求,世界各国开展了大量的新型升力式再入制导方法的研究和验证,比较典型的为改进的加速度再入制导算法,但该方法选用三段线性曲线组成的阻力加速度剖面限制了飞行器的再入飞行能力,同时由于三段曲线不连续,不利于飞行器跟踪参考阻力剖面,此外,用到的搜索算法比较消耗机时。
发明内容
本发明的技术解决问题是:克服现有技术的不足之处,提供一种实时解析构造的空天飞行器再入制导方法,该方法能够获得高的制导精度,提高了空天飞行器的纵横航程在线调整能力,同时计算简单,工程实现容易。
本发明的技术方案是:
一种实时解析构造的空天飞行器再入制导方法,包括步骤如下:
1)判断飞行器弹道倾角是否为零,若飞行器弹道倾角不为零,则进入步骤2),若飞行器弹道倾角为零,则进入步骤3);
2)将保持飞行器弹道倾角为零作为制导信息,将所述制导信息传递给飞行器控制***,直至飞行器弹道倾角零后,进入步骤3);
3)确定飞行器能量-阻力加速度走廊上界和走廊下界,根据所述飞行器能量-阻力加速度走廊上界和走廊下界,生成飞行器能量阻力加速度特征剖面;
4)修正飞行器能量阻力加速度特征剖面,获得飞行器能量阻力加速度特征剖面的修正量;
5)根据飞行器能量阻力加速度特征剖面的修正量,解析构造再入制导模型;
6)根据所述解析构造的再入制导模型,确定飞行器倾侧角指令;
7)将所述确定的倾侧角指令作为制导信息发送给飞行器控制***,完成飞行器再入制导工作。
所述步骤3)确定飞行器能量-阻力加速度走廊上界根据下列等式确定,具体为:
所述步骤3)飞行器能量-阻力加速度走廊下界根据下列等式确定,具体为:
其中,qmax为最大动压约束,nmax为最大过载约束,为最大热流密度约束,Dqmax为动压约束下的飞行器能量-阻力加速度走廊上界,Dnmax为过载约束下的飞行器能量-阻力加速度走廊上界,为热流密度约束约束下的飞行器能量-阻力加速度走廊上界,DEQ min为平衡滑翔约束下的飞行器能量-阻力加速度走廊下界,CL为气动升力系数,kQ为与飞行器外形相关的常值参数,r为飞行器地心距,Sref为飞行器参考面积,m为飞行器质量,CD为气动阻力系数,kQ为与飞行器外形相关的常值参数,v为飞行器速度,g为引力加速度。
所述步骤3)生成的飞行器能量阻力加速度特征剖面,具体为:
根据所述飞行器能量-阻力加速度走廊上界和走廊下界,以飞行器能量为横坐标,以飞行器阻力加速度为纵坐标,生成飞行器能量阻力加速度特征剖面;
A)当E≥Ec时,所述飞行器能量阻力加速度特征剖面包括n段二次曲线,n为正整数;其中,第i段二次曲线对应的飞行器能量阻力加速度特征剖面,具体为:
Di0=Ci1E2+Ci2E+Ci3,Eif<E≤Eib;
其中,i=1,2,3,...,n,Eif、Eib分别为第i条阻力加速度二次曲线的终端能量和起始能量,Ci1、Ci2和Ci3为第i条阻力加速度二次曲线的系数,Ec为二次曲线和一次曲线分界点能量,E为飞行器能量;
B)当E<Ec时,所述飞行器能量阻力加速度特征剖面为一次曲线,具体为;
其中,Dc为Ec点对应的飞行器阻力加速度,Df为Ef点对应的飞行器阻力加速度,Ef为再入终端的标称能量,vf为再入终端标称速度,hf为再入终端标称高度,gf为再入终端标称引力加速度,ρf为再入终端标称大气密度。
所述步骤4)修正飞行器能量阻力加速度特征剖面的方法,具体为:
41)根据所述飞行器能量阻力加速度特征剖面,确定飞行器能量阻力加速度特征剖面对应的飞行航程sD0;
42)确定飞行器初始待飞航程st0;
43)修正所述步骤3)生成的飞行器能量阻力加速度特征剖面,使得所述步骤41)确定的sD0与所述步骤42)确定的st0一致,获得修正后的所述飞行器能量阻力加速度特征剖面。
所述步骤41)飞行器能量阻力加速度特征剖面对应的飞行航程sD0,具体为:
所述步骤42)飞行器初始待飞航程st0,具体为:
st0=r0arccos(sinφ0sinφf+cosφ0cosφfcos(λf-λ0)),
其中,φf为再入终端标称纬度,λf为再入终端标称经度,φ0为弹道倾角变化为零时飞行器所处的纬度,λ0为弹道倾角变化为零时飞行器所处的经度,r0为弹道倾角变化为零时飞行器的地心距。
所述步骤5)解析构造再入制导模型,具体为:
sDt=sD-st,
st=rarccos(sinφsinφf+cosφcosφfcos(λf-λ)),
其中,θ为飞行器弹道倾角,u=cosσ,D0为当前能量对应的初始阻力加速度剖面下的阻力加速度,D为飞行器阻力加速度,ΔD为飞行器能量阻力加速度特征剖面的修正量,L为飞行器升力加速度,hs为大气密度特征参数,DN=Dl0+ΔD,r为当前飞行器地心距;φ为当前飞行器纬度;λ为当前飞行器经度;
EN确定方法如下
a)当E>Ec时:
EN=Ec,
b)当Ef<En≤Ec时:
EN=En,
所述步骤6),确定飞行器倾侧角指令具体为:
σ=σabssign(σ),
σabs=arccos(u),
其中,a0>0,a1>0,a2>0,ε0>0,其中G-1为G的逆矩阵,sign()为符号函数,σi-1为上一制导周期的倾侧角指令,Δψd为侧向方位边界,Δψ为侧向方位
本发明与现有技术相比的优点在于:
1)本发明采用阻力加速度特征剖面自主生成,以及基于此的阻力加速度剖面实时调整的航程误差调整方式,使得阻力加速度剖面曲线连续且易于跟踪,起到了航程偏差在线实时削除的作用,提高了制导精度;
2)本发明采用基于实时构造的再入制导方式,起到了提高了鲁棒性和适应性,避免了传统的标准弹道再入制导方法的固有缺陷,飞行器的纵横航程在线调整能力强和制导精度高;
3)本发明采用了解析计算实时制导参数的方式,减少了飞控机数据运算能力的要求,避免了数值预测校正计算量大,使得该方法易于工程实施。
附图说明
图1为方法流程图;
图2为本发明实施例的飞行器能量-阻力加速度走廊示意图;
图3为本发明实施例的侧向方位边界示意图;
图4为本发明实施例的速度曲线;
图5为本发明实施例的全程高度曲线;
图6为本发明实施例的全程地轨迹曲线;
图7为本发明实施例的全程弹道倾角曲线;
图8为本发明实施例的全程的倾侧角指令;
图9为本发明实施例的全程的阻力加速度曲线。
具体实施方式
从再入制导开始直到弹道倾角大小减小到零,为飞行器航迹拉起段;从飞行器航迹拉起段结束到飞行器能量减少至再入终端能量,为飞行器滑翔飞行段。
如图1所示,一种实时解析构造的空天飞行器再入制导方法,包括步骤如下:
1)判断飞行器弹道倾角是否为零,若飞行器弹道倾角不为零,则飞行器的飞行阶段为航迹拉起段,进入步骤2),若飞行器弹道倾角为零,则飞行器的飞行阶段为滑翔飞行段,进入步骤3);
2)将保持飞行器弹道倾角为零作为制导信息,将所述制导信息传递给飞行器控制***,直至飞行器弹道倾角零后,进入步骤3);
3)确定飞行器能量-阻力加速度走廊上界和走廊下界,根据所述飞行器能量-阻力加速度走廊上界和走廊下界,生成飞行器能量阻力加速度特征剖面;
4)修正飞行器能量阻力加速度特征剖面,获得飞行器能量阻力加速度特征剖面的修正量;
5)根据飞行器能量阻力加速度特征剖面的修正量,即根据获得的修正后当前飞行器能量阻力加速度特征剖面对航程偏差的影响,解析构造再入制导模型;
6)根据所述解析构造的再入制导模型,确定飞行器倾侧角指令;
7)将所述确定的倾侧角指令作为制导信息发送给飞行器控制***,完成飞行器再入制导工作。
步骤3)确定飞行器能量-阻力加速度走廊上界根据下列等式确定,具体为:
所述步骤3)飞行器能量-阻力加速度走廊下界根据下列等式确定,具体为:
其中,qmax为最大动压约束,nmax为最大过载约束,为最大热流密度约束,Dqmax为动压约束下的飞行器能量-阻力加速度走廊上界,Dnmax为过载约束下的飞行器能量-阻力加速度走廊上界,为热流密度约束约束下的飞行器能量-阻力加速度走廊上界,DEQ min为平衡滑翔约束下的飞行器能量-阻力加速度走廊下界,CL为气动升力系数,kQ为与飞行器外形相关的常值参数,r为飞行器地心距,Sref为飞行器参考面积,m为飞行器质量,CD为气动阻力系数,kQ为与飞行器外形相关的常值参数,v为飞行器速度,g为引力加速度。
飞行器能量阻力加速度特征剖面是以飞行器能量为横坐标,飞行器阻力加速度为纵坐标的一条连续光滑曲线,具体由多个光滑连续的二次曲线和一次曲线组成,光滑曲线位于步骤3确定飞行器能量-阻力加速度走廊上界和走廊下界范围内部。其中,v为飞行器速度,g为引力加速度,h为飞行器高度,大气密度ρ0为海平面大气密度,hs为大气特征常数,Sref为飞行器参考面积,m为飞行器质量,CD为气动阻力系数。
步骤3)生成的飞行器能量阻力加速度特征剖面,具体为:根据所述飞行器能量-阻力加速度走廊上界和走廊下界,以飞行器能量为横坐标,以飞行器阻力加速度为纵坐标,生成飞行器能量阻力加速度特征剖面;
A)当E≥Ec时,所述飞行器能量阻力加速度特征剖面包括n段二次曲线,n为正整数;其中,第i段二次曲线对应的飞行器能量阻力加速度特征剖面,具体为:
Di0=Ci1E2+Ci2E+Ci3,Eif<E≤Eib;
其中,i=1,2,3,...,n,Eif、Eib分别为第i条阻力加速度二次曲线的终端能量和起始能量,Ci1、Ci2和Ci3为第i条阻力加速度二次曲线的系数,Ec为二次曲线和一次曲线分界点能量,E为飞行器能量;Ci1、Ci2和Ci3由Eif对应的阻力加速度Dif、Eib对应的阻力加速度Dib以及Eim(Eif<Eim≤Eib)对应的阻力加速度Dim,这三个特征点确定。利用该阻力加速度曲线,可计算出对应的飞行航程为:
B)当E<Ec时,所述飞行器能量阻力加速度特征剖面为一次曲线,具体为;
其中,Dc为Ec点对应的飞行器阻力加速度,Df为Ef点对应的飞行器阻力加速度,Ef为再入终端的标称能量,vf为再入终端标称速度,hf为再入终端标称高度,gf为再入终端标称引力加速度,ρf为再入终端标称大气密度。利用该阻力加速度曲线,可计算出对应的飞行航程为:
步骤4)修正飞行器能量阻力加速度特征剖面的方法,具体为:
41)根据所述飞行器能量阻力加速度特征剖面,确定飞行器能量阻力加速度特征剖面对应的飞行航程sD0;
42)确定飞行器初始待飞航程st0;
43)修正所述步骤3)生成的飞行器能量阻力加速度特征剖面,使得所述步骤41)确定的sD0与所述步骤42)确定的st0一致,获得修正后的所述飞行器能量阻力加速度特征剖面。
步骤41)飞行器能量阻力加速度特征剖面对应的飞行航程sD0,具体为:
步骤42)飞行器初始待飞航程st0,具体为:
st0=r0arccos(sinφ0sinφf+cosφ0cosφfcos(λf-λ0)),
其中,φf为再入终端标称纬度,λf为再入终端标称经度,φ0为弹道倾角变化为零时飞行器所处的纬度,λ0为弹道倾角变化为零时飞行器所处的经度,r0为弹道倾角变化为零时飞行器的地心距。
步骤43)修正所述步骤3)生成的飞行器能量阻力加速度特征剖面,使得所述步骤41)确定的sD0与所述步骤42)确定的st0一致,获得修正后的所述飞行器能量阻力加速度特征剖面,具体为:
431)通过调整各条二次系数中Eif对应的阻力加速度Dif、Eib对应的阻力加速度Dib以及Eim(Eif<Eim≤Eib)对应的阻力加速度Dim,最终使得生成的飞行器能量阻力加速度特征剖面对应的飞行航程sD0与待飞航程st0一致,即sD0=st0。同时,应保证选取的特征点拟合出的阻力加速度曲线之间具有连续性;
432)得到航程匹配对应的Eif对应的阻力加速度Dif、Eib对应的阻力加速度Dib以及Eim(Eif<Eim≤Eib)对应的阻力加速度Dim,根据上述系数确定,Ci1、Ci2和Ci3为第i条阻力加速度二次曲线的系数。得到特征剖面的最终结果。
步骤5)根据飞行器能量阻力加速度特征剖面的修正量,解析构造再入制导模型,具体为:
51)当前飞行器能量阻力加速度特征剖面和再实时调整dDDC后的飞行器能量阻力加速度特征剖面公式为:
511)当E≥Ec时,当前阻力加速度剖面中第i条阻力加速度曲线为:
Di=Di0+ΔD,Eif<E≤Eib;
其中,ΔD为飞行器能量阻力加速度特征剖面的修正量,当前飞行器能量阻力加速度特征剖面再实时调整dDDC后的第i条阻力加速度剖面为:
DiC=Di0+ΔD+dDDC,Eif<E≤Eib;
512)当E<Ec时,当前阻力加速度剖面为:
其中,Ec为预先设定的二次曲线和一次曲线分界点能量,En为当前能量,当En≤Ef时,Dl=Df;当Ef<En≤Ec时,EN=En;En>Ec时,EN=Ec。ΔD为EN处飞行器能量阻力加速度特征剖面的修正量。在EN处再实时调整dDDC,并保持终端能量点处阻力加速度不变,实时调整后的一次阻力加速度曲线为:
513)当前阻力加速度剖面实时修正dDDC对航程偏差的影响为:
其中,sDt=sD-st,sD为当前阻力加速度对应的飞行航程,st为当前待飞航程,为导航***提供的阻力加速度变化率。
st=rarccos(sinφsinφf+cosφcosφfcos(λf-λ)),
其中,r为当前飞行器地心距;φ为当前飞行器纬度;λ为当前飞行器经度,当飞行器当前能量En≤Eif时,EiN=Eif;当Eif<En≤Eib时,EiN=En;En>Eib时,EiN=Eib。
52)根据实时调整阻力加速度剖面对飞行器航程偏差的影响,解析构造再入纵向制导模型,具体方法为:
sDt=sD-st,
st=rarccos(sinφsinφf+cosφcosφfcos(λf-λ)),
其中,θ为飞行器弹道倾角,u=cosσ,D0为当前能量对应的初始阻力加速度剖面下的阻力加速度,D为飞行器阻力加速度,ΔD为飞行器能量阻力加速度特征剖面的修正量,L为飞行器升力加速度,hs为大气密度特征参数,DN=Dl0+ΔD,r为当前飞行器地心距;φ为当前飞行器纬度;λ为当前飞行器经度;
EN确定方法如下:
a)当E>Ec时:
EN=Ec,
b)当Ef<En≤Ec时:
EN=En,
步骤6),确定飞行器倾侧角指令具体为:
σ=σabssign(σ),
σabs=arccos(u),
其中,a0>0,a1>0,a2>0,ε0>0,其中G-1为G的逆矩阵,sign()为符号函数,σi-1为上一制导周期的倾侧角指令,Δψd为侧向方位边界,Δψ为侧向方位。
实施例
初始条件为:初始高度120000m,初始速度7511.4m/s,初始弹道倾角-1.1°,初始航迹方位角46.3°,初始经度15.589°,初始纬度8.370°。
终端条件为:终端高度26740m,终端经度85.668075°,终端纬度42.194275°,终端地速765.8m/s,终端弹道倾角-8.5°,终端航迹方位角0°。
过程约束:最大动压约束12000Pa,最大过载约束2.5g,驻点热流密度约束800kW/m2。图5为本实施例全程高度曲线,图6为本实施例全程地轨迹曲线,图7为本实施例全程弹道倾角曲线,图8为本实施例的全程的航迹偏角曲线,图9为本实施例的全程的阻力加速度曲线。
本发明的具体步骤如下:
步骤1)飞行器能量-阻力加速度走廊的确定,具体为:
其中,qmax=12000Pa,nmax=2.5g,飞行器能量-阻力加速度走廊具体如图2所示。
步骤3)生成飞行器能量阻力加速度特征剖面简称阻力加速度剖面,具体方法为:
31)当E≥Ec(Ec=4.4×106)时,为二次曲线段,由多条二次曲线拼接而成,具体为:
利用该阻力加速度曲线,可计算出对应的飞行航程为:
32)当E<Ec时,为一次曲线段,由一条一次曲线组成,具体为:
其中,Dc=12.45,Df=7.3081,Ef=5.5306×105。利用该阻力加速度曲线,可计算出对应的飞行航程为:
sf0=3.6710×105,
步骤41)利用如下公式计算该阻力加速度曲线下的飞行航程:
sD0=3.47539×106,
步骤42)利用如下公式计算初始待飞航程st0:
st0=r0arccos(sinφ0sinφf+cosφ0cosφfcos(λf-λ0)),
其中:φf=42.194275°,λf=85.668075°,φ0=24.858637°,λ0=37.355623°,r0=6.457264956×106。
步骤43)通过调整各条二次系数中Eif对应的阻力加速度Dif、Eib对应的阻力加速度Dib以及Eim(Eif<Eim≤Eib)对应的阻力加速度Dim,最终使得生成的飞行器能量阻力加速度特征剖面对应的飞行航程sD0与待飞航程st0一致,即sD0=st0。同时,应保证选取的特征点拟合出的阻力加速度曲线之间具有连续性。
最终选取的再入阻力加速度特征点及飞行器能量阻力加速度特征剖面具体如图2所示。
所述步骤4)、步骤5)根据当前飞行器能量阻力加速度特征剖面实时修正对航程偏差的影响,实时调整飞行器能量阻力加速度特征剖面,解析构造再入纵向制导模型,具体方法为:
51)根据当前飞行器能量阻力加速度特征剖面对应的飞行航程与当前待飞航程的偏差的情况,实时调整飞行器能量阻力加速度特征剖面,并解算出对飞行器航程偏差的影响。当前飞行器能量阻力加速度特征剖面和再实时调整dDDC后的飞行器能量阻力加速度特征剖面公式为:
511)当E≥Ec时,当前阻力加速度剖面中第i条阻力加速度曲线为:
Di=Di0+ΔD,Eif<E≤Eib;
其中,ΔD为当前阻力加速度剖面相对于阻力加速度特征剖面的已修正量。当前飞行器能量阻力加速度特征剖面再实时调整dDDC后的第i条阻力加速度剖面为:
DiC=Di0+ΔD+dDDC,Eif<E≤Eib;
512)当E<Ec时,当前阻力加速度剖面为:
其中,当En≤Ef(En为当前能量)时,Dl=Df;当Ef<En≤Ec时,EN=En;En>Ec时,EN=Ec。ΔD为EN处当前阻力加速度相对于阻力加速度特征剖面的已修正量。在EN处再实时调整dDDC,并保持终端能量点处阻力加速度不变,实时调整后的一次阻力加速度曲线为:
当前阻力加速度剖面实时修正dDDC对航程偏差的影响为:
其中,sDt=sD-st,sD为当前阻力加速度对应的飞行航程,st为当前待飞航程。
st=rarccos(sinφsinφf+cosφcosφfcos(λf-λ)),
其中,r为当前飞行器地心距,φ为当前飞行器纬度,λ为当前飞行器经度,sin为能量大于等于Ec部分中的第i条阻力加速度剖面对应的飞行航程,具体为:
其中,当飞行器当前能量En≤Eif时,EiN=Eif;当Eif<En≤Eib时,EiN=En;En>Eib时,EiN=Eib。
sfn为能量小于Ec部分对应的一次曲线阻力加速度剖面对应的飞行航程,具体为:
根据实时调整阻力加速度剖面对飞行器航程偏差的影响,解析构造再入纵向制导模型,具体方法为:
根据实时调整阻力加速度剖面对飞行器航程偏差的影响,得出解析构造的再入制导模型为:
其中,G=μDG0,θ为飞行器弹道倾角, u=cosσ,D0为当前能量对应的初始阻力加速度剖面下的阻力加速度, 其中,DN=Dl0+ΔD,飞行器速度v曲线如图4所示。
而的解析求解公式分别为:
a)能量大于等于Ec部分中的第i条阻力加速度剖面曲线:
b)能量小于Ec部分:
步骤6)根据解析构造的再入纵向制导模型,确定飞行器纵向制导律,并解算出倾侧角指令大小,具体方法为:
61)利用解析构造的再入制导模型,得到的再入纵向制导律为:
其中,a0=0.0001,a1=0.1,a2=0.01,ε0=1。
62)根据u,计算倾侧角指令大小计算公式为:
σabs=arccos(u),
所述步骤6)使用侧向制导律,根据侧向方位误差,实时确定倾侧角指令的正负值,具体方法为:
其中,sign()为符号函数,σi-1为上一制导周期的倾侧角指令,Δψd为侧向方位边界具体如图3所示,Δψ为侧向方位。
所述步骤6)最后根据倾侧角指令大小和倾侧角指令的正负,得出倾侧角指令结果如图8所示,具体方法为:
σ=σabssign(σ)。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (8)
1.一种实时解析构造的空天飞行器再入制导方法,其特征在于,包括步骤如下:
1)判断飞行器弹道倾角是否为零,若飞行器弹道倾角不为零,则进入步骤2),若飞行器弹道倾角为零,则进入步骤3);
2)将保持飞行器弹道倾角为零作为制导信息,将所述制导信息传递给飞行器控制***,直至飞行器弹道倾角零后,进入步骤3);
3)确定飞行器能量-阻力加速度走廊上界和走廊下界,根据所述飞行器能量-阻力加速度走廊上界和走廊下界,生成飞行器能量阻力加速度特征剖面;
4)修正飞行器能量阻力加速度特征剖面,获得飞行器能量阻力加速度特征剖面的修正量;
5)根据飞行器能量阻力加速度特征剖面的修正量,解析构造再入制导模型;
6)根据所述解析构造的再入制导模型,确定飞行器倾侧角指令;
7)将所述确定的倾侧角指令作为制导信息发送给飞行器控制***,完成飞行器再入制导工作。
2.根据权利要求1所述的一种实时解析构造的空天飞行器再入制导方法,其特征在于,所述步骤3)确定飞行器能量-阻力加速度走廊上界根据下列等式确定,具体为:
所述步骤3)飞行器能量-阻力加速度走廊下界根据下列等式确定,具体为:
其中,qmax为最大动压约束,nmax为最大过载约束,为最大热流密度约束,Dqmax为动压约束下的飞行器能量-阻力加速度走廊上界,Dnmax为过载约束下的飞行器能量-阻力加速度走廊上界,为热流密度约束约束下的飞行器能量-阻力加速度走廊上界,DEQmin为平衡滑翔约束下的飞行器能量-阻力加速度走廊下界,CL为气动升力系数,kQ为与飞行器外形相关的常值参数,r为飞行器地心距,Sref为飞行器参考面积,m为飞行器质量,CD为气动阻力系数,kQ为与飞行器外形相关的常值参数,v为飞行器速度,g为引力加速度。
3.根据权利要求2所述的一种实时解析构造的空天飞行器再入制导方法,其特征在于,所述步骤3)生成的飞行器能量阻力加速度特征剖面,具体为:
根据所述飞行器能量-阻力加速度走廊上界和走廊下界,以飞行器能量为横坐标,以飞行器阻力加速度为纵坐标,生成飞行器能量阻力加速度特征剖面;
A)当E≥Ec时,所述飞行器能量阻力加速度特征剖面包括n段二次曲线,n为正整数;其中,第i段二次曲线对应的飞行器能量阻力加速度特征剖面,具体为:
Di0=Ci1E2+Ci2E+Ci3,Eif<E≤Eib;
其中,i=1,2,3,...,n,Eif、Eib分别为第i条阻力加速度二次曲线的终端能量和起始能量,Ci1、Ci2和Ci3为第i条阻力加速度二次曲线的系数,Ec为二次曲线和一次曲线分界点能量,E为飞行器能量;
B)当E<Ec时,所述飞行器能量阻力加速度特征剖面为一次曲线,具体为;
其中,Dc为Ec点对应的飞行器阻力加速度,Df为Ef点对应的飞行器阻力加速度,Ef为再入终端的标称能量,vf为再入终端标称速度,hf为再入终端标称高度,gf为再入终端标称引力加速度,ρf为再入终端标称大气密度。
4.根据权利要求3所述的一种实时解析构造的空天飞行器再入制导方法,其特征在于,所述步骤4)修正飞行器能量阻力加速度特征剖面的方法,具体为:
41)根据所述飞行器能量阻力加速度特征剖面,确定飞行器能量阻力加速度特征剖面对应的飞行航程sD0;
42)确定飞行器初始待飞航程st0;
43)修正所述步骤3)生成的飞行器能量阻力加速度特征剖面,使得所述步骤41)确定的sD0与所述步骤42)确定的st0一致,获得修正后的所述飞行器能量阻力加速度特征剖面。
5.根据权利要求4所述的一种实时解析构造的空天飞行器再入制导方法,其特征在于,所述步骤41)飞行器能量阻力加速度特征剖面对应的飞行航程sD0,具体为:
6.根据权利要求4或5之一所述的一种实时解析构造的空天飞行器再入制导方法,其特征在于,所述步骤42)飞行器初始待飞航程st0,具体为:
st0=r0arccos(sinφ0sinφf+cosφ0cosφfcos(λf-λ0)),
其中,φf为再入终端标称纬度,λf为再入终端标称经度,φ0为弹道倾角变化为零时飞行器所处的纬度,λ0为弹道倾角变化为零时飞行器所处的经度,r0为弹道倾角变化为零时飞行器的地心距。
7.根据权利要求6所述的一种实时解析构造的空天飞行器再入制导方法,其特征在于,所述步骤5)解析构造再入制导模型,具体为:
sDt=sD-st,
st=rarccos(sinφsinφf+cosφcosφfcos(λf-λ)),
其中,θ为飞行器弹道倾角,u=cosσ,D0为当前能量对应的初始阻力加速度剖面下的阻力加速度,D为飞行器阻力加速度,ΔD为飞行器能量阻力加速度特征剖面的修正量,L为飞行器升力加速度,hs为大气密度特征参数,DN=Dl0+ΔD,r为当前飞行器地心距;φ为当前飞行器纬度;λ为当前飞行器经度;
EN确定方法如下:
a)当E>Ec时:
EN=Ec,
b)当Ef<En≤Ec时:
EN=En,
8.根据权利要求7所述的一种实时解析构造的空天飞行器再入制导方法,其特征在于,所述步骤6),确定飞行器倾侧角指令具体为:
σ=σabssign(σ),
σabs=arccos(u),
其中,a0>0,a1>0,a2>0,ε0>0,其中G-1为G的逆矩阵,sign()为符号函数,σi-1为上一制导周期的倾侧角指令,Δψd为侧向方位边界,Δψ为侧向方位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811302375.8A CN109240323B (zh) | 2018-11-02 | 2018-11-02 | 一种实时解析构造的空天飞行器再入制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811302375.8A CN109240323B (zh) | 2018-11-02 | 2018-11-02 | 一种实时解析构造的空天飞行器再入制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109240323A true CN109240323A (zh) | 2019-01-18 |
CN109240323B CN109240323B (zh) | 2021-07-13 |
Family
ID=65076558
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811302375.8A Active CN109240323B (zh) | 2018-11-02 | 2018-11-02 | 一种实时解析构造的空天飞行器再入制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109240323B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110232215A (zh) * | 2019-05-09 | 2019-09-13 | 中国人民解放军国防科技大学 | 考虑机动任务需求的三维剖面分层迭代规划方法、***及介质 |
CN110262536A (zh) * | 2019-06-25 | 2019-09-20 | 湖北航天技术研究院总体设计所 | 无动力飞行器的纵向控制飞行能量管理方法及*** |
CN110309590A (zh) * | 2019-06-28 | 2019-10-08 | 北京理工大学 | 一种再入飞行器速度-高度再入走廊预测方法 |
CN111813146A (zh) * | 2020-07-01 | 2020-10-23 | 大连理工大学 | 基于bp神经网络预测航程的再入预测-校正制导方法 |
CN112256061A (zh) * | 2020-10-30 | 2021-01-22 | 北京航空航天大学 | 复杂环境及任务约束下的高超声速飞行器再入制导方法 |
CN113110558A (zh) * | 2021-05-12 | 2021-07-13 | 南京航空航天大学 | 一种混合推进无人机需求功率预测方法 |
CN116384600A (zh) * | 2023-06-06 | 2023-07-04 | 中国空气动力研究与发展中心超高速空气动力研究所 | 基于能量分析的航天器leo椭圆轨道衰降过程参数预报方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102139768A (zh) * | 2010-10-28 | 2011-08-03 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行的攻角制导方法 |
CN103486905A (zh) * | 2013-09-06 | 2014-01-01 | 中国运载火箭技术研究院 | 一种再入飞行器末制导交班条件确定方法 |
CN104656450A (zh) * | 2015-01-20 | 2015-05-27 | 北京航空航天大学 | 一种高超声速飞行器平稳滑翔再入弹道设计方法 |
US20160299233A1 (en) * | 2012-12-19 | 2016-10-13 | Elwha LLC, a limited liability corporation of the State of Delaware | Unoccupied flying vehicle (ufv) location assurance |
CN107941087A (zh) * | 2017-10-18 | 2018-04-20 | 北京航空航天大学 | 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法 |
CN107992074A (zh) * | 2017-12-07 | 2018-05-04 | 大连理工大学 | 一种基于飞行路径角规划的再入轨迹设计方法 |
-
2018
- 2018-11-02 CN CN201811302375.8A patent/CN109240323B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102139768A (zh) * | 2010-10-28 | 2011-08-03 | 中国科学院力学研究所 | 一种亚轨道飞行器再入飞行的攻角制导方法 |
US20160299233A1 (en) * | 2012-12-19 | 2016-10-13 | Elwha LLC, a limited liability corporation of the State of Delaware | Unoccupied flying vehicle (ufv) location assurance |
CN103486905A (zh) * | 2013-09-06 | 2014-01-01 | 中国运载火箭技术研究院 | 一种再入飞行器末制导交班条件确定方法 |
CN104656450A (zh) * | 2015-01-20 | 2015-05-27 | 北京航空航天大学 | 一种高超声速飞行器平稳滑翔再入弹道设计方法 |
CN107941087A (zh) * | 2017-10-18 | 2018-04-20 | 北京航空航天大学 | 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法 |
CN107992074A (zh) * | 2017-12-07 | 2018-05-04 | 大连理工大学 | 一种基于飞行路径角规划的再入轨迹设计方法 |
Non-Patent Citations (3)
Title |
---|
严晗 等: "有限时间收敛的阻力加速度跟踪进入制导律", 《空间控制技术与应用》 * |
唐军刚: "天基再入飞行器全程机动制导方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
王涛 等: "考虑阻力加速度的再入预测-校正制导算法", 《宇航学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110232215A (zh) * | 2019-05-09 | 2019-09-13 | 中国人民解放军国防科技大学 | 考虑机动任务需求的三维剖面分层迭代规划方法、***及介质 |
CN110232215B (zh) * | 2019-05-09 | 2023-03-03 | 中国人民解放军国防科技大学 | 考虑机动任务需求的三维剖面分层迭代规划方法、***及介质 |
CN110262536A (zh) * | 2019-06-25 | 2019-09-20 | 湖北航天技术研究院总体设计所 | 无动力飞行器的纵向控制飞行能量管理方法及*** |
CN110262536B (zh) * | 2019-06-25 | 2022-03-01 | 湖北航天技术研究院总体设计所 | 无动力飞行器的纵向控制飞行能量管理方法及*** |
CN110309590A (zh) * | 2019-06-28 | 2019-10-08 | 北京理工大学 | 一种再入飞行器速度-高度再入走廊预测方法 |
CN110309590B (zh) * | 2019-06-28 | 2021-01-19 | 北京理工大学 | 一种再入飞行器速度-高度再入走廊预测方法 |
CN111813146A (zh) * | 2020-07-01 | 2020-10-23 | 大连理工大学 | 基于bp神经网络预测航程的再入预测-校正制导方法 |
CN112256061A (zh) * | 2020-10-30 | 2021-01-22 | 北京航空航天大学 | 复杂环境及任务约束下的高超声速飞行器再入制导方法 |
CN113110558A (zh) * | 2021-05-12 | 2021-07-13 | 南京航空航天大学 | 一种混合推进无人机需求功率预测方法 |
CN116384600A (zh) * | 2023-06-06 | 2023-07-04 | 中国空气动力研究与发展中心超高速空气动力研究所 | 基于能量分析的航天器leo椭圆轨道衰降过程参数预报方法 |
CN116384600B (zh) * | 2023-06-06 | 2023-08-22 | 中国空气动力研究与发展中心超高速空气动力研究所 | 基于能量分析的航天器leo椭圆轨道衰降过程参数预报方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109240323B (zh) | 2021-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109240323A (zh) | 一种实时解析构造的空天飞行器再入制导方法 | |
CN110471450B (zh) | 在高度速度剖面内直接规划再入轨迹的方法 | |
CN110908396B (zh) | 可重复使用运载器全阶段再入返回制导方法 | |
CN111306989B (zh) | 一种基于平稳滑翔弹道解析解的高超声速再入制导方法 | |
CN104035335A (zh) | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导律 | |
Yu et al. | Analytical entry guidance for coordinated flight with multiple no-fly-zone constraints | |
CN106371312B (zh) | 基于模糊控制器的升力式再入预测-校正制导方法 | |
CN109240335B (zh) | 一种空天飞行器进场着陆制导方法 | |
CN101893892B (zh) | 一种无人机自动伞降回收控制方法 | |
CN107844128B (zh) | 一种基于复合比例导引的高超声速飞行器巡航段制导方法 | |
CN110015446B (zh) | 一种半解析的火星进入制导方法 | |
Slegers et al. | Terminal guidance of autonomous parafoils in high wind-to-airspeed ratios | |
US9684309B2 (en) | Aircraft guidance based on partial differential equation for miss distance | |
CN106444822A (zh) | 一种基于空间矢量场制导的平流层飞艇路径跟踪控制方法 | |
CN106249592A (zh) | 用于航空器的优化下降和进场剖面的自动确定方法和*** | |
De Ridder et al. | Optimal longitudinal trajectories for reusable space vehicles in the terminal area | |
CN108534785B (zh) | 一种大气进入制导轨迹自适应规划方法 | |
Scheuermann et al. | Combined lateral and longitudinal control of parafoils using upper-surface canopy spoilers | |
Gong et al. | Using flight manual data to derive aero-propulsive models for predicting aircraft trajectories | |
CN113835442A (zh) | 高超声速滑翔飞行器线性伪谱再入制导方法和*** | |
CN106444793A (zh) | 一种基于标称速度补偿思想的rlv进场着陆段速度控制方法 | |
CN114442673B (zh) | 一种基于混沌模型与dnn的飞行器智能轨迹规划方法 | |
Wang et al. | Entry guidance command generation for hypersonic glide vehicles under threats and multiple constraints | |
Miyazawa et al. | Flight control system for the automatic landing flight experiment | |
CN106444832B (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 |