CN112149234A - 一种基于俯仰角速率输入的飞行器质点运动模型设计方法 - Google Patents
一种基于俯仰角速率输入的飞行器质点运动模型设计方法 Download PDFInfo
- Publication number
- CN112149234A CN112149234A CN202011072038.1A CN202011072038A CN112149234A CN 112149234 A CN112149234 A CN 112149234A CN 202011072038 A CN202011072038 A CN 202011072038A CN 112149234 A CN112149234 A CN 112149234A
- Authority
- CN
- China
- Prior art keywords
- axis
- calculating
- angle
- elevator
- rate
- 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
- 239000002245 particle Substances 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 11
- 238000013461 design Methods 0.000 title abstract description 11
- 230000008859 change Effects 0.000 claims description 30
- 230000001133 acceleration Effects 0.000 claims description 16
- 238000013016 damping Methods 0.000 claims description 13
- 238000013178 mathematical model Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000012886 linear function Methods 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000005406 washing Methods 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract description 2
- 230000008878 coupling Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 238000005859 coupling reaction Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 206010034719 Personality change Diseases 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003028 elevating effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供了一种基于俯仰角速率输入的飞行器质点运动模型设计方法,利用现有的飞行器六自由度刚体运动方程,仅考虑俯仰运动,不考虑滚转与偏航运动,将原方程的输入由原来的升降舵变为俯仰角速率,将原微分方程中状态量俯仰角速率的特性用典型二阶系统来描述,并将飞行器的控制能力限制在角速率的微分方程中;本发明针对以迎角为输入的质点动力学方程的不足,建立一种基于俯仰角速率输入的质点运动模型,将飞行器控制能力的限制以及姿态运动的影响直接融入到质点动力学方程中,为轨迹规划和制导设计提供模型基础。
Description
技术领域
本发明涉及飞行器建模技术领域,主要涉及一种基于俯仰角速率输入的飞行器质点运动模型设计方法。
背景技术
随着航空航天技术的快速发展,飞行器的飞行速度越来越快、飞行高度越来越高、气动布局越来越先进、气动特性越来越复杂、飞行任务越来越多样,这都对控制系统的设计提出了严重的挑战。飞行器飞行过程中马赫数、迎角、高度、动压变化范围很大,不同状态下的气动特性差异大,飞行状态与动力系统耦合严重,姿态运动与质点运动耦合严重。研究姿态运动必须考虑质点运动的影响,研究质点运动也必须考虑姿态运动的影响,而姿态运动又取决于飞行器自身的控制能力。另一方面,受总体、结构、防热系统的限制,飞行器的控制能力有限,不同的飞行状态控制能力不同。在大迎角、大马赫数飞行阶段操纵能力明显不足,其质点运动受其控制能力的严格限制,同时姿态的变化对动力系统影响严重,轨迹规划及制导设计时,必须充分考虑其控制能力的限制,并将姿态变化的影响严格控制在允许的范围内。因此,需要建立考虑控制约束的质点运动模型,以适应姿态运动与质点运动之间的严重耦合。
目前,国内外轨迹规划及制导设计时普遍采用以迎角为输入的质点动力学方程,通过规划迎角规划轨迹剖面,并进行制导律的设计与仿真。文献“高超声速飞行器多约束再入轨迹快速优化”(2019,Vol.40(No.7):758~767)给出了以迎角和倾侧角为输入的高超声速飞行器无量纲三自由度运动方程,文献“Integration methods for aircraftscheduling and trajectory optimization at a busy terminal manoeuvring area”(2019,Vol.41(No.3):641~681)给出了以迎角和倾侧角为输入的无量纲三自由度运动方程。
以迎角为输入的质点动力学方程,虽然可以直接体现迎角的变化及变化率,对于无动力飞行器也可以间接反映对角速率的约束,但是在轨迹设计过程中,还存在一些问题。首先,对于带动力的飞行器,质点动力学方程中迎角的变化并不能直接体现姿态的变化,但是姿态的变化对飞行器的质点运动耦合严重;其次,质点动力学方程中没有体现对控制能力的约束,无法直接判断轨迹设计的合理性。因此,需要构建新的动力学方程,直接融入对姿态变化和控制能力的约束。
发明内容
发明目的:本发明提供了一种基于俯仰角速率输入的飞行器质点运动模型设计方法,针对以迎角为输入的质点动力学方程的不足,建立一种基于俯仰角速率输入的质点运动模型,适应面对称飞行器质点运动与姿态运动之间严重的耦合影响,将飞行器控制能力的限制以及姿态运动的影响直接融入到质点动力学方程中,为轨迹规划和制导设计提供模型基础。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种基于俯仰角速率输入的飞行器质点运动模型设计方法,其特征在于,包括以下步骤:
步骤S1、获取飞行器气动数据;所述气动数据包括机体轴x轴的力系数的基本项Cx0和升降舵产生的增量Cxc、机体轴z轴的力系数的基本项Cz0和升降舵产生的增量Czc、俯仰通道的稳定力矩系数Cm0和控制力矩系数Cmc、俯仰通道的俯仰阻尼导数Cmq和平尾下洗流时差阻尼导数
步骤S2、分别建立飞行器气动力系数和气动力矩系数的数学模型;
所述气动力系数基本项是马赫数Ma、迎角α和高度H的非线性函数,升降舵产生的气动力系数增量项是马赫数Ma、迎角α、高度H和升降舵δe的非线性函数,气动力系数表示如下:
Cx0=Cx0(Ma,α,H)
Cxc=Cxc(Ma,α,H,δe)
Cz0=Cz0(Ma,α,H)
Czc=Czc(Ma,α,H,δe)
所述气动力矩系数包括稳定力矩系数和控制力矩系数两部分,所述稳定力矩系数是马赫数Ma、迎角α和高度H的非线性函数;控制力矩系数是马赫数Ma、迎角α、高度H和升降舵δe的非线性函数;气动力矩系数表示如下:
Cm=Cm0(Ma,α,H)+Cmc(Ma,α,H,δe)
所述俯仰通道的阻尼导数是马赫数Ma和迎角α的非线性函数,表示如下:
步骤S3、获取飞行器推力数据,所述推力是时间t、马赫数Ma和高度H的非线性函数,具体表示如下:
T=T(t,Ma,H);
密度ρ和音速VS表示如下:
其中g为重力加速度,e为自然常数;
步骤S5、根据发动机推力方向与机体轴x轴的夹角η,计算推力在机体轴x轴和z轴上的分量Tx和Tz:
步骤S6、计算满足力矩平衡的升降舵配平舵面δe0;
在当前马赫数Ma、迎角α、高度H下,计算满足力矩平衡的升降舵配平舵面δe0,如下:
Cmc(Ma,α,H,δe0)=-Cm0(Ma,α,H);
步骤S7、计算机体系x轴和z轴方向的合力Fx和Fz如下:
其中,S为机翼参考面积;
根据升降舵的最大值δemax和最小值δemin,计算俯仰力矩最大值Mmax和最小值Mmin如下:
其中bA为机翼平均气动弦长;
其中Iyy为绕机体轴y轴的转动惯量;
其中,Δα为平衡状态下迎角的扰动量,Δδe为平衡状态下升降舵的增量;
步骤S10、计算俯仰角速率控制等效模型的频率ωn和阻尼ξ如下:
其中Kp为俯仰角速率反馈到升降舵的增益,Kα为迎角反馈到升降舵的增益;
其中Qc为等效模型的输入;
其中,R0为地球半径,ψ为偏航角,取固定值;
步骤S15、根据步骤S11-S14的计算结果,构造基于俯仰角速率输入的飞行器质点运动方程如下:
其中状态量x和输入量u分别为:
有益效果:本发明具备以下优点:
(1)以俯仰角速率为输入建立飞行器质点运动模型,将对角速率的约束直接融入到质点运动模型中,且通过对俯仰角速率的约束实现对迎角的约束。
(2)建立了俯仰角速率数学模型,将俯仰角速率的模型描述为典型的二阶系统,将控制能力的约束直接融入到质点运动模型中。
附图说明
图1是本发明提供的俯仰角速率控制结构图;
图2是本发明提供的飞行器质点运动模型设计流程图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
飞行器的刚体运动可以用微分方程描述如下:
其中,x,u分别表示系统的状态和输入。飞行器刚体运动方程分为运动学方程和动力学方程,运动学方程描述位置和速度的关系,动力学描述加速度和力/力矩的关系。
飞行器的纵向运动主要指俯仰运动和沿速度方向的线运动,在机体坐标系和地理坐标系下描述飞行器纵向运动的微分方程,飞行状态x包括俯仰角速率Q、俯仰角θ、偏航角ψ、机体系x轴的速度U、机体系z轴的速度W、高度H、经度l、纬度λ,输入u为升降舵δe。
下面通过线运动学方程、角运动学方程、线动力学方程和角动力学方程描述飞行器的纵向运动微分方程。
线运动学方程:
其中,R0为地球半径。
角运动学方程:
线动力学方程:
其中,Fx为机体系x轴方向的力,Fz为机体系z轴方向的力,m为飞行器质量,g为重力加速度。
角动力学方程:
其中,M为俯仰力矩,Iyy为绕机体轴y轴的转动惯量。
飞行器的合力是气动力和推力的总和,合力矩则为气动力矩和推力矩的总和:
其中,飞行器的气动力/力矩与当前飞行状态下的迎角α、马赫数Ma、高度H和升降舵δe有关。Tx为推力在机体轴x轴上的分量,Tz为推力在机体轴z轴上的分量。
仅考虑质点运动时,无法直接获取控制产生的升降舵,式(7)中俯仰力矩无法直接计算。因此本发明采用将输入量由升降舵改为俯仰角速率,俯仰角速率的动态特性用等效二阶系统描述,并将舵面带来的控制能力的约束体现在二阶系统的描述中。因为小扰动情况下,气动舵面变化对气动力的影响较小,直接用升降舵配平舵面计算气动力。
在当前飞行状态下进行配平和小扰动线性化,平衡状态满足:
在平衡状态下,进行小扰动线性化,得到线性化的动力学方程:
根据线性化方程可以得到传递函数:
考虑俯仰角速率控制律:
Δδe=KpΔQ+KαΔα (11)
图1给出了俯仰角速率控制结构图,其闭环系统传递函数可以描述为
其中,ξ和ωn分别为二阶环节的阻尼和频率。Mq为俯仰力矩对俯仰角速率Q的偏导数,为俯仰力矩对迎角变化率的偏导数,Mα为俯仰力矩对迎角α的偏导数,为俯仰力矩对升降舵δe的偏导数,Kp为俯仰角速率反馈到升降舵的增益,Kα为迎角反馈到升降舵的增益。
将闭环系统传递函数描述为微分方程的形式:
小扰动线性方程和配平状态叠加形成全量微分方程:
用式(14)替换式(6),式(3)、(4)、(5)、(14)和(16)共同构成基于俯仰角速率输入的飞行器质点运动方程。其中状态量和输入量分别为:
下面结合图2,以典型的面对称飞行器为例,阐述基于以俯仰角速率为输入的质点运动模型的实施方式。本发明中所提到的“升降舵”是指俯仰通道的所有控制舵面的统称,不是飞行器机体结构上单个的物理舵面,如左升降、右升降、左V尾、右V尾,左升降和右升降可以一起定义为“升降舵”,左V尾和右V尾也可以一起定义为“升降舵”。
步骤S1、获取飞行器气动数据;所述气动数据包括机体轴x轴的力系数的基本项Cx0和升降舵产生的增量Cxc、机体轴z轴的力系数的基本项Cz0和升降舵产生的增量Czc、俯仰通道的稳定力矩系数Cm0和控制力矩系数Cmc、俯仰通道的俯仰阻尼导数Cmq和平尾下洗流时差阻尼导数
步骤S2、分别建立飞行器气动力系数和气动力矩系数的数学模型;
所述气动力系数基本项是马赫数Ma、迎角α和高度H的非线性函数,升降舵产生的气动力系数增量项是马赫数Ma、迎角α、高度H和升降舵δe的非线性函数,气动力系数表示如下:
Cx0=Cx0(Ma,α,H)
Cxc=Cxc(Ma,α,H,δe)
Cz0=Cz0(Ma,α,H)
Czc=Czc(Ma,α,H,δe)
所述气动力矩系数包括稳定力矩系数和控制力矩系数两部分,所述稳定力矩系数是马赫数Ma、迎角α和高度H的非线性函数;控制力矩系数是马赫数Ma、迎角α、高度H和升降舵δe的非线性函数;气动力矩系数表示如下:
Cm=Cm0(Ma,α,H)+Cmc(Ma,α,H,δe)
所述俯仰通道的阻尼导数是马赫数Ma和迎角α的非线性函数,表示如下:
步骤S3、获取飞行器推力数据,所述推力是时间t、马赫数Ma和高度H的非线性函数,具体表示如下:
T=T(t,Ma,H);
密度ρ和音速VS表示如下:
其中g为重力加速度,e为自然常数;
步骤S5、根据发动机推力方向与机体轴x轴的夹角η,计算推力在机体轴x轴和z轴上的分量Tx和Tz:
步骤S6、计算满足力矩平衡的升降舵配平舵面δe0;
在当前马赫数Ma、迎角α、高度H下,计算满足力矩平衡的升降舵配平舵面δe0,如下:
Cmc(Ma,α,H,δe0)=-Cm0(Ma,α,H);
步骤S7、计算机体系x轴和z轴方向的合力Fx和Fz如下:
其中,S为机翼参考面积;
根据升降舵的最大值δemax和最小值δemin,计算俯仰力矩最大值Mmax和最小值Mmin如下:
其中bA为机翼平均气动弦长;
其中Iyy为绕机体轴y轴的转动惯量;
其中,Δα为平衡状态下迎角的扰动量,Δδe为平衡状态下升降舵的增量;
步骤S10、计算俯仰角速率控制等效模型的频率ωn和阻尼ξ如下:
其中Kp为俯仰角速率反馈到升降舵的增益,Kα为迎角反馈到升降舵的增益;
其中Qc为等效模型的输入;
其中,R0为地球半径,ψ为偏航角,取固定值;
步骤S15、根据步骤S11-S14的计算结果,构造基于俯仰角速率输入的飞行器质点运动方程如下:
其中状态量x和输入量u分别为:
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (1)
1.一种基于俯仰角速率输入的飞行器质点运动模型设计方法,其特征在于,包括以下步骤:
步骤S1、获取飞行器气动数据;所述气动数据包括机体轴x轴的力系数的基本项Cx0和升降舵产生的增量Cxc、机体轴z轴的力系数的基本项Cz0和升降舵产生的增量Czc、俯仰通道的稳定力矩系数Cm0和控制力矩系数Cmc、俯仰通道的俯仰阻尼导数Cmq和平尾下洗流时差阻尼导数
步骤S2、分别建立飞行器气动力系数和气动力矩系数的数学模型;
所述气动力系数基本项是马赫数Ma、迎角α和高度H的非线性函数,升降舵产生的气动力系数增量项是马赫数Ma、迎角α、高度H和升降舵δe的非线性函数,气动力系数表示如下:
Cx0=Cx0(Ma,α,H)
Cxc=Cxc(Ma,α,H,δe)
Cz0=Cz0(Ma,α,H)
Czc=Czc(Ma,α,H,δe)
所述气动力矩系数包括稳定力矩系数和控制力矩系数两部分,所述稳定力矩系数是马赫数Ma、迎角α和高度H的非线性函数;控制力矩系数是马赫数Ma、迎角α、高度H和升降舵δe的非线性函数;气动力矩系数表示如下:
Cm=Cm0(Ma,α,H)+Cmc(Ma,α,H,δe)
所述俯仰通道的阻尼导数是马赫数Ma和迎角α的非线性函数,表示如下:
步骤S3、获取飞行器推力数据,所述推力是时间t、马赫数Ma和高度H的非线性函数,具体表示如下:
T=T(t,Ma,H);
密度ρ和音速VS表示如下:
其中g为重力加速度,e为自然常数;
步骤S5、根据发动机推力方向与机体轴x轴的夹角η,计算推力在机体轴x轴和z轴上的分量Tx和Tz:
步骤S6、计算满足力矩平衡的升降舵配平舵面δe0;
在当前马赫数Ma、迎角α、高度H下,计算满足力矩平衡的升降舵配平舵面δe0,如下:
Cmc(Ma,α,H,δe0)=-Cm0(Ma,α,H);
步骤S7、计算机体系x轴和z轴方向的合力Fx和Fz如下:
其中,S为机翼参考面积;
根据升降舵的最大值δemax和最小值δemin,计算俯仰力矩最大值Mmax和最小值Mmin如下:
其中bA为机翼平均气动弦长;
其中Iyy为绕机体轴y轴的转动惯量;
其中,Δα为平衡状态下迎角的扰动量,Δδe为平衡状态下升降舵的增量;
步骤S10、计算俯仰角速率控制等效模型的频率ωn和阻尼ξ如下:
其中Kp为俯仰角速率反馈到升降舵的增益,Kα为迎角反馈到升降舵的增益;
其中Qc为等效模型的输入;
其中,R0为地球半径,ψ为偏航角,取固定值;
步骤S15、根据步骤S11-S14的计算结果,构造基于俯仰角速率输入的飞行器质点运动方程如下:
其中状态量x和输入量u分别为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011072038.1A CN112149234B (zh) | 2020-10-09 | 2020-10-09 | 一种基于俯仰角速率输入的飞行器质点运动模型设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011072038.1A CN112149234B (zh) | 2020-10-09 | 2020-10-09 | 一种基于俯仰角速率输入的飞行器质点运动模型设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112149234A true CN112149234A (zh) | 2020-12-29 |
CN112149234B CN112149234B (zh) | 2023-05-23 |
Family
ID=73952672
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011072038.1A Active CN112149234B (zh) | 2020-10-09 | 2020-10-09 | 一种基于俯仰角速率输入的飞行器质点运动模型设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112149234B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112965014A (zh) * | 2021-03-04 | 2021-06-15 | 哈尔滨工业大学 | 一种飞机机械结构变化引起磁干扰的补偿方法及装置 |
CN113607378A (zh) * | 2021-08-02 | 2021-11-05 | 厦门大学 | 一种绳系支撑飞行器模型强迫自由角运动模拟与抑制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1036721A (en) * | 1961-06-16 | 1966-07-20 | Smiths Industries Ltd | Improvements in or relating to automatic control systems for aircraft |
ITTO20080923A1 (it) * | 2008-12-11 | 2010-06-12 | Alenia Aeronautica Spa | Procedimento di stima dell'angolo di incidenza e dell'angolo di derapata di un aeromobile |
CN106707790A (zh) * | 2015-11-13 | 2017-05-24 | 成都飞机工业(集团)有限责任公司 | 一种无人机非线性数学模型建立方法 |
CN110263497A (zh) * | 2019-07-19 | 2019-09-20 | 南京航空航天大学 | 一种基于相对增益的气动耦合影响分析方法 |
CN110316358A (zh) * | 2019-03-29 | 2019-10-11 | 南京航空航天大学 | 基于动态逆的战斗机大迎角控制方法 |
CN111708377A (zh) * | 2020-06-21 | 2020-09-25 | 西北工业大学 | 基于惯导/飞控系统信息融合的飞行控制方法 |
-
2020
- 2020-10-09 CN CN202011072038.1A patent/CN112149234B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1036721A (en) * | 1961-06-16 | 1966-07-20 | Smiths Industries Ltd | Improvements in or relating to automatic control systems for aircraft |
ITTO20080923A1 (it) * | 2008-12-11 | 2010-06-12 | Alenia Aeronautica Spa | Procedimento di stima dell'angolo di incidenza e dell'angolo di derapata di un aeromobile |
CN106707790A (zh) * | 2015-11-13 | 2017-05-24 | 成都飞机工业(集团)有限责任公司 | 一种无人机非线性数学模型建立方法 |
CN110316358A (zh) * | 2019-03-29 | 2019-10-11 | 南京航空航天大学 | 基于动态逆的战斗机大迎角控制方法 |
CN110263497A (zh) * | 2019-07-19 | 2019-09-20 | 南京航空航天大学 | 一种基于相对增益的气动耦合影响分析方法 |
CN111708377A (zh) * | 2020-06-21 | 2020-09-25 | 西北工业大学 | 基于惯导/飞控系统信息融合的飞行控制方法 |
Non-Patent Citations (2)
Title |
---|
王建华;刘鲁华;汤国建;: "高超声速飞行器俯冲段制导与姿控系统设计" * |
高正红,焦天峰: "飞行器快速俯仰产生大迎角非定常气动力数学模型研究" * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112965014A (zh) * | 2021-03-04 | 2021-06-15 | 哈尔滨工业大学 | 一种飞机机械结构变化引起磁干扰的补偿方法及装置 |
CN113607378A (zh) * | 2021-08-02 | 2021-11-05 | 厦门大学 | 一种绳系支撑飞行器模型强迫自由角运动模拟与抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112149234B (zh) | 2023-05-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105159305B (zh) | 一种基于滑模变结构的四旋翼飞行控制方法 | |
CN102073755B (zh) | 近空间高超声速飞行器运动控制仿真方法 | |
CN112241125B (zh) | 一种基于微分平坦特性的无人机轨迹跟踪方法 | |
KR102294829B1 (ko) | 가변 하중을 가지는 무인비행체 비행제어 시스템 및 방법 | |
CN108090302B (zh) | 一种直升机飞行力学模拟方法及系统 | |
CN109703768B (zh) | 一种基于姿态/轨迹复合控制的软式空中加油对接方法 | |
CN111258216B (zh) | 一种适用于四旋翼飞行器的滑模重复控制器 | |
Wang et al. | Nonlinear aeroelastic control of very flexible aircraft using model updating | |
CN105759609B (zh) | 一种基于显式非线性模型预测控制的舰载机自主着舰方法 | |
CN109507890A (zh) | 一种基于eso的无人机动态逆广义预测控制器 | |
Pfeifle et al. | Cascaded incremental nonlinear dynamic inversion for three-dimensional spline-tracking with wind compensation | |
CN110316400B (zh) | 一种鸭翼布局固定翼无人机直接升力控制方法 | |
CN112149234A (zh) | 一种基于俯仰角速率输入的飞行器质点运动模型设计方法 | |
CN112631316B (zh) | 变负载四旋翼无人机的有限时间控制方法 | |
CN116300992A (zh) | 一种基于l1自适应动态逆的变体飞行器控制方法 | |
CN109582039A (zh) | 一种采用相对导航信息的j2摄动下最优队形重构方法 | |
Garcia et al. | Robust geometric navigation of a quadrotor UAV on SE (3) | |
Zhu et al. | Trajectory linearization control for a miniature unmanned helicopter | |
CN115329458A (zh) | 基于Fluent的海空多栖航行器空中飞行智能控制仿真方法 | |
CN113391555B (zh) | 基于非线性扰动观测器的四旋翼无人机降落控制方法 | |
Cheng et al. | A corridor-based flight mode transition strategy for agile ducted-fan tail-sitter UAV: Altitude-hold transition | |
CN116088549B (zh) | 一种尾座式垂直起降无人机姿态控制方法 | |
Zhong et al. | $ L_ {1} $ adaptive control of a dual-rotor tail-sitter unmanned aerial vehicle with input constraints during hover flight | |
CN115686036A (zh) | 一种基于预设性能的变外形飞行器多维复合控制方法 | |
CN112668092B (zh) | 一种耦合气动干扰的飞行器混合配平分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |