CN103838914B - 一种高超声速飞行器滑翔段弹道解析求解方法 - Google Patents
一种高超声速飞行器滑翔段弹道解析求解方法 Download PDFInfo
- Publication number
- CN103838914B CN103838914B CN201310744002.7A CN201310744002A CN103838914B CN 103838914 B CN103838914 B CN 103838914B CN 201310744002 A CN201310744002 A CN 201310744002A CN 103838914 B CN103838914 B CN 103838914B
- Authority
- CN
- China
- Prior art keywords
- cos
- gamma
- sin
- psi
- phi
- 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
Links
Landscapes
- Navigation (AREA)
Abstract
本发明公开了一种高超声速飞行器滑翔段弹道解析求解方法,包括以下几个步骤:步骤1:定义当地地心坐标系;步骤2:获取在当地地心坐标系下的运动方程;步骤3:获取滑翔弹道的解析解。本发明提出了升力系数坡度率调整量的概念,使得动力学方程能够将速度项分离,获得控制量与几何弹道的直接关系;本发明提出了再入弹道分段线性化的策略,通过引入当地地心坐标系,使得再入横向弹道能够线性化。
Description
技术领域
本发明涉及一种高超声速飞行器滑翔段弹道解析求解方法,属于航天技术、武器技术领域。
背景技术
高超声速飞行器是指飞行马赫数大于或者等于5的飞行器,它具有飞行速度快、突防能力强、全球到达、毁伤威力大等独特的优势,已经成为当今世界各国武器研制的热点和焦点。再入弹道规划则是高超声速飞行器的一项重要的关键技术,对于典型的高超声速飞行器再入弹道,滑翔段占据航程的绝大部分,因此,滑翔段的弹道规划就显得尤为重要。
在滑翔段的飞行过程中,飞行的速度和空域远大于常规导弹,将经历高度、马赫数的大范围变化,并且受到热流密度、动压、过载等诸多约束,使得滑翔段弹道规划成为一个高敏度问题。当前滑翔段弹道规划多采用数值优化算法来获得,如物理规划、高斯伪谱法、序列二次规划法、遗传算法+序列二次规划法、Legendre伪谱法等。然而这些方法都需要进行大量的计算,其中很大一部分是关于弹道积分的。因此,若能获得滑翔段弹道的解析解,将能够大大提高弹道规划的效率,同时,它还能够反应滑翔段飞行物理特性。因此滑翔段弹道解析解具有很高的理论意义和实际价值。
发明内容
本发明的目的是为了解决上述问题,提出一种高超声速飞行器滑翔段弹道解析求解方法,获得高度、弹道倾角、弹道偏角、经纬度、速度、攻角和倾侧角的解析解,为快速弹道规划提供理论支持。
一种高超声速飞行器滑翔段弹道解析求解方法,包括以下几个步骤:
步骤1:定义当地地心坐标系;
步骤2:获取在当地地心坐标系下的运动方程;
步骤3:获取滑翔弹道的解析解。
本发明的优点在于:
(1)提出了升力系数坡度率调整量的概念,使得动力学方程能够将速度项分离,获得控制量与几何弹道的直接关系;
(2)提出了再入弹道分段线性化的策略,通过引入当地地心坐标系,使得再入横向弹道能够线性化;
(3)提出了以CN1和CY为控制变量的再入弹道规划策略,使得横向和纵向弹道能够完全解耦;
(4)提出了考虑纵向弹道机动的能量-射程映射关系,从而获得更精确的速度解析解。
(5)在进行弹道规划时不需要进行弹道积分,大大加快规划运算速度;
(6)利用解析弹道横纵解耦的特点可使得弹道规划时可以分开考虑横向弹道和纵向弹道。
附图说明
图1是广义赤道示意图;
图2是当地地心坐标系示意图;
图3是地心赤道旋转坐标系与当地地心坐标系的关系图;
图4是攻角曲线精度验证;
图5是倾侧角曲线精度验证;
图6是高度曲线精度验证;
图7是横向弹道精度验证;
图8是弹道偏角曲线精度验证;
图9是弹道倾角曲线精度验证;
图10是速度曲线精度验;
图11(a)是不同初始弹道倾角下的攻角曲线;
图11(b)是不同初始弹道倾角下的倾侧角曲线;
图11(c)是不同初始弹道倾角下的高度曲线;
图11(d)是不同初始弹道倾角下的横向弹道曲线;
图11(e)是不同初始弹道倾角下的弹道倾角曲线族;
图11(f)是不同初始弹道倾角下的H-V走廊曲线;
图12(a)是攻角曲线族;
图12(b)是倾侧角曲线族;
图12(c)是高度曲线族;
图12(d)是横向弹道曲线族;
图12(e)是弹道偏角曲线族;
图12(f)是弹道倾角曲线族;
图12(g)是速度曲线族;
图12(h)是过载曲线族;
图12(i)是热流密度曲线族。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明在拟平衡滑翔的基础上,提出了纵向升力系数余量的概念,在此基础上推导出了速度无关运动方程;利用上述方程推导获得了纵向升力系数余量与飞行距离、飞行高度和飞行弹道倾角之间的关系;在此基础上,采用分段策略,并引入当地地心坐标系,并利用它求解横向机动解析解;最后利用射程与能量之间的关系,获得了速度的解析解。
本发明是一种高超声速飞行器滑翔段弹道解析求解方法,包括以下几个步骤:
步骤1:定义当地地心坐标系
为了能够更准确的求解飞行器横程,引入广义赤道的概念,如图1所示,从而使得飞行器的运动在广义赤道附近。为了描述广义赤道,引入当地地心坐标系,用符号Se1表示,如图2所示,坐标原点位于地心;xe1由地心指向飞行器;ye1在由xe1与V的平面内,且与V同一侧;ze1由右手定则确定(地心赤道旋转坐标系坐标原点同样为地心,ze由地心北极,xe处在赤道平面内指向本初子午线,ye由右手定则确定)。
进一步,地心赤道旋转坐标系与当地地心坐标系的关系,如图3所示。其中:i1为轨道倾角(所谓轨道,即指广义赤道),取值范围为0°至180°;Ω1为升交点赤经,取值范围为0°至360°;ω1为近地点幅角,取值范围为0°至360°。x1、y1、y2均为过程坐标轴,用于描述坐标转换的过程。
(1)求解转换矩阵
设飞行器的当前经纬度为[θ,φ],则它在地心赤道旋转坐标系中的单位位置矢量为,
其中:θ为在地心赤道旋转坐标系下的当地经度;φ为地心赤道旋转坐标系下的纬度;为飞行器当前位置矢量。
为了在地心赤道旋转坐标系中将单位速度矢量表示出,需要进行如下坐标转换,
其中:Se表示地心赤道旋转坐标系;Su为当地铅垂坐标系;Sk为航迹坐标系。
ψ为地心赤道旋转坐标系下的行向角;γ为地心赤道旋转坐标系下的弹道倾角,最终,可以将在[θ,φ]处且弹道倾角为ψ的单位速度矢量表述为,
其中:为在地心赤道旋转坐标系下的飞行器速度矢量。
利用式(1)和式(3),可得广义赤道平面的单位法向量为:
其中:为广义赤道平面单位法向量。而地心赤道旋转坐标系z轴的单位向量为, 由轨道倾角i1的定义可知结合式(4)可得i1的取值为,
i1=arccos(cosφsinψ) (5)
利用和可进一步获得当地地心坐标系的x轴单位向量
上式在ψ=π/2,φ=0时存在歧义,此时,
由于并且升交点赤经Ω1满足利用和可获得Ω1如下所示:
由于, 并且 利用和可得,
其中:ω1为近地点幅角,judge为判断近地点幅角正负号的特征量。
最终,获得Se1与Se之间的转换关系如下所示。
其中:表示当地地心坐标系与地心赤道旋转坐标系之间的转换矩阵;Lz(ω1)表示绕z轴旋转ω1对应的转换矩阵;Lx(i1)表示绕x轴旋转i1对应的转换矩阵;Lz(Ω1)表示绕z轴旋转Ω1对应的转换矩阵。由于线性化的动力学模型存在于当地地心坐标系中,而最终获得的弹道结果需要在地心赤道旋转坐标系中表示,因此,需要来完成两者之间的转换。
(2)坐标转换
设在当地地心坐标系中存在一点[θ1,φ1],并且该点的弹道偏角为ψ1,假设该点在地心赤道旋转坐标系中的坐标为[θ,φ],弹道偏角为ψ,则它们满足,
从而,获得在地心赤道旋转坐标系中的位置为,
其中:
l11=cosω1cosi1cosΩ1-sinω1sinΩ1
l21=-sinω1cosi1cosΩ1-cosω1sinΩ1
l31=sini1cosΩ1
l12=cosω1cosi1sinΩ1+sinω1cosΩ1
l22=-sinω1cosi1sinΩ1+cosω1cosΩ1
l32=sini1sinΩ1
kθ=l11cosφ1cosθ1+l21cosφ1sinθ1+l31sinφ1
单位速度矢量定义如下,
另外,在地心赤道旋转坐标系中可得,
由式(14)可以解得,
从而可得,
其中:vy为单位速度矢量y轴分量;vz为单位速度矢量z轴分量。
步骤2:获取在当地地心坐标系下的运动方程
无动力滑翔飞行器单位质量机械能表述为e=μ/r-V2/2,忽略地球自转,在当地地心坐标系中三自由度运动模型可表述为,
其中:R0=6378km,θ1和φ1为当地地心坐标系下的经度和纬度;V为飞行器与地球的相对速度;ψ1为当地地心坐标系下的弹道偏角;γ为地心赤道旋转坐标系下的弹道倾角(与当地地心左边系下的弹道倾角γ1相等);h为地心赤道旋转坐标系下的高度(与当地地心坐标系下的高度相等);s为滑翔段的飞行路程;和分别为高度、当地地心坐标系下的经度、当地地心坐标系下的纬度、负比能量、弹道倾角和当地地心坐标系下的弹道偏角及射程对时间的一阶导数。D和L分别为以过载形式表示的升力和阻力,即D=0.5ρV2SCD,L=0.5ρV2SCL,其中ρ为大气密度,S为气动参考面积,m为飞行器的质量;CL和CD分别为升力系数和阻力系数;g为当地重力加速度;σ为速度滚转角。
由于φ1≈0、γ≈0、ψ1≈π/2,并假设0.5ρV2SCN0=-m(V2/r-g)cosγ且CLcosσ=CN0+CN1,其中CN0为纵向平衡滑翔升力系数、CN1为升力系数坡度率调整量,式(17)可进一步化为:
其中:CN1为升力系数坡度率调整量,CY为升力系数横向调整量,满足CY=CLsinσ。
取ψ2=π/2-ψ1,并对上式进行如下的变换,
(h+R0)dθ1=Vdt (20)
(h+R0)ds=Vdt (25)
其中:dh、dθ1、dφ1、de、dγ、dψ1和ds分别为高度、当地地心坐标系下的经度、当地地心坐标系下的纬度、负比能量、弹道倾角和当地地心坐标系下的弹道偏角及射程的微分,dt为时间的微分。
步骤3:滑翔弹道解析解
设起滑点坐标为(θ0,φ0),起滑高度为h0,起滑弹道倾角为γ0,终端高度为hf,纵向升力系数余量为CN1,横向升力系数为CY。利用式(11)可获得该点处当地地心坐标系与地心赤道旋转坐标系之间的转换关系,并且该点在当地地心坐标系中满足θ10=0、φ10=0以及ψ10=π/2。则在当地地心坐标系下滑翔段弹道的解析解如下。
(1)滑翔弹道倾角解析解
由式(19)和式(22)可得,
对式(26)进行积分可得,
其中:ρf为预测终点大气密度;γf为预测终点弹道倾角。βh为指数大气模型常数。对式(27)进一步处理可得,
从而可以得到一个对滑翔段弹道影响非常大的常数,
上述的K*即为决定滑翔段纵向弹道弯曲程度的常数。
(2)飞行路程解析解
由式(20)和式(22)可得:
将式(28)和式(29)带入式(30)可得:
设为地心平均距离。对式(31)进行积分可得:
其中:sf为预测终点航程(初始航程为0)。
(3)滑翔经度解析解
由式(20)和式(22)可得,
式(33)与式(30)比较可得,
dθ1=ds (34)
从而可得,
θ1f=sf (35)
其中:θ1f为当地地心坐标系下的预测终点经度。
(4)滑翔纬度与弹道偏角解析解
由式(21)、式(22)和式(23)可得,
将式(28)和式(29)带入式(36),可得,
由于φ10=0,并且在整段弹道中φ1为小量,故对ψ2的影响不大,在这里取,
则,式(37)可化为,
由于φ10=0,ψ20=0,对上式进行拉式变换可得,
由上式可得,
对式(41)进行拉式反变换可得,
其中:φ1f为当地地心坐标系下的预测终点纬度。由上式还可得当地地心坐标系下的预测终点弹道偏角ψ1f为ψ1f=π/2-ψ2f。利用式(42)和式(12)给出的坐标转换关系,就可以获得在地心赤道旋转坐标系下预测终点经度θf预测终点纬度φf和预测终点弹道偏角ψf。
(5)滑翔速度解析解
假定在滑翔过程中,保持定常纵向升阻比KN(KN=CLcosσ/CD),则滑翔阻力可写为,
由于Lcosσ=0.5ρV2S(CN0+CN1),而CN1为一小量,忽略后可得,
其中:D为阻力;m为飞行器质量;V为飞行器相对地球的速度;CN1为升力系数坡度率调整量;g为当地重力加速度,还可以写成如下形式,
其中:μ为由引力常数和地球质量决定的常数。
式(44)中,阻力可分为两个部分,即保持平衡滑翔的升力产生的阻力以及用于纵向转弯的升力产生的阻力,其中第一项占主要部分,第二项仅为小的修正。忽略第二部分的影响,第一部分阻力满足,
由于h相对R0来说是小量,这里用它的算术平均值代替,则上式积分可得,
其中:e0为初始负比能量。从而可得,
其中:ef为预测终点负比能量。对于另一部分用于转弯的升力,其消耗的机械能为,
其中:为弹道倾角的导数。
取 则可得,
其中:Δe为负比能量修正量,综上可得,终端速度预测为,
其中:Vf为预测终点速度;为距离地心的平均距离;KN为升阻比。需要注意的是,在式(50)的计算中利用了平均速度,因此需要通过迭代来不断修正同时,还可以利用迭代修正对KN的估计(在实际弹道中KN不为常值,为了估算的精确一般采用KN的平均值,为了简单起见,可令其中KN0为起始点处的纵向升阻比,)采取迭代修正的方式获得。
(6)滑翔攻角和倾侧角解析解
式(28)、式(35)、式(42)、式(50)分别给出了预测终点弹道倾角γf、经度θf、纬度φf、弹道偏角ψf、速度Vf与高度hf之间的关系。利用平衡滑翔关系可知,
从而可得攻角和倾侧角分别满足,
其中:α为攻角;σ为速度滚转角;Ma为马赫数。
(7)最优初始下滑角度解析解
在上面的计算中,γ0为任意给定的初值。然而由弹道阻尼理论知,合理的初始下滑角γ0将会使得整个滑翔段的控制规律更加的平稳,具体推导如下。
由拟平衡滑翔条件知,
其中:CN0为拟平衡滑翔对应的纵向升力系数。假设当γ=γ*时,飞行器能够保持平衡滑翔状态,对上式求全微分可得,
将式(18)带入上式可得,
其中:γ*为拟平衡滑翔最优下滑弹道倾角;Lcosσ为升力纵向分量。由于(Lcosσ/m)V2βh/(2g)□(Lcosσ/m)[1+mV4/(2μLcosσ)],从而上式可进一步简化为,
其中:KN为纵向升阻比。式(27)、式(32)、式(35)、式(42)、式(50)、式(52)和式(56)共同构成了滑翔段弹道的解析解。
实施例
为了检验解析求解算法的精度,选用CAV作为计算模型,进行数值仿真效验。图4至图10给出了解析弹道精度的验证结果。由图4可以看出,纵向弹道解析解的经度非常高,不管是分段多还是少,得到的结果均与数值积分完全吻合,这是由于纵向弹道求解过程中没有设定任何假设,是精确的几何弹道。同理,图5中的弹道倾角也是和数值积分完全吻合。图6可以看出,在分5段的时候解析求解的速度与数值积分得到的数度曲线还是存在微小的差异,这是由于在实际飞行过程中不可能完全做到等纵向升阻比。由于攻角和倾侧角的精度与高度精度及速度精度相关,因此图7和图8的攻角和倾侧角曲线与速度类似(高度精度高,忽略不计)。由图9可以看出,横向弹道的精度相对较低,在分5段的时候横向弹道有明显的误差,在13段的时候误差变得较小,20段的时候精度已经较高。由图10可以看出,弹道偏角的误差较小,这说明在纬度较小情况下,纬度对弹道偏角影响较小的假设是对的。图11给出了不同初始弹道倾角情况下的解析求解弹道的对比图,由图可以看出,γ0越大,则射程越远,在H-V走廊中的位置越靠上,但同时也会使得控制变量的震荡越厉害。因此,为了使规划出来的弹道控制上更易于实现,一般选取γ0=γ*,此时对应的控制震荡最小(但还是会有震荡,因为实际的γ*是随高度和速度变化的)。图12给出了基于该解析解的弹道曲线族。由图12(a)和图12(b)可以看出,以升力系数坡度率调整量和升力系数横向分量为控制变量的滑翔段规划弹道对应的攻角和倾侧角均较为平滑;由图12(c)和图12(d)可以看出,该方案的规划弹道纵程和横程覆盖范围均较大;由图12(e)可以看出,速度变化也较为光滑;由图12(f)、图12(g)和图12(h)可以看出,过载、热流密度和动压均满足要求;由图12(i)可以看出,基于该解析解的弹道族在约束走廊中的位置几乎重合。
Claims (1)
1.一种高超声速飞行器滑翔段弹道解析求解方法,包括以下几个步骤:
步骤1:定义当地地心坐标系;
当地地心坐标系Se1的坐标原点位于地心,xe1由地心指向飞行器;ye1在由xe1与V的平面内,且与V同一侧,ze1由右手定则确定;
地心赤道旋转坐标系Se的坐标原点也为地心,ze由地心北极,xe处在赤道平面内指向本初子午线,ye由右手定则确定;
(1)获取转换矩阵;
当地地心坐标系Se1与地心赤道旋转坐标系Se之间的转换矩阵为:
其中:表示当地地心坐标系与地心赤道旋转坐标系之间的转换矩阵;Lz(ω1)表示绕z轴旋转ω1对应的转换矩阵;Lx(i1)表示绕x轴旋转i1对应的转换矩阵;Lz(Ω1)表示绕z轴旋转Ω1对应的转换矩阵;ω1为近地点幅角,Ω1为升交点赤经,i1为轨道倾角,具体为:
i1=arccos(cosφsinψ) (2)
θ为在地心赤道旋转坐标系下的当地经度;φ为地心赤道旋转坐标系下的纬度;ψ为地心赤道旋转坐标系下的行向角;judge为判断近地点幅角正负号的特征量;
(2)坐标转换
设在当地地心坐标系中存在一点[θ1,φ1],并且该点的弹道偏角为ψ1,假设该点在地心赤道旋转坐标系中的坐标为[θ,φ],弹道偏角为ψ,则它们满足,
从而,获得在地心赤道旋转坐标系中的位置为,
φ=arcsin(l11cosφ1cosθ1+l21cosφ1sinθ1+l31sinφ1)
其中:
l11=cosω1cosi1cosΩ1-sinω1sinΩ1
l21=-sinω1cosi1cosΩ1-cosω1sinΩ1
l31=sini1cosΩ1
l12=cosω1cosi1sinΩ1+sinω1cosΩ1
l22=-sinω1cosi1sinΩ1+cosω1cosΩ1
l32=sini1sinΩ1
kθ=l11cosφ1cosθ1+l21cosφ1sinθ1+l31sinφ1
单位速度矢量定义如下,
另外,在地心赤道旋转坐标系中,
由式(7)得,
从而可得,
其中:vy为单位速度矢量y轴分量;vz为单位速度矢量z轴分量;
步骤2:获取在当地地心坐标系下的运动方程
无动力滑翔飞行器单位质量机械能表述为e=μ/r-V2/2,忽略地球自转,在当地地心坐标系中三自由度运动模型可表述为,
其中:R0=6378km,θ1和φ1为当地地心坐标系下的经度和纬度;V为飞行器与地球的相对速度;ψ1为当地地心坐标系下的弹道偏角;γ为地心赤道旋转坐标系下的弹道倾角;h为地心赤道旋转坐标系下的高度;s为滑翔段的飞行路程;为高度对时间的一阶导数;为当地地心坐标系下的经度对时间的一阶导数;为当地地心坐标系下的纬度对时间的一阶导数;为负比能量对时间的一阶导数;为弹道倾角对时间的一阶导数;为当地地心坐标系下的弹道偏角对时间的一阶导数;为射程对时间的一阶导数;D和L分别为以过载形式表示的升力和阻力,即D=0.5ρV2SCD,L=0.5ρV2SCL,其中ρ为大气密度,S为气动参考面积,m为飞行器的质量;CL和CD分别为升力系数和阻力系数;g为当地重力加速度;σ为速度滚转角;
由于φ1≈0、γ≈0、ψ1≈π/2,并假设0.5ρV2SCN0=-m(V2/r-g)cosγ且CLcosσ=CN0+CN1,其中CN0为纵向平衡滑翔升力系数、CN1为升力系数坡度率调整量,式(10)可进一步化为:
其中:CN1为升力系数坡度率调整量,CY为升力系数横向调整量,满足CY=CLsinσ;
取ψ2=π/2-ψ1,并对上式进行如下的变换,
(h+R0)dθ1=Vdt (13)
(h+R0)ds=Vdt (18)
其中:dh、dθ1、dφ1、de、dγ、dψ1和ds分别为高度、当地地心坐标系下的经度、当地地心坐标系下的纬度、负比能量、弹道倾角和当地地心坐标系下的弹道偏角及射程的微分,dt为时间的微分;
步骤3:获取滑翔弹道的解析解
设起滑点坐标为(θ0,φ0),起滑高度为h0,起滑弹道倾角为γ0,终端高度为hf,纵向升力系数余量为CN1,横向升力系数为CY;利用式(1)可获得该点处当地地心坐标系与地心赤道旋转坐标系之间的转换关系,并且该点在当地地心坐标系中满足θ10=0、φ10=0以及ψ10=π/2;则在当地地心坐标系下滑翔段弹道的解析解如下;
(1)滑翔弹道倾角解析解
由式(12)和式(15)可得,
对式(19)进行积分可得,
其中:ρf为预测终点大气密度;γf为预测终点弹道倾角;βh为指数大气模型常数;对式(20)进一步处理可得,
从而得到常数,
上述的K*即为决定滑翔段纵向弹道弯曲程度的常数;
(2)飞行路程解析解
由式(15)和式(18)可得:
将式(21)和式(22)带入式(23)可得:
设为地心平均距离;对式(24)进行积分可得:
其中:sf为预测终点航程;
(3)滑翔经度解析解
由式(13)和式(15)可得,
式(26)与式(23)比较可得,
dθ1=ds (27)
从而可得,
θ1f=sf (28)
其中:θ1f为当地地心坐标系下的预测终点经度;
(4)滑翔纬度与弹道偏角解析解
由式(14)、式(15)和式(16)可得,
将式(21)和式(22)带入式(29),可得,
取,
则,式(30)可化为,
由于φ10=0,ψ20=0,对上式进行拉式变换可得,
sφ1(s)=K1ψ2(s)
由上式可得,
对式(34)进行拉式反变换可得,
其中:φ1f为当地地心坐标系下的预测终点纬度;由上式还可得当地地心坐标系下的预测终点弹道偏角ψ1f为ψ1f=π/2-ψ2f;利用式(35)和式(1)给出的坐标转换关系,就可以获得在地心赤道旋转坐标系下预测终点经度θf预测终点纬度φf和预测终点弹道偏角ψf;
(5)滑翔速度解析解
假定在滑翔过程中,保持定常纵向升阻比KN,KN=CLcosσ/CD,则滑翔阻力为,
由于Lcosσ=0.5ρV2S(CN0+CN1),忽略CN1后可得,
其中:D为阻力;m为飞行器质量;V为飞行器相对地球的速度;CN1为升力系数坡度率调整量;g为当地重力加速度,写成如下形式,
其中:μ为由地球引力常数;
式(37)中,阻力可分为两个部分,即保持平衡滑翔的升力产生的阻力以及用于纵向转弯的升力产生的阻力,忽略第二部分的影响,第一部分阻力满足,
由于h相对R0为小量,采用它的算术平均值代替,则上式积分可得,
其中:e0为初始负比能量;从而可得,
其中:ef为预测终点负比能量;对于另一部分用于转弯的升力,其消耗的机械能为,
其中:为弹道倾角的导数;
取则可得,
其中:Δe为负比能量修正量,综上可得,终端速度预测为:
其中:Vf为预测终点速度;为距离地心的平均距离;KN为升阻比;
在式(43)的计算中利用了平均速度,需要通过迭代来不断修正同时,利用迭代修正对KN的估计;
(6)滑翔攻角和倾侧角解析解
式(21)、式(28)、式(35)、式(43)分别给出了预测终点弹道倾角γf、经度θf、纬度φf、弹道偏角ψf、速度Vf与高度hf之间的关系;利用平衡滑翔关系可知,
从而可得攻角和倾侧角分别满足,
其中:α为攻角;σ为速度滚转角;Ma为马赫数;
(7)最优初始下滑角度解析解
具体推导如下;
由拟平衡滑翔条件知,
其中:CN0为拟平衡滑翔对应的纵向升力系数;假设当γ=γ*时,飞行器能够保持平衡滑翔状态,对上式求全微分可得,
将式(11)带入上式可得,
其中:γ*为拟平衡滑翔最优下滑弹道倾角;Lcosσ为升力纵向分量;由于(Lcosσ/m)V2βh/(2g)>>(Lcosσ/m)[1+mV4/(2μLcosσ)],从而上式可进一步简化为:
其中:KN为纵向升阻比;式(20)、式(25)、式(28)、式(35)、式(43)、式(45)和式(49)共同构成了滑翔段弹道的解析解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310744002.7A CN103838914B (zh) | 2013-12-30 | 2013-12-30 | 一种高超声速飞行器滑翔段弹道解析求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310744002.7A CN103838914B (zh) | 2013-12-30 | 2013-12-30 | 一种高超声速飞行器滑翔段弹道解析求解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103838914A CN103838914A (zh) | 2014-06-04 |
CN103838914B true CN103838914B (zh) | 2017-01-18 |
Family
ID=50802405
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310744002.7A Active CN103838914B (zh) | 2013-12-30 | 2013-12-30 | 一种高超声速飞行器滑翔段弹道解析求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103838914B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104392047B (zh) * | 2014-11-25 | 2017-05-10 | 北京航空航天大学 | 一种基于平稳滑翔弹道解析解的快速弹道规划方法 |
CN104571125B (zh) * | 2014-12-18 | 2016-01-27 | 北京控制工程研究所 | 一种利用一条标准弹道应对多种返回条件的控制方法 |
CN104850129B (zh) * | 2014-12-19 | 2017-10-24 | 北京控制工程研究所 | 一种跳跃式再入的射向预偏置横向制导方法 |
CN104656450B (zh) * | 2015-01-20 | 2017-04-19 | 北京航空航天大学 | 一种高超声速飞行器平稳滑翔再入弹道设计方法 |
CN105022858B (zh) * | 2015-05-08 | 2016-06-01 | 北京航天自动控制研究所 | 一种确定滑翔飞行器阻力加速度走廊边界的方法 |
CN106021628B (zh) * | 2015-07-03 | 2019-06-18 | 中国运载火箭技术研究院 | 一种运载火箭垂直返回弹道设计方法 |
CN105259904B (zh) * | 2015-10-15 | 2018-01-30 | 山东科技大学 | 基于模型预测控制的多操纵面无人机纵向解耦控制方法 |
CN105138808B (zh) * | 2015-10-19 | 2018-11-27 | 中国人民解放军国防科学技术大学 | 基于摄动理论的滑翔弹道误差传播分析方法 |
CN105354380B (zh) * | 2015-11-03 | 2018-12-28 | 中国人民解放军国防科学技术大学 | 面向摄动因素影响补偿的滑翔弹道快速修正方法 |
CN105550402B (zh) * | 2015-12-07 | 2018-08-31 | 北京航空航天大学 | 一种基于攻角或倾侧角变频的高超平稳机动滑翔弹道设计方法 |
CN105676638B (zh) * | 2016-01-11 | 2018-06-29 | 北京航空航天大学 | 平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法 |
CN105760573B (zh) * | 2016-01-21 | 2019-07-02 | 中国工程物理研究院总体工程研究所 | 沿临近空间大范围机动弹道的扰动引力延拓逼近方法 |
CN105718660B (zh) * | 2016-01-21 | 2019-03-01 | 中国工程物理研究院总体工程研究所 | 临近空间大范围机动弹道三维包络计算方法 |
CN105740506B (zh) * | 2016-01-21 | 2018-12-11 | 中国工程物理研究院总体工程研究所 | 沿临近空间大范围机动弹道空间包络的扰动引力逼近方法 |
CN105653827B (zh) * | 2016-03-17 | 2020-03-13 | 北京工业大学 | 高超声速飞行器Terminal滑模控制器设计方法 |
CN106227972A (zh) * | 2016-08-04 | 2016-12-14 | 北京航空航天大学 | 一种高超声速飞行器平稳滑翔弹道的优化方法 |
CN107341295B (zh) * | 2017-06-16 | 2020-09-01 | 湖北航天技术研究院总体设计所 | 具有终端角度和速度约束的下压段弹道设计方法 |
CN108549785B (zh) * | 2018-05-03 | 2021-09-24 | 中国人民解放军国防科技大学 | 一种基于三维飞行剖面的高超声速飞行器精准弹道快速预测方法 |
CN109145451B (zh) * | 2018-08-22 | 2022-11-29 | 哈尔滨工业大学 | 一种高速滑翔飞行器的运动行为识别与航迹估计方法 |
CN109190248B (zh) * | 2018-09-03 | 2023-07-18 | 中国运载火箭技术研究院 | 一种用于滑翔飞行器的滑翔射程解析方法及解析系统 |
CN111580547B (zh) * | 2020-04-15 | 2021-06-15 | 北京理工大学 | 一种高超声速飞行器编队控制方法 |
CN111580555B (zh) * | 2020-05-13 | 2022-04-08 | 北京控制工程研究所 | 一种高超声速飞行器上升段分段自适应预测校正制导方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103076015A (zh) * | 2013-01-04 | 2013-05-01 | 北京航空航天大学 | 一种基于全面最优校正的sins/cns组合导航系统及其导航方法 |
CN103090728A (zh) * | 2013-01-07 | 2013-05-08 | 北京理工大学 | 一种基于滑模控制的带末角约束制导方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6877691B2 (en) * | 2002-03-12 | 2005-04-12 | Bae Systems Information And Electronic Systems Integration Inc. | High altitude stripping for threat discrimination |
-
2013
- 2013-12-30 CN CN201310744002.7A patent/CN103838914B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103076015A (zh) * | 2013-01-04 | 2013-05-01 | 北京航空航天大学 | 一种基于全面最优校正的sins/cns组合导航系统及其导航方法 |
CN103090728A (zh) * | 2013-01-07 | 2013-05-08 | 北京理工大学 | 一种基于滑模控制的带末角约束制导方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103838914A (zh) | 2014-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103838914B (zh) | 一种高超声速飞行器滑翔段弹道解析求解方法 | |
CN104035335B (zh) | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导方法 | |
CN104392047B (zh) | 一种基于平稳滑翔弹道解析解的快速弹道规划方法 | |
CN104615144B (zh) | 基于目标规划的高超声速飞行器再入轨迹在线优化方法 | |
CN107966156A (zh) | 一种适用于运载火箭垂直回收段的制导律设计方法 | |
CN110015446B (zh) | 一种半解析的火星进入制导方法 | |
CN111306989A (zh) | 一种基于平稳滑翔弹道解析解的高超声速再入制导方法 | |
CN104215242B (zh) | 一种基于横向游移坐标系的极区惯性导航方法 | |
Zhu et al. | Impact time and angle control guidance independent of time-to-go prediction | |
CN106643341B (zh) | 基于准平衡滑翔原理的力热控制耦合设计方法 | |
Li et al. | Time-coordination entry guidance for multi-hypersonic vehicles | |
CN109657256A (zh) | 一种高精度弹道式再入标称返回轨道仿真方法 | |
CN106371312A (zh) | 基于模糊控制器的升力式再入预测‑校正制导方法 | |
CN107121929A (zh) | 基于线性协方差模型预测控制的鲁棒再入制导方法 | |
CN113900448B (zh) | 一种基于滑模干扰观测器的飞行器预测校正复合制导方法 | |
CN103453907B (zh) | 基于分层大气模型的行星进入段导航滤波方法 | |
CN104635495B (zh) | 一种含模型参数不确定性平流层卫星的经向偏移分级控制方法 | |
CN104019818A (zh) | 一种基于预测轨迹的行星导航轨道器布局优化方法 | |
CN110908407A (zh) | 一种rlv再入热流率跟踪的改进预测制导方法 | |
CN102508492B (zh) | 一种飞行器在等高航路点间的定高度大圆飞行实现方法 | |
Zhang et al. | Entry trajectory generation with complex constraints based on three-dimensional acceleration profile | |
Harlin et al. | Ballistic missile trajectory prediction using a state transition matrix | |
Yan et al. | Analysis of optimal initial glide conditions for hypersonic glide vehicles | |
CN105740506B (zh) | 沿临近空间大范围机动弹道空间包络的扰动引力逼近方法 | |
CN105354380B (zh) | 面向摄动因素影响补偿的滑翔弹道快速修正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |