CN102707629B - 基于飞行器切换模型的全维控制器区域设计方法 - Google Patents
基于飞行器切换模型的全维控制器区域设计方法 Download PDFInfo
- Publication number
- CN102707629B CN102707629B CN2012101769430A CN201210176943A CN102707629B CN 102707629 B CN102707629 B CN 102707629B CN 2012101769430 A CN2012101769430 A CN 2012101769430A CN 201210176943 A CN201210176943 A CN 201210176943A CN 102707629 B CN102707629 B CN 102707629B
- Authority
- CN
- China
- Prior art keywords
- beta
- alpha
- centerdot
- cos
- sin
- 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.)
- Expired - Fee Related
Links
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于飞行器切换模型的全维控制器区域设计方法,用于解决现有的控制器设计方法不能直接确定给定飞行区域整体稳定性的技术问题。该方法通过气动力、力矩方程得到给定控制目标高度和马赫数时的平衡点,采用相平面分析模型确定***的区域稳定性,在此基础上确定反馈控制器的参数,直接对飞行器三维大迎角运动进行控制,避免了在力矩方程中忽略气动力作用等不正确近似,使得控制器在整个设计区域都能保证飞行器的稳定性,减少甚至避免了分析模型导致的不稳定、不安全飞行等问题发生。
Description
技术领域
本发明涉及一种飞行器控制器设计方法,特别涉及基于飞行器切换模型的全维控制器区域设计方法。
背景技术
飞行控制的基本目的是改善飞机的稳定性和操纵性,从而提高执行任务的能力;最近几十年来,随着飞机性能的不断提高,飞行控制技术发生了很大的变化,出现了主动控制技术、综合控制技术、自主飞行控制技术等先进的飞行控制技术,飞行控制***与航电***出现了高度综合化的趋势。现代高性能飞机对飞行控制***提出了更高的要求,使用古典控制理论设计先进飞机的飞行控制***已越来越困难;为了获得更好的飞行品质,许多现代控制方法被应用到飞机飞行控制***的设计中,如线性二次型调节器/线性二次型高斯函数/回路传递恢复(LQR/LQG/LTR)方法、定量反馈方法、动态逆方法、反馈线性化方法、反步控制方法、滑模变结构控制方法等;这些方法,需要飞行器准确的数学模型,然而,飞行器模型是一个很复杂的非线性微分方程式,人们很难得到准确的数学模型;工程上,飞机模型都是在通过风洞实验和飞行试验得到的,实际飞行控制***设计中还要考虑以下问题:(1)在已经建立起数学模型的飞机参数发生变化或存在结构不确定时,飞行控制***应该具有小的灵敏度响应;(2)由于控制器频带比较宽,使得飞机性能受飞机结构和执行机构动态性能变化的影响比较有小的灵敏度响应比较大;(3)反馈控制器的设计虽然对飞行员指令会得到较理想的响应,但是对于外部干扰的响应可能会是破坏性的;(4)执行部件与控制元件存在制造容差,***运行过程中也存在老化、磨损及环境和运行条件恶化等现象;(5)在实际工程问题中,通常对数学模型要人为地进行简化,去掉一些复杂的因素;为此,非线性H∞和μ综合鲁棒控制等非线性设计方法也在飞行控制器设计中得到广泛关注;上述方法,能够得到仅适于某个给定飞行状态的控制律结构及参数,在此基础上,需要逐次对整个飞行包线内不同飞行状态下的控制律设计,得到适于不同飞行状态的控制律结构和参数,并利用不同的方法进行控制律参数及结构的调整参数规律进行设计,最后得到一个适合于整个包线的完整的飞行控制律;依赖以上控制器设计方法,设计人员不能直接确定在给定飞行区域的稳定性;文献“Hsien-Keng Chenand Ching-I Lee,Anti-control of chaos in rigid body motion,Chaos,Solitons&Fractals,2004,Vol.21(4):957-965”直接根据飞行器通用的气动力、力矩表达式进行了相平面分析,既不考虑飞行器机型、又不考虑气动导数;论文方法偏离实际太远,给出的结果不被人们认可。
发明内容
为了克服现有控制器设计方法不能直接确定给定飞行区域整体稳定性的不足,本发明提供一种基于飞行器切换模型的全维控制器区域设计方法,该方法通过气动力、力矩方程得到给定控制目标高度、马赫数时的飞行器平稳平飞气流迎角和配平舵面,引入气流迎角、侧滑角等状态反馈控制器,采用相平面分析模型确定***的区域稳定性,在此基础上确定反馈控制器的参数,直接对飞行器三维大迎角运动进行控制,避免了在力矩方程中忽略气动力作用等不正确近似,使得控制器在整个设计区域都能保证飞行器的稳定性,减少甚至避免了分析模型导致的不稳定、不安全飞行等问题发生。
本发明解决其技术问题所采用的技术方案:一种基于飞行器切换模型的全维控制器区域设计方法,其特点是包括以下步骤:
1、根据方程:
其中:u,v,w分别为沿飞行器机体轴系x,y,z轴的速度分量;nx,ny,nz分别为沿x,y,z轴的过载;p,q,r分别为滚转、俯仰、偏航角速度;欧拉角θ,ψ分别指滚转、俯仰、偏航角;h为高度;g为重力加速度,当θ在360°k+90°(k=0,1,2,…)附近时,取
当θ在360°k-90°(k=0,1,2,…)附近时,取
和气动力、力矩模型
其中:q为俯仰角速度,α为气流迎角,β为侧滑角,θ为俯仰角,为滚转角,ψ为偏航角,p为滚转角速度,r为偏航角速度,g为重力加速度,δ为包含方向舵、副翼、升降舵、油门开度、鸭翼等在内的输入向量,Ix为绕轴x的转动惯量,Iy为绕轴y的转动惯量,Iz为绕轴z的转动惯量,Izx=Ixz为乘积转动惯量,V0为空速,Mpα(α,β)、Mrα(α,β)、Me(α,β,δ)为有关纵向力矩函数表达式, Lqβ(α,β), Nqβ(α,β), 为有关力矩函数表达式,nx,ny,nz分别为沿飞行器机体轴系x,y,z轴的过载;δs,αs,βs分别为对应控制目标高度、马赫数时的配平舵面、气流迎角、给定转弯半径稳定盘旋的侧滑角;全申请书符号相同;
2、选取反馈控制器表达式为:
δ=δ0+k(α,β,p,r,q)
其中:δ0为舵面输入的常数值,k(α,β,p,r,q)为反馈控制函数;
3、在给定飞行区域内,采用以下相平面分析模型:
分析***收敛性,根据收敛性指标和平衡点条件:满足条件:p=0,r=0,q=0, α=αs,β=βs时,δ=δs共同确定反馈控制器的参数;其中:定义如步骤1;aij(i=,p,r,q,j=l,m)为有关气动力系数。
本发明的有益效果是:通过气动力、力矩方程得到给定控制目标高度和马赫数时的平衡点,采用相平面分析模型确定***的区域稳定性,在此基础上确定反馈控制器的参数,直接对飞行器三维大迎角运动进行控制,避免了在力矩方程中忽略气动力作用等不正确近似,使得控制器在整个设计区域都能保证飞行器的稳定性,减少甚至避免了分析模型导致的不稳定、不安全飞行等问题发生。
下面结合实施例对本发明作详细说明。
具体实施方式
以某飞行器三维模型为例。
1、该飞行器三维模型中的气动力、力矩为:
gny/V0=-0.01r-0.40226β+0.0236β2-0.010221β3-0.035δr
gnz/V0=-0.877α+0.47α2+3.846α3-0.215δe
gnx/V0=-0.01265α+0.0047α3
-0.01859521β+0.002145291β3-0.2232δx=0
-0.0323β-0.1335δr=0
-4.208α-0.47α2-3.564α3-20.967δe+6.265α2δe=0;
2、选取反馈控制器表达式为:
δ=δ0+k(α,β,p,r,q)
3、在给定飞行区域-10°≤θ≤90°,-30°≤α≤80°,采用以下相平面分析模型:
分析***收敛性,根据收敛性指标和平衡点条件:满足条件:p=0,r=0,q=0, α=αs,β=βs时,δ=δs共同确定反馈控制器的参数为:δx=0.0961β3,δr=0,δe=-α3/(5.883-1.758α2);
Claims (1)
1.一种基于飞行器切换模型的全维控制器区域设计方法,其特征在于包括以下步骤:
(a)根据方程:
其中:u,v,w分别为沿飞行器机体轴系x,y,z轴的速度分量;nx,ny,nz分别为沿x,y,z轴的过载;p,q,r分别为滚转、俯仰、偏航角速度;欧拉角分别指滚转、俯仰、偏航角;h为高度;g为重力加速度,当θ在360°k+90°附近时,k=0,1,2,;…;取
并根据气动力、力矩模型
其中:q为俯仰角速度,α为气流迎角,β为侧滑角,θ为俯仰角,为滚转角,ψ为偏航角,p为滚转角速度,r为偏航角速度,g为重力加速度,δ为包含方向舵、副翼、升降舵、油门开度、鸭翼等在内的输入向量,Ix为绕轴x的转动惯量,Iy为绕轴y的转动惯量,Iz为绕轴z的转动惯量,Izx=Ixz为乘积转动惯量,V0为空速,Mpα(α,β)、Mrα(α,β)、Me(α,β,δ)为有关纵向力矩函数表达式, 为有关力矩函数表达式,nx,ny,nz分别为沿飞行器机体轴系x,y,z轴的过载;δs,αs,βs分别为对应控制目标高度、马赫数时的配平舵面、气流迎角、给定转弯半径稳定盘旋的侧滑角;
(b)选取反馈控制器表达式为:
δ=δ0+k(α,β,p,r,q)
其中:δ0为舵面输入的常数值,k(α,β,p,r,q)为反馈控制函数;
(c)在给定飞行区域内,采用以下相平面分析模型:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012101769430A CN102707629B (zh) | 2012-05-31 | 2012-05-31 | 基于飞行器切换模型的全维控制器区域设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012101769430A CN102707629B (zh) | 2012-05-31 | 2012-05-31 | 基于飞行器切换模型的全维控制器区域设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102707629A CN102707629A (zh) | 2012-10-03 |
CN102707629B true CN102707629B (zh) | 2013-12-04 |
Family
ID=46900525
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012101769430A Expired - Fee Related CN102707629B (zh) | 2012-05-31 | 2012-05-31 | 基于飞行器切换模型的全维控制器区域设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102707629B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103149930B (zh) * | 2013-03-24 | 2015-04-08 | 西安费斯达自动化工程有限公司 | 飞行器大迎角运动切换模型的故障诊断和容错控制方法 |
CN103197560A (zh) * | 2013-04-06 | 2013-07-10 | 西安费斯达自动化工程有限公司 | 飞行器三维飞行区域控制器宽适应性设计方法 |
CN105843079B (zh) * | 2016-05-30 | 2020-01-07 | 中国科学院光电技术研究所 | 一种多阶目标运动信息的估计方法 |
CN108303897B (zh) * | 2018-03-02 | 2020-12-29 | 西安费斯达自动化工程有限公司 | 飞行器颤振分析网格模型拉盖尔建模方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101625571A (zh) * | 2009-07-25 | 2010-01-13 | 大连理工大学 | 一种模拟自旋式飞行器六自由度运动的方法 |
CN102073755A (zh) * | 2010-11-10 | 2011-05-25 | 南京航空航天大学 | 近空间高超声速飞行器运动控制仿真方法 |
CN102176160A (zh) * | 2010-12-31 | 2011-09-07 | 北京控制工程研究所 | 基于偏量特征模型的高超声速飞行器自适应学习控制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7686253B2 (en) * | 2006-08-10 | 2010-03-30 | The Boeing Company | Systems and methods for tracing aircraft vortices |
US20090265049A1 (en) * | 2008-04-22 | 2009-10-22 | Honeywell International, Inc. | Aircraft system emissions and noise estimation mechanism |
-
2012
- 2012-05-31 CN CN2012101769430A patent/CN102707629B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101625571A (zh) * | 2009-07-25 | 2010-01-13 | 大连理工大学 | 一种模拟自旋式飞行器六自由度运动的方法 |
CN102073755A (zh) * | 2010-11-10 | 2011-05-25 | 南京航空航天大学 | 近空间高超声速飞行器运动控制仿真方法 |
CN102176160A (zh) * | 2010-12-31 | 2011-09-07 | 北京控制工程研究所 | 基于偏量特征模型的高超声速飞行器自适应学习控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102707629A (zh) | 2012-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102707624B (zh) | 基于飞行器常规模型的纵向控制器区域设计方法 | |
CN102163059B (zh) | 推力变向无人机的姿态控制***及控制方法 | |
CN102707723B (zh) | 基于飞行器常规模型的横航向控制器区域设计方法 | |
Huang et al. | Fixed-time autonomous shipboard landing control of a helicopter with external disturbances | |
CN104950898B (zh) | 一种再入飞行器全阶非奇异终端滑模姿态控制方法 | |
CN106444799B (zh) | 基于模糊扩张状态观测器和自适应滑模的四旋翼无人机控制方法 | |
Heidlauf et al. | Verification Challenges in F-16 Ground Collision Avoidance and Other Automated Maneuvers. | |
CN104950901B (zh) | 无人直升机姿态误差有限时间收敛非线性鲁棒控制方法 | |
CN104460681B (zh) | 倾转旋翼无人直升机过渡段的飞行控制方法 | |
CN102707616B (zh) | 基于飞行器三角模型的控制器区域设计方法 | |
CN102692928B (zh) | 基于飞行器四元数模型的控制器区域设计方法 | |
Chen et al. | Design of Flight Control System for a Novel Tilt‐Rotor UAV | |
Zhang et al. | Adaptive constrained backstepping controller with prescribed performance methodology for carrier-based UAV | |
CN105629734B (zh) | 一种近空间飞行器的轨迹跟踪控制方法 | |
CN102707722B (zh) | 基于飞行器常规模型的全维控制器区域设计方法 | |
CN102707629B (zh) | 基于飞行器切换模型的全维控制器区域设计方法 | |
CN104298109A (zh) | 基于多控制器融合的无尾飞行器协调转弯控制方法 | |
CN110320927A (zh) | 智能变形飞行器的飞行控制方法及*** | |
CN108595790A (zh) | 基于动态方程的高超声速飞行器耦合分析方法 | |
CN104881035A (zh) | 飞行器操纵耦合补偿方法、姿态运动控制方法及*** | |
CN104965515A (zh) | 一种基于自适应控制的螺旋桨无人靶机的姿态控制方法 | |
Wang et al. | A neural network based MRAC scheme with application to an autonomous nonlinear rotorcraft in the presence of input saturation | |
Li et al. | Simulation method for wind tunnel based virtual flight testing | |
CN104166348B (zh) | 一种单滑块变质心控制飞行器的动态稳定性判定方法 | |
Yang et al. | Non-linear position control for hover and automatic landing of unmanned aerial vehicles |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131204 Termination date: 20210531 |