CN102073755A - 近空间高超声速飞行器运动控制仿真方法 - Google Patents

近空间高超声速飞行器运动控制仿真方法 Download PDF

Info

Publication number
CN102073755A
CN102073755A CN 201010537312 CN201010537312A CN102073755A CN 102073755 A CN102073755 A CN 102073755A CN 201010537312 CN201010537312 CN 201010537312 CN 201010537312 A CN201010537312 A CN 201010537312A CN 102073755 A CN102073755 A CN 102073755A
Authority
CN
China
Prior art keywords
aircraft
mergeformat
thrust
mtplaceref
mteqn
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
Application number
CN 201010537312
Other languages
English (en)
Other versions
CN102073755B (zh
Inventor
程路
都延丽
姜长生
薛雅丽
蒲明
傅健
文杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN 201010537312 priority Critical patent/CN102073755B/zh
Publication of CN102073755A publication Critical patent/CN102073755A/zh
Application granted granted Critical
Publication of CN102073755B publication Critical patent/CN102073755B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明公开一种近空间高超声速飞行器的运动控制仿真方法，包括如下步骤：(1)建立飞行器在变化风场条件下的运动模型；(2)建立所述近空间高超声速飞行器的发动机推力模型；(3)建立与运动模型和推力模型相关的气动参数库和发动机参数库；(4)根据不同的任务要求，基于运动模型设计姿态稳定跟踪控制律及重心运动控制律；(5)实时计算当前状态下的大气数据和风场数据，并据此以及所设计的控制律计算飞行器的各个状态，重复此步骤直至仿真结束。此种仿真方法在充分剖析近空间高超声速飞行器的动力学和运动学规律的基础上，结合发动机推力模型和大气风场环境的实际情况，提出一整套在变化风场环境下的近空间高超声速飞行器的运动模型。本发明还公开一种基于该飞行器的控制仿真平台。

Description

近空间高超声速飞行器运动控制仿真方法
技术领域
 本发明涉及一种航空航天领域的建模和控制仿真问题，具体地说，是飞行器的运动建模与控制仿真方法，可应用于近空间高超声速飞行器飞行稳定性分析和自动控制方案的设计中。
背景技术
 由于近空间环境和高超声速飞行的特殊性，近空间高超声速飞行器的稳定智能自主自适应控制面临着一些重要的技术难题，如近空间环境的精确数学描述、高超声速飞行器的气动布局设计和运动建模、气动参数数据库的实验和计算获取、所采用发动机的推力模型和推力方案设计以及与研究密切相关的飞控系统仿真平台的设计等。
目前，近空间高超声速飞行器主要采用以下四种构型的运动模型：乘波体、升力体、翼身融合体和轴对称体。乘波体构型在高超声速范围内具有较高的升阻比，Rasmussen等人(Rasmussen M L. Waverider configurations derived from inclined circular and elliptic cones. Journal of Spacecraft and Rockets, 1980,17(6): 537-545)提出了通过圆锥流场生成乘波体的方法，得到了广泛的理论和实验研究。升力体作为一种航天飞机和可重复使用运载器普遍使用的构型主要依靠机身形成的升力面而稳定飞行，可消除机身部件所产生的附加阻力和翼身之间的干扰。翼身融合体的气动布局是将机翼和机身两个部件相互融合没有明显界限，具有结构重量轻、气动阻力小等特点，如B2轰炸机。轴对称体构型具有尖头部、后掠角大、展弦比小等特点，适合高超声速飞行器一体化设计。
美国国家宇航局兰利研究中心针对高超声速飞行器提出了一种带翼锥形体构型的气动模型。该构型的前体是轴对称圆锥体，喷管部分是锥平截头体, 发动机模块环绕整个机体。其贡献在于，采用计算流体力学软件和部分风洞试验验证相结合的手段，提供了一套比较完整的气动数据，为飞行控制研究提供了参考对象。随后的研究者利用数值计算的方法拟合了气动数据，给出了气动参数同飞行器迎角和马赫数的非线性表达式，也提供了包括气动力、气动力矩和刚体转动惯量的数学表达。而由公开文献所代表的此类研究工作并未考虑飞行器在复杂的近空间环境中和特定的发动机推力模型下的建模问题。
因此，需要一种可以综合考虑近空间风场环境、发动机推力模型的近空间高超声速飞行器运动建模方法，并以此为基础设计一种较为通用的仿真平台结构。
发明内容
 本发明所要解决的技术问题是提供一种近空间高超声速飞行器的运动建模方法及基于该飞行器的控制仿真平台，其在充分剖析近空间高超声速飞行器的动力学和运动学规律的基础上，结合发动机推力模型和大气风场环境的实际情况，提出近空间高超声速飞行器在不同坐标系的运动方程，并以此为基础，研究飞行器控制仿真的结构。
本发明为解决上述技术问题，采用的技术方案是：
一种近空间高超声速飞行器的控制仿真方法，包括如下步骤：
(1)建立飞行器在变化风场条件下的运动模型；
Figure 872647DEST_PATH_IMAGE001
                                       (1)
Figure 727470DEST_PATH_IMAGE002
                                          (2)
Figure 444890DEST_PATH_IMAGE003
                                          (3)
Figure 777783DEST_PATH_IMAGE004
                                      (4)
Figure 529838DEST_PATH_IMAGE005
                                      (5)
Figure 352301DEST_PATH_IMAGE006
                                    (6)
Figure 554087DEST_PATH_IMAGE007
               (7)
Figure 690670DEST_PATH_IMAGE008
(8)
Figure 359549DEST_PATH_IMAGE009
 (9)
 (10)
Figure 982609DEST_PATH_IMAGE011
(11)
Figure 657304DEST_PATH_IMAGE012
(12)
Figure 180689DEST_PATH_IMAGE013
                                 (13)
Figure 548216DEST_PATH_IMAGE014
                               (14)
Figure 727525DEST_PATH_IMAGE015
                               (15)
其中，
Figure 940331DEST_PATH_IMAGE016
、
Figure 518556DEST_PATH_IMAGE017
和
Figure 791405DEST_PATH_IMAGE018
分别为飞行器在地面坐标系中三个方向的位置，
Figure 458010DEST_PATH_IMAGE019
为飞行器对空速度，
Figure 474508DEST_PATH_IMAGE020
为航迹倾斜角，
Figure 972485DEST_PATH_IMAGE021
为航迹方位角，
Figure 681815DEST_PATH_IMAGE022
为迎角，
Figure DEST_PATH_IMAGE023
为侧滑角，为航迹滚转角，
Figure DEST_PATH_IMAGE025
、
Figure 328008DEST_PATH_IMAGE026
和为飞行器相对惯性系(地面坐标系)的转动角速度矢量在机体坐标系上的分量，、和
Figure 433160DEST_PATH_IMAGE030
分别为推力
Figure DEST_PATH_IMAGE031
在机体坐标系三个方向上的分量，
Figure 808778DEST_PATH_IMAGE032
为阻力，为侧力，
Figure 104761DEST_PATH_IMAGE034
为升力，
Figure DEST_PATH_IMAGE035
、
Figure 515014DEST_PATH_IMAGE036
和
Figure DEST_PATH_IMAGE037
为气动力矩，分别为滚转力矩、俯仰力矩和偏航力矩，
Figure 500900DEST_PATH_IMAGE038
、
Figure DEST_PATH_IMAGE039
和
Figure 363813DEST_PATH_IMAGE040
为操纵面动作所产生的控制力矩，
Figure DEST_PATH_IMAGE041
、
Figure 463488DEST_PATH_IMAGE042
和
Figure DEST_PATH_IMAGE043
是由风场引起的不确定力矩，
Figure 462668DEST_PATH_IMAGE044
、
Figure DEST_PATH_IMAGE045
和
Figure 684702DEST_PATH_IMAGE046
为飞行器的惯性矩，为飞行器总质量，
Figure 972595DEST_PATH_IMAGE048
为重力加速度，
Figure DEST_PATH_IMAGE049
、
Figure 935347DEST_PATH_IMAGE050
和
Figure DEST_PATH_IMAGE051
风速矢量
Figure 992296DEST_PATH_IMAGE052
在地面坐标系上的分量；
(2)建立与所述近空间高超声速飞行器匹配的发动机推力模型，单台发动机估算模型为
Figure DEST_PATH_IMAGE053
                                          (16)
Figure 119652DEST_PATH_IMAGE054
                                              (17)
其中，
Figure DEST_PATH_IMAGE055
为动压，
Figure 894841DEST_PATH_IMAGE056
为推力系数，
Figure DEST_PATH_IMAGE057
为发动机中燃料流速，推力
Figure 398634DEST_PATH_IMAGE031
和比冲量
Figure 572739DEST_PATH_IMAGE058
是与马赫数、燃料当量比
Figure 870997DEST_PATH_IMAGE060
相关的函数；此外，发动机系统采用具有推力矢量控制的变推力组合发动机系统，采用“X”型发动机配置，即同时安装4台相同的摆动发动机，它们可分别沿切线方向摆动，定义发动机推力矢量等效偏转角
Figure DEST_PATH_IMAGE061
、
Figure 133482DEST_PATH_IMAGE062
、，沿机体坐标系三轴方向的推力分量和推力力矩可简化为
                                           (18)
Figure DEST_PATH_IMAGE065
                                          (19)
Figure 206928DEST_PATH_IMAGE066
                                          (20)
Figure DEST_PATH_IMAGE067
                                        (21)
Figure 941666DEST_PATH_IMAGE068
                                   (22)
Figure DEST_PATH_IMAGE069
                                  (23)
其中，
Figure 485255DEST_PATH_IMAGE070
为单台发动机推力，
Figure DEST_PATH_IMAGE071
为发动机推力作用线与机体纵轴距离，为质心到前缘的距离，
Figure DEST_PATH_IMAGE073
为主发动机推力中心到前缘距离；
(3)建立与运动模型和推力模型相关的气动参数库和发动机参数库；
(4)根据不同的任务要求，基于运动模型设计姿态稳定跟踪控制律及重心运动控制律；
(5)实时计算当前状态下的大气数据和风场数据，并据此以及所设计的控制律计算飞行器的各个状态，重复此步骤直至仿真结束。
采用上述方案后，本发明与现有技术相比的优点在于：
1)通过对变化风场扰动下近空间高超声速飞行器飞行运动的机理分析，建立了在风场中的飞行器运动全维状态方程模型，使得飞行控制的研究更具有针对性且更切合物理实际；
2)本发明提出了一种用于近空间高超声速飞行器的吸气式发动机推力估算模型和推力矢量配置方案，解决了对发动机推力进行闭环控制的问题，有利于进行机体/发动机一体化控制方案的设计；
3)本发明提出了面向控制的模型处理方式，给出了控制器设计的基本思路，便于展开针对在不确定条件下近空间高超声速飞行器的非线性鲁棒自适应控制方案的设计；
4)本发明提出了一种用于在变化风场和不确定环境下进行近空间高超声速飞行器控制仿真的软件平台，使得控制方案的仿真验证更加方便和明了。
 附图说明
 图1(a)是近空间高超声速飞行器气动构型的俯视图；
图1(b)是近空间高超声速飞行器气动构型的侧视图；
图2是近空间高超声速飞行器相关坐标系之间关系的纵向剖视图；
图3是发动机推力系数在给定试验点数据的网格图；
图4是发动机推力系数的全局拟合曲面；
图5是发动机比冲量在给定试验点数据的网格图；
图6是发动机比冲量的全局拟合曲面；
图7是推力矢量发动机的配置图；
图8是近空间高超声速飞行器飞控系统结构图；
图9是飞行控制仿真平台流程图；
图10是基于MATLAB\SIMULINK飞控仿真平台的软件界面；
图11是大气风场生成器的仿真结构。
具体实施方式
 本发明提供一种近空间高超声速飞行器的控制仿真方法，包括如下步骤：
(1)建立飞行器在变化风场条件下的运动模型；
以下将详细说明近空间高超声速飞行器的建模过程。
本发明提供的近空间高超声速飞行器是基于美国国家宇航局兰利研究中心的技术报告(Shaughnessy J D, Pinckney S Z, McMinn J D, Cruz C I, Kelley M-L. Hypersonic Vehicle Simulation Model: Winged-cone Configuration. NASA TM-102610, 1990)提出的有翼锥形体结构，具有三角形机翼，具有单垂尾方向舵、可独立工作的左右升降副翼舵，其气动布局如图1(a)和图1(b)所示，该飞行器的相关几何参数可参见此技术报告，在此不再赘述。
由于近空间飞行环境相对复杂，这给近空间高超声速飞行器的控制系统设计带来了一定的难度，因此在复杂大气环境下的建模就显得至关重要。大气环境的描述一般包括气温、空气密度以及风速等与高度的关系，本发明主要研究在变化风场下的飞行器建模问题。
为了后续说明方便，首先由图2给出描述飞行器运动方程的不同坐标系在纵向剖面上相互之间的关系，其中包括地面坐标系
Figure 420161DEST_PATH_IMAGE074
、机体坐标系、气流坐标系
Figure 325801DEST_PATH_IMAGE076
以及航迹坐标系
Figure DEST_PATH_IMAGE077
。
Figure 297299DEST_PATH_IMAGE078
、
Figure DEST_PATH_IMAGE079
、
Figure 779594DEST_PATH_IMAGE080
及
Figure DEST_PATH_IMAGE081
分别为相应的两个坐标系之间的转换矩阵，各自对应的逆转换矩阵分别为
Figure 520148DEST_PATH_IMAGE081
、、
Figure DEST_PATH_IMAGE083
及
Figure 117799DEST_PATH_IMAGE084
，此为现有知识，详细表达式可见相关参考资料(如肖业伦、金长江所著“大气扰动中的飞行原理”，国防工业出版社，1993)。
一般情况下，在飞行过程中，飞行器的对地速度矢量、对空速度矢量
Figure 242881DEST_PATH_IMAGE086
与即时风速矢量
Figure DEST_PATH_IMAGE087
构成矢量三角形关系
                       
Figure 165838DEST_PATH_IMAGE088
                      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (1)
其中各矢量的元素分别为在机体坐标系三个坐标轴上的分量。考虑到目前对风速的测量一般采用在地面坐标系上的分量
Figure DEST_PATH_IMAGE089
，因此，下文的推导中在必要处作了相应的坐标转换。
根据牛顿第二定律，飞行器的质心运动方程可写作
                       
Figure 348033DEST_PATH_IMAGE090
                      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (2)
其中，
Figure 356440DEST_PATH_IMAGE047
为飞行器总质量，
Figure DEST_PATH_IMAGE091
为飞行器所受的合外力矢量。由此可以得到飞行器在不同坐标系下的质心运动方程。
首先，不考虑地球曲率和旋转，将地面坐标系作为惯性系，可得
       
Figure 347530DEST_PATH_IMAGE092
     MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (3)
其中，
Figure DEST_PATH_IMAGE093
为惯性系中飞行器的位置矢量，将上式化简后可得飞行器的导航方程组
                                   MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (4)
其中，
Figure DEST_PATH_IMAGE095
为飞行器对空速度矢量
Figure 481019DEST_PATH_IMAGE096
的模，
Figure DEST_PATH_IMAGE097
为航迹倾斜角，
Figure 711144DEST_PATH_IMAGE098
为航迹方位角。
其次，由机体坐标系下的质心运动方程可得
            
