CN110888451A - 一种多旋翼无人机容错控制方法及*** - Google Patents

一种多旋翼无人机容错控制方法及*** Download PDF

Info

Publication number
CN110888451A
CN110888451A CN201911322352.8A CN201911322352A CN110888451A CN 110888451 A CN110888451 A CN 110888451A CN 201911322352 A CN201911322352 A CN 201911322352A CN 110888451 A CN110888451 A CN 110888451A
Authority
CN
China
Prior art keywords
unmanned aerial
aerial vehicle
fault
control
matrix
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
CN201911322352.8A
Other languages
English (en)
Other versions
CN110888451B (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.)
Electric Power Research Institute of State Grid Shanxi Electric Power Co Ltd
Original Assignee
Electric Power Research Institute of State Grid Shanxi Electric Power 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 Electric Power Research Institute of State Grid Shanxi Electric Power Co Ltd filed Critical Electric Power Research Institute of State Grid Shanxi Electric Power Co Ltd
Priority to CN201911322352.8A priority Critical patent/CN110888451B/zh
Publication of CN110888451A publication Critical patent/CN110888451A/zh
Application granted granted Critical
Publication of CN110888451B publication Critical patent/CN110888451B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
    • G05D1/0816Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
    • G05D1/0825Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using mathematical models
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft

Landscapes

  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明涉及一种多旋翼无人机容错控制方法及***,采用改进的线性自抗扰控制的姿态控制算法控制所述无人机的飞行,以保证所述无人机在飞行过程的鲁棒性;当检测到所述无人机的部分电机异常时,构建故障矩阵Ri;基于所述故障矩阵Ri在线建立故障模型;基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息;采用所述改进的线性自抗扰控制的姿态控制算法控制所述故障模型下的无人机的飞行,并按照所述电机的控制分配信息控制所述无人机,以达到所需的姿态和高度。该控制方法及***提高了多旋翼无人机的容错能力,保证了多旋翼无人机具有更大的负载能力和更高的稳定性。

Description

一种多旋翼无人机容错控制方法及***
技术领域
本发明涉及无人机技术领域,尤其涉及一种多旋翼无人机容错控制方法及***。
背景技术
近年来,随着科技的不断进步,无人飞行器,尤其是在四旋翼以上的多旋翼无人飞行器领域取得了快速发展。多旋翼无人机是一种配备了数据处理与传输***、传感器、自动控制***和通信***等机载设备的飞行器,能够进行一定的稳态控制和飞行,而且具备一定的自主飞行能力。当前,多旋翼飞行器现在已广泛应用于农林植保、电力巡检、物流运输等领域,极大程度地方便了人民的生产生活。
当多旋翼无人飞行器发生故障时,飞行状态会发生突变,从而造成不可估计的后果,因此,需要设计出一种容错控制方法来提高多旋翼无人机的容错能力,从而保证多旋翼无人机具有更大的负载能力和更高的稳定性。
发明内容
本发明所要解决的技术问题是提供一种多旋翼无人机容错控制方法及***,提高了多旋翼无人机的容错能力,保证了多旋翼无人机具有更大的负载能力和更高的稳定性。
为解决上述问题,本发明所采取的技术方案是:
一方面,提供了一种多旋翼无人机容错控制方法,其包括:
采用改进的线性自抗扰控制的姿态控制算法控制所述无人机的飞行,以保证所述无人机在飞行过程的鲁棒性;
当检测到所述无人机的部分电机异常时,构建故障矩阵Ri,i为大于等于0小于等于所述无人机中所有电机个数的整数;
基于所述故障矩阵Ri在线建立故障模型;
基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息;
采用所述改进的线性自抗扰控制的姿态控制算法控制所述故障模型下的无人机的飞行,并按照所述电机的控制分配信息控制所述无人机,以达到所需的姿态和高度。
作为本发明的进一步改进,所述改进的线性自抗扰控制的姿态控制算法包括:
安排过渡过程:采用以下公式,通过二阶环节将输入的突变信号转化为缓变信号,然后使输出信号达到期望的输入信号:
Figure BDA0002327484090000011
其中,其中,G(s)代表二阶环节的传递函数,T代表二阶环节的时间常数,s代表代表传递函数中的变量符号;
线性扩张状态观测器:采用以下状态空间方程和公式,实现对模型中各变量进行实时跟踪:
Figure BDA0002327484090000021
其中,x1,x2,x3分别代表所***的状态变量,
Figure BDA0002327484090000022
b0代表估计的控制增益,w代表外部扰动,y代表所述模型的输出,u代表所述模型的输入;
Figure BDA0002327484090000023
其中,z1,z2,z3分别代表所述线性扩张状态观测器的***状态变量,β123分别代表所述线性扩张状态观测器的增益。
作为本发明的进一步改进,所述当检测到所述无人机的部分电机异常时,构建故障矩阵Ri,包括:
实时检测所述无人机的所有电机;
当检测到所述无人机的部分电机异常时,计算故障电机的输出与无故障时的输出的比值,根据所述比值构建故障矩阵Ri
作为本发明的进一步改进,所述基于所述故障矩阵Ri在线建立故障模型,包括:
采用以下公式,在线建立故障模型:
Figure BDA0002327484090000024
其中
Figure BDA0002327484090000025
分别代表大地坐标系下的位置加速度,
Figure BDA0002327484090000026
分别代表在大地坐标系下所述无人机的飞行器姿态角的角加速度,
Figure BDA0002327484090000027
θ,ψ分别代表横滚角、俯仰角和偏航角,Ix,Iy,Iz分别代表所述无人机机身在三个方向的转动惯量,m代表所述无人机的质量,g代表重力加速度,UR,UP,UY,UT分别代表所述无人机的电机均无故障时的横滚力矩、俯仰力矩、偏航力矩以及升力,fp,fq,fr,fz分别表示横滚力矩误差、俯仰力矩误差、偏航力矩误差以及升力误差。
作为本发明的进一步改进,所述基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息,包括:
采用以下公式,得到优化后的分配矩阵Nf,将所述优化后的分配矩阵Nf作为所述无人机上所有电机的控制分配信息:
Nf=Af-1
Nf=AfT(Af·AfT)-1;其中,Af代表部分电机故障后的控制效率矩阵,AfT代表Af的转置;
Af=ARi;其中,A代表故障前的控制效率矩阵;
Figure BDA0002327484090000031
其中,b为升力系数,l为所述无人机的轴距,d为反扭矩系数。
另一方面,提供了一种多旋翼无人机容错控制***,其包括:
第一控制模块,用于采用改进的线性自抗扰控制的姿态控制算法控制所述无人机的飞行,以保证所述无人机在飞行过程的鲁棒性;
故障矩阵构建模块,用于当检测到所述无人机的部分电机异常时,构建故障矩阵Ri,i为大于等于0小于等于所述无人机中所有电机个数的整数;
故障模型建立模块,用于基于所述故障矩阵Ri在线建立故障模型;
分配信息获取模块,用于基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息;
第二控制模块,采用所述改进的线性自抗扰控制的姿态控制算法控制所述故障模型下的无人机的飞行,并按照所述电机的控制分配信息控制所述无人机,以达到所需的姿态和高度。
作为本发明的进一步改进,所述故障矩阵构建模块包括:
安排过渡过程单元,用于采用以下公式,通过二阶环节将输入的突变信号转化为缓变信号,然后使输出信号达到期望的输入信号:
Figure BDA0002327484090000032
其中,其中,G(s)代表二阶环节的传递函数,T代表二阶环节的时间常数,s代表代表传递函数中的变量符号;
线性扩张状态观测器单元,用于采用以下状态空间方程和公式,实现对模型中各变量进行实时跟踪:
Figure BDA0002327484090000041
其中,x1,x2,x3分别代表所***的状态变量,
Figure BDA0002327484090000042
b0代表估计的控制增益,w代表外部扰动,y代表所述模型的输出,u代表所述模型的输入;
Figure BDA0002327484090000043
其中,z1,z2,z3分别代表所述线性扩张状态观测器的***状态变量,β123分别代表所述线性扩张状态观测器的增益。
作为本发明的进一步改进,所述故障模型建立模块包括:
检测单元,用于实时检测所述无人机的所有电机;
故障模型建立单元,用于当检测到所述无人机的部分电机异常时,计算故障电机的输出与无故障时的输出的比值,根据所述比值构建故障矩阵Ri
作为本发明的进一步改进,所述故障模型建立模块包括:
故障模型建立单元,用于采用以下公式,在线建立故障模型:
Figure BDA0002327484090000044
其中
Figure BDA0002327484090000045
分别代表大地坐标系下的位置加速度,
Figure BDA0002327484090000046
分别代表在大地坐标系下所述无人机的飞行器姿态角的角加速度,
Figure BDA0002327484090000047
θ,ψ分别代表横滚角、俯仰角和偏航角,Ix,Iy,Iz分别代表所述无人机机身在三个方向的转动惯量,m代表所述无人机的质量,g代表重力加速度,UR,UP,UY,UT分别代表所述无人机的电机均无故障时的横滚力矩、俯仰力矩、偏航力矩以及升力,fp,fq,fr,fz分别表示横滚力矩误差、俯仰力矩误差、偏航力矩误差以及升力误差。
作为本发明的进一步改进,所述分配信息获取模块包括:
分配信息获取单元,用于采用以下公式,采用以下公式,得到优化后的分配矩阵Nf,将所述优化后的分配矩阵Nf作为所述无人机上所有电机的控制分配信息:
Nf=Af-1
Nf=AfT(Af·AfT)-1;其中,Af代表部分电机故障后的控制效率矩阵,AfT代表Af的转置;
Af=ARi;其中,A代表故障前的控制效率矩阵;
Figure BDA0002327484090000051
其中,b为升力系数,l为所述无人机的轴距,d为反扭矩系数。
采用上述技术方案所产生的有益效果在于:
本发明所提供的一种本发明实施例提供的一种多旋翼无人机容错控制方法及***,采用改进的线性自抗扰控制的姿态控制算法控制所述无人机的飞行,以保证所述无人机在飞行过程的鲁棒性;当检测到所述无人机的部分电机异常时,构建故障矩阵Ri;基于所述故障矩阵Ri在线建立故障模型;基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息;采用所述改进的线性自抗扰控制的姿态控制算法控制所述故障模型下的无人机的飞行,并按照所述电机的控制分配信息控制所述无人机,以达到所需的姿态和高度。在无人机的电机未发生故障时,选用的基本控制律为改进的线性自抗扰控制的姿态控制算法,该控制算法对干扰具有较强的鲁棒能力。而且,当无人机的部分电机发生故障时,还能基于发生故障后的故障矩阵,得到新的飞行模型-故障模型,以及所有电机的控制分配信息,从而保证无人机在该故障模式下飞行的同时,还能减少故障电机的使用,从而使无人机达到平稳飞行状态。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1是本发明实施例提供的一种多旋翼无人机容错控制方法的流程图。
图2是本发明实施例提供的一种发生第一种故障时***未进行分配优化的横滚角的响应曲线图。
图3是本发明实施例提供的一种发生第一种故障时***未进行分配优化的俯仰角的响应曲线图。
图4是本发明实施例提供的一种发生第一种故障时***未进行分配优化的偏航角的响应曲线图。
图5是本发明实施例提供的一种发生第一种故障时***进行分配优化与未发生故障时的横滚角对比图。
图6是本发明实施例提供的一种发生第一种故障时***进行分配优化与未发生故障时的俯仰角对比图。
图7是本发明实施例提供的一种发生第一种故障时***进行分配优化与未发生故障时的偏航角对比图。
图8是本发明实施例提供的一种发生第二种故障时***未进行分配优化的横滚角的响应曲线图。
图9是本发明实施例提供的一种发生第二种故障时***未进行分配优化的俯仰角的响应曲线图。
图10是本发明实施例提供的一种发生第二种故障时***未进行分配优化的偏航角的响应曲线图。
图11是本发明实施例提供的一种发生第二种故障时***进行分配优化与未发生故障时的横滚角对比图。
图12是本发明实施例提供的一种发生第二种故障时***进行分配优化与未发生故障时的俯仰角对比图。
图13是本发明实施例提供的一种发生第二种故障时***进行分配优化与未发生故障时的偏航角对比图。
图14是本发明实施例提供的一种多旋翼无人机容错控制***的结构图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合附图和具体实施例对发明进行清楚、完整的描述。
当前,多旋翼无人机主要分为:四旋翼无人机和四旋翼以上的多旋翼无人机,四旋翼无人机缺少旋翼组件的冗余,一旦发送故障,飞行姿态发生突变,在某些应用领域会造成不可估计的后果。而对于四旋翼以上的多旋翼无人机来说,可以通过优化控制算法,使多旋翼无人机具有很好的容错能力、更大的负载能力和更高的稳定性,以保证无人机在遇到强外力干扰或部分电机受损时仍具有良好的稳定性和安全性,从而可以携带更多的任务设备,完成更复杂的任务。
基于此,本发明涉及的多旋翼无人机指代的是四旋翼以上的多旋翼无人机。
图1和图2是本发明提供的多旋翼无人机容错控制方法的流程图,如图1和图2所示,该容错控制方法包括:
S101:采用改进的线性自抗扰控制的姿态控制算法控制无人机的飞行,以保证无人机在飞行过程的鲁棒性。
其中,改进的线性自抗扰控制的姿态控制算法包括:
(1)安排过渡过程:采用以下公式,通过二阶环节将输入的突变信号转化为缓变信号,然后使输出信号达到期望的输入信号:
Figure BDA0002327484090000071
其中,G(s)代表二阶环节的传递函数,T代表二阶环节的时间常数,T可取期望过渡时间的
Figure BDA0002327484090000072
s代表传递函数中的变量符号,是一个复参数。
(2)线性扩张状态观测器:采用以下状态空间方程和公式,实现对模型中各变量进行实时跟踪:
Figure BDA0002327484090000073
其中,x1,x2,x3分别代表***的状态变量,
Figure BDA0002327484090000074
b0代表估计的控制增益,w代表外部扰动,y代表模型的输出,u代表模型的输入;
Figure BDA0002327484090000075
其中,z1,z2,z3分别代表线性扩张状态观测器的***状态变量,β123分别代表线性扩张状态观测器的增益。
需要说明的是,多旋翼无人机的主要通过改变每个旋翼对应的电机的转速来改变该无人机飞行的姿态,电机的转速是通过改变PWM信号的占空比改变的,电机转速发生变化后,电机产生的扬力和转矩发生变化,同时根据多旋翼无人机自身电机的位置分布,决定了横滚方向、俯仰方向的转矩和偏航方向的反转矩。横滚方向和俯仰方向的变化产生了飞机X,Y轴方向的线速度,偏航角的变化为飞机航向的变化,扬力的变化产生了飞机Z轴即高度的变化。
其中,机体坐标系定义为:原点取在无人机的质心,坐标系与机体固连;X轴与机身设计的纵轴平行,且处于无人机的对称平面内,指向前方;Y轴垂直于无人机对称平面指向右方;Z轴在无人机对称平面内,且垂直于X轴指向下方。整个坐标系符合欧拉坐标系右手定则。地面坐标系即惯性坐标系定义为:采用北东地坐标系,XE轴指向北面,YE轴指向东面,而ZE轴指向地心方向。地面坐标系为仿真实验的环境中的坐标系。
无人机根据用户操作所产生的虚拟控制指令经控制分配信息分配为每个电机的实际控制指令。在四个电机正常的情况下,控制分配信息为固定值,不发生改变。但当无人机的部分电机发生故障时,电机无法正确响应分配的实际控制指令,即电机对于相同的控制指令无法调整为对应的转速,不能满足控制的要求,因此需对控制分配信息进行优化再分配。
S102:当检测到无人机的部分电机异常时,构建故障矩阵Ri,i为大于等于0小于等于无人机中所有电机个数的整数。
其中,该步骤包括:
S1021:实时检测无人机的所有电机。
S1022:当检测到无人机的部分电机异常时,计算故障电机的输出与无故障时的输出的比值,根据比值构建故障矩阵。
关于判断部分电机异常的方式,在一种可能的实现方式中,在无人机飞行过程中,用户可通过推动遥控器的摇杆以向无人机飞行控制***发送控制指令,无人机飞行控制***接收到控制指令后,通过内部的姿态解算环节和姿态控制环节对每个电调输出一定的控制量,之后每个电调输出与该控制量对应的PWM值以控制与该电调连接的电机的转速,即,每个电调输出的PWM值与控制指令相对应。因此,当某一个电调输出的PWM值不与该控制指令相对应时,则确定该电调对应的电机发生故障。
当检测到无人机的部分电机异常时,无人机飞行控制***获取所有电机的输出值,即,电机的效率。其中,本发明实施例中,将未发生故障的电机的效率设为1;而发生故障的电机效率为故障电机的输出相对于正常时的输出的量化值。将每个电机的效率分别以k1,k2…ki表示,则构建得到的故障矩阵Ri为:
Ri=diag[k1,k2…ki]。
例如:该多旋翼无人机为六旋翼无人机,将该六旋翼无人机上的六个电机分别标号为1-6,当1号电机和2号电机由于发生故障导致其电机的效率分别为k1,k2时,该六旋翼无人机的故障矩阵为:
Figure BDA0002327484090000091
S103:基于故障矩阵Ri在线建立故障模型。
采用以下公式,在线建立故障模型:
Figure BDA0002327484090000092
其中,
Figure BDA0002327484090000093
分别代表大地坐标系下的位置加速度,
Figure BDA0002327484090000094
分别代表在大地坐标系下无人机的飞行器姿态角的角加速度,
Figure BDA0002327484090000095
θ,ψ分别代表横滚角、俯仰角和偏航角,Ix,Iy,Iz分别代表无人机机身在三个方向的转动惯量,m代表无人机的质量,g代表重力加速度,UR,UP,UY,UT分别代表无人机的电机均无故障时的横滚力矩、俯仰力矩、偏航力矩以及升力,fp,fq,fr,fz分别表示横滚力矩误差、俯仰力矩误差、偏航力矩误差以及升力误差。
上述
Figure BDA0002327484090000096
可通过传感器得到,传感器包括测量加速度的三轴加速度计和测量角速度的三轴陀螺仪。
fp,fq,fr,fz可由故障矩阵Ri得到,需要说明的是,fp,fq,fr,fz与无人机的旋翼数量及旋翼的布局类型相关,一般情况下,一旦确定无人机的旋翼数量和旋翼的类型,即可得到fp,fq,fr,fz
在一种可能的实现方式中,该多旋无人机为“X”型六旋翼无人机时,将该六旋翼无人机上的六个电机分别标号为1-6,各轴之间夹角为60°,可采用以下方式,得到无人机的电机均无故障时的横滚力矩UR、俯仰力矩UP、偏航力矩UY以及升力UT
Figure BDA0002327484090000101
则当1号电机发生故障后其电机效率为k1,其他电机未发生故障时,故障矩阵R6为:
Figure BDA0002327484090000102
可采用以下方式得到无人机的部分电机出现故障时的的横滚力矩
Figure BDA0002327484090000103
俯仰力矩
Figure BDA0002327484090000104
偏航力矩
Figure BDA0002327484090000105
以及升力
Figure BDA0002327484090000106
Figure BDA0002327484090000107
其中,上述
Figure BDA0002327484090000108
分别代表1-6号电机未出故障时每个电机所对应的转速。
之后,采用以下公式,求得横滚力矩误差fp、俯仰力矩误差fq、偏航力矩误差fr以及升力误差fz
Figure BDA0002327484090000109
S104:基于故障矩阵Ri得到无人机上所有电机的控制分配信息。
采用以下公式,得到优化后的分配矩阵Nf,将优化后的分配矩阵Nf作为无人机上所有电机的控制分配信息:
Nf=Af-1
Nf=AfT(Af·AfT)-1;其中,Af代表部分电机故障后的控制效率矩阵,AfT代表Af的转置;
Af=ARi;其中,A代表故障前的控制效率矩阵;
Figure BDA0002327484090000111
其中,b为升力系数,l为无人机的轴距,d为反扭矩系数。
S105:采用改进的线性自抗扰控制的姿态控制算法控制故障模型下的无人机的飞行,并按照电机的控制分配信息控制无人机,以达到所需的姿态和高度。
根据实际控制指令调整每个电机的转速,进而使每个电机驱动无人机达到所需的姿态和高度。姿态包括俯仰、滚转和偏航。
例如,可采用以下公式调整每个电机的转速:
τf=[UR UP UY UT]T;其中,τf代表力矩矩阵;
Figure BDA0002327484090000112
其中,
Figure BDA0002327484090000113
表示当无人机部分电机异常时,每个电机所对应的转速。
本发明实施例中,无人机飞行控制***根据用户操作分布虚拟控制指令给执行机构,执行机构中的电机的转速发生变化,则无人机飞行控制***在采用改进的线性自抗扰控制的姿态控制算法控制无人机飞行过程中,基于新的控制分配信息对每个电机的控制分配信息进行优化;还会按照建立的故障模型,补偿由于电机故障而产生的横滚力矩误差fp、俯仰力矩误差fq、偏航力矩误差fr以及升力误差fz,从而使得无人机的高度和姿态发生变化,达到无人机所要达到的高度和姿态。
另外,本发明实施例采用六旋翼无人机为对象,对容错控制方法进行了仿真验证,如图2-图13所示,且图2-图13均为给定期望的横滚角、俯仰角和偏航角为15°,给定期望的高度为1m的情况下得到的曲线图。
其中,图2-图4分别为1号电机的效率为1/5时,***未进行分配优化的横滚角、俯仰角和偏航角响应曲线图。由图2-图4可知,当1号电机力效为1/5时且不采用本发明实施例提供的无人机容错控制方法时,三个姿态角无法跟随输入的期望姿态角达到稳定。
图5-图7分别为1号电机的效率为1/5时,***进行分配优化与未发生故障时的横滚角、俯仰角和偏航角的对比图,其中,(a)线表示***进行分配优化的三个姿态角,(b)线代表未发生故障时的三个姿态角。由图5-图7可知,当采用了本发明实施例提供的无人机容错控制方法后,与无故障时的飞行状况相比,虽然没稳定时间较长,但是依然可以很好的跟踪期望值并达到最终稳定。
图8-图10分别为1号电机的效率为1/2时,***未进行分配优化的横滚角、俯仰角和偏航角响应曲线图。由图8-图10可知,当1号电机力效为1/2时且不采用本发明实施例提供的无人机容错控制方法时,三个姿态角无法跟随输入的期望姿态角达到稳定且震荡明显。
图11-图13分别为1号电机的效率为1/2时,***进行分配优化与未发生故障时的横滚角、俯仰角和偏航角的对比图,其中,(a)线表示***进行分配优化的三个姿态角,(b)线代表未发生故障时的三个姿态角。由图11-图13可知,当采用了本发明实施例提供的无人机容错控制方法后,与无故障时的飞行状况相比,虽然没稳定时间较长,但是依然可以很好的跟踪期望值并达到最终稳定。
本发明实施例提供的一种多旋翼无人机容错控制方法,采用改进的线性自抗扰控制的姿态控制算法控制无人机的飞行,以保证无人机在飞行过程的鲁棒性;当检测到无人机的部分电机异常时,构建故障矩阵Ri;基于故障矩阵Ri在线建立故障模型;基于故障矩阵Ri得到无人机上所有电机的控制分配信息;采用改进的线性自抗扰控制的姿态控制算法控制故障模型下的无人机的飞行,并按照电机的控制分配信息控制无人机,以达到所需的姿态和高度。在无人机的电机未发生故障时,选用的基本控制律为改进的线性自抗扰控制的姿态控制算法,该控制算法对干扰具有较强的鲁棒能力。而且,当无人机的部分电机发生故障时,还能基于发生故障后的故障矩阵,得到新的飞行模型-故障模型,以及所有电机的控制分配信息,从而保证无人机在该故障模式下飞行的同时,还能减少故障电机的使用,从而使无人机达到平稳飞行状态。
图14是本发明实施例提供的一种多旋翼无人机容错控制***的结构图,如图14所示,其包括:
第一控制模块1401,用于采用改进的线性自抗扰控制的姿态控制算法控制无人机的飞行,以保证无人机在飞行过程的鲁棒性;
故障矩阵构建模块1402,用于当检测到无人机的部分电机异常时,构建故障矩阵Ri,i为大于等于0小于等于无人机中所有电机个数的整数;
故障模型建立模块1403,用于基于故障矩阵Ri在线建立故障模型;
分配信息获取模块1404,用于基于故障矩阵Ri得到无人机上所有电机的控制分配信息;
第二控制模块1405,采用改进的线性自抗扰控制的姿态控制算法控制故障模型下的无人机的飞行,并按照电机的控制分配信息控制无人机,以达到所需的姿态和高度。
其中,第一控制模块1401和第二控制模块1405可以为同一个控制模块。
在一种可能的实现方式中,故障矩阵构建模块1402包括:
安排过渡过程单元,用于采用以下公式,通过二阶环节将输入的突变信号转化为缓变信号,然后使输出信号达到期望的输入信号:
Figure BDA0002327484090000131
其中,其中,G(s)代表二阶环节的传递函数,T代表二阶环节的时间常数,s代表传递函数中的变量符号;
线性扩张状态观测器单元,用于采用以下状态空间方程和公式,实现对模型中各变量进行实时跟踪:
Figure BDA0002327484090000132
其中,x1,x2,x3分别代表所***的状态变量,
Figure BDA0002327484090000133
b0代表估计的控制增益,w代表外部扰动,y代表模型的输出,u代表模型的输入;
Figure BDA0002327484090000134
其中,z1,z2,z3分别代表线性扩张状态观测器的***状态变量,β123分别代表线性扩张状态观测器的增益。
在一种可能的实现方式中,故障模型建立模块1403包括:
检测单元,用于实时检测无人机的所有电机;
故障模型建立单元,用于当检测到无人机的部分电机异常时,计算故障电机的输出与无故障时的输出的比值,根据比值构建故障矩阵Ri
在一种可能的实现方式中,故障模型建立模块1403包括:
故障模型建立单元,用于采用以下公式,在线建立故障模型:
Figure BDA0002327484090000135
其中
Figure BDA0002327484090000136
分别代表大地坐标系下的位置加速度,
Figure BDA0002327484090000137
分别代表在大地坐标系下无人机的飞行器姿态角的角加速度,
Figure BDA0002327484090000138
θ,ψ分别代表横滚角、俯仰角和偏航角,Ix,Iy,Iz分别代表无人机机身在三个方向的转动惯量,m代表无人机的质量,g代表重力加速度,UR,UP,UY,UT分别代表无人机的电机均无故障时的横滚力矩、俯仰力矩、偏航力矩以及升力,fp,fq,fr,fz分别表示横滚力矩误差、俯仰力矩误差、偏航力矩误差以及升力误差。
在一种可能的实现方式中,分配信息获取模块1404包括:
分配信息获取单元,用于采用以下公式,得到优化后的分配矩阵Nf,将优化后的分配矩阵Nf作为无人机上所有电机的控制分配信息:
Nf=Af-1
Nf=AfT(Af·AfT)-1;其中,Af代表部分电机故障后的控制效率矩阵,AfT代表Af的转置;
Af=ARi;其中,A代表故障前的控制效率矩阵;
Figure BDA0002327484090000141
其中,b为升力系数,l为无人机的轴距,d为反扭矩系数。
本发明实施例提供的一种多旋翼无人机容错控制***,采用改进的线性自抗扰控制的姿态控制算法控制无人机的飞行,以保证无人机在飞行过程的鲁棒性;当检测到无人机的部分电机异常时,构建故障矩阵Ri;基于故障矩阵Ri在线建立故障模型;基于故障矩阵Ri得到无人机上所有电机的控制分配信息;采用改进的线性自抗扰控制的姿态控制算法控制故障模型下的无人机的飞行,并按照电机的控制分配信息控制无人机,以达到所需的姿态和高度。在无人机的电机未发生故障时,选用的基本控制律为改进的线性自抗扰控制的姿态控制算法,该控制算法对干扰具有较强的鲁棒能力。而且,当无人机的部分电机发生故障时,还能基于发生故障后的故障矩阵,得到新的飞行模型-故障模型,以及所有电机的控制分配信息,从而保证无人机在该故障模式下飞行的同时,还能减少故障电机的使用,从而使无人机达到平稳飞行状态。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围。

Claims (10)

1.一种多旋翼无人机容错控制方法,其特征在于,其包括:
采用改进的线性自抗扰控制的姿态控制算法控制所述无人机的飞行,以保证所述无人机在飞行过程的鲁棒性;
当检测到所述无人机的部分电机异常时,构建故障矩阵Ri,i为大于等于0小于等于所述无人机中所有电机个数的整数;
基于所述故障矩阵Ri在线建立故障模型;
基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息;
采用所述改进的线性自抗扰控制的姿态控制算法控制所述故障模型下的无人机的飞行,并按照所述电机的控制分配信息控制所述无人机,以达到所需的姿态和高度。
2.根据权利要求1所述的一种多旋翼无人机容错控制方法,其特征在于,所述改进的线性自抗扰控制的姿态控制算法包括:
安排过渡过程:采用以下公式,通过二阶环节将输入的突变信号转化为缓变信号,然后使输出信号达到期望的输入信号:
Figure FDA0002327484080000011
其中,其中,G(s)代表二阶环节的传递函数,T代表二阶环节的时间常数,s代表传递函数中的变量符号;
线性扩张状态观测器:采用以下状态空间方程和公式,实现对模型中各变量进行实时跟踪:
Figure FDA0002327484080000012
其中,x1,x2,x3分别代表所***的状态变量,
Figure FDA0002327484080000013
b0代表估计的控制增益,w代表外部扰动,y代表所述模型的输出,u代表所述模型的输入;
Figure FDA0002327484080000014
其中,z1,z2,z3分别代表所述线性扩张状态观测器的***状态变量,β123分别代表所述线性扩张状态观测器的增益。
3.根据权利要求1所述的一种多旋翼无人机容错控制方法,其特征在于,所述当检测到所述无人机的部分电机异常时,构建故障矩阵Ri,包括:
实时检测所述无人机的所有电机;
当检测到所述无人机的部分电机异常时,计算故障电机的输出与无故障时的输出的比值,根据所述比值构建故障矩阵Ri
4.根据权利要求1所述的一种多旋翼无人机容错控制方法,其特征在于,所述基于所述故障矩阵Ri在线建立故障模型,包括:
采用以下公式,在线建立故障模型:
Figure FDA0002327484080000021
其中
Figure FDA0002327484080000022
分别代表大地坐标系下的位置加速度,
Figure FDA0002327484080000023
分别代表在大地坐标系下所述无人机的飞行器姿态角的角加速度,
Figure FDA0002327484080000024
θ,ψ分别代表横滚角、俯仰角和偏航角,Ix,Iy,Iz分别代表所述无人机机身在三个方向的转动惯量,m代表所述无人机的质量,g代表重力加速度,UR,UP,UY,UT分别代表所述无人机的电机均无故障时的横滚力矩、俯仰力矩、偏航力矩以及升力,fp,fq,fr,fz分别表示横滚力矩误差、俯仰力矩误差、偏航力矩误差以及升力误差。
5.根据权利要求1所述的一种多旋翼无人机容错控制方法,其特征在于,所述基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息,包括:
采用以下公式,得到优化后的分配矩阵Nf,将所述优化后的分配矩阵Nf作为所述无人机上所有电机的控制分配信息:
Nf=Af-1
Nf=AfT(Af·AfT)-1;其中,Af代表部分电机故障后的控制效率矩阵,AfT代表Af的转置;
Af=ARi;其中,A代表故障前的控制效率矩阵;
Figure FDA0002327484080000025
其中,b为升力系数,l为所述无人机的轴距,d为反扭矩系数。
6.一种多旋翼无人机容错控制***,其特征在于,其包括:
第一控制模块,用于采用改进的线性自抗扰控制的姿态控制算法控制所述无人机的飞行,以保证所述无人机在飞行过程的鲁棒性;
故障矩阵构建模块,用于当检测到所述无人机的部分电机异常时,构建故障矩阵Ri,i为大于等于0小于等于所述无人机中所有电机个数的整数;
故障模型建立模块,用于基于所述故障矩阵Ri在线建立故障模型;
分配信息获取模块,用于基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息;
第二控制模块,采用所述改进的线性自抗扰控制的姿态控制算法控制所述故障模型下的无人机的飞行,并按照所述电机的控制分配信息控制所述无人机,以达到所需的姿态和高度。
7.根据权利要求6所述的一种多旋翼无人机容错控制***,其特征在于,所述故障矩阵构建模块包括:
安排过渡过程,用于采用以下公式,通过二阶环节将输入的突变信号转化为缓变信号,然后使输出信号达到期望的输入信号:
Figure FDA0002327484080000031
其中,其中,G(s)代表二阶环节的传递函数,T代表二阶环节的时间常数,s代表传递函数中的变量符号;
线性扩张状态观测器单元,用于采用以下状态空间方程和公式,实现对模型中各变量进行实时跟踪:
Figure FDA0002327484080000032
其中,x1,x2,x3分别代表所***的状态变量,
Figure FDA0002327484080000033
b0代表估计的控制增益,w代表外部扰动,y代表所述模型的输出,u代表所述模型的输入;
Figure FDA0002327484080000034
其中,z1,z2,z3分别代表所述线性扩张状态观测器的***状态变量,β123分别代表所述线性扩张状态观测器的增益。
8.根据权利要求6所述的一种多旋翼无人机容错控制***,其特征在于,所述故障模型建立模块包括:
检测单元,用于实时检测所述无人机的所有电机;
故障模型建立单元,用于当检测到所述无人机的部分电机异常时,计算故障电机的输出与无故障时的输出的比值,根据所述比值构建故障矩阵Ri
9.根据权利要求6所述的一种多旋翼无人机容错控制***,其特征在于,所述故障模型建立模块包括:
故障模型建立单元,用于采用以下公式,在线建立故障模型:
Figure FDA0002327484080000041
其中
Figure FDA0002327484080000042
分别代表大地坐标系下的位置加速度,
Figure FDA0002327484080000043
分别代表在大地坐标系下所述无人机的飞行器姿态角的角加速度,
Figure FDA0002327484080000044
θ,ψ分别代表横滚角、俯仰角和偏航角,Ix,Iy,Iz分别代表所述无人机机身在三个方向的转动惯量,m代表所述无人机的质量,g代表重力加速度,UR,UP,UY,UT分别代表所述无人机的电机均无故障时的横滚力矩、俯仰力矩、偏航力矩以及升力,fp,fq,fr,fz分别表示横滚力矩误差、俯仰力矩误差、偏航力矩误差以及升力误差。
10.根据权利要求6所述的一种多旋翼无人机容错控制***,其特征在于,所述分配信息获取模块包括:
分配信息获取单元,用于采用以下公式,得到优化后的分配矩阵Nf,将所述优化后的分配矩阵Nf作为所述无人机上所有电机的控制分配信息:
Nf=Af-1
Nf=AfT(Af·AfT)-1;其中,Af代表部分电机故障后的控制效率矩阵,AfT代表Af的转置;
Af=ARi;其中,A代表故障前的控制效率矩阵;
Figure FDA0002327484080000051
其中,b为升力系数,l为所述无人机的轴距,d为反扭矩系数。
CN201911322352.8A 2019-12-20 2019-12-20 一种多旋翼无人机容错控制方法及*** Active CN110888451B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911322352.8A CN110888451B (zh) 2019-12-20 2019-12-20 一种多旋翼无人机容错控制方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911322352.8A CN110888451B (zh) 2019-12-20 2019-12-20 一种多旋翼无人机容错控制方法及***

