发明名称 平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法
摘要 一种平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法,它包括以下步骤:1:再入滑翔纵向动力学建模;2:跳跃滑翔弹道高度振幅解析解;3:含高度振幅反馈的弹道阻尼控制方法设计;4:平稳滑翔参考攻角和倾侧角设计;5:准自然频率跳跃滑翔模态切换方案设计;最终获得了平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法。本发明获得了解析的再入弹道高度振幅解析表达式,具有很高的求解精度,实现了对再入弹道高度振幅的精确控制;控制方法简单效率高。
申请公布号 CN105676638A 申请公布日期 2016.06.15
申请号 CN201610016198.1 申请日期 2016.01.11
申请人 北京航空航天大学 发明人 陈万春;胡锦川
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 一种平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法,其特征在于:它包括以下几个步骤:步骤1:再入滑翔纵向动力学建模忽略地球自转,高超声速飞行器再入滑翔弹道的纵向运动方程,如下所示:<maths num="0001"><math><![CDATA[<mrow><mover><mi>h</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>V</mi><mi> </mi><mi>s</mi><mi>i</mi><mi>n</mi><mi>&gamma;</mi></mrow>]]></math><img file="FDA0000904234810000011.GIF" wi="228" he="75" /></maths><maths num="0002"><math><![CDATA[<mrow><mover><mi>V</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mo>-</mo><msub><mi>a</mi><mi>D</mi></msub><mo>-</mo><mi>g</mi><mi> </mi><mi>s</mi><mi>i</mi><mi>n</mi><mi>&gamma;</mi></mrow>]]></math><img file="FDA0000904234810000012.GIF" wi="356" he="77" /></maths><maths num="0003"><math><![CDATA[<mrow><mover><mi>&gamma;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>V</mi></mfrac><mo>&lsqb;</mo><msub><mi>a</mi><mrow><mi>L</mi><mn>1</mn></mrow></msub><mo>+</mo><mrow><mo>(</mo><mrow><msup><mi>V</mi><mn>2</mn></msup><mo>/</mo><mi>r</mi><mo>-</mo><mi>g</mi></mrow><mo>)</mo></mrow><mi>c</mi><mi>o</mi><mi>s</mi><mi>&gamma;</mi><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000904234810000013.GIF" wi="613" he="126" /></maths>上式中,h、V、γ分别为再入飞行器的高度、速度、弹道倾角;<img file="FDA0000904234810000014.GIF" wi="238" he="78" />分别为再入飞行器的高度、速度和弹道倾角对时间的导数;r为从地心至飞行器的径向距离,与高度的关系为h=r‑R<sub>0</sub>,其中R<sub>0</sub>为地球半径;g为重力加速度;步骤2:跳跃滑翔弹道高度振幅解析解定义平稳滑翔高度增量和弹道倾角增量分别如下:Δh=h‑h<sub>sg</sub>···············Δγ=γ‑γ<sub>sg</sub>···············上式中,Δh和Δγ分别为平稳滑翔高度增量和弹道倾角增量;h<sub>sg</sub>和γ<sub>sg</sub>分别为平稳滑翔高度和弹道倾角;由于再入过程中弹道倾角包括平稳滑翔弹道倾角为小量,因此假设sinγ<sub>sg</sub>≈γ<sub>sg</sub>,sinγ≈γ,cosγ<sub>sg</sub>≈cosγ≈1,从而得平稳滑翔高度增量和弹道倾角增量的关系如下:<maths num="0004"><math><![CDATA[<mrow><mfrac><mrow><mi>d</mi><mi>&Delta;</mi><mi>h</mi></mrow><mrow><mi>d</mi><mi>&Delta;</mi><mi>&gamma;</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>&Delta;</mi><mi>&gamma;</mi></mrow><mrow><msub><mi>k</mi><mrow><mi>L</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&beta;</mi><mi>h</mi></msub><mi>&Delta;</mi><mi>h</mi></mrow></msup><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mn>.............</mn></mrow>]]></math><img file="FDA0000904234810000015.GIF" wi="1125" he="159" /></maths>上式中,k<sub>L1</sub>为平稳滑翔高度相关的量;β<sub>h</sub>为指数大气模型常数;对上式进行积分得:<maths num="0005"><math><![CDATA[<mrow><mo>-</mo><mfrac><msub><mi>k</mi><mrow><mi>L</mi><mn>1</mn></mrow></msub><msub><mi>&beta;</mi><mi>h</mi></msub></mfrac><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&beta;</mi><mi>h</mi></msub><mi>&Delta;</mi><mi>h</mi></mrow></msup><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&beta;</mi><mi>h</mi></msub><msub><mi>&Delta;h</mi><mn>0</mn></msub></mrow></msup><mo>)</mo></mrow><mo>-</mo><msub><mi>k</mi><mrow><mi>L</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>&Delta;</mi><mi>h</mi><mo>-</mo><msub><mi>&Delta;h</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mfrac><mrow><msup><mi>&Delta;&gamma;</mi><mn>2</mn></msup></mrow><mn>2</mn></mfrac><mo>-</mo><mfrac><mrow><msubsup><mi>&Delta;&gamma;</mi><mn>0</mn><mn>2</mn></msubsup></mrow><mn>2</mn></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000904234810000016.GIF" wi="1100" he="159" /></maths>上式中,Δh<sub>0</sub>和Δγ<sub>0</sub>分别为定义平稳滑翔高度增量和弹道倾角增量的初值;上式中若取Δγ=0,则求得高度增量的最大值Δh<sub>max</sub>和最小值Δh<sub>min</sub>,从而获得跳跃滑翔弹道高度振幅解析解;步骤3:含高度振幅反馈的弹道阻尼控制方法设计在获得高度振幅的解析解之后,下一步的工作则是对高度振幅进行精确控制;由阻尼影响分析可知,当ζ<sub>c</sub>>‑0.028时,高度振幅随时间逐渐减小;当ζ<sub>c</sub><‑0.028时,高度振幅随时间逐渐增大;当ζ<sub>c</sub>=‑0.028时,高度振幅几乎不变,因此,能构造如下反馈:ζ<sub>c</sub>=‑0.028+k<sub>h</sub>(|Δh<sub>min</sub>|‑A<sub>need</sub>)···········上式中,ζ<sub>c</sub>为由高度振幅反馈获得的滑翔弹道阻尼;A<sub>need</sub>为期望高度振幅;|Δh<sub>min</sub>|为实际预测振幅;k<sub>h</sub>为高度相关常数,单位为m<sup>‑1</sup>;为了防止振幅剧烈变化,同时避免过阻尼,还需将上获得的阻尼进行限幅,当ζ<sub>c</sub>大于0.8时,取ζ<sub>c</sub>=0.8;当ζ<sub>c</sub>小于‑1时,取ζ<sub>c</sub>=‑1;将上带入改进弹道阻尼控制方法后得高度振幅反馈的弹道阻尼控制方法如下:<maths num="0006"><math><![CDATA[<mrow><msub><mi>C</mi><mrow><mi>L</mi><mn>1</mn></mrow></msub><mo>=</mo><msubsup><mi>C</mi><mrow><mi>L</mi><mn>1</mn></mrow><mo>*</mo></msubsup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>2</mn><msub><mi>&zeta;</mi><mi>c</mi></msub><msqrt><mrow><mfrac><mrow><msub><mi>&beta;</mi><mi>h</mi></msub><mi>r</mi></mrow><mrow><mi>g</mi><mi>r</mi><mo>/</mo><msup><mi>V</mi><mn>2</mn></msup><mo>-</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mn>1</mn><msub><mi>K</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub></mfrac></mrow></msqrt><mo>)</mo></mrow><msup><mi>e</mi><mrow><msub><mi>&beta;</mi><mi>h</mi></msub><mi>&Delta;</mi><mi>h</mi></mrow></msup><mi>&Delta;</mi><mi>&gamma;</mi><mn>........</mn></mrow>]]></math><img file="FDA0000904234810000021.GIF" wi="1370" he="167" /></maths>上式中,<img file="FDA0000904234810000022.GIF" wi="80" he="71" />为参考纵向升力系数分量;C<sub>L1</sub>为纵向升力系数分量;K<sub>ref</sub>为由参考攻角和参考倾侧角决定的纵向升阻比;步骤4:平稳滑翔参考攻角和倾侧角设计为了使得再入滑翔弹道满足再入过程约束与终端约束,首先需要设计平稳滑翔所需的参考攻角曲线α<sub>ref</sub>和倾侧角曲线σ<sub>ref</sub>;参考攻角曲线的设计原则是既能满足飞行任务所需的射程覆盖范围,又留有启动准自然频率跳跃滑翔所需的攻角调整余量;在获得参考攻角曲线α<sub>ref</sub>后,利用平衡滑翔条件获得的倾侧角走廊的上边界,只要所选择的参考倾侧角曲线位于上述走廊内,则所获得的滑翔弹道满足过程约束;调整参考倾侧角绝对值的大小能改变再入滑翔飞行距离,而调整倾侧角的反转点t<sub>c1</sub>和t<sub>c2</sub>则使得滑翔弹道满足终端经纬度约束;步骤5:准自然频率跳跃滑翔模态切换方案设计当飞行器遭遇不可从顶上越过、且无法侧向绕飞的禁飞区时,需要采用准自然频率跳跃滑翔进行机动突防;若滑翔弹道与禁飞区中心的距离S<sub>radius</sub>小于禁飞区半径R<sub>nofly</sub>,且飞行器与禁飞区边界的距离S<sub>dist</sub>小于临界距离S<sub>lim</sub>,则需切换到准自然频率跳跃滑翔模态;此时,跳跃滑翔参考攻角为平稳滑翔参考攻角曲线α<sub>ref</sub>的基础上额外增加Δα,该Δα为启动准自然频率跳跃滑翔的攻角增量,由期望高度振幅确定;此外,在平稳滑翔模态下,飞行器的滑翔弹道阻尼取0.8,能较好的抑制弹道振荡;而在准自然频率跳跃滑翔模态下,飞行器的滑翔弹道阻尼高度振幅反馈确定;通过上述5个步骤,最终获得了本发明所述的一种平稳滑翔/准自然频率跳跃滑翔组合机动突防弹道规划方法。
地址 100191 北京市海淀区学院路37号