Figure DEST_PATH_IMAGE099
          MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (5)
其中，定义为飞行器相对惯性系(地面坐标系)的转动角速度矢量在机体坐标系上的分量，
Figure DEST_PATH_IMAGE101
、
Figure 72647DEST_PATH_IMAGE102
和分别为推力在机体坐标系三个方向上的分量，
Figure 644891DEST_PATH_IMAGE104
为阻力，
Figure 977783DEST_PATH_IMAGE033
为侧力，
Figure 464260DEST_PATH_IMAGE034
为升力。将式(1)代入式(5)整理后，可得到飞行器的力方程组
        
Figure DEST_PATH_IMAGE105
      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (6)
其中，
Figure 489984DEST_PATH_IMAGE022
为迎角，
Figure 691771DEST_PATH_IMAGE023
为侧滑角。
再者，由航迹坐标系下的质心运动方程可得
          
Figure 828354DEST_PATH_IMAGE106
         MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (7)
其中，
Figure DEST_PATH_IMAGE107
为飞行器对地速度矢量
Figure 169337DEST_PATH_IMAGE108
的模。同样将式(1)代入式(7)，并考虑到
Figure DEST_PATH_IMAGE109
，可整理得到另一组力方程组
          
Figure 303646DEST_PATH_IMAGE110
         MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (8)
此外，根据关系式
Figure DEST_PATH_IMAGE111
可计算出航迹滚转角的状态方程。
在上述公式(4)、(6)、(8)中，其中的最后一项均是关于风切变的描述，具体表达如下
                       
Figure 995659DEST_PATH_IMAGE112
                      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (9)
其中，
Figure DEST_PATH_IMAGE113
为地面坐标系下风场矢量对时间的导数；
                                      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (10)
为风场的梯度矩阵，其中
Figure DEST_PATH_IMAGE115
，
Figure 66175DEST_PATH_IMAGE116
，，其它元素的表达类同；此外，
Figure 409749DEST_PATH_IMAGE118
为地面坐标系下飞行器对地速度矢量。显然地
                    MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (11)
将式(11)代入式(9)，即可得到风场矢量的导数。
进一步，作为刚体的飞行器的转动矢量方程为
                        
Figure 469289DEST_PATH_IMAGE120
                      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (12)
其中，
Figure 742139DEST_PATH_IMAGE121
为飞行器对质心的动量矩，取决于惯性矩和角速度，
Figure 471060DEST_PATH_IMAGE122
是作用在飞行器上的所有外力对质心的总力矩矢量。因此，在机体坐标系下的转动矢量方程
                   
Figure 487558DEST_PATH_IMAGE123
                 MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (13)
式中， 
Figure 923218DEST_PATH_IMAGE124
、
Figure 429286DEST_PATH_IMAGE125
和
Figure 583187DEST_PATH_IMAGE126
是
Figure 400445DEST_PATH_IMAGE122
在机体坐标系下的分量，分别称滚转力矩、俯仰力矩和偏航力矩，且有
                      
Figure 690612DEST_PATH_IMAGE127
                    MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (14)
其中，为气动力矩矢量，
Figure 8778DEST_PATH_IMAGE129
为操作面动作产生的控制力矩矢量；此外，
Figure 367079DEST_PATH_IMAGE130
、和
Figure 625202DEST_PATH_IMAGE132
为在机体坐标系下的分量，并有
                  
Figure 650106DEST_PATH_IMAGE133
                MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (15)
Figure 914866DEST_PATH_IMAGE134
、
Figure 133970DEST_PATH_IMAGE135
和为飞行器的惯性矩，
Figure 246599DEST_PATH_IMAGE137
、
Figure 365865DEST_PATH_IMAGE138
和
Figure 758800DEST_PATH_IMAGE139
为惯性积。考虑到飞行器为对称体，则有
Figure 658623DEST_PATH_IMAGE140
，将上式代入式(13)整理化简后可得力矩方程组
             
Figure 162417DEST_PATH_IMAGE141
           MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (16)
这里，气动力矩分量
Figure 136189DEST_PATH_IMAGE142
、
Figure 762343DEST_PATH_IMAGE143
和
Figure 87145DEST_PATH_IMAGE144
是关于、
Figure 219978DEST_PATH_IMAGE146
和
Figure 17033DEST_PATH_IMAGE147
的函数，具体表达式可见下文，而
Figure 563552DEST_PATH_IMAGE148
是飞行器相对于周围空气场的转动角速度矢量，现有传感技术无法对其精确测量，因此式(16)中的气动力矩的数值无法获取，所以本发明做如下处理。根据式(1)所表达的三角关系，有如下关系成立
                       
Figure 674727DEST_PATH_IMAGE149
                     MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (17)
其中，
Figure 623092DEST_PATH_IMAGE150
为空气介质相对于惯性场的转动角速度矢量，它是由风梯度造成的。将此式展开，有
                
Figure 591048DEST_PATH_IMAGE151
              MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (18)
其中最后一项中的各变量为风梯度，具体定义可见式(10)。这里，
Figure 624863DEST_PATH_IMAGE100
的数值可由陀螺仪测量得到，为此我们将气动力矩分量
Figure 8571DEST_PATH_IMAGE142
、和的计算分解成两部分：可由
Figure 736990DEST_PATH_IMAGE154
、和
Figure 930170DEST_PATH_IMAGE156
的测量值计算的部分，
Figure 177611DEST_PATH_IMAGE157
、和
Figure 239425DEST_PATH_IMAGE159
；以及由存在风梯度而产生的
Figure 16889DEST_PATH_IMAGE160
、和所引起的不确定项
Figure 522453DEST_PATH_IMAGE163
、
Figure 154423DEST_PATH_IMAGE164
和
Figure DEST_PATH_IMAGE165
。因此，式(16)可改写为
           
Figure 743667DEST_PATH_IMAGE166
         MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (19)
将以上式(4)、(6)、(8)、(9)和(19)整理展开，得到近空间高超声速飞行器在风场下完整的运动学方程组
                                      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (20)
                    
Figure 661420DEST_PATH_IMAGE168
                  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (21)
                     
Figure DEST_PATH_IMAGE169
                   MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (22)
                    
Figure 994313DEST_PATH_IMAGE170
                 MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (23)
                    
Figure 480789DEST_PATH_IMAGE005
                 MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (24)
                   
Figure 506513DEST_PATH_IMAGE006
                 MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (25)
        
Figure 711230DEST_PATH_IMAGE007
      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (26)
Figure 644551DEST_PATH_IMAGE008
 MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (27)
Figure 251113DEST_PATH_IMAGE009
  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (28)
Figure 447739DEST_PATH_IMAGE010
  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (29)
Figure 936489DEST_PATH_IMAGE011
  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (30)
Figure 873833DEST_PATH_IMAGE012
  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (31)
                  
Figure 334902DEST_PATH_IMAGE013
                MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (32)
                 
Figure 436850DEST_PATH_IMAGE014
               MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (33)
                 
Figure 678475DEST_PATH_IMAGE015
               MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (34)
其中，各个气动力的计算如下
                        
Figure DEST_PATH_IMAGE171
                      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (35)
                        
Figure 94544DEST_PATH_IMAGE172
                      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (36)
                        
Figure DEST_PATH_IMAGE173
                      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (37)
                        
Figure 410119DEST_PATH_IMAGE174
                      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (38)
                      
Figure DEST_PATH_IMAGE175
                    MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (39)
                      
Figure 886231DEST_PATH_IMAGE176
                    MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (40)
式中
                       
Figure DEST_PATH_IMAGE177
                     MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (41)
                       
Figure 284327DEST_PATH_IMAGE178
                     MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (42)
                     
Figure DEST_PATH_IMAGE179
                  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (43)
                   
Figure 504087DEST_PATH_IMAGE180
                 MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (44)
                              MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (45)
                
Figure 939747DEST_PATH_IMAGE182
              MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (46)
           
Figure DEST_PATH_IMAGE183
         MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (47)
              
Figure 321181DEST_PATH_IMAGE184
           MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (48)
           
Figure 537399DEST_PATH_IMAGE185
         MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (49)
其中，
Figure 357587DEST_PATH_IMAGE186
为机翼参考面积， 
Figure 457874DEST_PATH_IMAGE187
为平均气动翼弦，
Figure 72526DEST_PATH_IMAGE188
为翼展，
Figure 510461DEST_PATH_IMAGE189
为质心到参考力矩中心的距离，
Figure 134340DEST_PATH_IMAGE190
为动压，为空气密度，
Figure 330146DEST_PATH_IMAGE192
为相应的气动导数，它们均是迎角和马赫数的函数；
Figure 255377DEST_PATH_IMAGE193
和
Figure 417368DEST_PATH_IMAGE194
分别为左、右升降副翼舵偏角，
Figure 679197DEST_PATH_IMAGE195
为方向舵偏角。
(2)建立所述近空间高超声速飞行器的发动机推力模型；
以下，本发明将根据公开文献所提供的数据，分析设计一种适合近空间高超声速飞行器的发动机推力工程估算模型。本发明提供了一种通过动压来估算发动机推力及比冲量的关系式
                      
Figure 901231DEST_PATH_IMAGE196
                    MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (50)
                        
Figure 48179DEST_PATH_IMAGE197
                      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (51)
                       
Figure 13861DEST_PATH_IMAGE198
                     MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (52)
其中，
Figure 133127DEST_PATH_IMAGE031
为发动机推力，单位
Figure 526062DEST_PATH_IMAGE199
；
Figure 425885DEST_PATH_IMAGE200
为动压，单位
Figure 929678DEST_PATH_IMAGE201
；
Figure 903451DEST_PATH_IMAGE202
为推力系数，单位
Figure 467287DEST_PATH_IMAGE203
；
Figure 588827DEST_PATH_IMAGE059
为马赫数，无量纲；
Figure 893382DEST_PATH_IMAGE204
为燃料当量比，无量纲；
Figure 987240DEST_PATH_IMAGE205
为比冲量，无量纲；
Figure DEST_PATH_IMAGE206
、
Figure 721978DEST_PATH_IMAGE207
分别为发动机中燃料和空气流速，单位。
由上面的模型可知，推力系数和比冲量均是与马赫数及燃料当量比相关的函数，根据相关参考资料(Irene M. Gregory, Rajiv S. Chowdhry, John D. McMinn, John D. Shaughnessy. Hypersonic Vehicle Model and Control Law Development Using 
Figure 206180DEST_PATH_IMAGE209
 and 
Figure DEST_PATH_IMAGE210
 Synthesis. NASA TM-4562. 1994)提供的实验和理论数据，本发明对某些给定的参考点的推力系数和比冲量进行了预测计算，具体数值由表1和表2给出，其中，表1表示推力系数与马赫数、燃料当量比的关系，表2表示比冲量与马赫数、燃料当量比的关系。
 
表1
Figure 51776DEST_PATH_IMAGE211
表2
Figure DEST_PATH_IMAGE212
在给定的马赫数情况下，推力系数与燃料当量比近似为线性关系，而当燃料当量比小于0.3之后，推力系数将随其减小而急剧下降，这正体现了发动机的熄火过程。
为了全方位估算在高超声速阶段不同马赫数和不同燃料当量比的情况下的推力系数和发动机比冲量，本发明根据以上数据，采用双立方插值算法进行二元拟合，从而实现高超声速阶段发动机推力参数的全局估算。双立方插值算法作为一种有效的线性插值算法，具有较高的预测精度。图3至图6以三维坐标图的形式分别给出了推力系数数据、推力系数全局拟合结果、比冲量数据以及比冲量全局拟合结果。从图中也可看出双立方插值算法较好的逼近效果。
发动机作为近空间高超声速飞行器模型的主要组成部分之一，在飞行控制中起着不可或缺的作用。在控制仿真模型中，发动机系统采用具有推力矢量控制的变推力组合发动机系统，采用“X”型发动机配置，即同时安装4台相同的摆动发动机，它们可分别沿切线方向摆动, 其编号以及相应的摆动角
Figure 937823DEST_PATH_IMAGE213
、
Figure DEST_PATH_IMAGE214
、
Figure 840533DEST_PATH_IMAGE215
及
Figure DEST_PATH_IMAGE216
的正方向定义如图7所示。不难计算出4台摆动发动机总推力沿机体轴三个方向的分量大小：
                
Figure 874348DEST_PATH_IMAGE217
              MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (53)
                            MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (54)
               
Figure 461318DEST_PATH_IMAGE219
             MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (55)
相应的推力力矩分别为
               
Figure DEST_PATH_IMAGE220
             MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (56)
                     
Figure 529768DEST_PATH_IMAGE221
                   MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (57)
                                        MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (58)
其中，
Figure 543992DEST_PATH_IMAGE070
为单台发动机推力，
Figure 799524DEST_PATH_IMAGE223
为质心到前缘的距离，
Figure DEST_PATH_IMAGE224
为发动机推力作用线与机体纵轴距离，
Figure 273010DEST_PATH_IMAGE225
为主发动机推力中心到前缘距离，为参考力矩中心到前缘距离。
为了便于控制系统的设计，本发明考虑采用等效偏转角来代替四台发动机的摆动角。考虑到摆动发动机的摆动角均限制在
Figure 195966DEST_PATH_IMAGE227
的范围内，因此可近似表达为
                 
Figure DEST_PATH_IMAGE228
               MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (59)
定义发动机推力矢量等效偏转角为
                    
Figure 381091DEST_PATH_IMAGE229
                  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (60)
                    
Figure DEST_PATH_IMAGE230
                  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (61)
                                      MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (62)
那么，沿机体坐标系三轴方向的推力分量和推力力矩可简化为
                       
Figure 115009DEST_PATH_IMAGE064
                     MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (63)
                                           MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (64)
                       
Figure 310815DEST_PATH_IMAGE066
                     MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (65)
                      
Figure 868836DEST_PATH_IMAGE067
                   MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (66)
                   
Figure 395107DEST_PATH_IMAGE068
                 MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (67)
                   
Figure 292656DEST_PATH_IMAGE069
                 MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (68)
据此，我们可以通过对发动机燃料当量比和推力矢量等效偏转角的控制来实现发动机推力大小和方向的变化。
至此，包含变化风场及发动机推力模型在内的近空间高超声速飞行器的运动建模已基本完成。
(3)建立与运动模型和推力模型相关的气动参数库和发动机参数库；
以下，根据气动力和气动力矩的计算方式(44)-(49)以及单台发动机的模型数据(见表1和表2)，参考NASA技术报告(Shaughnessy J D, Pinckney S Z, McMinn J D, Cruz C I, Kelley M-L. Hypersonic Vehicle Simulation Model: Winged-cone Configuration. NASA TM-102610, 1990)所提供的数据，基于计算流体力学和插值拟合算法建立近空间高超声速飞行器气动参数库和发动机参数库。
(4)根据不同的任务要求，基于运动模型设计姿态稳定跟踪控制律及重心运动控制律；
以下，对控制系统的结构和设计思路做简要说明。
飞行控制系统的任务主要包括：飞行轨迹跟踪和飞行姿态稳定两个部分。飞机的操纵机构包括：左右升降副翼舵和方向舵的舵偏角指令
Figure DEST_PATH_IMAGE232
、
Figure 881900DEST_PATH_IMAGE233
、
Figure DEST_PATH_IMAGE234
；发动机燃料当量比；推力矢量等效偏转角
Figure 932213DEST_PATH_IMAGE235
、
Figure DEST_PATH_IMAGE236
、
Figure 621952DEST_PATH_IMAGE237
。根据不同的飞行任务将采用不同的控制器结构和操纵机构。近空间高超声速飞行器飞行控制系统的一般结构如图8所示。
一个合理可行的飞控系统的方案设计必然是在保证姿态稳定的前提下根据不同的飞行任务分别设计针对性的控制器，以达到相应的任务要求。此外，由于近空间飞行环境复杂且高超声速飞行动态变化快，控制系统难免存在外干扰、不确定和建模误差(这里统称复合干扰)，因此有必要对复合干扰作观测估计并加以补偿控制，以提高系统的鲁棒性。控制器的具体设计思路说明如下：
其一、在姿态稳定控制中，将发动机控制作为开环控制量，仅考虑其推力矢量控制，飞控系统的设计就是设计控制力矩
Figure DEST_PATH_IMAGE238
，并根据一定的算法分配成舵偏角指令
Figure 382097DEST_PATH_IMAGE239
、
Figure DEST_PATH_IMAGE240
、和推力矢量等效偏转角
Figure DEST_PATH_IMAGE242
、、
Figure DEST_PATH_IMAGE244
，最终使飞行器的姿态稳定跟踪给定的指令。本发明采用的控制分配关系为
                
Figure 202395DEST_PATH_IMAGE245
              MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (69)
其中，分配矩阵
Figure DEST_PATH_IMAGE246
的具体表达式可见文献(朱亮，空天飞行器不确定非线性鲁棒自适应控制，南京航空航天大学博士论文，2006)。
为此我们将与之相关的运动方程写成仿射非线性形式，并将其分成慢、快两个回路，分别进行控制律设计。为了实现对气流角及航迹角状态
Figure 71125DEST_PATH_IMAGE247
的稳定跟踪，姿态角慢回路(29)～(31)和角速率快回路(32)～(34)的控制方程可分别写为
                   
Figure DEST_PATH_IMAGE248
                 MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (70)
                  
Figure 697891DEST_PATH_IMAGE249
                MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT (71)
式中，
Figure DEST_PATH_IMAGE250
和
Figure 638165DEST_PATH_IMAGE251
分别为快、慢回路的控制量，
Figure DEST_PATH_IMAGE252
、
Figure 99234DEST_PATH_IMAGE253
分别为各回路的复合干扰，各变量的具体表达式在此不做展开。研究工作者可以通过这一组仿射非线性方程，采用不同的控制方法进行控制律的设计。
其二、在重心运动控制中，是在姿态稳定的基础上，考虑到发动机推力及推力矢量控制的作用，飞控系统的设计就是要设计包括发动机推力、推力矢量等效偏转角和气动舵面偏转角在内的控制律，使得飞行器在保证姿态平滑变化的飞行品质下实现重心运动的稳定跟踪。
控制器采用内外环递阶控制的结构，外环是关于导航控制回路，基于非仿射非线性方程(26)、(28)设计相应的气动力和推力
Figure 380490DEST_PATH_IMAGE031
控制以满足导航任务的要求，并通过一定的分配算法将其分配成发动机燃料当量比
Figure 858876DEST_PATH_IMAGE060
和相应的姿态角
Figure 236768DEST_PATH_IMAGE247
。然后再根据前面所提的姿态稳定控制方法，最终设计出推力矢量等效偏转角
Figure 775197DEST_PATH_IMAGE235
、
Figure DEST_PATH_IMAGE254
、
Figure 110976DEST_PATH_IMAGE255
和气动舵面偏角
Figure 455369DEST_PATH_IMAGE232
、
Figure 625451DEST_PATH_IMAGE233
、
Figure 334781DEST_PATH_IMAGE234
等控制量。
(5)实时计算当前状态下的大气数据和风场数据，并据此以及所设计的控制律计算飞行器的各个状态，重复此步骤直至仿真结束。
以下，本发明在运动建模和控制系统设计的前提下，通过所设计的仿真平台对相关飞行环境数据和飞行状态进行实时计算和输出。
为了方便飞行控制系统的设计以及对其中各种算法的测试和验证，有必要研究一套基于近空间高超声速飞行器运动模型的数字仿真平台，利用其完成各种飞行控制方案的设计、验证和评估。本发明基于MATLAB等仿真软件，设计了针对近空间高超声速飞行器飞行控制的仿真平台。该软件平台的流程如图9所示，图10给出了基于MATLAB\SIMULINK飞控仿真平台的软件界面。软件结构主要由支持层和执行层组成。支持层主要为仿真提供相应的数据支持和计算，主要包括大气环境数据库、大气风场生成器、控制方案设计及控制参数调节模块、气动参数库、发动机参数库、飞行器运动模型等。执行层主要根据相应的算法安排，进行一系列必要的仿真计算，并实时输出飞行器的飞行状态。
仿真的一般流程如下：
(1)初始化相关参数，调入飞行器初始状态信息和指令信息；
(2)确定本机信息；
(3)根据飞行器实时信息和大气环境数据库，调节环境参数；
(4)进行当前状态下的大气数据计算；
(5)根据风场生成器所提供信息，进行风速、风梯度等风场数据计算；
(6)综合各种参数和状态，做相应的任务管理；
(7)对预先设计的控制器进行参数调节；
(8)根据控制器输出，分别进行推进计算和气动计算；
(9)根据计算结果，基于飞行器运动方程计算飞行器各个状态；
(10)实时反馈飞行器状态，并判断仿真是否结束，若未结束，则返回步骤(2)；
(11)仿真结束。
本发明近空间高超声速飞行器运动控制仿真平台的具体实施步骤如下：
(1)确定飞行器的气动数据库和发动机参数库完整，建立描述大气环境的数据库，以支持飞行器运动方程解算；
(2)根据任务信息，判断所采用的控制器类型，并根据相应的运动方程设计飞行控制律，并针对复合干扰设计复合干扰观测器、设计补偿控制律，最终确定控制器结构。
(3)通过仿真平台，设计风场模型的仿真模拟，图11给出了大气风场生成器的仿真结构。这里总共需要3个独立的噪音源，即通过白噪声生成器(WNG)来获得大气风场的随机信号，并将其作用于近空间高超声速飞行器运动模型中。图中各方框中的传递函数分别为针对生成三个方向的风速分量及相应的风梯度而简化的成形滤波器传递函数。其参数确定方法可参见文献(肖业伦,金长江.大气扰动中的飞行原理.国防工业出版社,1993)。
(4)将控制器结构、风场生成器等信息导入仿真平台，根据相应的控制仿真任务，设定任务指令和初始状态信息，进行实时的仿真验证。
以上实施方式仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。

Claims (1)

1.一种近空间高超声速飞行器的控制仿真方法，其特征在于包括如下步骤：
(1)建立飞行器在变化风场条件下的运动模型：
Figure DEST_PATH_IMAGE002
                                       (1)
Figure DEST_PATH_IMAGE004
                                           (2)
                                           (3)
Figure DEST_PATH_IMAGE008
                                       (4)
Figure DEST_PATH_IMAGE010
                                       (5)
Figure DEST_PATH_IMAGE012
                                      (6)
Figure DEST_PATH_IMAGE014
               (7)
Figure DEST_PATH_IMAGE016
 (8)
Figure DEST_PATH_IMAGE018
 (9)
Figure DEST_PATH_IMAGE020
 (10)
