CN105676638A - 平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法 - Google Patents
平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法 Download PDFInfo
- Publication number
- CN105676638A CN105676638A CN201610016198.1A CN201610016198A CN105676638A CN 105676638 A CN105676638 A CN 105676638A CN 201610016198 A CN201610016198 A CN 201610016198A CN 105676638 A CN105676638 A CN 105676638A
- Authority
- CN
- China
- Prior art keywords
- glide
- angle
- height
- trajectory
- gliding
- 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
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Feedback Control In General (AREA)
Abstract
一种平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法,它包括以下步骤:1:再入滑翔纵向动力学建模;2:跳跃滑翔弹道高度振幅解析解;3:含高度振幅反馈的弹道阻尼控制方法设计;4:平稳滑翔参考攻角和倾侧角设计;5:准自然频率跳跃滑翔模态切换方案设计;最终获得了平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法。本发明获得了解析的再入弹道高度振幅解析表达式,具有很高的求解精度,实现了对再入弹道高度振幅的精确控制;控制方法简单效率高。
Description
技术领域
本发明提供一种平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法,属于航天技术、武器技术领域。
背景技术
近年来,随着反导技术的不断进步,机动突防设计成为任务规划中不可或缺的一环。相对于传统的弹道导弹来说,助推/滑翔式高超声速飞行器拥有更强的机动能力,通过倾侧反转能够进行大范围的横向机动;而纵向机动则主要有平稳机动滑翔和跳跃机动滑翔两种方案。尽管平稳机动滑翔具有高度变化平缓、热流密度和动压峰值小、攻角和倾侧角曲线光滑等优点,且便于确定再入弹道在约束走廊中的位置,但跳跃滑翔弹道具有更大的纵向机动幅度,若能解决跳跃机动过程中的过程约束问题,则其在机动突防时更具优势。
发明内容
(一)目的:本发明的目的是提供一种平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法,它利用再入滑翔弹道高度跳跃滑翔弹道的高度振幅预测,通过阻尼反馈对弹道高度振幅进行精确控制,最后给出了一种平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道快速规划方法。
(二)技术方案:本发明一种平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法,包括以下几个步骤:
步骤1:再入滑翔纵向动力学建模
忽略地球自转,高超声速飞行器再入滑翔弹道的纵向运动方程,如下所示:
上式中,h、V、γ分别为再入飞行器的高度、速度、弹道倾角;分别为再入飞行器的高度、速度和弹道倾角对时间的导数;r为从地心至飞行器的径向距离,与高度的关系为h=r-R0,其中R0为地球半径;g为重力加速度。
步骤2:跳跃滑翔弹道高度振幅解析解
定义平稳滑翔高度增量和弹道倾角增量分别如下:
Δh=h-hsg
Δγ=γ-γsg
上式中,Δh和Δγ分别为平稳滑翔高度增量和弹道倾角增量;hsg和γsg分别为平稳滑翔高度和弹道倾角;由于再入过程中弹道倾角(包括平稳滑翔弹道倾角)为小量,因此可假设sinγsg≈γsg,sinγ≈γ,cosγsg≈cosγ≈1,从而可得平稳滑翔高度增量和弹道倾角增量的关系如下:
上式中,kL1为平稳滑翔高度相关的量;βh为指数大气模型常数。对上式进行积分可得:
上式中,Δh0和Δγ0分别为定义平稳滑翔高度增量和弹道倾角增量的初值。上式中若取Δγ=0,则可求得高度增量的最大值Δhmax和最小值Δhmin,从而获得跳跃滑翔弹道高度振幅解析解。
步骤3:含高度振幅反馈的弹道阻尼控制方法设计
在获得高度振幅的解析解之后,下一步的工作则是对高度振幅进行精确控制。由阻尼影响分析可知,当ζc>-0.028时,高度振幅随时间逐渐减小;当ζc<-0.028时,高度振幅随时间逐渐增大;当ζc=-0.028时,高度振幅几乎不变。因此,可构造如下反馈:
ζc=-0.028+kh(Δhmin|-Aneed)
上式中,ζc为由高度振幅反馈获得的滑翔弹道阻尼;Aneed为期望高度振幅;|Δhmin|为实际预测振幅;kh为高度相关常数,单位为m-1。为了防止振幅剧烈变化,同时避免过阻尼,还需将上获得的阻尼进行限幅,当ζc大于0.8时,取ζc=0.8;当ζc小于-1时,取ζc=-1。将上带入改进弹道阻尼控制方法后可得高度振幅反馈的弹道阻尼控制方法如下:
上式中,为参考纵向升力系数分量;CL1为纵向升力系数分量;Kref为由参考攻角和参考倾侧角决定的纵向升阻比。
步骤4:平稳滑翔参考攻角和倾侧角设计
为了使得再入滑翔弹道满足再入过程约束与终端约束,首先需要设计平稳滑翔所需的参考攻角曲线αref和倾侧角曲线σref;参考攻角曲线的设计原则是既能满足飞行任务所需的射程覆盖范围,又留有启动准自然频率跳跃滑翔所需的攻角调整余量;在获得参考攻角曲线αref后,利用平衡滑翔条件获得的倾侧角走廊的上边界(由最大热流密度约束、最大动压约束和最大过载约束决定),只要所选择的参考倾侧角曲线位于上述走廊内,则所获得的滑翔弹道满足过程约束;调整参考倾侧角绝对值的大小可改变再入滑翔飞行距离,而调整倾侧角的反转点tc1和tc2则可使得滑翔弹道满足终端经纬度约束。
步骤5:准自然频率跳跃滑翔模态切换方案设计
当飞行器遭遇图5所示的不可从顶上越过、且无法侧向绕飞的禁飞区时,需要采用准自然频率跳跃滑翔进行机动突防;若滑翔弹道与禁飞区中心的距离Sradius小于禁飞区半径Rnofly,且飞行器与禁飞区边界的距离Sdist小于临界距离Slim,则需切换到准自然频率跳跃滑翔模态。此时,跳跃滑翔参考攻角为平稳滑翔参考攻角曲线αref的基础上额外增加Δα,该Δα为启动准自然频率跳跃滑翔的攻角增量,由期望高度振幅确定。此外,在平稳滑翔模态下,飞行器的滑翔弹道阻尼取0.8,可较好的抑制弹道振荡;而在准自然频率跳跃滑翔模态下,飞行器的滑翔弹道阻尼高度振幅反馈确定。
通过上述5个步骤,按照图1给出的流程,最终获得了本发明所述的一种平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法。
(三)本发明的优点及功效在于:
(1)根据给定攻角和倾侧角的平稳滑翔高度增量与弹道倾角增量的微分方程,获得了解析的再入弹道高度振幅解析表达式,具有很高的求解精度;
(2)含高度振幅反馈的弹道阻尼控制方法设计,实现了对再入弹道高度振幅的精确控制;
(3)给出了一种平稳滑翔/准自然频率跳跃滑翔机动突防弹道快速规划方法,具有规划效率高、控制方法简单等特点。
附图说明
图1本发明所述方法流程示意图。
图2a是自然振荡纵向弹道与平稳滑翔纵向弹道示意图。
图2b是自然阻尼振荡弹道高度增量振幅预测示意图。
图3a是含高度振幅反馈的弹道阻尼控制方法获得的纵向弹道曲线示意图。
图3b是含高度振幅反馈的弹道阻尼控制方法获得的高度增量曲线示意图。
图3c是含高度振幅反馈的弹道阻尼控制方法获得的攻角曲线示意图。
图3d是含高度振幅反馈的弹道阻尼控制方法获得的道阻尼曲线示意图。
图4a是不同期望高度振幅下的高度增量曲线示意图。
图4b是不同期望高度振幅下的高度增量曲线局部放大示意图。
图4c是不同期望高度振幅下的纵向弹道曲线示意图。
图4d是不同期望高度振幅下的攻角曲线示意图。
图5是不可从顶上跃过禁飞区示意图。
图6是滑翔弹道与禁飞区关系示意图。
图7是平稳滑翔/准自然频率跳跃滑组合机动突防弹道纵向弹道曲线示意图。
图8是平稳滑翔/准自然频率跳跃滑组合机动突防弹道横向弹道曲线示意图。
图9是平稳滑翔/准自然频率跳跃滑组合机动突防弹道速度曲线示意图。
图10是平稳滑翔/准自然频率跳跃滑组合机动突防弹道动压曲线示意图。
图11是平稳滑翔/准自然频率跳跃滑组合机动突防弹道热流密度曲线示意图。
图12是平稳滑翔/准自然频率跳跃滑组合机动突防弹道过载曲线示意图。
图13是平稳滑翔/准自然频率跳跃滑组合机动突防弹道倾侧角曲线示意图。
图14是平稳滑翔/准自然频率跳跃滑组合机动突防弹道攻角曲线示意图。
图15是平稳滑翔/准自然频率跳跃滑组合机动突防弹道阻尼曲线示意图。
图中符号、代号说明如下:tc1和tc2分别为倾侧角反转时刻;αref和σref分别为平稳滑翔的参考攻角和参考倾侧角;Δα为启动准自然频率跳跃滑翔的攻角增量;Δhmax和Δhmin分别为预测高度增量的最大值和最小值;ζc为实际弹道阻尼;Sradius为滑翔弹道与禁飞区中心的距离;Rnofly为禁飞区半径;θnofly和φnofly分别为禁飞区中心的经度和纬度。
具体实施方式
下面将结合附图和实施案例对本发明作进一步的详细说明。
针对平稳滑翔/准自然频率跳跃滑翔组合机动突防的弹道规划问题,本发明首先建立再入滑翔纵向动力学模型;然后对跳跃滑翔弹道的高度振幅进行解析预测;并设计了高度振幅反馈的弹道阻尼控制方法,实现了对滑翔弹道高度振幅的精确控制;最后,给出了平稳滑翔参考攻角和倾侧角设计方法和准自然频率跳跃滑翔模态切换方案。
本发明一种平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法,见图1所示,包括以下几个步骤:
步骤1:再入滑翔纵向动力学建模。
忽略地球自转,高超声速飞行器再入滑翔弹道的纵向运动方程如下所示:
式中,h、V、γ分别为再入飞行器的高度、速度、弹道倾角;分别为再入飞行器的高度、速度和弹道倾角对时间的导数;r为从地心至飞行器的径向距离,与高度的关系为h=r-R0,其中R0为地球半径;g为重力加速度。aL1和aD分别为升力加速度纵向分量和阻力加速度,它们的表达式为,
式中,CL1为纵向升力系数分量,CL1=CLcosσ;CL和CD分别为升力系数和阻力系数;σ为倾侧角;m为飞行器质量;S为气动参考面积;ρ为大气密度,通常采用指数大气模型,如下所示,
式中,ρsea为海平面大气密度;βh为指数大气模型常数,通常取βh=1/7200。
步骤2:跳跃滑翔弹道高度振幅解析预测。
定义平稳滑翔高度增量和弹道倾角增量分别如下:
Δh=h-hsg(6)
Δγ=γ-γsg(7)
上式中,Δh和Δγ分别为平稳滑翔高度增量和弹道倾角增量;hsg和γsg分别为平稳滑翔高度和弹道倾角,满足如下关系:
上式中,和分别为平稳滑翔高度和弹道倾角对时间的导数;为参考升力系数纵向分量,由参考攻角和倾侧角决定。
再入过程中,弹道倾角(包括平稳滑翔弹道倾角)为小量,因此可假设sinγsg≈γsg,sinγ≈γ,cosγsg≈cosγ≈1,从而由式(1)、(3)、(6)、(7)、(8)和(9)可得,
上式中,kL1为平稳滑翔高度相关的量,表达式为
由式(10)和(11)可得,
对上式进行积分,
可得,
上式中,Δh0和Δγ0分别为定义平稳滑翔高度增量和弹道倾角增量的初值。式(14)给出了振荡弹道的高度增量和弹道倾角增量之间的关系,若取Δh=0,可求得弹道倾角增量的最大值Δγmax;若取Δγ=0,则可求得高度增量的最大值Δhmax和最小值Δhmin。需要指出的是,通常情况下,Δhmax的绝对值并不等于Δhmin的绝对值,考虑到热流密度或动压的峰值通常出现在高度增量取最小值的时候,因此定义|Δhmin|为高度振幅。
为了检验式(14)的精度,取αref=15deg、σref=0deg,初始高度增量Δh0=5km,初始弹道倾角增量Δγ0=0deg,阻尼分别取自然阻尼和-0.03,得到的结果如图2a-图2b所示。由图可以看出,采用式(14)能够较为准确的预测的高度振幅,但由于线性化误差及速度衰减的影响,整个滑翔段预测的高度振幅曲线存在着小幅度的扰动,扰动的幅度约为100m。
步骤3:含高度振幅反馈的弹道阻尼控制方法
在获得高度振幅的解析解之后,下一步的工作则是对高度振幅进行精确控制。由阻尼影响分析可知,当ζc>-0.028时,高度振幅随时间逐渐减小;当ζc<-0.028时,高度振幅随时间逐渐增大;当ζc=-0.028时,高度振幅几乎不变。因此,可构造如下反馈:
ζc=-0.028+kh(Δhmin|-Aneed)(15)
上式中,ζc为由高度振幅反馈获得的滑翔弹道阻尼;Aneed为期望高度振幅;|Δhmin|为实际预测振幅;kh为高度相关常数,单位为m-1。为了防止振幅剧烈变化,同时避免过阻尼,还需将式(15)获得的阻尼进行限幅,当ζc大于0.8时,取ζc=0.8;当ζc小于-1时,取ζc=-1。将式(15)带入改进弹道阻尼控制方法后可得高度振幅反馈的弹道阻尼控制方法如下:
上式中,为参考纵向升力系数;CL1为实际纵向升力系数;Kref为由参考攻角和参考倾侧角决定的纵向升阻比。
图3a-图3d给出了不同初始弹道倾角增量下采用含高度振幅反馈的弹道阻尼控制方法获得的滑翔弹道,图中参考攻角和倾侧角分别取15deg和0deg,期望高度振幅为2km,kh=0.001m-1。由图可以看出,不论初始状态如何,上述方法均能很快使得滑翔弹道的高度振幅收敛到2km,并且仅在初始动态调整段需要较大的攻角调整量。图4a-图4d则给出了采用含高度振幅反馈的弹道阻尼控制方法获得不同期望高度振幅的结果,可以看出,采用该方法获得的滑翔弹道高度振幅误差在2%以内。
步骤4:平稳滑翔参考攻角和倾侧角设计
为了使得再入滑翔弹道满足再入过程约束与终端约束,首先需要设计平稳滑翔所需的参考攻角曲线αref和倾侧角曲线σref。参考攻角曲线的设计原则是既能满足飞行任务所需的射程覆盖范围,又留有启动准自然频率跳跃滑翔所需的攻角调整余量。在获得参考攻角曲线αref后,利用平衡滑翔条件,可将再入滑翔过程中的最大热流密度约束、最大动压约束和最大过载约束转化为倾侧角走廊的上边界,只要所选择的参考倾侧角曲线位于上述走廊内,则所获得的滑翔弹道满足过程约束。调整参考倾侧角绝对值的大小可改变再入滑翔飞行距离,而调整倾侧角的反转点tc1和tc2则可使得滑翔弹道满足终端经纬度约束。
步骤5:准自然频率跳跃滑翔模态切换方案设计
当飞行器遭遇图5所示的不可从顶上越过、且无法侧向绕飞的禁飞区时,需要采用准自然频率跳跃滑翔进行机动突防。规划过程中判断滑翔弹道是否穿过禁飞区的标准为滑翔弹道与禁飞区中心的距离Sradius是否小于禁飞区半径Rnofly,其中Sradius可由下式预测:
Sradius=Srangesin(ψnofly-ψ)(17)
上式中,ψ为飞行器的实际弹道偏角;Sradius为滑翔弹道与禁飞区中心的距离;ψnofly为禁飞区中心的方向角,Srange为飞行器与禁飞区中心的球面距离,两者定义如图6所示,表达式如下:
Srange=arccos[sinφsinφnofly+cosφcosφnoflycos(θ-θnofly)](18)
上式中,θ和φ分别为飞行器的当前经度和纬度;θnofly和φnofly分别为禁飞区中心的经度和纬度。若Sradius>Rnofly(Rnofly为禁飞区半径),表示再入滑翔弹道不会进入禁飞区,无需进行准自然频率跳跃机动突防;Sradius≤Rnofly则表示再入弹道将进入禁飞区,需要适时启动准自然频率跳跃滑翔模态,启动的标准为飞行器与禁飞区边界的距离Sdist小于临界距离Slim。
为了使飞行器快速从平稳滑翔模态转换到给定高度振幅的准自然频率跳跃滑翔模态,可在平稳滑翔参考攻角曲线αref的基础上额外增加Δα,其中Δα为期望高度振幅所对应的攻角增量,Δα+αref可由下面的升力系数求解获得。
上式中,和CLz分别为αref和Δα+αref所对应的升力系数。
此外,在平稳滑翔模态下,飞行器的滑翔弹道阻尼取0.8,可较好的抑制弹道振荡;而在准自然频率跳跃滑翔模态下,飞行器的滑翔弹道阻尼则由式(15)确定。
通过上述5个步骤,按照图1给出的流程,最终获得了种平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法。
实施案例
设初始弹道偏角、初始经度、初始纬度、初始速度和终端速度分别取90deg、0deg、0deg、7000m/s和2000m/s,目标经纬度如表1所示,禁飞区1和禁飞区2的中心经纬度分别为(30deg,0deg)和(90deg,0deg),半径均为1000km。
表1平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道目标经纬度设定
则采用上述规划流程获得的结果如图7至图15所示。由图7至图8可以看出,规划获得的弹道满足终端位置约束,并且在禁飞区内保持振荡突防形态。由图9至图12可以看出,规划获得的弹道满足终端速度约束、动压约束、热流密度约束、过载约束。由图13和图14可以看出,规划获得的弹道倾侧角曲线较为光滑,而攻角曲线则变化较大,这是由于本文中采用攻角来调整纵向升力系数分量。若采用倾侧角来调整纵向升力系数分量,则攻角曲线会更加光滑,而倾侧角曲线则变化较大。由图15则可以看出本文提出的规划方法设计的阻尼曲线比较简单。
Claims (1)
1.一种平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法,其特征在于:它包括以下几个步骤:
步骤1:再入滑翔纵向动力学建模
忽略地球自转,高超声速飞行器再入滑翔弹道的纵向运动方程,如下所示:
上式中,h、V、γ分别为再入飞行器的高度、速度、弹道倾角;分别为再入飞行器的高度、速度和弹道倾角对时间的导数;r为从地心至飞行器的径向距离,与高度的关系为h=r-R0,其中R0为地球半径;g为重力加速度;
步骤2:跳跃滑翔弹道高度振幅解析解
定义平稳滑翔高度增量和弹道倾角增量分别如下:
Δh=h-hsg···············
Δγ=γ-γsg···············
上式中,Δh和Δγ分别为平稳滑翔高度增量和弹道倾角增量;hsg和γsg分别为平稳滑翔高度和弹道倾角;由于再入过程中弹道倾角包括平稳滑翔弹道倾角为小量,因此假设sinγsg≈γsg,sinγ≈γ,cosγsg≈cosγ≈1,从而得平稳滑翔高度增量和弹道倾角增量的关系如下:
上式中,kL1为平稳滑翔高度相关的量;βh为指数大气模型常数;对上式进行积分得:
上式中,Δh0和Δγ0分别为定义平稳滑翔高度增量和弹道倾角增量的初值;上式中若取Δγ=0,则求得高度增量的最大值Δhmax和最小值Δhmin,从而获得跳跃滑翔弹道高度振幅解析解;
步骤3:含高度振幅反馈的弹道阻尼控制方法设计
在获得高度振幅的解析解之后,下一步的工作则是对高度振幅进行精确控制;由阻尼影响分析可知,当ζc>-0.028时,高度振幅随时间逐渐减小;当ζc<-0.028时,高度振幅随时间逐渐增大;当ζc=-0.028时,高度振幅几乎不变,因此,能构造如下反馈:
ζc=-0.028+kh(|Δhmin|-Aneed)···········
上式中,ζc为由高度振幅反馈获得的滑翔弹道阻尼;Aneed为期望高度振幅;|Δhmin|为实际预测振幅;kh为高度相关常数,单位为m-1;为了防止振幅剧烈变化,同时避免过阻尼,还需将上获得的阻尼进行限幅,当ζc大于0.8时,取ζc=0.8;当ζc小于-1时,取ζc=-1;将上带入改进弹道阻尼控制方法后得高度振幅反馈的弹道阻尼控制方法如下:
上式中,为参考纵向升力系数分量;CL1为纵向升力系数分量;Kref为由参考攻角和参考倾侧角决定的纵向升阻比;
步骤4:平稳滑翔参考攻角和倾侧角设计
为了使得再入滑翔弹道满足再入过程约束与终端约束,首先需要设计平稳滑翔所需的参考攻角曲线αref和倾侧角曲线σref;参考攻角曲线的设计原则是既能满足飞行任务所需的射程覆盖范围,又留有启动准自然频率跳跃滑翔所需的攻角调整余量;在获得参考攻角曲线αref后,利用平衡滑翔条件获得的倾侧角走廊的上边界,只要所选择的参考倾侧角曲线位于上述走廊内,则所获得的滑翔弹道满足过程约束;调整参考倾侧角绝对值的大小能改变再入滑翔飞行距离,而调整倾侧角的反转点tc1和tc2则使得滑翔弹道满足终端经纬度约束;
步骤5:准自然频率跳跃滑翔模态切换方案设计
当飞行器遭遇不可从顶上越过、且无法侧向绕飞的禁飞区时,需要采用准自然频率跳跃滑翔进行机动突防;若滑翔弹道与禁飞区中心的距离Sradius小于禁飞区半径Rnofly,且飞行器与禁飞区边界的距离Sdist小于临界距离Slim,则需切换到准自然频率跳跃滑翔模态;此时,跳跃滑翔参考攻角为平稳滑翔参考攻角曲线αref的基础上额外增加Δα,该Δα为启动准自然频率跳跃滑翔的攻角增量,由期望高度振幅确定;此外,在平稳滑翔模态下,飞行器的滑翔弹道阻尼取0.8,能较好的抑制弹道振荡;而在准自然频率跳跃滑翔模态下,飞行器的滑翔弹道阻尼高度振幅反馈确定;
通过上述5个步骤,最终获得了本发明所述的一种平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610016198.1A CN105676638B (zh) | 2016-01-11 | 2016-01-11 | 平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610016198.1A CN105676638B (zh) | 2016-01-11 | 2016-01-11 | 平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105676638A true CN105676638A (zh) | 2016-06-15 |
CN105676638B CN105676638B (zh) | 2018-06-29 |
Family
ID=56300022
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610016198.1A Active CN105676638B (zh) | 2016-01-11 | 2016-01-11 | 平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105676638B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106294280A (zh) * | 2016-08-22 | 2017-01-04 | 麻毅威 | 一种弹道规划方法 |
CN107621198A (zh) * | 2017-08-28 | 2018-01-23 | 北京航空航天大学 | 一种绕多禁飞区的多步决策弹道规划方法 |
CN109145451A (zh) * | 2018-08-22 | 2019-01-04 | 哈尔滨工业大学 | 一种高速滑翔飞行器的运动行为识别与航迹估计方法 |
CN109190248A (zh) * | 2018-09-03 | 2019-01-11 | 中国运载火箭技术研究院 | 一种用于滑翔飞行器的滑翔射程解析方法及解析系统 |
CN110147521A (zh) * | 2019-04-25 | 2019-08-20 | 北京航空航天大学 | 一种高超声速飞行器跳跃滑翔弹道解析求解方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5402965A (en) * | 1993-09-20 | 1995-04-04 | Rockwell International Corporation | Reusable flyback satellite |
US20090112535A1 (en) * | 2007-10-26 | 2009-04-30 | Sensis Corporation | Method of integrating point mass equations to include vertical and horizontal profiles |
CN103176476A (zh) * | 2013-03-08 | 2013-06-26 | 北京航空航天大学 | 一种滑翔式无人机自主进场航路规划方法 |
CN103838914A (zh) * | 2013-12-30 | 2014-06-04 | 北京航空航天大学 | 一种高超声速飞行器滑翔段弹道解析求解方法 |
CN104176268A (zh) * | 2014-08-06 | 2014-12-03 | 北京航空航天大学 | 一种滑翔飞行弹道阻尼控制方法 |
CN104392047A (zh) * | 2014-11-25 | 2015-03-04 | 北京航空航天大学 | 一种基于平稳滑翔弹道解析解的快速弹道规划方法 |
CN104656450A (zh) * | 2015-01-20 | 2015-05-27 | 北京航空航天大学 | 一种高超声速飞行器平稳滑翔再入弹道设计方法 |
CN104731104A (zh) * | 2015-03-09 | 2015-06-24 | 北京航天自动控制研究所 | 一种高超声速飞行器滑翔飞行段的纵向制导方法 |
CN104778376A (zh) * | 2015-05-04 | 2015-07-15 | 哈尔滨工业大学 | 一种临近空间高超声速滑翔弹头跳跃弹道预测方法 |
CN104965418A (zh) * | 2015-06-01 | 2015-10-07 | 北京航空航天大学 | 一种基于弹道阻尼控制和热流解析预测的引入段制导方法 |
CN105022403A (zh) * | 2015-05-11 | 2015-11-04 | 北京航天自动控制研究所 | 滑翔飞行器的纵向轨迹控制增益的确定方法 |
-
2016
- 2016-01-11 CN CN201610016198.1A patent/CN105676638B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5402965A (en) * | 1993-09-20 | 1995-04-04 | Rockwell International Corporation | Reusable flyback satellite |
US20090112535A1 (en) * | 2007-10-26 | 2009-04-30 | Sensis Corporation | Method of integrating point mass equations to include vertical and horizontal profiles |
CN103176476A (zh) * | 2013-03-08 | 2013-06-26 | 北京航空航天大学 | 一种滑翔式无人机自主进场航路规划方法 |
CN103838914A (zh) * | 2013-12-30 | 2014-06-04 | 北京航空航天大学 | 一种高超声速飞行器滑翔段弹道解析求解方法 |
CN104176268A (zh) * | 2014-08-06 | 2014-12-03 | 北京航空航天大学 | 一种滑翔飞行弹道阻尼控制方法 |
CN104392047A (zh) * | 2014-11-25 | 2015-03-04 | 北京航空航天大学 | 一种基于平稳滑翔弹道解析解的快速弹道规划方法 |
CN104656450A (zh) * | 2015-01-20 | 2015-05-27 | 北京航空航天大学 | 一种高超声速飞行器平稳滑翔再入弹道设计方法 |
CN104731104A (zh) * | 2015-03-09 | 2015-06-24 | 北京航天自动控制研究所 | 一种高超声速飞行器滑翔飞行段的纵向制导方法 |
CN104778376A (zh) * | 2015-05-04 | 2015-07-15 | 哈尔滨工业大学 | 一种临近空间高超声速滑翔弹头跳跃弹道预测方法 |
CN105022403A (zh) * | 2015-05-11 | 2015-11-04 | 北京航天自动控制研究所 | 滑翔飞行器的纵向轨迹控制增益的确定方法 |
CN104965418A (zh) * | 2015-06-01 | 2015-10-07 | 北京航空航天大学 | 一种基于弹道阻尼控制和热流解析预测的引入段制导方法 |
Non-Patent Citations (2)
Title |
---|
WENBIN YU 等: "Trajectory-Shaping Guidance with final speed and load factor constraints", 《ISA TRANSACTIONS》 * |
胡锦川 等: "高超声速飞行器平稳滑翔弹道设计方法", 《北京航空航天大学学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106294280A (zh) * | 2016-08-22 | 2017-01-04 | 麻毅威 | 一种弹道规划方法 |
CN106294280B (zh) * | 2016-08-22 | 2018-11-20 | 麻毅威 | 一种弹道规划方法 |
CN107621198A (zh) * | 2017-08-28 | 2018-01-23 | 北京航空航天大学 | 一种绕多禁飞区的多步决策弹道规划方法 |
CN107621198B (zh) * | 2017-08-28 | 2019-04-12 | 北京航空航天大学 | 一种绕多禁飞区的多步决策弹道规划方法 |
CN109145451A (zh) * | 2018-08-22 | 2019-01-04 | 哈尔滨工业大学 | 一种高速滑翔飞行器的运动行为识别与航迹估计方法 |
CN109190248A (zh) * | 2018-09-03 | 2019-01-11 | 中国运载火箭技术研究院 | 一种用于滑翔飞行器的滑翔射程解析方法及解析系统 |
CN109190248B (zh) * | 2018-09-03 | 2023-07-18 | 中国运载火箭技术研究院 | 一种用于滑翔飞行器的滑翔射程解析方法及解析系统 |
CN110147521A (zh) * | 2019-04-25 | 2019-08-20 | 北京航空航天大学 | 一种高超声速飞行器跳跃滑翔弹道解析求解方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105676638B (zh) | 2018-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105676638A (zh) | 平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法 | |
CN108036676B (zh) | 一种基于三维再入弹道解析解的全射向自主再入制导方法 | |
CN109508030B (zh) | 一种考虑多禁飞区约束的协同解析再入制导方法 | |
CN107941087A (zh) | 一种基于阻力剖面的高升阻比高超平稳滑翔再入制导方法 | |
US11286065B2 (en) | Method for designing reentry trajectory based on flight path angle planning | |
CN108387140B (zh) | 一种考虑多个禁飞区约束的解析再入制导方法 | |
CN109740198B (zh) | 一种基于解析预测的滑翔飞行器三维再入制导方法 | |
CN104656450B (zh) | 一种高超声速飞行器平稳滑翔再入弹道设计方法 | |
CN104648695B (zh) | 一种基于倾侧角可用性的再入走廊最优规划方法 | |
CN111306989A (zh) | 一种基于平稳滑翔弹道解析解的高超声速再入制导方法 | |
CN106371312B (zh) | 基于模糊控制器的升力式再入预测-校正制导方法 | |
CN106643341B (zh) | 基于准平衡滑翔原理的力热控制耦合设计方法 | |
CN106054604B (zh) | 基于模型预测控制理论的再入飞行器鲁棒最优制导方法 | |
CN112947573B (zh) | 终端时间约束下的高超声速飞行器再入制导方法 | |
CN102139768A (zh) | 一种亚轨道飞行器再入飞行的攻角制导方法 | |
CN105550402B (zh) | 一种基于攻角或倾侧角变频的高超平稳机动滑翔弹道设计方法 | |
Zhu et al. | Impact time and angle control guidance independent of time-to-go prediction | |
US11845573B2 (en) | Configuration construction and attitude control method for pyramid deorbit sail | |
Zimmermann et al. | Optimization of the tether-assisted return mission of a guided re-entry capsule | |
CN110733670B (zh) | 一种短航程低过载的再入轨迹设计方法 | |
CN107121929A (zh) | 基于线性协方差模型预测控制的鲁棒再入制导方法 | |
CN110147521A (zh) | 一种高超声速飞行器跳跃滑翔弹道解析求解方法 | |
CN104965418B (zh) | 一种基于弹道阻尼控制和热流解析预测的引入段制导方法 | |
CN104850129B (zh) | 一种跳跃式再入的射向预偏置横向制导方法 | |
CN109240335A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |