Next Article in Journal
Development of POLON—A Green Microsatellite Propulsion Module Utilizing 98% Hydrogen Peroxide
Previous Article in Journal
Structural Damage Assessment Using Multiple-Stage Dynamic Flexibility Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Control of a Single Tail Tilt Servomotor Tri-Rotor Ducted VTOL-UAV

1
School of Automation and Electrical Engineering, University of Science and Technology Beijing, Beijing 100083, China
2
Institute of Mechanics, Chinese Academy of Sciences, Beijing 100083, China
3
Institute of Engineering Thermophysics, Chinese Academy of Sciences, Beijing 100083, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Aerospace 2022, 9(6), 296; https://doi.org/10.3390/aerospace9060296
Submission received: 22 March 2022 / Revised: 21 April 2022 / Accepted: 25 April 2022 / Published: 31 May 2022
(This article belongs to the Section Aeronautics)

Abstract

:
This paper explores a novel single tail tilt servomotor tri-rotor ducted vertical takeoff and landing unmanned aerial vehicle system (VTOL-UAV), and proposes a compound control method combining disturbance observer, model predictive control (MPC) and sliding mode nonlinear dynamic inversion (NDI), and realizes the robust tracking control of the VTOL-UAV trajectory under compound disturbance. Specifically, the inner loop adopts nonlinear dynamic inverse which improved by sliding mode to realize the pseudo linearization of the system. The outer loop adopts the model predictive control based on the E-SSPC (State Space Predictive Controller based on the Error model) method, on this basis, the sliding mode disturbance observer based on fast Super-twisting algorithm is introduced into the position loop to observe and compensate the disturbance in real time, which improves the robustness of the outer loop system. Numerical simulation experiments verify the effectiveness and robustness of the control method. Finally, the flight test of the VTOL-UAV is carried out.

1. Introduction

Compound helicopter combines the advantages of vertical take-off, landing and hovering of helicopter, as well as the rapidity of fixed wing aircraft. It has broad application prospect and is one of the most important development directions of the helicopter. The Lockheed Martin Sikorsky-Boeing SB > 1 DEFIANT X helicopter exceeded 426 km/h (230 knots) in flight testing [1]. The Eurocopter X3 hybrid fixed wing helicopter opens the frontiers of aviation by attaining a speed milestone of 472 km/h (255 knots) in level flight, and reaches a speed of 487 km/h (263 knots) during a descent [2].
The compound helicopter has the characteristics of complex aerodynamic characteristics and multiple flight modes, so it is difficult to achieve the precision flight. There are many nonlinear control methods to solve the problem, such as gain scheduling, nonlinear dynamic inversion, backstepping, active disturbance rejection control, etc., each has its advantages and disadvantages. Different nonlinear control methods are combined to get better control effect than a single method. Reference [3] deduced the model of a compound helicopter, and predicted the maximum L/D (lift drag ratio) relation between the lift and propulsive force shares in different flight speeds and take-off weights, as well as rotor speed and blade twist. Reference [4] adopted the adaptive control to estimate and compensate the external disturbances to reduce the negative effects derived from actuator faults and input saturation, and a robust inner-outer loop adaptive FTC (Fault Tolerant Control) scheme improved by the backstepping was proposed to guarantee the tracking accuracy. Reference [5] divided the helicopter’s dynamic mode into the longitudinal-lateral and the heading-heave subsystem, and employed sliding mode control technique to compensate for the immeasurable flapping angles’ dynamic effects and external disturbances. In order to process the uncertainty and external disturbance of the longitudinal model of hypersonic vehicles, Reference [6] designed feedback linearization to realize the decoupling of the height and velocity channels, and applied nonlinear dynamic inverse control with the pole configuration method. Finally, sliding mode disturbance observer was used to estimate and compensate the equivalent disturbance. Reference [7] introduced extended state estimation algorithm with saturation constrainted observations to estimate the system disturbance and the system state. Via compensating the external disturbance and the coupling relationship among the system states, the active disturbance rejection controller was designed. Reference [8] proposed an improved piecewise constant adaptive NDI (Nonlinear Dynamic Inversion) control method suitable for general flight control systems, and gave the angular rates control strategy.
A novel tri-rotor ducted vertical takeoff and landing composite unmanned aerial vehicle system (Tri-Rotor Ducted VTOL-UAV) is developed by the authors, and it is equipped with a single tail tilt servomotor to delete the inverse torque and enhance controllability of yaw moment. The tri-rotor ducted VTOL-UAV is shown in the Figure 1. The two fronts of rotors rotate in opposite directions, which created reaction torque almost zero. The duct is applied to improve the efficiency of hovering and vertical flight, which is blended designed with wing body to reduce the resistance in forward flight. The working mode of the tail rotor tilt tri-rotor VTOL-UAV is as follows: it operates in helicopter mode when taking off and landing vertically, then the aircraft tilts forward and enters the low-speed forward flight state, when the flight speed increases to the transition threshold, the VTOL-UAV enters the conversion mode. When flying in airplane mode, the rotor system is locked out.
For the tri-rotor ducted VTOL-UAV, the aerodynamic characteristics are more complex than the compound UAV due to the effect of tilting tail, aerodynamic effects of ducts and rotors, three rotors instead of an even number of rotor, thus it is difficult to build an accurate model of the aircraft and obtain satisfactory flight accuracy through the conventional control method. Based on the relative independence of the inner and outer loops and the different requirements of response speed, the inner loop attitude controller and the outer loop position tracking controller are designed respectively. The dynamic inverse which can realize the pseudo linearization of nonlinear system by employing nonlinear cancellation has strong nonlinear processing ability [9]. Model predictive control (MPC) is suitable for slow varying systems, and has strong robustness and good control effect for position and speed control that has the characteristics of time delay, nonlinearity and uncertainty [10,11]. Therefore, a dual loop hierarchical compound control scheme which combines improved model predictive control and sliding mode nonlinear dynamic inversion (SMNDI) is proposed. The improved model predictive controller is designed for position control to provide the desired attitude command. For the attitude control, full state feedback on the angular motion model is carried out by using nonlinear dynamic inverse to counteract the nonlinear term, and a pseudo linear system with linear input and output relationship is obtained. However, there are some problems when combining model prediction with dynamic inversion, that is the dynamic inversion requires high-precision model, which also leads to insufficient robustness. The robustness of model predictive controller will be greatly challenged in the face of large disturbances. Because of the above problems and analysis, this paper introduces a sliding mode disturbance observer to observe and compensate the disturbance in real time, so as to strengthen the anti-interference of the position tracking controller [12,13].
The rest of this paper is organized as follows: In Section 1, the flight dynamic model of the tri-rotor ducted VTOL-UAV is established; In Section 2, the nonlinear control algorithm is studied; Section 3 carries out simulation test and flight test; Section 4 sets out the relevant conclusions.

2. Flight Dynamic Model

The motion equation of the tri-rotor ducted VTOL-UAV can be written as follows
P ˙ = V , V ˙ = g e 1 R F / m + d 1 , Θ ˙ = H ω , ω ˙ = J 1 ( ω × J ω ) + J 1 τ + d 2 ,
where, P = [ x , y , z ] T and V = [ u , v , w ] T are the position and velocity in inertial frame respectively; m is the mass; e 1 = [ 0 , 0 , 1 ] T ; g is the acceleration of gravity; R is the transformation matrix from body coordinate frame to inertial frame,
R = c o s θ c o s ψ c o s ψ s i n θ s i n ϕ c o s ϕ s i n ψ c o s ϕ c o s ψ s i n θ + s i n ϕ s i n ψ c o s θ s i n ψ s i n ψ s i n θ s i n ϕ + c o s ϕ c o s ψ c o s ϕ s i n ψ s i n θ s i n ϕ c o s ψ s i n θ c o s θ s i n ϕ c o s θ c o s ϕ ,
d 1 is acceleration disturbance, which is caused by external force. d 2 is angular acceleration disturbance, which is caused by external torque. Both d 1 and d 2 are the unknown external bounded disturbances. Θ = [ ϕ , θ , ψ ] T is the attitude angle including roll angle, pitch angle and yaw angle. and ω = [ p , q , r ] T is the attitude rate including roll rate, pitch rate and yaw rate. H is the attitude vector transformation matrix,
H = 1 s i n ϕ t a n θ c o s ϕ t a n θ 0 c o s ϕ s i n ϕ 0 s i n ϕ / c o s θ c o s ϕ / c o s θ ,
F = [ X , Y , Z ] T is the total aerodynamic force acting on the aircraft, and τ = [ L , M , N ] T is the resultant moment in body coordinate frame.
After obtaining the above equations, it is necessary to calculate the forces and moments that acting on the VTOL-UAV, which are shown in Figure 2. The total aerodynamic force acting on the aircraft including the F r m r and F l m r generated by the left and right rotors, the F t r generated by the tilting tail rotor, the aerodynamic force acting on the fuselage F f u s , aerodynamic force on horizontal stabilizer F h f and vertical fin F v f . In addition, there is gravity on the aircraft m g .
For the ducted rotor system, the influence of the duct on rotors is reflected as follows [14]: the duct makes the air flow through the rotor more concentrated and stable, and the vortex loss at the blade tip is less and thus the efficiency of providing lift is higher. The duct is similar to an annular wing while vertical motion, which makes the internal air velocity greater than the external air velocity and generates additional lift. The calculation of the rotors with culvert is based on the vortex theory of the rotor in the axial flow created by Rukovsky [14,15]. Due to the mature calculation method of the rotor pull and torque, the calculation results are directly given,
T = ( 1 + k ) A { 0.5 ρ [ ( ( 1 ζ χ 2 ) v 0 ) 2 + ( v 0 s i n α + v ) 2 ] + 0.5 ρ [ ( ( 1 ζ χ 1 ) v 0 ) 2 + ( v 0 s i n α + v i ) 2 ] 0.5 ρ ( 1 ζ v i ) v 0 2 }
Q = [ C T ( λ 0 μ z ) + C D 0 σ / 8 ( 1 + 7 / 3 μ 2 ) ] ρ ( Ω R ) 2 π R 3
where, T is the pull force of the rotor system. Q is the torque of the rotor system. A is rotor disk area. ρ is Earth’s atmosphere density. α is the attack angle of blade element. k = 0.2 is the correction coefficient of the pull increment caused by the culvert [16,17]. C T = T / ρ ( Ω R ) 2 π R 2 represents the pull coefficient. C D 0 is the blade drag coefficient. v 0 , v 1 and v 2 are the upstream infinite air flow velocity, the induced velocity and the downstream infinite air flow velocity respectively. v is the increase of the vertical velocity at the downstream infinity compared with that at the upstream infinity. ζ 1 and ζ 2 are the horizontal velocity blocking coefficient of the air flow in the rotor plane due to the culvert and the retardation coefficient at the downstream infinity. ζ v i is the local loss coefficient caused by the bending of the fluid at the culvert inlet. μ z represents the dimensionless form of the vertical velocity component. μ represents the induced velocity at the rotor disc. λ 0 is the rotor inflow ratio. σ is rotor solidity. R is rotor disc radius. Ω is the blade rotation rate.
The influence of the wake generated by rotor blade rotation on the thrust propeller mainly depends on whether it is in the rotor downwash area [18]. Considering the influence of the culvert and the relative position of the rotor and propeller, there is the following conclusion: the aerodynamic interference between the culvert rotor and the propulsion propeller is small and can be ignored. Referring to the calculation method of rotor tension, the thrust and torque of propeller are respectively,
T T = π ρ R p 2 ( V p + v p ) v p
C T p = 0.5 N p ρ Ω p 2 c p 0 R p C L p ϕ p r 3 d r
where, V p and v p are the flight velocity and the induced speed at the propeller disc; V p + v p is the air velocity flowing through the propulsion propeller, which is superimposed by the flow velocity and induced velocity [19], and T p is the rotor pull.
The aerodynamic force of the fuselage is calculated by CFD method. The lift coefficient and drag coefficient with the attack angle are shown in Figure 3 and Figure 4 respectively. It can be seen that the lift coefficient has linear characteristics in the range of attack angle 0 ° to 10 ° , and the lift line slope decreases when the attack angle is greater than 10 ° . In order to ensure the safe flight of the aircraft, the attack angle should be controlled within 0 ° to 10 ° .
From the Figure 4, it can be seen that the drag coefficient has linear characteristics in the range of attack angle 0 ° to 12 ° . Besides, the drag coefficient is a little large which leads to small lift drag ratio, because the aircraft drag is mainly generated by the culvert and the head of the aircraft, as shown in Figure 5, the points with large pressure of the whole aircraft are concentrated in the culvert and nose [20]. The aircraft resistance is small when it flies at a small rake angle, and resistance increases rapidly with the increase of flight speed. The aerodynamic shape of the aircraft is optimized to increase lift and reduce drag. The research result will be published in another paper.

3. Controller Design

In this section, preliminaries of nonlinear dynamic inversion and model predictive control are introduced, and they will be further used in the analysis and design of control algorithm.
The tri-rotor ducted VTOL-UAV has obvious nonlinear characteristics caused by rotors and ducts, and the longitudinal-heading coupling is serious due to the left-right tilt of tail rotor. In addition, the influence of external disturbance needs to be considered. These factors bring great difficulties to the controller design. The combination of model predictive control and sliding mode dynamic inversion are used in this paper. The inner loop adopts the nonlinear dynamic inversion improved with the sliding mode idea to cancel the nonlinearity of the attitude and realize the pseudo linearization of the system, and the outer loop adopts the model predictive control based on the E-SSPC method. On this basis, the sliding mode disturbance observer based on fast super-twisting algorithm is introduced to observe the compound disturbance in the position loop, and the observed value compensates the original model predictive control law in real time, so as to improve the robustness and strengthen the anti-interference of the outer loop position control system. The complete UAV trajectory tracking hierarchical compound controller is shown in Figure 6 [20].

3.1. Attitude Dynamic Inverse Controller Design

In order to carry out the dynamic inverse control design, the attitude control system dynamic model in Equation (1) can be transformed into the following affine form [20]
Θ ˙ = f 1 + g 1 ω
ω ˙ = f 2 + g 2 u + d 1
where, u = [ L , M , N ] T , f 1 and f 2 are state matrices of attitude and attitude angle rate, f 1 = 0 3 × 1 , d 1 is the external unknown bounded disturbance.
g 1 = 1 s i n ϕ t a n θ c o s ϕ t a n θ 0 c o s ϕ s i n ϕ 0 s i n ϕ / c o s θ c o s ϕ / c o s θ
f 2 = [ J z ( J y J z ) J x z 2 ] q r J x J z J x z 2 + J x z ( J x J y + J x z ) p q J x J z J x z 2 1 J y [ ( J z J x ) p r J x z ( p 2 r 2 ) ] J x z ( J y J z ) J x ) q r J x J z J x z 2 + [ J x ( J x J y ) + J x z 2 ] p q J x J z J x z 2
g 2 = J z J x J z J x z 2 0 J x z J x J z J x z 2 0 1 J y 0 J x z J x J z J x z 2 0 J x J x J z J x z 2
where, J x , J y , J z are the moment of inertia about roll-axis, yaw-axis and pitch-axis respectively, and J x z is cross moment of inertia about roll axis and pitch axis. Other physical meanings are the same as above.
According to the difference of response speed, the attitude subsystem can be divided into angle rate loop and attitude angle loop [21], and the control laws can be designed respectively.

