首页 > 专利 > 滨州学院 > 针对无人机舵机故障的姿态鲁棒自适应容错控制方法专利详情

针对无人机舵机故障的姿态鲁棒自适应容错控制方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2019-12-04
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2020-04-10
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2022-04-29
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2039-12-04
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201911224613.2 申请日 2019-12-04
公开/公告号 CN110888324B 公开/公告日 2022-04-29
授权日 2022-04-29 预估到期日 2039-12-04
申请年 2019年 公开/公告年 2022年
缴费截止日
分类号 G05B13/04 主分类号 G05B13/04
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 4
权利要求数量 5 非专利引证数量 1
引用专利数量 1 被引证专利数量 0
非专利引证 1、2019.03.28王志方 等.“电力巡检无人机的鲁棒自适应容错控制”《.武汉大学学报(工学版)》.2018,第51卷(第7期),第646-653页.;
引用专利 WO2019057634A 被引证专利
专利权维持 3 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 滨州学院 当前专利权人 滨州学院
发明人 郝伟 第一发明人 郝伟
地址 山东省滨州市滨城区黄河五路391号滨州学院 邮编 256600
申请人数量 1 发明人数量 1
申请人所在省 山东省 申请人所在市 山东省滨州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
济南泉城专利商标事务所 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
张贵宾
摘要
本发明提供了一种针对无人机舵机故障的姿态鲁棒自适应容错控制方法,包括:根据当前电机升力向量F、当前姿态角速度Ω以及故障估计值和按照预先设计的滑模自适应观测器对故障进行观测,确定并更新姿态角速度估计值根据姿态角速度估计值和当前姿态角速度Ω,计算并更新姿态角速度估计误差eΩ;根据当前姿态四元数q和目标姿态四元数qd,确定姿态误差四元数eq;按照预先设计的容错控制器,计算电机升力向量F,以按照电机升力向量控制调整无人机姿态;根据当前电机升力向量F、姿态角速度估计误差eΩ和姿态误差四元数eq,按照故障估计值和自适应律和估算更新故障估计值和通过本发明的技术方案,能够有效地克服无人机故障的影响,保持姿态稳定,大大地减少了计算量,提高了控制效率。
  • 摘要附图
    针对无人机舵机故障的姿态鲁棒自适应容错控制方法
  • 说明书附图:图1
    针对无人机舵机故障的姿态鲁棒自适应容错控制方法
  • 说明书附图:图2
    针对无人机舵机故障的姿态鲁棒自适应容错控制方法
  • 说明书附图:图3
    针对无人机舵机故障的姿态鲁棒自适应容错控制方法
  • 说明书附图:图4
    针对无人机舵机故障的姿态鲁棒自适应容错控制方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2022-04-29 授权
2 2020-04-10 实质审查的生效 IPC(主分类): G05B 13/04 专利申请号: 201911224613.2 申请日: 2019.12.04
3 2020-03-17 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.一种针对无人机舵机故障的姿态鲁棒自适应容错控制方法,适用于三旋翼无人机,其特征在于,包括:
当无人机舵机故障时,确定当前电机升力向量F、当前姿态四元数q、当前姿态角速度Ω以及目标姿态四元数qd;
根据所述当前电机升力向量F、所述当前姿态角速度Ω以及故障估计值 和 按照预先设计的滑模自适应观测器对故障进行观测,确定并更新姿态角速度估计值
根据所述姿态角速度估计值 和所述当前姿态角速度Ω,计算并更新姿态角速度估计误差eΩ;
根据所述当前姿态四元数q和所述目标姿态四元数qd,确定姿态误差四元数eq;
根据所述姿态角速度估计误差eΩ、所述姿态误差四元数eq,按照预先设计的容错控制器,计算电机升力向量,以按照所述电机升力向量控制调整无人机姿态;
根据所述当前电机升力向量F、所述姿态角速度估计误差eΩ和姿态误差四元数eq,按照预先设计的故障估计值 和 自适应律,估算更新故障估计值 和
定义惯性坐标系I、机体坐标系B和目标坐标系Bd,预先构建无人机舵机故障时的第一非线性动力学模型,
T 3×1
其中,Ω=[Ω1 Ω2 Ω3] ∈R ,表征为机体坐标系B相对于惯性坐标系I的姿态角速T
度,Ω1,Ω2,Ω3分别表征为滚转角速度、俯仰角速度和偏航角速度,[·]表示矩阵的转置,
3×1
∈表示集合间的“属于”关系,R 表示3行1列的实数向量, 表示求取Ω的一阶时间导数;
3×3
J=[J1 J2 J3]∈R 表征为转动惯量矩阵,J1,J2,J3分别表征为三旋翼无人机绕三个轴的转动惯量;×表示叉乘,Ω×JΩ表示求取Ω和JΩ的向量积;
为一系数矩阵,其中l表示前面某一电机中
心到无人机轴心的距离,l3表示舵机中心到无人机轴心的距离,α表示前面两个电机连线与某一电机和无人机轴心连线之间的夹角,k为电机的升力系数,δ为舵机发生堵塞时的偏转T 3×1
角度,l,l3,α,k均为已知常数,δ为未知常数;F=[F1 F2 F3]∈R 表征为电机升力向量,F1,F2,F3分别表征为故障发生后三个电机产生的升力,·表示点乘,Δ(δ)·F表示求取Δ(δ)和T 3×1
F的内积;D=[d1 d2 d3]∈R 为角速度变化引起的扰动系数;
定义故障变量β1,β2,β1=l3cosδ,β2=kcosδ+l3sinδ,变换所述第一非线性动力学模型的变量,得到含有所述故障变量β1,β2的第二非线性动力学模型,
根据所述第二非线性动力学模型,预先设计所述滑模自适应观测器,以对故障进行观测,
所述滑模自适应观测器为:
其中, 表征为姿态角速度估计值,也即表示对Ω的估计值, 表示求取 的一阶
T 3×1
时间导数,v=[v1 v2 v3]∈R ,∫vdt表示求取v关于时间的积分,定义Ω的估计误差,也即姿态角速度估计误差 故障估计值 分别表示故障变
1/2
量β1,β2的估计值, SIG1=[k11|eΩ1| sign(eΩ1)k12|
1/2 1/2 T T
eΩ2| sign(eΩ2)k13|eΩ3| sign(eΩ3)] ,v=[k21sign(eΩ1)k22sign(eΩ2)k23sign(eΩ3)] ,其中k11,k12,k13,k21,k22,k23均为正常数,eΩ1,eΩ2,eΩ3为姿态角速度估计误差eΩ的三个元素,
1/2
sign表示符号函数,|·| 表示求取绝对值的 次方,
所述根据所述当前姿态四元数q和目标姿态四元数qd,确定姿态误差四元数eq,包括:
采用基于单位四元数的姿态表示方法,机体坐标系B在惯性坐标系I下的表达用“等效轴
3×1
角坐标系”方法,将B和I重合,将B绕矢量k0∈R 按右手定则旋转 角,得到当前姿态四元数其中, 且满足 k0
3×1
∈R 为定义在惯性坐标系I中的任意单位矢量, 为机体坐标系B绕矢量k0旋转的任意角度;
由机体坐标系B到惯性坐标系I的坐标变换矩阵用单位四元数表示为
其中,I3为3×3的单位矩阵,S(qv)表示求取qv对
应的反对称矩阵;
采用基于单位四元数的姿态表示方法,目标坐标系Bd在惯性坐标系I下的表达用“等效
3×1
轴角坐标系”方法,将Bd和I重合,将Bd绕矢量kd∈R 按右手定则旋转 角,得到目标姿态四元数 其中 且满足
3×1
kd∈R 为定义在惯性坐标系I中的任意单位矢量, 为目标坐标系Bd绕矢
量kd旋转的任意角度;
由目标坐标系Bd到惯性坐标系I的坐标变换矩阵用四元数表示为
S(qvd)表示求取qvd对应的反对称矩阵;
定义姿态误差四元数 其中,e0和
ev满足
由目标坐标系Bd到机体坐标系B的坐标变换矩阵用四元数表示为
S(ev)表示求取ev对应的反对称矩阵,定义角速度
3×1
跟踪误差 其中Ωd∈R 表示目标坐标系Bd相对于惯性坐标系I的目标
姿态角速度,
还包括:
定义误差信号变量s,r,预先设计容错控制器,
3×3
s=eΩ+ksev, 其中,ks∈R , 均为常系数矩阵,
所述容错控制器为:
其中, 为控制器增益矩阵,diag{[σ1 σ2 σ3]}表示求取以σ1,σ2,
3×3
σ3为对角线元素的对角矩阵;kf∈R 为常系数矩阵,I3为3行3列的单位矩阵,r(t)为r在t时刻的取值,r(0)为r在0时刻的参数, 表示对(kf+I3)r(τ)+σ
sign(s(τ))关于τ从0到t积分; 表示求取矩阵 的逆矩阵。

2.根据权利要求1所述的针对无人机舵机故障的姿态鲁棒自适应容错控制方法,其特征在于,还包括:
预先设计故障估计值 和 自适应律,故障估计值 和 自适应律满足:
其中 和 分别表示 和 的自适应律,γ1,γ2,p1,p2均为
正常数, 分别表示F2,F3的导数。

3.根据权利要求1所述的针对无人机舵机故障的姿态鲁棒自适应容错控制方法,其特征在于,
所选控制器增益矩阵σ满足:

4.一种计算机设备,其特征在于,计算机设备包括处理器,处理器用于执行存储器中存储的计算机程序时实现如权利要求1至3中任一项所述的针对无人机舵机故障的姿态鲁棒自适应容错控制方法的步骤。

5.一种计算机可读存储介质,其特征在于,其上存储有计算机程序,计算机程序被处理器执行时实现如权利要求1至3中任一项所述的针对无人机舵机故障的姿态鲁棒自适应容错控制方法的步骤。
说明书

技术领域

[0001] 本发明涉及无人机控制技术领域,具体而言,涉及一种针对无人机舵机故障的姿态鲁棒自适应容错控制方法。

背景技术

[0002] 近年来,多旋翼无人机在高空摄影、灾后救援、环境监测等军事和民用领域得到越来越广泛的应用。与传统四旋翼无人机、六旋翼无人机不同,三旋翼无人机通常由三个电机和一个舵机组成,结构更简单、成本更低、能耗更小、机动性更强。三旋翼无人机依靠三个电机的转动及舵机的偏转实现俯仰、滚转、偏航等动作,受无人机飞行稳定性及自身工艺等影响,舵机极易发生堵塞故障,对无人机的飞行性能产生严重影响。
[0003] 目前国内外很多研究机构已经开始致力于三旋翼无人机动力学建模以及飞行控制的研究,但是针对舵机发生堵塞故障时的飞行控制研究尚没有相关文献涉及。三旋翼无人机作为一个四输入六输出的欠驱动系统,当舵机发生堵塞故障时,输入量减少一个,使得其欠驱动特性更为明显。如何针对三旋翼无人机发生舵机堵塞故障时进行姿态控制成为亟待解决的技术问题。

发明内容

