发明名称 近空间高超声速飞行器运动控制仿真方法
摘要 本发明公开一种近空间高超声速飞行器的控制仿真方法,包括如下步骤:(1)建立飞行器在变化风场条件下的运动模型;(2)建立所述近空间高超声速飞行器的发动机推力模型;(3)建立与运动模型和推力模型相关的气动参数库和发动机参数库;(4)根据不同的任务要求,基于运动模型设计姿态稳定跟踪控制律及重心运动控制律;(5)实时计算当前状态下的大气数据和风场数据,并据此以及所设计的控制律计算飞行器的各个状态,重复此步骤直至仿真结束。此种仿真方法在充分剖析近空间高超声速飞行器的动力学和运动学规律的基础上,结合发动机推力模型和大气风场环境的实际情况,提出一整套在变化风场环境下的近空间高超声速飞行器的运动模型。本发明还公开一种基于该飞行器的控制仿真平台。
申请公布号 CN102073755B 申请公布日期 2013.04.10
申请号 CN201010537312.8 申请日期 2010.11.10
申请人 南京航空航天大学 发明人 程路;都延丽;姜长生;薛雅丽;蒲明;傅健;文杰
分类号 G06F17/50(2006.01)I;B64G7/00(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 许方
主权项 1.一种近空间高超声速飞行器运动控制仿真方法,其特征在于包括如下步骤:(1)建立飞行器在变化风场条件下的运动模型:<img file="FDA00002074562900011.GIF" wi="1645" he="64" /><img file="FDA00002074562900012.GIF" wi="1644" he="63" /><img file="FDA00002074562900013.GIF" wi="1645" he="64" /><img file="FDA00002074562900014.GIF" wi="1646" he="64" /><img file="FDA00002074562900015.GIF" wi="1646" he="58" /><img file="FDA00002074562900016.GIF" wi="1645" he="64" /><img file="FDA00002074562900017.GIF" wi="1646" he="93" /><img file="FDA00002074562900018.GIF" wi="760" he="47" /><img file="FDA00002074562900019.GIF" wi="1540" he="100" /><img file="FDA000020745629000110.GIF" wi="1581" he="100" /><img file="FDA000020745629000111.GIF" wi="965" he="61" /><img file="FDA000020745629000112.GIF" wi="1348" he="99" /><img file="FDA000020745629000113.GIF" wi="1575" he="100" /><img file="FDA000020745629000114.GIF" wi="1247" he="61" /><img file="FDA000020745629000115.GIF" wi="1496" he="100" /><img file="FDA000020745629000116.GIF" wi="1593" he="100" /><img file="FDA000020745629000117.GIF" wi="996" he="60" /><img file="FDA000020745629000118.GIF" wi="1536" he="100" /><img file="FDA000020745629000119.GIF" wi="1600" he="100" /><img file="FDA000020745629000120.GIF" wi="312" he="61" /><img file="FDA00002074562900021.GIF" wi="1466" he="100" /><img file="FDA00002074562900022.GIF" wi="1263" he="51" /><img file="FDA00002074562900023.GIF" wi="1570" he="82" /><img file="FDA00002074562900024.GIF" wi="1362" he="99" /><img file="FDA00002074562900025.GIF" wi="1025" he="52" /><img file="FDA00002074562900026.GIF" wi="1454" he="61" /><img file="FDA00002074562900027.GIF" wi="1617" he="119" /><img file="FDA00002074562900028.GIF" wi="1615" he="114" /><img file="FDA00002074562900029.GIF" wi="1616" he="120" />其中,x<sub>g</sub>、y<sub>g</sub>和z<sub>g</sub>分别为飞行器在地面坐标系中三个方向的位置,V<sub>a</sub>为飞行器对空速度,γ<sub>a</sub>为航迹倾斜角,χ<sub>a</sub>为航迹方位角,α为迎角,β为侧滑角,μ<sub>a</sub>为航迹滚转角,p<sub>k</sub>、q<sub>k</sub>和r<sub>k</sub>为飞行器相对地面坐标系的转动角速度矢量在机体坐标系上的分量,T<sub>xb</sub>、T<sub>yb</sub>和T<sub>zb</sub>分别为推力T在机体坐标系三个方向上的分量,D为阻力,Y为侧力,L为升力,l<sub>A0</sub>、m<sub>A0</sub>和n<sub>A0</sub>为气动力矩,分别为滚转力矩、俯仰力矩和偏航力矩,l<sub>C</sub>、m<sub>C</sub>和n<sub>C</sub>为操纵面动作所产生的控制力矩,Δl<sub>A</sub>、Δm<sub>A</sub>和Δn<sub>A</sub>是由风场引起的不确定力矩,I<sub>x</sub>、I<sub>y</sub>和I<sub>z</sub>为飞行器的惯性矩,M为飞行器总质量,g为重力加速度,u<sub>wg</sub>、v<sub>wg</sub>和w<sub>wg</sub>是风速矢量V<sub>w</sub>在地面坐标系上的分量;(2)建立与所述近空间高超声速飞行器匹配的发动机推力模型,单台发动机估算模型为<img file="FDA000020745629000210.GIF" wi="1619" he="58" /><img file="FDA000020745629000211.GIF" wi="1618" he="108" />其中,<img file="FDA000020745629000212.GIF" wi="42" he="84" />为动压,C<sub>T</sub>为推力系数,<img file="FDA000020745629000213.GIF" wi="50" he="46" />为发动机中燃料流速,推力T和比冲量I<sub>sp</sub>是与马赫数Ma、燃料当量比η相关的函数;此外,发动机系统采用具有推力矢量控制的变推力组合发动机系统,采用“X”型发动机配置,即同时安装4台相同的摆动发动机,它们可分别沿切线方向摆动,定义发动机推力矢量等效偏转角δ<sub>x</sub>、δ<sub>y</sub>、δ<sub>z</sub>,沿机体坐标系三轴方向的推力分量和推力力矩可简化为T<sub>xb</sub>=4F<sub>c</sub>=T             (18)<img file="FDA000020745629000214.GIF" wi="1581" he="101" /><img file="FDA000020745629000215.GIF" wi="1578" he="103" /><img file="FDA00002074562900031.GIF" wi="1591" he="81" /><img file="FDA00002074562900032.GIF" wi="1587" he="101" /><img file="FDA00002074562900033.GIF" wi="1583" he="103" />其中,F<sub>c</sub>为单台发动机推力,X<sub>rc</sub>为发动机推力作用线与机体纵轴距离,X<sub>m</sub>为质心到前缘的距离,X<sub>T</sub>为主发动机推力中心到前缘距离;(3)建立与运动模型和推力模型相关的气动参数库和发动机参数库;(4)根据不同的任务要求,基于运动模型设计姿态稳定跟踪控制律及重心运动控制律;(5)实时计算当前状态下的大气数据和风场数据,并据此以及所设计的控制律计算飞行器的各个状态,重复此步骤直至仿真结束。 
地址 210016 江苏省南京市白下区御道街29号