3.1.1. Attitude Angle Rate Loop Controller Design

For the angle rate of UAV, the following form of sliding mode dynamic inverse control law can be designed
s 2 = ω ω d μ 2 = k ω s 2 ϵ 2 s i g n ( s 2 ) + ω d ˙ u = g 2 1 ( μ 2 f 2 )
where, s 2 is an auxiliary variable; k ω and ϵ 2 are positive definite diagonal matrix; ω d is the desired attitude angle rate, and it is given by the attitude angle sliding mode dynamic inverse controller; u is control input of attitude angle rate control.
Theorem 1.
For attitude angle system (7), under the action of control law (10), the actual attitude angle rate ω will converge to the desired attitude angle rate ω d asymptotically.
Proof. 
Substituting the control law (10) to the system (7) to obtain
s 2 ˙ = k ω s 2 ϵ 2 s i g n ( s 2 ) + d 1
Selecting the Lyapunov candidate L = 1 2 s 2 T s 2 , both sides of L ˙ = s 2 T s ˙ 2 is derivated, and substitute Equation (11) into it, there is
L ˙ = s 2 T k ω s 2 ϵ 2 s i g n s 2 + d 1 s 2 T k ω s 2 ϵ 2 s i g n s 2 + Φ
When selecting appropriate control parameter ϵ 2 , it is easy to obtain L ˙ s 2 T k ω s 2 . According to the BarbaLat lemma, because k ω and ϵ 2 are positive definite matrix, when t, there is s 2 → 0, according to the definition of auxiliary variables s 2 , there is an ω convergence to ω d asymptotically.
It can be seen from (12) that k ω s 2 and ϵ 2 s i g n ( s 2 ) have different functions, the k ω s 2 is used to realize the asymptotic tracking of s 2 , while the ϵ 2 s i g n ( s 2 ) is used for dealing with the uncertainty d 1 , their joint action ensures the balance of performance and robustness.

3.1.2. Attitude Angle Loop Controller Design

The attitude angle of the VTOL-UAV is an independent slowly varying state variable, the following form of sliding mode dynamic inverse controller can be designed
s 1 = Θ Θ d μ 1 = k Θ s 1 ϵ 1 s i g n ( s 1 ) + Θ d ˙ u = g 1 1 ( μ 1 f 1 )
where, s 1 is an auxiliary variable; k Θ and ϵ 1 are positive definite diagonal matrix; Θ d is the desired attitude angle, and it is given by the outer loop MPC position controller; u is control input of attitude control.

3.2. Model Predictive Position Controller Design

3.2.1. Prediction Model

The prediction model is applied to the position control system. Considering the external bounded disturbance, the position subsystem can be written as follows
P ˙ = V
V ˙ = F / m + d 2
in which, d 2 is the external unknown bounded disturbance, meanings of the other variables are the same as above.
Now, suppose ξ ( t ) = [ x ( t ) , y ( t ) , z ( t ) , u ( t ) , v ( t ) , w ( t ) ] T and μ = [ μ 1 ( t ) , μ 2 ( t ) , μ 3 ( t ) ] T as the new state varieties and control inputs, the (14) and (15) can be written as
ξ ˙ ( t ) = A ( t ) ξ ( t ) + B ( t ) μ ( t ) + Δ
The goal of position controller is to make the aircraft flight along the desired trajectory with as little error as possible under the action of control μ . The desired trajectory is designed as ξ r ( t ) = [ x r ( t ) , y r ( t ) , z r ( t ) , u r ( t ) , v r ( t ) , w r ( t ) ] T , and the error dynamics of the position subsystem can be obtained as
ξ ˜ ˙ ( t ) = A ( t ) ξ ˜ ( t ) + B ( t ) μ ˜ ( t ) + Δ
in which, ξ ˜ ( t ) = ξ ( t ) ξ r ( t ) , μ ( t ) = μ ( t ) μ r ( t ) ,
A ( t ) = 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ,   B ( t ) = 0 0 0 0 0 0 0 0 0 1 / m 0 0 0 1 / m 0 0 0 1 / m ,   Δ ( t ) = 0 0 0 d 1 ( t ) d 2 ( t ) d 3 ( t ) .
At first, the uncertainty Δ is not considered temporarily, and the backward difference is used to discretize the above Equation (17), supposing the discrete time step is Δ t , there is
ξ ˜ ( k + 1 ) = A ( k ) ξ ˜ ( k ) + B ( k ) μ ˜ ( k )
where,
A ( k ) = 1 Δ t 0 1 0 0 0 1 0 0 1 0 0 0 1 Δ t 0 1 0 0 0 1 0 0 0 0 0 0 1 Δ t 0 0 0 0 0 1 ,   B ( k ) = 0 0 0 Δ t / m 0 0 0 0 0 0 Δ t / m 0 0 0 0 0 0 Δ t / m .
Let N p and N u are the predicted time length and the control time length respectively, and there is N p N u . Using Equation (17), through iteration and sorting, the expression of the system prediction model can be deduced as follows
ξ ˜ ^ = P ( t ) ξ ˜ ( t ) + H ( t ) u ˜
where,
ξ ˜ ^ ξ ˜ ( t + Δ t ) ξ ˜ ( t + N p Δ t ) ,   u ˜ u ˜ ( t ) ξ ˜ ( t + ( N u 1 ) Δ t ) ,
ξ r ˜ ^ ξ r ˜ ( t + Δ t ) ξ r ( t ) ξ r ˜ ( t + N p Δ t ) ξ r ( t ) ,   u r ˜ ^ μ r ( t ) μ r ( t Δ t ) μ r ˜ ( t + ( N u 1 ) Δ t ) μ r ( t Δ t ) ,
P ( t ) A ( t ) A ( t ) A ( t + Δ t ) α ( t , 0 , 2 ) α ( t , 0 , 1 ) ,
in which,
α ( t , j , l ) i = j N p l A ( t + i Δ t ) ,
H ( t ) B ( t ) 0 0 A ( t + Δ t ) B ( t ) B ( t + Δ t ) 0 α ( t , 1 , 2 ) B ( t ) α ( t , 2 , 2 ) B ( t + Δ t ) 0 α ( t , 1 , 1 ) B ( t ) α ( t , 2 , 1 ) B ( t + Δ t ) B ( t + ( N u 1 ) Δ t ) .

3.2.2. Optimal Control Law

Model predictive control is essentially an optimal control based on finite time domain. To obtain the optimal control law, the penalty function about the terminal state ξ ˜ ( t + Δ t ) based on the traditional linear quadratic cost function is introduced
Ω = [ ξ ˜ ^ ( t + N p Δ t ) ξ ˜ r ( t + N p Δ t ) ] T G [ ξ ˜ ^ ( t + N p Δ t ) ξ ˜ r ( t + N p Δ t ) ]
and the terminal state ξ ˜ ( t + N p Δ t ) can be written as
ξ ˜ ( t + Δ t ) = α ( t , 0 , 1 ) ξ ˜ ( t ) + M u ˜
where, M = [ α ( t , 1 , 1 ) B ( t ) α ( t , 2 , 1 ) B ( t + Δ t ) B ( t + ( N u 1 ) Δ t ) ] , and it is the last row of the coefficient matrix H(t).
Therefore, the final objective cost function is defined as
J = [ ξ ˜ ^ ξ ˜ r ] T Q [ ξ ˜ ^ ξ ˜ r ] + [ u ˜ u ˜ r ] T R [ u ˜ u ˜ r ] + Ω
where, ξ ˜ ^ is the predictive output of the vector ξ ˜ ; Q and R are the state weight matrix and the control weight matrix respectively. The optimal control law can be obtained by calculating the partial derivative of Equation (22) with respect to the control u ˜ and making
J u ˜ = 0
The optimal control law is obtained by solving the Equation (23),
u ˜ = ( H T Q H + R + M T G M ) 1 H T Q [ ξ ˜ r P ξ ˜ ] + R u ˜ r + M T G [ ξ ˜ r ( t + N p Δ t ) α ( t , 0 , 1 ) ξ ˜ ( t ) ]
When the control sequence is obtained, the observed value Δ ^ of the uncertainty Δ in (17) estimated by the sliding mode disturbance observer is introduced for compensation, and the final control law of the position subsystem is obtained
μ ˜ = u ˜ ( t ) B 1 Δ ^
From the definition of control error μ ˜ = μ μ r , the desired control input can be obtained as
μ d = μ ˜ + μ r
From the relationship between the control inputs and the forces, the desired ϕ d , θ d and ψ d can be obtained, which together serve as the desired command of attitude control and provide a reference command for the sliding mode dynamic inverse controller.

3.2.3. Sliding Mode Disturbance Observer Design

In order to effectively deal with the problem of insufficient robustness caused by external disturbance, and improve the control accuracy, the fast super-twisting algorithm is used to design the sliding mode disturbance observer to observe and estimate the external uncertainty Δ continuously in real time [22,23,24,25,26]. The design of the sliding mode disturbance observer is as follows
s 1 = ξ ˜ z 1 z ˙ 1 = A ξ ˜ + B μ ˜ + Δ ^ Δ ^ ˙ = ω ¯ 1 d i a g ( | s 1 | 0.5 ) s i g n s 1 + ω ¯ 2 s 1 + 0 t [ ω ¯ 3 s i g n s 1 ( τ ) + ω ¯ 4 s 1 ( τ ) ] d τ
where, s 1 = [ s 11 , s 12 , s 13 , s 14 , s 15 , s 16 ] T is the auxiliary sliding mode vector, Δ ^ R 6 represents the sliding mode control input, and it is also the observation output of the disturbance, ω ¯ i = d i a g ω ¯ i 1 , ω ¯ i 2 , ω ¯ i 3 , ω ¯ i 4 , ω ¯ i 5 , ω ¯ i 6 , i = 1 , 2 , 3 , 4 , in which, each component is greater than zero, | s 1 | 0.5 s i g n ( s 1 ) = [ | s 11 | 0.5 s i g n ( s 11 ) , , | s 16 | 0.5 s i g n ( s 16 ) ] T .

4. Numerical Simulation and Flight Test

4.1. Simulation

To verify the effectiveness of the control method, the numerical simulation is designed to test the trajectory tracking performance of the controller and its robustness to compound disturbances. Set the total simulation time as 30 s, and the reference trajectory expression for design tracking is
x r ( t ) = 7 s i n ( π t / 7 ) y r ( t ) = 14 s i n ( π t / 14 + 1 / 36 ) z r ( t ) = 130
The initial position of the UAV is P 0 = [ 0 , 0 , 57 ] m, the initial velocity is V 0 = [ 0.0 , 0.0 , 3.3 ] m/s, other variables are set to zero. At the same time, compound disturbance is added into the velocity loop and angle rate loop respectively from the beginning of flight, and the disturbance expression is as follows
d 1 = [ 1.2 s i n ( 3 t ) + 0.5 c o s ( t ) , s i n ( 3 t ) + 0.5 c o s ( t ) , 0.2 s i n ( t ) + 0.5 c o s ( 3 t ) ] T
d 2 = [ 0.3 c o s ( 5 t ) + 0.15 , 0.3 s i n ( 5 t ) + 0.2 , 0.3 s i n ( 5 t ) + 0.12 ] T
where, the units of d 1 and d 2 are m/s 2 and rad/s 2 respectively. The parameters of model prediction controller are as follows, prediction step N 2 = 10 , control step N u = 5 ; The state weight matrix and control weight matrix are Q = d i a g 160 , 8 , 160 , 8 , 140 , 4 , R = d i a g 0.004 , 0.002 , 0.002 , the parameters of disturbance observer are ω ¯ 1 = d i a g { 0.84 , 0.84 , 0.84 , 0.336 , 0.336 , 0.336 } , ω ¯ 2 = d i a g { 9.4 , 9.4 , 9.4 , 3.2 , 3.2 , 3.2 } , ω ¯ 3 = d i a g { 0.75 , 0.75 , 0.75 , 0.325 , 0.325 , 0.325 } , ω ¯ 4 = d i a g { 0.48 , 0.48 , 0.48 , 0.208 , 0.208 , 0.208 } . The parameters of sliding mode dynamic inverse controller are k Θ = 8.4 , ξ = 0.01 , k ω = 8.4 , ξ = 0.02 . The simulation result is shown in Figure 7, Figure 8, Figure 9 and Figure 10.
The Figure 7 shows the tracking of the desired trajectory in three-dimensional space. It can be seen from the figure that the compound method proposed can realize the accurate tracking of the desired trajectory.
The Figure 8 shows the attitude angle and their tracking of the UAV when horizontal flight. The angle command is reasonable and tracking effect is ideal, the tracking error is small, which indicates that the sliding mode dynamic inverse method has strong robustness to compound disturbance.
Figure 9 and Figure 10 are the estimate of the composite disturbance d 1 and d 2 obtained by the sliding mode disturbance observer. It can be seen that the the desired and actual curves coincide basically, and the designed sliding mode disturbance observer has high tracking accuracy for composite disturbance. The observation and estimation of external disturbance has laid a foundation for disturbance compensation.
Figure 11, Figure 12, Figure 13 and Figure 14 are the control inputs of the VTOL-UAV, the amplitude of each control input is kept within a reasonable range, and the control is adjusted repeatedly due to the existence of composite disturbance. From the above analysis, the compound control method proposed in this paper has good anti-interference, and realizes the trajectory tracking control of the VTOL-UAV under the combined action of force and torque disturbance. At the same time, the controller has high tracking accuracy, rapid response, small overshoot, and the control amplitude value is kept within a reasonable range.

4.2. Flight Test

The flight test of 8-word route tracking is carried out. The conditions are set as follows: the initial height is 57 m, the desired height is 130 m, and the desired horizontal velocity is 20 m/s. The test results are shown in Figure 15, Figure 16, Figure 17 and Figure 18.
From the flight test results, the VTOL-UAV can complete the 8-word route flight with high consistency. Specifically, the roll control has high precision and fast response. There is obvious time delay in pitch angle and heading angle control. The altitude control accuracy is within ±4 m, which is due to the large aircraft mass, slow system response and large control force. The longitudinal control and heading control are mainly realized through the tail rotor, with small rotor disc area and small control force. At the same time, the pitch channel is coupled with the heading channel, resulting in height and heading oscillation. This phenomenon can be improved by optimizing the control strategy in the transition process and limiting the change rate of pitch angle command.

5. Conclusions

In this paper, a compound control method combining model prediction control improved with disturbance observer and sliding mode nonlinear dynamic inversion is proposed for the three culvert rotor composite VTOL-UAV with tiltable tail rotor. The robust trajectory tracking control with the composite disturbance is realized. First, the nonlinear dynamic inverse is improved with the sliding mode for the inner loop, and the controller is designed. For the outer loop, the E-SSPC method is used to deduce the prediction model, and the model predictive controller is designed by using nominal error dynamic system. The sliding mode disturbance observer based on super-twisting algorithm is introduced into the position controller to estimate and compensate the composite disturbance, which enhances the robustness of the model predictive controller. The effectiveness and robustness of the compound control scheme are verified by numerical simulation experiments. Finally, the flight test is carried out.

