Next Article in Journal
A Variation-Aware Design Methodology for Distributed Arithmetic
Next Article in Special Issue
Bio-Inspired Autonomous Visual Vertical and Horizontal Control of a Quadrotor Unmanned Aerial Vehicle
Previous Article in Journal
Drift-Diffusion Simulation of High-Speed Optoelectronic Devices
Previous Article in Special Issue
Preliminary Design of an Unmanned Aircraft System for Aircraft General Visual Inspection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Longitudinal Attitude Control Decoupling Algorithm Based on the Fuzzy Sliding Mode of a Coaxial-Rotor UAV

1
School of Mechatronical Engineering, Beijing Institute of Technology, Beijing 100081, China
2
China North Industries Corp., Beijing 100053, China
*
Author to whom correspondence should be addressed.
Electronics 2019, 8(1), 107; https://doi.org/10.3390/electronics8010107
Submission received: 21 December 2018 / Revised: 14 January 2019 / Accepted: 14 January 2019 / Published: 18 January 2019
(This article belongs to the Special Issue Autonomous Control of Unmanned Aerial Vehicles)

Abstract

:
A longitudinal attitude decoupling algorithm based on the fuzzy sliding mode control for a small coaxial rotor unmanned aerial vehicle (UAV) is presented in this paper. The attitude system of a small coaxial rotor UAV is characterized by nonlinearity, strong coupling and uncertainty, which causes difficulties pertaining to its flight control. According to its six-degree-of-freedom model and structural characteristics, the dynamic model was established, and a longitudinal attitude decoupling algorithm was proposed. A fuzzy sliding mode control was used to design the controller to adapt to the underactuated system. Compared with the uncoupled fuzzy sliding mode control, simulation results indicated that the proposed method could improve the stability of the system, presented with a better adapting ability, and could effectively suppress the modeling error and external interference of the coaxial rotor aircraft attitude system. The proposed method also has the advantages of high accuracy, good stability, and the ease of implementation.

1. Introduction

In recent years, due to small unmanned aerial vehicles’ (UAV) characteristics regarding maneuverability, flexibility and location difficulties, research on this type of UAV has drawn wide attention. With the unprecedented development of small aircrafts, the autonomous flight control of UAVs has become a research priority in the field of aviation [1]. Compared with fixed-wing aircrafts, the coaxial rotor uses a pair of coaxial reversing rotors which compensate for each other’s torque, instead of balancing the yaw moment of the aircraft without the tail rotor [2]. Therefore, the aircraft has a compact structure, a small radial size, and a higher power efficiency. The data indicate that it is 35–40% smaller than the single rotor structure with a tail rotor, and in the same hovering conditions, the coaxial-rotor consumes 5% less energy than the single rotor [3]. In addition, with the reduction of the radial size of the aircraft along the rotor, the inertia of the aircraft decreases and its controllability and maneuverability are enhanced. The design without the tail rotor has also eliminated some hidden problems [4]. Research on coaxial-rotor helicopters has already had significant achievements, but the small coaxial-rotor UAV has received special attention in recent years. The small size of the aircraft and the different maneuvering modes brings about differences in control methods. The operation mode of the coaxial vehicle is different from that of an ordinary vehicle which is also a typical underactuated system [5]. For the small coaxial-rotor UAV, the six-degree-of-freedom non-linear coupling problem is prominent, and decoupling is important for stability and control of the vehicle.
The general attitude control system of the UAV coordinates and controls the longitudinal, lateral and heading channels. The design of the longitudinal channel controller is the most critical and complex of the three channels, and its control rate significantly affects the UAV’s flight performance [6,7]. In the literature [8], the decoupling control method is used to design the aircraft control system. Both the adaptability and control effect of the system, however, need improvement. In Reference [9], by combining the advantages of feedback linearization and variable structure control, the attitude controller of the aircraft was designed. However, it was unable to effectively weaken the sliding mode chattering of the system, and the controller parameters could not be adjusted in real time according to the disturbance, which caused poor control performance. Reference [7] proposed a control law which was designed using the adaptive backstepping method and which did not require any knowledge of aircraft aerodynamics. Simulation results showed good performance of the feedback law, but the actual implementation was complicated and difficult to achieve. In Reference [10], a fuzzy logic control of the longitudinal motion of an aircraft based on the Takagi–Sugeno modeling approach was presented, and while the stability and tracking effect were good, the problem of system coupling had not been solved well and the control precision needed improvement. There are many studies of the decoupling controls of aircraft, but few are focused specifically on coaxial aircrafts [11].
The main role of this paper is to propose a decoupling algorithm that improves the reliability of the attitude control for the longitudinal motion stability of the coaxial rotor UAV. In order to satisfy the stability requirements of a coaxial-rotor UAV’s longitudinal motion [12], a suitable control algorithm and controller needed to be designed. Before this, we required a dynamic model which featured a qualified and effective vehicle longitudinal motion [13]. In accordance to the lab-developed coaxial rotor UAV, a rigorous and effective non-linear mathematical model of longitudinal motion was established, and an under-actuated controller was designed using the fuzzy sliding mode. Simulation results showed that the position control performance of the aircraft was improved when the decoupling algorithm was applied to the coaxial rotor longitudinal motion control system. The position and attitude were significantly improved [14] and the method was simple and effective.
This paper is organized as follows: In Section 2, according to the self-developed coaxial vehicle, the modeling and derivation processes are given. In Section 3, the decoupling algorithm design is introduced. The controller design and stability analysis based on the fuzzy sliding mode control are described in Section 4. Finally, the simulation results and comparison with the decoupling algorithm are shown in Section 5.

2. Aircraft Longitudinal Flight Model

There is a large degree of coupling among the control inputs of the aircraft. The general method is to regard these coupling quantities as external disturbances, but this method introduces large errors. To solve this problem, we used the method of controlling the correlation coefficient of the input, by selecting the appropriate correlation coefficient so that coupling among the control inputs would be handled better.

2.1. Rotor System Modeling

In order to establish a simplified model that could both reflect the aerodynamic characteristics of a coaxial-rotor and be suitable for controller design, we first made the following assumptions: The blade was rigid, the blade root truncation effect was ignored, and the tip loss and the flapping hinge extension were assumed without considering the unsteady effect [15]. The structure design and force analysis of the coaxial-rotor are shown in Figure 1.
Based on the blade-element theory, the integral expression of the rotor pulling force and torque could be obtained:
T = N b 2 π 0 2 π 0 R e ( l c v d s v ) d r d ψ f Q = N b 2 π 0 2 π 0 R e ( l c v d s v ) r d r d ψ f
where N b is the number of blades, R is the radius of the rotor, e is the amount of hinge extension of blade swing, ψ f is the local azimuth, r is the local radial coordinates, l and d respectively are local lift and resistance, and c v and s v are the correction terms related to the aircraft and the flying environment [16]. The approximate inflow ratio of a blade-element was:
λ = λ i n + λ f s = v i n Ω R v b Ω R
where v i n is the induced velocity, Ω is the rotor speed, v b is the rotor speed, is the body speed, and λ i n and λ f s are the inflow ratio corresponding to the induced velocity and body velocity, respectively.
To solve the problem of interference between coaxial rotors, the Pitt–Peters dynamic inflow method was used to model the induced velocity. The specific method was to connect the dynamic variation of the induced velocity with the variation of aerodynamic parameters through a first-order linear differential equation. We could then clarify the relationship between the pull coefficient and the induced velocity through integral calculation. The method was simple and in good agreement with the experimental data. The proposed model could be well applied to the simulation and controller designs. According to the dynamic inflow model, the relationship between the induced velocity and the pull coefficient was as follows:
M λ ˙ i n + V L 1 λ i n = C
where M , V and L are the parameter matrix of inflow dynamics, respectively, λ i n = ( λ 0 + λ s + λ c ) T represents the time-averaged, first-order horizontal, and vertical components of the induced inflow ratio. C = ( C T + C l + λ m ) T represents the pull torque, roll torque and pitch torque coefficients of the rotors. The interaction of the induced velocity is expressed as:
λ i = λ i n , i + K j i λ i n , j + λ f s e 1
where i , j are the upper and lower rotors and K j i is the parameter matrix, indicating the influence of the induced velocity between the upper and lower rotors related to the distance between the rotors, airfoil, and flight state, while e 1 = ( 1 0 0 ) T .
Considering that the interaction of the induced velocity mainly affected the channels of total distance and heading, the induced velocity in the plane of the propeller disk was almost unaffected. Assuming that the induced velocity was uniformly distributed in the plane of the propeller disk, the average inflow ratio and the differential inflow ratio were defined respectively as:
λ a = 1 2 ( λ u + λ l ) λ r = 1 2 ( λ u λ l )
From the formulas above, the pull and torque coefficients could be obtained using integral calculation ( τ is the rotor flapping time constant):
C T i = σ a 0 2 ( 1 3 θ i 1 2 ) λ i C Q i = λ i C T i + σ 8 C D
So, the single rotor thrust and torque were presented:
T i = ρ A ( Ω i R ) 2 C T i Q i = ρ A Ω i 2 R 3 C Q i
where subscript i represents the upper rotor (u) or the lower rotor (l), θ i is the pitch, c is the chord length of the blade, A is the paddle area, a 0 is the slope of the lift line of the airfoil, C D is the airfoil drag coefficient, ρ is the air density, and σ = ( N b c ) / ( π R ) is the real degree of the paddle.