Publications (2)

Publication Number Publication Date
CN110888451A true CN110888451A (zh) 2020-03-17
CN110888451B CN110888451B (zh) 2022-12-06

Family

ID=69752610

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911322352.8A Active CN110888451B (zh) 2019-12-20 2019-12-20 一种多旋翼无人机容错控制方法及***

Country Status (1)

Country Link
CN (1) CN110888451B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111290424A (zh) * 2020-03-26 2020-06-16 南方医科大学南方医院 用于医院血液样本运输的无人机姿态控制方法及无人机
CN112327896A (zh) * 2020-10-29 2021-02-05 东北大学 旋翼容错控制方法、装置、计算机存储介质及计算机设备
CN112373677A (zh) * 2020-11-18 2021-02-19 三生万物(北京)人工智能技术有限公司 一种六旋翼无人机动力缺失保护***、保护方法
CN112947532A (zh) * 2021-04-13 2021-06-11 中国计量大学 一种基于药液测量的植保四旋翼无人机高度控制方法
CN113039502A (zh) * 2020-05-07 2021-06-25 深圳市大疆创新科技有限公司 多旋翼无人机及其控制方法、控制装置和计算机可读存储介质
CN113467244A (zh) * 2021-07-15 2021-10-01 之江实验室 基于控制矩阵重分配的多旋翼无人机容错控制方法及装置
WO2021223136A1 (zh) * 2020-05-07 2021-11-11 深圳市大疆创新科技有限公司 控制方法、数据处理方法、设备、无人机及存储介质
WO2022193075A1 (zh) * 2021-03-15 2022-09-22 深圳市大疆创新科技有限公司 无人飞行器控制方法、无人飞行器和存储介质
CN115407802A (zh) * 2022-10-26 2022-11-29 广东汇天航空航天科技有限公司 多旋翼飞行器及其控制方法、设备、计算机可读存储介质
CN115857309A (zh) * 2023-02-27 2023-03-28 华东交通大学 一种可倾斜六旋翼无人机的容错控制方法及***
CN115877717A (zh) * 2022-12-26 2023-03-31 南通大学 一种基于自抗扰控制的飞行器容错控制结构与控制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102902277A (zh) * 2012-10-19 2013-01-30 北京航空航天大学 一种六旋翼飞行器有一个电机停机时的安全着陆方法
CN106444812A (zh) * 2016-10-26 2017-02-22 华南智能机器人创新研究院 一种基于四旋翼无人机的姿态控制的方法及其***
CN108614573A (zh) * 2018-05-15 2018-10-02 上海扩博智能技术有限公司 六旋翼无人机的自动容错姿态控制方法
CN109308064A (zh) * 2017-07-28 2019-02-05 深圳禾苗通信科技有限公司 一种四旋翼无人机的故障容错控制方法及***

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102902277A (zh) * 2012-10-19 2013-01-30 北京航空航天大学 一种六旋翼飞行器有一个电机停机时的安全着陆方法
CN106444812A (zh) * 2016-10-26 2017-02-22 华南智能机器人创新研究院 一种基于四旋翼无人机的姿态控制的方法及其***
CN109308064A (zh) * 2017-07-28 2019-02-05 深圳禾苗通信科技有限公司 一种四旋翼无人机的故障容错控制方法及***
CN108614573A (zh) * 2018-05-15 2018-10-02 上海扩博智能技术有限公司 六旋翼无人机的自动容错姿态控制方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JIMOH O. PEDRO: "PI-Based Fault Tolerant Control For Fixed-Wing UAVs Using Control Allocation", 《IFAC-PAPERSONLINE》 *
SONGMING JIAO: "Fault tolerant control Algorithm of hexarotor UAV", 《JOURNAL OF ROBOTICS》 *
刘登攀: "六旋翼无人机故障容错及精准定高策略研究", 《工程科技Ⅱ辑》 *
袁东 等: "二阶***线性自抗扰控制器频带特性与参数配置研究", 《控制理论与应用》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111290424A (zh) * 2020-03-26 2020-06-16 南方医科大学南方医院 用于医院血液样本运输的无人机姿态控制方法及无人机
WO2021223136A1 (zh) * 2020-05-07 2021-11-11 深圳市大疆创新科技有限公司 控制方法、数据处理方法、设备、无人机及存储介质
CN113039502B (zh) * 2020-05-07 2024-04-12 深圳市大疆创新科技有限公司 多旋翼无人机及其控制方法、控制装置和计算机可读存储介质
CN113039502A (zh) * 2020-05-07 2021-06-25 深圳市大疆创新科技有限公司 多旋翼无人机及其控制方法、控制装置和计算机可读存储介质
CN112327896A (zh) * 2020-10-29 2021-02-05 东北大学 旋翼容错控制方法、装置、计算机存储介质及计算机设备
CN112373677A (zh) * 2020-11-18 2021-02-19 三生万物(北京)人工智能技术有限公司 一种六旋翼无人机动力缺失保护***、保护方法
WO2022193075A1 (zh) * 2021-03-15 2022-09-22 深圳市大疆创新科技有限公司 无人飞行器控制方法、无人飞行器和存储介质
CN112947532B (zh) * 2021-04-13 2022-05-03 中国计量大学 一种基于药液测量的植保四旋翼无人机高度控制方法
CN112947532A (zh) * 2021-04-13 2021-06-11 中国计量大学 一种基于药液测量的植保四旋翼无人机高度控制方法
CN113467244A (zh) * 2021-07-15 2021-10-01 之江实验室 基于控制矩阵重分配的多旋翼无人机容错控制方法及装置
CN113467244B (zh) * 2021-07-15 2024-03-26 之江实验室 基于控制矩阵重分配的多旋翼无人机容错控制方法及装置
CN115407802A (zh) * 2022-10-26 2022-11-29 广东汇天航空航天科技有限公司 多旋翼飞行器及其控制方法、设备、计算机可读存储介质
WO2024087271A1 (zh) * 2022-10-26 2024-05-02 广东汇天航空航天科技有限公司 多旋翼飞行器及其控制方法、设备、计算机可读存储介质
CN115877717A (zh) * 2022-12-26 2023-03-31 南通大学 一种基于自抗扰控制的飞行器容错控制结构与控制方法
CN115877717B (zh) * 2022-12-26 2023-06-13 南通大学 一种基于自抗扰控制的飞行器容错控制结构与控制方法
CN115857309A (zh) * 2023-02-27 2023-03-28 华东交通大学 一种可倾斜六旋翼无人机的容错控制方法及***

