CN107941087A - 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法 - Google Patents
一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法 Download PDFInfo
- Publication number
- CN107941087A CN107941087A CN201710970667.8A CN201710970667A CN107941087A CN 107941087 A CN107941087 A CN 107941087A CN 201710970667 A CN201710970667 A CN 201710970667A CN 107941087 A CN107941087 A CN 107941087A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- mtd
- msup
- 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
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F41—WEAPONS
- F41G—WEAPON SIGHTS; AIMING
- F41G3/00—Aiming or laying means
- F41G3/22—Aiming or laying means for vehicle-borne armament, e.g. on aircraft
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/242—Orbits and trajectories
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Radar, Positioning & Navigation (AREA)
- Combustion & Propulsion (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Chemical & Material Sciences (AREA)
- Secondary Cells (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法,它包括以下几个步骤:步骤一:高超声速飞行器建模;步骤二:滑翔距离解析解推导;步骤三:再入轨迹规划与跟踪制导;步骤四:基于地球自转模型的平稳滑翔弹道倾角推导。本发明的优点在于:(1)获得了精度更高的基于能量的滑翔距离解析解,形式简单,便于应用。(2)以能量为自变量规划阻力剖面,在旋转地球背景下将阻力剖面转换为纵向升阻比,然后调整倾侧角在最大升阻比下进行跟踪,横向运动由航向误差门限控制。(3)针对再入弹道的振荡问题,推导基于地球自转模型的平稳滑翔弹道倾角,并采用以此为反馈项的弹道阻尼技术抑制振荡。
Description
技术领域
本发明涉及了一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法,属于航天技术、武器技术领域。
背景技术
高超声速滑翔飞行器一般由助推火箭送入预定高度,然后无动力从轨道或亚轨道再入大气层。为了获得远程飞行能力和强横向机动能力,飞行器通常采用大升阻比的升力体结构。
再入飞行是一个多约束问题,具有任务多变、飞行时间长、速度和位置变化大、飞行环境变化剧烈等特点。这给滑翔式飞行器再入制导设计带来了困难和挑战。标准轨道制导律通过跟踪标准轨迹进行制导,其特点是计算量小,易于实现,通常以阻力加速度剖面为参考剖面,通过滑翔距离解析解求解名义弹道预测剩余航程,调整PID增益控制倾侧角跟踪参考剖面。但多数滑翔距离解析解由于忽略了地球自转引起的惯性力,在旋转地球背景下精度不高,存在跟踪误差。
发明内容
本发明的目的在于提供一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法,针对高超声速飞行器平稳滑翔上述问题,在旋转地球背景下,证明了线化地球自转引起的惯性力方法的精确性。在此基础上,推导出精度更高的滑翔距离解析解,并根据基于能量的解析解规划参考阻力剖面进行跟踪。针对再入弹道存在的长周期弱阻尼问题,采用考虑地球自转的弹道阻尼技术进行抑制。
本发明一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法,其特征在于:它包括以下几个步骤(如图12所示):
步骤一:高超声速飞行器建模
在旋转地球背景下,采用标准大气模型和平方反比引力场模型,建立高超声速飞行器的六自由度动力学方程如下:
其中,d/dt表示对时间求导,t是时间,h是海拔高度,R是射程,λ是经度,φ是纬度,V是飞行器相对于地球的速率,γ是弹道倾角,ψ是飞行器航向角,以当地北向为基准,σ是倾侧角,Re是地球半径,取值6378.137km,ωe是地球自转角速度,L=ρV2SCL/2m为升力加速度,D=ρV2SCD/2m为阻力加速度。
步骤二:滑翔距离解析解推导;
1.积分形式的射程预测公式
在旋转地球背景下,考虑由地球自转引起的惯性力,基于步骤一中的动力学模型推导积分形式的射程预测公式。
(1)基于速度的射程预测公式
由射程导数公式(2)与式(5),可得到射程对速度的导数。假设cosγ≈1,sinγ≈0,得到考虑地球自转的基于速度的射程预测公式
(2)基于能量的射程预测公式
由于飞行器在再入过程中是无动力的,能量逐渐衰减,且终端条件也可以由终端速度或终端能量定义,这里使用能量作为独立变量,其好处是不用考虑再入过程开始、结束时间的差异,且可以统一合理地描述不同轨迹段。定义能量为
E=gh+1/2V2 (9)
将上式对高度求导可得
根据式(1),(2),(5)及式(10),经过一系列代数化简,假设cosγ≈1,sinγ≈0,以及量级很小,可得基于能量的考虑地球自转的射程积分公式
2.惯性力线化及证明
由于地球自转引起的牵连加速度和科氏加速度均是时变的,且具有较强的非线性,无法直接积分,故对其进行线化处理。由于平稳滑翔段在再入飞行过程中占很大比例,主要对平稳滑翔段的惯性力进行线化处理。
在平稳滑翔段,弹道倾角很小,令式(6)为0,cosγ≈1,得到平稳滑翔段阻力满足:
上式结合式(8),式(11),考虑可近似为常数,故仅需要线化cos2φ,sinφcosφcosψ和cosφsinψ项。cos2φ与sinφcosφcosψ前面的系数(即)很小,约为10-4级别,所以在此粗略的用能量的线性函数来拟合,拟合公式如下:
sinφcosφcosψ=k11E+k10 (13)
cos2φ=k21E+k20 (14)
不用拟合cosφsinψ是因为其值近似为常数,严谨的证明过程如下:
因为cosφsinψ与速度无关,所以可忽略速度影响,视飞行器仅受重力作用做圆周运动。由于重力对ze轴力矩为0,故可得
在NED坐标系下,由
转到地心赤道旋转坐标系中可得
所以
又因为飞行器仅受重力做匀速圆周运动,且Re+h≈Re为定值,cosγ≈1,可得
cosφsinψ≈const (20)
故cosφsinψ可视为常值,代入初始值φ0和ψ0求解此常值。
3.平稳滑翔段滑翔距离解析解推导
在推导基于能量的滑翔距离解析解时,式(12)中的速度一次项难以处理,故在此将改速度项线化为能量的一次函数如下:
在此基础上,经过一系列代数变换,得到常值纵向升阻比的平稳滑翔段滑翔距离解析解如下:
(1)基于速度的滑翔距离解析解
另
若Δ≥0
若Δ<0
其中,
a2=(k11+k21)ωe 2(Re+h)2/2+1 (24)
(2)基于能量的滑翔距离解析解
其中,
b1=2+2kV1ωe(Re+h)cosφsinψ+k21ωe 2(Re+h)2+k11ωe 2(Re+h) (29)
对比分析可知,基于能量的滑翔距离解析解精度更高,形式更为简洁,故以能量为自变量进行再入轨迹规划。
步骤三:再入轨迹规划与跟踪制导;
1.基准攻角剖面规划
对于再入机动飞行器而言,射程和热防护是攻角剖面设计的主要考虑因素。针对CAV-H(High performance Common Aero Vehicle)模型,这里设计攻角曲线为能量分段函数:
其中,Es1和Es2表示攻角分段函数的分界点。在再入初始下降阶段,为尽量拉高弹道,减小热流密度,采用较大攻角α1飞行在平稳滑翔段,为了发挥最大机动能力,飞行器以最大升阻比攻角滑翔,即α2=(L/D)max。出于控制系统跟踪性和气动稳定性的考虑,设计Es1和Es2之间的线性段来保证攻角剖面平稳过渡。
2.阻力-能量剖面规划
由于初始段阻力加速度走廊狭窄,且走廊上边界可以近似为二次曲线,故将初始段的参考D(E)曲线设计成两个关于能量E的二次多项式,形式如下:
其中,A11,A12,A13,A21,A22,A23为常值。
再入走廊中部范围较大,考虑平稳滑翔条件,设计滑翔段参考D(E)曲线为一条直线。根据式(6),得:
即为关于能量E的一次曲线,比较剩余射程和预测射程大小,确定其与二次多项式阶段交点。其中,L1=Lcosσ,为升力在铅垂平面的分量,L1/D为纵向升阻比,可由步骤二中平稳滑翔段滑翔距离解析解得到。
由于再入走廊末端较为宽阔,在线性段末端规划一次曲线段,以满足过程约束和终端约束。
D=k3E+b3(E<E3) (34)
其中,k3,b3为常值。
3.参考弹道跟踪指令
通过对公式(1)和阻力公式求导,进行一系列代数化简,忽略阻力系数CD的一阶、二阶导数,得到如下一个公式,其可以将参考阻力转换为纵向升阻比L1/D。
进而通过调节基准倾侧角σbsl来跟踪纵向升阻比曲线,实现倾侧反转。
飞行器横向运动通过航向误差门限进行控制。
步骤四:基于地球自转模型的平稳滑翔弹道倾角推导;
弹道阻尼技术以平稳滑翔弹道倾角作为反馈项进行控制,可以成功抑制再入弹道存在的长周期弱阻尼振荡。在平稳滑翔段,相比于传统的PID反馈控制,弹道阻尼控制技术跟踪精确性更高,形式简单,便于应用。在旋转地球背景下,根据式(6)推导考虑地球自转的平稳滑翔弹道倾角γ'SG,满足
d2γ'SG 2+d1γ'SG+d0=0 (37)
各项系数为:
代入弹道阻尼技术中作为反馈项添加到基准攻角和倾侧角上进行控制:
αcmd=αbsl+cosσbslKγ(γ'SG-γ) (41)
用上述弹道阻尼控制技术可精确跟踪平稳滑翔弹道,有效抑制弹道振荡。
通过上述四个步骤,可使弹道稳定在平稳滑翔弹道附近,最终得到了一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法。
本发明一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法,优点在于:
(1)线化由地球自转引起的惯性力,并给出了严格证明。获得了精度更高的基于能量的滑翔距离解析解,形式简单,便于应用。
(2)以能量为自变量规划阻力剖面,在旋转地球背景下将阻力剖面转换为纵向升阻比,然后调整倾侧角在最大升阻比下进行跟踪,横向运动由航向误差门限控制。
(3)针对再入弹道的振荡问题,推导基于地球自转模型的平稳滑翔弹道倾角,并采用以此为反馈项的弹道阻尼技术抑制振荡。
附图说明
图1是地心赤道旋转坐标系(GER)和当地北东下坐标系(NED)示意图。
图2是滑翔距离解析解对比曲线。
图3(a)是再入高度走廊。
图3(b)是再入阻力加速度走廊。
图4是规划的阻力加速度剖面。
图5是航向误差门限。
图6(a)是反馈控制对比高度曲线。
图6(b)是反馈控制对比阻力剖面跟踪曲线。
图7是七种目标状态的经纬度曲线。
图8是七种目标状态高度射程曲线。
图9是七种目标状态速度时间曲线。
图10是七种目标状态弹道倾角曲线。
图11(a)是目标东南方向高度曲线。
图11(b)是目标东南方向速度时间曲线。
图11(c)是目标东南方向弹道倾角时间曲线。
图11(d)是目标东南方向倾侧角时间曲线。
图11(e)是目标东南方向航向角曲线。
图11(f)是目标东南方向航向角误差曲线。
图12为本发明方法流程图。
上述图中,涉及到的符号、代号说明如下:
图1中,表示地心赤道旋转坐标系(GER),oxyz表示当地北东下坐标系(NED)。h是海拔高度,R是射程,λ是经度,φ是纬度,V是飞行器相对于地球的速率,γ是弹道倾角,ψ是飞行器航向角,以当地北向为基准,σ是倾侧角。图1中,Trajectory Simulation表示弹道仿真结果,AGRFREE表示基于能量考虑地球自转的滑翔距离解析解,AGRFRES表示基于速度考虑地球自转的滑翔距离解析解,AGRFNREE表示基于能量无地球自转的滑翔距离解析解,AGRFNRES表示基于速度无地球自转的滑翔距离解析解。图4中,E1,E2,E3表示规划的阻力加速度剖面分段点。图5中,Without Control表示不加反馈控制,PID Control表示添加PID控制,TDCT表示弹道阻尼控制。图9中Trajectory Simulation表示弹道仿真结果,ReferenceTrajectory表示参考剖面,Entry Corridor表示高度走廊,LOS表示视线航向角。所有图例中,Height、Altitude表示高度,Range表示射程,Longitude表示经度,Latitude表示纬度,Velocity表示速度,Energy表示能量,Time表示时间,Heading Angle表示航向角,HeadingError表示航向误差,Flight-Path Angle表示弹道倾角,Downrange表示纵程,DragAcceleration表示阻力加速度。
具体实施方式
下面将结合附图和实施案例对本发明作进一步的详细说明。
针对高超声速飞行器平稳滑翔上述问题,在旋转地球背景下,证明了线化地球自转引起的惯性力方法的精确性。在此基础上,推导出精度更高的滑翔距离解析解,并根据基于能量的解析解规划参考阻力剖面进行跟踪。针对再入弹道存在的长周期弱阻尼问题,采用考虑地球自转的弹道阻尼技术进行抑制。
本发明一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法,其特征在于:它包括以下几个步骤:
步骤一:再入模型建立;
1.再入运动方程
如图1所示,这里在地心赤道旋转坐标系(GER)和当地北东下坐标系(NED)下观察飞行器运动。
其中,d/dt表示对时间求导,t是时间,h是海拔高度,R是射程,λ是经度,φ是纬度,V是飞行器相对于地球的速率,γ是弹道倾角,ψ是飞行器航向角,以当地北向为基准,σ是倾侧角,Re是地球半径,取值6378.137km,ωe是地球自转角速度,L=ρV2SCL/2m为升力加速度,D=ρV2SCD/2m为阻力加速度。
2.过程约束
飞行器滑翔飞行中所要满足的热流、动压、过载与平衡滑翔条件等约束共同限制了滑翔弹道的可行范围,即滑翔飞行再入走廊。高超声速飞行器滑翔段的过程约束如下:
其中,KQ是热流密度相关常数,是最大热流密度,nmax是飞行器允许得最大过载,qmax为最大来流动压。
3.终端约束
终端约束主要为终端速度约束,终端高度约束,终端能量约束,弹目距离约束,,航向误差约束,终端倾侧角约束以及约束弹道倾角为一0°附近的小量。
步骤二:滑翔距离解析解;
1.积分形式的射程预测公式
(1)基于速度的射程预测公式
由射程导数公式
结合式(4),可得到射程对速度的导数
由于弹道倾角接近于0,假设cosγ≈1,sinγ≈0,则上述公式可以简化为
(2)基于能量的射程预测公式
由于飞行器在再入过程中是无动力的,能量逐渐衰减,且终端条件也可以由终端速度或终端能量定义,这里使用能量作为独立变量,其好处是不用考虑再入过程开始、结束时间的差异,且可以统一合理地描述不同轨迹段。定义能量为
E=gh+1/2V2 (13)
将上式对高度求导可得
结合式(1)和式(4)可得
将式(15)代入式(14)中,g项被抵消,可得
结合式(1)和式(10)可得
将式(17)代入式(16)中可得
得
考虑弹道倾角很小,即cosγ≈1,sinγ≈0,且量级很小,射程可由下式估算:
2.惯性力线化
由于地球自转引起的牵连加速度和科氏加速度均是时变的,且具有较强的非线性,无法直接积分。故对其进行线化处理。由于平稳滑翔段在再入飞行过程中占很大比例,主要对平稳滑翔段的惯性力进行线化处理。
在平稳滑翔段,弹道倾角很小,令式(5)为0,cosγ≈1,得到平稳滑翔段阻力满足:
上式结合式(12)和式(20),考虑可近似为常数,故仅需要线化cos2φ,sinφcosφcosψ和cosφsinψ项。cos2φ与sinφcosφcosψ前面的系数(即)很小,约为10-4级别,所以在此粗略的用能量的线性函数来拟合,拟合公式如下:
sinφcosφcosψ=k11E+k10 (22)
cos2φ=k21E+k20 (23)
不用拟合cosφsinψ是因为其值近似为常数,严谨的证明过程如下:
因为cosφsinψ与速度无关,所以可忽略速度影响,视飞行器仅受重力作用做圆周运动。由于重力对ze轴力矩为0,故可得
在NED坐标系下,由
转到地心赤道旋转坐标系中可得
所以
又因为飞行器仅受重力做匀速圆周运动,且Re+h≈Re为定值,cosγ≈1,可得
cosφsinψ≈const (29)
故cosφsinψ可视为常值,代入初始值φ0和ψ0求解此常值。
3.平稳滑翔段滑翔距离解析解推导
下面试推导平稳滑翔段滑翔距离解析解。在推导基于能量的滑翔距离解析解时,式(21)中的速度一次项难以处理,故在此将改速度项线化为能量的一次函数如下:
在此基础上,经过一系列代数变换,得到常值纵向升阻比的平稳滑翔段滑翔距离解析解如下表1。
表1
在旋转地球环境下通过仿真对比以上四种滑翔距离解析解的精度。
以下表1为解析解与弹道仿真结果比较
表2
由图2和以下表1可以看出,忽略地球自转影响会大幅降低滑翔距离解析解的精度,因此,推导滑翔距离解析解时考虑地球自转是十分必要的。尽管粗略地线化地球自转引起的惯性力造成一定的误差,但仿真结果表明,线化后得到的基于能量的解析解精度仍然很高。此外,仿真结果表明,无论是否考虑地球自转,基于能量的滑翔距离解析解精度都高于基于速度的滑翔距离解析解。对比式(11)与式(12),发现二者的区别在于基于速度的滑翔距离解析解简化时忽略了项,而基于能量的滑翔距离解析解忽略了项。可见前者忽略的势能的影响,并且易知前者的数量级远大于后者。在弹道倾角量级较大,即sinγ不严格等于0的情况下,以能量为自变量来预测射程相较以速度为自变量更为准确。同时,对比式(11)与式(12),易看出基于能量的滑翔距离解析解形式更为简单,故选取基于能量的滑翔距离解析解来进行射程预测。
步骤三:再入轨迹规划与跟踪制导;
1.基准攻角剖面规划
对于再入机动飞行器而言,射程和热防护是攻角剖面设计的主要考虑因素。针对CAV-H模型,这里设计攻角曲线为能量分段函数:
在再入初始下降阶段,为尽量拉高弹道,减小热流密度,采用较大攻角α1飞行在平稳滑翔段,为了发挥最大机动能力,飞行器以最大升阻比攻角滑翔,即α2=(L/D)max。出于控制系统跟踪性和气动稳定性的考虑,设计Es1和Es2之间的线性段来保证攻角剖面平稳过渡。
2.再入走廊设计
再入走廊上边界由式(5)在和σ=0的情况下求解得到。再入走廊下边界通过求解过程约束(7)-(9)得到,如图3所示。
3.阻力-能量剖面规划
由于初始段阻力加速度走廊狭窄,且走廊上边界可以近似为二次曲线,故将初始段的参考D(E)曲线设计成两个关于能量E的二次多项式,形式如下:
其中,A11,A12,A13,A21,A22,A23为常值。
再入走廊中部范围较大,考虑平稳滑翔条件,设计滑翔段参考D(E)曲线为一条直线。根据式(5),得:
即为关于能量E的一次曲线,比较剩余射程和预测射程大小,确定其与二次多项式阶段交点。其中,L1=Lcosσ,为升力在铅垂平面的分量,L1/D为纵向升阻比,可由平稳滑翔段滑翔距离解析解得到。
由于再入走廊末端较为宽阔,在线性段末端规划一次曲线段,以满足过程约束和终端约束。
D=k3E+b3(E<E3) (0.34)
其中,k3,b3为常值。规划的阻力加速度剖面如图4所示。
4.参考弹道跟踪指令
通过对公式(1)和阻力公式求导,进行一系列代数化简,忽略阻力系数CD的一阶、二阶导数,得到如下一个公式,其可以将参考阻力转换为纵向升阻比L1/D。
进而通过调节基准倾侧角σbsl来跟踪纵向升阻比曲线,实现倾侧反转。
5.倾侧角反转
侧向制导方法利用预先设计的航向误差门限控制倾侧角反转,将飞行器的侧向运动限定再预先设定的走廊内,当航向误差到达走廊边界时,令倾侧角反向,使横向运动向航向角误差减小的方向进行。
定义航向误差为
Δψ=ψ-ψLOS (0.37)
其中,ψLOS为视线方位角,即为在飞行器当前铅垂面内,飞行器指向目标的向量与当地北向的夹角。
为了保证Δψ能够进入航向误差走廊,再入初始阶段飞行器应朝终端目标点飞行。进入再入阶段后,倾侧角符号由下式确定:
其中,sign[φ(tk-1)]表示上一时刻由侧向制导算法给出的倾侧角符号,ΔψH(E)为航向误差门限。通常根据横程大小、倾侧角反转次数和终端误差等多方面约束将ΔψH(E)设计为能量的分段线性函数。按照终端航向角误差要求,将航向误差门限设计为如图5所示的漏斗形。
步骤四:基于地球自转模型的平稳滑翔弹道倾角推导;
采用基准攻角和基准倾侧角跟踪参考阻力-能量剖面,对比传统的PID反馈控制和Yu提出的平稳滑翔弹道倾角作为反馈项进行控制的弹道阻尼控制方法。所得到的弹道如图6所示,可以看出,弹道阻尼控制技术和PID控制都可以成功抑制弹道长周期弱阻尼振荡。但图6(b)显示,在平稳滑翔段,PID控制的阻力加速度跟踪曲线存在较大的跟踪误差,而弹道阻尼控制技术则可以较为准确的跟踪参考剖面,优势明显,故选用弹道阻尼技术抑制弹道振荡。
在旋转地球背景下,根据式(5)推导考虑地球自转的平稳滑翔弹道倾角γ'SG,满足
d2γ'SG 2+d1γ'SG+d0=0 (39)
各项系数为:
代入弹道阻尼技术中作为反馈项添加到基准攻角和倾侧角上进行控制:
αcmd=αbsl+cosσbslKγ(γSG-γ) (43)
用上述弹道阻尼控制技术可精确跟踪平稳滑翔弹道,有效抑制弹道振荡。
实施案例:
为了验证算法的精度,以CAV-H作为再入飞行器模型在七种不同目标状态下进行仿真。再入初始状态均为h0=80km,V0=7000m/s,λ0=0°,θ0=45°,γ0=0°。ψ0依据目标方向取值。弹道阻尼反馈项系数Kγ=5。终端约束为Vf=2000m/s,hf=25km,弹目距离Rtm≤5km。终端航向角误差|Δψf|≤5°。仿真结果如图7-图10。
表3
图7中,按照逆时针方向从左至右,分别是七种不同的目标位置,即从西北方位到东北方位。图8-图10和表3(终端状态量)显示了弹道阻尼控制技术能控制飞行器在不同目标状态下成功抑制弹道振荡,精准命中目标。以东南方向的目标位置为例具体分析参考剖面跟踪情况。
由图11(a)可明显看出,所得弹道处于狭窄的高度走廊内,满足过程约束。同时,相较于原先再入弹道,弹道振荡的情况得到了极大的缓解,弹道的平稳性大大增强。但在初始下降阶段,弹道仍存在一个波峰的振荡,可通过增大弹道阻尼系数来消除。图11(b)中,实际速度基本跟踪上了参考速度,达到了末端速度要求。图11(c)中弹道倾角精确跟踪参考弹道倾角,基本稳定在0附近,满足要求。图11(d)为倾侧角时间曲线,结合图11(e)和图11(f),可看出在仿真过程中,倾侧角翻转了6次,末端横向误差|Δψf|≤5°,基本命中目标点。
综上所述,通过上述步骤,推导出了本发明方法,即一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法,案例仿真结果表明本发明方法能够精准跟踪平稳滑翔弹道,成功抑制弹道振荡,提高导弹的命中精度,综合性能优异。
Claims (1)
1.一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法,其特征在于:它包括以下几个步骤:
步骤一:高超声速飞行器建模
在旋转地球背景下,采用标准大气模型和平方反比引力场模型,建立高超声速飞行器的六自由度动力学方程如下:
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>h</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>V</mi>
<mi> </mi>
<mi>sin</mi>
<mi>&gamma;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>R</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>V</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&gamma;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>&lambda;</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>V</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&gamma;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&psi;</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&phi;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>&phi;</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>V</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&gamma;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&psi;</mi>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>V</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mi>D</mi>
<mo>-</mo>
<mi>g</mi>
<mi> </mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&gamma;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mi>&phi;</mi>
<mi>cos</mi>
<mi>&phi;</mi>
<mi>cos</mi>
<mi>&gamma;</mi>
<mi>cos</mi>
<mi>&psi;</mi>
<mo>+</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
<mi>sin</mi>
<mi>&gamma;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>&gamma;</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mi>L</mi>
<mi>V</mi>
</mfrac>
<mi>cos</mi>
<mi>&sigma;</mi>
<mo>-</mo>
<mfrac>
<mrow>
<mi>g</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&gamma;</mi>
</mrow>
<mi>V</mi>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>V</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&gamma;</mi>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mi>V</mi>
</mfrac>
<mo>&lsqb;</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
<mi>cos</mi>
<mi>&gamma;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mi>&phi;</mi>
<mi>cos</mi>
<mi>&phi;</mi>
<mi>sin</mi>
<mi>&gamma;</mi>
<mi>cos</mi>
<mi>&psi;</mi>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>V&omega;</mi>
<mi>e</mi>
</msub>
<mi>cos</mi>
<mi>&phi;</mi>
<mi>sin</mi>
<mi>&psi;</mi>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>&psi;</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>L</mi>
<mi> </mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&sigma;</mi>
</mrow>
<mrow>
<mi>V</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&gamma;</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>V</mi>
<mi> </mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&phi;</mi>
<mi>cos</mi>
<mi>&gamma;</mi>
<mi>sin</mi>
<mi>&psi;</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
<mi>cos</mi>
<mi>&phi;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&phi;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&phi;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&psi;</mi>
</mrow>
<mrow>
<mi>V</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&gamma;</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<mi>cos</mi>
<mi>&gamma;</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mi>cos</mi>
<mi>&gamma;</mi>
<mo>-</mo>
<mi>cos</mi>
<mi>&phi;</mi>
<mi>sin</mi>
<mi>&gamma;</mi>
<mi>cos</mi>
<mi>&psi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,d/dt表示对时间求导,t是时间,h是海拔高度,R是射程,λ是经度,φ是纬度,V是飞行器相对于地球的速率,γ是弹道倾角,ψ是飞行器航向角,以当地北向为基准,σ是倾侧角,Re是地球半径,取值6378.137km,ωe是地球自转角速度,L=ρV2SCL/2m为升力加速度,D=ρV2SCD/2m为阻力加速度;
步骤二:滑翔距离解析解推导
1.积分形式的射程预测公式
在旋转地球背景下,考虑由地球自转引起的惯性力,基于步骤一中的动力学模型推导积分形式的射程预测公式;
(1)基于速度的射程预测公式
由射程导数公式(2)与式(5),可得到射程对速度的导数;假设cosγ≈1,sinγ≈0,得到考虑地球自转的基于速度的射程预测公式
<mrow>
<mi>R</mi>
<mo>=</mo>
<mo>-</mo>
<mo>&Integral;</mo>
<mfrac>
<mrow>
<mi>V</mi>
<mi>d</mi>
<mi>V</mi>
</mrow>
<mrow>
<mi>D</mi>
<mo>+</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mi>&phi;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&phi;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&psi;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
(2)基于能量的射程预测公式
由于飞行器在再入过程中是无动力的,能量逐渐衰减,且终端条件也可以由终端速度或终端能量定义,这里使用能量作为独立变量,其好处是不用考虑再入过程开始、结束时间的差异,且可以统一合理地描述不同轨迹段;定义能量为
E=gh+1/2V2 (9)
将上式对高度求导可得
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>E</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>h</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>g</mi>
<mo>+</mo>
<mi>V</mi>
<mfrac>
<mrow>
<mi>d</mi>
<mi>V</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>h</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
根据式(1),(2),(5)及式(10),经过一系列代数化简,假设cosγ≈1,sinγ≈0,以及量级很小,可得基于能量的考虑地球自转的射程积分公式
<mrow>
<mi>R</mi>
<mo>=</mo>
<mo>-</mo>
<mo>&Integral;</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>E</mi>
</mrow>
<mrow>
<mi>D</mi>
<mo>+</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&phi;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&phi;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&psi;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
2.惯性力线化及证明
由于地球自转引起的牵连加速度和科氏加速度均是时变的,且具有较强的非线性,无法直接积分,故对其进行线化处理;由于平稳滑翔段在再入飞行过程中占很大比例,主要对平稳滑翔段的惯性力进行线化处理;
在平稳滑翔段,弹道倾角很小,令式(6)为0,cosγ≈1,得到平稳滑翔段阻力满足:
<mrow>
<mi>D</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>g</mi>
<mo>-</mo>
<mfrac>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
<mrow>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
</mrow>
</mfrac>
<mo>-</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>V&omega;</mi>
<mi>e</mi>
</msub>
<mi>cos</mi>
<mi>&phi;</mi>
<mi>sin</mi>
<mi>&psi;</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
上式结合式(8),式(11),考虑可近似为常数,故仅需要线化cos2φ,sinφcosφcosψ和cosφsinψ项;cos2φ与sinφcosφcosψ前面的系数(即)很小,约为10-4级别,所以在此粗略的用能量的线性函数来拟合,拟合公式如下:
sinφcosφcosψ=k11E+k10 (13)
cos2φ=k21E+k20 (14)
不用拟合cosφsinψ是因为其值近似为常数,严谨的证明过程如下:
因为cosφsinψ与速度无关,所以可忽略速度影响,视飞行器仅受重力作用做圆周运动;由于重力对ze轴力矩为0,故可得
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<msub>
<mover>
<mi>H</mi>
<mo>&RightArrow;</mo>
</mover>
<msub>
<mi>z</mi>
<mi>e</mi>
</msub>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
在NED坐标系下,由
<mrow>
<msup>
<mover>
<mi>V</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>E</mi>
<mi>D</mi>
</mrow>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>V</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&gamma;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&psi;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>V</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&gamma;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&psi;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>V</mi>
<mi> </mi>
<mi>sin</mi>
<mi>&gamma;</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mover>
<mi>R</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>E</mi>
<mi>D</mi>
</mrow>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
转到地心赤道旋转坐标系中可得
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mover>
<mi>H</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>G</mi>
<mi>E</mi>
<mi>R</mi>
</mrow>
</msup>
<mo>=</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>N</mi>
<mi>E</mi>
<mi>D</mi>
</mrow>
<mrow>
<mi>G</mi>
<mi>E</mi>
<mi>R</mi>
</mrow>
</msubsup>
<msup>
<mover>
<mi>H</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>E</mi>
<mi>D</mi>
</mrow>
</msup>
<mo>=</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>N</mi>
<mi>E</mi>
<mi>D</mi>
</mrow>
<mrow>
<mi>G</mi>
<mi>E</mi>
<mi>R</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msup>
<mover>
<mi>R</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>E</mi>
<mi>D</mi>
</mrow>
</msup>
<mo>&times;</mo>
<mi>m</mi>
<msup>
<mover>
<mi>V</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>N</mi>
<mi>E</mi>
<mi>D</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>m</mi>
<mi>V</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&gamma;</mi>
<mi>sin</mi>
<mi>&psi;</mi>
<mi>cos</mi>
<mi>&lambda;</mi>
<mi>sin</mi>
<mi>&phi;</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>m</mi>
<mi>V</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&gamma;</mi>
<mi>cos</mi>
<mi>&psi;</mi>
<mi>sin</mi>
<mi>&lambda;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>m</mi>
<mi>V</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&gamma;</mi>
<mi>sin</mi>
<mi>&psi;</mi>
<mi>sin</mi>
<mi>&lambda;</mi>
<mi>sin</mi>
<mi>&phi;</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>m</mi>
<mi>V</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&gamma;</mi>
<mi>cos</mi>
<mi>&psi;</mi>
<mi>cos</mi>
<mi>&lambda;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
<mi>m</mi>
<mi>V</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&gamma;</mi>
<mi>cos</mi>
<mi>&phi;</mi>
<mi>sin</mi>
<mi>&psi;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
所以
<mrow>
<msub>
<mover>
<mi>H</mi>
<mo>&RightArrow;</mo>
</mover>
<msub>
<mi>z</mi>
<mi>e</mi>
</msub>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>m</mi>
<mi>V</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&gamma;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&phi;</mi>
<mi>sin</mi>
<mi>&psi;</mi>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
<mi>s</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
又因为飞行器仅受重力做匀速圆周运动,且Re+h≈Re为定值,cosγ≈1,可得
cosφsinψ≈const (20)
故cosφsinψ可视为常值,代入初始值φ0和ψ0求解此常值;
3.平稳滑翔段滑翔距离解析解推导
在推导基于能量的滑翔距离解析解时,式(12)中的速度一次项难以处理,故在此将改速度项线化为能量的一次函数如下:
<mrow>
<mi>V</mi>
<mo>=</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>V</mi>
<mn>1</mn>
</mrow>
</msub>
<mi>E</mi>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>V</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mi>E</mi>
<mo>+</mo>
<msub>
<mi>V</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<msub>
<mi>E</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>21</mn>
<mo>)</mo>
</mrow>
</mrow>
在此基础上,经过一系列代数变换,得到常值纵向升阻比的平稳滑翔段滑翔距离解析解如下:
(1)基于速度的滑翔距离解析解
另
若Δ≥0
<mrow>
<mi>R</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<mi>D</mi>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mi>l</mi>
<mi>n</mi>
<mfrac>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>-</mo>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>-</mo>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
<mi>l</mi>
<mi>n</mi>
<mfrac>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>-</mo>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>22</mn>
<mo>)</mo>
</mrow>
</mrow>
若Δ<0
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>R</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<mi>D</mi>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mi>ln</mi>
<mfrac>
<mrow>
<mo>(</mo>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mi>V</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msubsup>
<mi>V</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>V</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msqrt>
<mrow>
<mo>-</mo>
<mi>&Delta;</mi>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mi>arctan</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>V</mi>
<mo>+</mo>
<mfrac>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mfrac>
</mrow>
<msqrt>
<mrow>
<mo>-</mo>
<mi>&Delta;</mi>
</mrow>
</msqrt>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>arctan</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mfrac>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mfrac>
</mrow>
<msqrt>
<mrow>
<mo>-</mo>
<mi>&Delta;</mi>
</mrow>
</msqrt>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>23</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,
a2=(k11+k21)ωe 2(Re+h)2/2+1 (24)
<mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mi>&phi;</mi>
<mi>sin</mi>
<mi>&psi;</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mn>11</mn>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>21</mn>
</msub>
<mo>)</mo>
<msup>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>/</mo>
<mn>2</mn>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>25</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>-</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mn>21</mn>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>11</mn>
</msub>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
<mi>g</mi>
<mi>h</mi>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>20</mn>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>k</mi>
<mn>10</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mn>11</mn>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>21</mn>
</msub>
<mo>)</mo>
<msup>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>/</mo>
<mn>2</mn>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>26</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msqrt>
<mi>&Delta;</mi>
</msqrt>
</mrow>
<mn>2</mn>
</mfrac>
<mo>,</mo>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msqrt>
<mi>&Delta;</mi>
</msqrt>
</mrow>
<mn>2</mn>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>27</mn>
<mo>)</mo>
</mrow>
</mrow>
(2)基于能量的滑翔距离解析解
<mrow>
<mi>R</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mi>E</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<msub>
<mi>E</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>28</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,
b1=2+2kV1ωe(Re+h)cosφsinψ+k21ωe 2(Re+h)2+k11ωe 2(Re+h) (29)
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mn>2</mn>
<msub>
<mi>k</mi>
<mrow>
<mi>V</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mi>&phi;</mi>
<mi>sin</mi>
<mi>&psi;</mi>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>20</mn>
</msub>
<msup>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>10</mn>
</msub>
<msup>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mn>3</mn>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>30</mn>
<mo>)</mo>
</mrow>
</mrow>
对比分析可知,基于能量的滑翔距离解析解精度更高,形式更为简洁,故以能量为自变量进行再入轨迹规划;
步骤三:再入轨迹规划与跟踪制导
1.基准攻角剖面规划
对于再入机动飞行器而言,射程和热防护是攻角剖面设计的主要考虑因素;针对CAV-H模型,这里设计攻角曲线为能量分段函数:
<mrow>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>b</mi>
<mi>s</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>E</mi>
<mo>></mo>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mfrac>
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mi>&alpha;</mi>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo><</mo>
<mi>E</mi>
<mo>&le;</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>E</mi>
<mo>&le;</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>31</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Es1和Es2表示攻角分段函数的分界点;在再入初始下降阶段,为尽量拉高弹道,减小热流密度,采用较大攻角α1飞行在平稳滑翔段,为了发挥最大机动能力,飞行器以最大升阻比攻角滑翔,即α2=(L/D)max;出于控制系统跟踪性和气动稳定性的考虑,设计Es1和Es2之间的线性段来保证攻角剖面平稳过渡;
2.阻力-能量剖面规划
由于初始段阻力加速度走廊狭窄,且走廊上边界可以近似为二次曲线,故将初始段的参考D(E)曲线设计成两个关于能量E的二次多项式,形式如下:
<mrow>
<mi>D</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mn>11</mn>
</msub>
<msup>
<mi>E</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>12</mn>
</msub>
<mi>E</mi>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>13</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>E</mi>
<mo>></mo>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mn>21</mn>
</msub>
<msup>
<mi>E</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>22</mn>
</msub>
<mi>E</mi>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<mo><</mo>
<mi>E</mi>
<mo>&le;</mo>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>32</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,A11,A12,A13,A21,A22,A23为常值;
再入走廊中部范围较大,考虑平稳滑翔条件,设计滑翔段参考D(E)曲线为一条直线;根据式(6),得:
<mrow>
<mi>D</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>g</mi>
<mo>-</mo>
<mfrac>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
<mrow>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
</mrow>
</mfrac>
<mo>-</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>V&omega;</mi>
<mi>e</mi>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&phi;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&psi;</mi>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<mi>D</mi>
<mo>)</mo>
</mrow>
</mfrac>
<mo>,</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mn>3</mn>
</msub>
<mo>&le;</mo>
<mi>E</mi>
<mo>&le;</mo>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>33</mn>
<mo>)</mo>
</mrow>
</mrow>
即为关于能量E的一次曲线,比较剩余射程和预测射程大小,确定其与二次多项式阶段交点;其中,L1=Lcosσ,为升力在铅垂平面的分量,L1/D为纵向升阻比,可由步骤二中平稳滑翔段滑翔距离解析解得到;
由于再入走廊末端较为宽阔,在线性段末端规划一次曲线段,以满足过程约束和终端约束;
D=k3E+b3 (E<E3) (34)
其中,k3,b3为常值;
3.参考弹道跟踪指令
通过对公式(1)和阻力公式求导,进行一系列代数化简,忽略阻力系数CD的一阶、二阶导数,得到如下一个公式,其可以将参考阻力转换为纵向升阻比L1/D;
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mi>D</mi>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>h</mi>
<mi>S</mi>
</msub>
<msup>
<mi>D</mi>
<mn>2</mn>
</msup>
</mfrac>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mover>
<mi>D</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>-</mo>
<mover>
<mi>D</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mover>
<mi>D</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>D</mi>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mi>D</mi>
</mrow>
<mi>V</mi>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mn>4</mn>
<msup>
<mi>D</mi>
<mn>3</mn>
</msup>
</mrow>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>+</mo>
<mfrac>
<mi>D</mi>
<msub>
<mi>h</mi>
<mi>S</mi>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
<mrow>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>g</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>D</mi>
</mfrac>
<mo>&lsqb;</mo>
<msup>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>V&omega;</mi>
<mi>e</mi>
</msub>
<mi>cos</mi>
<mi>&phi;</mi>
<mi>sin</mi>
<mi>&psi;</mi>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<msub>
<mi>h</mi>
<mi>S</mi>
</msub>
<mrow>
<mi>D</mi>
<mi>V</mi>
</mrow>
</mfrac>
<mo>(</mo>
<mrow>
<mfrac>
<mover>
<mi>D</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>D</mi>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>D</mi>
</mrow>
<mi>V</mi>
</mfrac>
</mrow>
<mo>)</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mi>&phi;</mi>
<mi>cos</mi>
<mi>&phi;</mi>
<mi>cos</mi>
<mi>&psi;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>35</mn>
<mo>)</mo>
</mrow>
</mrow>
进而通过调节基准倾侧角σbsl来跟踪纵向升阻比曲线,实现倾侧反转;
<mrow>
<mo>|</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>b</mi>
<mi>s</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>=</mo>
<msup>
<mi>cos</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<mi>D</mi>
</mrow>
<mrow>
<mi>L</mi>
<mo>/</mo>
<mi>D</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>36</mn>
<mo>)</mo>
</mrow>
</mrow>
飞行器横向运动通过航向误差门限进行控制;
步骤四:基于地球自转模型的平稳滑翔弹道倾角推导
弹道阻尼技术以平稳滑翔弹道倾角作为反馈项进行控制,可以成功抑制再入弹道存在的长周期弱阻尼振荡;在平稳滑翔段,相比于传统的PID反馈控制,弹道阻尼控制技术跟踪精确性更高,形式简单,便于应用;在旋转地球背景下,根据式(6)推导考虑地球自转的平稳滑翔弹道倾角γ'SG,满足
<mrow>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
<msup>
<msubsup>
<mi>&gamma;</mi>
<mrow>
<mi>S</mi>
<mi>G</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>&gamma;</mi>
<mrow>
<mi>S</mi>
<mi>G</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>37</mn>
<mo>)</mo>
</mrow>
</mrow>
各项系数为:
<mrow>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&phi;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&phi;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&psi;</mi>
<mi>V</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>38</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>L</mi>
<mi> </mi>
<msub>
<mi>cos&sigma;</mi>
<mrow>
<mi>b</mi>
<mi>s</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>S</mi>
</msub>
</mfrac>
<mi>V</mi>
<mo>-</mo>
<mfrac>
<msup>
<mi>V</mi>
<mn>3</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mo>-</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>L</mi>
</mrow>
<mi>V</mi>
</mfrac>
<msub>
<mi>cos&sigma;</mi>
<mrow>
<mi>b</mi>
<mi>s</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>V</mi>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&phi;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&psi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&phi;</mi>
<mi>V</mi>
<mo>+</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mn>2</mn>
<mi>&phi;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&psi;</mi>
<mover>
<mi>&phi;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>-</mo>
<mi>sin</mi>
<mi>&phi;</mi>
<mi>cos</mi>
<mi>&phi;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&psi;</mi>
<mover>
<mi>&psi;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>39</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>L</mi>
</mrow>
<mi>V</mi>
</mfrac>
<msub>
<mi>cos&sigma;</mi>
<mrow>
<mi>b</mi>
<mi>s</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>V</mi>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&phi;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&psi;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mi>&rho;V</mi>
<mn>2</mn>
</msup>
<mi>s</mi>
</mrow>
<mrow>
<mn>2</mn>
<mi>m</mi>
</mrow>
</mfrac>
<msub>
<mover>
<mi>C</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>L</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mover>
<mi>&phi;</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mn>2</mn>
<mi>&phi;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>40</mn>
<mo>)</mo>
</mrow>
</mrow>
代入弹道阻尼技术中作为反馈项添加到基准攻角和倾侧角上进行控制:
αcmd=αbsl+cosσbslKγ(γ'SG-γ) (41)
<mrow>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>c</mi>
<mi>m</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>b</mi>
<mi>s</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>b</mi>
<mi>s</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mi>&gamma;</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&gamma;</mi>
<mrow>
<mi>S</mi>
<mi>G</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<mi>&gamma;</mi>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>b</mi>
<mi>s</mi>
<mi>l</mi>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>42</mn>
<mo>)</mo>
</mrow>
</mrow>
用上述弹道阻尼控制技术可精确跟踪平稳滑翔弹道,有效抑制弹道振荡;
通过上述四个步骤,可使弹道稳定在平稳滑翔弹道附近,最终得到了一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710970667.8A CN107941087B (zh) | 2017-10-18 | 2017-10-18 | 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710970667.8A CN107941087B (zh) | 2017-10-18 | 2017-10-18 | 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107941087A true CN107941087A (zh) | 2018-04-20 |
CN107941087B CN107941087B (zh) | 2019-08-06 |
Family
ID=61936132
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710970667.8A Active CN107941087B (zh) | 2017-10-18 | 2017-10-18 | 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107941087B (zh) |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108549785A (zh) * | 2018-05-03 | 2018-09-18 | 中国人民解放军国防科技大学 | 一种基于三维飞行剖面的高超声速飞行器精准弹道快速预测方法 |
CN109190248A (zh) * | 2018-09-03 | 2019-01-11 | 中国运载火箭技术研究院 | 一种用于滑翔飞行器的滑翔射程解析方法及解析系统 |
CN109240335A (zh) * | 2018-10-29 | 2019-01-18 | 北京控制工程研究所 | 一种空天飞行器进场着陆制导方法 |
CN109240323A (zh) * | 2018-11-02 | 2019-01-18 | 北京控制工程研究所 | 一种实时解析构造的空天飞行器再入制导方法 |
CN109446582A (zh) * | 2018-09-29 | 2019-03-08 | 北京航空航天大学 | 一种考虑地球自转的高精度降阶平稳滑翔动力学建模方法 |
CN109740198A (zh) * | 2018-12-14 | 2019-05-10 | 中国人民解放军国防科技大学 | 一种基于解析预测的滑翔飞行器三维再入制导方法 |
CN109977543A (zh) * | 2019-03-26 | 2019-07-05 | 中国人民解放军国防科技大学 | 基于侧向优先的三维剖面覆盖区域计算方法、系统及介质 |
CN110309590A (zh) * | 2019-06-28 | 2019-10-08 | 北京理工大学 | 一种再入飞行器速度-高度再入走廊预测方法 |
CN110471450A (zh) * | 2019-08-29 | 2019-11-19 | 大连理工大学 | 在高度速度剖面内直接规划再入轨迹的方法 |
CN111366044A (zh) * | 2019-12-29 | 2020-07-03 | 湖北航天飞行器研究所 | 一种平飞过渡段制导控制方法 |
CN111488646A (zh) * | 2020-03-02 | 2020-08-04 | 北京航空航天大学 | 一种旋转地球下高超声速平稳滑翔弹道的解析求解方法 |
CN111580535A (zh) * | 2020-05-18 | 2020-08-25 | 中国人民解放军国防科技大学 | 一种基于凸优化的再入轨迹三维剖面规划方法及系统 |
CN111580555A (zh) * | 2020-05-13 | 2020-08-25 | 北京控制工程研究所 | 一种高超声速飞行器上升段分段自适应预测校正制导方法 |
CN112034879A (zh) * | 2019-12-16 | 2020-12-04 | 哈尔滨工业大学 | 一种基于高度-射程比的标准轨迹跟踪制导方法 |
CN112278336A (zh) * | 2020-11-03 | 2021-01-29 | 北京空间飞行器总体设计部 | 高速再入热流飞行试验验证方法 |
CN112507467A (zh) * | 2020-12-21 | 2021-03-16 | 中国人民解放军96901部队23分队 | 基于阻力和纵向升阻比的滑翔弹道随速度变化降阶解 |
CN112507465A (zh) * | 2020-12-21 | 2021-03-16 | 中国人民解放军96901部队23分队 | 基于阻力和纵向升阻比的滑翔弹道随能量变化降阶解 |
CN112861251A (zh) * | 2020-12-21 | 2021-05-28 | 中国人民解放军96901部队23分队 | 考虑地球自转速率的滑翔弹道随速度变化降阶解 |
CN112861249A (zh) * | 2020-12-21 | 2021-05-28 | 中国人民解放军96901部队23分队 | 一种基于攻角和阻力加速度的滑翔弹道随速度变化降阶解 |
CN112861250A (zh) * | 2020-12-21 | 2021-05-28 | 中国人民解放军96901部队23分队 | 一种基于攻角和阻力加速度的滑翔弹道随能量变化降阶解 |
CN113093790A (zh) * | 2021-03-22 | 2021-07-09 | 北京航空航天大学 | 一种基于解析模型的飞行器再入滑翔轨迹规划方法 |
CN114167888A (zh) * | 2021-11-19 | 2022-03-11 | 湖北航天技术研究院总体设计所 | 一种滑翔高超声速飞行器末端位置和速度控制方法 |
CN114355777A (zh) * | 2022-01-06 | 2022-04-15 | 浙江大学 | 一种基于分布式压力传感器与分段姿态控制的动态滑翔方法及系统 |
CN114895700A (zh) * | 2022-04-06 | 2022-08-12 | 北京理工大学 | 带终端时间约束的进入飞行器阻力剖面规划方法 |
CN116991171A (zh) * | 2023-07-12 | 2023-11-03 | 中国科学院力学研究所 | 基于能量的无动力滑降无人机着陆下滑道动态生成方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100250031A1 (en) * | 2009-03-27 | 2010-09-30 | The Charles Stark Draper Laboratory, Inc. | Propulsive guidance for atmospheric skip entry trajectories |
RU2408847C1 (ru) * | 2009-12-24 | 2011-01-10 | Открытое акционерное общество "Концерн радиостроения "Вега" | Способ самонаведения летательных аппаратов на гиперзвуковые цели |
CN104731104A (zh) * | 2015-03-09 | 2015-06-24 | 北京航天自动控制研究所 | 一种高超声速飞行器滑翔飞行段的纵向制导方法 |
CN105354380A (zh) * | 2015-11-03 | 2016-02-24 | 中国人民解放军国防科学技术大学 | 面向摄动因素影响补偿的滑翔弹道快速修正方法 |
CN105550402A (zh) * | 2015-12-07 | 2016-05-04 | 北京航空航天大学 | 一种基于攻角或倾侧角变频的高超平稳机动滑翔弹道设计方法 |
JP2016182874A (ja) * | 2015-03-26 | 2016-10-20 | 三菱電機株式会社 | リエントリ宇宙機誘導制御システム |
CN106054604A (zh) * | 2016-06-02 | 2016-10-26 | 西北工业大学 | 基于模型预测控制理论的再入飞行器鲁棒最优制导方法 |
CN106547991A (zh) * | 2016-11-25 | 2017-03-29 | 中国工程物理研究院总体工程研究所 | 沿滑翔弹道的扰动引力重构模型优化方法 |
-
2017
- 2017-10-18 CN CN201710970667.8A patent/CN107941087B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100250031A1 (en) * | 2009-03-27 | 2010-09-30 | The Charles Stark Draper Laboratory, Inc. | Propulsive guidance for atmospheric skip entry trajectories |
RU2408847C1 (ru) * | 2009-12-24 | 2011-01-10 | Открытое акционерное общество "Концерн радиостроения "Вега" | Способ самонаведения летательных аппаратов на гиперзвуковые цели |
CN104731104A (zh) * | 2015-03-09 | 2015-06-24 | 北京航天自动控制研究所 | 一种高超声速飞行器滑翔飞行段的纵向制导方法 |
JP2016182874A (ja) * | 2015-03-26 | 2016-10-20 | 三菱電機株式会社 | リエントリ宇宙機誘導制御システム |
CN105354380A (zh) * | 2015-11-03 | 2016-02-24 | 中国人民解放军国防科学技术大学 | 面向摄动因素影响补偿的滑翔弹道快速修正方法 |
CN105550402A (zh) * | 2015-12-07 | 2016-05-04 | 北京航空航天大学 | 一种基于攻角或倾侧角变频的高超平稳机动滑翔弹道设计方法 |
CN106054604A (zh) * | 2016-06-02 | 2016-10-26 | 西北工业大学 | 基于模型预测控制理论的再入飞行器鲁棒最优制导方法 |
CN106547991A (zh) * | 2016-11-25 | 2017-03-29 | 中国工程物理研究院总体工程研究所 | 沿滑翔弹道的扰动引力重构模型优化方法 |
Cited By (44)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108549785A (zh) * | 2018-05-03 | 2018-09-18 | 中国人民解放军国防科技大学 | 一种基于三维飞行剖面的高超声速飞行器精准弹道快速预测方法 |
CN108549785B (zh) * | 2018-05-03 | 2021-09-24 | 中国人民解放军国防科技大学 | 一种基于三维飞行剖面的高超声速飞行器精准弹道快速预测方法 |
CN109190248B (zh) * | 2018-09-03 | 2023-07-18 | 中国运载火箭技术研究院 | 一种用于滑翔飞行器的滑翔射程解析方法及解析系统 |
CN109190248A (zh) * | 2018-09-03 | 2019-01-11 | 中国运载火箭技术研究院 | 一种用于滑翔飞行器的滑翔射程解析方法及解析系统 |
CN109446582A (zh) * | 2018-09-29 | 2019-03-08 | 北京航空航天大学 | 一种考虑地球自转的高精度降阶平稳滑翔动力学建模方法 |
CN109446582B (zh) * | 2018-09-29 | 2023-04-07 | 北京航空航天大学 | 一种考虑地球自转的高精度降阶平稳滑翔动力学建模方法 |
CN109240335A (zh) * | 2018-10-29 | 2019-01-18 | 北京控制工程研究所 | 一种空天飞行器进场着陆制导方法 |
CN109240323A (zh) * | 2018-11-02 | 2019-01-18 | 北京控制工程研究所 | 一种实时解析构造的空天飞行器再入制导方法 |
CN109240323B (zh) * | 2018-11-02 | 2021-07-13 | 北京控制工程研究所 | 一种实时解析构造的空天飞行器再入制导方法 |
CN109740198B (zh) * | 2018-12-14 | 2022-12-23 | 中国人民解放军国防科技大学 | 一种基于解析预测的滑翔飞行器三维再入制导方法 |
CN109740198A (zh) * | 2018-12-14 | 2019-05-10 | 中国人民解放军国防科技大学 | 一种基于解析预测的滑翔飞行器三维再入制导方法 |
CN109977543B (zh) * | 2019-03-26 | 2023-01-06 | 中国人民解放军国防科技大学 | 基于侧向优先的三维剖面覆盖区域计算方法、系统及介质 |
CN109977543A (zh) * | 2019-03-26 | 2019-07-05 | 中国人民解放军国防科技大学 | 基于侧向优先的三维剖面覆盖区域计算方法、系统及介质 |
CN110309590A (zh) * | 2019-06-28 | 2019-10-08 | 北京理工大学 | 一种再入飞行器速度-高度再入走廊预测方法 |
WO2021036778A1 (zh) * | 2019-08-29 | 2021-03-04 | 大连理工大学 | 在高度速度剖面内直接规划再入轨迹的方法 |
CN110471450A (zh) * | 2019-08-29 | 2019-11-19 | 大连理工大学 | 在高度速度剖面内直接规划再入轨迹的方法 |
US11079239B2 (en) | 2019-08-29 | 2021-08-03 | Dalian University Of Technology | Method for directly planning reentry trajectory in height-velocity profile |
CN112034879A (zh) * | 2019-12-16 | 2020-12-04 | 哈尔滨工业大学 | 一种基于高度-射程比的标准轨迹跟踪制导方法 |
CN111366044A (zh) * | 2019-12-29 | 2020-07-03 | 湖北航天飞行器研究所 | 一种平飞过渡段制导控制方法 |
CN111488646A (zh) * | 2020-03-02 | 2020-08-04 | 北京航空航天大学 | 一种旋转地球下高超声速平稳滑翔弹道的解析求解方法 |
CN111488646B (zh) * | 2020-03-02 | 2022-04-01 | 北京航空航天大学 | 一种旋转地球下高超声速平稳滑翔弹道的解析求解方法 |
CN111580555A (zh) * | 2020-05-13 | 2020-08-25 | 北京控制工程研究所 | 一种高超声速飞行器上升段分段自适应预测校正制导方法 |
CN111580535A (zh) * | 2020-05-18 | 2020-08-25 | 中国人民解放军国防科技大学 | 一种基于凸优化的再入轨迹三维剖面规划方法及系统 |
CN111580535B (zh) * | 2020-05-18 | 2023-06-06 | 中国人民解放军国防科技大学 | 一种基于凸优化的再入轨迹三维剖面规划方法及系统 |
CN112278336B (zh) * | 2020-11-03 | 2022-03-22 | 北京空间飞行器总体设计部 | 高速再入热流飞行试验验证方法 |
CN112278336A (zh) * | 2020-11-03 | 2021-01-29 | 北京空间飞行器总体设计部 | 高速再入热流飞行试验验证方法 |
CN112861250B (zh) * | 2020-12-21 | 2023-03-28 | 中国人民解放军96901部队23分队 | 基于攻角和阻力的滑翔弹道随能量变化降阶解计算方法 |
CN112507467A (zh) * | 2020-12-21 | 2021-03-16 | 中国人民解放军96901部队23分队 | 基于阻力和纵向升阻比的滑翔弹道随速度变化降阶解 |
CN112861250A (zh) * | 2020-12-21 | 2021-05-28 | 中国人民解放军96901部队23分队 | 一种基于攻角和阻力加速度的滑翔弹道随能量变化降阶解 |
CN112507467B (zh) * | 2020-12-21 | 2023-05-12 | 中国人民解放军96901部队23分队 | 基于阻力和升阻比的滑翔弹道随速度变化降阶解计算方法 |
CN112861251A (zh) * | 2020-12-21 | 2021-05-28 | 中国人民解放军96901部队23分队 | 考虑地球自转速率的滑翔弹道随速度变化降阶解 |
CN112507465A (zh) * | 2020-12-21 | 2021-03-16 | 中国人民解放军96901部队23分队 | 基于阻力和纵向升阻比的滑翔弹道随能量变化降阶解 |
CN112861251B (zh) * | 2020-12-21 | 2023-03-28 | 中国人民解放军96901部队23分队 | 考虑地球自转速率的滑翔弹道随速度变化降阶解计算方法 |
CN112507465B (zh) * | 2020-12-21 | 2023-04-14 | 中国人民解放军96901部队23分队 | 基于阻力和升阻比的滑翔弹道随能量变化降阶解计算方法 |
CN112861249B (zh) * | 2020-12-21 | 2023-03-31 | 中国人民解放军96901部队23分队 | 基于攻角和阻力的滑翔弹道随速度变化降阶解计算方法 |
CN112861249A (zh) * | 2020-12-21 | 2021-05-28 | 中国人民解放军96901部队23分队 | 一种基于攻角和阻力加速度的滑翔弹道随速度变化降阶解 |
CN113093790A (zh) * | 2021-03-22 | 2021-07-09 | 北京航空航天大学 | 一种基于解析模型的飞行器再入滑翔轨迹规划方法 |
CN114167888A (zh) * | 2021-11-19 | 2022-03-11 | 湖北航天技术研究院总体设计所 | 一种滑翔高超声速飞行器末端位置和速度控制方法 |
CN114167888B (zh) * | 2021-11-19 | 2023-06-20 | 湖北航天技术研究院总体设计所 | 一种滑翔高超声速飞行器末端位置和速度控制方法 |
CN114355777A (zh) * | 2022-01-06 | 2022-04-15 | 浙江大学 | 一种基于分布式压力传感器与分段姿态控制的动态滑翔方法及系统 |
CN114355777B (zh) * | 2022-01-06 | 2023-10-10 | 浙江大学 | 一种基于分布式压力传感器与分段姿态控制的动态滑翔方法及系统 |
CN114895700A (zh) * | 2022-04-06 | 2022-08-12 | 北京理工大学 | 带终端时间约束的进入飞行器阻力剖面规划方法 |
CN116991171A (zh) * | 2023-07-12 | 2023-11-03 | 中国科学院力学研究所 | 基于能量的无动力滑降无人机着陆下滑道动态生成方法 |
CN116991171B (zh) * | 2023-07-12 | 2024-04-16 | 中国科学院力学研究所 | 基于能量的无动力滑降无人机着陆下滑道动态生成方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107941087B (zh) | 2019-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107941087B (zh) | 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法 | |
CN111306989B (zh) | 一种基于平稳滑翔弹道解析解的高超声速再入制导方法 | |
CN107966156B (zh) | 一种适用于运载火箭垂直回收段的制导律设计方法 | |
Lee et al. | Polynomial guidance laws considering terminal impact angle and acceleration constraints | |
Lu et al. | Adaptive terminal guidance for hypervelocity impact in specified direction | |
CN103090728B (zh) | 一种基于滑模控制的带末角约束制导方法 | |
CN109740198B (zh) | 一种基于解析预测的滑翔飞行器三维再入制导方法 | |
CN103245257B (zh) | 基于Bezier曲线的多约束飞行器导引方法 | |
CN104656450B (zh) | 一种高超声速飞行器平稳滑翔再入弹道设计方法 | |
CN106227972A (zh) | 一种高超声速飞行器平稳滑翔弹道的优化方法 | |
CN103486905B (zh) | 一种再入飞行器末制导交班条件确定方法 | |
CN104199303B (zh) | 一种基于矢量场制导的平流层卫星平面路径跟踪控制方法 | |
CN106444430A (zh) | 运载火箭一子级再入控制系统及方法、仿真系统及方法 | |
Wang et al. | Terminal guidance for a hypersonic vehicle with impact time control | |
CN106444822A (zh) | 一种基于空间矢量场制导的平流层飞艇路径跟踪控制方法 | |
CN113031642B (zh) | 动态禁飞区约束的高超声速飞行器轨迹规划方法和系统 | |
CN110764523B (zh) | 基于抗饱和平滑变换的比例积分前置导引攻击目标的方法 | |
CN112256061A (zh) | 复杂环境及任务约束下的高超声速飞行器再入制导方法 | |
CN109240323A (zh) | 一种实时解析构造的空天飞行器再入制导方法 | |
Leavitt et al. | Performance of evolved acceleration guidance logic for entry (EAGLE) | |
US10480904B2 (en) | Gbias for rate based autopilot | |
CN114035616B (zh) | 一种飞行器对移动目标打击控制方法及系统 | |
Chen et al. | Impact time and angle constrained guidance via range‐based line‐of‐sight shaping | |
Yang | Analysis of optimal midcourse guidance law | |
Dhabale et al. | Impact angle constraint guidance law using cubic splines for intercepting stationary targets |
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 |