CN110618694A - 一种空间拦截器鲁棒姿态控制方法 - Google Patents

一种空间拦截器鲁棒姿态控制方法 Download PDF

Info

Publication number
CN110618694A
CN110618694A CN201910853364.7A CN201910853364A CN110618694A CN 110618694 A CN110618694 A CN 110618694A CN 201910853364 A CN201910853364 A CN 201910853364A CN 110618694 A CN110618694 A CN 110618694A
Authority
CN
China
Prior art keywords
attitude
robust
control
control method
aircraft
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
CN201910853364.7A
Other languages
English (en)
Other versions
CN110618694B (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.)
Northern Information Control Research Institute Group Co Ltd
Original Assignee
Northern Information Control Research Institute Group Co Ltd
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 Northern Information Control Research Institute Group Co Ltd filed Critical Northern Information Control Research Institute Group Co Ltd
Priority to CN201910853364.7A priority Critical patent/CN110618694B/zh
Publication of CN110618694A publication Critical patent/CN110618694A/zh
Application granted granted Critical
Publication of CN110618694B publication Critical patent/CN110618694B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明提供了一种空间拦截器鲁棒姿态控制方法,包括以下步骤:步骤1,建立飞行器姿态模型步骤2,设计鲁棒非线性姿态稳定控制器。

Description

一种空间拦截器鲁棒姿态控制方法
技术领域
本发明涉及一种飞行器控制技术,特别是一种空间拦截器鲁棒姿态控制方法。
背景技术
由三个欧拉角的最少元素表示法,一个空间拦截器精确的姿态运动描述表示为一组多输入多输出、交叉耦合的非线性方程组。对于使用欧拉角描述法表示的空间拦截器姿态控制问题,现有的技术是采用backsteppng控制方法设计了非线性姿态控制律,然而并没有考虑在外部扰动变化缓慢的条件下姿态控制的鲁棒性,也没有考虑在外部扰动变化缓慢的条件下姿态稳定控制的稳态精度问题。
发明内容
本发明的目的在于提供一种空间拦截器鲁棒姿态控制方法。
实现本发明目的的技术方案为:一种空间拦截器鲁棒姿态控制方法,包括以下步骤:
步骤1,建立飞行器姿态模型
步骤2,设计鲁棒非线性姿态稳定控制器。
进一步地,步骤1中所述的飞行器姿态模型为欧拉角描述的姿态运动学模型
其中,θ表示俯仰角,ψ表示偏航角,γ表示倾斜角,ω=[ωx ωy ωz]T是体坐标系下角速度向量。
进一步地,步骤2中所述的鲁棒非线性姿态稳定控制器为
其中,J表示体坐标系下的名义转动惯量矩阵,α2是虚拟控制量,c3和c4是正常数,
x2=[θ ψ γ]T
x3=[ωx ωy ωz]T
z2=x21
z3=x32
α1、α2是虚拟控制量,c3和c4是正常数,Γ是一个正定对角矩阵,L为线性操作符。
本发明于backsteppng系统化控制技术,通过对空间飞行器存在变化缓慢外部扰动进行自适应估计,使得对空间飞行器进行姿态稳定控制时对变化缓慢外部扰动具有鲁棒性,由于在设计中引入了姿态欧拉角的积分,相对于以往的控制算法,可以提高在外部扰动变化缓慢的条件下姿态稳定控制的稳态精度,使用自适应控制方案对飞行器惯性参数进行估计,适用于空间飞行器惯性参数未知的情况。
下面结合说明书附图对本发明做进一步描述。
附图说明
图1为欧拉角稳定示意图,其中a)为俯仰角稳定示意图,b)为偏航角稳定示意图,c)为滚动角稳定示意图。
图2为角速度示意图,其中a)为角速度ωx示意图,b)为角速度ωy示意图,c)为角速度ωz示意图。
图3为控制力拒示意图,其中a)为控制力矩Tx示意图,b)为控制力矩Ty,c)为控制力矩Tz示意图。
图4为本发明的方法流程示意图。
具体实施方式
一种空间拦截器鲁棒姿态控制方法,包括以下步骤:
步骤1,建立飞行器姿态模型
步骤2,生化危机鲁棒非线性姿态稳定控制器。
具体的
步骤1具体过程如下:
步骤1.1,姿态动力学模型
考虑外部扰动力矩,一个刚体飞行器的姿态动力学模型由下式描述
式中,J=diag[Jx Jy Jz]表示体坐标系下的名义转动惯量矩阵;ω=[ωx ωyωz]T是体坐标系下角速度向量;u=[TxTyTz]T是控制力矩,可由反作用控制推进器产生;F(t)是外部扰动力矩向量,假设在本专利的设计中F(t)变化缓慢;操作符号ω×作用在ω=[ωx ωy ωz]T上形成反对称阵:
步骤1.2,欧拉角描述的姿态运动学模型
飞行器的姿态运动采用欧拉角描述的数学模型。用θ表示俯仰角,ψ表示偏航角,γ表示倾斜角。采用ψ→θ→γ(2→3→1)转序得到欧拉角描述的姿态运动学方程
这些用欧拉角描述的方程的的优点是物理意义明确,并且是飞行器姿态的最少元素表示法。但是,三个欧拉角的姿态表示有奇异性的问题,方程式和方程式显示奇异点在θ=±90°。因此,在某些实际情况下,某一个特定的旋转顺序可能是首选。在许多工程姿态控制问题中,这个姿态描述可以被采用,因为工作的俯仰角θ不接近±90°。从形式上看,这些方程对于ωx、ωy和ωz来说都是线性的,但是对于欧拉角来说是非线性的。
步骤2的具体过程为:
定义向量xθ=[θ ψ γ]T,x3=[ωx ωy ωz]T,x1=∫x2dt∈R3
总的姿态稳定系统表示成
式中
定义
z1=x1,z2=x21,z3=x32
式中,α1和α2是虚拟控制量,可以得到
α1=-c1x1
式中,c1是一个正常数。可以得到
式中,c2是一个正常数。可以得到
为F的估计值,F的估计误差为假设外部扰动力矩F(t)变化缓慢,取为系统构造一个Lyapunov函数
对上式沿着系统求导可得
设计鲁棒非线性姿态稳定控制器
式中,c3和c4是正常数。
从而可以得到
从图1可见,鲁棒非线性姿态稳定控制器提供了较好的瞬态过程,正如在图1(a)、图1(b)和图1(c)中所见,瞬态过程过后,飞行器的姿态稳定到0,稳定控制过程中欧拉角有较小的稳定超调。
图1表明鲁棒非线性姿态稳定控制器可以减轻变化缓慢外部干扰的影响,实现了较好的稳定控制性能。鲁棒非线性稳定控制律作用下的角速度和控制力矩分别绘制在图2和图3,由图可见,暂态过程过后,角速度趋近于0,说明姿态欧拉角稳定下来,变化很小,控制力矩变成小量,说明需要输出较小的力矩就能达到控制效果。
实施例
使用技术方案中系统设计的鲁棒非线性姿态稳定控制器,可以对空间拦截器进行自动控制。这部分将通过数值仿真分析说明具体实施方式和验证所提出控制算法的有效性。假设飞行器的转动惯量为Jx=0.3kgm2,Jy=Jz=2kgm2。飞行器的初始姿态角为θ(0)=40°,ψ(0)=-40°,γ(0)=10°,初始角速度为ωx(0)=0°/s,ωy(0)=0°/s,ωz(0)=0°/s。执行机构的饱和限制为5Nm。假设飞行器进行从一个静态到另一个静态的姿态机动,终端姿态角为0。由于在本专利的设计中假设F(t)变化缓慢,可以假设
F(t)=[0.15sin(0.1t)0.2sin(0.1t)0.1sin(0.1t)]TNm (9)
控制器中的参数取为c1=0.2,c2=11,c3=11,选取c4=30,滤波器的时间常数为τi=0.01s,i=1,2,3。
为了改善暂态过程,仿真中对欧拉角指令进行了规划,把暂态过程设计为斜坡信号(见图1)。