Figure DEST_PATH_IMAGE022
(11)
Figure DEST_PATH_IMAGE024
(12)
Figure DEST_PATH_IMAGE026
                                 (13)
Figure DEST_PATH_IMAGE028
                               (14)
Figure DEST_PATH_IMAGE030
                                (15)
其中，
Figure DEST_PATH_IMAGE032
、
Figure DEST_PATH_IMAGE034
和
Figure DEST_PATH_IMAGE036
分别为飞行器在地面坐标系中三个方向的位置，
Figure DEST_PATH_IMAGE038
为飞行器对空速度，
Figure DEST_PATH_IMAGE040
为航迹倾斜角，
Figure DEST_PATH_IMAGE042
为航迹方位角，
Figure DEST_PATH_IMAGE044
为迎角，
Figure DEST_PATH_IMAGE046
为侧滑角，
Figure DEST_PATH_IMAGE048
为航迹滚转角，
Figure DEST_PATH_IMAGE050
、
Figure DEST_PATH_IMAGE052
和
Figure DEST_PATH_IMAGE054
为飞行器相对惯性系(地面坐标系)的转动角速度矢量在机体坐标系上的分量，、和
Figure DEST_PATH_IMAGE060
分别为推力
Figure DEST_PATH_IMAGE062
在机体坐标系三个方向上的分量，
Figure DEST_PATH_IMAGE064
为阻力，
Figure DEST_PATH_IMAGE066
为侧力，
Figure DEST_PATH_IMAGE068
为升力，
Figure DEST_PATH_IMAGE070
、
Figure DEST_PATH_IMAGE072
和
Figure DEST_PATH_IMAGE074
为气动力矩，分别为滚转力矩、俯仰力矩和偏航力矩，
Figure DEST_PATH_IMAGE076
、
Figure DEST_PATH_IMAGE078
和
Figure DEST_PATH_IMAGE080
为操纵面动作所产生的控制力矩，
Figure DEST_PATH_IMAGE082
、
Figure DEST_PATH_IMAGE084
和
Figure DEST_PATH_IMAGE086
是由风场引起的不确定力矩，
Figure DEST_PATH_IMAGE088
、
Figure DEST_PATH_IMAGE090
和
Figure DEST_PATH_IMAGE092
为飞行器的惯性矩，
Figure DEST_PATH_IMAGE094
为飞行器总质量，
Figure DEST_PATH_IMAGE096
为重力加速度，
Figure DEST_PATH_IMAGE098
、和
Figure DEST_PATH_IMAGE102
风速矢量
Figure DEST_PATH_IMAGE104
在地面坐标系上的分量；
(2)建立与所述近空间高超声速飞行器匹配的发动机推力模型，单台发动机估算模型为
Figure DEST_PATH_IMAGE106
                                          (16)
Figure DEST_PATH_IMAGE108
                                              (17)
其中，
Figure DEST_PATH_IMAGE110
为动压，
Figure DEST_PATH_IMAGE112
为推力系数，
Figure DEST_PATH_IMAGE114
为发动机中燃料流速，推力
Figure 617916DEST_PATH_IMAGE062
和比冲量
Figure DEST_PATH_IMAGE116
是与马赫数
Figure DEST_PATH_IMAGE118
、燃料当量比
Figure DEST_PATH_IMAGE120
相关的函数；此外，发动机系统采用具有推力矢量控制的变推力组合发动机系统，采用“X”型发动机配置，即同时安装4台相同的摆动发动机，它们可分别沿切线方向摆动，定义发动机推力矢量等效偏转角、
Figure DEST_PATH_IMAGE124
、
Figure DEST_PATH_IMAGE126
，沿机体坐标系三轴方向的推力分量和推力力矩可简化为
Figure DEST_PATH_IMAGE128
                                           (18)
Figure DEST_PATH_IMAGE130
                                          (19)
Figure DEST_PATH_IMAGE132
                                          (20)
Figure DEST_PATH_IMAGE134
                                        (21)
Figure DEST_PATH_IMAGE136
                                   (22)
Figure DEST_PATH_IMAGE138
                                  (23)
其中，
Figure DEST_PATH_IMAGE140
为单台发动机推力，为发动机推力作用线与机体纵轴距离，
Figure DEST_PATH_IMAGE144
为质心到前缘的距离，为主发动机推力中心到前缘距离；
(3)建立与运动模型和推力模型相关的气动参数库和发动机参数库；
(4)根据不同的任务要求，基于运动模型设计姿态稳定跟踪控制律及重心运动控制律；
(5)实时计算当前状态下的大气数据和风场数据，并据此以及所设计的控制律计算飞行器的各个状态，重复此步骤直至仿真结束。
CN 201010537312 2010-11-10 2010-11-10 近空间高超声速飞行器运动控制仿真方法 Expired - Fee Related CN102073755B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010537312 CN102073755B (zh) 2010-11-10 2010-11-10 近空间高超声速飞行器运动控制仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010537312 CN102073755B (zh) 2010-11-10 2010-11-10 近空间高超声速飞行器运动控制仿真方法

Publications (2)

Publication Number Publication Date
CN102073755A true CN102073755A (zh) 2011-05-25
CN102073755B CN102073755B (zh) 2013-04-10

Family

ID=44032293

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010537312 Expired - Fee Related CN102073755B (zh) 2010-11-10 2010-11-10 近空间高超声速飞行器运动控制仿真方法

Country Status (1)

Country Link
CN (1) CN102073755B (zh)

Cited By (38)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101937233A (zh) * 2010-08-10 2011-01-05 南京航空航天大学 近空间高超声速飞行器非线性自适应控制方法
CN102320378A (zh) * 2011-06-20 2012-01-18 北京航空航天大学 多操纵面飞机的一种均衡操纵分配方法
CN102566446A (zh) * 2012-01-16 2012-07-11 南京航空航天大学 基于线性模型组的无人直升机全包线数学模型构建方法
CN102566427A (zh) * 2012-01-09 2012-07-11 南京航空航天大学 一种飞行器鲁棒控制方法
CN102707722A (zh) * 2012-05-31 2012-10-03 西北工业大学 基于飞行器常规模型的全维控制器区域设计方法
CN102707723A (zh) * 2012-05-31 2012-10-03 西北工业大学 基于飞行器常规模型的横航向控制器区域设计方法
CN102707629A (zh) * 2012-05-31 2012-10-03 西北工业大学 基于飞行器切换模型的全维控制器区域设计方法
CN103425135A (zh) * 2013-07-30 2013-12-04 南京航空航天大学 一种具有输入饱和的近空间飞行器鲁棒控制方法
CN103454927A (zh) * 2013-08-22 2013-12-18 哈尔滨工业大学 飞行器分布式网络化全物理地面仿真装置及方法
CN103473470A (zh) * 2013-09-24 2013-12-25 中国航天空气动力技术研究院 一种地面效应风洞试验数据处理方法
CN103837151A (zh) * 2014-03-05 2014-06-04 南京航空航天大学 一种四旋翼飞行器的气动模型辅助导航方法
CN104408243A (zh) * 2014-11-19 2015-03-11 中国航空工业集团公司沈阳飞机设计研究所 一种发动机模型安装推力修正方法
CN104634183A (zh) * 2014-12-18 2015-05-20 北京控制工程研究所 一种基于升阻比实时估计的自适应制导方法
CN105388763A (zh) * 2015-12-11 2016-03-09 北京航空航天大学 一种对流层间歇滑翔飞行控制方法
CN105653827A (zh) * 2016-03-17 2016-06-08 北京工业大学 高超声速飞行器Terminal滑模控制器设计方法
CN106019930A (zh) * 2016-08-03 2016-10-12 中国人民解放军63821部队 飞行器机动过程的气动/控制一体化耦合模拟技术
CN103970957B (zh) * 2014-05-20 2017-04-12 北京航空航天大学 一种弹性乘波体高超声速飞行器仿真方法
CN104331084B (zh) * 2014-09-30 2017-05-03 中国运载火箭技术研究院 一种基于方向舵控滚转策略的气动舵偏范围计算方法
CN106707790A (zh) * 2015-11-13 2017-05-24 成都飞机工业（集团）有限责任公司 一种无人机非线性数学模型建立方法
CN106802669A (zh) * 2017-03-16 2017-06-06 广西师范大学 一种基于速度矢量合成的航迹规划算法及其飞行器
CN106842912A (zh) * 2016-11-30 2017-06-13 南京航空航天大学 高超声速机动飞行抗舵面饱和鲁棒控制方法
CN108139425A (zh) * 2015-10-28 2018-06-08 欧姆龙株式会社 空速测量系统
CN108459611A (zh) * 2018-05-21 2018-08-28 济南大学 一种近空间飞行器的姿态跟踪控制方法
CN108958038A (zh) * 2018-08-16 2018-12-07 上海航天控制技术研究所 一种适应飞行器推力离散特性的控制参数调整方法
CN109145451A (zh) * 2018-08-22 2019-01-04 哈尔滨工业大学 一种高速滑翔飞行器的运动行为识别与航迹估计方法
CN109189087A (zh) * 2018-08-20 2019-01-11 哈尔滨工业大学 一种垂直起降重复使用运载器的自适应容错控制方法
CN109541963A (zh) * 2018-11-12 2019-03-29 北京应用气象研究所 一种基于侧滑角信息的无人机测风建模技术
CN109558660A (zh) * 2018-11-21 2019-04-02 中国航天空气动力技术研究院 一种航天器碎片陨落落区预报方法
CN109720535A (zh) * 2017-10-30 2019-05-07 成都飞机工业（集团）有限责任公司 一种翼身融合飞机
CN110456781A (zh) * 2019-09-16 2019-11-15 桂林航天工业学院 一种飞行器控制系统的空间稳定性分析方法
CN110471313A (zh) * 2019-08-26 2019-11-19 中仿智能科技（上海）股份有限公司 一种模拟飞行器的飞行仿真分系统
CN111290299A (zh) * 2020-03-20 2020-06-16 湖南云顶智能科技有限公司 一种高超声速试飞器助推段控制系统仿真平台
CN111752145A (zh) * 2019-03-26 2020-10-09 索尼公司 设计uav控制器模型的方法和系统、存储介质以及无人机
CN113252287A (zh) * 2021-07-19 2021-08-13 中国空气动力研究与发展中心高速空气动力研究所 一种加速度复合控制风洞半实物仿真试验方法
CN113342028A (zh) * 2021-07-26 2021-09-03 中国人民解放军63629部队 大型飞行器跟踪数字化仿真平台及其仿真方法
CN113721663A (zh) * 2021-10-29 2021-11-30 北京航空航天大学 一种柔性飞行器的起降轨迹规划方法
CN113985916A (zh) * 2021-10-26 2022-01-28 西安电子科技大学 基于压强闭环控制的飞行器变推力发动机控制分配方法、系统、装置及存储介质
CN118192303A (zh) * 2024-05-16 2024-06-14 成都流体动力创新中心 一种半实物仿真方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080073609A1 (en) * 2006-09-23 2008-03-27 Eldert Akkermann Compressed-air needle valve for controlling an air flow for driving engine simulators in aircraft models for wind tunnel experiments
CN101730437A (zh) * 2008-10-24 2010-06-09 南京理工大学 用于临近空间飞行器电子设备温控的辐射散热系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080073609A1 (en) * 2006-09-23 2008-03-27 Eldert Akkermann Compressed-air needle valve for controlling an air flow for driving engine simulators in aircraft models for wind tunnel experiments
CN101730437A (zh) * 2008-10-24 2010-06-09 南京理工大学 用于临近空间飞行器电子设备温控的辐射散热系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《计算机仿真》 20100131 张胜涛等 近空间高超飞行器气动热红外特性数值仿真 第114-118页 1 第27卷, 第1期 *

Cited By (61)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101937233A (zh) * 2010-08-10 2011-01-05 南京航空航天大学 近空间高超声速飞行器非线性自适应控制方法
CN101937233B (zh) * 2010-08-10 2012-05-30 南京航空航天大学 近空间高超声速飞行器非线性自适应控制方法
CN102320378A (zh) * 2011-06-20 2012-01-18 北京航空航天大学 多操纵面飞机的一种均衡操纵分配方法
CN102320378B (zh) * 2011-06-20 2013-07-24 北京航空航天大学 多操纵面飞机的一种均衡操纵分配方法
CN102566427A (zh) * 2012-01-09 2012-07-11 南京航空航天大学 一种飞行器鲁棒控制方法
CN102566446A (zh) * 2012-01-16 2012-07-11 南京航空航天大学 基于线性模型组的无人直升机全包线数学模型构建方法
CN102566446B (zh) * 2012-01-16 2014-05-14 南京航空航天大学 基于线性模型组的无人直升机全包线数学模型构建方法
CN102707629A (zh) * 2012-05-31 2012-10-03 西北工业大学 基于飞行器切换模型的全维控制器区域设计方法
CN102707722B (zh) * 2012-05-31 2014-04-16 西北工业大学 基于飞行器常规模型的全维控制器区域设计方法
CN102707629B (zh) * 2012-05-31 2013-12-04 西北工业大学 基于飞行器切换模型的全维控制器区域设计方法
CN102707723A (zh) * 2012-05-31 2012-10-03 西北工业大学 基于飞行器常规模型的横航向控制器区域设计方法
CN102707722A (zh) * 2012-05-31 2012-10-03 西北工业大学 基于飞行器常规模型的全维控制器区域设计方法
CN103425135A (zh) * 2013-07-30 2013-12-04 南京航空航天大学 一种具有输入饱和的近空间飞行器鲁棒控制方法
CN103454927A (zh) * 2013-08-22 2013-12-18 哈尔滨工业大学 飞行器分布式网络化全物理地面仿真装置及方法
CN103454927B (zh) * 2013-08-22 2016-06-08 哈尔滨工业大学 飞行器分布式网络化全物理地面仿真装置
CN103473470A (zh) * 2013-09-24 2013-12-25 中国航天空气动力技术研究院 一种地面效应风洞试验数据处理方法
CN103473470B (zh) * 2013-09-24 2016-08-03 中国航天空气动力技术研究院 一种地面效应风洞试验数据处理方法
CN103837151A (zh) * 2014-03-05 2014-06-04 南京航空航天大学 一种四旋翼飞行器的气动模型辅助导航方法
CN103837151B (zh) * 2014-03-05 2016-08-24 南京航空航天大学 一种四旋翼飞行器的气动模型辅助导航方法
CN103970957B (zh) * 2014-05-20 2017-04-12 北京航空航天大学 一种弹性乘波体高超声速飞行器仿真方法
CN104331084B (zh) * 2014-09-30 2017-05-03 中国运载火箭技术研究院 一种基于方向舵控滚转策略的气动舵偏范围计算方法
CN104408243B (zh) * 2014-11-19 2019-02-12 中国航空工业集团公司沈阳飞机设计研究所 一种发动机模型安装推力修正方法
CN104408243A (zh) * 2014-11-19 2015-03-11 中国航空工业集团公司沈阳飞机设计研究所 一种发动机模型安装推力修正方法
CN104634183A (zh) * 2014-12-18 2015-05-20 北京控制工程研究所 一种基于升阻比实时估计的自适应制导方法
CN108139425B (zh) * 2015-10-28 2020-11-20 欧姆龙株式会社 空速测量系统
US11008118B2 (en) 2015-10-28 2021-05-18 Omron Corporation Airspeed measurement system
CN108139425A (zh) * 2015-10-28 2018-06-08 欧姆龙株式会社 空速测量系统
CN106707790A (zh) * 2015-11-13 2017-05-24 成都飞机工业（集团）有限责任公司 一种无人机非线性数学模型建立方法
CN105388763B (zh) * 2015-12-11 2018-03-06 北京航空航天大学 一种对流层间歇滑翔飞行控制方法
CN105388763A (zh) * 2015-12-11 2016-03-09 北京航空航天大学 一种对流层间歇滑翔飞行控制方法
CN105653827B (zh) * 2016-03-17 2020-03-13 北京工业大学 高超声速飞行器Terminal滑模控制器设计方法
CN105653827A (zh) * 2016-03-17 2016-06-08 北京工业大学 高超声速飞行器Terminal滑模控制器设计方法
CN106019930A (zh) * 2016-08-03 2016-10-12 中国人民解放军63821部队 飞行器机动过程的气动/控制一体化耦合模拟技术
CN106842912A (zh) * 2016-11-30 2017-06-13 南京航空航天大学 高超声速机动飞行抗舵面饱和鲁棒控制方法
CN106842912B (zh) * 2016-11-30 2020-08-14 南京航空航天大学 高超声速机动飞行抗舵面饱和鲁棒控制方法
CN106802669A (zh) * 2017-03-16 2017-06-06 广西师范大学 一种基于速度矢量合成的航迹规划算法及其飞行器
CN109720535A (zh) * 2017-10-30 2019-05-07 成都飞机工业（集团）有限责任公司 一种翼身融合飞机
CN108459611B (zh) * 2018-05-21 2021-08-17 济南大学 一种近空间飞行器的姿态跟踪控制方法
CN108459611A (zh) * 2018-05-21 2018-08-28 济南大学 一种近空间飞行器的姿态跟踪控制方法
CN108958038A (zh) * 2018-08-16 2018-12-07 上海航天控制技术研究所 一种适应飞行器推力离散特性的控制参数调整方法
CN108958038B (zh) * 2018-08-16 2021-04-23 上海航天控制技术研究所 一种适应飞行器推力离散特性的控制参数调整方法
CN109189087A (zh) * 2018-08-20 2019-01-11 哈尔滨工业大学 一种垂直起降重复使用运载器的自适应容错控制方法
CN109145451A (zh) * 2018-08-22 2019-01-04 哈尔滨工业大学 一种高速滑翔飞行器的运动行为识别与航迹估计方法
CN109541963A (zh) * 2018-11-12 2019-03-29 北京应用气象研究所 一种基于侧滑角信息的无人机测风建模技术
CN109541963B (zh) * 2018-11-12 2021-06-25 北京应用气象研究所 一种基于侧滑角信息的无人机测风建模方法
CN109558660A (zh) * 2018-11-21 2019-04-02 中国航天空气动力技术研究院 一种航天器碎片陨落落区预报方法
CN109558660B (zh) * 2018-11-21 2023-10-03 中国航天空气动力技术研究院 一种航天器碎片陨落落区预报方法
CN111752145A (zh) * 2019-03-26 2020-10-09 索尼公司 设计uav控制器模型的方法和系统、存储介质以及无人机
US11984038B2 (en) 2019-03-26 2024-05-14 Sony Corporation Concept for designing and using an UAV controller model for controlling an UAV
CN110471313B (zh) * 2019-08-26 2022-07-22 中仿智能科技（上海）股份有限公司 一种模拟飞行器的飞行仿真分系统
CN110471313A (zh) * 2019-08-26 2019-11-19 中仿智能科技（上海）股份有限公司 一种模拟飞行器的飞行仿真分系统
CN110456781A (zh) * 2019-09-16 2019-11-15 桂林航天工业学院 一种飞行器控制系统的空间稳定性分析方法
CN110456781B (zh) * 2019-09-16 2021-08-24 桂林航天工业学院 一种飞行器控制系统的空间稳定性分析方法
CN111290299A (zh) * 2020-03-20 2020-06-16 湖南云顶智能科技有限公司 一种高超声速试飞器助推段控制系统仿真平台
CN113252287A (zh) * 2021-07-19 2021-08-13 中国空气动力研究与发展中心高速空气动力研究所 一种加速度复合控制风洞半实物仿真试验方法
CN113342028B (zh) * 2021-07-26 2022-08-09 中国人民解放军63629部队 大型飞行器跟踪数字化仿真平台及其仿真方法
CN113342028A (zh) * 2021-07-26 2021-09-03 中国人民解放军63629部队 大型飞行器跟踪数字化仿真平台及其仿真方法
CN113985916A (zh) * 2021-10-26 2022-01-28 西安电子科技大学 基于压强闭环控制的飞行器变推力发动机控制分配方法、系统、装置及存储介质
CN113985916B (zh) * 2021-10-26 2024-04-05 西安电子科技大学 基于压强闭环控制的飞行器变推力发动机控制分配方法、系统、装置及存储介质
CN113721663A (zh) * 2021-10-29 2021-11-30 北京航空航天大学 一种柔性飞行器的起降轨迹规划方法
CN118192303A (zh) * 2024-05-16 2024-06-14 成都流体动力创新中心 一种半实物仿真方法及系统

Also Published As

Publication number Publication date
CN102073755B (zh) 2013-04-10

Similar Documents

Publication Publication Date Title
CN102073755A (zh) 近空间高超声速飞行器运动控制仿真方法
CN103488814B (zh) 一种适用于再入飞行器姿态控制的闭环仿真系统
Lu et al. Real-time simulation system for UAV based on Matlab/Simulink
CN111399531B (zh) 高超声速飞行器滑翔段制导与姿态控制一体化设计方法
CN108090302B (zh) 一种直升机飞行力学模拟方法及系统
CN104007665A (zh) 一种固液动力飞行器飞行仿真测试系统
Sun et al. Design and implementation of a real-time hardware-in-the-loop testing platform for a dual-rotor tail-sitter unmanned aerial vehicle
Hua et al. Effect of elastic deformation on flight dynamics of projectiles with large slenderness ratio
CN107102565A (zh) 无人机集群软件在环仿真系统
Oruc et al. Towards real-time pilot-in-the-loop CFD simulations of helicopter/ship dynamic interface
CN114721261A (zh) 一种火箭子级姿态翻转着陆在线制导方法
Dalle et al. Flight envelope calculation of a hypersonic vehicle using a first principles-derived model
CN114706413A (zh) 近地轨道微纳卫星变质心姿态控制方法及系统
Cheng et al. A corridor-based flight mode transition strategy for agile ducted-fan tail-sitter UAV: Altitude-hold transition
Tunik et al. LMI-based synthesis of quadrotor guidance and control system
MacLean Modeling and simulation of a sounding rocket active stabilization system
Gonzalez et al. Linear control of highly flexible aircraft based on loop separation
Wuilbercq et al. Influence of boundary layer transition on the trajectory optimisation of a reusable launch vehicle
Ward et al. Performance characteristics of an autonomous airdrop system in realistic wind environments
Joseph et al. Towards Development of a Dynamic-Soaring, Morphing-Wing UAV: Flight Simulations
Zheng et al. Control-oriented modeling and analysis of the longitudinal dynamics of variable-sweep aircraft
Ou et al. Integrated flight dynamics modelling for unmanned aerial vehicles
Hall et al. An unsteady simulation technique for missile guidance and control applications
Guo et al. Time-accurate simulation of longitudinal flight mechanics with control by CFD/RBD coupling
Ritter et al. Comparison of Flight Test Data and Aeroelastic Simulation Results for DLR's Research Aircraft ATRA

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130410

Termination date: 20131110