[0004] 本发明旨在至少解决现有技术或相关技术中存在的技术问题之一。
[0005] 为此,本发明的一个目的在于提供一种针对无人机舵机故障的姿态鲁棒自适应容错控制方法。
[0006] 本发明的另一个目的在于提供一种计算机设备。
[0007] 本发明的再一个目的在于提供一种计算机可读存储介质。
[0008] 为了实现上述目的,本发明的第一方面的技术方案提供了一种针对无人机舵机故障的姿态鲁棒自适应容错控制方法,适用于三旋翼无人机,包括:当无人机舵机故障时,确定当前电机升力向量F、当前姿态四元数q、当前姿态角速度Ω以及目标姿态四元数qd;根据当前电机升力向量F、当前姿态角速度Ω以及故障估计值 和 按照预先设计的滑模自适应观测器对故障进行观测,确定并更新姿态角速度估计值 根据姿态角速度估计值和当前姿态角速度Ω,计算并更新姿态角速度估计误差eΩ;根据当前姿态四元数q和目标姿态四元数qd,确定姿态误差四元数eq;根据姿态角速度估计误差eΩ、姿态误差四元数eq,按照预先设计的容错控制器,计算电机升力向量,以按照电机升力向量控制调整无人机姿态;根据当前电机升力向量F、姿态角速度估计误差eΩ和姿态误差四元数eq,按照预先设计的故障估计值 和 自适应律,估算更新故障估计值 和
[0009] 在该技术方案中,通过预先设计的滑模自适应观测器对故障进行观测,有利于对三旋翼无人机舵机堵塞故障进行更有针对性的容错控制,采用基于单位四元数的姿态表示方法,用“等效轴角坐标系”方法等,有效避免了姿态表示奇异性问题,按照预先设计的容错控制器,计算电机升力向量,之后按照电机升力向量控制调整无人机姿态,既能对故障进行有效地抑制,而且又不需要主动容错控制所需要的故障隔离,大大地减少了计算量,提高了控制效率。经过仿真实验证明,本发明提出的针对无人机舵机故障的姿态鲁棒自适应容错控制方法具有较好的鲁棒性,当三旋翼无人机舵机发生堵塞故障时,无人机能够较快地克服故障影响,保持姿态稳定。
[0010] 需要说明的是,姿态角速度估计值 和故障估计值 都需要预设初始值,在滑模自适应观测器对故障进行观测,在控制过程中进行计算更新迭代,形成循环,当前电机升力向量F、当前姿态角速度Ω、当前姿态四元数q实时获取。
[0011] 在上述技术方案中,优选地,还包括:定义惯性坐标系{I}、机体坐标系{B}和目标坐标系{Bd},预先构建无人机舵机故障时的第一非线性动力学模型,T 3×1
其中,Ω=[Ω1 Ω2 Ω3]∈R ,表征为机体坐标系{B}相对于惯性坐标系{I}的姿态角速度,T
Ω1,Ω2,Ω3分别表征为滚转角速度、俯仰角速度和偏航角速度,[·]表示矩阵的转置,∈表示
3×1
集合间的“属于”关系,R 表示3行1列的实数向量, 表示求取Ω的一阶时间导数;J=[J1 J 2
3×3
J3]∈R 表征为转动惯量矩阵,J1,J2,J3分别表征为三旋翼无人机绕三个轴的转动惯量;×表示叉乘,Ω×JΩ表示求取Ω和JΩ的向量积; 为
一系数矩阵,其中l表示前面某一电机中心到无人机轴心的距离,l3表示舵机中心到无人机轴心的距离,α表示前面两个电机连线与某一电机和无人机轴心连线之间的夹角,k为电机的升力系数,δ为舵机发生堵塞时的偏转角度,l,l3,α,k均为已知常数,δ为未知常数;
[0012] F=[F1 F2 F3]T∈R3×1表征为电机升力向量,F1,F2,F3分别表征为故障发生后三个T 3×1电机产生的升力,·表示点乘,Δ(δ)·F表示求取Δ(δ)和F的内积;D=[d1 d2 d3]∈R 为角速度变化引起的扰动系数;定义故障变量β1,β2,β1=l3cosδ,β2=kcosδ+l3sinδ,变换第一非线性动力学模型的变量,得到含有故障变量β1,β2的第二非线性动力学模型,根据第二非线性动力学模型,预先设计滑模自适应观测
器,以对故障进行观测。
[0013] 在上述任一项技术方案中,优选地,滑模自适应观测器为:
[0014] 其中,其中 表征为姿态角速度T 3×1
估计值,也即表示对Ω的估计值, 表示求取 的一阶时间导数,v=[v1 v2 v3] ∈R ,∫vdt表示求取v关于时间的积分,定义Ω的估计误差,也即姿态角速度估计误差
故障估计值 分别表示故障变量β1,β2的估计值,
1/2 1/2
SIG1=[k11|eΩ1| sign(eΩ1) k12|eΩ2| sign(eΩ2) 
1/2 T T
k13|eΩ3| sign(eΩ3)] ,v=[k21sign(eΩ1) k22sign(eΩ2) k23sign(eΩ3)] ,其中k11,k12,k13,k21,k22,k23均为正常数,eΩ1,eΩ2,eΩ3为姿态角速度估计误差eΩ的三个元素,sign表示符号函
1/2
数,|·| 表示求取绝对值的 次方。
[0015] 在该技术方案中,通过分析舵机对三旋翼无人机的作用原理,并考虑角速度变化带来的扰动对其动力学特性的影响,得到三旋翼无人机执行器发生故障时的非线性动力学模型,来设计滑模自适应观测器,以对故障进行观测,滑模自适应观测器的观测性能更加,进一步有利于对三旋翼无人机舵机堵塞故障进行更有针对性的容错控制。
[0016] 在上述任一项技术方案中,优选地,根据当前姿态四元数q和目标姿态四元数qd,确定姿态误差四元数eq,包括:采用基于单位四元数的姿态表示方法,机体坐标系{B}在惯3×1
性坐标系{I}下的表达用“等效轴角坐标系”方法,将{B}和{I}重合,将{B}绕矢量k0∈R 按右手定则旋转 角,得到当前姿态四元数 其中,
3×1
且满足 k0∈R 为定义在惯性坐标系{I}
中的任意单位矢量, 为机体坐标系{B}绕矢量k0旋转的任意角度;由机体坐标系{B}到惯性坐标系{I}的坐标变换矩阵用单位四元数表示为
其中,I3为3×3的单位矩阵,S(qv)表示求取qv对应的反对称矩阵;采用基于单位四元数的姿态表示方法,目标坐标系{Bd}在惯性坐标系{I}下的表达用“等效轴角坐标系”方法,将{Bd}
3×1
和{I}重合,将{Bd}绕矢量kd∈R 按右手定则旋转 角,得到目标姿态四元数
其中 且满足
3×1
kd∈R 为定义在惯性坐标系{I}中的任意单位矢量, 为目标坐标系{Bd}绕矢量kd旋转的任意角度;由目标坐标系{Bd}到惯性坐标系{I}的坐标变换矩阵用四元数表示为
S(qvd)表示求取qvd对应的反对称矩阵;定义
姿态误差四元数 其中,e0和ev满足
由目标坐标系{Bd}到机体坐标系{B}的坐标变换矩阵用四元数表示为
S(ev)表示求取ev对应的反对称矩阵,定义角速度
3×1
跟踪误差 其中Ωd∈R 表示目标坐标系{Bd}相对于惯性坐标系{I}的
目标姿态角速度。
[0017] 在该技术方案中,采用基于单位四元数的姿态表示方法,用“等效轴角坐标系”方法等,有效避免了姿态表示奇异性问题,定义了姿态误差四元数eq,便于描述表达三旋翼无人机当前姿态与目标姿态之间的差异,便于定义误差信号变量,设计构建容错控制器,以对故障进行有效地抑制。
[0018] 在上述任一项技术方案中,优选地,还包括:定义误差信号变量s,r,预先设计容错3×3
控制器,s=eΩ+ksev, 其中,ks∈R , 均为常系数矩阵,
容错控制器为:
[0019]
[0020] 其中, 为控制器增益矩阵,diag{[σ1σ2σ3]}表示求取以σ1,3×3
σ2,σ3为对角线元素的对角矩阵;kf∈R 为常系数矩阵,I3为3行3列的单位矩阵,r(t)为r在t时刻的取值,r(0)为r在0时刻的参数, 表示对(kf+I3)r(τ)
+σsign(s(τ))关于τ从0到t积分; 表示求取矩阵 的逆矩阵。
[0021] 在该技术方案中,设计的容错控制器,采用的参数更少,大大减少了计算量,既能对故障进行有效地抑制,又不需要主动容错控制所需要的故障隔离,提高了控制效率。
[0022] 在上述任一项技术方案中,优选地,还包括:预先设计故障估计值 和 自适应律,故障估计值 和 自适应律满足:
[0023] 其中 和 分别表示 和 的自适应律,γ1,γ2,p1,p2均为正常数, 分别表示F2,F3的导数。
[0024] 在该技术方案中,设计的故障估计值 和 自适应律,自适应性较强,稳定性较强,有利于无人机克服故障影响,保持姿态稳定。
[0025] 在上述任一项技术方案中,优选地,所选控制器增益矩阵σ满足:
[0026]
[0027] 在该技术方案中,采用该控制器增益矩阵,使得所设计的容错控制器能够让姿态跟踪误差四元数和角速度跟踪误差半全局渐近收敛到0,进一步有利于对故障进行有效地抑制,使得无人机能够较快地克服故障影响,保持姿态稳定。
[0028] 本发明的第二方面的技术方案提出了一种计算机设备,计算机设备包括处理器,处理器用于执行存储器中存储的计算机程序时实现如上述本发明的第一方面的技术方案提出的任一项的针对无人机舵机故障的姿态鲁棒自适应容错控制方法的步骤。
[0029] 在该技术方案中,计算机设备包括处理器,处理器用于执行存储器中存储的计算机程序时实现如上述本发明的第一方面的技术方案提出的任一项的针对无人机舵机故障的姿态鲁棒自适应容错控制方法的步骤,因此具有上述本发明的第一方面的技术方案提出的任一项的针对无人机舵机故障的姿态鲁棒自适应容错控制方法的全部有益效果,在此不再赘述。
[0030] 本发明的第三方面的技术方案提出了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现本发明的第一方面的技术方案提出的任一项的针对无人机舵机故障的姿态鲁棒自适应容错控制方法的步骤。
[0031] 在该技术方案中,计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现本发明的第一方面的技术方案提出的任一项的针对无人机舵机故障的姿态鲁棒自适应容错控制方法的步骤,因此具有上述本发明的第一方面的技术方案提出的任一项的针对无人机舵机故障的姿态鲁棒自适应容错控制方法的全部有益效果,在此不再赘述。
[0032] 通过以上技术方案,针对三旋翼无人机发生舵机堵塞故障时的姿态控制问题采用基于观测器技术的方法,该容错控制方法既能对故障进行有效地抑制,而且又不需要主动容错控制所需要的故障隔离,大大地减少了计算量,提高了控制效率。经过仿真实验证明,本发明提出的针对无人机舵机故障的姿态鲁棒自适应容错控制方法具有较好的鲁棒性,当三旋翼无人机舵机发生堵塞故障时,无人机能够较快地克服故障影响,保持姿态稳定。
[0033] 本发明的附加方面和优点将在下面的描述部分中给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。

实施方案

[0039] 为了能够更清楚地理解本发明的上述目的、特征和优点,下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是,在不冲突的情况下,本申请的实施例及实施例中的特征可以相互组合。
[0040] 在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用其他不同于在此描述的其他方式来实施,因此,本发明的保护范围并不受下面公开的具体实施例的限制。
[0041] 根据本发明的实施例的针对无人机舵机故障的姿态鲁棒自适应容错控制方法。
[0042] 首先,定义惯性坐标系{I}、机体坐标系{B}和目标坐标系{Bd},通过分析舵机对三旋翼无人机的作用原理,并考虑角速度变化带来的对其动力学特性的影响,得到三旋翼无人机执行器发生故障时的非线性动力学模型:
[0043]
[0044] 式(1)中各变量定义如下:Ω=[Ω1 Ω2 Ω3]T∈R3×1表示机体坐标系{B}相对于惯性坐标系{I}的姿态角速度,Ω1,Ω2,Ω3分别表示滚转角速度、俯仰角速度和偏航角速度,T 3×1[·]表示矩阵的转置,∈表示集合间的“属于”关系,R 表示3行1列的实数向量, 表示求
3×3
取Ω的一阶时间导数;J=[J1 J2 J3]∈R 为转动惯量矩阵,J1,J2,J3分别为三旋翼无人机绕三个轴的转动惯量;×表示叉乘,Ω×JΩ表示求取Ω和JΩ的向量积;
为一系数矩阵,其中l3表示前面某一电机中心
到无人机轴心的距离,α表示舵机中心到无人机轴心的距离,表示前面两个电机连线与某一电机和无人机轴心连线之间的夹角,k为电机的升力系数,δ为舵机发生堵塞时的偏转角度,l,l3,α,k均为已知常数,δ为未知常数,sin和cos分别表示正弦函数和余弦函数;F=[F1 F2 T 3×1
F3]∈R 表示升力向量,F1,F2,F3分别表示故障发生后三个电机产生的升力;·表示点乘,T 3×1
Δ(δ)·F表示求取Δ(δ)和F的内积;D=[d1 d2 d3]∈R 为角速度变化引起的扰动系数。
[0045] 其次,定义变量β1,β2分别为:
[0046] β1=l3cosδ,   (2)
[0047] β2=kcosδ+l3sinδ,   (3)
[0048] 由于δ为未知常数,则β1,β2为未知常数,Δ(δ)可写为式(1)可表示为
[0049]
[0050] 为避免姿态表示奇异性问题,采用基于单位四元数的姿态表示方法,机体坐标系{B}在惯性坐标系{I}下的表达用“等效轴角坐标系”方法,将{B}和{I}重合,将{B}绕矢量k03×1
∈R 按右手定则旋转 角,得到当前姿态四元数 其中,
3×1
且满足 k0∈R 为定义在坐标系{I}中
的任意单位矢量, 为坐标系{B}绕矢量k0旋转的任意角度;由机体坐标系{B}到惯性坐标系{I}的坐标变换矩阵用四元数表示为 I3为3×3的单位矩
阵,S(qv)表示求取qv对应的反对称矩阵。同理,目标坐标系{Bd}在惯性坐标系{I}下的表达也用
3×1
“等效轴角坐标系”方法,将{Bd}和{I}重合,将{Bd}绕矢量kd∈R 按右手定则旋转 角,得到目标姿态四元数 其中
3×1
且满足 kd∈R 同样为定义在坐标系{I}中的任意单位矢量, 为坐标系{Bd}绕
矢量kd旋转的任意角度;由目标坐标系{Bd}到惯性坐标系{I}的坐标变换矩阵用四元数表示为S(qvd)表示求取qvd对应的反对称矩阵,为了描
述三旋翼无人机当前姿态与目标姿态之间的差异,定义姿态误差四元数
[0051]
[0052] 其中e0和ev同样满足 由目标坐标系{Bd}到机体坐标系{B}的坐标变换矩阵示为 S(ev)表示求取ev对应的反对称矩阵,定
3×1
义角速度跟踪误差 其中Ωd∈R 表示目标坐标系{Bd}相对于惯性坐
标系{I}的姿态角速度;
[0053] 为了对三旋翼无人机舵机堵塞故障进行更有针对性的容错控制,设计滑模自适应观测器对故障进行观测:
[0054]T
[0055] 其中 表示对Ω的估计值, 表示求取 的一阶时间导数,v=[v1 v2 v3] ∈3×1
R ,∫vdt表示求取v关于时间的积分,定义Ω的估计误差为
1/
分别表示β1,β2的估计值, SIG1=[k11|eΩ1|
2 1/2 1/2 T
sign(eΩ1) k12|eΩ2| sign(eΩ2) k13|eΩ3| sign(eΩ3)] ,v=[k21sign(eΩ1) k22sign(eΩ2) T
k23sign(eΩ3)] ,其中k11,k12,k13,k21,k22,k23均为正常数,eΩ1,eΩ2,eΩ3为姿态角速度估计误
1/2
差eΩ的三个元素,sign表示符号函数,|·| 表示求取绝对值的 次方。
[0056] 然后,定义新的误差信号变量:
[0057] s=eΩ+ksev   (7)
[0058]3×3
[0059] 其中ks∈R , 均为常系数矩阵。
[0060] 最后,设计容错控制器为:
[0061]
[0062] 其中 为控制器增益矩阵,diag{[σ1 σ2 σ3]}表示求取以σ1,3×3
σ2,σ3为对角线元素的对角矩阵;kf∈R 为常系数矩阵,I3为3行3列的单位矩阵,r(t)为r在t时刻的取值,r(0)为r在0时刻的参数, 表示对(kf+I3)r(τ)
+σsign(s(τ))关于τ从0到t积分。 表示求取矩阵 的逆矩阵, 和 的自
适应律满足:
[0063]
[0064] 其中 和 分别表示 和 的自适应律,γ1,γ2,p1,p2均为正常数, 分别表示F2,F3的导数。
[0065] 如果所选控制器增益σ满足:
[0066]
[0067] 那么,所设计的控制器能够使得姿态跟踪误差四元数和角速度跟踪误差半全局渐近收敛到0。
[0068] 验证步骤具体是,采用基于Lyapunov李亚普诺夫稳定性的分析方法可以证明当时T间趋于无穷时,eΩ和ev分别渐近收敛到[0 0 0]。
[0069] 当无人机舵机故障时,确定当前电机升力向量F、当前姿态四元数q、当前姿态角速度Ω以及目标姿态四元数qd;根据当前电机升力向量F、当前姿态角速度Ω以及故障估计值和 按照预先设计的滑模自适应观测器对故障进行观测,确定并更新姿态角速度估计值 根据姿态角速度估计值 和当前姿态角速度Ω,计算并更新姿态角速度估计误差eΩ;根据当前姿态四元数q和目标姿态四元数qd,确定姿态误差四元数eq;根据姿态角速度估计误差eΩ、姿态误差四元数eq,按照预先设计的容错控制器,计算电机升力向量,以按照电机升力向量控制调整无人机姿态;根据当前电机升力向量F、姿态角速度估计误差eΩ和姿态误差四元数eq,按照预先设计的故障估计值 和 自适应律,估算更新故障估计值 和[0070] 需要说明的是,姿态角速度估计值 和故障估计值 都需要预设初始值,在滑模自适应观测器对故障进行观测,在控制过程中进行计算更新迭代,形成循环,当前电机升力向量F、当前姿态角速度Ω、当前姿态四元数q实时获取。
[0071] 采用本发明提出的针对无人机舵机故障的姿态鲁棒自适应容错控制方法进行仿真试验,涉及的各参数取值如下:
[0072] J=diag{[0.01 0.01 0.02]T}kg·m2,l=0.16m,l3=0.2m,α=26.56°,k=0.005,T T TD=diag{[0.1 0.1 0.1] },ks=diag{[3 5 5]},kr=diag{[0.1 0.1 0.1]},kf=diagT T
{[0.15 0.85 0.95]},σ=diag{[1 1 4]},k11=0.0009,k12=0.1,k13=0.1,k21=0.0005,k22=0.1,k23=0.2,γ1=0.2,γ2=0.1,p1=0.005,p2=1000。
[0073] 初始姿态四元数和角速度分别为q0=[0.9961 0.04 ‑0.05 0.06]T,ω0=[2.11 T2.3 2.2]deg/s。三旋翼无人机的舵机堵塞角约为8°。仿真结果分别如图1至图4所示。
[0074] 如图1所示,从三旋翼无人机姿态误差四元数的变化曲线,可见,舵机发生堵塞故障时,无人机姿态发生变化,并在5s内迅速收敛至0。如图2所示,从姿态角速度估计误差的变化曲线,可见,故障发生后,其在3s之内迅速收敛到0。由此可见,控制目标得到很好的实现,该方法既能对故障进行有效地抑制,无人机能够较快地克服故障影响,保持姿态稳定。如图3所示,从控制输入变化曲线可见,升力在合理变化范围内。如图4所示,从舵机堵塞故障估计值变化曲线,可见,为稳定状态,与理论计算结果相符。综上可见,本发明提出的针对无人机舵机故障的姿态鲁棒自适应容错控制方法具有有效性。
[0075] 以上结合附图详细说明了本发明的技术方案,本发明提出了一种针对无人机舵机故障的姿态鲁棒自适应容错控制方法,针对三旋翼无人机发生舵机堵塞故障时的姿态控制问题采用基于观测器技术的方法,该容错控制方法既能对故障进行有效地抑制,而且又不需要主动容错控制所需要的故障隔离,大大地减少了计算量,提高了控制效率。经过仿真实验证明,本发明提出的针对无人机舵机故障的姿态鲁棒自适应容错控制方法具有较好的鲁棒性,当三旋翼无人机舵机发生堵塞故障时,无人机能够较快地克服故障影响,保持姿态稳定。
[0076] 本发明方法中的步骤可根据实际需要进行顺序调整、合并和删减。
[0077] 本发明装置中的单元可根据实际需要进行合并、划分和删减。
[0078] 本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存储介质中,存储介质包括只读存储器(Read‑Only Memory,ROM)、随机存储器(Random Access Memory,RAM)、可编程只读存储器(Programmable Read‑only Memory,PROM)、可擦除可编程只读存储器(Erasable Programmable Read Only Memory,EPROM)、一次可编程只读存储器(One‑time Programmable Read‑Only Memory,OTPROM)、电子抹除式可复写只读存储器(Electrically‑Erasable Programmable Read‑Only Memory,EEPROM)、只读光盘(Compact Disc Read‑Only Memory,CD‑ROM)或其他光盘存储器、磁盘存储器、磁带存储器、或者能够用于携带或存储数据的计算机可读的任何其他介质。
[0079] 以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

附图说明

[0034] 本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
[0035] 图1示出了三旋翼无人机姿态误差四元数的变化曲线图;
[0036] 图2示出了姿态角速度估计误差的变化曲线图;
[0037] 图3示出了控制输入变化曲线图;
[0038] 图4示出了舵机堵塞故障估计值变化曲线图。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号