一种高超声速飞行器弹道的仿真方法
一、相关研究
临近空间高超声速飞行器(Near Space Hypersonic Vehicle, NSHV)飞行速度快、突防能力强、打击距离远,主要有 2 类: 一类是由火箭助推至大气层高空后,无动力再入大气层而长距离滑翔的升力式再入滑翔弹,典型的有 HTV-2 和AHW; 一类是机载平台投放后,由吸气式超燃冲压发动机连续启停助推的乘波体滑翔弹,典型的有X-43A、X-51A。
对高超声速飞行器的弹道特性进行仿真,可为高超声速飞行器的探测与拦截问题研究提供支撑;也可用于各种模拟系统,为其增加仿真对象和仿真样式。
文献[1]对X-51的试验过程与弹道特性进行了分析,但并未给出相应的计算模型。文献[2]将飞行器的飞行过程划分为高空巡航段和俯冲攻击段,在高空巡航段设计了常值稳态式、非周期性滑跃式和周期性跳跃式3种控制方案,在俯冲攻击段则采用比例导引法进行控制。文献[3]指出,高超声速滑翔飞行器具有两种典型的飞行弹道模式:一种称为平衡滑翔,最早由德国科学家 Saenger 于 20 世纪 30年代在研究高速升力式飞行时提出,其弹道特点是几乎没有波动的平坦滑翔下降弹道,近年来研究中把倾侧角时变条件的平衡滑翔称为准平衡滑翔;另一种是 Saenger 所设计的“银鸟”飞行器采用的滑翔跳跃弹道,即具有一定跳跃波动幅度的滑翔轨迹。该文献对纵向弹道特性进行了较为深入的分析。文献[4] 以升力式再入滑翔弹为对象,研究了控制量和初始运动状态对弹道轨迹的具体影响。文献[5]设计了以高度为节点约束、攻角为控制量的周期性跳跃巡航式制导飞行方案并生成弹道轨迹。文献[6]给出了常攻角飞行、常升阻比飞行、最优升阻比飞行、平衡滑翔飞行和横向机动的动力学模型。文献[7]给出了基于椭球的弹道仿真动力学模型。文献[8]则专门针对导弹爬升段建立了仿真模型。在飞行过程中为实现机动突防、过航路点或躲避禁飞区的目的,往往会选择适当的机动模式。高超声速目标典型的机动模式有纵向跳跃滑翔机动、侧向摆动突防机动以及螺旋机动等。文献[9]探讨了飞行器弹道中的机动特性,主要是目标的纵向和侧向机动能力。
二、基本思路
结合利用多个文献中给定模型,设计并实现了一个高超声速飞行器弹道仿真算法,并应用到基于OSGEarth的态势中。
在美飞行试验中,X-51A 由 B-52 轰炸机携带升空,自母机投放后经火箭推进至超燃冲压发动机的工作高度及飞行马赫数,然后超燃冲压发动机点火,将飞行器由马赫数 4.5 加速到马赫数 6 的巡航速度。
因此,主要仿真加速到6马赫之后的跳跃式飞行然后俯冲至目标的过程。即模型的输入输出如下
输入:发射点经度、纬度(高度限定为固定值,如20km);目标点的经度、纬度、高度;
输出:飞行器弹道数据(采样点的位置矢量和速度矢量)。
滑翔和俯冲过程中需要各自独立的计算模型,但都需要求解微分方程,因此采用龙格库塔法求解。为提高方法的通用性,定义方程右值类为
class EquationRightCalculator{
int getNumOfEquation();
void getCurrentValue(double& t,double* y);
virtual void setCurrentValue(double t,double* y);
virtual void calculate(double t,double* y,double* d)=0;
int _numOfEquation;
double* _currentY;
double _currentT;};
定义龙格库塔法类为
class RunghKutta{
void rk1
5971




被折叠的 条评论
为什么被折叠?