Claims (3)

1.一种空间拦截器鲁棒姿态控制方法,其特征在于,包括以下步骤:
步骤1,建立飞行器姿态模型
步骤2,设计鲁棒非线性姿态稳定控制器。
2.根据权利要求1所述的方法,其特征在于,步骤1中所述的飞行器姿态模型为欧拉角描述的姿态运动学模型
其中,θ表示俯仰角,ψ表示偏航角,γ表示倾斜角,ω=[ωx ωy ωz]T是体坐标系下角速度向量。
3.根据权利要求1所述的方法,其特征在于,步骤2中所述的鲁棒非线性姿态稳定控制器为
其中,J表示体坐标系下的名义转动惯量矩阵,α2是虚拟控制量,c3和c4是正常数,
x2=[θ ψ γ]T
x3=[ωx ωy ωz]T
z2=x21
z3=x32
α1、α2是虚拟控制量,c3和c4是正常数,Γ是一个正定对角矩阵,L为线性操作符。
CN201910853364.7A 2019-09-10 2019-09-10 一种空间拦截器鲁棒姿态控制方法 Active CN110618694B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910853364.7A CN110618694B (zh) 2019-09-10 2019-09-10 一种空间拦截器鲁棒姿态控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910853364.7A CN110618694B (zh) 2019-09-10 2019-09-10 一种空间拦截器鲁棒姿态控制方法

Publications (2)

Publication Number Publication Date
CN110618694A true CN110618694A (zh) 2019-12-27
CN110618694B CN110618694B (zh) 2022-08-12

Family

ID=68923203

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910853364.7A Active CN110618694B (zh) 2019-09-10 2019-09-10 一种空间拦截器鲁棒姿态控制方法

Country Status (1)

Country Link
CN (1) CN110618694B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6061611A (en) * 1998-01-06 2000-05-09 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Closed-form integrator for the quaternion (euler angle) kinematics equations
EP1901153A1 (en) * 2006-09-12 2008-03-19 OFFIS e.V. Control system for unmanned 4-rotor-helicopter
CN105468007A (zh) * 2015-12-02 2016-04-06 哈尔滨工业大学 一种基于干扰观测器的挠性卫星轨迹线性化姿态控制方法
CN107831658A (zh) * 2017-11-01 2018-03-23 北京航空航天大学 一种无人飞行器的鲁棒最优姿态控制方法
CN107992084A (zh) * 2017-12-27 2018-05-04 北京航空航天大学 不依靠角速度反馈的无人机鲁棒姿态控制方法和装置
CN108181807A (zh) * 2017-12-06 2018-06-19 北京航空航天大学 一种卫星初态阶段自适应容错姿态控制方法
CN108803648A (zh) * 2018-08-17 2018-11-13 北京航空航天大学 无人飞行器鲁棒姿态控制方法、装置及电子设备
CN109521786A (zh) * 2019-01-23 2019-03-26 南京邮电大学 基于比例积分观测器的四旋翼无人机滑模姿态控制方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6061611A (en) * 1998-01-06 2000-05-09 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Closed-form integrator for the quaternion (euler angle) kinematics equations
EP1901153A1 (en) * 2006-09-12 2008-03-19 OFFIS e.V. Control system for unmanned 4-rotor-helicopter
CN105468007A (zh) * 2015-12-02 2016-04-06 哈尔滨工业大学 一种基于干扰观测器的挠性卫星轨迹线性化姿态控制方法
CN107831658A (zh) * 2017-11-01 2018-03-23 北京航空航天大学 一种无人飞行器的鲁棒最优姿态控制方法
CN108181807A (zh) * 2017-12-06 2018-06-19 北京航空航天大学 一种卫星初态阶段自适应容错姿态控制方法
CN107992084A (zh) * 2017-12-27 2018-05-04 北京航空航天大学 不依靠角速度反馈的无人机鲁棒姿态控制方法和装置
CN108803648A (zh) * 2018-08-17 2018-11-13 北京航空航天大学 无人飞行器鲁棒姿态控制方法、装置及电子设备
CN109521786A (zh) * 2019-01-23 2019-03-26 南京邮电大学 基于比例积分观测器的四旋翼无人机滑模姿态控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周成宝: "飞行器姿态非线性控制方法研究", 《中国博士学位论文全文数据库》 *
周成宝: "飞行器姿态非线性控制方法研究", 《中国博士学位论文全文数据库》, no. 02, 31 December 2017 (2017-12-31), pages 1 - 108 *

Also Published As

Publication number Publication date
CN110618694B (zh) 2022-08-12

Similar Documents

Publication Publication Date Title
CN106773679B (zh) 一种基于角速度观测器的航天器容错控制方法
CN107121961B (zh) 一种基于迭代学习干扰观测器的航天器姿态容错控制方法
CN108319140B (zh) 一种重定义输出式无模型自适应航向控制方法及系统
CN109765920B (zh) 集成故障观测器和控制分配策略的航天器姿态容错的控制方法
Wu et al. Quantized fixed-time fault-tolerant attitude control for hypersonic reentry vehicles
Jasim et al. Integral backstepping controller for quadrotor path tracking
CN110442020A (zh) 一种新型的基于鲸鱼优化算法的容错控制方法
CN109709978A (zh) 一种高超声速飞行器制导控制一体化设计方法
Wu et al. Robust attitude maneuver control of spacecraft with reaction wheel low-speed friction compensation
CN115686048B (zh) 执行器受限航天器交会系统的动态触发有限时间控制方法
CN110618694B (zh) 一种空间拦截器鲁棒姿态控制方法
CN116923730B (zh) 一种具有自调节预设性能约束的航天器姿态主动容错控制方法
CN104460678B (zh) 一种基于Sigmoid型非线性滑模函数的航天器姿态控制方法
CN116466732B (zh) 一种抗振荡的模型参考自适应的飞行器滚转角控制方法
CN112305917A (zh) 充液航天器的固定时间终端滑模鲁棒容错控制方法及装置
CN110667890B (zh) 一种空间飞行器非线性姿态稳定方法
Xu et al. Singularity-free adaptive control of MIMO discrete-time nonlinear systems with general vector relative degrees
CN109032163B (zh) 一种考虑电机电压的四旋翼飞行器姿态控制方法
Jang et al. Ares i flight control system design
Khelfi et al. Robust control with sliding mode for a quadrotor unmanned aerial vehicle
CN113848958B (zh) 基于四元数的全驱动抗退绕水下机器人有限时间容错轨迹跟踪控制方法
CN115793447A (zh) 一种气动大不确定情况下的运载火箭分离体姿态控制方法
CN109799704A (zh) 一种挠性航天器姿态与结构振动的耦合控制方法
CN105912012A (zh) 一种部分执行机构故障下航天器混沌姿态控制方法
US6477433B1 (en) Control of velocity limited systems

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