Next Article in Journal
Developing an Image Manipulation Detection Algorithm Based on Edge Detection and Faster R-CNN
Previous Article in Journal
Greedy Method for Boarding a Partially Occupied Airplane Using Apron Buses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vector-Coupled Flight Controller Design Based on Multivariable Backstepping Sliding Mode

College of Intelligent Science and Technology, National University of Defense Technology, Changsha 410073, China
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(10), 1225; https://doi.org/10.3390/sym11101225
Submission received: 25 August 2019 / Revised: 23 September 2019 / Accepted: 24 September 2019 / Published: 1 October 2019

Abstract

:
This study aims to address the issues of the simultaneous control of the angle of attack, sideslip angle, and airspeed of a flight in a vectorial form. A vector controller, with a symmetry structure, is developed to transform the attitude and speed control problem into a space-vector tracking problem. We first establish flight vector-coupled dynamics, i.e., describe velocity and angular-velocity vectors in a body-fixed frame, and then propose a multivariable backstepping sliding mode control algorithm along with nonlinear disturbance observers for the vectorial dynamics. The theoretical analysis ensures that the states of the system can be enforced to reach a small neighborhood of the desired sliding manifold. The results of the numerical simulation illustrate the effectiveness and robustness of the combined vector-control scheme.

1. Introduction

There are many interesting and challenging problems with respect to the control of modern aircraft. These include, but are not limited to, the cross-coupling of dynamics [1], non-minimum phase behavior [2], constraints on the control inputs as a result of position and rate limits [3], and high nonlinearity and uncertainty [4] in flight equations. Moreover, the performance of a flight control system is expected to be very high, and the demands for flight safety have dictated the development of the robust control designs [5].
Traditionally, flight control systems were developed based on the “small perturbation theory”. This theory follows the assumption that flight dynamics are time-invariant and linear around the operating points, while the longitudinal dynamics are independent of the lateral ones. Then, the equations of motion can be decoupled and addressed separately with linear control methods such as proportion integration differentiation (PID) control [6] and linear quadratic control [7,8]. However, the performance of flight control systems deteriorates owing to the un-modeled terms of the nonlinearities intrinsic to flight dynamics, and the cross-coupling between the longitudinal and lateral motions, which become crucial under an ultra condition, such as a large angle of attack (AOA) or speedy rolling [9].
To overcome the disadvantages of the linear methods, nonlinear control laws have been developed rapidly in the past decades. Feedback linearization is one of the most widely used amongst such methods [10,11]. However, the primary drawback of these methods is that all nonlinearities of a control system must be known for feedback linearization [12]. However, the effective and robust nonlinear control methods that do not require all nonlinearities to be accurate have been used in flight control. Sliding mode control [13,14,15], disturbance observer-enhanced control [16,17], model predictive control [18], adaptive backstepping [19], and dynamic surface control [20] are some nonlinearities that have been considered.
The design of controls of flights have achieved various degrees of success. However, there are several challenges that deserve further research:
  • Cross-coupling that is inherent to flight dynamics. In previous studies [13,14,16,17,18,19,20,21], the models used to design flight control were decoupled into a scalar form, i.e., state equations with velocity, AOA, sideslip angle, roll rate, pitch rate, and yaw rate, and multiple single-channel controllers were designed. During the decoupling process, some coupling components terms are first calculated/estimated and then compensated for, while the other coupling terms may be disregarded owing to their negligible influence. However, this is a conservative method because the coupling terms may include the control input.
  • Design of control laws that ensure the stability of the whole control system. In some studies [22,23], the design of the control systems is based on the assumption with “timescale separation”, wherein the slow attitude dynamics are separated from the fast angular-rate dynamics. The outer and inner controllers correspond to the slow and fast subsystems, respectively. They can be designed individually to simplify the complexity of the control system. However, the inherent weakness of this framework is that the overall system’s stability cannot be achieved theoretically.
  • Uncertainties regarding various disturbances, including external disturbances, un-modeled dynamics, and parameter uncertainties. Balancing the robustness and the control performance of an aircraft control system during the design process is important. The robustness and control performance of a control system have an inverse relationship, i.e., a more robust controller implies more certain attenuation; however, the robustness is obtained by sacrificing the nominal control performance to a certain extent.
Motivated by the aforementioned challenge, and for a multivariable backstepping sliding mode approach through disturbance estimation for the flight control systems within multivariable design, this study aimed to make the following contributions:
  • We model the flight dynamics using a body-fixed frame in the vectorial form and consider both matched and unmatched disturbances. The attitude control of an aircraft, along with the total velocity, can be converted into a space-vector tracking equation. The triplet airspeed ( V T ), AOA ( α ) and sideslip angle ( β ) are controlled simultaneously in a vectorial manner. This also allows us to deal with the cross-term actively and correctly. A key aspect of this study is the active use of the cross-coupling in terms of flight dynamics instead of decoupling and passive suppression and compensation.
  • The control-oriented model has a lower-triangular form. Then, the structure of a Lyapunov-based backstepping approach [24,25] is used in this study to ensure the stability of the closed-loop system theoretically. This work establishes a combined multivariable backstepping sliding mode controller, along with nonlinear disturbance observers. The disturbance observers enhance the scheme of construction of the control law by combining the backstepping sliding mode control feedback with disturbance estimation, based on straightforward and feedforward compensation. Unlike similar works [12,16,17,18,19], the developed controller not only avoids solving complex matrix equations and treating inverse matrixes, but also fully realizes vector-coupled control.
  • The developed control structure is concise and aesthetically appealing compared with traditional control structures that use a decoupled collection of single variables. The combined control scheme has a symmetry structure and each term is meaningful, and this feature is significant in that the control parameters can be adjusted in each term directively.
The remainder of this paper is organized as follows. The control-oriented model of the flight dynamics, established in a vectorial manner, is given in Section 2. The multivariable backstepping sliding mode controller, along with the nonlinear disturbance observers, is proposed based on the vectorial scheme in Section 3. The performance of the developed control algorithm is illustrated through simulations in Section 4, and conclusions are drawn in Section 5.

2. Modelling

To describe the flight motion in this study, three coordinates need to be defined as follows:
  • Body-fixed Frame: Σ b , the reference frame with the origin at the gravity center and axes pointing forward, over the right wing, and down (relative to the pilot).
  • Inertial Frame: Σ i , the reference frame with a specific ground origin and axes pointing the North, East, and down to the Earth center.
  • Wind Frame: Σ w , the reference frame with the origin at the gravity center and the x-axis pointing to the velocity direction of the flight. The orientation of this frame relative to the body-fixed frame is determined by AOA ( α ) and sideslip angle ( β ). The lift, drag, and side forces are defined naturally in this reference frame, respectively.
Figure 1 depicts an aircraft model [26] with the defined frames.The controlled aircraft is assumed to be a rigid body with mass m, and the motion in a body-fixed coordinate system with the origin at the center of mass, which can be described as:
F b = d d t m V b + W × m V ,
where V is the velocity vector and W is the angular-velocity vector. We can express the velocity vector as V = V ( cos α cos β sin β sin α cos β ) T = V V 0 , where V = V is the amplitude of the velocity V , and V 0 = ( cos α cos β sin β sin α cos β ) T is the unit velocity vector, α and β represent AOA and sideslip angle, respectively. F b = G + R + T b , where G , R , and T represent the gravity vector, the aerodynamics-force vector, and the thrust vector, respectively, they can be found in the Appendix A.
Moreover, the moment equation in the body-fixed coordinates can be described as
J W ˙ = W × J W + M + M d ,
where M = m x m y m z T represents the control torque vector, W = p q r T is the angular-velocity vector in the body-fixed coordinates, J is the inertia moment, and M d is the disturbance-torque vector.
By combining Equation (1) with Equation (2) and considering the uncertainties, we obtain
V ˙ = W × V + 1 m F a + u V V 0 + D 1 , W ˙ = J 1 W × JW + J 1 M + D 2 ,
where F a = G + R , and the lumped uncertainties, D 1 and D 2 , represent the discrepancies between the actual aircraft and its mathematical model used for the design of controllers. These discrepancies arise from unknown external disturbances, aircraft parameters, and unmodeled dynamics.

3. Vector-Coupled Flight Controller Design

In this section, a multivariable backstepping sliding mode controller together with a nonlinear disturbance observer is proposed for the vector-coupled flight model (Eqaution (3)) so that the velocity vector V tracks the desired vector, V d . First, the control objective is established, and then the controller is derived. Figure 2 shows a schematic of the control structure.

3.1. Control Objective

Our aim is to establish a multivariable control scheme with the following properties of the flight control system.
(1) The attitude commands, α d , β d , and the airspeed command of V d are tracked simultaneously by the system’s output vector V because these commands can be transformed into a vectorial command V d = V d ( cos α d cos β d sin β d sin α d cos β d ) T .
(2) The performance of the nominal closed-loop control system is recovered as if the lumped disturbances, D 1 and D 2 , have been known exactly.
(3) The stability of the overall system can be guaranteed through a theoretical perspective.

3.2. Control Law Design

For convenience in deducing the control algorithm, an angular-acceleration vector is introduced as
u T = J 1 ( M W × J W ) ,
and u T is the control input of the system. Then, define F a V = F a / m and the dynamics can be rewritten by
V ˙ = W × V + F a ( V ) + u V V 0 + D 1 , W ˙ = u T + D 2 .
Remark 1.
The devices onboard to control a flight include the deflections of the aerodynamic control surfaces δ and the engine thrust T . The vector δ contains the positions of canard wings δ c , right elevon δ r e , left elevon δ l e , and rudder δ r , respectively. Furthermore, the engine thrust produces an acceleration of the flight along the x-axis. In this study, for the convenience of developing the vectorial controller, the control-oriental model is established and control inputs are the angular-acceleration vector u T = J 1 ( M W × ( JW ) ) and the engine thrust term u V , which is assumed to be along with the velocity direction. In the angular-acceleration term, the controlled angular acceleration M = q p S b C l c C M b C n T is obtained by the deflections of the aerodynamic control surfaces, where q p = 0.5 ρ V T 2 is the dynamic pressure, ρ is the air density, V T is the velocity, S is the wing surface area, b and c are the wingspan and the mean aerodynamic chord, respectively. C l , C M , and C n represent the dimensionless coefficients in the moment expressions. If we can obtain the control torque vector M = J u T + W × JW , δ can be obtained by δ = B 1 J u T + W × JW , where M = B δ , and B is the distribution matrix.
Assumption 1.
Consider the lumped disturbances of acceleration and angular acceleration, D 1 and D 2 , in Equation (5), their derivatives are bounded by δ 1 = sup t > 0 D ˙ 1 and δ 2 = sup t > 0 D ˙ 2 , where δ 1 and δ 2 are positive constants [4,14,27,28,29].
Remark 2.
The lumped disturbances D 1 and D 2 in Equation (5) arise from the external disturbances and un-modeled dynamics. However, they do not only depend on time, but also on the states of the flight. Considering that the flight conditions and the states may not change rapidly, we propose the above assumption [4,14,27,28,29] for the controller design in this study. In practice, because the derivatives of these uncertainties are difficult to acquire beforehand, their upper bounds, δ 1 and δ 2 , are set relatively large to adapt with extreme flying conditions. Our future research will focus on designing an adaptive-gain controller to avoid such overshooting phenomenon.

3.2.1. Nonlinear Disturbance Observer

The uncertainties, D 1 and D 2 , are first estimated through nonlinear disturbance observers.
Inspired by the work of Chen [30], a vectorial disturbance observer D ^ 1 is designed to estimate the uncertain term D 1 , which is given by
D 1 = z 1 + λ 1 e 1 , D ^ 1 = z ^ 1 + λ 1 e 1 , z ^ ˙ 1 = λ 1 ( W × V + z ^ 1 + λ 1 e 1 + u V V 0 ) . .
Then, the time derivative of D ^ 1 is
D ^ ˙ 1 = z ^ ˙ 1 + λ 1 e ˙ 1 = λ 1 ( W × V + z ^ 1 + λ 1 e 1 + u V V 0 ) + λ 1 ( W × V + D 1 + u V V 0 ) = λ 1 ( D ^ 1 + D 1 ) = λ 1 D ˜ 1 .
Remark 3.
The nonlinear disturbance observer’s error dynamics can be described as D ˜ ˙ 1 = D ˙ 1 D ^ ˙ 1 = D ˙ 1 λ 1 D ˜ 1 , which implies D ˜ 1 s = 1 s + λ 1 D ˙ 1 s , where s is the Laplace operator. Hence, the nonlinear disturbance observer’s error dynamics can be approximated by a first-order system, which can be seen as a first-order low-pass filter when the parameter λ 1 is selected appropriately. The nonlinear disturbance observer can handle the high-frequency noise of the disturbance differential D ˙ 1 to a certain extent. This property can be verified through Case 2 in the simulation section.
Proof. 
The stability of the above nonlinear-disturbance-observer system is analyzed by introducing a Lyapunov candidate:
L D 1 = 0.5 D ˜ 1 2 , D ˜ 1 = D 1 D ^ 1 ,
and differentiating the expression in Equation (8) yields
L ˙ D 1 = D ˜ 1 T D ˜ ˙ 1 = D ˜ 1 T ( D ˙ 1 D ^ ˙ 1 ) = D ˜ 1 T ( D ˙ 1 λ 1 D ˜ 1 ) 0.5 D ˜ 1 2 + 0.5 D ˙ 1 2 λ 1 D ˜ 1 2 ( λ 1 0.5 ) D ˜ 1 2 + 0.5 δ 1 2 .
It can be found, under the condition
λ 1 > 0.5 , D ˜ 1 > δ 1 2 ( λ 1 0.5 ) ,
that
L ˙ D 1 0 .
Therefore, the estimation error of the disturbance observer in Equation (8) is bounded by the following region
S ( D ˜ 1 ) = D ˜ 1 δ 1 2 ( λ 1 0.5 ) .
Similarly, D ^ 2 is the disturbance-observer vector of D 2 which is expressed by
D 2 = z 2 + λ 2 e 2 , D ^ 2 = z ^ 2 + λ 2 e 2 , z ^ ˙ 2 = λ 2 ( u T + z ^ 2 + λ 2 e 2 ) .
Moreover, the proof procedure of the above nonlinear disturbance observation of D 2 is akin to that of D 1 . For the sake of space, the details are omitted herein. □

3.2.2. Multivariable Backstepping Sliding Mode Controller

Theorem 1.
Considering the control-oriented model of Equation (5), then multivariable backstepping sliding mode controller, Equations (14)–(16), can make the velocity-vector tracking error of Equation (17) be semi-globally uniformly ultimately bounded:
u V = k 1 ( V d V ) T V 0 V 0 T D ^ 1 a V 0 T s ( e 1 ) V 0 T F a ( V ) ,
W d = k 2 ( V d × V ) + 1 V V 0 × D ^ 1 + a V V 0 × s ( e 1 ) + 1 V V 0 × F a ( V ) ,
u T = ( V d × V ) D ^ 2 b s ( e 2 ) k 3 ( W W d ) + d w ˙ ,
where
e 1 = V V d ,
e 2 = W W d ,
s ( e 1 ) = e 1 e 1 + ε 1 ,
s ( e 2 ) = e 2 e 2 + ε 2 ,
a > δ 1 2 ( λ 1 0.5 ) ,
b > δ 2 2 ( λ 2 0.5 ) ,
τ d ˙ w ˙ + d w ˙ = W d , τ > 0 .
D ^ 1 and D ^ 2 are the observers of D 1 and D 2 , respectively, which are depicted by
D 1 = z 1 + λ 1 e 1 , D ^ 1 = z ^ 1 + λ 1 e 1 , z ^ ˙ 1 = λ 1 ( W × V + z ^ 1 + λ 1 e 1 + u V V 0 ) ,
D 2 = z 2 + λ 2 e 2 , D ^ 2 = z ^ 2 + λ 2 e 2 , z ^ ˙ 2 = λ 2 ( z ^ 2 + λ 2 e 2 + u T ) ,
λ 1 > 0.5 , λ 2 > 0.5 .
Proof of Theorem 1. 
Theorem 1 can be proved by the following two steps.
Step 1. Firstly, the angular velocity W and u V are viewed as control input variables to control the velocity vector V tracking the given command V d .
The error e 1 is defined as the expected sliding mode surface, and a Lyapunov candidate is selected as:
L 1 = 0.5 e 1 2 , e 1 = V V d ,
and introduce the virtual control term W d as the expectation of the angular velocity W .
Differentiating the first equation of Equation (27) yields
L ˙ 1 = ( V V d ) T ( W × V + F a ( V ) + u V V 0 + D 1 ) .
We construct the control input u V and virtual angular-velocity control term W d as
u V = u ¯ V V 0 T D ^ 1 a V 0 T s ( e 1 ) V 0 T F a ( V ) , W d = W ¯ + 1 V V 0 × D ^ 1 + a V V 0 × s ( e 1 ) + 1 V V 0 × F a ( V ) ,
where u ¯ V and W ¯ are parts of u V and W d , respectively, and they will be designed in the following part.
To eliminate the chattering caused by the sign function, a continuous saturation function is employed to replace it, which can be described as
s ( e 1 ) = e 1 e 1 + ε 1 ,
where ε 1 is a small positive constant.
For an arbitrary vector χ and a unit vector V 0 , the relationship below exists
χ = V 0 T · χ · V 0 + V 0 × χ × V 0 .
Then,
D ^ 1 = [ V 0 T D ^ 1 ] V 0 + 1 V [ V 0 × D ^ 1 ] × V , s ( e 1 ) = [ V 0 T s ( e 1 ) ] V 0 + 1 V [ V 0 × s ( e 1 ) ] × V , F a ( V ) = [ V 0 T F a ( V ) ] V 0 + 1 V [ V 0 × F a ( V ) ] × V .
Substituting Equations (29), (30) and (32) into Equation (28) leads to
L ˙ 1 = ( V V d ) T { W ¯ × V + F a ( V ) + D 1 + u ¯ V V 0 ( V 0 T D ^ 1 ) V ^ ( V 0 T F a ( V ) ) V 0 1 V [ V 0 × D ^ 1 ] × V 1 V [ V 0 × F a ( V ) ] × V a [ ( V 0 T s ( e 1 ) ) V ^ ] a 1 V [ V 0 × s ( e 1 ) ] × V } = ( V V d ) T ( W ¯ × V + u ¯ V V 0 + D 1 D ^ 1 a s ( e 1 ) ) = ( V V d ) T ( W ¯ × V + u ¯ V V 0 + D ˜ 1 a s ( e 1 ) ) .
We then select W ¯ and u ¯ V to make L ˙ 1 negative
L ˙ 1 = ( V V d ) T ( W ¯ × V + u ¯ V V 0 + D ˜ 1 a s ( e 1 ) ) = ( V V d ) T ( W ¯ × V + u ¯ V V 0 ) + e 1 T D ˜ 1 a e 1 T s ( e 1 ) V d T ( W ¯ × V ) + u ¯ V ( V V d ) T V 0 + δ 1 2 ( λ 1 0.5 ) i = 1 3 e 1 i a e 1 T s ( e 1 ) .
It can be found, under the condition
u ¯ V = k 1 ( V d V ) T V 0 , W ¯ = k 2 ( V d × V ) ,
that
L ˙ 1 k 1 ( ( V V d ) T V ^ ) 2 k 2 V d × V 2 + δ 1 2 ( λ 1 0.5 ) e 1 a e 1 T e 1 e 1 + ε 1 2 k 1 L 1 + δ 1 2 ( λ 1 0.5 ) e 1 a e 1 T e 1 e 1 + ε 1 + a ε 1 2 e 1 + ε 1 = 2 k 1 L 1 + δ 1 2 ( λ 1 0.5 ) e 1 a e 1 + a ε 1 = 2 k 1 L 1 + ( δ 1 2 ( λ 1 0.5 ) a ) e 1 + a ε 1 .
As long as
a > δ 1 2 ( λ 1 0.5 ) ,
one can obtain
L ˙ 1 2 k 1 L 1 + a ε 1 .
By solving the differential inequality (Equation (38)), we have
0 L 1 a ε 1 2 k 1 + ( L 1 ( 0 ) a ε 1 2 k 1 ) e 2 k 1 t , t 0 .
We can conclude that L 1 is bounded by a ε 1 2 k 1 , i.e., 0 L 1 a ε 1 2 k 1 , t .
In summary, the control inputs
u V = k 1 ( V d V ) T V 0 V 0 T D ^ 1 a V 0 T s ( e 1 ) V 0 T F a ( V ) , W d = k 2 ( V d × V ) + 1 V V 0 × D ^ 1 + a V V 0 × s ( e 1 ) + 1 V V 0 × F a ( V ) ,
along with the Lyapunov function L 1 , thus ensuring that V converges to an arbitrary small neighborhood around the desired V d if u V and W can achieve the control values described in Equation (40).
Step 2. Second, the angular acceleration u T is the control input to maintain the angular-velocity vector W tracking the desired virtual angular-velocity vector W d .
We define the error e 2 as the expectedu sliding-mode-surface vector, which represents the tracking error of the angular velocity
e 2 = W W d .
Then, the dynamics can be rewritten as
V ˙ = W d × V + F a ( V ) + u V V 0 e 2 × V + D 1 .
The time derivative of e 2 is
e ˙ 2 = u T d w ˙ + D 2 ,
where
τ d ˙ w ˙ + d w ˙ = W d , τ > 0 .
The control algorithm is proposed by
u T = ( V d × V ) D ^ 2 b s ( e 2 ) k 3 ( W W d ) + d w ˙ ,
where
s ( e 2 ) = e 2 e 2 + ε 2 ,
and ε 2 is a small positive value.
From the above, a Lyapunov candidate is chosen by:
L = L 1 + L 2 + L D 1 + L D 2 = 1 2 e 1 2 + 1 2 e 2 2 + 1 2 D ˜ 1 2 + 1 2 D ˜ 2 2 .
Then, the time derivative of (Equation (47)) is
L ˙ = e 1 T e ˙ 1 + e 2 T e ˙ 2 + D ˜ 1 T D ˜ ˙ 1 + D ˜ 2 T D ˜ ˙ 2 k 1 ( ( V V d ) T V ^ ) 2 k 2 V d × V 2 a δ 1 2 ( λ 1 0.5 ) i = 1 3 e 1 i + δ 2 2 ( λ 2 0.5 ) b i = 1 3 e 2 i k 3 e 2 2 ( λ 1 0.5 ) D ˜ 1 2 + 0.5 δ 1 2 ( λ 2 0.5 ) D ˜ 2 2 + 0.5 δ 2 2 + a ε 1 + b ε 2 k 1 e 1 2 k 3 e 2 2 ( λ 1 0.5 ) D ˜ 1 2 ( λ 2 0.5 ) D ˜ 2 2 + ( 0.5 δ 1 2 + 0.5 δ 2 2 + a ε 1 + b ε 2 ) 2 ξ ( 1 2 e 1 2 + 1 2 e 2 2 + 1 2 D ˜ 1 2 + 1 2 D ˜ 2 2 ) + η = 2 ξ L + η ,
where
ξ = min k 1 , k 3 , λ 1 0.5 , λ 2 0.5 ,
η = 0.5 ( δ 1 2 + δ 2 2 ) + a ε 1 + b ε 2 .
By solving the differential equation of (Equation (48)), we can obtain
0 L η 2 ξ + ( L ( 0 ) η 2 ξ ) e 2 ξ t , t 0 .
We can conclude that L is bounded by η 2 ξ , i.e., 0 L η 2 ξ , t .
Then, it can be deduced that e 1 , e 2 , D ˜ 1 , and D ˜ 2 are all semi-globally uniformly ultimately bounded:
e 1 η ξ + 2 L ( 0 ) , t 0
and
e 1 η ξ , t .
We can obtain
e 1 i η ξ + 2 L ( 0 ) , t 0 , i = 1 , 2 , 3
and
e 1 i η ξ , t , i = 1 , 2 , 3 .
From inequalities Equation (52) to Equation (55), it can be concluded that the respective tracking errors of airspeed ( V T ), AOA ( α ), and sideslip angle ( β ) can converge to arbitrary small domains of the origin. In addition, the designed controller can guarantee that the vectorial sliding manifold practically converges to the origin. □
Remark 4.
d w ˙ is the derivative vector of the virtual angular-velocity controller W d . However, W d may not be practically differentiable owing to undifferentiable noise or disturbance. To solve this problem, this study is inspired by the work of [31,32]; we adopt the low-pass filter technique in a vectorial manner.
Remark 5.
Flight control law interpretation: the different terms in the control algorithm, Equations (14)–(16) can be interpreted as follows:
u V = k 1 ( V d V ) T V 0 c o n t r o l o f V V 0 T D ^ 1 a V 0 T s ( e 1 ) c a n c e l o f D 1 V 0 T F a ( V ) c a n c e l o f F a ( V ) ,
W d = k 2 ( V d × V ) c o n t r o l o f V s d i r e c t i o n + 1 V V 0 × D ^ 1 + a V V 0 × s ( e 1 ) c a n c e l o f D 1 + 1 V V 0 × F a ( V ) c a n c e l o f F a ( V ) ,
u T = ( V d × V ) c o n t r o l o f V s d i r e c t i o n k 3 ( W W d ) c o n t r o l o f W D ^ 2 b s ( e 2 ) c a n c e l o f D 2 + d w ˙ c a n c e l o f d W d d t .
From the interpretation of Equations (56)–(58), the developed combined control scheme has a symmetry structure between u V and W d , which is concise and aesthetically appealing compared with the traditional control structures that decoupled the dynamic equations and used a collection of single variables [12,16,17,18,19]. Moreover, each term is meaningful, and this feature is significant as the control parameters can be adjusted in each term directively.
Remark 6.
The developed multivariable backstepping and sliding mode controller with nonlinear disturbance observer is deduced directly on the vector-coupled dynamics. In addition, this study deals with the cross-coupling by vector operations, including inner product algorithm and cross multiply algorithm, which takes the interaction between the state variables automatically and directly into account. Unlike similar works [12,16,17,18,19], the developed controller not only avoids solving complex matrix equations and treating inverse matrixes, but also fully realizes vector-coupled control.

4. Simulations

The effectiveness and robustness of the proposed combined multivariable backstepping sliding mode controller is verified by using the ADMIRE aircraft model [33]. The deflections of the aerodynamic control surfaces δ in the ADMIRE aircraft model include the canard wings, the right and left elevons, and the rudder. These control inputs are limited by δ c 55 , 25 , δ r e 25 , 25 , δ l e 25 , 25 , and δ r 30 , 30 , respectively. They are listed on Table 1. Furthermore, m = 9100 kg and maximum engine thrust is set with 100 kN, which means the engine thrust control input term u V [ 0 , 10.99 ] m/s 2 . The initial values of the simulation are chosen as V 0 = 80 m/s, H 0 = 3000 m, α 0 = 0 deg, β 0 = 0 deg, p 0 = q 0 = r 0 = 0 rad/s. The gains of the controller are selected as k 1 = 0.2 , k 2 = 0.005 , k 3 = 0.01 , a = 10 , b = 20 , λ 1 = 5 and λ 2 = 20 . The command signal is given by
V d = V T r e f cos α r e f cos β r e f , sin β r e f , sin α r e f cos β r e f T .
The disturbances D 1 and D 2 are selected as
D 1 = 2.1 0.1 cos 0.5 t 0.1 sin t 1.2 0.3 sin 0.5 t + 0.1 sin t 1.1 + 0.2 cos 0.5 t + 0.1 sin t m / s 2 ,
D 2 = 0.1 cos 0.5 t 0.1 sin t 0.1 0.3 sin 0.5 t 0.1 sin t 0.1 cos 2 t 0.1 sin t rad / s 2 .
To illustrate the effectiveness and robustness of the proposed vectorial controller, three cases are selected. The sampling step is 0.001 s, and the simulation results are shown in Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12. To avoid the abrupt control input change, the command passes through the filter with desired flying qualities:
x r s x c s = w n 2 s 2 + 2 ξ n w n + w n 2 ,
where w n = 3 rad / s , and ξ n = 1 .

4.1. Effectiveness of the Proposed Controller

Case 1.
To illustrate the effectiveness of the proposed controller, the commands were given as follows:
V T r e f = 80 m / s , α r e f = 10 , β r e f = 0 , 0 < t 10 s , V T r e f = 90 m / s , α r e f = 2 , β r e f = 4 , 10 < t 20 s .
The simulation results in Case 1 are shown in Figure 3, Figure 4, Figure 5 and Figure 6. From Figure 3, it can be seen that, when 0 < t 10 s, only the AOA is controlled and, for 10 < t 20 s, the total velocity V T , AOA ( α ), and sideslip angle ( β ) are all controlled simultaneously. The controlled variables can converge to a small neighborhood around the given values within 4 s in the presence of uncertainties. The tracking errors of total velocity, AOA, and sideslip angle are 7.5 × 10 4 m/s, 6.4 × 10 5 deg, and 4.6 × 10 6 deg, respectively. Figure 4 illustrates the fluctuation of the control variables in the simulation, including the force control u V , virtual control of angular velocity W d , torque vector u T , and deflection vector δ (canard wings δ c , right elevon δ r e , left elevon δ l e , and rudder δ r ). The control variables can actuate promptly in response to the external disturbances. The actuators are appropriately managed to maintain high tracking performance. Based on Figure 5 and Figure 6, the nonlinear disturbance observers give better estimations of the disturbances. The estimation errors are with the maximum magnitude of 10 3 . The simulation results in Case 1 show that the triplet states can be controlled separately and simultaneously in the form of vectors, which illustrate the effectiveness of the developed controller.

4.2. Robustness of the Proposed Nonlinear Disturbance Observer Enhanced Controller

Case 2.
To illustrate the robustness of the proposed nonlinear disturbance observer enhanced controller, the total velocity ( V T ) , AOA (α) and sideslip angle (β) are controlled under noisy disturbances. The simulation conditions are similar to Case 1, and zero-mean white Gaussian noises are added on the disturbances D 1 and D 2 in Equations (60) and (61), respectively. The simulation results in Case 2 are shown from Figure 7, Figure 8, Figure 9 and Figure 10.
Remark 7.
In fact, the lumped disturbances, D 1 and D 2 , are correlated in a realistic physics environment. However, the correlation is difficult to model and simulate. Then, inspired by the simulations of [4,17,20,27], the uncertainties, D 1 and D 2 , are selected with white Gaussian noises for illustrating the robustness and filter’s character of the developed nonlinear disturbance observer. Our future research will focus on considering the correlated D 1 and D 2 which are state-dependent and develop an adaptive-gain controller under the proposed scheme.
The performance of the developed flight controller against noisy disturbances is evaluated in Case 2. From Figure 7, the reference signals can be quickly tracked even when the noise and the model uncertainties exist simultaneously. The states of the aircraft can be steered to the intended values in a timely manner. The tracking errors of total velocity, AOA, and sideslip angle are 1.7 × 10 5 m/s, 4.4 × 10 5 deg, and 7.8 × 10 5 deg, respectively. Figure 8 illustrates the fluctuation of the control variables in the simulation. The control variables can actuate promptly in response to the noisy disturbances and the chattering phenomenon exists, but it is not serious in the control inputs. In Figure 9 and Figure 10, the performance of the nonlinear disturbance observers is illustrated. The magnitude of chartering in the observer is much smaller than that in the noise. The disturbance observers can handle noisy disturbances, which can be seen as low-pass filters to some degree, and this reflects the theoretical analysis in Remark 3. In summary, the robustness and applicability of the developed multivariable backstepping sliding mode controller is further verified through the simulation studies in Case 2.

4.3. Comparison with a Decoupled Controller

Case 3.
To verify the superiority of the proposed vector-coupled multivariable backstepping sliding mode controller (VCBSC), over our previous work, it is compared with a nonhomogeneous nonlinear disturbance observer-based dynamic surface controller (NNDODSC) presented in [20] through simulations under the same conditions. The NNDODSC is proposed based on the traditional flight control model, which was described as univariate form, i.e., state equations with AOA, sideslip angle, roll angle, roll rate, pitch rate, and yaw rate, and multiple scalar equations were considered and designed using a novel nonhomogeneous nonlinear disturbance observer and dynamic surface control theory. The simulation results in Case 3 are shown in Figure 11 and Figure 12.
As shown in Figure 11, the AOA and sideslip angle stably follow their respective reference commands, and better tracking performance is achieved with faster convergent rate by using the proposed VCBSC scheme than the traditional decoupled method of the NNDODSC. Moreover, the control inputs of the two control algorithms are shown in Figure 12. In addition, we can detect from Figure 11 and Figure 12, though the control surface deflections under the decoupled controller are larger than those of the proposed vector-coupled controller when t [0,4 s], the tracking performance of the decoupled controller is worse and the settling time is longer. Larger control deflections are required for NNDODSC than VCBSC, especially in the initial stage, which illustrates that the proposed combined vector controller, which deals with the inherent cross-coupling in the flight dynamics, can improve the efficiency of the control inputs.
Remark 8.
To testify the superiority of the proposed vector-coupled multivariable backstepping sliding mode controller, we have spent time and effort on regulating the parameters for the decoupled controller by trial and error to make the performance of the closed-loop system as good as possible under the simulation scenario, and make sure that the control surface deflections of the two simulations are nearly the same. Then, it can be guaranteed that the comparison between the two controllers is fair.
Based on the above simulation results, we can conclude that the proposed vectorial controller is effective and robust, the triplet airspeed ( V T ), AOA ( α ), and sideslip angle ( β ) can converge to small regions around the desired commands, and the design objective is achieved: a high control precision can be obtained even in the case of noisy disturbances. This also validates the theoretical analysis results.

5. Conclusions

We developed a vector-coupled control scheme for multivariable flight control based on the Lyapunov stability theory. The controller had a combination of vector sliding mode control, vector backstepping control, and nonlinear distance observers. Our key findings and the main contributions of our investigation are summarized as follows:
  • The developed control scheme allowed the conversion of attitude and airspeed control of an aircraft into a space-vector tracking problem. The results showed that the triplet airspeed ( V T ), AOA ( α ), and sideslip angle ( β ) could be controlled separately and simultaneously in the form of vectors. A key feature of this study is the active use of cross-coupling in the flight dynamics instead of decoupling and passive suppression and compensation.
  • The use of the Lyapunov stability theory enabled the development of a flight control system that is semi-globally uniformly ultimately bounded. The simulation results were comparable to the theoretical results, and it was shown that the developed controller is effective and robust.
  • The developed control scheme is concise and aesthetically more appealing compared with traditional control structures, which uses decoupled collections of single variables. The combined control scheme has a symmetry structure and each term is meaningful. The feature that the control parameters can be considered and adjusted in each term directively is significant.
Our future work will concentrate on adaptive flight control design under the proposed vector-coupled control scheme considering control input saturation.

Author Contributions

Conceptualization, Y.H. and P.L.; methodology, Y.H.; software, Y.H.; validation, J.M. and P.L.; formal analysis, J.M.; writing–original draft preparation, Y.H.; writing–review and editing, J.M. and P.L.; project administration, P.L.; funding acquisition, P.L.

Funding

This work was supported by Research Fund of State Key Laboratory of High Performance Computing (NUDT) under Grant 201613-02.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations and nomenclatures are used in this manuscript:
Error signal
Estimation signal
α Angle of attack (AOA) (rad)
β Sideslip angle (rad)
γ Roll angle (rad)
ψ Yaw angle (rad)
ρ Atmospheric density (kg/ m 3 )
φ Pitch angle (rad)
gAcceleration of gravity ( m / s 2 )
mFlight mass (kg)
p , q , r Roll, pitch and yaw rates (rad/s)
q p Dynamic pressure 0.5 ρ V T 2 ( p a )
δ Deflection vector δ c δ r e δ l e δ r T (rad)
B Distribution matrix
F External force vector (N)
G Gravity vector (N)
J Inertia matrix d i a g J x , J y , J z ( kg · m 2 )
R Aerodynamics force vector (N)
T Thrust force vector (N)
u V Thrust control term ( m / s 2 )
V 0 Unit vector of flight velocity
V Flight velocity vector v x v y v z T (m/s)
W Flight angular velocity vector p q r T (rad/s)
W d Virtual angular velocity control term p d q d r d T (rad/s)
M Control torque vector m x m y m z T ( N · m )
b Body-fixed frame coordinates O x b y b z b
i Inertial frame coordinates O x i y i z i
w Wind frame coordinates O x w y w z w
V T Magnitude of velocity (m/s)
x i , y i , z i Flight position in the inertial frame (m)
D 1 Lumped acceleration disturbances vector D 11 D 12 D 13 T ( m / s 2 )
D 2 Lumped angular acceleration disturbances vector D 21 D 22 D 23 T ( deg / s 2 )

Appendix A

The gravity vector, the aerodynamics-force vector, and the thrust vector can be calculated through
R = q S R 1 T α , β C D , C Y , C L T ,
G = m R 2 T φ , ψ , γ g ,
T = m u V V 0 ,
where C D , C Y , and C L are the dimensionless coefficients, and R 1 α , β S O 3 maps the body-fixed frame coordinates Σ b to the wind frame Σ w which is given by:
R 1 α , β = cos α c o s β sin β sin α cos β cos α sin β cos β sin α sin β sin α 0 cos α ,
and rotation matrix R 2 φ , ψ , γ S O 3 maps the body-fixed frame Σ b to inertial frame Σ i :
R 2 φ , ψ , γ = c ψ c φ c ψ s φ s γ s ψ c γ c ψ s φ c γ s ψ s γ s ψ c φ c ψ c γ + s ψ s φ s γ s ψ s φ c γ c ψ s γ s φ c φ s γ c γ c φ ,
where s x and c y stand for the sin x and cos y functions with their corresponding arguments, respectively. The Euler angles φ , ψ , and γ stand for pitch angle, yaw angle, and roll angle, respectively. Then, they can be obtained through
γ ˙ φ ˙ ψ ˙ = 1 s γ t φ c γ t φ 0 c γ s γ 0 s γ / c φ c γ / c φ W ,
where t z stands for the tan z function.

References

  1. Zhu, J.H. A Survey of Advanced Flight Control Theory and Application. In Proceedings of the Conference on Computational Engineering in Systems Applications, Beijing, China, 4–6 October 2007; pp. 655–658. [Google Scholar]
  2. Shkolnikov, A.; Shtessel, Y. Aircraft Nonminimum Phase Control in Dynamic Sliding Manifolds. J. Guid. Control Dyn. 2001, 24, 566–571. [Google Scholar] [CrossRef]
  3. Zong, Q.; Wang, F.; Tian, B. Nonlinear adaptive filter backstepping flight control for reentry vehicle with input constraint and external disturbances. J. Aerosp. Eng. 2013, 28, 889–907. [Google Scholar] [CrossRef]
  4. You, M.; Zong, Q.; Tian, B.; Zhao, X.; Zeng, F. The comprehensive design of uniform robust exact disturbance observer and fixed-time controller for reusable launch vehicles. IET Control Threory Appl. 2018, 12, 638–648. [Google Scholar] [CrossRef]
  5. Li, P.; Yu, X.; Zhang, Y. Adaptive Multivariable Integral TSMC of a Hypersonic Gliding Vehicle With Actuator Faults and Model Uncertainties. IEEE/ASME Trans. Mechatron. 2017, 22, 2723–2735. [Google Scholar] [CrossRef]
  6. Zhang, L.; Bi, S.; Yang, H. Fuzzy-PID Control Algorithm of the Helicopter Model Flight Attitude Control. In Proceedings of the 2010 Chinese Control and Decision Conference, Xuzhou, China, 26–28 May 2010; pp. 1438–1443. [Google Scholar]
  7. Amato, F.; Mattei, M.; Scala, S.; Verde, L. Robust Flight Control Design for the HIRM based on Linear Quadratic Control. Aerosp. Sci. Technol. 2000, 4, 423–438. [Google Scholar] [CrossRef]
  8. Yoshimasa, O. Design of a Flight Controller for Hypersonic Flight Experiment Vehicle. Asian J. Control 2004, 6, 353–361. [Google Scholar]
  9. Huang, Y.; Xu, K.; Han, J.; Lam, J. Flight Control Design Using Extended State Observer and Non-smooth Feedback. In Proceedings of the 40th IEEE Conference on Decision and Control, Orlando, FL, USA, 4–7 December 2001; pp. 223–228. [Google Scholar]
  10. Snell, A.; Enns, D.; Garrard, W. Nonlinear Inversion Flight Control for a Supermaneuverable Aircraft. J. Guid. Control Dyn. 1992, 15, 976–984. [Google Scholar] [CrossRef]
  11. Da Costa, R.; Chu, Q.; Mulder, J. Reentry Flight Controller Design Using Nonlinear Dynamic Inversion. J. Guid. Control Dyn. 2003, 40, 64–71. [Google Scholar] [CrossRef]
  12. Sonneveldt, L.; Chu, Q.; Mulder, J. Nonlinear Flight Control Design Using Constrained Adaptive Backstepping. J. Guid. Control Dyn. 2007, 30, 322–336. [Google Scholar] [CrossRef]
  13. Hao, A.; Fidan, B.; Wu, Q.; Wang, C.; Cao, X. Sliding Mode Differentiator Based Tracking Control of Uncertain Nonlinear Systems with Application to Hypersonic Flight. Asian J. Control 2019, 21, 143–155. [Google Scholar]
  14. Castañeda, H.; Salas-Peña, O.; Leon-Morales, J. Robust flight control for a fixed-wing unmanned aerial vehicle using adaptive super-twisting approach. J. Aerosp. Eng. 2014, 228, 2310–2322. [Google Scholar]
  15. Nagesh, I.; Edwards, C. A multivariable super-twisting sliding mode approach. Automatica 2014, 50, 984–988. [Google Scholar] [CrossRef] [Green Version]
  16. 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. J. Aerosp. Eng. 2017, 232, 344–361. [Google Scholar] [CrossRef]
  17. Chen, M.; Yu, J. Disturbance Observer-Based Adaptive Sliding Mode Control For Near-Space Vehicles. Nonlinear Dyn. 2015, 15, 1671–1682. [Google Scholar] [CrossRef]
  18. Vansoest, W.R.; Chu, Q.P.; Mulder, J.A. Combined feedback linearization and constrained model predictive control for entry flight. J. Guid. Control Dyn. 2006, 29, 427–434. [Google Scholar] [CrossRef]
  19. Li, C.Y.; Jing, W.X.; Gao, C.S. Adaptive backstepping-based flight control system using integral filters. Aerosp. Sci. Technol. 2009, 13, 105–113. [Google Scholar] [CrossRef]
  20. Ma, J.; Li, P.; Zheng, Z. Disturbance observer based dynamic surface flight control for an uncertain aircraft. J. Aerosp. Eng. 2017, 232, 729–744. [Google Scholar] [CrossRef]
  21. Yu, P.; Shtessel, Y. Continuous Higher Order Sliding Mode Control with Adaption of Air Breathing Hypersonic Missile. Int. J. Adapt. Control Signal Process. 2016, 30, 1099–1117. [Google Scholar] [CrossRef]
  22. Shtessel, Y. Multiple Time Scale Flight Control Using Re-configurable Sliding Modes. In Proceedings of the AIAA Guidance, Navigation and Control Conference and Exhibit, Hilton Head, SC, USA, 20–23 August 2007; pp. 1–10. [Google Scholar]
  23. Dong, Q.; Zong, Q.; Tian, B.; Wang, F. Adaptive-gain multivariable super-twisting sliding mode control for reentry RLV with torque perturbation. Int. J. Robust Nonlinear Control 2016, 27, 620–638. [Google Scholar] [CrossRef]
  24. Härkegård, O.; Glad, T. Vector backstepping design for flight control. In Proceedings of the AIAA Guidance, Navigation and Control Conference and Exhibit, Hilton Head, SC, USA, 20–23 August 2007; pp. 1–10. [Google Scholar]
  25. Fossena, T.; Berg, S. Nonlinear Vectorial Backstepping Design for Global Exponential Backing of Marine Vessels in the Presence of Actuator Dynamics. In Proceedings of the 36th Conference on Decision and Control, San Diego, CA, USA, 12 December 2003; pp. 4237–4242. [Google Scholar]
  26. Bates, D.; Hagström, M. Nonlinear Analysis and Synthesis Techniques for Aircraft Control; Springer: Berlin, Geramany, 2007. [Google Scholar]
  27. Tian, B.; Yin, L.; Wang, H. Finite Time Reentry Attitude Control Based on Adaptive Multivariable Disturbance Compensation. IEEE Trans. Ind. Electron. 2015, 62, 5889–5898. [Google Scholar] [CrossRef]
  28. 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. 2017, 65, 4135–4144. [Google Scholar] [CrossRef]
  29. Lie, F.; Pradipta, A.; Gebre-Egziabher, D. Synthetic Air Data System. J. Aircr. 2013, 50, 1234–1249. [Google Scholar] [CrossRef]
  30. Chen, W. Disturbance observer based control for nonlinear systems. IEEE/ASME Trans. Mechatron. 2004, 9, 706–710. [Google Scholar] [CrossRef]
  31. Swaroop, D. Dynamic Surface Control of Nonlinear Systems. In Proceedings of the American Control Conference Albuquerque, Albuquerque, NM, USA, 6 June 1997; pp. 3028–3034. [Google Scholar]
  32. Swaroop, D. Dynamic Surface Control for Nonlinear Systems. IEEE Trans. Autom. Control 2000, 45, 1893–1899. [Google Scholar] [CrossRef]
  33. Forssell, L.; Nilsson, U. ADMIRE the Aero-Data Model in a Research Environment Version 4.0, Model Description; Technical Report, FOI-R-1624-SE; FOI: Stockholm, Sweden, 2005. [Google Scholar]
Figure 1. Flight model with inertial frame Σ i , body-fixed frame Σ b , and wind frame Σ w .
Figure 1. Flight model with inertial frame Σ i , body-fixed frame Σ b , and wind frame Σ w .
Symmetry 11 01225 g001
Figure 2. Schematic of the control structure.
Figure 2. Schematic of the control structure.
Symmetry 11 01225 g002
Figure 3. Case 1: flight states: (a) airspeed V T ; (c) angle of attack α ; (e) angle of sideslip β ; and their tracking errors: (b) tracking error of V T ; (d) tracking error of α ; (f) tracking error of β .
Figure 3. Case 1: flight states: (a) airspeed V T ; (c) angle of attack α ; (e) angle of sideslip β ; and their tracking errors: (b) tracking error of V T ; (d) tracking error of α ; (f) tracking error of β .
Symmetry 11 01225 g003
Figure 4. Case 1: control variables: (a) force control u V ; (b) virtual angular velocity W d ; (c) angular acceleration control u T ; and deflection vector δ : (d) canard wings δ c ; (e) right elevon δ r e ; (f) left elevon δ l e , and (g) rudder δ r .
Figure 4. Case 1: control variables: (a) force control u V ; (b) virtual angular velocity W d ; (c) angular acceleration control u T ; and deflection vector δ : (d) canard wings δ c ; (e) right elevon δ r e ; (f) left elevon δ l e , and (g) rudder δ r .
Symmetry 11 01225 g004
Figure 5. Case 1: (a) the disturbance D 11 and its observer D ^ 11 ; (b) the estimation error D ˜ 11 ; (c) the disturbance D 12 and its observer D ^ 12 ; (d) the estimation error D ˜ 12 ; (e) the disturbance D 13 and its observer D ^ 13 ; (f) the estimation error D ˜ 13 .
Figure 5. Case 1: (a) the disturbance D 11 and its observer D ^ 11 ; (b) the estimation error D ˜ 11 ; (c) the disturbance D 12 and its observer D ^ 12 ; (d) the estimation error D ˜ 12 ; (e) the disturbance D 13 and its observer D ^ 13 ; (f) the estimation error D ˜ 13 .
Symmetry 11 01225 g005
Figure 6. Case 1: (a) the disturbance D 21 and its observer D ^ 21 ; (b) the estimation error D ˜ 21 ; (c) the disturbance D 22 and its observer D ^ 22 ; (d) the estimation error D ˜ 22 ; (e) the disturbance D 23 and its observer D ^ 23 ; (f) the estimation error D ˜ 23 .
Figure 6. Case 1: (a) the disturbance D 21 and its observer D ^ 21 ; (b) the estimation error D ˜ 21 ; (c) the disturbance D 22 and its observer D ^ 22 ; (d) the estimation error D ˜ 22 ; (e) the disturbance D 23 and its observer D ^ 23 ; (f) the estimation error D ˜ 23 .
Symmetry 11 01225 g006
Figure 7. Case 2: flight states: (a) airspeed V T ; (c) angle of attack α ; (e) angle of sideslip β ; and their tracking errors: (b) tracking error of V T ; (d) tracking error of α ; (f) tracking error of β .
Figure 7. Case 2: flight states: (a) airspeed V T ; (c) angle of attack α ; (e) angle of sideslip β ; and their tracking errors: (b) tracking error of V T ; (d) tracking error of α ; (f) tracking error of β .
Symmetry 11 01225 g007
Figure 8. Case 2: control variables: (a) force control u V ; (b) virtual angular velocity W d ; (c) angular acceleration control u T ; and deflection vector δ : (d) canard wings δ c ; (e) right elevon δ r e ; (f) left elevon δ l e , and (g) rudder δ r .
Figure 8. Case 2: control variables: (a) force control u V ; (b) virtual angular velocity W d ; (c) angular acceleration control u T ; and deflection vector δ : (d) canard wings δ c ; (e) right elevon δ r e ; (f) left elevon δ l e , and (g) rudder δ r .
Symmetry 11 01225 g008
Figure 9. Case 2: (a) the disturbance D 11 and its observer D ^ 11 ; (b) the estimation error D ˜ 11 ; (c) the disturbance D 12 and its observer D ^ 12 ; (d) the estimation error D ˜ 12 ; (e) the disturbance D 13 and its observer D ^ 13 ; (f) the estimation error D ˜ 13 .
Figure 9. Case 2: (a) the disturbance D 11 and its observer D ^ 11 ; (b) the estimation error D ˜ 11 ; (c) the disturbance D 12 and its observer D ^ 12 ; (d) the estimation error D ˜ 12 ; (e) the disturbance D 13 and its observer D ^ 13 ; (f) the estimation error D ˜ 13 .
Symmetry 11 01225 g009
Figure 10. Case 2: (a) the disturbance D 21 and its observer D ^ 21 ; (b) the estimation error D ˜ 21 ; (c) the disturbance D 22 and its observer D ^ 22 ; (d) the estimation error D ˜ 22 ; (e) the disturbance D 23 and its observer D ^ 23 ; (f) the estimation error D ˜ 23 .
Figure 10. Case 2: (a) the disturbance D 21 and its observer D ^ 21 ; (b) the estimation error D ˜ 21 ; (c) the disturbance D 22 and its observer D ^ 22 ; (d) the estimation error D ˜ 22 ; (e) the disturbance D 23 and its observer D ^ 23 ; (f) the estimation error D ˜ 23 .
Symmetry 11 01225 g010
Figure 11. Case 3: flight states: (a) angle of attack α ; (c) angle of sideslip β ; and their tracking errors: (b) tracking error of α ; (d) tracking error of β .
Figure 11. Case 3: flight states: (a) angle of attack α ; (c) angle of sideslip β ; and their tracking errors: (b) tracking error of α ; (d) tracking error of β .
Symmetry 11 01225 g011
Figure 12. Case 3: (a) canard wings δ c ; (b) right elevon δ r e ; (c) left elevon δ l e ; and (d) rudder δ r .
Figure 12. Case 3: (a) canard wings δ c ; (b) right elevon δ r e ; (c) left elevon δ l e ; and (d) rudder δ r .
Symmetry 11 01225 g012
Table 1. Control input units and maximum values.
Table 1. Control input units and maximum values.
ControlUnitMin.Max.Rate Limit
Canard wingsdeg−5525±50 deg/s
Right elevondeg−2525±50 deg/s
Left elevondeg−2525±50 deg/s
Rudderdeg−3030±50 deg/s

Share and Cite

MDPI and ACS Style

Han, Y.; Li, P.; Ma, J. Vector-Coupled Flight Controller Design Based on Multivariable Backstepping Sliding Mode. Symmetry 2019, 11, 1225. https://doi.org/10.3390/sym11101225

AMA Style

Han Y, Li P, Ma J. Vector-Coupled Flight Controller Design Based on Multivariable Backstepping Sliding Mode. Symmetry. 2019; 11(10):1225. https://doi.org/10.3390/sym11101225

Chicago/Turabian Style

Han, Yang, Peng Li, and Jianjun Ma. 2019. "Vector-Coupled Flight Controller Design Based on Multivariable Backstepping Sliding Mode" Symmetry 11, no. 10: 1225. https://doi.org/10.3390/sym11101225

APA Style

Han, Y., Li, P., & Ma, J. (2019). Vector-Coupled Flight Controller Design Based on Multivariable Backstepping Sliding Mode. Symmetry, 11(10), 1225. https://doi.org/10.3390/sym11101225

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