Also Published As

Publication number Publication date
CN110888451B (zh) 2022-12-06

Similar Documents

Publication Publication Date Title
CN110888451B (zh) 一种多旋翼无人机容错控制方法及***
Song et al. Neuroadaptive fault-tolerant control of quadrotor UAVs: A more affordable solution
CN109991991B (zh) 一种无人直升机鲁棒容错跟踪方法
Chen et al. Robust backstepping sliding-mode control and observer-based fault estimation for a quadrotor UAV
EP3424820B1 (en) Controlled flight of a multicopter experiencing a failure affecting an effector
Ranjbaran et al. Fault recovery of an under-actuated quadrotor aerial vehicle
CN111766899A (zh) 一种基于干扰观测器的四旋翼无人机集群抗干扰编队控制方法
US20220036739A1 (en) Method for Controlling Steady Flight of Unmanned Aircraft
CN106647781A (zh) 基于重复控制补偿神经模糊pid四旋翼飞行器的控制方法
US20210107626A1 (en) Method of controlling an actuator system and aircraft using same
CN111880410B (zh) 一种针对电机故障的四旋翼无人机容错控制方法
CN111273688A (zh) 基于事件触发的四旋翼无人机一致性编队控制方法
CN112684705B (zh) 一种四旋翼飞行器编队追踪控制方法
CN115657730B (zh) 一种大规模多旋翼无人机鲁棒分簇编队控制方法
Sun et al. Nonlinear robust compensation method for trajectory tracking control of quadrotors
CN112327896A (zh) 旋翼容错控制方法、装置、计算机存储介质及计算机设备
Sahul et al. A novel method on Disturbance Rejection PID Controller for Quadcopter based on Optimization algorithm
Khebbache et al. Robust stabilization of a quadrotor UAV in presence of actuator and sensor faults
Khebbache Robust control algorithm considering the actuator faults for attitude tracking of an uav quadrotor aircraft
CN114326819A (zh) 基于耦合力场无人机建模与结构损伤自适应容错控制方法
CN116848485A (zh) 无人飞行器控制方法、无人飞行器和存储介质
CN116009606A (zh) 无人机的动力分配方法、装置、无人机及存储介质
Heidarian et al. Attitude control of VTOL-UAVs
Li et al. Leader-follower formation control of mutilple UAVs based on ADRC: experiment research
CN116414148B (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