2.2. External Force Modeling

In the steady state of hovering, the external force of the aircraft was determined by the lift F 1 and F 2 , the gravity G , and the aerodynamic resistance, F D :
F = F 1 + F 2 + G + F D

2.2.1. The Lift of the Rotor System

The lift of the designed aircraft was provided by the upper and lower rotors. Since the swash plate related to the upper rotor, the upper rotor provided both lift and lateral force, while the lower rotor only provided lift. The following could be obtained from Formulas (2) and (7):
F 1 = ρ A ( Ω u R ) 2 C T u cos δ F 2 = ρ A ( Ω l R ) 2 C T l

2.2.2. The Gravity of Aircraft

The mass of aircraft is m , and the body gravity was expressed by the body coordinate system data:
G = m g [ sin θ cos θ sin ϕ cos θ cos ϕ ]

2.2.3. The Aerodynamic Resistance

According to the empirical formula of aerodynamics, the resistance of the fuselage in the hovering state could be expressed as:
F D = 1 2 ρ V 2 A f u s C D f u s
where V is the relative fly-forward speed of the aircraft, A f u s is the equivalent cross-sectional area of the body, and C D f u s is the resistance coefficient of the whole body.

2.3. Establishing the Longitudinal Posture Model

{ m x ¨ = ( F 1 cos δ + F 2 F D ) sin ϕ + ε 0 F 1 d sin δ cos ϕ m z ¨ = ( F 1 cos δ + F 2 F D ) cos ϕ + ε 0 F 1 d sin δ sin ϕ m g I x x ϕ ¨ = F 1 d sin δ
Take x 1 = x , x 2 = x ˙ , z 1 = z , z 2 = z ˙ , u 1 = ( F 1 cos δ + F 2 F D ) / m , u 2 = ( F 1 d sin δ ) / I x x ,
ε = ε 0 I x x m . Formula (12) [17,18] could be written as:
{ x ˙ 1 = x 2 x ˙ 2 = u 1 sin ϕ + ε u 2 cos ϕ z ˙ 1 = z 2 z ˙ 2 = u 1 cos ϕ + ε u 2 sin ϕ g ϕ ˙ = φ φ ˙ = u 2
where x 1 ( t ) , z 1 ( t ) and ϕ ( t ) are the mass center and deflection angles of aircraft, u 1 and u 2 are the control inputs, g is the gravitational acceleration, and ε 0 represents the parasitic force of lateral displacement generated by rolling torque and ε 0 0 [19]. The system’s output was:
y ( t ) = [ x 1 ( t ) z 1 ( t ) ϕ ( t ) ] T

3. Design of the Decoupling Algorithm

Formula (13) was a strong non-linear coupling model [20]. In order to transform the model into the under-driven standard form, the following method was used to decouple the model [21,22].

3.1. Eliminate the Control Coupling of x ˙ 2 and z ˙ 2

In order to eliminate the control coupling of x ˙ 2 and z ˙ 2 , the following decoupling algorithm could be designed:
[ u 1 u 2 ] = [ sin ϕ ε cos ϕ cos ϕ ε sin ϕ ] 1 [ u m 1 u m 2 + g ]
where u m 1 and u m 2 are the control items to be designed. Next, Formula (13) could be changed:
{ x ¨ 1 = u m 1 z ¨ 1 = u m 2 ε ϕ ¨ = u m 1 cos ϕ + ( u m 2 + g ) sin ϕ

3.2. Eliminate the Coupling of u m 1 and u m 2

In order to eliminate the coupling of u m 1 and u m 2 in ε ϕ ¨ , the following decoupling algorithm could be designed:
{ x m = x 1 ε sin ϕ z m = z 1 + ε cos ϕ u m 1 = ( u m 3 ε ϕ ˙ 2 ) sin ϕ + ε u m 4 cos ϕ u m 2 = ( u m 3 ε ϕ ˙ 2 ) cos ϕ + ε u m 4 sin ϕ g
where u m 3 and u m 4 are the control items to be designed, then Formula (16) could be changed:
{ x ¨ m = u m 3 sin ϕ z ¨ m = u m 3 cos ϕ g ϕ ¨ = u m 4

3.3. Eliminate the Coupling of sin ϕ and cos ϕ

In order to eliminate the coupling of sin ϕ and cos ϕ , the following decoupling algorithm was designed:
We took q 1 = x m , q 2 = q ˙ 1 , q 3 = z m , q 4 = q ˙ 3 , q 5 = tan ϕ , and q 6 = q ˙ 5 , and when ϕ was very small, we took sin ϕ ϕ . Then we could obtain the following:
q ¨ 5 = u m 4 cos 2 ϕ + 2 ϕ ˙ 2 cos 3 ϕ sin ϕ cos 4 ϕ
We took h 1 = u m 3 cos ϕ , h 2 = tan ϕ . Then we could get u m 4 = h 2 cos 2 ϕ 2 ϕ ˙ 2 tan ϕ . Now, h 1 and h 2 were the control items to be designed. Then, Formula (18) could be changed as:
{ q ˙ 1 = q 2 q ˙ 2 = q 5 h 1 q ˙ 3 = q 4 { q ˙ 4 = h 1 g q ˙ 5 = q 6 q ˙ 6 = h 2

3.4. Transform the Control Model into an Under-Driven Standard Form

The under-driven standard form could be obtained:
{ y ˙ 1 = y 2 y ˙ 2 = f 1 ( y 1 y 2 y 3 ) y ˙ 3 = y 4 y ˙ 4 = f 2 ( y 1 y 2 y 3 ) + b h + d
where f 2 ( y 1 y 2 y 3 ) = [ g 0 ] , b = [ 1 0 0 1 ] , h = [ h 1 h 2 ] , and d is the control disturbance. In order to make f 1 ( y 1 y 2 y 3 ) y 3 invertible, which would be helpful for the design of the control law, we took y 2 = [ q 2 + q 4 q 5 t 0 q 3 d t ] , y 1 = [ q 1 + t 0 q 4 q 5 d t t 0 ( t 0 q 3 d t ) d t ] , y 3 = [ q 3 q 5 ] , y 4 = [ q 4 q 6 ] , and f 1 ( y 1 y 2 y 3 ) = [ g y 3 ( 2 ) y 3 ( 1 ) ] .

4. Design of the Controller

The structure of the controller [16] is shown in Figure 2.

4.1. Design of the Control Law

In order to design the control law for Formula (21), we took y 1 d as the reference instruction of y 1 , and the error variables were as follows [19]:
e 1 = y 1 y 1 d        ,                    e 2 = y 2 y ˙ 1 d e 3 = f 1 y ¨ 1 d        ,                     e 4 = f 1 y 3 y 4 y 1 d ( 3 )
Next, we designed the sliding surface:
s = c 1 e 1 + c 2 e 2 + c 3 e 3 + e 4
where c i > 0 and i = 1 , 2 , 3 .
When s ˙ = 0 , we could obtain h = u e q and c 1 e ˙ 1 + c 2 e ˙ 2 + c 3 e ˙ 3 + e ˙ 4 = 0 .
The equivalent switch control items could thus be obtained [23]:
u e q = ( f 1 y 3 b ) 1 ( c 1 y 2 c 1 y ˙ 1 d                             + c 2 f 1 c 2 y ¨ 1 d    + c 3 f 1 y 3 y 4                              c 3 y 1 d ( 3 ) + f 1 y 3 f 2 y 1 d ( 4 ) )
where y ˙ 1 = y 2 , y ¨ 1 = f 1 , y 1 ( 3 ) = f ˙ 1 = f 1 y 3 y 4 and y 1 ( 4 ) = f ¨ 1 = f 1 y 3 y ˙ 4
The switching control item was then designed:
u s w 1 = ( f 1 y 3 b ) 1 [ K ( t ) sgn ( s ) + λ s + E 1 ( t ) ]
where E 1 ( t ) is unknown interference. Both E 1 ( t ) and K ( t ) will be described in more detail in Section 4.2. The control law could be designed as follows:
h = u e q + u s w

4.2. Stability Analysis of the Control System

Taking Formulas (23)–(25) into s ˙ , the following could be obtained:
s ˙ = c 1 e ˙ 1 + c 2 e ˙ 2 + c 3 e ˙ 3 + e ˙ 4         = c 1 ( y 2 y ˙ 1 d ) + c 2 ( f 1 y ¨ 1 d ) + c 3 ( f 1 y 3 y 4 y 1 d ( 3 ) )                 + d d t [ f 1 y 3 ] y 4 + f 1 y 3 f 2 y 1 d ( 4 )         = K ( t ) sgn ( s ) λ s + f 1 y 3 d
We took K ( t ) = α d ¯ + ρ , where d ¯ ( 1 ) | d ( 1 ) | , d ¯ ( 2 ) | d ( 2 ) | , ρ ( 1 ) > 0 , ρ ( 2 ) > 0 , and α > 0 .
We took the Lyapunov function as V = 1 2 s T s , so:
V ˙ = s T s ˙ = s T [ ( α d ¯ + ρ ) sgn ( s ) λ s + f 1 y 3 d ]    = ( α d ¯ + ρ ) s λ s 2 + s T E ( t )    ρ s λ s 2 0
where the gain of the switching K ( t ) was the cause of chattering, and the control disturbance can be expressed as E ( t ) = f 1 y 3 d , which was used for ensuring that the necessary sliding mode presence conditions were met. When s = 0 , we could obtain e 4 = c 1 e 1 + c 2 e 2 + c 3 e 3 . We took the following:
A = [ 0 1 0 0 0 1 c 1 c 2 c 3 ]
A is the Hurwitz function, and λ represents the Eigenvalues of A , λ > 0 .
Taking E 1 = [ e 1            e 2 e 3 ] T , the error equation of the state could be written as E ˙ 1 = A E 1 .
Taking Q = Q T > 0 , we could get the Lyapunov equation A T P + P A = Q . The solution was P = P T > 0 . We took the Lyapunov function as the following:
V ˙ 1 = E ˙ 1 T P E 1 + E 1 T P E ˙ 1 = ( A E 1 ) T P E 1 + E 1 T P ( A E 1 )              = E 1 T A T P E 1 + E 1 T P A E 1 = E 1 T ( A T P + P A ) E 1              = E 1 T Q E 1 λ min ( Q ) | | E 1 | | 2 2 0
where λ min ( Q ) is the minimum eigenvalue of positive definite matrix, Q .
From V ˙ 1 0 , we could obtain: e 1 0 , e 2 0 , e 2 0 , then y 1 y 1 d , y 2 y 2 d , and y 3 y 3 d . From the stability of the sliding mode, we could obtain y 4 y 4 d . In the end, x x d , z z d , and ϕ ϕ d .

4.3. Establish the Fuzzy System

The condition for the existence of sliding mode was s T s ˙ < 0 , and when the system reached the sliding surface, it would remain on the sliding surface [24]. From Formula (26), we could see that in order to ensure that the system movement reached the gain of the sliding surface, K ( t ) needed to be sufficient to eliminate the impact of uncertainty. Then we could ensure the existence of the sliding condition.
The idea of the fuzzy rules was represented as follows:
If s T s ˙ > 0 , then K ( t ) should increase;
if s T s ˙ < 0 , then K ( t ) should be reduced.
From the two types above, we could design the fuzzy system using s T s ˙ and Δ K ( t ) . In this system, s T s ˙ is the input, and Δ K ( t ) is the output. The fuzzy sets of the system were defined as follows:
s T s ˙ = { N B N M Z O P M P B } Δ K = { N B N M Z O P M P B }
where N B represents the negative big, N M is the negative middle, Z O is the zero, P M is the positive middle, and P B is the positive big. The input and output membership functions of the fuzzy system are shown in Figure 3 and Figure 4. The upper bound of K ^ ( t ) was estimated using the integral method:
K ^ ( t ) = G 1 0 t Δ K d t
where G 1 is the proportion coefficient, determined according to experience. The control law was:
u s w 2 = ( f 1 y 3 b ) 1 [ K ^ ( t ) sgn ( s ) + λ s + E ( t ) ]

5. Simulation Analysis

The physical parameters of the system model were obtained through the three-dimensional model established in the CAD software Inventor, and the model reference coefficients were calculated in accordance to the physical parameters and the modeling results. The parameters are shown in Table 1.
For the controlled Formula (13), we took ε = 10 and g = 9.8 , and set a predetermined track as x d = t ,     z d = sin t and ϕ d = 0 .
In order to make A become the Hurwitz function, we took the control law parameters c 1 = 27 ,     c 2 = 27 ,     c 3 = 9 and λ = 0.10 . The initial state of the controlled system was taken as [ 5 0 0.5 0 0.1 0 ] . We used the control law (Formula (25)) and saturation function method, and took the thickness of the boundary layer Δ to be 0.10 .
According to the structural characteristics of the coaxial-rotor UAV, a dynamic model of the longitudinal motion was established. The dynamic model of the aircraft was then decoupled, the fuzzy control and sliding mode controls were combined, and then a fuzzy sliding mode control based on the decoupling algorithm was designed for the coaxial-rotor. The control method was then simulated by MATLAB/Simulink. The results showed that the control method could track the command signal more quickly and efficiently compared to the method of the traditional sliding mode control. It could quickly reduce the yaw attitude angle deviation and the steady-state error could reach almost zero, and with a strong self-adaptive ability, it could achieve a better control effect. The response speed, tracking accuracy, and efficiency of the system were significantly improved.
The proposed control method could improve the stability of the system, which could effectively restrain the modeling errors and external disturbances of the aircraft’s attitude system. This method had the advantages of high control precision, strong robustness, and ease of implementation in engineering. In future studies, we will focus on the design of the decoupling algorithm under the influence of more inputs and interferences, and will apply this algorithm to specific engineering practices.
Figure 5 and Figure 6 show the performance of position tracking in the horizontal direction while the two control methods were used. The instruction given along x was a straight-line motion. The former figure indicates the position tracking with the decoupling algorithm and fuzzy control. The latter indicates the fuzzy control without the decoupling algorithm. From these figures, we could see that the performance of the control method with the decoupling algorithm and fuzzy control was faster, more accurate, and more stable than the general sliding mode control, ensuring that the aircraft would be more stable in actual movement and in improving the flight.
Figure 7 and Figure 8 show the performance of the position tracking in the vertical direction while the two control methods were used. The instruction given along z was a sinusoidal motion. The former figure indicates the position tracking with the decoupling algorithm and fuzzy control, and the latter indicates the general sliding mode control. From these figures we could see that the time required for the two methods to track from the initial position to the specified trajectory was almost the same, but the performance of the control method with the decoupling algorithm and fuzzy control was smoother, and the tracking error was also smaller. Thus, the flight of the aircraft would be more stable.
The angle tracking with the decoupling algorithm and fuzzy control are shown in Figure 9, and when compared with the angle tracking (Figure 10) without the decoupling algorithm and fuzzy control, it could be seen that in the former the tracking errors decreased while the response times were basically the same, and the movement accuracy and resistance to disturbances of the system improved. The system had better tracking performance. The impact of these disturbances, controller outputs chattering, external disturbances, and noise of the measurement were significantly reduced after the decoupling algorithm and fuzzy controller are added.
Figure 11, Figure 12 and Figure 13 show the speed tracking and angular speed tracking with the decoupling algorithm and fuzzy control. It could be seen that the change trajectory was smooth, and the time required to reach a stable state from the initial state was very short. Figure 14 and Figure 15 show the control input of the system. We could see that the control input curves were smooth and no chattering occurred. Theoretical analysis and experimental simulation results showed that the fuzzy sliding mode control based on the decoupling algorithm could improve the stability of the system, had a better self-adaptive ability, and effectively restrained the modeling errors and external disturbances of the co-rotating twin-rotor aircraft attitude system.

6. Conclusions

According to the structural characteristics of the coaxial-rotor UAV, a dynamic model of longitudinal motion was established. The dynamic model of the aircraft was then decoupled, the fuzzy control and sliding mode controls were combined, and a fuzzy sliding mode control based on the decoupling algorithm was designed for the coaxial-rotor. The control method was then simulated by MATLAB/Simulink. The results showed that the control method could track the command signal more quickly and efficiently compared to the method of the traditional sliding mode control. It could quickly reduce the yaw attitude angle deviation and the steady-state error could reach almost zero. With a strong self-adaptive ability, it could achieve a better control effect. The response speed, tracking accuracy, and efficiency of the system had been significantly improved.
The proposed control method could improve the stability of the system, which could effectively restrain the modeling errors and external disturbances of the aircraft’s attitude system. This method had the advantages of high control precision, strong robustness, and ease of implementation in engineering. In future studies, we will focus on the design of the decoupling algorithm under the influence of more inputs and interferences, and will apply this algorithm to specific engineering practices.

Author Contributions

K.L. designed the main parts of the method and was also responsible for writing the paper. Y.W. built the Simulation program. C.W. and H.D. provided some technical comments.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kannan, N.; Bhat, M.S. Longitudinal H Infinity Stability Augmentation System for a Thrust Vectored Unmanned Aircraft. J. Guid. Control Dyn. 2005, 28, 1240–1250. [Google Scholar] [CrossRef]
  2. Krämer, P.; Gimonet, B.; v Grünhagen, W. A systematic approach to nonlinear rotorcraft model identification. Aerosp. Sci. Technol. 2002, 6, 579–590. [Google Scholar] [CrossRef]
  3. Wang, X.; Li, K.; Zhao, N.; Deng, H. Nonlinear dynamics modeling and simulation of cylindrical coaxial UAV. In Proceedings of the IEEE International Conference on Real-Time Computing & Robotics, Okinawa, Japan, 14–18 July 2017; pp. 505–510. [Google Scholar]
  4. Niu, S.; Li, J.; Shen, Y. Design, modeling and disturbance rejection control of a bio-inspired coaxial helicopter MAV in Atmospheric Boundary Layer. In Proceedings of the 2015 IEEE International Conference on Robotics and Biomimetics (ROBIO), Zhuhai, China, 6–9 December 2016; pp. 1272–1277. [Google Scholar]
  5. Mettler, B.; Tischler, M.B.; Kanade, T. System Identification Modeling of a Small-Scale Unmanned Rotorcraft for Flight Control Design. J. Am. Helicopter Soc. 2002, 47, 50–63. [Google Scholar] [CrossRef]
  6. Khandani, K.; Majd, V.J.; Darestani, M.R.; Talebi, H. A stochastic sliding mode scheme for longitudinal control of an aircraft model. In Proceedings of the 25th Iranian Conference on Electrical Engineering, Tehran, Iran, 2–4 May 2017; pp. 795–800. [Google Scholar]
  7. Gavilan, F.; Vazquez, R.; Acosta, J.Á. Adaptive control for aircraft longitudinal dynamics with thrust saturation. J. Guid. Control Dyn. 2015, 38, 651–661. [Google Scholar] [CrossRef]
  8. Duan, F.H.; Han, C.Z. Application of Nonlinear System Decoupling Control Theory in Flight Control. Electron. Opt. Control 2001, 84, 30–34. [Google Scholar]
  9. Ao, B.Q.; Li, J.L. Variable-structure control method of aircraft attitude based on feedback linearization. Mod. Def. Technol. 2003, 31, 41–44. [Google Scholar]
  10. Petr, H.; Kashyapa, N. Aircraft longitudinal motion control based on Takagi–Sugeno fuzzy model. Appl. Soft Comput. 2016, 49, 269–278. [Google Scholar]
  11. BC Campos, L.M. Nonlinear Longitudinal Stability of a Symmetric Aircraft. J. Aircr. 1997, 34, 360–369. [Google Scholar] [CrossRef]
  12. McLean, D.; Zaludin, Z.A. Stabilization of longitudinal motion of a hypersonic transport aircraft. Trans. Inst. Meas. Control 1999, 21, 99–105. [Google Scholar] [CrossRef]
  13. Jyothi, J.; Bindu, G.R.; Jayakumar, M. Robust Longitudinal Controller Design for an Unmanned Tailless Aircraft. In Proceedings of the 6th AIAA Aviation Technology, Integration and Operations Conference, Wichita, KS, USA, 25–27 September 2006; pp. 1–7. [Google Scholar]
  14. Yue, T.; Wang, L.X.; Ai, J.Q. Longitudinal Linear Parameter Varying Modeling and Simulation of Morphing Aircraft. J. Aircr. 2013, 50, 1673–1681. [Google Scholar] [CrossRef]
  15. Adam, M.W.; Ephrahim, G. Longitudinal Dynamics of a Perching Aircraft. Proc. SPIE Int. Soc. Opt. Eng. 2012, 43, 1386–1392. [Google Scholar]
  16. Yuan, X.M.; Zhu, J.H.; Mao, M. Modeling and robust tracking control for coaxial unmanned helicopter. Control Theory Appl. 2014, 31, 1286–1294. [Google Scholar]
  17. Yang, X.M.; Li, W.J. Four rotor aircraft control based on sliding mode controller. J. Univ. Technol. Nat. Sci. 2016, 39, 924–928. [Google Scholar]
  18. Shiyang, R.; Dominick, A. Longitudinal Flying Qualities Prediction for Nonlinear Aircraft. J. Guid. Control Dyn. 2003, 26, 474–482. [Google Scholar]
  19. Al-Hiddabi, S.A.; McClamroch, N.H. Aggressive Longitudinal Aircraft Path Tracking Using Nonlinear Control. Asian J. Control 2010, 3, 280–288. [Google Scholar] [CrossRef]
  20. Speyer, J.; White, J.; Douglas, R.; Hull, D. MIMO controller design for longitudinal decoupled aircraft motion. In Proceedings of the Guidance and Control Conference, Boston, MA, USA, 19–22 August 2013; pp. 729–737. [Google Scholar]
  21. Shan, S.Q.; Hou, Z.X.; Wang, W.K. Aircraft longitudinal decoupling based on a singular perturbation approach. Adv. Mech. Eng. 2017, 9, 1–8. [Google Scholar] [CrossRef]
  22. Phillips, W.F.; Santana, B.W. Aircraft Small-Disturbance Theory with Longitudinal-Lateral Coupling. J. Aircr. 2002, 39, 973–980. [Google Scholar] [CrossRef]
  23. Xu, R.; Ozguner, U. Sliding mode control of a class of underactuated systems. Automatica 2008, 44, 233–241. [Google Scholar] [CrossRef]
  24. Wang, Y.C.; Sun, H. Sliding Mode Controller Design for VTOL Aircraft. Comput. Meas. Control 2016, 6, 102–105. [Google Scholar]
Figure 1. (a) The structure design and (b) force analysis of the coaxial-rotor Rotorcraft.
Figure 1. (a) The structure design and (b) force analysis of the coaxial-rotor Rotorcraft.
Electronics 08 00107 g001
Figure 2. The structure of the controller.
Figure 2. The structure of the controller.
Electronics 08 00107 g002
Figure 3. The input membership function.
Figure 3. The input membership function.
Electronics 08 00107 g003
Figure 4. The output membership function.
Figure 4. The output membership function.
Electronics 08 00107 g004
Figure 5. Position tracking with the decoupling algorithm of the x-axis.
Figure 5. Position tracking with the decoupling algorithm of the x-axis.
Electronics 08 00107 g005
Figure 6. Position tracking without the decoupling algorithm of the x-axis.
Figure 6. Position tracking without the decoupling algorithm of the x-axis.
Electronics 08 00107 g006
Figure 7. Position tracking with the decoupling algorithm of the z-axis.
Figure 7. Position tracking with the decoupling algorithm of the z-axis.
Electronics 08 00107 g007
Figure 8. Position tracking without the decoupling algorithm of the z-axis.
Figure 8. Position tracking without the decoupling algorithm of the z-axis.
Electronics 08 00107 g008
Figure 9. Angle tracking with the decoupling algorithm.
Figure 9. Angle tracking with the decoupling algorithm.
Electronics 08 00107 g009
Figure 10. Angle tracking without the decoupling algorithm.
Figure 10. Angle tracking without the decoupling algorithm.
Electronics 08 00107 g010
Figure 11. Speed tracking with the decoupling algorithm.
Figure 11. Speed tracking with the decoupling algorithm.
Electronics 08 00107 g011
Figure 12. Speed tracking with the decoupling algorithm.
Figure 12. Speed tracking with the decoupling algorithm.
Electronics 08 00107 g012
Figure 13. Position tracking with the decoupling algorithm.
Figure 13. Position tracking with the decoupling algorithm.
Electronics 08 00107 g013
Figure 14. Control input u 1 with the decoupling algorithm.
Figure 14. Control input u 1 with the decoupling algorithm.
Electronics 08 00107 g014
Figure 15. Control input u 2 with the decoupling algorithm.
Figure 15. Control input u 2 with the decoupling algorithm.
Electronics 08 00107 g015
Table 1. The main symbols and parameters.
Table 1. The main symbols and parameters.
ParameterDescriptionValueUnit
ρ Local air density 1.14 kg / cm 2
g Local gravitational acceleration 9.804 m / s 2
R Rotor radius 0.20 m
A Area of rotor plane 0.126 m 2
m Weight 1 kg
d Distance between upper rotor and the gravity center 0.20 m
e Amount of hinge extension of blade flapping 0.05 m
I x x Inertia about x-axis 9.16 × 10 4 kg / m 2
C T u Trust coefficient of upper rotor 9.42 × 10 3 rad 1
C T l Trust coefficient of lower rotor 6.77 × 10 3 rad 1
C Q u Torque coefficient of upper rotor 6.14 × 10 4 rad 1
C Q l Torque coefficient of lower rotor 6.01 × 10 4 rad 1

Share and Cite

MDPI and ACS Style

Li, K.; Wei, Y.; Wang, C.; Deng, H. Longitudinal Attitude Control Decoupling Algorithm Based on the Fuzzy Sliding Mode of a Coaxial-Rotor UAV. Electronics 2019, 8, 107. https://doi.org/10.3390/electronics8010107

AMA Style

Li K, Wei Y, Wang C, Deng H. Longitudinal Attitude Control Decoupling Algorithm Based on the Fuzzy Sliding Mode of a Coaxial-Rotor UAV. Electronics. 2019; 8(1):107. https://doi.org/10.3390/electronics8010107

Chicago/Turabian Style

Li, Kewei, Yiran Wei, Chao Wang, and Hongbin Deng. 2019. "Longitudinal Attitude Control Decoupling Algorithm Based on the Fuzzy Sliding Mode of a Coaxial-Rotor UAV" Electronics 8, no. 1: 107. https://doi.org/10.3390/electronics8010107

APA Style

Li, K., Wei, Y., Wang, C., & Deng, H. (2019). Longitudinal Attitude Control Decoupling Algorithm Based on the Fuzzy Sliding Mode of a Coaxial-Rotor UAV. Electronics, 8(1), 107. https://doi.org/10.3390/electronics8010107

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