Author Contributions

Conceptualization, Y.H., J.G. and P.Y.; methodology, Y.H. and J.G.; software, J.G.; data curation, P.Y. and G.Z.; writing—original draft preparation, Y.H.; writing—review and editing, Y.H.; supervision, N.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported in part by the Beijing Natural Science Foundation (4222050) and in part by the National Natural Science Foundation of China (62173030).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors gratefully acknowledge the reviewers for their careful work and thoughtful suggestions that have helped improve this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Future Long Range Assault Aircraft. Available online: https://www.boeing.com/defense/future-long-range-assault-aircraft/index.page{#}/videos/defiant-x-exceeds-230-knots (accessed on 21 March 2022).
  2. Helicopters History. Available online: https://www.airbus.com/en/who-we-are/our-history/helicopters-history/x3 (accessed on 21 March 2022).
  3. Yang, K.; Han, D.; Shi, Q. Study on the lift and propulsive force shares to improve the flight performance of a compound helicopter. Chin. J. Aeronaut. 2021, 35, 365–375. [Google Scholar] [CrossRef]
  4. Yan, K.; Chen, M.; Wu, Q.; Zhu, R. Robust adaptive compensation control for unmanned autonomous helicopter with input saturation and actuator faults. Chin. J. Aeronaut. 2019, 32, 2299–2310. [Google Scholar] [CrossRef]
  5. Xian, B.; Guo, J.; Zhang, Y.; Zhao, B. Sliding mode tracking control for miniature unmanned helicopters. Chin. J. Aeronaut. 2015, 28, 277–284. [Google Scholar] [CrossRef] [Green Version]
  6. Wang, C.; Wang, X. Design of a controller for hypersonic vehicle based on sliding mode disturbance observer. Electron. Opt. Control. 2020, 27, 90–94. [Google Scholar]
  7. Tan, S.; Guo, J.; Zhao, Y.; Zhang, J. Adaptive control with saturation-constrainted observations for drag-free satellites—A set-valued identification approach. Sci. China Inf. Sci. 2021, 64, 182–193. [Google Scholar] [CrossRef]
  8. Li, Y.; Liu, X.; He, Q. An adaptive dynamic inversion control method based on improved piecewise constant for flight control system. J. Northwestern Polytech. Univ. 2021, 39, 167–174. [Google Scholar] [CrossRef]
  9. Lu, B. Research on High-Bandwidth/High-Precision Longitudinal Control for a Hypersonic Vehicle; Nanjing University of Aeronautics and Astronautics: Nanjing, China, 2008. [Google Scholar]
  10. Zhou, H. Small-Scale Unmanned Helicopter Modeling and Controller Design; South China University of Technology: Guangzhou, China, 2011. [Google Scholar]
  11. Wang, X.; Meng, X.; Li, C. Design of trajectory tracking controller for UAV based on MPC. Syst. Eng. Electron. 2021, 43, 191–198. [Google Scholar]
  12. Chen, J.; Sun, R.; Zhu, B. Disturbance observer-based control for small nonlinear UAV systems with transient performance constraint. Aerosp. Sci. Technol. 2020, 105, 106028. [Google Scholar] [CrossRef]
  13. Li, Y.; Chen, M.; Ge, S.S.; Li, D. Anti-disturbance control for attitude and altitude systems of the helicopter under random disturbances. Aerosp. Sci. Technol. 2019, 96, 105561. [Google Scholar] [CrossRef]
  14. Wang, Q. A Numerical Analysis of the Aerodynamic Performance of the Ducted-Fan UAV; National University of Defense Technology: Changsha, China, 2008. [Google Scholar]
  15. Wang, S. Theory and Simulation on Overall Design and Control System of the Saucer Ducted Unmanned Aerial Vehicle; Fudan University: Shanghai, China, 2011. [Google Scholar]
  16. Kong, W. Research on some Key Technical Issues of Composed High Speed Helicopter; Nanjing University of Aeronautics and Astronautics: Nanjing, China, 2011. [Google Scholar]
  17. Hu, Y. Research on System-Identification-Based Dynamic Modelling and Control Technology for High-Speed Helicopter; Beihang University: Beijing, China, 2020. [Google Scholar]
  18. Shen, H. Research in the Interaction Flow Field of Compound High-Speed Helicopter Rotor/Fuselage/Auxiliary Propulsor; Nanjing University of Aeronautics and Astronautics: Nanjing, China, 2015. [Google Scholar]
  19. Wang, Q.; Chen, M.; Xu, G. Control strategy for compound coaxial helicopter under conversion mode. J. Aerosp. Power 2014, 29, 458–466. [Google Scholar]
  20. Chen, N. Research on Nonlinear Robust Control Technology for Unmanned Helicopter; Beihang University: Beijing, China, 2018. [Google Scholar]
  21. Bai, Y.; Shi, X.; Wang, S. Auv attitude control system design with robust nonlinear dynamic inversion. Tactical Missile Technol. 2016, 2, 93–97. [Google Scholar]
  22. Dávila, A.; Moreno, J.A.; Fridman, L. Variable gains super-twisting algorithm: A lyapunov based design. In Proceedings of the 2010 American Control Conference, Baltimore, MD, USA, 30 June–2 July 2010; IEEE: Piscataway, NJ, USA, 2010; pp. 968–973. [Google Scholar]
  23. Evangelista, C.; Puleston, P.; Valenciaga, F. Lyapunov-designed super-twisting sliding mode control for wind energy conversion optimization. IEEE Trans. Ind. Electron. 2013, 60, 538–545. [Google Scholar] [CrossRef]
  24. Moreno, J.A.; Osorio, M. Strict lyapunov functions for the super-twisting algorithm. IEEE Trans. Autom. Control. 2012, 57, 1035–1040. [Google Scholar] [CrossRef]
  25. Tan, J.; Zhou, Z.; Zhu, X.; Xu, X. Fast super twisting algorithm and its application to attitude control of flying wing UAV. Control. Decis. 2016, 31, 143–148. [Google Scholar]
  26. Yang, F.; Tan, S.; Xue, W.; Guo, J.; Zhao, Y. Extended state filtering with saturation-constrainted observations and active disturbance rejection control of position and attitude for drag-free satellites. Acta Autom. Sin. 2020, 46, 2337–2349. [Google Scholar]
Figure 1. The tail rotor tilting three ducted fan VTOL-UAV.
Figure 1. The tail rotor tilting three ducted fan VTOL-UAV.
Aerospace 09 00296 g001
Figure 2. The forces on the tri-rotor ducted VTOL-UAV.
Figure 2. The forces on the tri-rotor ducted VTOL-UAV.
Aerospace 09 00296 g002
Figure 3. The lift coefficient with attack angle of the fuselage.
Figure 3. The lift coefficient with attack angle of the fuselage.
Aerospace 09 00296 g003
Figure 4. The drag coefficient with attack angle of the fuselage.
Figure 4. The drag coefficient with attack angle of the fuselage.
Aerospace 09 00296 g004
Figure 5. The fuselage drag distribution of the tri-rotor ducted VTOL-UAV.
Figure 5. The fuselage drag distribution of the tri-rotor ducted VTOL-UAV.
Aerospace 09 00296 g005
Figure 6. The diagram of layered compound controller.
Figure 6. The diagram of layered compound controller.
Aerospace 09 00296 g006
Figure 7. The desired trajectory tracking.
Figure 7. The desired trajectory tracking.
Aerospace 09 00296 g007
Figure 8. The attitude angle response.
Figure 8. The attitude angle response.
Aerospace 09 00296 g008
Figure 9. The estimate of the composite disturbance d1.
Figure 9. The estimate of the composite disturbance d1.
Aerospace 09 00296 g009
Figure 10. The estimate of the composite disturbance d2.
Figure 10. The estimate of the composite disturbance d2.
Aerospace 09 00296 g010
Figure 11. The left main rotor control input of the VTOL-UAV.
Figure 11. The left main rotor control input of the VTOL-UAV.
Aerospace 09 00296 g011
Figure 12. The right main rotor control input of the VTOL-UAV.
Figure 12. The right main rotor control input of the VTOL-UAV.
Aerospace 09 00296 g012
Figure 13. The tail rotor speed control input of the VTOL-UAV.
Figure 13. The tail rotor speed control input of the VTOL-UAV.
Aerospace 09 00296 g013
Figure 14. The tail rotor tilt angle control input of the VTOL-UAV.
Figure 14. The tail rotor tilt angle control input of the VTOL-UAV.
Aerospace 09 00296 g014
Figure 15. The flight test.
Figure 15. The flight test.
Aerospace 09 00296 g015
Figure 16. The flight trajectory.
Figure 16. The flight trajectory.
Aerospace 09 00296 g016
Figure 17. The attitude angle response of the tri-rotor ducted VTOL-UAV.
Figure 17. The attitude angle response of the tri-rotor ducted VTOL-UAV.
Aerospace 09 00296 g017
Figure 18. The height response of the tri-rotor ducted VTOL-UAV.
Figure 18. The height response of the tri-rotor ducted VTOL-UAV.
Aerospace 09 00296 g018
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hu, Y.; Guo, J.; Ying, P.; Zeng, G.; Chen, N. Nonlinear Control of a Single Tail Tilt Servomotor Tri-Rotor Ducted VTOL-UAV. Aerospace 2022, 9, 296. https://doi.org/10.3390/aerospace9060296

AMA Style

Hu Y, Guo J, Ying P, Zeng G, Chen N. Nonlinear Control of a Single Tail Tilt Servomotor Tri-Rotor Ducted VTOL-UAV. Aerospace. 2022; 9(6):296. https://doi.org/10.3390/aerospace9060296

Chicago/Turabian Style

Hu, Yanpeng, Jin Guo, Pei Ying, Guannan Zeng, and Nanyu Chen. 2022. "Nonlinear Control of a Single Tail Tilt Servomotor Tri-Rotor Ducted VTOL-UAV" Aerospace 9, no. 6: 296. https://doi.org/10.3390/aerospace9060296

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