[go: up one dir, main page]

CN110888451B - A fault-tolerant control method and system for a multi-rotor UAV - Google Patents

A fault-tolerant control method and system for a multi-rotor UAV Download PDF

Info

Publication number
CN110888451B
CN110888451B CN201911322352.8A CN201911322352A CN110888451B CN 110888451 B CN110888451 B CN 110888451B CN 201911322352 A CN201911322352 A CN 201911322352A CN 110888451 B CN110888451 B CN 110888451B
Authority
CN
China
Prior art keywords
fault
control
uav
matrix
motors
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
Application number
CN201911322352.8A
Other languages
Chinese (zh)
Other versions
CN110888451A (en
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.)
State Grid Electric Power Research Institute Of Sepc
Original Assignee
State Grid Electric Power Research Institute Of Sepc
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 State Grid Electric Power Research Institute Of Sepc filed Critical State Grid Electric Power Research Institute Of Sepc
Priority to CN201911322352.8A priority Critical patent/CN110888451B/en
Publication of CN110888451A publication Critical patent/CN110888451A/en
Application granted granted Critical
Publication of CN110888451B publication Critical patent/CN110888451B/en
Expired - Fee Related 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

The invention relates to a fault-tolerant control method and a fault-tolerant control system for a multi-rotor unmanned aerial vehicle, wherein the flight of the unmanned aerial vehicle is controlled by adopting an improved attitude control algorithm of linear active disturbance rejection control so as to ensure the robustness of the unmanned aerial vehicle in the flight process; when detecting that partial motor of unmanned aerial vehicle is unusual, construct trouble matrix R i (ii) a Based on the fault matrix R i Establishing a fault model on line; based on the fault matrix R i Obtaining control distribution information of all motors on the unmanned aerial vehicle; and controlling the flight of the unmanned aerial vehicle under the fault model by adopting the attitude control algorithm of the improved linear active disturbance rejection control, and controlling the unmanned aerial vehicle according to the control distribution information of the motor so as to achieve the required attitude and height. The control method and the control system improve the fault-tolerant capability of the multi-rotor unmanned aerial vehicle, and ensure that the multi-rotor unmanned aerial vehicle has larger load capacity and higher stability.

Description

一种多旋翼无人机容错控制方法及系统A fault-tolerant control method and system for a multi-rotor UAV

技术领域technical field

本发明涉及无人机技术领域,尤其涉及一种多旋翼无人机容错控制方法及系统。The invention relates to the technical field of unmanned aerial vehicles, in particular to a method and system for fault-tolerant control of a multi-rotor unmanned aerial vehicle.

背景技术Background technique

近年来,随着科技的不断进步,无人飞行器,尤其是在四旋翼以上的多旋翼无人飞行器领域取得了快速发展。多旋翼无人机是一种配备了数据处理与传输系统、传感器、自动控制系统和通信系统等机载设备的飞行器,能够进行一定的稳态控制和飞行,而且具备一定的自主飞行能力。当前,多旋翼飞行器现在已广泛应用于农林植保、电力巡检、物流运输等领域,极大程度地方便了人民的生产生活。In recent years, with the continuous advancement of science and technology, unmanned aerial vehicles, especially in the field of multi-rotor unmanned aerial vehicles with more than four rotors, have achieved rapid development. A multi-rotor UAV is an aircraft equipped with airborne equipment such as data processing and transmission systems, sensors, automatic control systems, and communication systems. It can perform certain steady-state control and flight, and has certain autonomous flight capabilities. At present, multi-rotor aircraft has been widely used in fields such as agriculture, forestry and plant protection, power inspection, logistics and transportation, which greatly facilitates people's production and life.

当多旋翼无人飞行器发生故障时,飞行状态会发生突变,从而造成不可估计的后果,因此,需要设计出一种容错控制方法来提高多旋翼无人机的容错能力,从而保证多旋翼无人机具有更大的负载能力和更高的稳定性。When the multi-rotor unmanned aerial vehicle fails, the flight state will change suddenly, resulting in unpredictable consequences. Therefore, it is necessary to design a fault-tolerant control method to improve the fault-tolerant ability of the multi-rotor UAV, so as to ensure that the multi-rotor unmanned The machine has greater load capacity and higher stability.

发明内容Contents of the invention

本发明所要解决的技术问题是提供一种多旋翼无人机容错控制方法及系统,提高了多旋翼无人机的容错能力,保证了多旋翼无人机具有更大的负载能力和更高的稳定性。The technical problem to be solved by the present invention is to provide a multi-rotor UAV fault-tolerant control method and system, which improves the fault-tolerant capability of the multi-rotor UAV and ensures that the multi-rotor UAV has a larger load capacity and a higher load capacity. stability.

为解决上述问题,本发明所采取的技术方案是:In order to solve the problems referred to above, the technical scheme that the present invention takes is:

一方面,提供了一种多旋翼无人机容错控制方法,其包括:On the one hand, a fault-tolerant control method for a multi-rotor UAV is provided, which includes:

采用改进的线性自抗扰控制的姿态控制算法控制所述无人机的飞行,以保证所述无人机在飞行过程的鲁棒性;The attitude control algorithm of the improved linear active disturbance rejection control is used to control the flight of the UAV, so as to ensure the robustness of the UAV during flight;

当检测到所述无人机的部分电机异常时,构建故障矩阵Ri,i为大于等于0小于等于所述无人机中所有电机个数的整数;When some motors of the UAV are detected to be abnormal, a fault matrix R i is constructed, where i is an integer greater than or equal to 0 and less than or equal to the number of all motors in the UAV;

基于所述故障矩阵Ri在线建立故障模型;Establishing a fault model online based on the fault matrix R i ;

基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息;Obtain the control allocation information of all motors on the UAV based on the fault matrix R i ;

采用所述改进的线性自抗扰控制的姿态控制算法控制所述故障模型下的无人机的飞行,并按照所述电机的控制分配信息控制所述无人机,以达到所需的姿态和高度。Adopt the attitude control algorithm of the improved linear ADRC to control the flight of the unmanned aerial vehicle under the fault model, and control the unmanned aerial vehicle according to the control distribution information of the motor to achieve the required attitude and high.

作为本发明的进一步改进,所述改进的线性自抗扰控制的姿态控制算法包括:As a further improvement of the present invention, the attitude control algorithm of the improved linear active disturbance rejection control includes:

安排过渡过程:采用以下公式,通过二阶环节将输入的突变信号转化为缓变信号,然后使输出信号达到期望的输入信号:Arrange the transition process: use the following formula to convert the input sudden change signal into a slow change signal through the second-order link, and then make the output signal reach the desired input signal:

Figure BDA0002327484090000011
其中,其中,G(s)代表二阶环节的传递函数,T代表二阶环节的时间常数,s代表代表传递函数中的变量符号;
Figure BDA0002327484090000011
Wherein, G(s) represents the transfer function of the second-order link, T represents the time constant of the second-order link, and s represents the variable symbol in the transfer function;

线性扩张状态观测器:采用以下状态空间方程和公式,实现对模型中各变量进行实时跟踪:Linear expansion state observer: use the following state space equations and formulas to realize real-time tracking of variables in the model:

Figure BDA0002327484090000021
其中,x1,x2,x3分别代表所系统的状态变量,
Figure BDA0002327484090000022
b0代表估计的控制增益,w代表外部扰动,y代表所述模型的输出,u代表所述模型的输入;
Figure BDA0002327484090000021
Among them, x 1 , x 2 , x 3 respectively represent the state variables of the system,
Figure BDA0002327484090000022
b 0 represents the estimated control gain, w represents the external disturbance, y represents the output of the model, and u represents the input of the model;

Figure BDA0002327484090000023
其中,z1,z2,z3分别代表所述线性扩张状态观测器的系统状态变量,β123分别代表所述线性扩张状态观测器的增益。
Figure BDA0002327484090000023
Wherein, z 1 , z 2 , z 3 respectively represent the system state variables of the linearly extended state observer, and β 1 , β 2 , β 3 represent the gain of the linearly extended state observer respectively.

作为本发明的进一步改进,所述当检测到所述无人机的部分电机异常时,构建故障矩阵Ri,包括:As a further improvement of the present invention, when the part of the motor of the drone is detected to be abnormal, the fault matrix R i is constructed, including:

实时检测所述无人机的所有电机;Real-time detection of all motors of the drone;

当检测到所述无人机的部分电机异常时,计算故障电机的输出与无故障时的输出的比值,根据所述比值构建故障矩阵RiWhen some motors of the UAV are detected to be abnormal, the ratio of the output of the faulty motor to the output without fault is calculated, and a fault matrix R i is constructed according to the ratio.

作为本发明的进一步改进,所述基于所述故障矩阵Ri在线建立故障模型,包括:As a further improvement of the present invention, the online establishment of a fault model based on the fault matrix R i includes:

采用以下公式,在线建立故障模型:The fault model is established online using the following formula:

Figure BDA0002327484090000024
Figure BDA0002327484090000024

其中

Figure BDA0002327484090000025
分别代表大地坐标系下的位置加速度,
Figure BDA0002327484090000026
分别代表在大地坐标系下所述无人机的飞行器姿态角的角加速度,
Figure BDA0002327484090000027
θ,ψ分别代表横滚角、俯仰角和偏航角,Ix,Iy,Iz分别代表所述无人机机身在三个方向的转动惯量,m代表所述无人机的质量,g代表重力加速度,UR,UP,UY,UT分别代表所述无人机的电机均无故障时的横滚力矩、俯仰力矩、偏航力矩以及升力,fp,fq,fr,fz分别表示横滚力矩误差、俯仰力矩误差、偏航力矩误差以及升力误差。in
Figure BDA0002327484090000025
Represent the position acceleration in the earth coordinate system,
Figure BDA0002327484090000026
represent the angular acceleration of the aircraft attitude angle of the unmanned aerial vehicle under the earth coordinate system respectively,
Figure BDA0002327484090000027
θ, ψ represent the roll angle, pitch angle and yaw angle respectively, I x , I y , I z represent the moments of inertia of the UAV fuselage in three directions respectively, and m represents the mass of the UAV , g represents the acceleration of gravity, U R , U P , U Y , U T respectively represent the rolling moment, pitching moment, yaw moment and lift when the motors of the UAV are not faulty, f p , f q , f r , f z represent roll moment error, pitch moment error, yaw moment error and lift error respectively.

作为本发明的进一步改进,所述基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息,包括:As a further improvement of the present invention, the control distribution information of all motors on the drone is obtained based on the fault matrix R i , including:

采用以下公式,得到优化后的分配矩阵Nf,将所述优化后的分配矩阵Nf作为所述无人机上所有电机的控制分配信息:The following formula is used to obtain the optimized allocation matrix N f , and the optimized allocation matrix N f is used as the control allocation information of all motors on the drone:

Nf=Af-1Nf = Af -1 ;

Nf=AfT(Af·AfT)-1;其中,Af代表部分电机故障后的控制效率矩阵,AfT代表Af的转置;N f = A fT (A f A fT ) -1 ; wherein, A f represents the control efficiency matrix after a partial motor failure, and A fT represents the transposition of A f ;

Af=ARi;其中,A代表故障前的控制效率矩阵;A f = AR i ; where, A represents the control efficiency matrix before failure;

Figure BDA0002327484090000031
其中,b为升力系数,l为所述无人机的轴距,d为反扭矩系数。
Figure BDA0002327484090000031
Wherein, b is the lift coefficient, l is the wheelbase of the UAV, and d is the counter torque coefficient.

另一方面,提供了一种多旋翼无人机容错控制系统,其包括:On the other hand, a fault-tolerant control system for a multi-rotor UAV is provided, which includes:

第一控制模块,用于采用改进的线性自抗扰控制的姿态控制算法控制所述无人机的飞行,以保证所述无人机在飞行过程的鲁棒性;The first control module is used to control the flight of the unmanned aerial vehicle by adopting the attitude control algorithm of the improved linear active disturbance rejection control, so as to ensure the robustness of the unmanned aerial vehicle during flight;

故障矩阵构建模块,用于当检测到所述无人机的部分电机异常时,构建故障矩阵Ri,i为大于等于0小于等于所述无人机中所有电机个数的整数;A fault matrix construction module, used to construct a fault matrix R i when detecting abnormality of some motors of the drone, where i is an integer greater than or equal to 0 and less than or equal to the number of all motors in the drone;

故障模型建立模块,用于基于所述故障矩阵Ri在线建立故障模型;A fault model building module, configured to build a fault model online based on the fault matrix R i ;

分配信息获取模块,用于基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息;An allocation information acquisition module, configured to obtain control allocation information of all motors on the UAV based on the fault matrix R i ;

第二控制模块,采用所述改进的线性自抗扰控制的姿态控制算法控制所述故障模型下的无人机的飞行,并按照所述电机的控制分配信息控制所述无人机,以达到所需的姿态和高度。The second control module adopts the attitude control algorithm of the improved linear active disturbance rejection control to control the flight of the UAV under the fault model, and controls the UAV according to the control distribution information of the motor, so as to achieve Desired attitude and height.

作为本发明的进一步改进,所述故障矩阵构建模块包括:As a further improvement of the present invention, the fault matrix building block includes:

安排过渡过程单元,用于采用以下公式,通过二阶环节将输入的突变信号转化为缓变信号,然后使输出信号达到期望的输入信号:The transition process unit is arranged to convert the input sudden change signal into a slow change signal through the second-order link by using the following formula, and then make the output signal reach the desired input signal:

Figure BDA0002327484090000032
其中,其中,G(s)代表二阶环节的传递函数,T代表二阶环节的时间常数,s代表代表传递函数中的变量符号;
Figure BDA0002327484090000032
Wherein, G(s) represents the transfer function of the second-order link, T represents the time constant of the second-order link, and s represents the variable symbol in the transfer function;

线性扩张状态观测器单元,用于采用以下状态空间方程和公式,实现对模型中各变量进行实时跟踪:The Linear Extended State Observer unit is used to implement real-time tracking of variables in the model using the following state space equations and formulas:

Figure BDA0002327484090000041
其中,x1,x2,x3分别代表所系统的状态变量,
Figure BDA0002327484090000042
b0代表估计的控制增益,w代表外部扰动,y代表所述模型的输出,u代表所述模型的输入;
Figure BDA0002327484090000041
Among them, x 1 , x 2 , x 3 respectively represent the state variables of the system,
Figure BDA0002327484090000042
b 0 represents the estimated control gain, w represents the external disturbance, y represents the output of the model, and u represents the input of the model;

Figure BDA0002327484090000043
其中,z1,z2,z3分别代表所述线性扩张状态观测器的系统状态变量,β123分别代表所述线性扩张状态观测器的增益。
Figure BDA0002327484090000043
Wherein, z 1 , z 2 , z 3 respectively represent the system state variables of the linearly extended state observer, and β 1 , β 2 , β 3 represent the gain of the linearly extended state observer respectively.

作为本发明的进一步改进,所述故障模型建立模块包括:As a further improvement of the present invention, the fault model building module includes:

检测单元,用于实时检测所述无人机的所有电机;A detection unit is used to detect all motors of the drone in real time;

故障模型建立单元,用于当检测到所述无人机的部分电机异常时,计算故障电机的输出与无故障时的输出的比值,根据所述比值构建故障矩阵RiThe fault model building unit is used to calculate the ratio of the output of the faulty motor to the output when there is no fault when some motors of the drone are detected to be abnormal, and construct a fault matrix R i according to the ratio.

作为本发明的进一步改进,所述故障模型建立模块包括:As a further improvement of the present invention, the fault model building module includes:

故障模型建立单元,用于采用以下公式,在线建立故障模型:The fault model establishing unit is used to establish the fault model online by adopting the following formula:

Figure BDA0002327484090000044
Figure BDA0002327484090000044

其中

Figure BDA0002327484090000045
分别代表大地坐标系下的位置加速度,
Figure BDA0002327484090000046
分别代表在大地坐标系下所述无人机的飞行器姿态角的角加速度,
Figure BDA0002327484090000047
θ,ψ分别代表横滚角、俯仰角和偏航角,Ix,Iy,Iz分别代表所述无人机机身在三个方向的转动惯量,m代表所述无人机的质量,g代表重力加速度,UR,UP,UY,UT分别代表所述无人机的电机均无故障时的横滚力矩、俯仰力矩、偏航力矩以及升力,fp,fq,fr,fz分别表示横滚力矩误差、俯仰力矩误差、偏航力矩误差以及升力误差。in
Figure BDA0002327484090000045
Represent the position acceleration in the earth coordinate system,
Figure BDA0002327484090000046
represent the angular acceleration of the aircraft attitude angle of the unmanned aerial vehicle under the earth coordinate system respectively,
Figure BDA0002327484090000047
θ, ψ represent the roll angle, pitch angle and yaw angle respectively, I x , I y , I z represent the moments of inertia of the UAV fuselage in three directions respectively, and m represents the mass of the UAV , g represents the acceleration of gravity, U R , U P , U Y , U T respectively represent the rolling moment, pitching moment, yaw moment and lift when the motors of the UAV are not faulty, f p , f q , f r , f z represent roll moment error, pitch moment error, yaw moment error and lift error respectively.

作为本发明的进一步改进,所述分配信息获取模块包括:As a further improvement of the present invention, the allocation information acquisition module includes:

分配信息获取单元,用于采用以下公式,采用以下公式,得到优化后的分配矩阵Nf,将所述优化后的分配矩阵Nf作为所述无人机上所有电机的控制分配信息:The allocation information acquisition unit is used to adopt the following formula to obtain the optimized allocation matrix Nf , and use the optimized allocation matrix Nf as the control allocation information of all motors on the drone:

Nf=Af-1Nf = Af -1 ;

Nf=AfT(Af·AfT)-1;其中,Af代表部分电机故障后的控制效率矩阵,AfT代表Af的转置;N f = A fT (A f A fT ) -1 ; wherein, A f represents the control efficiency matrix after a partial motor failure, and A fT represents the transposition of A f ;

Af=ARi;其中,A代表故障前的控制效率矩阵;A f = AR i ; where, A represents the control efficiency matrix before failure;

Figure BDA0002327484090000051
其中,b为升力系数,l为所述无人机的轴距,d为反扭矩系数。
Figure BDA0002327484090000051
Wherein, b is the lift coefficient, l is the wheelbase of the UAV, and d is the counter torque coefficient.

采用上述技术方案所产生的有益效果在于:The beneficial effects produced by adopting the above-mentioned technical scheme are:

本发明所提供的一种本发明实施例提供的一种多旋翼无人机容错控制方法及系统,采用改进的线性自抗扰控制的姿态控制算法控制所述无人机的飞行,以保证所述无人机在飞行过程的鲁棒性;当检测到所述无人机的部分电机异常时,构建故障矩阵Ri;基于所述故障矩阵Ri在线建立故障模型;基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息;采用所述改进的线性自抗扰控制的姿态控制算法控制所述故障模型下的无人机的飞行,并按照所述电机的控制分配信息控制所述无人机,以达到所需的姿态和高度。在无人机的电机未发生故障时,选用的基本控制律为改进的线性自抗扰控制的姿态控制算法,该控制算法对干扰具有较强的鲁棒能力。而且,当无人机的部分电机发生故障时,还能基于发生故障后的故障矩阵,得到新的飞行模型-故障模型,以及所有电机的控制分配信息,从而保证无人机在该故障模式下飞行的同时,还能减少故障电机的使用,从而使无人机达到平稳飞行状态。The present invention provides a multi-rotor unmanned aerial vehicle fault-tolerant control method and system provided by an embodiment of the present invention. The attitude control algorithm of the improved linear active disturbance rejection control is used to control the flight of the unmanned aerial vehicle to ensure that all Describe the robustness of the unmanned aerial vehicle in the flight process; when detecting the abnormal part of the motor of the unmanned aerial vehicle, construct the fault matrix R i ; establish the fault model online based on the fault matrix R i ; based on the fault matrix R i obtain the control distribution information of all motors on the UAV; adopt the attitude control algorithm of the improved linear active disturbance rejection control to control the flight of the UAV under the fault model, and distribute information according to the control distribution of the motors Control the drone to achieve the desired attitude and altitude. When the motor of the UAV does not fail, the basic control law selected is the attitude control algorithm of the improved linear active disturbance rejection control, which has strong robustness to disturbance. Moreover, when some motors of the UAV fail, based on the failure matrix after the failure, a new flight model-fault model and the control distribution information of all motors can be obtained, so as to ensure that the UAV can operate in this failure mode. While flying, it can also reduce the use of faulty motors, so that the UAV can achieve a stable flight state.

附图说明Description of drawings

为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。In order to more clearly illustrate the specific implementation of the present invention or the technical solutions in the prior art, the following will briefly introduce the accompanying drawings that need to be used in the specific implementation or description of the prior art. Obviously, the accompanying drawings in the following description The drawings show some implementations of the present invention, and those skilled in the art can obtain other drawings based on these drawings without creative work.

图1是本发明实施例提供的一种多旋翼无人机容错控制方法的流程图。Fig. 1 is a flowchart of a fault-tolerant control method for a multi-rotor UAV provided by an embodiment of the present invention.

图2是本发明实施例提供的一种发生第一种故障时系统未进行分配优化的横滚角的响应曲线图。Fig. 2 is a response curve diagram of the roll angle for which the system does not perform distribution optimization when the first type of fault occurs according to an embodiment of the present invention.

图3是本发明实施例提供的一种发生第一种故障时系统未进行分配优化的俯仰角的响应曲线图。Fig. 3 is a response curve diagram of the pitch angle at which the system does not perform distribution optimization when the first type of fault occurs according to an embodiment of the present invention.

图4是本发明实施例提供的一种发生第一种故障时系统未进行分配优化的偏航角的响应曲线图。Fig. 4 is a response curve diagram of the yaw angle without allocation optimization of the system when the first type of fault occurs according to an embodiment of the present invention.

图5是本发明实施例提供的一种发生第一种故障时系统进行分配优化与未发生故障时的横滚角对比图。Fig. 5 is a comparison diagram of the roll angle when the system performs allocation optimization when the first type of failure occurs and when no failure occurs according to an embodiment of the present invention.

图6是本发明实施例提供的一种发生第一种故障时系统进行分配优化与未发生故障时的俯仰角对比图。Fig. 6 is a comparison diagram of pitch angles when the system performs allocation optimization when the first type of failure occurs and when no failure occurs according to an embodiment of the present invention.

图7是本发明实施例提供的一种发生第一种故障时系统进行分配优化与未发生故障时的偏航角对比图。Fig. 7 is a comparison diagram of the yaw angle when the system performs allocation optimization when the first type of failure occurs and when no failure occurs according to an embodiment of the present invention.

图8是本发明实施例提供的一种发生第二种故障时系统未进行分配优化的横滚角的响应曲线图。Fig. 8 is a response curve diagram of the roll angle for which the system does not perform distribution optimization when the second type of fault occurs according to an embodiment of the present invention.

图9是本发明实施例提供的一种发生第二种故障时系统未进行分配优化的俯仰角的响应曲线图。Fig. 9 is a response curve diagram of the pitch angle for which the system does not perform allocation optimization when the second type of fault occurs according to an embodiment of the present invention.

图10是本发明实施例提供的一种发生第二种故障时系统未进行分配优化的偏航角的响应曲线图。Fig. 10 is a response curve diagram of the yaw angle without allocation optimization of the system when the second type of fault occurs according to an embodiment of the present invention.

图11是本发明实施例提供的一种发生第二种故障时系统进行分配优化与未发生故障时的横滚角对比图。Fig. 11 is a comparison diagram of the roll angle when the system performs allocation optimization when the second type of failure occurs and when no failure occurs according to an embodiment of the present invention.

图12是本发明实施例提供的一种发生第二种故障时系统进行分配优化与未发生故障时的俯仰角对比图。Fig. 12 is a comparison diagram of pitch angles when the system performs allocation optimization when the second type of failure occurs and when no failure occurs according to an embodiment of the present invention.

图13是本发明实施例提供的一种发生第二种故障时系统进行分配优化与未发生故障时的偏航角对比图。Fig. 13 is a comparison diagram of the yaw angle when the system performs allocation optimization when the second type of failure occurs and when no failure occurs according to an embodiment of the present invention.

图14是本发明实施例提供的一种多旋翼无人机容错控制系统的结构图。Fig. 14 is a structural diagram of a fault-tolerant control system for a multi-rotor UAV provided by an embodiment of the present invention.

具体实施方式detailed description

为使本发明的目的、技术方案和优点更加清楚,下面结合附图和具体实施例对发明进行清楚、完整的描述。In order to make the purpose, technical solution and advantages of the present invention clearer, the invention will be clearly and completely described below in conjunction with the accompanying drawings and specific embodiments.

当前,多旋翼无人机主要分为:四旋翼无人机和四旋翼以上的多旋翼无人机,四旋翼无人机缺少旋翼组件的冗余,一旦发送故障,飞行姿态发生突变,在某些应用领域会造成不可估计的后果。而对于四旋翼以上的多旋翼无人机来说,可以通过优化控制算法,使多旋翼无人机具有很好的容错能力、更大的负载能力和更高的稳定性,以保证无人机在遇到强外力干扰或部分电机受损时仍具有良好的稳定性和安全性,从而可以携带更多的任务设备,完成更复杂的任务。At present, multi-rotor UAVs are mainly divided into: quad-rotor UAVs and multi-rotor UAVs with more than four rotors. Quad-rotor UAVs lack the redundancy of rotor components. Once the transmission fails, the flight attitude will change suddenly. Some application areas can cause incalculable consequences. For multi-rotor UAVs with four or more rotors, the control algorithm can be optimized to make the multi-rotor UAV have good fault tolerance, greater load capacity and higher stability, so as to ensure the stability of the UAV. It still has good stability and safety when encountering strong external force interference or part of the motor is damaged, so that it can carry more mission equipment and complete more complex missions.

基于此,本发明涉及的多旋翼无人机指代的是四旋翼以上的多旋翼无人机。Based on this, the multi-rotor UAV involved in the present invention refers to a multi-rotor UAV with more than four rotors.

图1和图2是本发明提供的多旋翼无人机容错控制方法的流程图,如图1和图2所示,该容错控制方法包括:Fig. 1 and Fig. 2 are the flowcharts of the fault-tolerant control method of the multi-rotor UAV provided by the present invention, as shown in Fig. 1 and Fig. 2, the fault-tolerant control method comprises:

S101:采用改进的线性自抗扰控制的姿态控制算法控制无人机的飞行,以保证无人机在飞行过程的鲁棒性。S101: Control the flight of the UAV by using the attitude control algorithm of the improved linear ADRC to ensure the robustness of the UAV during the flight.

其中,改进的线性自抗扰控制的姿态控制算法包括:Among them, the attitude control algorithm of the improved linear ADRC includes:

(1)安排过渡过程:采用以下公式,通过二阶环节将输入的突变信号转化为缓变信号,然后使输出信号达到期望的输入信号:(1) Arrange the transition process: use the following formula to convert the input mutation signal into a slow-change signal through the second-order link, and then make the output signal reach the desired input signal:

Figure BDA0002327484090000071
其中,G(s)代表二阶环节的传递函数,T代表二阶环节的时间常数,T可取期望过渡时间的
Figure BDA0002327484090000072
s代表传递函数中的变量符号,是一个复参数。
Figure BDA0002327484090000071
Among them, G(s) represents the transfer function of the second-order link, T represents the time constant of the second-order link, and T can be taken as the expected transition time
Figure BDA0002327484090000072
s represents the variable symbol in the transfer function and is a complex parameter.

(2)线性扩张状态观测器:采用以下状态空间方程和公式,实现对模型中各变量进行实时跟踪:(2) Linear expansion state observer: use the following state space equations and formulas to realize real-time tracking of each variable in the model:

Figure BDA0002327484090000073
其中,x1,x2,x3分别代表系统的状态变量,
Figure BDA0002327484090000074
b0代表估计的控制增益,w代表外部扰动,y代表模型的输出,u代表模型的输入;
Figure BDA0002327484090000073
Among them, x 1 , x 2 , x 3 represent the state variables of the system respectively,
Figure BDA0002327484090000074
b 0 represents the estimated control gain, w represents the external disturbance, y represents the output of the model, and u represents the input of the model;

Figure BDA0002327484090000075
其中,z1,z2,z3分别代表线性扩张状态观测器的系统状态变量,β123分别代表线性扩张状态观测器的增益。
Figure BDA0002327484090000075
Among them, z 1 , z 2 , z 3 represent the system state variables of the linearly extended state observer respectively, and β 1 , β 2 , β 3 represent the gain of the linearly extended state observer respectively.

需要说明的是,多旋翼无人机的主要通过改变每个旋翼对应的电机的转速来改变该无人机飞行的姿态,电机的转速是通过改变PWM信号的占空比改变的,电机转速发生变化后,电机产生的扬力和转矩发生变化,同时根据多旋翼无人机自身电机的位置分布,决定了横滚方向、俯仰方向的转矩和偏航方向的反转矩。横滚方向和俯仰方向的变化产生了飞机X,Y轴方向的线速度,偏航角的变化为飞机航向的变化,扬力的变化产生了飞机Z轴即高度的变化。It should be noted that the multi-rotor UAV mainly changes the flying attitude of the UAV by changing the speed of the motor corresponding to each rotor. The speed of the motor is changed by changing the duty cycle of the PWM signal. After the change, the lifting force and torque generated by the motor change, and at the same time, according to the position distribution of the multi-rotor UAV's own motor, the torque in the roll direction, the pitch direction and the counter torque in the yaw direction are determined. The change of roll direction and pitch direction produces the linear velocity of the aircraft in the X and Y axis directions, the change of yaw angle is the change of the aircraft heading, and the change of lift force produces the change of the aircraft Z axis, that is, the height.

其中,机体坐标系定义为:原点取在无人机的质心,坐标系与机体固连;X轴与机身设计的纵轴平行,且处于无人机的对称平面内,指向前方;Y轴垂直于无人机对称平面指向右方;Z轴在无人机对称平面内,且垂直于X轴指向下方。整个坐标系符合欧拉坐标系右手定则。地面坐标系即惯性坐标系定义为:采用北东地坐标系,XE轴指向北面,YE轴指向东面,而ZE轴指向地心方向。地面坐标系为仿真实验的环境中的坐标系。Among them, the body coordinate system is defined as: the origin is taken at the center of mass of the drone, and the coordinate system is fixedly connected to the body; the X-axis is parallel to the longitudinal axis of the fuselage design, and is in the plane of symmetry of the drone, pointing forward; the Y-axis Pointing to the right perpendicular to the plane of symmetry of the drone; the Z-axis is within the plane of symmetry of the drone and pointing downwards perpendicular to the X-axis. The entire coordinate system conforms to the right-hand rule of the Euler coordinate system. The ground coordinate system, that is, the inertial coordinate system, is defined as follows: using the northeast ground coordinate system, the XE axis points to the north, the YE axis points to the east, and the ZE axis points to the center of the earth. The ground coordinate system is the coordinate system in the environment of the simulation experiment.

无人机根据用户操作所产生的虚拟控制指令经控制分配信息分配为每个电机的实际控制指令。在四个电机正常的情况下,控制分配信息为固定值,不发生改变。但当无人机的部分电机发生故障时,电机无法正确响应分配的实际控制指令,即电机对于相同的控制指令无法调整为对应的转速,不能满足控制的要求,因此需对控制分配信息进行优化再分配。The virtual control command generated by the UAV according to the user's operation is assigned as the actual control command of each motor through the control distribution information. When the four motors are normal, the control assignment information is a fixed value and does not change. However, when some motors of the UAV fail, the motors cannot correctly respond to the assigned actual control commands, that is, the motors cannot be adjusted to the corresponding speed for the same control commands, and cannot meet the control requirements. Therefore, it is necessary to optimize the control distribution information redistribute.

S102:当检测到无人机的部分电机异常时,构建故障矩阵Ri,i为大于等于0小于等于无人机中所有电机个数的整数。S102: When some motors of the UAV are detected to be abnormal, construct a fault matrix Ri, where i is an integer greater than or equal to 0 and less than or equal to the number of all motors in the UAV.

其中,该步骤包括:Among them, this step includes:

S1021:实时检测无人机的所有电机。S1021: Detect all motors of the drone in real time.

S1022:当检测到无人机的部分电机异常时,计算故障电机的输出与无故障时的输出的比值,根据比值构建故障矩阵。S1022: When some motors of the UAV are detected to be abnormal, calculate the ratio of the output of the faulty motor to the output when there is no fault, and construct a fault matrix according to the ratio.

关于判断部分电机异常的方式,在一种可能的实现方式中,在无人机飞行过程中,用户可通过推动遥控器的摇杆以向无人机飞行控制系统发送控制指令,无人机飞行控制系统接收到控制指令后,通过内部的姿态解算环节和姿态控制环节对每个电调输出一定的控制量,之后每个电调输出与该控制量对应的PWM值以控制与该电调连接的电机的转速,即,每个电调输出的PWM值与控制指令相对应。因此,当某一个电调输出的PWM值不与该控制指令相对应时,则确定该电调对应的电机发生故障。Regarding the method of judging the abnormality of some motors, in one possible implementation, during the flight of the drone, the user can push the joystick of the remote control to send control instructions to the flight control system of the drone, and the drone will fly After the control system receives the control command, it outputs a certain amount of control to each ESC through the internal attitude calculation link and attitude control link, and then each ESC outputs a PWM value corresponding to the control amount to control the ESC. The speed of the connected motor, that is, the PWM value output by each ESC corresponds to the control command. Therefore, when the PWM value output by a certain ESC does not correspond to the control command, it is determined that the motor corresponding to the ESC is faulty.

当检测到无人机的部分电机异常时,无人机飞行控制系统获取所有电机的输出值,即,电机的效率。其中,本发明实施例中,将未发生故障的电机的效率设为1;而发生故障的电机效率为故障电机的输出相对于正常时的输出的量化值。将每个电机的效率分别以k1,k2…ki表示,则构建得到的故障矩阵Ri为:When some motors of the drone are detected to be abnormal, the flight control system of the drone acquires the output values of all the motors, that is, the efficiency of the motors. Wherein, in the embodiment of the present invention, the efficiency of a motor without a fault is set as 1; and the efficiency of a motor with a fault is a quantized value of the output of the faulty motor relative to the normal output. The efficiency of each motor is represented by k 1 , k 2 ...k i respectively, then the constructed fault matrix R i is:

Ri=diag[k1,k2…ki]。R i =diag[k 1 , k 2 . . . k i ].

例如:该多旋翼无人机为六旋翼无人机,将该六旋翼无人机上的六个电机分别标号为1-6,当1号电机和2号电机由于发生故障导致其电机的效率分别为k1,k2时,该六旋翼无人机的故障矩阵为:For example: the multi-rotor UAV is a six-rotor UAV, and the six motors on the six-rotor UAV are respectively marked as 1-6. When k 1 and k 2 , the fault matrix of the hexacopter UAV is:

Figure BDA0002327484090000091
Figure BDA0002327484090000091

S103:基于故障矩阵Ri在线建立故障模型。S103: Establish a fault model online based on the fault matrix R i .

采用以下公式,在线建立故障模型:The fault model is established online using the following formula:

Figure BDA0002327484090000092
Figure BDA0002327484090000092

其中,

Figure BDA0002327484090000093
分别代表大地坐标系下的位置加速度,
Figure BDA0002327484090000094
分别代表在大地坐标系下无人机的飞行器姿态角的角加速度,
Figure BDA0002327484090000095
θ,ψ分别代表横滚角、俯仰角和偏航角,Ix,Iy,Iz分别代表无人机机身在三个方向的转动惯量,m代表无人机的质量,g代表重力加速度,UR,UP,UY,UT分别代表无人机的电机均无故障时的横滚力矩、俯仰力矩、偏航力矩以及升力,fp,fq,fr,fz分别表示横滚力矩误差、俯仰力矩误差、偏航力矩误差以及升力误差。in,
Figure BDA0002327484090000093
Represent the position acceleration in the earth coordinate system,
Figure BDA0002327484090000094
Represent the angular acceleration of the UAV's aircraft attitude angle in the earth coordinate system,
Figure BDA0002327484090000095
θ, ψ represent the roll angle, pitch angle and yaw angle respectively, I x , I y , I z represent the moments of inertia of the UAV fuselage in three directions, m represents the mass of the UAV, and g represents the gravity Acceleration, U R , U P , U Y , U T respectively represent the rolling moment, pitching moment, yaw moment and lift when the motors of the UAV are not faulty, f p , f q , f r , f z respectively Indicates roll moment error, pitch moment error, yaw moment error and lift error.

上述

Figure BDA0002327484090000096
可通过传感器得到,传感器包括测量加速度的三轴加速度计和测量角速度的三轴陀螺仪。the above
Figure BDA0002327484090000096
Available through sensors that include a three-axis accelerometer that measures acceleration and a three-axis gyroscope that measures angular velocity.

fp,fq,fr,fz可由故障矩阵Ri得到,需要说明的是,fp,fq,fr,fz与无人机的旋翼数量及旋翼的布局类型相关,一般情况下,一旦确定无人机的旋翼数量和旋翼的类型,即可得到fp,fq,fr,fzf p , f q , f r , f z can be obtained from the fault matrix R i . It should be noted that f p , f q , f r , f z are related to the number of rotors of the UAV and the layout type of the rotors. In general Next, once the number of rotors and the type of rotors are determined, f p , f q , f r , f z can be obtained.

在一种可能的实现方式中,该多旋无人机为“X”型六旋翼无人机时,将该六旋翼无人机上的六个电机分别标号为1-6,各轴之间夹角为60°,可采用以下方式,得到无人机的电机均无故障时的横滚力矩UR、俯仰力矩UP、偏航力矩UY以及升力UTIn a possible implementation, when the multi-rotor UAV is an "X" type six-rotor UAV, the six motors on the six-rotor UAV are respectively marked as 1-6, and each axis is clamped The angle is 60°, and the following methods can be used to obtain the rolling moment U R , pitching moment U P , yaw moment U Y and lift U T when the motors of the UAV are not faulty:

Figure BDA0002327484090000101
Figure BDA0002327484090000101

则当1号电机发生故障后其电机效率为k1,其他电机未发生故障时,故障矩阵R6为:

Figure BDA0002327484090000102
可采用以下方式得到无人机的部分电机出现故障时的的横滚力矩
Figure BDA0002327484090000103
俯仰力矩
Figure BDA0002327484090000104
偏航力矩
Figure BDA0002327484090000105
以及升力
Figure BDA0002327484090000106
Then when the No. 1 motor fails, its motor efficiency is k 1 , and other motors do not fail, the fault matrix R 6 is:
Figure BDA0002327484090000102
The rolling moment when some motors of the UAV fail can be obtained in the following way
Figure BDA0002327484090000103
pitching moment
Figure BDA0002327484090000104
Yaw moment
Figure BDA0002327484090000105
and lift
Figure BDA0002327484090000106

Figure BDA0002327484090000107
Figure BDA0002327484090000107

其中,上述

Figure BDA0002327484090000108
分别代表1-6号电机未出故障时每个电机所对应的转速。Among them, the above
Figure BDA0002327484090000108
Respectively represent the corresponding speed of each motor when the No. 1-6 motors are not faulty.

之后,采用以下公式,求得横滚力矩误差fp、俯仰力矩误差fq、偏航力矩误差fr以及升力误差fzAfterwards, use the following formulas to obtain roll moment error f p , pitch moment error f q , yaw moment error f r and lift error f z :

Figure BDA0002327484090000109
Figure BDA0002327484090000109

S104:基于故障矩阵Ri得到无人机上所有电机的控制分配信息。S104: Obtain control distribution information of all motors on the UAV based on the fault matrix R i .

采用以下公式,得到优化后的分配矩阵Nf,将优化后的分配矩阵Nf作为无人机上所有电机的控制分配信息:Using the following formula, the optimized distribution matrix N f is obtained, and the optimized distribution matrix N f is used as the control distribution information of all motors on the UAV:

Nf=Af-1Nf = Af -1 ;

Nf=AfT(Af·AfT)-1;其中,Af代表部分电机故障后的控制效率矩阵,AfT代表Af的转置;N f = A fT (A f A fT ) -1 ; wherein, A f represents the control efficiency matrix after a partial motor failure, and A fT represents the transposition of A f ;

Af=ARi;其中,A代表故障前的控制效率矩阵;A f = AR i ; where, A represents the control efficiency matrix before failure;

Figure BDA0002327484090000111
其中,b为升力系数,l为无人机的轴距,d为反扭矩系数。
Figure BDA0002327484090000111
Among them, b is the lift coefficient, l is the wheelbase of the UAV, and d is the counter torque coefficient.

S105:采用改进的线性自抗扰控制的姿态控制算法控制故障模型下的无人机的飞行,并按照电机的控制分配信息控制无人机,以达到所需的姿态和高度。S105: Use the attitude control algorithm of the improved linear ADRC to control the flight of the UAV under the fault model, and control the UAV according to the control distribution information of the motor to achieve the required attitude and height.

根据实际控制指令调整每个电机的转速,进而使每个电机驱动无人机达到所需的姿态和高度。姿态包括俯仰、滚转和偏航。Adjust the speed of each motor according to the actual control instructions, and then make each motor drive the UAV to achieve the required attitude and height. Attitude includes pitch, roll and yaw.

例如,可采用以下公式调整每个电机的转速:For example, the following formula can be used to adjust the speed of each motor:

τf=[UR UP UY UT]T;其中,τf代表力矩矩阵;τ f =[U R U P U Y U T ] T ; where, τ f represents the moment matrix;

Figure BDA0002327484090000112
其中,
Figure BDA0002327484090000113
表示当无人机部分电机异常时,每个电机所对应的转速。
Figure BDA0002327484090000112
in,
Figure BDA0002327484090000113
Indicates the corresponding speed of each motor when some motors of the UAV are abnormal.

本发明实施例中,无人机飞行控制系统根据用户操作分布虚拟控制指令给执行机构,执行机构中的电机的转速发生变化,则无人机飞行控制系统在采用改进的线性自抗扰控制的姿态控制算法控制无人机飞行过程中,基于新的控制分配信息对每个电机的控制分配信息进行优化;还会按照建立的故障模型,补偿由于电机故障而产生的横滚力矩误差fp、俯仰力矩误差fq、偏航力矩误差fr以及升力误差fz,从而使得无人机的高度和姿态发生变化,达到无人机所要达到的高度和姿态。In the embodiment of the present invention, the UAV flight control system distributes virtual control instructions to the actuator according to the user operation, and the speed of the motor in the actuator changes, the UAV flight control system adopts the improved linear active disturbance rejection control. During the flight process of the UAV controlled by the attitude control algorithm, the control distribution information of each motor is optimized based on the new control distribution information; and the rolling moment error f p , The pitch moment error f q , the yaw moment error f r and the lift error f z make the height and attitude of the UAV change, reaching the desired height and attitude of the UAV.

另外,本发明实施例采用六旋翼无人机为对象,对容错控制方法进行了仿真验证,如图2-图13所示,且图2-图13均为给定期望的横滚角、俯仰角和偏航角为15°,给定期望的高度为1m的情况下得到的曲线图。In addition, the embodiment of the present invention uses the six-rotor UAV as the object, and simulates and verifies the fault-tolerant control method, as shown in Figures 2-13, and Figures 2-13 are all given the desired roll angle and pitch The resulting plots are given a desired altitude of 1m for a 15° angle and a yaw angle of 15°.

其中,图2-图4分别为1号电机的效率为1/5时,系统未进行分配优化的横滚角、俯仰角和偏航角响应曲线图。由图2-图4可知,当1号电机力效为1/5时且不采用本发明实施例提供的无人机容错控制方法时,三个姿态角无法跟随输入的期望姿态角达到稳定。Among them, Fig. 2-Fig. 4 respectively show the response curves of roll angle, pitch angle and yaw angle when the efficiency of No. 1 motor is 1/5, and the system is not optimized for distribution. It can be seen from Figures 2 to 4 that when the No. 1 motor has a power effect of 1/5 and the UAV fault-tolerant control method provided by the embodiment of the present invention is not used, the three attitude angles cannot follow the input expected attitude angles to achieve stability.

图5-图7分别为1号电机的效率为1/5时,系统进行分配优化与未发生故障时的横滚角、俯仰角和偏航角的对比图,其中,(a)线表示系统进行分配优化的三个姿态角,(b)线代表未发生故障时的三个姿态角。由图5-图7可知,当采用了本发明实施例提供的无人机容错控制方法后,与无故障时的飞行状况相比,虽然没稳定时间较长,但是依然可以很好的跟踪期望值并达到最终稳定。Figures 5-7 are the comparison diagrams of the roll angle, pitch angle and yaw angle when the efficiency of the No. 1 motor is 1/5, the distribution optimization of the system and the failure-free time, where (a) line represents the system The three attitude angles for allocation optimization, the line (b) represents the three attitude angles without failure. It can be seen from Fig. 5-Fig. 7 that when the UAV fault-tolerant control method provided by the embodiment of the present invention is adopted, compared with the flight condition when there is no fault, although the unsteady time is longer, the expected value can still be tracked very well. and reach final stability.

图8-图10分别为1号电机的效率为1/2时,系统未进行分配优化的横滚角、俯仰角和偏航角响应曲线图。由图8-图10可知,当1号电机力效为1/2时且不采用本发明实施例提供的无人机容错控制方法时,三个姿态角无法跟随输入的期望姿态角达到稳定且震荡明显。Figures 8-10 are the response curves of roll angle, pitch angle, and yaw angle when the efficiency of the No. 1 motor is 1/2, and the system is not optimized for distribution. From Figures 8 to 10, it can be seen that when the power effect of the No. 1 motor is 1/2 and the UAV fault-tolerant control method provided by the embodiment of the present invention is not used, the three attitude angles cannot follow the input expected attitude angles to reach stability and Vibration is obvious.

图11-图13分别为1号电机的效率为1/2时,系统进行分配优化与未发生故障时的横滚角、俯仰角和偏航角的对比图,其中,(a)线表示系统进行分配优化的三个姿态角,(b)线代表未发生故障时的三个姿态角。由图11-图13可知,当采用了本发明实施例提供的无人机容错控制方法后,与无故障时的飞行状况相比,虽然没稳定时间较长,但是依然可以很好的跟踪期望值并达到最终稳定。Figures 11-13 are the comparison diagrams of the roll angle, pitch angle and yaw angle when the efficiency of the No. 1 motor is 1/2, the distribution optimization of the system and the failure-free time, where the (a) line represents the system The three attitude angles for allocation optimization, the line (b) represents the three attitude angles without failure. From Figures 11 to 13, it can be seen that when the UAV fault-tolerant control method provided by the embodiment of the present invention is adopted, compared with the flight condition without failure, although the unsteady time is longer, the expected value can still be tracked very well and reach final stability.

本发明实施例提供的一种多旋翼无人机容错控制方法,采用改进的线性自抗扰控制的姿态控制算法控制无人机的飞行,以保证无人机在飞行过程的鲁棒性;当检测到无人机的部分电机异常时,构建故障矩阵Ri;基于故障矩阵Ri在线建立故障模型;基于故障矩阵Ri得到无人机上所有电机的控制分配信息;采用改进的线性自抗扰控制的姿态控制算法控制故障模型下的无人机的飞行,并按照电机的控制分配信息控制无人机,以达到所需的姿态和高度。在无人机的电机未发生故障时,选用的基本控制律为改进的线性自抗扰控制的姿态控制算法,该控制算法对干扰具有较强的鲁棒能力。而且,当无人机的部分电机发生故障时,还能基于发生故障后的故障矩阵,得到新的飞行模型-故障模型,以及所有电机的控制分配信息,从而保证无人机在该故障模式下飞行的同时,还能减少故障电机的使用,从而使无人机达到平稳飞行状态。A fault-tolerant control method for a multi-rotor UAV provided by an embodiment of the present invention uses an improved linear active disturbance rejection control attitude control algorithm to control the flight of the UAV to ensure the robustness of the UAV during flight; When some motors of the UAV are detected to be abnormal, the fault matrix R i is constructed; the fault model is established online based on the fault matrix R i ; the control distribution information of all motors on the UAV is obtained based on the fault matrix R i ; the improved linear active disturbance rejection The attitude control algorithm of the control controls the flight of the UAV under the fault model, and controls the UAV according to the control distribution information of the motor to achieve the required attitude and height. When the motor of the UAV does not fail, the basic control law selected is the attitude control algorithm of the improved linear active disturbance rejection control, which has strong robustness to disturbance. Moreover, when some motors of the UAV fail, based on the failure matrix after the failure, a new flight model-fault model and the control distribution information of all motors can be obtained, so as to ensure that the UAV can operate in this failure mode. While flying, it can also reduce the use of faulty motors, so that the UAV can achieve a stable flight state.

图14是本发明实施例提供的一种多旋翼无人机容错控制系统的结构图,如图14所示,其包括:Figure 14 is a structural diagram of a multi-rotor UAV fault-tolerant control system provided by an embodiment of the present invention, as shown in Figure 14, which includes:

第一控制模块1401,用于采用改进的线性自抗扰控制的姿态控制算法控制无人机的飞行,以保证无人机在飞行过程的鲁棒性;The first control module 1401 is used to control the flight of the UAV by adopting the attitude control algorithm of the improved linear ADRC to ensure the robustness of the UAV during the flight;

故障矩阵构建模块1402,用于当检测到无人机的部分电机异常时,构建故障矩阵Ri,i为大于等于0小于等于无人机中所有电机个数的整数;The fault matrix construction module 1402 is used to construct a fault matrix R i when some motors of the drone are abnormal, and i is an integer greater than or equal to 0 and less than or equal to the number of all motors in the drone;

故障模型建立模块1403,用于基于故障矩阵Ri在线建立故障模型;A fault model building module 1403, configured to build a fault model online based on the fault matrix R i ;

分配信息获取模块1404,用于基于故障矩阵Ri得到无人机上所有电机的控制分配信息;Distribution information acquisition module 1404, for obtaining the control distribution information of all motors on the UAV based on the fault matrix R i ;

第二控制模块1405,采用改进的线性自抗扰控制的姿态控制算法控制故障模型下的无人机的飞行,并按照电机的控制分配信息控制无人机,以达到所需的姿态和高度。The second control module 1405 uses the attitude control algorithm of the improved linear ADRC to control the flight of the UAV under the fault model, and controls the UAV according to the control distribution information of the motor to achieve the required attitude and altitude.

其中,第一控制模块1401和第二控制模块1405可以为同一个控制模块。Wherein, the first control module 1401 and the second control module 1405 may be the same control module.

在一种可能的实现方式中,故障矩阵构建模块1402包括:In a possible implementation, the fault matrix construction module 1402 includes:

安排过渡过程单元,用于采用以下公式,通过二阶环节将输入的突变信号转化为缓变信号,然后使输出信号达到期望的输入信号:The transition process unit is arranged to convert the input sudden change signal into a slow change signal through the second-order link by using the following formula, and then make the output signal reach the desired input signal:

Figure BDA0002327484090000131
其中,其中,G(s)代表二阶环节的传递函数,T代表二阶环节的时间常数,s代表传递函数中的变量符号;
Figure BDA0002327484090000131
Wherein, G(s) represents the transfer function of the second-order link, T represents the time constant of the second-order link, and s represents the variable symbol in the transfer function;

线性扩张状态观测器单元,用于采用以下状态空间方程和公式,实现对模型中各变量进行实时跟踪:The Linear Extended State Observer unit is used to implement real-time tracking of variables in the model using the following state space equations and formulas:

Figure BDA0002327484090000132
其中,x1,x2,x3分别代表所系统的状态变量,
Figure BDA0002327484090000133
b0代表估计的控制增益,w代表外部扰动,y代表模型的输出,u代表模型的输入;
Figure BDA0002327484090000132
Among them, x 1 , x 2 , x 3 respectively represent the state variables of the system,
Figure BDA0002327484090000133
b 0 represents the estimated control gain, w represents the external disturbance, y represents the output of the model, and u represents the input of the model;

Figure BDA0002327484090000134
其中,z1,z2,z3分别代表线性扩张状态观测器的系统状态变量,β123分别代表线性扩张状态观测器的增益。
Figure BDA0002327484090000134
Among them, z 1 , z 2 , z 3 represent the system state variables of the linearly extended state observer respectively, and β 1 , β 2 , β 3 represent the gain of the linearly extended state observer respectively.

在一种可能的实现方式中,故障模型建立模块1403包括:In a possible implementation, the fault model building module 1403 includes:

检测单元,用于实时检测无人机的所有电机;A detection unit for real-time detection of all motors of the drone;

故障模型建立单元,用于当检测到无人机的部分电机异常时,计算故障电机的输出与无故障时的输出的比值,根据比值构建故障矩阵RiThe fault model building unit is used to calculate the ratio of the output of the faulty motor to the output when there is no fault when some motors of the UAV are detected to be abnormal, and construct a fault matrix R i according to the ratio.

在一种可能的实现方式中,故障模型建立模块1403包括:In a possible implementation, the fault model building module 1403 includes:

故障模型建立单元,用于采用以下公式,在线建立故障模型:The fault model establishing unit is used to establish the fault model online by adopting the following formula:

Figure BDA0002327484090000135
Figure BDA0002327484090000135

其中

Figure BDA0002327484090000136
分别代表大地坐标系下的位置加速度,
Figure BDA0002327484090000137
分别代表在大地坐标系下无人机的飞行器姿态角的角加速度,
Figure BDA0002327484090000138
θ,ψ分别代表横滚角、俯仰角和偏航角,Ix,Iy,Iz分别代表无人机机身在三个方向的转动惯量,m代表无人机的质量,g代表重力加速度,UR,UP,UY,UT分别代表无人机的电机均无故障时的横滚力矩、俯仰力矩、偏航力矩以及升力,fp,fq,fr,fz分别表示横滚力矩误差、俯仰力矩误差、偏航力矩误差以及升力误差。in
Figure BDA0002327484090000136
Represent the position acceleration in the earth coordinate system,
Figure BDA0002327484090000137
Represent the angular acceleration of the UAV's aircraft attitude angle in the earth coordinate system,
Figure BDA0002327484090000138
θ, ψ represent the roll angle, pitch angle and yaw angle respectively, I x , I y , I z represent the moments of inertia of the UAV fuselage in three directions, m represents the mass of the UAV, and g represents the gravity Acceleration, U R , U P , U Y , U T respectively represent the rolling moment, pitching moment, yaw moment and lift when the motors of the UAV are not faulty, f p , f q , f r , f z respectively Indicates roll moment error, pitch moment error, yaw moment error and lift error.

在一种可能的实现方式中,分配信息获取模块1404包括:In a possible implementation manner, the allocation information obtaining module 1404 includes:

分配信息获取单元,用于采用以下公式,得到优化后的分配矩阵Nf,将优化后的分配矩阵Nf作为无人机上所有电机的控制分配信息:The allocation information acquisition unit is used to obtain the optimized allocation matrix Nf by using the following formula, and use the optimized allocation matrix Nf as the control allocation information of all motors on the drone:

Nf=Af-1Nf = Af -1 ;

Nf=AfT(Af·AfT)-1;其中,Af代表部分电机故障后的控制效率矩阵,AfT代表Af的转置;N f = A fT (A f A fT ) -1 ; wherein, A f represents the control efficiency matrix after a partial motor failure, and A fT represents the transposition of A f ;

Af=ARi;其中,A代表故障前的控制效率矩阵;A f = AR i ; where, A represents the control efficiency matrix before failure;

Figure BDA0002327484090000141
其中,b为升力系数,l为无人机的轴距,d为反扭矩系数。
Figure BDA0002327484090000141
Among them, b is the lift coefficient, l is the wheelbase of the UAV, and d is the counter torque coefficient.

本发明实施例提供的一种多旋翼无人机容错控制系统,采用改进的线性自抗扰控制的姿态控制算法控制无人机的飞行,以保证无人机在飞行过程的鲁棒性;当检测到无人机的部分电机异常时,构建故障矩阵Ri;基于故障矩阵Ri在线建立故障模型;基于故障矩阵Ri得到无人机上所有电机的控制分配信息;采用改进的线性自抗扰控制的姿态控制算法控制故障模型下的无人机的飞行,并按照电机的控制分配信息控制无人机,以达到所需的姿态和高度。在无人机的电机未发生故障时,选用的基本控制律为改进的线性自抗扰控制的姿态控制算法,该控制算法对干扰具有较强的鲁棒能力。而且,当无人机的部分电机发生故障时,还能基于发生故障后的故障矩阵,得到新的飞行模型-故障模型,以及所有电机的控制分配信息,从而保证无人机在该故障模式下飞行的同时,还能减少故障电机的使用,从而使无人机达到平稳飞行状态。A fault-tolerant control system for a multi-rotor UAV provided by an embodiment of the present invention uses an improved linear active disturbance rejection control attitude control algorithm to control the flight of the UAV to ensure the robustness of the UAV during flight; When some motors of the UAV are detected to be abnormal, the fault matrix R i is constructed; the fault model is established online based on the fault matrix R i ; the control distribution information of all motors on the UAV is obtained based on the fault matrix R i ; the improved linear active disturbance rejection The attitude control algorithm of the control controls the flight of the UAV under the fault model, and controls the UAV according to the control distribution information of the motor to achieve the required attitude and altitude. When the motor of the UAV does not fail, the basic control law selected is the attitude control algorithm of the improved linear active disturbance rejection control, which has strong robustness to disturbance. Moreover, when some motors of the UAV fail, based on the failure matrix after the failure, a new flight model-fault model and the control distribution information of all motors can be obtained, so as to ensure that the UAV can operate in this failure mode. While flying, it can also reduce the use of faulty motors, so that the UAV can achieve a stable flight state.

最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围。Finally, it should be noted that: the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: it can still be Modifications are made to the technical solutions described in the foregoing embodiments, or equivalent replacements are made to some of the technical features; these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims (8)

1.一种多旋翼无人机容错控制方法,其特征在于,其包括:1. A multi-rotor unmanned aerial vehicle fault-tolerant control method, is characterized in that, it comprises: 采用改进的线性自抗扰控制的姿态控制算法控制所述无人机的飞行,以保证所述无人机在飞行过程的鲁棒性;The attitude control algorithm of the improved linear active disturbance rejection control is used to control the flight of the UAV, so as to ensure the robustness of the UAV during flight; 当检测到所述无人机的部分电机异常时,构建故障矩阵Ri,i为大于等于0小于等于所述无人机中所有电机个数的整数;When some motors of the UAV are detected to be abnormal, a fault matrix R i is constructed, where i is an integer greater than or equal to 0 and less than or equal to the number of all motors in the UAV; 基于所述故障矩阵Ri在线建立故障模型;Establishing a fault model online based on the fault matrix R i ; 基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息;Obtain the control allocation information of all motors on the UAV based on the fault matrix R i ; 采用所述改进的线性自抗扰控制的姿态控制算法控制所述故障模型下的无人机的飞行,并按照所述电机的控制分配信息控制所述无人机,以达到所需的姿态和高度;Adopt the attitude control algorithm of the improved linear ADRC to control the flight of the unmanned aerial vehicle under the fault model, and control the unmanned aerial vehicle according to the control distribution information of the motor to achieve the required attitude and high; 所述基于所述故障矩阵Ri在线建立故障模型,包括:The online establishment of a fault model based on the fault matrix R i includes: 采用以下公式,在线建立故障模型:The fault model is established online using the following formula:
Figure FDA0003886206230000011
Figure FDA0003886206230000011
其中
Figure FDA0003886206230000012
分别代表大地坐标系下的位置加速度,
Figure FDA0003886206230000013
分别代表在大地坐标系下所述无人机的飞行器姿态角的角加速度,
Figure FDA0003886206230000014
θ,ψ分别代表横滚角、俯仰角和偏航角,Ix,Iy,Iz分别代表所述无人机机身在三个方向的转动惯量,m代表所述无人机的质量,g代表重力加速度,UR,UP,UY,UT分别代表所述无人机的电机均无故障时的横滚力矩、俯仰力矩、偏航力矩以及升力,fp,fq,fr,fz分别表示横滚力矩误差、俯仰力矩误差、偏航力矩误差以及升力误差。
in
Figure FDA0003886206230000012
Represent the position acceleration in the earth coordinate system,
Figure FDA0003886206230000013
represent the angular acceleration of the aircraft attitude angle of the unmanned aerial vehicle under the earth coordinate system respectively,
Figure FDA0003886206230000014
θ, ψ represent the roll angle, pitch angle and yaw angle respectively, I x , I y , I z represent the moments of inertia of the UAV fuselage in three directions respectively, and m represents the mass of the UAV , g represents the acceleration of gravity, U R , U P , U Y , U T respectively represent the rolling moment, pitching moment, yaw moment and lift when the motors of the UAV are not faulty, f p , f q , f r , f z represent roll moment error, pitch moment error, yaw moment error and lift error respectively.
2.根据权利要求1所述的一种多旋翼无人机容错控制方法,其特征在于,所述改进的线性自抗扰控制的姿态控制算法包括:2. a kind of multi-rotor UAV fault-tolerant control method according to claim 1, is characterized in that, the attitude control algorithm of described improved linear active disturbance rejection control comprises: 安排过渡过程:采用以下公式,通过二阶环节将输入的突变信号转化为缓变信号,然后使输出信号达到期望的输入信号:Arrange the transition process: use the following formula to convert the input sudden change signal into a slow change signal through the second-order link, and then make the output signal reach the desired input signal:
Figure FDA0003886206230000015
其中,其中,G(s)代表二阶环节的传递函数,T代表二阶环节的时间常数,s代表传递函数中的变量符号;
Figure FDA0003886206230000015
Wherein, G(s) represents the transfer function of the second-order link, T represents the time constant of the second-order link, and s represents the variable symbol in the transfer function;
线性扩张状态观测器:采用以下状态空间方程和公式,实现对模型中各变量进行实时跟踪:Linear expansion state observer: use the following state space equations and formulas to realize real-time tracking of variables in the model:
Figure FDA0003886206230000021
其中,x1,x2,x3分别代表所系统的状态变量,
Figure FDA0003886206230000022
b0代表估计的控制增益,w代表外部扰动,y代表所述模型的输出,u代表所述模型的输入;
Figure FDA0003886206230000021
Among them, x 1 , x 2 , x 3 respectively represent the state variables of the system,
Figure FDA0003886206230000022
b 0 represents the estimated control gain, w represents the external disturbance, y represents the output of the model, and u represents the input of the model;
Figure FDA0003886206230000023
其中,z1,z2,z3分别代表所述线性扩张状态观测器的系统状态变量,β123分别代表所述线性扩张状态观测器的增益。
Figure FDA0003886206230000023
Wherein, z 1 , z 2 , z 3 respectively represent the system state variables of the linearly extended state observer, and β 1 , β 2 , β 3 represent the gain of the linearly extended state observer respectively.
3.根据权利要求1所述的一种多旋翼无人机容错控制方法,其特征在于,所述当检测到所述无人机的部分电机异常时,构建故障矩阵Ri,包括:3. A kind of multi-rotor unmanned aerial vehicle fault-tolerant control method according to claim 1, is characterized in that, described when detecting the part motor abnormality of described unmanned aerial vehicle, constructing fault matrix R i , comprising: 实时检测所述无人机的所有电机;Real-time detection of all motors of the drone; 当检测到所述无人机的部分电机异常时,计算故障电机的输出与无故障时的输出的比值,根据所述比值构建故障矩阵RiWhen some motors of the UAV are detected to be abnormal, the ratio of the output of the faulty motor to the output without fault is calculated, and a fault matrix R i is constructed according to the ratio. 4.根据权利要求1所述的一种多旋翼无人机容错控制方法,其特征在于,所述基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息,包括:4. a kind of multi-rotor unmanned aerial vehicle fault-tolerant control method according to claim 1, is characterized in that, described based on described failure matrix R i obtains the control distribution information of all motors on the described unmanned aerial vehicle, comprising: 采用以下公式,得到优化后的分配矩阵Nf,将所述优化后的分配矩阵Nf作为所述无人机上所有电机的控制分配信息:The following formula is used to obtain the optimized allocation matrix N f , and the optimized allocation matrix N f is used as the control allocation information of all motors on the drone: Nf=Af-1Nf = Af -1 ; Nf=AfT(Af·AfT)-1;其中,Af代表部分电机故障后的控制效率矩阵,AfT代表Af的转置;N f = A fT (A f A fT ) -1 ; wherein, A f represents the control efficiency matrix after a partial motor failure, and A fT represents the transposition of A f ; Af=ARi;其中,A代表故障前的控制效率矩阵;A f = AR i ; where, A represents the control efficiency matrix before failure;
Figure FDA0003886206230000024
其中,b为升力系数,l为所述无人机的轴距,d为反扭矩系数。
Figure FDA0003886206230000024
Wherein, b is the lift coefficient, l is the wheelbase of the UAV, and d is the counter torque coefficient.
5.一种多旋翼无人机容错控制系统,其特征在于,其包括:5. A multi-rotor unmanned aerial vehicle fault-tolerant control system, is characterized in that, it comprises: 第一控制模块,用于采用改进的线性自抗扰控制的姿态控制算法控制所述无人机的飞行,以保证所述无人机在飞行过程的鲁棒性;The first control module is used to control the flight of the unmanned aerial vehicle by adopting the attitude control algorithm of the improved linear active disturbance rejection control, so as to ensure the robustness of the unmanned aerial vehicle during flight; 故障矩阵构建模块,用于当检测到所述无人机的部分电机异常时,构建故障矩阵Ri,i为大于等于0小于等于所述无人机中所有电机个数的整数;A fault matrix construction module, used to construct a fault matrix R i when detecting abnormality of some motors of the drone, where i is an integer greater than or equal to 0 and less than or equal to the number of all motors in the drone; 故障模型建立模块,用于基于所述故障矩阵Ri在线建立故障模型;A fault model building module, configured to build a fault model online based on the fault matrix R i ; 分配信息获取模块,用于基于所述故障矩阵Ri得到所述无人机上所有电机的控制分配信息;An allocation information acquisition module, configured to obtain control allocation information of all motors on the UAV based on the fault matrix R i ; 第二控制模块,采用所述改进的线性自抗扰控制的姿态控制算法控制所述故障模型下的无人机的飞行,并按照所述电机的控制分配信息控制所述无人机,以达到所需的姿态和高度;The second control module adopts the attitude control algorithm of the improved linear active disturbance rejection control to control the flight of the UAV under the fault model, and controls the UAV according to the control distribution information of the motor, so as to achieve desired attitude and altitude; 所述故障模型建立模块包括:故障模型建立单元,用于采用以下公式,在线建立故障模型:The fault model establishment module includes: a fault model establishment unit, which is used to establish a fault model online by using the following formula:
Figure FDA0003886206230000031
Figure FDA0003886206230000031
其中
Figure FDA0003886206230000032
分别代表大地坐标系下的位置加速度,
Figure FDA0003886206230000033
分别代表在大地坐标系下所述无人机的飞行器姿态角的角加速度,
Figure FDA0003886206230000034
θ,ψ分别代表横滚角、俯仰角和偏航角,Ix,Iy,Iz分别代表所述无人机机身在三个方向的转动惯量,m代表所述无人机的质量,g代表重力加速度,UR,UP,UY,UT分别代表所述无人机的电机均无故障时的横滚力矩、俯仰力矩、偏航力矩以及升力,fp,fq,fr,fz分别表示横滚力矩误差、俯仰力矩误差、偏航力矩误差以及升力误差。
in
Figure FDA0003886206230000032
Represent the position acceleration in the earth coordinate system,
Figure FDA0003886206230000033
represent the angular acceleration of the aircraft attitude angle of the unmanned aerial vehicle under the earth coordinate system respectively,
Figure FDA0003886206230000034
θ, ψ represent the roll angle, pitch angle and yaw angle respectively, I x , I y , I z represent the moments of inertia of the UAV fuselage in three directions respectively, and m represents the mass of the UAV , g represents the acceleration of gravity, U R , U P , U Y , U T respectively represent the rolling moment, pitching moment, yaw moment and lift when the motors of the UAV are not faulty, f p , f q , f r , f z represent roll moment error, pitch moment error, yaw moment error and lift error respectively.
6.根据权利要求5所述的一种多旋翼无人机容错控制系统,其特征在于,所述故障矩阵构建模块包括:6. a kind of multi-rotor unmanned aerial vehicle fault-tolerant control system according to claim 5, is characterized in that, described failure matrix building block comprises: 安排过渡过程,用于采用以下公式,通过二阶环节将输入的突变信号转化为缓变信号,然后使输出信号达到期望的输入信号:The transition process is arranged to convert the input sudden change signal into a slow change signal through the second-order link by using the following formula, and then make the output signal reach the desired input signal:
Figure FDA0003886206230000041
其中,其中,G(s)代表二阶环节的传递函数,T代表二阶环节的时间常数,s代表传递函数中的变量符号;
Figure FDA0003886206230000041
Wherein, G(s) represents the transfer function of the second-order link, T represents the time constant of the second-order link, and s represents the variable symbol in the transfer function;
线性扩张状态观测器单元,用于采用以下状态空间方程和公式,实现对模型中各变量进行实时跟踪:The Linear Extended State Observer unit is used to implement real-time tracking of variables in the model using the following state space equations and formulas:
Figure FDA0003886206230000042
其中,x1,x2,x3分别代表所系统的状态变量,
Figure FDA0003886206230000043
b0代表估计的控制增益,w代表外部扰动,y代表所述模型的输出,u代表所述模型的输入;
Figure FDA0003886206230000042
Among them, x 1 , x 2 , x 3 respectively represent the state variables of the system,
Figure FDA0003886206230000043
b 0 represents the estimated control gain, w represents the external disturbance, y represents the output of the model, and u represents the input of the model;
Figure FDA0003886206230000044
其中,z1,z2,z3分别代表所述线性扩张状态观测器的系统状态变量,β123分别代表所述线性扩张状态观测器的增益。
Figure FDA0003886206230000044
Wherein, z 1 , z 2 , z 3 respectively represent the system state variables of the linearly extended state observer, and β 1 , β 2 , β 3 represent the gain of the linearly extended state observer respectively.
7.根据权利要求5所述的一种多旋翼无人机容错控制系统,其特征在于,所述故障模型建立模块包括:7. a kind of multi-rotor unmanned aerial vehicle fault-tolerant control system according to claim 5, is characterized in that, described failure model building module comprises: 检测单元,用于实时检测所述无人机的所有电机;A detection unit is used to detect all motors of the drone in real time; 故障模型建立单元,用于当检测到所述无人机的部分电机异常时,计算故障电机的输出与无故障时的输出的比值,根据所述比值构建故障矩阵RiThe fault model building unit is used to calculate the ratio of the output of the faulty motor to the output when there is no fault when some motors of the drone are detected to be abnormal, and construct a fault matrix R i according to the ratio. 8.根据权利要求5所述的一种多旋翼无人机容错控制系统,其特征在于,所述分配信息获取模块包括:8. A kind of multi-rotor UAV fault-tolerant control system according to claim 5, is characterized in that, described distribution information acquisition module comprises: 分配信息获取单元,用于采用以下公式,得到优化后的分配矩阵Nf,将所述优化后的分配矩阵Nf作为所述无人机上所有电机的控制分配信息:The allocation information acquisition unit is used to obtain the optimized allocation matrix N f by using the following formula, and use the optimized allocation matrix N f as the control allocation information of all motors on the drone: Nf=Af-1Nf = Af -1 ; Nf=AfT(Af·AfT)-1;其中,Af代表部分电机故障后的控制效率矩阵,AfT代表Af的转置;N f = A fT (A f A fT ) -1 ; wherein, A f represents the control efficiency matrix after a partial motor failure, and A fT represents the transposition of A f ; Af=ARi;其中,A代表故障前的控制效率矩阵;A f = AR i ; where, A represents the control efficiency matrix before failure;
Figure FDA0003886206230000051
其中,b为升力系数,l为所述无人机的轴距,d为反扭矩系数。
Figure FDA0003886206230000051
Wherein, b is the lift coefficient, l is the wheelbase of the UAV, and d is the counter torque coefficient.
CN201911322352.8A 2019-12-20 2019-12-20 A fault-tolerant control method and system for a multi-rotor UAV Expired - Fee Related CN110888451B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911322352.8A CN110888451B (en) 2019-12-20 2019-12-20 A fault-tolerant control method and system for a multi-rotor UAV

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911322352.8A CN110888451B (en) 2019-12-20 2019-12-20 A fault-tolerant control method and system for a multi-rotor UAV

Publications (2)

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

Family

ID=69752610

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911322352.8A Expired - Fee Related CN110888451B (en) 2019-12-20 2019-12-20 A fault-tolerant control method and system for a multi-rotor UAV

Country Status (1)

Country Link
CN (1) CN110888451B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111290424B (en) * 2020-03-26 2020-11-06 南方医科大学南方医院 Unmanned aerial vehicle attitude control method for hospital blood sample transportation and unmanned aerial vehicle
WO2021223173A1 (en) 2020-05-07 2021-11-11 深圳市大疆创新科技有限公司 Multi-rotor unmanned aerial vehicle and control method therefor, control device, and computer-readable storage medium
WO2021223136A1 (en) * 2020-05-07 2021-11-11 深圳市大疆创新科技有限公司 Control method, data processing method, device, unmanned aerial vehicle, and storage medium
CN112327896A (en) * 2020-10-29 2021-02-05 东北大学 Rotor fault-tolerant control method and device, computer storage medium and computer equipment
CN112373677B (en) * 2020-11-18 2021-06-22 三生万物(北京)人工智能技术有限公司 Power loss protection system and protection method for six-rotor unmanned aerial vehicle
WO2022193075A1 (en) * 2021-03-15 2022-09-22 深圳市大疆创新科技有限公司 Unmanned aerial vehicle control method, unmanned aerial vehicle, and storage medium
CN112947532B (en) * 2021-04-13 2022-05-03 中国计量大学 Height control method of plant protection quad-rotor unmanned aerial vehicle based on liquid medicine measurement
CN113467244B (en) * 2021-07-15 2024-03-26 之江实验室 Multi-rotor unmanned aerial vehicle fault-tolerant control method and device based on control matrix redistribution
CN115407802B (en) * 2022-10-26 2023-03-24 广东汇天航空航天科技有限公司 Multi-rotor aircraft, control method and device thereof, and computer-readable storage medium
CN115877717B (en) * 2022-12-26 2023-06-13 南通大学 Aircraft fault-tolerant control structure and control method based on active disturbance rejection control
CN115857309A (en) * 2023-02-27 2023-03-28 华东交通大学 Fault-tolerant control method and system for tiltable six-rotor unmanned aerial vehicle

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102902277A (en) * 2012-10-19 2013-01-30 北京航空航天大学 Safe landing method of six-rotor aircraft during one motor shutdown
CN106444812A (en) * 2016-10-26 2017-02-22 华南智能机器人创新研究院 Attitude control method and system based on quad-rotor unmanned helicopter
CN108614573A (en) * 2018-05-15 2018-10-02 上海扩博智能技术有限公司 The automatic fault tolerant attitude control method of six rotor wing unmanned aerial vehicles
CN109308064A (en) * 2017-07-28 2019-02-05 深圳禾苗通信科技有限公司 A kind of the failure tolerant control method and system of quadrotor drone

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102902277A (en) * 2012-10-19 2013-01-30 北京航空航天大学 Safe landing method of six-rotor aircraft during one motor shutdown
CN106444812A (en) * 2016-10-26 2017-02-22 华南智能机器人创新研究院 Attitude control method and system based on quad-rotor unmanned helicopter
CN109308064A (en) * 2017-07-28 2019-02-05 深圳禾苗通信科技有限公司 A kind of the failure tolerant control method and system of quadrotor drone
CN108614573A (en) * 2018-05-15 2018-10-02 上海扩博智能技术有限公司 The automatic fault tolerant attitude control method of six rotor wing unmanned aerial vehicles

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Fault tolerant control Algorithm of hexarotor UAV;Songming Jiao;《Journal of Robotics》;20200804;全文 *
PI-Based Fault Tolerant Control For Fixed-Wing UAVs Using Control Allocation;Jimoh O. Pedro;《IFAC-PapersOnLine》;20171231;全文 *
二阶系统线性自抗扰控制器频带特性与参数配置研究;袁东 等;《控制理论与应用》;20131215;第1630-1640页 *
六旋翼无人机故障容错及精准定高策略研究;刘登攀;《工程科技Ⅱ辑》;20220315;全文 *

Also Published As

Publication number Publication date
CN110888451A (en) 2020-03-17

Similar Documents

Publication Publication Date Title
CN110888451B (en) A fault-tolerant control method and system for a multi-rotor UAV
Labbadi et al. Adaptive fractional-order nonsingular fast terminal sliding mode based robust tracking control of quadrotor UAV with Gaussian random disturbances and uncertainties
Raza et al. Intelligent flight control of an autonomous quadrotor
CN109991991B (en) Robust fault-tolerant tracking method for unmanned helicopter
CN109856972B (en) Robust fault-tolerant tracking control method for unmanned helicopter
CN112684705B (en) Four-rotor aircraft formation tracking control method
Yu et al. Formation control for quadrotor swarm systems: Algorithms and experiments
CN113296507B (en) Multi-power positioning ship cooperative formation control method based on space-time decoupling
CN115639830B (en) An air-ground agent cooperative formation control system and formation control method thereof
CN119088082B (en) A fault-tolerant formation control method for quadrotor drones with anti-interference capability
CN116483124A (en) Anti-interference four-rotor unmanned aerial vehicle formation control method for wireless speed measurement
CN114610072A (en) Distributed time-varying formation tracking control method and system for unmanned aerial vehicle cluster system
CN112882484A (en) Anti-interference fault-tolerant control method for unmanned aerial vehicle
CN117193345A (en) Autonomous landing control method of aircraft-ship collaborative platform based on dynamic event triggering
CN114020038B (en) A multi-quadrotor UAV adaptive sliding mode attitude tracking control method and system
Jiao et al. Fault tolerant control algorithm of hexarotor uav
CN116301041B (en) Method, device and system for controlling formation flight of fixed-wing slave following collar machine
CN114995521B (en) Multi-unmanned aerial vehicle distributed formation control method and device and electronic equipment
CN114545977B (en) An impedance control method for cooperative formation of UAVs with robotic arms
CN117289604A (en) Optimal control system and method for multi-quad-rotor unmanned aerial vehicle based on reinforcement learning
CN116610106A (en) A fault-tolerant tracking control method for multiple UAVs
CN116301026A (en) Large maneuvering agile flight method of four-rotor unmanned aerial vehicle in complex environment
Ban et al. Nonlinear disturbance observer based dynamic surface control for trajectory tracking of a quadrotor uav
Yecheskel et al. Modeling and Control of a Hexacopter with a Rotating Seesaw
Di et al. Tracking Control of the Six-rotor UAV with Input Dead-zone and Saturation

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
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: 20221206