Next Article in Journal
Estimation Error Based Disturbance Observer Design for Flexible Loop Shaping
Next Article in Special Issue
Preliminary Design of an Unmanned Aircraft System for Aircraft General Visual Inspection
Previous Article in Journal
Multiple-Mode Wideband Bandpass Filter Using Split Ring Resonators in a Rectangular Waveguide Cavity
Previous Article in Special Issue
An Efficient SC-FDM Modulation Technique for a UAV Communication Link
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Order Sliding Mode-Based Fixed-Time Active Disturbance Rejection Control for Quadrotor Attitude System

School of Astronautics, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Electronics 2018, 7(12), 357; https://doi.org/10.3390/electronics7120357
Submission received: 4 November 2018 / Revised: 20 November 2018 / Accepted: 21 November 2018 / Published: 26 November 2018
(This article belongs to the Special Issue Autonomous Control of Unmanned Aerial Vehicles)

Abstract

:
This article presents a fixed-time active disturbance rejection control approach for the attitude control problem of quadrotor unmanned aerial vehicle in the presence of dynamic wind, mass eccentricity and an actuator fault. The control scheme applies the feedback linearization technique and enhances the performance of the traditional active disturbance rejection control (ADRC) based on the fixed-time high-order sliding mode method. A switching-type uniformly convergent differentiator is used to improve the extended state observer for estimating and attenuating the lumped disturbance more accurately. A multivariable high-order sliding mode feedback law is derived to achieve fixed time convergence. The timely convergence of the designed extended state observer and the feedback law is proved theoretically. Mathematical simulations with detailed actuator models and real time experiments are performed to demonstrate the robustness and practicability of the proposed control scheme.

Graphical Abstract

1. Introduction

Unmanned aerial vehicles (UAVs) are useful for tasks that are dangerous or unaccessible for human operation and popular in military and civilian applications such as investigation, inspection and surveillance with the advantage of moving in three-dimensional space flexibly. A quadrotor is a kind of vertical takeoff and landing (VTOL) aircraft lifted and controlled by four rotors [1]. Being simpler in structure, less sensitive to damage, easier to handle and more cost-effective, quadrotors have gained more attention in small unmanned aerial system (UAS) research than traditional helicopters. Therefore, plenty of remarkable research achievements related to quadrotors have been made in recent years [2,3].
In practical terms, the position control method is designed according to specific mission requirements and could be implemented onboard or remotely, while the attitude control is performed via an onboard processor to stabilize the attitude reliably [4]. Traditional efficient proportional–derivative (PD) control method has been used in the attitude control of quadrotor [5,6] and achieves good practical performance. Feedback linearization is a commonly used technique to control non-linear coupled systems and are adaptable to the controller design of quadrotor [7]. In [8], the attitude states are transformed into a new state space with a nonlinear transformation to obtain a linear system on which linear error-feedback control can be applied. The idea to design a controller on a transformed linear system instead of the original non-linear system is valuable, however the results in [8] show that the linear error-feedback controller is not robust in terms of uncertainties and measurement noises.
The small low-cost quadrotor is a non-redundant aircraft which suffers from external and internal disturbances such as dynamic wind and actuator faults. In order to design the attitude controller to be robust against unknown disturbances, the observer-based online disturbance estimation and attenuation strategies have been widely studied in recent years. A disturbance observer based on Q-filter is applied in [9,10] for robust hovering control. In [11], a neural network is used to learn the uncertain terms in UAV dynamics. Aboudonia integrated a disturbance observer with feedback linearization-based control in [12] for robustness, and combined a sliding mode controller with a disturbance observer in [13]. The control schemes are verified with the Dryden turbulence model and Von-Karman wind model, respectively. Disturbance observer-based trajectory tracking methods were studied in [14,15]. The sliding mode-based observer and high-order sliding mode-based observer are also extensively studied to enhance the robustness of feedback linearization-based controller [16], sliding mode controller [17] and back-stepping controller [18]. Shi designed extended state observers in [19,20] for the attitude control problem of quadrotors.
Active disturbance rejection control (ADRC) is proposed by Han in 2009 [21]. ADRC inherits from the classical proportional–integral–derivative (PID) controller and executes in an error-driven manner. By augmenting the system uncertainty into the state vector and constructing an extended state observer (ESO), the uncertainty is observed in real-time. Aiming at replacing PID in industrial applications, the ADRC scheme is developed experimentally in the first place. Recently, Guo [22] has proved the convergence of ADRC for non-linear systems mathematically. As a practical control scheme, ADRC has been applied to many engineering aspects, such as motor control, power plant control, ship control, etc [23] (pp. 6–9). ADRC has been adopted for attitude control of a quadrotor directly in [23,24,25,26].
In order to avoid the chattering effect in the sliding mode method, Levant [27,28] designed a general-form arbitrary-order exact robust differentiator based on high-order sliding mode (HOSM) algorithm. The proposed differentiator was proved to converge in finite time with bounded uncertainties. In [29], Levant’s robust differentiator was improved via an extra exponential function with an exponent greater than 1 in order to achieve fixed-time convergence regardless of the initial deviation. In [30], a kind of fixed-time observer is derived and the constraints of the corrective term is formulated. In fact, the observer in [29] is an instance of the formula of the fixed-time observer in [30]. Angulo designed a fixed-time observer operated in a switching manner in [31]. The observer is switched from a uniform differentiator to Levant’s finite-time differentiator. Compared with the methods used in [29,30], Angulo’s switching-type observer performs only one exponential function at each step. Therefore, it is more practical on micro-processors with limited computing power. The switching-type observer has been applied to hypersonic vehicle flight control system [32] and the Brunovsky system [33]. An adaptive super-twisting-based controller is designed for a hypersonic vehicle in [34,35], in which a fixed-time observer is applied. In [36], a fixed-time observer and an integral terminal sliding mode method are studied for fault-tolerant control of a hypersonic vehicle.
In this paper, motivated by the ADRC structure and fixed-time observer, a fixed-time active disturbance rejection control (FTADRC) scheme based on the high-order sliding mode is proposed for the quadrotor attitude system in the presence of unknown disturbances including model uncertainty, dynamic wind, actuator fault and mass eccentricity. The proposed control scheme is analyzed theoretically and verified experimentally. The main contributions of this paper are as follows.
  • The ESO in ADRC is improved via robust uniform high-order sliding mode differentiator to achieve fixed time convergence given bounded differential of lumped disturbance.
  • A non-linear feedback control law combining a high-order sliding mode with feedback linearization is applied in the improved ADRC scheme. In this way, the attitude controller provides fixed-time stability.
The remainder of this paper is organized as follows. The mathematical models are presented in Section 2. The classical ADRC structure is described in Section 3. The proposed high-order sliding mode-based FTADRC is detailed in Section 4. Simulations with a Dryden wind model and experiments with simulated disturbance are carried out in Section 5. Finally, the discussions of the experiments and the conclusions are presented in Section 6 and Section 7.

2. Mathematical Models

2.1. Rigid Body Dynamics

In this section, the quadrotor body is considered as a symmetrical rigid body attached with four sets of actuators. And the center of mass coincides exactly with the center of the body. Figure 1 shows a simplified model of an “X” type quadrotor and the coordinate frames used in this paper. The body coordinate frame Obxbybzb is fixed with the quadrotor body, and the Earth coordinate frame OeNED is fixed with the ground. As the flight time is short, the rotation of the earth is neglected such that the Earth coordinate frame is static in inertial space.
We use modified Rodrigues parameters (MRPs) to represent rigid rotation to avoid the singularity problem of the Euler angles and additional normalization requirement of the unit quaternions. MRPs are expressed in the form of three-dimensional vectors ρ = [ ρ 1 ρ 2 ρ 3 ] T . ρ can be expressed in the form of unit quaternion q = [ q 0 q v T ] T and angle ϕ and axis n of rotation as following.
ρ = q v 1 + q 0 = tan ϕ 4 n
The rotation matrix represented in form of MRP is expressed as follows:
R ( ρ ) = [ 1 8 ( ρ 2 2 + ρ 3 2 ) ( | ρ | 2 + 1 ) 2 4 ( | ρ | 2 ρ 3 2 ρ 1 ρ 2 ρ 3 ) ( | ρ | 2 + 1 ) 2 4 ( | ρ | 2 ρ 2 + 2 ρ 1 ρ 3 ρ 2 ) ( | ρ | 2 + 1 ) 2 4 ( | ρ | 2 ρ 3 + 2 ρ 1 ρ 2 ρ 3 ) ( | ρ | 2 + 1 ) 2 1 8 ( ρ 1 2 + ρ 3 2 ) ( | ρ | 2 + 1 ) 2 4 ( | ρ | 2 ρ 1 2 ρ 2 ρ 3 ρ 1 ) ( | ρ | 2 + 1 ) 2 4 ( | ρ | 2 ρ 2 2 ρ 1 ρ 3 ρ 2 ) ( | ρ | 2 + 1 ) 2 4 ( | ρ | 2 ρ 1 + 2 ρ 2 ρ 3 ρ 1 ) ( | ρ | 2 + 1 ) 2 1 8 ( ρ 1 2 + ρ 2 2 ) ( | ρ | 2 + 1 ) 2 ]
The mathematical model of the quadrotor attitude described by MRPs is defined as follows:
{ ρ ˙ 1 = 1 4 ω x ( ρ 1 2 ρ 2 2 ρ 3 2 + 1 ) + 1 2 ω y ( ρ 1 ρ 2 ρ 3 ) + 1 2 ω z ( ρ 1 ρ 3 + ρ 2 ) ρ ˙ 2 = 1 4 ω y ( ρ 2 2 ρ 1 2 ρ 3 2 + 1 ) + 1 2 ω x ( ρ 1 ρ 2 + ρ 3 ) + 1 2 ω z ( ρ 2 ρ 3 ρ 1 ) ρ ˙ 3 = 1 4 ω z ( ρ 3 2 ρ 1 2 ρ 2 2 + 1 ) + 1 2 ω x ( ρ 1 ρ 3 ρ 2 ) + 1 2 ω y ( ρ 2 ρ 3 + ρ 1 ) ω ˙ x = I y I z I x ω y ω z + τ x I x + J r Ω r ω y I x ω ˙ y = I z I x I y ω x ω z + τ y I y J r Ω r ω x I y ω ˙ z = I x I y I z ω x ω y + τ z I z + J r Ω ˙ r I z
where ω = [ ω x ω y ω z ] T is the angular rate along body coordinate frame. I x , I y and I z are the diagonal elements of inertial matrix. The products of inertia are assumed to be zero. τ = [ τ x τ y τ z ] T is external torque. J r is the rotational inertia of each propeller. Ω r is linear combination of the rotor speeds Ω i ( i = 1 , 2 , 3 , 4 ) and calculated as follows:
Ω r = Ω 1 Ω 2 + Ω 3 Ω 4  
Each MRP ρ has a shadow MRP ρ S that represents the same orientation as ρ . Using the combination of MRPs and their shadow set, any rotation could be described without singularity. It is easy to be understand that rotating angle ϕ around the axis n is exactly the same as rotating angle 360 ϕ around the axis n , as shown in Figure 2. According to [37], the shadow point of MRP ρ is formulated as follows.
ρ S = ρ ρ 2
It is worth noting that ρ is also the shadow point of ρ S .

2.2. Actuator Mode

Each of the four sets of actuators consists of an electron speed regulator, a brushless direct-current (DC) motor (BLDM) and a propeller. The electron speed regulator is responsive to input signals. Thus, it is modelled as a simple proportional component and neglected in following sections.

2.2.1. Motor Model

The circuit of the BLDM can be simplified into a series circuit composed of an equivalent resistance R M , an equivalent inductance L M and the induction power supply e M produced by rotor rotation. The input voltage u satisfies the following voltage equation of the motor circuit:
u = R M i + L M d i d t + Ω K v
where i is the current, Ω is the rotational speed of the rotor, K v is the motor velocity constant. The motor torque is proportional to the current with coefficient K T .
Q M = K T i
Because the propeller is mechanically attached with a motor rotor, we ignore the rotational inertia of the motor rotor and treat it as part of the rotational inertia of the propeller J r . The rotation dynamics equation of the propeller is formulated as follows:
J r Ω ˙ = Q M Q
where Q represents the aerodynamic drag torque, which will be discussed later in the propeller aerodynamic model. The equivalent inductance L M of small BLDM is negligible. Therefore, the motor is modelled by the following first-order differential equation with input u and output Ω .
J r Ω ˙ = K T [ 1 R M ( u Ω K v ) i 0 ] Q

2.2.2. Propeller Aerodynamic Model

We build the propeller aerodynamic model for quadrotors mainly by consulting the aerodynamics for a helicopter rotor in [38]. The model combines momentum theory and blade element analysis, and achieves the following equation of the inflow ratio λ :
4 F λ 2 + μ 2 ( λ λ c ) r 1 4 σ C l α ( θ μ 2 + 2 θ r 2 2 λ r ) = 0
where μ is the rotor advance ratio; F is a function of λ named by Prandtl correction factor; λ c is the climb inflow ratio; σ is the rotor solidity; C l α is the lift-curve slope; θ is the blade pitch angle; r is the non-dimensional radial distance. The readers can refer to [38] for the detailed description and derivation of the above equation.
λ cannot be directly solved from the above equation as F is a non-linear function of λ . Instead, the equation could be solved in a nested iterative way as following.
  • Set F = 1 .
  • Solve the equation through the Newton method with initial value in the case of μ = 0 :
    λ = λ 0 ( r , λ c ) = ( σ C l α 16 F λ c 2 ) 2 + σ C l α 8 F θ r ( σ C l α 16 F λ c 2 )
  • Calculate F ( λ ) , go to step 2 and start the next iteration.
Generally, convergence is obtained after about 3 times of iteration (outer cycle).
Divide the blade into N elements with length of Δ r . The rotor thrust coefficient C T and torque coefficient C Q are approximated with the sum of corresponding coefficients of each annulus with the width of Δ r as shown in Figure 3.
C T = n = 1 N Δ C T n C Q = n = 1 N Δ C Q n Δ C T n = 4 F ( r n ) λ ( r n ) 2 + μ 2 ( λ ( r n ) λ c ) r n Δ r Δ C Q n = λ ( r n ) Δ C T n + 1 4 σ ( r n ) C d ( 2 r n 2 + μ 2 ) Δ r
where C d is the drag coefficient, n = 1 , , N . The trust and torque generated by the rotor could be calculated with the coefficients.
{ T = C T ρ π R 2 ( Ω R ) 2 Q = C Q ρ π R 3 ( Ω R ) 2
where ρ is the air density; R is the radius of the rotor disk.

2.3. Example of Measuring and Calculating Propeller Aerodynamic Model Parameters

Multiple parameters need to be determined before applying the above propeller aerodynamic model in simulations. We summarize the parameters in Table 1, and present the method for determination with an example.
As Table 1 shows, there are four kinds of parameters required in propeller aerodynamic model. The atmospheric parameter is directly determined with the Earth’s atmospheric parameters. Propeller parameters are manually measured. Velocity parameters are calculated online according to wind speed and quadrotor motion in simulations. Aerodynamic coefficients are usually determined with the aid of an accurate structural model, which is difficult to construct. We use an experimental method to estimate the coefficients through measuring rotor trust, torque and rotation speed.
We use 16″ carbon propellers as example for measuring and estimating. The propeller radius is 203.2 mm. The number of blades per propeller is 2. Blade chord and blade pitch angle vary with the radial position. We select 12 measurement points on the blade and measure the chord and pitch angle manually on the measurement points. The measured results are shown in Table 2.
We use a tachometer (Hobbywing, Shenzhen, China), a tensionmeter (Jnsensor, Bengbu, China) and a torquemeter (Jnsensor, Bengbu, China) to measure the thrust-speed and torque-speed curves of 16″ carbon propellers in static air. The three sensors are shown in Figure 4. In the measurement process, the tensionmeter and the torquemeter are fixed on a solid stable structure, while the motor with the propeller is mounted on the tensionmeter or the torquemeter with a flange. The tachometer is connected with arbitrary two phases of the three-phase motor and measures the rotation speed by counting the phase changes. The tensionmeter and the torquemeter output analog signals. The signals are converted into digital measurements with an analog-to-digital (A/D) converter and sent to the computer through a serial port.
According to Equation (13), the rotor thrust and torque increase with the square of the rotation speed with ratio of C T ρ π R 4 and C Q ρ π R 5 respectively. We take measurements in static air and fit the thrust-speed and torque-speed curves with 2-order polynomial function. The thrust coefficient and torque coefficient are recovered from the fitting parameters. We merge the measurements from four sets of motors and propellers to get the average coefficients. The resulting curves are shown in Figure 5. It is shown that the resulting 2-order polynomial functions fit the measurements well.
According to the fitting parameters, C T and C Q are 1.36 × 10−2 and 2.028 × 10−3 in the case of λ c = 0 and μ = 0 . The whole calculation process of C T described before is treated as a non-linear function of C l α , i.e.,
C T = f C T ( C l α )
It can be easily solved by the classical dichotomy method. With the above results, we obtain C l α = 23.91 . The solution of C d is 0.0384 based on Equation (12) of C Q .

3. Active Disturbance Rejection Control (ADRC) Method

The ADRC method applied to a second-order system is shown in Figure 6. The method consists of three parts: tracking differentiator (TD) is used to generate the desired states given the desired outputs; extended state observer (ESO) is used to estimate the system uncertainty; non-linear state error feedback (NLSEF) is used to feedback state errors effectively.
The traditional ESO is a model-independent disturbance observer and can be employed universally in the non-linear system control problem. However, the convergence rate of the traditional ESO is not sufficient for a quadrotor attitude control problem. Therefore, this paper focuses on improving the ESO with a high-order slide mode algorithm for quadrotor attitude control problem. Additionally, the non-linear state error feedback control law is improved with a multivariable high-order slide mode algorithm.

4. High-Order Sliding Mode-Based Fixed-Time Active Disturbance Rejection Control (FTADRC)

The FTADRC scheme is designed within the traditional ADRC structure based on the feedback-linearization technique. A schematic block diagram of the proposed FTADRC scheme for quadrotor attitude control is shown in Figure 7. The attitude control scheme consists of five parts:
  • Feedback linearization for regularizing the attitude dynamic model;
  • Fixed-time extended state observer (FTESO) for observing the unknown disturbances accurately;
  • MRP-TD for tracking the differential of input attitude described by MRP;
  • Non-linear feedback control law for driving the orientation of quadrotor to track the desired attitude timely;
  • Control allocation for generating pulse-width modulation (PWM) signals for motors.

4.1. Feedback Linearization

We remove the gyroscopic effect and reactionary torque related to the rotor speed in the attitude dynamic equations and treat them as part of the unknown disturbance. A six-dimensional system is defined with state vector x = [ ρ 1 , ρ 2 , ρ 3 , ω x , ω y , ω z ] T , input vector u = [ τ x , τ y , τ z ] T and measurable output vector y = [ y 1 , y 2 , y 3 ] T . A non-linear system is achieved according to the attitude model in Equation (3).
[ ρ ˙ 1 ρ ˙ 2 ρ ˙ 3 ω ˙ x ω ˙ y ω ˙ z ] x ˙ = [ 1 4 ω x ( ρ 1 2 ρ 2 2 ρ 3 2 + 1 ) + 1 2 ω y ( ρ 1 ρ 2 ρ 3 ) + 1 2 ω z ( ρ 1 ρ 3 + ρ 2 ) 1 4 ω y ( ρ 2 2 ρ 1 2 ρ 3 2 + 1 ) + 1 2 ω x ( ρ 1 ρ 2 + ρ 3 ) + 1 2 ω z ( ρ 2 ρ 3 ρ 1 ) 1 4 ω z ( ρ 3 2 ρ 1 2 ρ 2 2 + 1 ) + 1 2 ω x ( ρ 1 ρ 3 ρ 2 ) + 1 2 ω y ( ρ 2 ρ 3 + ρ 1 ) ( I y I z ) ω y ω z / I x ( I z I x ) ω x ω z / I y ( I x I y ) ω x ω y / I z ] f ( x ) + [ 0 0 0 0 0 0 0 0 0 1 I x 0 0 0 1 I y 0 0 0 1 I z ] g [ τ x τ y τ z ] u [ y 1 y 2 y 3 ] y = [ ρ 1 ρ 2 ρ 3 ] h ( x )
The Lie derivative and k-th Lie derivative of function h i ( x ) with respect to a vector-valued function f ( x ) are defined as follows:
L f h i ( x ) = j = 1 6 h i ( x ) x j f j ( x ) L f k h i ( x ) = L f ( L f k 1 h i ( x ) )
A 3 × 3 matrix E ( x ) is constructed with the Lie derivatives L g j L f h i ( x ) , ( i [ 1 , 3 ] , j [ 1 , 3 ] ) :
E ( x ) = [ L g 1 L f h 1 ( x ) L g 2 L f h 1 ( x ) L g 3 L f h 1 ( x ) L g 1 L f h 2 ( x ) L g 2 L f h 2 ( x ) L g 3 L f h 2 ( x ) L g 1 L f h 3 ( x ) L g 2 L f h 3 ( x ) L g 3 L f h 3 ( x ) ]
Substituting functions of (15) into E ( x ) yields:
E ( x ) = [ ρ 1 2 ρ 2 2 ρ 3 2 + 1 4 I x ρ 1 ρ 2 ρ 3 2 I y ρ 1 ρ 3 + ρ 2 2 I z ρ 1 ρ 2 + ρ 3 2 I x ρ 2 2 ρ 1 2 ρ 3 2 + 1 4 I y ρ 2 ρ 3 ρ 1 2 I z ρ 1 ρ 3 ρ 2 2 I x ρ 2 ρ 3 + ρ 1 2 I y ρ 3 2 ρ 1 2 ρ 2 2 + 1 4 I z ]
with a determinant of:
| E ( x ) | = ( ρ 1 2 + ρ 2 2 + ρ 3 2 + 1 ) 3 64 I x I y I z
Obviously | E ( x ) | > 0 holds for all ρ 3 . Therefore, E ( x ) is non-singular. The relative degree vector of system (15) is r = [ 2 , 2 , 2 ] T . System (15) can be transformed into a regular form.
[ y ¨ 1 y ¨ 2 y ¨ 3 ] = [ L f 2 h 1 ( x ) L f 2 h 2 ( x ) L f 2 h 3 ( x ) ] Ф ( x ) + E ( x ) u
The formula of a three-dimensional vector Ф ( x ) is in Appendix A. Defining a new six-dimensional state vector ξ = [ ξ 1 T , ξ 2 T ] T = [ ρ 1 , ρ 2 , ρ 3 , ρ ˙ 1 , ρ ˙ 2 , ρ ˙ 3 ] T , a new state equation is derived from (20).
{ ξ ˙ 1 = ξ 2 ξ ˙ 2 = Ф + E u
where Ф and E are determined by the system model parameters. We divide them into the nominal parts Ф ¯ , E ¯ and the uncertain parts Δ Ф , Δ E . The input torque u is divided into the control torque u c and the unknown disturbance torque Δ u . Ф , E and u are rewritten in the following format.
{ Ф = Ф ¯ + Δ Ф E = E ¯ + Δ E u = u c + Δ u
Substituting (22) into (21) yields:
{ ξ ˙ 1 = ξ 2 ξ ˙ 2 = ( Ф ¯ + E ¯ u c ) + ( Δ Ф + Δ E u c + E Δ u )
Define an auxiliary input vector μ and a lumped disturbance vector Δ μ as follows:
{ μ = Ф ¯ + E ¯ u c Δ μ = Δ Ф + Δ E u c + E Δ u
Substituting (24) into (23) yields:
{ ξ ˙ 1 = ξ 2 ξ ˙ 2 = μ + Δ μ

4.2. Fixed-Time Third-Order Sliding Mode Observer-Based Extended State Observer (ESO)

According to (24), the equivalent disturbance torque Δ τ is related to Δ μ in the format:
Δ τ = E ¯ 1 Δ μ
Since E ¯ is a non-singular matrix, the estimation of Δ τ is equivalent to the estimation of Δ μ . Augment Δ μ as an extended state ξ 3 = Δ μ , the reconstructed linear system is rewritten as:
{ ξ ˙ 1 = ξ 2 ξ ˙ 2 = ξ 3 + μ ξ ˙ 3 = ζ
where ζ is the derivative of ξ 3 . Assume that each element of ζ is bounded.
ζ i [ L , L ] , i = 1 , 2 , 3
The elements in a 3-dimensional vector ξ 1 , ξ 2 and ξ 3 are independent of each other. Motivated by [31], a fixed-time convergent extended state observer (FTESO) for each element i (I = 1,2,3) is designed as:
ξ ^ ˙ 1 , i = θ i α 1 , i ξ ^ 1 , i ξ 1 , i δ 2 3 ( 1 θ i ) β 1 , i ξ ^ 1 , i ξ 1 , i δ 1 + ε + ξ ^ 2 , i ξ ^ ˙ 2 , i = θ i α 2 , i δ ξ ^ 1 , i ξ 1 , i δ 1 3 ( 1 θ i ) β 2 , i δ ξ ^ 1 , i ξ 1 , i δ 1 + 2 ε + ξ ^ 3 , i + μ i ξ ^ ˙ 3 , i = θ i α 3 , i δ 2 ξ ^ 1 , i ξ 1 , i δ 0 ( 1 θ i ) β 3 , i δ 2 ξ ^ 1 , i ξ 1 , i δ 1 + 3 ε θ i = { 0 , ξ ^ 1 , i ξ 1 , i > δ 1 , ξ ^ 1 , i ξ 1 , i δ
where ξ ^ 1 , i , ξ ^ 2 , i and ξ ^ 3 , i are the estimation of ξ 1 , i . ξ 2 , i and ξ 3 , i respectively. δ > 0 is used to scale the estimation error. Function x k = | x | k sign ( x ) , sign ( · ) represents the signum function. ε > 0 is chosen small enough. θ i is used to switch between two different exponential functions. α j , i ( j = 1 , 2 , 3 ) are selected based on the boundary value L using the formulas for the HOSM differentiator in [28]. β j , i ( j = 1 , 2 , 3 ) are selected such that the following matrix is Hurwitz as suggested by [31].
A = [ β 1 , i 1 0 β 2 , i 0 1 β 3 , i 0 0 ]
There are two major improvements between the above observer with the differentiators used in [31,32].
1. The switching is conducted according to the estimation error instead of an arbitrary fixed time.
2. The estimation error ξ ^ 1 , i ξ 1 , i is scaled with 1 / δ .
The switching structure of the observer takes advantage of dealing with different estimation errors with different exponential functions. The function with exponential larger than 1 is more efficient with larger errors, while the function with exponential smaller than 1 is more efficient with smaller errors. So, we use a switch mechanism controlled by estimation errors to make full use of different exponential functions. The scaling factor is used to adjust flexibly the estimation errors for better performance.
Theorem 1.
Suppose that ξ 1 and μ are available in real time. FTESO (29) converges to the true extended states ξ 1 , i , ξ 2 , i and ξ 3 , i in a fixed time with the parameters chosen according to preceding rules.
Proof of Theorem 1.
Define the following notations for observation errors.
{ Δ ξ 1 , i = ξ ^ 1 , i ξ 1 , i Δ ξ 2 , i = ξ ^ 2 , i ξ 2 , i Δ ξ 3 , i = ξ ^ 3 , i ξ 3 , i
Differentiating Δ ξ 1 , i , Δ ξ 2 , i and Δ ξ 3 , i , one can obtain:
{ Δ ξ ˙ 1 , i = θ i α 1 , i ξ ^ 1 , i ξ 1 , i δ 2 3 ( 1 θ i ) β 1 , i ξ ^ 1 , i ξ 1 , i δ 1 + ε + Δ ξ ^ 2 , i Δ ξ ˙ 2 , i = θ i α 2 , i δ ξ ^ 1 , i ξ 1 , i δ 1 3 ( 1 θ i ) β 2 , i δ ξ ^ 1 , i ξ 1 , i δ 1 + 2 ε + Δ ξ ^ 3 , i Δ ξ ˙ 3 , i = θ i α 3 , i δ 2 ξ ^ 1 , i ξ 1 , i δ 0 ( 1 θ i ) β 3 , i δ 2 ξ ^ 1 , i ξ 1 , i δ 1 + 3 ε ζ i
Define:
{ η 1 , i = Δ ξ 1 , i ( δ t ) / δ η 2 , i = Δ ξ 2 , i ( δ t ) η 3 , i = δ Δ ξ 3 , i ( δ t )
Substituting (33) into (32) yields:
{ η ˙ 1 , i = θ i α 1 , i η 1 , i 2 3 ( 1 θ i ) β 1 , i η 1 , i 1 + ε + Δ ξ ^ 2 , i η ˙ 2 , i = θ i α 2 , i η 1 , i 1 3 ( 1 θ i ) β 2 , i ξ 1 , i 1 + 2 ε + Δ ξ ^ 3 , i ξ ˙ 3 , i = θ i α 3 , i ξ 1 , i 0 ( 1 θ i ) β 3 , i ξ 1 , i 1 + 3 ε δ 2 ζ i
According to Theorem 1 in [31], (34) converges to a 3-sliding surface S = { η 1 , i = 0 , η 2 , i = 0 , η 3 , i = 0 } in a fixed time T . Therefore, (29) converges to the true extended states ξ 1 , i , ξ 2 , i and ξ 3 , i in a fixed time T / δ . This complete the proof. □

4.3. Tracking Differentiator

The tracking differentiator for the desired attitude ρ d = [ ρ 1 , d , ρ 2 , d , ρ 3 , d ] T based on the above feedback linearization model is obtained by directly applying the TD in the classical ADRC structure.
{ ξ ˙ 1 , d = ξ 2 , d ξ ˙ 2 , d = f h a n 3 ( ξ 1 , d ρ d , ξ 2 , d , r , h 0 )
where function f h a n 3 is the vector version of tracking function fhan in [21]. It performs function fhan for each element of the 3-D vectors. The above tracking differentiator outputs the desired linearized states ξ 1 , d and ξ 2 , d . Considering the multi-solution characteristic of MRP in describing rotation, the shadow tracking differentiator for ρ d S should also be constructed as follows.
{ ξ ˙ 1 , d S = ξ 2 , d S ξ ˙ 2 , d S = f h a n 3 ( ξ 1 , d S ρ d S , ξ 2 , d S , r , h 0 )

4.4. Multivariable High-Order Sliding Mode (HOSM)-Based Fixed-Time Non-Linear Feedback Law

Define the sliding surface and shadow sliding surface:
σ = C ( ξ 1 ξ 1 , d ) + ( ξ 2 ξ 2 , d ) σ S = C ( ξ 1 ξ 1 , d S ) + ( ξ 2 ξ 2 , d S )
where C is a 3 × 3 diagonal matrix with positive elements of main diagonal. For tracking the desired attitude described by the MRP set { ρ d , ρ d S } , there are two candidate sliding mode surfaces to approach. We force the quadrotor moving towards the “nearby candidate” in { σ , σ S } with minimal 2-norm. Without loss of generality, the fixed-time feedback law is designed and analyzed for σ in the following.
Differentiating σ and substituting (25) into the expression of σ ˙ , one can obtain:
σ ˙ = C ( ξ 2 ξ 2 , d ) + μ + Δ μ
Inspired by [39], we use 2-norm of the 3-dimension vector σ to apply the high-order sliding mode algorithm on the above multivariable system. A fixed-time second-order sliding mode control law with disturbance rejection is designed as:
μ = k 1 σ σ p 1 k 2 σ σ p 2 C ( ξ 2 ξ 2 , d ) ξ ^ 3
where k 1 > 0 , k 2 > 0 , 0 < p 1 < 1 , 0 < p 2 < 1 , and the notation σ represents 2-norm of vector σ . According to Theorem 1, ξ ^ 3 becomes equal to Δ μ in a fixed time. Substituting (39) into (38), the close-loop dynamic equation of σ can be obtained.
σ ˙ = k 1 σ σ p 1 k 2 σ σ p 2
Theorem 2.
Consider the closed-loop system (40), the sliding variable σ and its derivative σ ˙ reach the origin in fixed time.
Proof of Theorem 2.
Given (40), the 3-dimension vectors σ ˙ and σ are parallel in opposite directions. That means the direction of σ does not change with time, and the following equation holds.
σ ˙ = d σ d t
Substituting (40) into (41) yields:
d σ d t = ( k 1 σ p 1 + k 2 σ p 2 ) σ = k 1 σ 1 p 1 k 2 σ 1 + p 2
When the initial value of the system satisfies the inequation σ ( 0 ) > 1 , according to the above equation, one can obtain:
d σ d t < k 2 σ 1 + p 2
Inequation (43) can be rewritten as:
d σ σ 1 + p 2 < k 2 d t
Denote the time consumed for σ ( t ) reducing from σ ( 0 ) to 1 to T 1 , i.e., σ ( T 1 ) = 1 . Integrate Equation (44), and it can be found that T 1 is bounded by a fixed value that is independent of initial value σ ( 0 ) .
T 1 < 1 k 2 p 2
After the system state reaching σ = 1 in time T 1 , according to (42), one can obtain:
d σ d t < k 1 σ 1 p 1
Denote the time consumed for σ ( t ) reducing from 1 to 0 to T 2 , i.e., σ ( T 1 + T 2 ) = 0 . Integrate Equation (46), and it can be concluded that if σ ( T 1 ) is in range (0,1], T 2 is bounded by a fixed value.
T 2 < 1 k 1 p 1
For a 3-dimensional vector σ , σ = 0 σ = 0 . According to Equation (40), σ ˙ = 0 when σ = 0 . This completes the proof. □
The control torque is calculated according to the first equation in (24).
u c = E ¯ 1 μ Ф ¯
The feedback law (39) performs exponential functions on the 2-norm of σ instead of three Euler angles used in [15]. In this way, the exponential functions-related computation is reduced. This characteristic is helpful when using a low-cost advanced RISC (reduced instruction set computing) machine (ARM) microprocessor that is ineffective to nonlinear functions such as exponential, trigonometric and inverse trigonometric functions.
By means of input–output feedback linearization, the non-linear attitude system is transformed into a linear system and the parameters in the observer and controller are adjusted in a standard way [28]. More importantly, the matrix B for tuning the disturbance estimation and compensation in traditional ADRC [21] is a model-related parameter even though the whole method is independent of the type of system function. This means the parameter is hard to tune in attitude control of quadrotor that the control torque is not generated directly by the controller output. The HOSM-based FTESO and non-linear feedback law are constructed based on the equivalent linear model, such that the observer and the controller are completely independent without common parameters and can be adjusted individually. In conclusion, the proposed FTADRC method is tuned more easily than the traditional ADRC method.

4.5. Non-Linear Control Allocation

The control torque is allocated to four actuators according to the actuator models. After the motors reach the steady state, the PWM signals outputted by the controller are proportional to the rotor speeds. Given the quadratic dependence of thrust and torque upon rotor speed, a simple nonlinear control allocation method is used in this paper.
{ r 1 = 1 2 2 τ x 2 l λ T + 2 τ y 2 l λ T + τ z λ Q + F λ T r 2 = 1 2 2 τ x 2 l λ T + 2 τ y 2 l λ T τ z λ Q + F λ T r 3 = 1 2 2 τ x 2 l λ T 2 τ y 2 l λ T + τ z λ Q + F λ T r 4 = 1 2 2 τ x 2 l λ T 2 τ y 2 l λ T τ z λ Q + F λ T  
where r i , ( i = 1 , 2 , 3 , 4 ) represents the duty circle of PWM signal for each actuator, l is the distance between the rotor and the center of mass of the quadcopter, F is the whole thrust provided by remote controller, λ T and λ Q are the coefficients related to actuator models which can be approximately calculated with maximum thrust T max , maximum torque Q max and maximum duty circle r max .
λ T = T max r max 2 λ Q = Q max r max 2

5. Simulation and Experimental Results

5.1. Simulation Results

A simulation model for quadrotor attitude control is built according to the mathematical models and data acquired in Section 2. The structure of the simulation model is shown in Figure 8. The controller module includes the control method to be verified. The output signals are passed into a zero-order holder before the mathematical models in order to simulate the discrete calculation on the digital processor. The model part includes the motor model, propeller aerodynamic model, and rigid body dynamic model.
The unknown wind is added into the model as external disturbance. The rotor thrust and torque are influenced by the wind speed near the propeller. As the quadrotor body is acting as an obstruction to the wind, the wind attenuates while it passes through the body. As a result, the wind speed at the upwind side is larger than the wind speed at the downwind side. The wind attenuation is related to the interaction between the quadrotor body and the wind. In this paper, the wind speed is resolved into a component w parallel to the rotor disk and a component u perpendicular to the rotor disk. The attenuation of u is neglected so that u is constant near the quadrotor body. w decreases proportionally with the distance traveled along the quadrotor body. Figure 9 is a graphical representation of w decreasing with its travel along the quadrotor body. The green dashed line is the component w passing the body centroid. p i ( i = 1 , 2 , 3 , 4 ) represents four rotor centers. After projecting p i on the gray line, the projection point closest to the wind direction is selected as the windward point p w . The wind speed near the rotor corresponding to p w is supposed to be w . The wind speed near other three rotors is calculated as follows:
w i = max ( 1 γ w d i , 0 ) w
where γ w is the attenuation factor, d i is the distance between rotor i with p w along vector w .
The above wind attenuation leads to different thrust at the upwind side and downwind side and in turn generates a wind-related overturning torque which acts as the disturbance in attitude control system.
The commonly used Dryden turbulence model is used to generate random wind. The fundamental of the Dryden model is constructing a transfer function based on spectrum functions of atmospheric turbulence velocity and its gradient to convert the white noise signal to a colored noise signal. The Dryden module in aerospace toolbox of Simulink (Version 8.9, MathWorks., Natick, MA, United States, 2017) is used in this paper to generate random wind speed. The low-altitude intensity in the model is set to 10m/s. The other parameters are decided according to MIL-F-8785C standard. The wind speed along three axes of Earth coordinates framed with an altitude of 6m is shown in Figure 10.
We present the mathematical simulations of proposed FTADRC attitude control method on a quadrotor model with the parameters in Table 3. The other parameters related to propeller shape and aerodynamics are chosen as Section 2.3. The first four parameters in Table 3 are used in the controller. Practically, l is easy to measure, while the inertias are hard to obtain. So, the nominal parameters are selected as the last column.
In the simulations, the position of the quadrotor is assumed to be free from the forces that include the propeller thrust and gravity. A square wave signal with amplitude of 20°, period of 8s is used as the desired roll angle input to analyze the step response. The random wind generated by the Dryden model is added into the simulation since the 15th second in order to compare the proposed FTADRC with the traditional ADRC in the absence and presence of dynamic wind disturbance. The simulation results of two different attitude control methods are shown in Figure 11, Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16.
As illustrated in Figure 11, quick convergence of the roll angle is achieved by both methods before 15 s with no disturbances because we tune the parameters of the two different controllers to achieve high-gain feedback control. With the presence of simulated dynamic wind, FTADRC tracks the desired attitude more accurately than traditional ADRC. Figure 14, Figure 15 and Figure 16 indicate that the disturbance estimation of FTESO is more precise than traditional ESO.
The sharp peaks in Figure 14, Figure 15 and Figure 16 are caused by motor inertia which is generally neglected in previous researches. The motors are not able to provide the rapidly changing control torque to follow the step input. The deviation of control torque acts as a rapidly changing disturbance which is difficult to track for both FTESO and ESO. However, the sudden disturbance vanishes quickly as the motor tracks the desired speed, and the attitude of the quadrotor is not seriously affected.

5.2. Experimental Results

In order to verify the effectiveness of the proposed FTADRC attitude control method in practical applications, we have developed a flight control unit (FCU) mainly using an STM32F103 micro-processor (STMicroelectronics, Geneva, Switzerland) and two ICM-20689 inertial measurement units (IMUs, InvenSense, San Jose, CA, United States). The angular speed and acceleration information provided by the two IMUs is averaged to generate more accurate measurements. An external HMC5983 magnetometer (Honeywell, Morris Plains, NJ, United States) is used for assisting the attitude determination process. The attitude determination algorithm in [40] is applied to the FCU providing attitude measurements. The core components including micro-processor and sensors cost less than 40 US dollars. In the experiments, the proposed fixed-time disturbance rejection control method consumes less than 1 ms on the 72 MHz ARM processor.
The experiments are conducted with a self-assembled QR450 quadrotor. The experimental setup consists of a remote controller, a ground control station (GCS) and the quadrotor with FCU onboard is shown in Figure 17. A pair of Xbee modules is used for communication between the FCU and GCS.
The quadrotor frame is made by aluminium alloy and carbon fiber to lower costs, and it has been deformed after being used for more than a year in plenty of flight tests, as shown in Figure 18. The deformations result in additional disturbances as the parameters are set ignoring them. The parameters are chosen as Table 4, in which I x , I y and I z are estimated intuitively without precise measurement.
The experiments are conducted indoors without Global Positioning System (GPS) signals. To assure flight security, the quadrotor is controlled remotely. The external unknown disturbance is simulated in three different ways. The performance of rejecting disturbance is verified respectively in all cases.
Case 1.Eccentric mass
In this case, an iron mass block of 194 g acting as the eccentric mass is hung from the right-hand front corner of the quadrotor frame, as shown in Figure 19. The original quadrotor without the eccentric mass is 1195 g. The eccentric mass is equivalent to 16% of the original quadrotor and 64% of individual motor thrust during hovering in weight. In addition, the block waves during flight and generates time-varying disturbance. The iron block hung from a corner of the frame mainly acts as an unknown constant torque in the attitude control system. In a simple error-driven control method, the constant disturbance torque will cause an attitude bias.
The curves of the attitude angles are shown in Figure 20, Figure 21 and Figure 22, and the curves of the estimated disturbances are shown in Figure 23. The roll and pitch angles are affected by the eccentric mass and the control errors increase to about 2° in 1.5 s after the quadrotor takes off from the ground. The FTESO tracks the unknown disturbances in about 3 s with an initial value of 0. The fast-changing desired roll and pitch angles are produced by remote controller to test the controller’s performance. By estimating and attenuating the estimated disturbances, the attitude angles are tracked accurately.
Case 2.Sudden fault of a single motor
In this case, we simulate the actuator fault of motor 1 (right hand front motor) with a programed effectiveness loss of 30%, which means the output duty circle of motor 1 is multiplied by 0.7. The time of fault and recovery is controlled by a digital switch on the remote controller. As a non-redundant system, the quadrotor is sensitive to motor faults. The attitude model is uncontrollable if any of the four motors is completely disabled. So, we employ a programed motor effectiveness loss instead of a complete motor failure to model an in-flight actuator fault.
In the experiment, the fault occurs at 6.4 s and the motor recovers at 32.4 s. The curves of attitude angles are shown in Figure 24, Figure 25 and Figure 26, and the curves of estimated disturbances are shown in Figure 27. The roll angle quickly changes to 19.8°, and the pitch angle quickly changes to −14.8° as a result of the motor fault. FTESO converges to the disturbances caused by the fault and the attitude angles track the desired values after 16 s.
The recovery of the motor is also treated as an external disturbance. It can be easily known from the motor model that the acceleration of the motor is slower than the deceleration. Thus, the recovery process has less impact on the quadrotor with a bounded disturbance differential-based FTESO. The roll angle changes to −16.4°, and the pitch angle changes to 14.4° 0.4 s after the motor recovers from the fault. The disturbance estimations of FTESO approximately turn back to the values before the motor fault 5.6 s after the recovery.
Case 3.Wind and damaged propeller
In this case, we install a damaged propeller on motor 2 (left hand front motor) and use two electrical fans to generate the external wind as shown in Figure 28. We make the quadrotor hover in the wind by remote control. The average wind speed generated by the two electrical fans is around 4.5 m/s. The main purpose of this case is to show the performance of FTADRC method with common lumped disturbances in quadrotor attitude control such as wind disturbance and propeller damage. Figure 29, Figure 30, Figure 31 and Figure 32 present the experimental results of the FTADRC and conventional ADRC schemes. In each figure, panel (a) indicates the result of FTADRC, and panel (b) indicates the result of conventional ADRC.
Figure 29, Figure 30 and Figure 31 present the tracking errors of attitude angles of FTADRC and conventional ADRC control schemes. The tracking errors of roll and pitch angles using FTADRC control scheme keep below 2° with average values around 0°, while the tracking errors of roll and pitch angles using the ADRC control scheme have peak-to-peak values larger than 15° and obviously non-zero average values. Figure 32 shows disturbance estimation with FTADRC and traditional ADRC schemes. The disturbance introduced by the wind is hard to evaluate accurately because of the complicated wind dynamics related to the electrical fans. However, it still could be seen from Figure 32 that the disturbance estimation using FTADRC changes faster than that using conventional ADRC. Along with the tracking errors of attitude angles, it can be demonstrated that FTADRC responds more quickly to the external disturbances and tracks the input attitude more accurately than conventional ADRC.

6. Discussion

In the experiment with eccentric mass, the quadrotor takes off from the ground slowly so that the eccentric mass could be regard as a disturbance accumulating during take-off. Thus, its influence on attitude angles is not obvious and FTESO converges to the lumped disturbances quickly. During the flight, the attitude of the quadrotor accurately tracks the desired attitude with disturbances rejected by FTESO. This experiment together with the simulations demonstrate the stability and robustness of FTADRC in the presence of slowly varying disturbance such as dynamic wind and eccentric mass.
In the experiment with a motor fault, the disturbance caused by the fault takes effect rapidly. Strictly speaking, the differential of disturbance exceeds the bound set in FTESO, and the proposed FTADRC is not suitable for such a case theoretically. However, the results show that the attitude system of quadrotor is stabled by FTADRC except that FTESO consumes more time for convergence. Therefore, the proposed robust control scheme could be applied to the practical situations beyond the theoretical constraints.
Since the motor merely fails in flight and the centroid of mass is generally adjusted before flight, the conditions in the previous two cases are rigorous, and not normal in the practical operation of quadrotors. The last case with damaged propeller and unknown wind represents a more common operation situation of a low-cost quadrotor. In this experiment, the proposed FTADRC and conventional ADRC attitude control schemes are compared. The performance of FTADRC clearly exceeds the conventional ADRC according to the results.

7. Conclusions

In this paper, the FTADRC scheme is proposed based on a HOSM method for the attitude control of a quadrotor with unknown disturbance. The control scheme employs a fixed-time ESO for observing the lumped disturbance, and uses a multivariable fixed-time second-order sliding mode method to fast approach the sliding surface for a feedback-linearized system. Detailed mathematical models are built for simulation analysis. According to the comparative simulation and experiment results, the attitude control accuracy is improved compared with the traditional ADRC because FTESO achieves higher precision than the traditional ESO, and the control scheme is effective on a low-cost flight control unit. It can be concluded that the proposed FTADRC is robust toward dynamic wind, mass eccentricity, and motor fault. Additionally, FTADRC is practical for a low-cost quadrotor because of its high time-efficiency.

Supplementary Materials

The video of the experiment with eccentric mass is available online at https://youtu.be/hf8jx9WN6Ec. The video of the experiment with a motor fault is available online at https://youtu.be/HxRsjkjuvY4. The code on STM32F103 used in the experiments is available online at https://1drv.ms/f/s!ApzqIuEnpaPHhSaF6lRLL5F_qh0Y.

Author Contributions

Funding acquisition, N.C.; Investigation, F.Y.; Methodology, C.S. and F.Y.; Project administration, C.W.; Resources, N.C.; Software, C.S.; Supervision, C.W.; Writing – original draft, C.S.

Funding

This work is supported by the National Nature Science Fund of China (Grant No. 61403100), the Fundamental Research Funds for the Central Universities (Grant No. HIT.NSRIF.2015037) and the open National Defense Key Disciplines Laboratory of Exploration of Deep Space Landing and Return Control Technology, Harbin Institute of Technology (Grant No. HIT.KLOF.2013.079).

Acknowledgments

We gratefully acknowledge the assistance of Mo He in assembling the quadrotor used in our experiments. We also thank Jinbo Wang for processing videos and images.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Define:
a = ρ 1 2 + ρ 2 2 + ρ 3 2 1 2 = q 0 1 + q 0
The three elements of Ф ( x ) are:
L f 2 h 1 ( x ) = 1 4 ( ω z + ρ 1 ω y ρ 2 ω x ) ( ρ 3 ω x ρ 1 ω z + ρ 2 2 ω y a ω y + ρ 1 ρ 2 ω x + ρ 2 ρ 3 ω z ) 1 4 ( ω y ρ 1 ω z + ρ 3 ω x ) ( ρ 1 ω y ρ 2 ω x + ρ 3 2 ω z a ω z + ρ 1 ρ 3 ω x + ρ 2 ρ 3 ω y ) + 1 4 ( ρ 1 ω x + ρ 2 ω y + ρ 3 ω z ) ( ρ 2 ω z ρ 3 ω y + ρ 1 2 ω x a ω x + ρ 1 ρ 2 ω y + ρ 1 ρ 3 ω z ) + 1 4 I x ( I y I z ) ω y ω z ( ρ 1 2 ρ 2 2 ρ 3 2 + 1 ) + 1 2 I z ( I x I y ) ω x ω y ( ρ 1 ρ 3 + ρ 2 ) + 1 2 I y ( I z I x ) ω x ω z ( ρ 1 ρ 2 ρ 3 ) L f 2 h 2 ( x ) = 1 4 ( ω x + ρ 2 ω z ρ 3 ω y ) ( ρ 1 ω y ρ 2 ω x + ρ 3 2 ω z a ω z + ρ 1 ρ 3 ω x + ρ 2 ρ 3 ω y ) 1 4 ( ω z + ρ 1 ω y ρ 2 ω x ) ( ρ 2 ω z ρ 3 ω y + ρ 1 2 ω x a ω x + ρ 1 ρ 2 ω y + ρ 1 ρ 3 ω z ) + 1 4 ( ρ 1 ω x + ρ 2 ω y + ρ 3 ω z ) ( ρ 3 ω x ρ 1 ω z + ρ 2 2 ω y a ω y + ρ 1 ρ 2 ω x + ρ 2 ρ 3 ω z ) + 1 4 I y ( I z I x ) ω x ω z ( ρ 2 2 ρ 1 2 ρ 3 2 + 1 ) + 1 2 I z ( I x I y ) ω x ω y ( ρ 2 ρ 3 ρ 1 ) + 1 2 I x ( I y I z ) ω y ω z ( ρ 1 ρ 2 + ρ 3 ) L f 2 h 3 ( x ) = 1 4 ( ω y ρ 1 ω z + ρ 3 ω x ) ( ρ 2 ω z ρ 3 ω y + ρ 1 2 ω x a ω x + ρ 1 ρ 2 ω y + ρ 1 ρ 3 ω z ) 1 4 ( ω x + ρ 2 ω z ρ 3 ω y ) ( ρ 3 ω x ρ 1 ω z + ρ 2 2 ω y a ω y + ρ 1 ρ 2 ω x + ρ 2 ρ 3 ω z ) 1 4 ( ρ 1 ω x + ρ 2 ω y + ρ 3 ω z ) ( ρ 1 ω y ρ 2 ω x + ρ 3 2 ω z a ω z + ρ 1 ρ 3 ω x + ρ 2 ρ 3 ω y ) + 1 4 I z ( I x I y ) ω x ω y ( ρ 3 2 ρ 1 2 ρ 2 2 + 1 ) + 1 2 I y ( I z I x ) ω x ω z ( ρ 2 ρ 3 + ρ 1 ) + 1 2 I x ( I y I z ) ω y ω z ( ρ 1 ρ 3 ρ 2 )

References

  1. Bouabdallah, S.; Murrieri, P.; Siegwart, R. Design and control of an indoor micro quadrotor. In Proceedings of the IEEE International Conference on Robotics and Automation, New Orleans, LA, USA, 26 April–1 May; pp. 4393–4398.
  2. Özbek, N.S.; Önkol, M.; Efe, M.Ö. Feedback control strategies for quadrotor-type aerial robots: A survey. Trans. Inst. Meas. Control 2016, 38, 529–554. [Google Scholar] [CrossRef]
  3. Lee, H.; Kim, H.J. Trajectory tracking control of multirotors from modelling to experiments: A survey. Int. J. Control Autom. Syst. 2017, 15, 281–292. [Google Scholar] [CrossRef]
  4. Tayebi, A.; McGilvray, S. Attitude stabilization of a vtol quadrotor aircraft. IEEE Trans. Control Syst. Technol. 2006, 14, 562–571. [Google Scholar] [CrossRef]
  5. Mahony, R.; Kumar, V.; Corke, P. Multirotor aerial vehicles modeling, estimation, and control of quadrotor. IEEE Robot. Autom. Mag. 2012, 19, 20–32. [Google Scholar] [CrossRef]
  6. Pounds, P.; Mahony, R.; Corke, P. Modelling and control of a large quadrotor robot. Control Eng. Practice 2010, 18, 691–699. [Google Scholar] [CrossRef] [Green Version]
  7. Kendoul, F. Survey of advances in guidance, navigation, and control of unmanned rotorcraft systems. J. Field Robot. 2012, 29, 315–378. [Google Scholar] [CrossRef]
  8. Lee, D.; Kim, H.J.; Sastry, S. Feedback linearization vs. Adaptive sliding mode control for a quadrotor helicopter. Int. J. Control Autom. Syst. 2009, 7, 419–428. [Google Scholar] [CrossRef]
  9. Kim, J.; Kang, M.S.; Park, S. Accurate modeling and robust hovering control for a quad-rotor vtol aircraft. In Proceedings of the 2nd International Symposium on UAVs, Reno, NV, USA, June 2009; pp. 9–26. [Google Scholar]
  10. Lyu, X.M.; Zhou, J.N.; Gu, H.W.; Li, Z.X.; Shen, S.J.; Zhang, F. Disturbance observer based hovering control of quadrotor tail-sitter vtol uavs using h-infinity synthesis. IEEE Robot. Autom. Lett. 2018, 3, 2910–2917. [Google Scholar] [CrossRef]
  11. Dierks, T.; Jagannathan, S. Output feedback control of a quadrotor uav using neural networks. IEEE Trans. Neural Netw. 2010, 21, 50–66. [Google Scholar] [CrossRef] [PubMed]
  12. Aboudonia, A.; El-Badawy, A.; Rashad, R. Disturbance observer-based feedback linearization control of an unmanned quadrotor helicopter. Proc. Inst. Mech. Eng. Part I: J. Syst. Control Eng. 2016, 230, 877–891. [Google Scholar] [CrossRef]
  13. Aboudonia, A.; Rashad, R.; El-Badawy, A. Composite hierarchical anti-disturbance control of a quadrotor uav in the presence of matched and mismatched disturbances. J. Intell. Robot. Syst. 2018, 90, 201–216. [Google Scholar] [CrossRef]
  14. Xiao, B.; Yin, S. A new disturbance attenuation control scheme for quadrotor unmanned aerial vehicles. IEEE Trans. Ind. Inform. 2017, 13, 2922–2932. [Google Scholar] [CrossRef]
  15. Ahmed, N.; Chen, M. Sliding mode control for quadrotor with disturbance observer. Adv. Mech. Eng. 2018, 10, 16. [Google Scholar] [CrossRef]
  16. Benallegue, A.; Mokhtari, A.; Fridman, L. High-order sliding-mode observer for a quadrotor uav. Int. J. Robust Nonlinear Control 2008, 18, 427–440. [Google Scholar] [CrossRef]
  17. Besnard, L.; Shtessel, Y.B.; Landrum, B. Quadrotor vehicle control via sliding mode controller driven by sliding mode disturbance observer. J. Frankl. Inst. 2012, 349, 658–684. [Google Scholar] [CrossRef]
  18. Fethalla, N.; Saad, M.; Michalska, H.; Ghommam, J. Robust observer-based dynamic sliding mode conroller for a quadrotor uav. IEEE Access 2018, 6, 45846–45859. [Google Scholar] [CrossRef]
  19. Shi, D.; Wu, Z.; Chou, W. Super-twisting extended state observer and sliding mode controller for quadrotor uav attitude system in presence of wind gust and actuator faults. Electronics 2018, 7, 128. [Google Scholar] [CrossRef]
  20. Shi, D.; Wu, Z.; Chou, W. Generalized extended state observer based high precision attitude control of quadrotor vehicles subject to wind disturbance. IEEE Access 2018, 6. [Google Scholar] [CrossRef]
  21. Han, J. From pid to active disturbance rejection control. IEEE Trans. Ind. Electron. 2009, 56, 900–906. [Google Scholar] [CrossRef]
  22. Wu, Z.-H.; Guo, B.-Z. On convergence of active disturbance rejection control for a class of uncertain stochastic nonlinear systems. Int. J. Control 2017, 1–14. [Google Scholar] [CrossRef]
  23. Guo, B.Z.; Zhao, Z.L. Active Disturbance Rejection Control for Nonlinear Systems: An Introduction, 1st ed.; John Wiley & Sons: Hoboken, NJ, USA, 2017; pp. 6–9. [Google Scholar]
  24. Yang, H.J.; Cheng, L.; Xia, Y.Q.; Yuan, Y. Active disturbance rejection attitude control for a dual closed-loop quadrotor under gust wind. IEEE Trans. Control Syst. Technol. 2018, 26, 1400–1405. [Google Scholar] [CrossRef]
  25. Ma, D.L.; Xia, Y.Q.; Li, T.Y.; Chang, K. Active disturbance rejection and predictive control strategy for a quadrotor helicopter. IET Contr. Theory Appl. 2016, 10, 2213–2222. [Google Scholar] [CrossRef]
  26. Guo, Y.; Jiang, B.; Zhang, Y. A novel robust attitude control for quadrotor aircraft subject to actuator faults and wind gusts. IEEE/CAA J. Autom. Sinica 2018, 5, 292–300. [Google Scholar] [CrossRef]
  27. Levant, A. Robust exact differentiation via sliding mode technique. Automatica 1998, 34, 379–384. [Google Scholar] [CrossRef]
  28. Levant, A. Higher-order sliding modes, differentiation and output-feedback control. Int. J. Control 2003, 76, 924–941. [Google Scholar] [CrossRef]
  29. Basin, M.; Yu, P.; Shtessel, Y. Finite-and fixed-time differentiators utilising hosm techniques. IET Control Theory Appl. 2016, 11, 1144–1152. [Google Scholar] [CrossRef]
  30. Menard, T.; Moulay, E.; Perruquetti, W. Fixed-time observer with simple gains for uncertain systems. Automatica 2017, 81, 438–446. [Google Scholar] [CrossRef] [Green Version]
  31. Angulo, M.T.; Moreno, J.A.; Fridman, L. Robust exact uniformly convergent arbitrary order differentiator. Automatica 2013, 49, 2489–2495. [Google Scholar] [CrossRef]
  32. Li, P.; Ma, J.; Zheng, Z. Disturbance-observer-based fixed-time second-order sliding mode control of an air-breathing hypersonic vehicle with actuator faults. Proc. Inst. Mech. Eng. Part G: J. Aerosp. Eng. 2018, 232, 344–361. [Google Scholar] [CrossRef]
  33. Ni, J.; Liu, L.; Chen, M.; Liu, C. Fixed-time disturbance observer design for brunovsky system. IEEE Trans. Circuit Syst. Part 2: Express Birefs 2018, 65, 341–345. [Google Scholar] [CrossRef]
  34. Yu, P.; Shtessel, Y.; Edwards, C. Continuous higher order sliding mode control with adaptation of air breathing hypersonic missile. Int. J. Adapt. Control Signal Process. 2016, 30, 1099–1117. [Google Scholar] [CrossRef] [Green Version]
  35. Basin, M.V.; Yu, P.; Shtessel, Y.B. Hypersonic missile adaptive sliding mode control using finite-and fixed-time observers. IEEE Trans. Ind. Electron. 2018, 65, 930–941. [Google Scholar] [CrossRef]
  36. Yu, X.; Li, P.; Zhang, Y. The design of fixed-time observer and finite-time fault-tolerant control for hypersonic gliding vehicles. IEEE Trans. Ind. Electron. 2018, 65, 4135–4144. [Google Scholar] [CrossRef]
  37. Schaub, H.; Junkins, J.L. Stereographic orientation parameters for attitude dynamics: A generalization of the rodrigues parameters. J. Astron. Sci. 1996, 44, 1–19. [Google Scholar]
  38. Leishman, G.J. Principles of Helicopter Aerodynamics, 2nd ed.; Cambridge university press: Cambridge, UK, 2006; pp. 55–167. [Google Scholar]
  39. Basin, M.; Panathula, C.B.; Shtessel, Y. Multivariable continuous fixed-time second-order sliding mode control: Design and convergence time estimation. IET Control Theory Appl. 2016, 11, 1104–1111. [Google Scholar] [CrossRef]
  40. Madgwick, S. An efficient orientation filter for inertial and inertial/magnetic sensor arrays. Report x-io Univ. Bristol 2010, 25, 113–118. [Google Scholar]
Figure 1. Definition of coordinate frames.
Figure 1. Definition of coordinate frames.
Electronics 07 00357 g001
Figure 2. Graphic representation of original modified Rodrigues parameter (MRP) and shadow MRP.
Figure 2. Graphic representation of original modified Rodrigues parameter (MRP) and shadow MRP.
Electronics 07 00357 g002
Figure 3. Blade element.
Figure 3. Blade element.
Electronics 07 00357 g003
Figure 4. Tensionmeter, torquemeter and tachometer.
Figure 4. Tensionmeter, torquemeter and tachometer.
Electronics 07 00357 g004
Figure 5. Measurements and fitting curves. (a) Trust curve; (b) torque curve.
Figure 5. Measurements and fitting curves. (a) Trust curve; (b) torque curve.
Electronics 07 00357 g005
Figure 6. Structure of active disturbance rejection controller for second-order system.
Figure 6. Structure of active disturbance rejection controller for second-order system.
Electronics 07 00357 g006
Figure 7. Block diagram of the proposed fixed-time active disturbance rejection control (FTADRC) attitude control scheme.
Figure 7. Block diagram of the proposed fixed-time active disturbance rejection control (FTADRC) attitude control scheme.
Electronics 07 00357 g007
Figure 8. Structure of simulation model.
Figure 8. Structure of simulation model.
Electronics 07 00357 g008
Figure 9. Chart of calculating the wind speed around each propeller.
Figure 9. Chart of calculating the wind speed around each propeller.
Electronics 07 00357 g009
Figure 10. Random wind speed generated by Dryden model.
Figure 10. Random wind speed generated by Dryden model.
Electronics 07 00357 g010
Figure 11. Curves of roll angle.
Figure 11. Curves of roll angle.
Electronics 07 00357 g011
Figure 12. Curves of pitch angle.
Figure 12. Curves of pitch angle.
Electronics 07 00357 g012
Figure 13. Curves of yaw angle.
Figure 13. Curves of yaw angle.
Electronics 07 00357 g013
Figure 14. Curves of observation error for Δ τ x .
Figure 14. Curves of observation error for Δ τ x .
Electronics 07 00357 g014
Figure 15. Curves of observation error for Δ τ y .
Figure 15. Curves of observation error for Δ τ y .
Electronics 07 00357 g015
Figure 16. Curves of observation error for Δ τ z .
Figure 16. Curves of observation error for Δ τ z .
Electronics 07 00357 g016
Figure 17. Experimental system.
Figure 17. Experimental system.
Electronics 07 00357 g017
Figure 18. Deformation of quadrotor frame used in the experiments.
Figure 18. Deformation of quadrotor frame used in the experiments.
Electronics 07 00357 g018
Figure 19. A mass block hung from the right hand front corner is used as the eccentric mass.
Figure 19. A mass block hung from the right hand front corner is used as the eccentric mass.
Electronics 07 00357 g019
Figure 20. Curves of roll angle.
Figure 20. Curves of roll angle.
Electronics 07 00357 g020
Figure 21. Curves of pitch angle.
Figure 21. Curves of pitch angle.
Electronics 07 00357 g021
Figure 22. Curves of yaw angle.
Figure 22. Curves of yaw angle.
Electronics 07 00357 g022
Figure 23. Curves of disturbance estimation.
Figure 23. Curves of disturbance estimation.
Electronics 07 00357 g023
Figure 24. Curves of roll angle.
Figure 24. Curves of roll angle.
Electronics 07 00357 g024
Figure 25. Curves of pitch angle.
Figure 25. Curves of pitch angle.
Electronics 07 00357 g025
Figure 26. Curves of yaw angle.
Figure 26. Curves of yaw angle.
Electronics 07 00357 g026
Figure 27. Curves of disturbance estimation.
Figure 27. Curves of disturbance estimation.
Electronics 07 00357 g027
Figure 28. Experimental environment with wind and damaged propeller.
Figure 28. Experimental environment with wind and damaged propeller.
Electronics 07 00357 g028
Figure 29. Tracking errors of roll angle using different control schemes: (a) FTADRC; (b) conventional ADRC.
Figure 29. Tracking errors of roll angle using different control schemes: (a) FTADRC; (b) conventional ADRC.
Electronics 07 00357 g029
Figure 30. Tracking errors of pitch angle using different control schemes: (a) FTADRC; (b) conventional ADRC.
Figure 30. Tracking errors of pitch angle using different control schemes: (a) FTADRC; (b) conventional ADRC.
Electronics 07 00357 g030
Figure 31. Tracking errors of yaw angle using different control schemes: (a) FTADRC; (b) conventional ADRC.
Figure 31. Tracking errors of yaw angle using different control schemes: (a) FTADRC; (b) conventional ADRC.
Electronics 07 00357 g031
Figure 32. Curves of disturbance estimation using different control schemes: (a) FTADRC; (b) conventional ADRC.
Figure 32. Curves of disturbance estimation using different control schemes: (a) FTADRC; (b) conventional ADRC.
Electronics 07 00357 g032
Table 1. Parameters used in the aerodynamic model of propeller.
Table 1. Parameters used in the aerodynamic model of propeller.
ParametersDescriptionDetermination method
Atmospheric parameter ρ Air density1.2kg/m3 in low altitude
Propeller
parameters
R Propeller radiusMeasuring directly
c ( r n ) Blade chordMeasuring directly
θ ( r n ) Blade pitch angleMeasuring directly
N b Number of rotor bladesMeasuring directly
Velocity
parameters
λ c Climb inflow ratioCalculated according to velocity of quadrotor and wind speed
μ Advance ratioCalculated according to velocity of quadrotor and wind speed
Aerodynamic coefficients C l α Lift-curve slopeEstimated with trust-rotation speed curve
C d Drag coefficientEstimated with torque-rotation speed curve
Table 2. Chord length and pitch angle of the blade.
Table 2. Chord length and pitch angle of the blade.
Radial Position (mm)Chord (mm)Pitch Angle (°)
322935.2
484028
644821.4
805217.6
965514.2
1125612.7
1285511.3
1445110.6
160476.6
176386.0
192286.0
203.206.0
Table 3. Parameters of the quadrotor in simulations.
Table 3. Parameters of the quadrotor in simulations.
ParameterDescriptionTrue ValueNominal Value
IxInertia along xb-axis0.1 kgm20.05 kgm2
IyInertia along yb-axis0.1 kgm20.05 kgm2
IzInertia along zb-axis0.22 kgm20.5 kgm2
lDistance between rotor and centroid0.4 m0.4 m
mMass8 kg--
KvMotor velocity constant325 rpm/V--
R M Equivalent resistance of motor0.26 Ω--
Table 4. Control parameters used in the experiments.
Table 4. Control parameters used in the experiments.
ParametersNominal Value
I x 0.01 kgm2
I y 0.01 kgm2
I z 0.02 kgm2
C diag(5,5,30)
k 1 70
k 2 10
p 1 0.5
p 2 0.5

Share and Cite

MDPI and ACS Style

Song, C.; Wei, C.; Yang, F.; Cui, N. High-Order Sliding Mode-Based Fixed-Time Active Disturbance Rejection Control for Quadrotor Attitude System. Electronics 2018, 7, 357. https://doi.org/10.3390/electronics7120357

AMA Style

Song C, Wei C, Yang F, Cui N. High-Order Sliding Mode-Based Fixed-Time Active Disturbance Rejection Control for Quadrotor Attitude System. Electronics. 2018; 7(12):357. https://doi.org/10.3390/electronics7120357

Chicago/Turabian Style

Song, Chunlin, Changzhu Wei, Feng Yang, and Naigang Cui. 2018. "High-Order Sliding Mode-Based Fixed-Time Active Disturbance Rejection Control for Quadrotor Attitude System" Electronics 7, no. 12: 357. https://doi.org/10.3390/electronics7120357

APA Style

Song, C., Wei, C., Yang, F., & Cui, N. (2018). High-Order Sliding Mode-Based Fixed-Time Active Disturbance Rejection Control for Quadrotor Attitude System. Electronics, 7(12), 357. https://doi.org/10.3390/electronics7120357

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop