CN106383523B - 一种飞行器非线性姿态控制系统稳定性分析方法 - Google Patents
一种飞行器非线性姿态控制系统稳定性分析方法 Download PDFInfo
- Publication number
- CN106383523B CN106383523B CN201610874908.4A CN201610874908A CN106383523B CN 106383523 B CN106383523 B CN 106383523B CN 201610874908 A CN201610874908 A CN 201610874908A CN 106383523 B CN106383523 B CN 106383523B
- Authority
- CN
- China
- Prior art keywords
- aircraft
- attitude control
- nonlinear
- threshold value
- dead zone
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
- G05D1/0816—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
Abstract
本发明公开了一种飞行器非线性姿态控制系统稳定性分析方法,属于飞行器姿态控制领域。该方法首先计算飞行器飞行状态对应的动力系数,用死区阈值h表达其中的对应项;之后绘制飞行器线性部分传递函数乃奎斯特曲线和非线性部分负倒描述函数曲线,若两曲线没有交汇点,则判断非线性系统稳定,设计俯仰通道死区阈值为h的取值,否则判断非线性系统不稳定;该方法给出了不同阈值条件下闭环系统的稳定性判别方法,从理论上提出了死区阈值对应的闭环回路稳定性判定依据。
Description
技术领域
本发明属于飞行器姿态控制领域,更具体地,涉及一种飞行器非线性姿态控制系统稳定性分析方法。
背景技术
飞行器姿态控制回路包含有死区的继电特性环节,是典型的非线性反馈控制系统。死区阈值h的选取直接影响控制回路的稳定性和动态跟踪性能。从提高姿态控制回路跟踪精度的角度考虑,要求选取较小的阈值,但是过小的阈值会导致控制回路不稳定;阈值选取过大会降低姿态控制动态跟踪性能及精度。
在实际工程设计过程中,广泛采用数学仿真调试设计的方法,多次设置不同的死区阈值分别进行仿真,分析仿真姿态控制结果,最终选取合理的死区阈值设置。通过这种仿真调试的方法能设计出满足姿态控制要求的死区阈值h,使得姿态控制回路具有比较理想的动态跟踪性能。但是这种方法需要进行多次仿真进行结果对比,是一种基于外部结果的设计方法,没有从系统内部分析出所设计闭环反馈系统的稳定性及稳定裕度。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种飞行器非线性姿态控制系统稳定性分析方法,其目的在于用死区阈值h表达动力系数,分析非线性闭环回路的稳定性,由此解决现有技术中需多次仿真及结果对比的问题。
为实现上述目的,按照本发明的一个方面,提供了一种飞行器非线性姿态控制系统稳定性分析方法,该方法包括以下步骤:
(1)采集飞行器飞行状态下的动压q、参考面积s、参考长度l、质量m、马赫数Ma、速度V、侧滑角β、攻角α、弹道倾角θ;
(2)计算飞行器飞行状态对应的俯仰通道动力系数b1、b2、c1和c2:
式中,Jz1为弹体绕OZ1轴的转动惯量;为俯仰通道阻尼力矩系数;为俯仰力矩系数对攻角导数;为升力系数对攻角导数;g为重力加速度;
由于非线性姿态控制被控对象Fctr_pLctr_p/hctr_p、Fctr_pcosα/hctr_p分别与线性被控对象具有一致的物理意义和数学量纲,均表示单位校正网络输出控制量产生的控制力矩,因此b3和c3可用死区阈值h进行定义
式中,Fctr_p为俯仰通道单个姿控喷管推力大小;Lctr_p为俯仰通道单个姿控喷管推力作用力臂;
(3)预选一个死区阈值h,代入到动力系数中,绘制飞行器线性部分传递函数乃奎斯特曲线和非线性部分负倒描述函数曲线,若两曲线没有交汇点,则判断非线性系统稳定,设计俯仰通道死区阈值为h的预选值;否则判断非线性系统不稳定。
进一步地,飞行器线性部分传递函数为
其中,Gjz(s)为校正网络传递函数;为机体姿态控制扰动量Δδp到俯仰姿态角扰动量的传递函数,表达式为
由此可以绘制出飞行器线性部分传递函数乃奎斯特曲线。
进一步地,非线性部分的输入输出关系如下:
根据非线性特性描述函数法,上述死区继电特性描述函数如下:
式中,自变量A为正弦输入信号幅值。基于线性部分动力系数b3、c3的定义,非线性特性描述函数中M=h;由此可以绘制出飞行器非线性部分负倒描述函数-1/N(A)的曲线。由极值条件dN(A)/dA=0,求解得到极大值点:
由闭环特征方程:
1+N(A)G1(jω)=0
看出引入非线性环节后,相当于在线性传递函数基础上乘以一个增益N(A)环节,其中N(A)为变增益实数,具有极大值N(Am)=2/π<1,因此,引入非线性环节后,相比原线性反馈回路会增大闭环回路的稳定裕度。
进一步地,该方法同样适用于偏航通道和滚动通道姿态控制系统稳定性分析。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下技术特征及有益效果:
(1)本发明方法给出了一种新的飞行器非线性姿态控制稳定性判定方法,进一步提出满足非线性姿态控制系统闭环稳定的死区阈值h设定方法;
(2)本发明方法给出了死区阈值h对应的非线性姿态控制系统稳定性判定理论性依据。
附图说明
图1是本发明方法流程图;
图2是h=3/180*π的乃奎斯特曲线和非线性部分负倒描述函数曲线图;
图3是h=3/180*π的跟踪误差响应图;
图4是h=0.5/180*π乃奎斯特曲线和非线性部分负倒描述函数曲线;
图5是h=0.5/180*π跟踪误差响应图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示,本发明包括以下步骤:
(1)采集飞行器飞行状态下的动压q、参考面积s、参考长度l、质量m、马赫数Ma、速度V、侧滑角β、攻角α、弹道倾角θ。
(2)计算飞行器飞行状态对应的俯仰通道动力系数b1、b2、c1和c2:
式中,Jz1为弹体绕OZ1轴的转动惯量;为俯仰通道阻尼力矩系数;为俯仰力矩系数对攻角导数;为升力系数对攻角导数;g为重力加速度;
由于非线性姿态控制被控对象Fctr_pLctr_p/hctr_p、Fctr_pcosα/hctr_p分别与线性被控对象具有一致的物理意义和数学量纲,均表示单位校正网络输出控制量产生的控制力矩,因此b3和c3可用死区阈值h的表示:
式中,Fctr_p为俯仰通道单个姿控喷管推力大小;Lctr_p为俯仰通道单个姿控喷管推力作用力臂。
(3)预选一个死区阈值h,代入到动力系数中,绘制飞行器线性部分传递函数乃奎斯特曲线和非线性部分负倒描述函数曲线,若两曲线没有交汇点,则判断非线性系统稳定,设计俯仰通道死区阈值为h的预选值;否则判断非线性系统不稳定。
其中,飞行器线性部分传递函数为
其中,Gjz(s)为校正网络传递函数;为机体姿态控制扰动量Δδp到俯仰姿态角扰动量的传递函数,表达式为
由此可以绘制出飞行器线性部分传递函数乃奎斯特曲线;
其中,非线性部分的输入输出关系如下:
根据非线性特性描述函数法,上述死区继电特性描述函数如下:
式中,自变量A为正弦输入信号幅值。基于线性部分动力系数b3、c3的定义,非线性特性描述函数中M=h;由此可以绘制出飞行器非线性部分负倒描述函数-1/N(A)的曲线。
由极值条件dN(A)/dA=0,求解得到极大值点:
由闭环特征方程:
1+N(A)G1(jω)=0
看出引入非线性环节后,相当于在线性传递函数基础上乘以一个增益N(A)环节,其中N(A)为变增益实数,具有极大值N(Am)=2/π<1,因此,引入非线性环节后,相比原线性反馈回路会增大闭环回路的稳定裕度。
实施例:
某飞行器一种飞行状态对应的俯仰通道动力系数数值如下:
b1=0;b2=0.6244;c1=0.0021;c2=-0.0015;
取h=3/180*π,绘制线性部分G1(s)乃奎斯特曲线和非线性部分负倒描述函数曲线-1/N(A)(A≥h),见图2,乃奎斯特曲线与非线性部分负倒描述函数曲线-1/N(A)没有交点,非线性系统稳定,振幅衰减。姿态控制系统对幅值为阶跃指令信号的跟踪误差响应见图3。
取h=0.5/180*π,绘制线性部分G1(s)乃奎斯特曲线和非线性部分负倒描述函数曲线-1/N(A)(A≥h),见图4,乃奎斯特曲线与非线性部分负倒描述函数曲线-1/N(A)之间存在交点,交点对应一个不稳定周期运动和一个稳定周期运动,当初始条件或外扰动使得振幅大于不稳定周期运动对应的振幅时,闭环系统将会产生自振荡,姿态控制系统对幅值为阶跃指令信号的跟踪误差响应见图5。
取h=3/180*π,姿态反馈控制系统不存在自振,闭环特征方程为
1+N(A)G1(jω)=0,
其中,N(A)取极大值N(Am)=2/π,计算得到反馈回路幅值裕度Gm=13.5dB,Pm=52deg,稳定裕度合理。进行数学仿真,动态跟踪特性理想并且跟踪精度满足设计要求,因此设计俯仰通道死区阈值为h=3/180*π,完成姿态反馈控制系统参数的设计和稳定性分析。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种飞行器非线性姿态控制系统稳定性分析方法,其特征在于,该方法包括以下步骤:
(1)采集飞行器飞行状态下的动压q、参考面积s、参考长度l、质量m、马赫数Ma、速度V、侧滑角β、攻角α、弹道倾角θ;
(2)计算飞行器飞行状态对应的俯仰通道动力系数b1、b2、c1和c2:
式中,Jz1为弹体绕OZ1轴的转动惯量;为俯仰通道阻尼力矩系数;为俯仰力矩系数对攻角导数;为升力系数对攻角导数;g为重力加速度;
b3和c3表示为含有死区阈值h的表达式:
式中,Fctr_p为俯仰通道单个姿控喷管推力大小;Lctr_p为俯仰通道单个姿控喷管推力作用力臂;
(3)预选一个死区阈值h,代入到动力系数中,绘制飞行器线性部分传递函数乃奎斯特曲线和非线性部分负倒描述函数曲线,若两曲线没有交汇点,则判断对应非线性系统稳定,设计俯仰通道死区阈值为h的预选值;否则判断对应非线性系统不稳定。
2.根据权利要求1所述的一种飞行器非线性姿态控制系统稳定性分析方法,其特征在于,该方法步骤(3)中飞行器线性部分传递函数为
其中,Gjz(s)为校正网络传递函数;为机体姿态控制扰动量Δδp到俯仰姿态角扰动量的传递函数,表达式为
由此可以绘制出飞行器线性部分传递函数乃奎斯特曲线。
3.根据权利要求1所述的一种飞行器非线性姿态控制系统稳定性分析方法,其特征在于,该方法步骤(3)中飞行器非线性部分负倒描述函数为-1/N(A),
其中,式中,M=h;自变量A为正弦输入信号幅值;由此可以绘制出飞行器非线性部分负倒描述函数曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610874908.4A CN106383523B (zh) | 2016-09-30 | 2016-09-30 | 一种飞行器非线性姿态控制系统稳定性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610874908.4A CN106383523B (zh) | 2016-09-30 | 2016-09-30 | 一种飞行器非线性姿态控制系统稳定性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106383523A CN106383523A (zh) | 2017-02-08 |
CN106383523B true CN106383523B (zh) | 2019-03-08 |
Family
ID=57936128
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610874908.4A Active CN106383523B (zh) | 2016-09-30 | 2016-09-30 | 一种飞行器非线性姿态控制系统稳定性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106383523B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108549219A (zh) * | 2018-03-28 | 2018-09-18 | 北京控制工程研究所 | 一种航天器喷气控制器参数辅助设计方法及系统 |
CN113641181B (zh) * | 2020-04-27 | 2023-05-26 | 华中科技大学 | 基于在线性能评估飞行器增益自适应姿态控制方法和系统 |
CN111813134B (zh) * | 2020-06-16 | 2024-02-02 | 湖北航天技术研究院总体设计所 | 飞行器控制系统的稳定性判别方法和判别系统 |
CN112325710B (zh) * | 2020-09-24 | 2023-03-31 | 北京航天自动控制研究所 | 一种运载火箭大推力直接入轨高精度姿态控制方法和系统 |
CN113467234B (zh) * | 2021-06-01 | 2024-04-09 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种确定飞机极限环振荡幅值与频率的方法及装置 |
CN114397074B (zh) * | 2022-01-20 | 2024-02-13 | 中山大学·深圳 | 一种正弦振动台控制方法、系统及装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000142360A (ja) * | 1998-11-05 | 2000-05-23 | Mazda Motor Corp | 車両の挙動制御装置 |
US6990396B2 (en) * | 2002-12-13 | 2006-01-24 | The Boeing Company | Method and apparatus for reaction wheel dynamic compensation in long-duration deployment of a large reflector |
CN100451898C (zh) * | 2005-12-14 | 2009-01-14 | 上海微小卫星工程中心 | 微小卫星的姿态控制方法及系统 |
CN101937233B (zh) * | 2010-08-10 | 2012-05-30 | 南京航空航天大学 | 近空间高超声速飞行器非线性自适应控制方法 |
US20150203215A1 (en) * | 2014-01-17 | 2015-07-23 | Eric T. Falangas | Early performance evaluation of conceptual flight and space vehicles |
CN105005196B (zh) * | 2015-05-14 | 2017-07-28 | 南京农业大学 | 农用车辆自主导航转向控制方法 |
CN105159311B (zh) * | 2015-09-14 | 2018-01-05 | 湖北航天技术研究院总体设计所 | 一种用于捷联导引头的自动驾驶仪的设计方法 |
-
2016
- 2016-09-30 CN CN201610874908.4A patent/CN106383523B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106383523A (zh) | 2017-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106383523B (zh) | 一种飞行器非线性姿态控制系统稳定性分析方法 | |
Theis et al. | Robust control design for active flutter suppression | |
CN109901606A (zh) | 一种用于四旋翼精确轨迹跟踪的混合有限时间控制方法 | |
CN106078742B (zh) | 一种针对带有输出约束的柔性机械臂的振动控制方法 | |
CN107065554B (zh) | 近空间可变翼飞行器自适应小翼切换控制系统及工作方法 | |
CN107491081A (zh) | 一种抗干扰四旋翼无人机姿态控制方法 | |
CN109062241B (zh) | 基于线性伪谱模型预测控制的自主全射向再入制导方法 | |
CN105182742A (zh) | 一种弹性体飞行器自适应受限跟踪控制间接法 | |
CN108241292B (zh) | 一种基于扩张状态观测器的水下机器人滑模控制方法 | |
CN108490788A (zh) | 一种基于双干扰观测的飞行器俯仰通道反演控制方法 | |
CN108459505B (zh) | 一种适合控制迭代设计的非常规布局飞行器快速建模方法 | |
CN109709978B (zh) | 一种高超声速飞行器制导控制一体化设计方法 | |
CN108983825A (zh) | 一种无人机时变编队的跟踪方法及系统 | |
Singh et al. | State feedback control of an aeroelastic system with structural nonlinearity | |
CN106708082B (zh) | 基于模糊控制的飞行器俯仰通道姿态指令快速跟踪方法 | |
CN105607472A (zh) | 非线性二元机翼的自适应反演滑模控制方法及装置 | |
CN106444807A (zh) | 一种栅格舵与侧喷流的复合姿态控制方法 | |
CN107263483B (zh) | 二自由度关节机器人轨迹的协调控制方法 | |
CN109703768B (zh) | 一种基于姿态/轨迹复合控制的软式空中加油对接方法 | |
CN109062055A (zh) | 一种基于Back-stepping鲁棒自适应动态面的近空间飞行器控制系统 | |
CN104122794A (zh) | 微陀螺仪的自适应模糊神经补偿非奇异终端滑模控制方法 | |
CN108180910A (zh) | 一种基于气动参数不确定的飞行器快速高精度制导方法 | |
CN110865539A (zh) | 一种随机干扰下的无人艇跟踪误差约束控制方法 | |
CN104597911A (zh) | 空中加油受油机自适应最优对接轨迹跟踪飞行控制方法 | |
CN106197173A (zh) | 基于扰动估计和补偿的战术导弹鲁棒姿态控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |