Next Article in Journal
On the Feasibility of a Launcher-Deployable High-Altitude Airship: Effects of Design Constraints in an Optimal Sizing Framework
Previous Article in Journal
Performance Analysis of Empennage Configurations on a Surveillance and Monitoring Mission of a VTOL-Plane UAV Using a Computational Fluid Dynamics Simulation
Previous Article in Special Issue
Studying and Optimizing the Take-Off Performance of Three-Surface Aircraft
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Backstepping Nonsingular Terminal Sliding-Mode Attitude Control of Flexible Airships with Actuator Faults

1
School of Aeronautics and Astronautics, Shanghai Jiao Tong University, Shanghai 200240, China
2
School of Aerospace, Transport and Manufacturing, Cranfield University, Cranfield MK43 0AL, UK
3
AVIC Xi’an Flight Automatic Control Research Institute, Xi’an 710065, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(4), 209; https://doi.org/10.3390/aerospace9040209
Submission received: 27 December 2021 / Revised: 27 February 2022 / Accepted: 15 March 2022 / Published: 11 April 2022
(This article belongs to the Special Issue Aircraft Modeling for Design, Simulation and Control)

Abstract

:
This paper studies the attitude tracking control of a flexible airship subjected to wind disturbances, actuator saturation and control surface faults. Efficient flexible airship models, including elastic deformation, rigid body motions, and their coupling, are established via Lagrange theory. A fast-nonsingular terminal sliding-mode (NTSM) combined with a backstepping control is proposed for the problem. The benefits of this approach are NTSM merits of high robustness, fast transient response, and finite time convergence, as well as the backstepping control in terms of globally asymptotic stability. However, the major limitation of the backstepping NTSM is that its design procedure is dependent on the prior knowledge of the bound values of the disturbance and faults. To overcome this limitation, a wind observer is designed to compensate for the effect of the wind disturbances, an anti-windup compensator is designed to compensate for actuator saturation, and an adaptive fault estimator is designed to estimate the faults of the control surfaces. Globally exponential stability of the closed-loop control system is guaranteed by using the Lyapunov stability theory. Finally, simulation results demonstrate effectiveness and advantages of the proposed control for the Skyship-500 flexible airship, even in the presence of unknown wind disturbances, control surface faults, and different stiffness variants.

1. Introduction

Airships have many potential applications for transport, environment surveillance, communication relay, and aerial photography, etc. These applications have attracted much research interest [1]. As one of the lighter-than-air unmanned vehicles, reliable autonomous control is highly desired for different airship missions. Thus, the key objective is to realize accurate trajectory-tracking control for the unmanned airship. However, due to the inherent system nonlinearity and underactuation of the airship under diverse environment disturbances, this makes the trajectory-tracking control of the airship quite challenging and difficult.
To achieve high-precision tracking performance, many control schemes of the airship have been proposed to solve the trajectory tracking problem, which includes gain scheduling control [2], backstepping design [3,4], sliding-mode control [5,6], neural network control [7,8], and fuzzy control [9]. Among these methods, the sliding-mode control and backstepping design are potentially useful approaches to overcome disturbance and faults of the airship. The backstepping approach has been proposed for the robust control of nonlinear systems in a strict feedback form through the construction of control Lyapunov functions (CLFs) [10]. The nonlinear backstepping design has been widely studied for various aerospace and underwater vehicle control problems [11,12]. Meanwhile, the sliding-mode control, due to its insensitivity to disturbances and model uncertainties, has been applied to the trajectory-tracking control for the airship with unknown disturbances and winds [5]. To make the output response quickly converge to the desired values, the terminal sliding-mode control (SMC) is proposed for robust control within a finite time [7,13]. Mofid et al. proposed a super-twisting terminal SMC and barrier function terminal SMC to deal with the quad-rotor unmanned aerial vehicles (UAVs) with delay input [14] and drive the error dynamics to converge on a region near the origin within a finite time [15]. The event-triggered fractional-order SMC technique was proposed to stabilize the quadrotor UAV with external random/time-varying disturbances [16]. In addition, recently, some researchers developed the backstepping sliding-mode control design to use combination merits of SMC and the backstepping design [17,18].
On the other hand, the system faults and disturbances of the airship are unknown and nonlinear. Thus, it is necessary for disturbance rejection and fault compensation to achieve the desired robust performance. Guo et al. [19] investigated a distributed fault-tolerant sliding-mode control for 2-D plane vehicular platoon systems, where a neuro-adaptive fault-tolerant control scheme was implemented by combining a sliding-mode control technique with a radial basis function neural network (RBFNN) to guarantee the finite-time stability. Xiao et al. [20] proposed an adaptive integral sliding-mode approach to realize fault-tolerant tracking control for a multi-vectored thrust ellipsoidal airship. Liu and Whidborne [8] proposed a neural network adaptive backstepping fault tolerant control for the unmanned airship with multi-vectored thrusters. These methods not only guarantee the system with faults of robust stability, but also tolerate the unknown actuator faults, even in the disturbance environment.
However, all the above nonlinear designs are focused on rigid-body airship models, which cannot guarantee the performance for flexible cases. To realize high altitude and long endurance missions, the stratospheric airship often uses the flexible inflated membrane material as the airship hull; thus, the weight is lighter and endurance is longer. Unlike conventional rigid airships, the mutual interaction between flexibility, aerodynamics, and flight dynamics of flexible airships should be modeled [21,22], because this interaction will improve the model precision and control performance [23]. For the trajectory-tracking control of the flexible airship, Li et al. [21] proposed an integrated model, including the flight dynamics, structural dynamics, aerostatics, and aerodynamics, but they only analyzed the linear mode characteristics of the flexible airship without consideration of nonlinear control problems. The main difficulties for the flexible airship are nonlinear and coupling control between the rigid and flexible body. Bennaceur and Azouz [24] proposed a Lagrangian method to build up an efficient model of flexible airships. In addition, Han et al. [25] proposed a command filtered backstepping approach to design the trajectory controller for a flexible airship. However, they did not consider disturbances and fault effects on the flexible airship, both of which can result in system instability. In fact, from long endurance-hovering in the air and working in diverse environments day and night, the sensors and actuators of the flexible stratospheric airship are easily impaired. Thus, fault tolerant control of the airship with actuator faults is studied here.
Motivated by Li’s [21] and Han’s work [25], the main work in this paper is a novel backstepping nonsingular terminal sliding-mode (BNTSM) control scheme proposed for attitude tracking of the airship with unknown faults, disturbances, and actuator saturation. In the proposed BNTSM control scheme, the nonsingular terminal sliding-mode technique and backstepping technique are integrated. The proposed BNTSM control has advantages of high robustness, fast transient response, and finite time convergence, and the globally asymptotic stability of the closed-loop system can be guaranteed. Furthermore, a wind observer and an adaptive fault estimator are designed to observe the wind disturbances and estimate the control surface faults; thus, the limitation of the BNTSM control dependent on prior knowledge is overcome. The main contributions of the proposed scheme are summarized as follows.
  • A backstepping nonsingular terminal sliding-mode control scheme is proposed for attitude tracking control of the flexible airship with actuator faults, actuator saturation, and uncertainties of stiffness.
  • A wind observer with an adaptive disturbance observer is designed to reject variable external bounded disturbances and cope with model parameter uncertainties.
  • An anti-windup compensator based on proportion to saturation errors is used to compensate actuator saturation.
  • An adaptive fault estimator is incorporated into the BNTSM control to implement fault estimation and fault tolerant control.
This paper is organized as follows. The nonlinear dynamics-modeling of the flexible airship and the robust control problem are described in Section 2. The BNTSM control is proposed in Section 3, and stability is proven for the corresponding closed-loop tracking system. In Section 4, three scenarios are simulated to demonstrate the BNTSM control performance. Section 5 gives the conclusions.

2. Airship Modeling and Problem Formulation

The structure of the flexible airship is shown in Figure 1. There are four control surfaces, including two rudders and two elevators on the tail fins, and two propellers are set up on each side of the hull. The gondola under the airship envelope houses the avionics system and the flight control system and other payloads.

2.1. Airship Kinematics Model

The equations of motion of a flexible airship are expressed in a body frame Fb{Obxbybzb} fixed on the undeformed airship with its origin, Ob, at the center of volume (CV), see Figure 1. The motion of the airship is described as the translation and rotation of this body frame with respect to the inertial frame FI{OIxIyIzI}, plus the deformation of the material points on the body relative to the body frame [21]. A centerline frame Fp{Opxpypzp} is introduced to calculate the aerodynamics forces of the flexible airship. The elastic displacement u of a material point on the airship can be written by a summation of the shape functions, as
u = q i ( t ) Φ i ( r ) ,   ( i = 1 ,   2 ,   N )
where q i are the time-dependent generalized coordinates, Φ i = [ Φ x i Φ y i Φ z i ] T are the structural mode shape vector, and r is the position vector of the material point from the origin, Ob. The associated velocity distribution over the elastic body is:
v = υ a + Ω × r + Ω × u + u ˙
where υ a = [u, v, w]T is the translational velocity vector in Fb, Ω = [ p q r ] T is the angular rate vector, and the × denotes cross-product operation. The first two terms of Equation (2) are related to the rigid body translation and rotation, while the last two terms reflect the influence of the elasticity on the local velocity.
Assume the airship is modeled as a free-free Euler–Bernoulli beam only undergoing bending deformation; hence, if the internal pressure is high enough to prevent wrinkling, then the hull bending can be described as the fundamental mode shape of the airship. The shape functions Φ are chosen as the natural vibration mode shapes of the beam, and 2N shape functions are employed to describe the airship’s deflection, that is,
{ Φ i = [ 0 Φ i 0 ] T ,   the   bending   in   the   Oxy   plane   Φ N + i = [ 0 0 Φ i ] T ,   the   bending   in   the   Oxz   plane ( i = 1 ,   2 ,   N )
where Φ i is the ith natural mode shape of the free-free Euler–Bernoulli beam. The kinematics equations of the position and attitude of the airship are [3]:
{ ξ ˙ = R ( η ) υ a + υ w η ˙ = J ( η ) Ω
where ξ = [x, y, z]T is the position vector, η = [ ϕ θ ψ ] T is the attitude vector, υ a = [u, v, w]T is the translational velocity vector in Fb, and υ w is the wind velocity in fixed frame. From Equation (4), it can be obtained that the local velocity is:
υ = υ a + R T ( η ) υ w  
where R(η) denotes the direction cosine matrix [3], and the transformation matrix J is:
J ( η ) = ( 1 s ( ϕ ) tan θ c ( ϕ ) tan θ 0 c ( ϕ ) s ( ϕ ) 0 s ( ϕ ) sec θ c ( ϕ ) sec θ )
where |θ| < π/2 is assumed to avoid the matrix singularity, and s(·) and c(·) denote sine and cosine functions, respectively.
The motions of the flexible airship are described as Lagrange equations [24], and the flexible dynamics model is obtained as follows [21]:
M s y s [ υ ˙ a Ω ˙ q ¨ ] = [ 0 0 S E ] + [ F I M I Q I ] + [ F G M G Q G ] + [ F A S M A S Q A S ] + [ F A D M A D Q A D ] + [ F C M C Q C ]
where M s y s = [ M 11 M 12 M 13 M 21 M 22 M 23 M 31 M 32 M 33 ] is the airship total mass, the virtual mass terms are included, FI denotes the inertia force vector, FG denotes the gravity vector, FAS denotes the aerostatic force vector, FAD denotes the aerodynamics force vector, and FC denotes the control force vector; M i and Q i denote the associated moments and general elastic forces of the flexible airship, i = I, G, AS, AD and C. The local velocity is as Equation (5); q = [ q 1 , q 2 , , q 2 N ] T denotes the generalized elastic coordinate vector. S E denotes the internal elastic force vector,
S E = K q
where the stiffness matrix K meets
K i , j = K N + i , N + j = L E I Φ i Φ j d x , ( i , j = 1 , 2 , N )  
where Φ i = d 2 Φ i d x 2 . EI denotes bending stiffness with EI = πR3 ET0, E is the elastic modulus of the hull envelope, T0 is its thickness, and R is the hull radius. Equation (7) can be rewritten as
X ˙ ( t ) = f ( X ) + g ( X ) u ( t )
where X = [ υ a T Ω T q ˙ T ] T is the augmented state vector, u = [ δ T δ e δ r ] T denote control inputs of the thruster, elevator, and rudder, and
f = 1 M s y s ( [ 0 0 S E ] + [ F I M I Q I ] + [ F G M G Q G ] + [ F A S M A S Q A S ] + [ F A D M A D Q A D ] )
g = [ g 1 g 2 g 3 ] 9 × 3
g 1 = 1 M s y s [ F ¯ k T ( S T + i = 1 2 N q i S Φ i ) F ¯ k T F ¯ k T T [ Φ x i Φ y i Φ z i ] T ] 9 × 1 , g 2 = 1 M s y s [ F ¯ δ e ( S f + i = 1 2 N q i S Φ i ) F ¯ δ e F ¯ δ e T [ Φ x i Φ y i Φ z i ] T ] 9 × 1 ,   g 3 = 1 M s y s [ F ¯ δ r ( S f + i = 1 2 N q i S Φ i ) F ¯ δ r F ¯ δ r T [ Φ x i Φ y i Φ z i ] T ] 9 × 1 ,
S T = [ 0 z T y T z T 0 x T y T x T 0 ] , S Φ i = [ 0 Φ z i Φ y i Φ z i 0 Φ x i Φ y i Φ x i 0 ] , S f = [ 0 z r y r z r 0 x r x m y r x r + x m 0 ] ,
F ¯ k T = [ k T 0 0 ] T , F ¯ δ e = q ¯ S r e f [ 0 0 C z δ e ] T , F ¯ δ r = q ¯ S r e f [ 0 C y δ r 0 ] T ,
where (xT, yT, zT) and (xr, yr, zr) denote positions of the thrust and control surface in the body frame, Fb, respectively; xm is the position of the aerostatic center along the x-axis, kT is control gain for thruster. C z δ e and C y δ r denote the aerodynamic derivative coefficients of the elevator and rudder, respectively, q ¯ = 1 2 ρ V 2 denotes the dynamic pressure, V is the airspeed, and ρ is the air density.
Next, we detailed and present each force and moment on the right-hand side of Equation (11). The inertia force and moment vector of the flexible airship in Equation (11) is [21]:
[ F I M I Q I ] = [ m Ω × Ω × c t o t a l × 2 P q ˙ c t o t a l × Ω × Ω × J t o t a l 2 H q ˙ , t o t a l P q ˙ T H q ˙ , t o t a l T 2 M e q ˙ ] [ υ a Ω q ˙ ]
where
c t o t a l × = M 21 = ( m ( r + u ) d m ) × = ( m r G + m ( q i ( t ) Φ i ( r ) ) d m ) × , J t o t a l = M 22 = m ( r + u ) × ( r + u ) × d m , P q ˙ = Ω × P = Ω × M 13 , H q ˙ = [ J r u , 1 Ω J r u , 2 Ω J r u , 2 N Ω ] , J r u , i = m r × Φ i × d m , H q ˙ , t o t a l = H q ˙ + H q ˙ u , H q ˙ u = [ q j ( t ) J u u , 1 j Ω q j ( t ) J u u , 2 j Ω q j ( t ) J u u , 2 N j Ω ] ( i ,   j = 1 ,   2 ,   ,   2 N ) ,   M e q ˙ = m Φ T Ω × Φ d m .
The superscript denotes the skew-symmetric matrix form of a vector (corresponding to a cross-product operation). Consider the gravity, aerostatic force, and moment vector of the flexible airship in Equation (11). Gravity, the resulted moment, and its generalized one are described as follows:
F G = m g g ^
M G = ( r G + u G ) × F G = r G × F G + 1 m i = 1 2 N q i ( t ) p i × F G
Q G i = m g g ^ T Φ i d m = g g ^ T p i , i = 1 , 2 , , 2 N
where g denotes that the acceleration of gravity and unit direction vector of gravity is:
g ^ = [ sin θ cos θ sin ϕ cos θ cos ϕ ] T
u G = m u d m / m = 1 m i = 1 2 N q i ( t ) p i , p i = m Φ i d m
The aerostatic force in the body frame in Equation (11) is:
F A S = ρ g V B g ^
The resulted moment and its generalized moment are presented as follows:
M A S = u V × F A S = ( V B u d V / V B ) × F A S = q i ( t ) I A S , i × F A S
Q A S , i = L ρ g ^ T Φ i S d x = F A S T I A S , i , i = 1 , 2 , , 2 N
where VB denotes the volume of the airship hull and S is the hull cross-sectional area:
u V = 1 V B V B u d V , I A S , i × = 1 V B V B Φ i d V = 1 V B L Φ i S ( x ) d x
The aerodynamics force and moment vector of the flexible airship in Equation (11) is:
[ F A D M A D Q A D ] = [ F p M p Q p ] + [ F v M v Q v ] + [ F F M F Q F ] + [ F H ( F ) M H ( F ) Q H ( F ) ] + [ F axial M axial Q axial ]
where subscripts p, v, F, H(F), and the axial on the right-hand side of Equation (27) denote the forces and moments due to the potential-flow (related to the added mass), viscous effect on the hull, aerodynamic force acting on fins and acting on the hull due to the fins, and axial drag, respectively. The calculation of the aerodynamics forces is the same as in [21].
The control force and moment vector F C , M C , Q C of the flexible airship includes those of propellers and control surface deflections. The thrust in the body frame is F T = [ F T x F T y F T z ] T ; the corresponding moment and generalized moment are as follows:
M T = r T × F T + q i ( t ) Φ i ( r T ) × F T
Q T i = F T Φ i ( r T )
where Φ i ( r T ) is the mode shape function at the propeller-mounted position r T . Forces due to the control surface deflection are presented in the body frame, as
F δ T = [ Δ C D Δ C L ( δ r ) Δ C L ( δ e ) ] q ¯ S f
where S f denotes the area of control surface, and Δ C L and Δ C D are the lift and drag resulted from the control surface deflection. The corresponding moment and generalized moment are:
M δ = r δ × F δ + q i ( t ) Φ i ( r δ ) × F δ
Q δ , i = F δ T Φ i ( r δ )
where Φ i ( r δ ) is the shape function at the control surface-mounted position of r δ .
Considering the model uncertainties or disturbances often occurring in the dynamics system of Equation (10), assume the atmospheric parameters are stable and disturbances are bounded; thus, a model uncertainty, d, and the actuator saturation are introduced, thus the system (10) can be modified, as
X ˙ ( t ) = f ( X ) + g ( X ) sat ( u a ( t ) ) + Γ d ( t )
where d denotes disturbances except for winds and model uncertainties, such as aerodynamic coefficients and structural stiffness; sat ( ) denotes saturation function. As the airship has a large inertia and its motion is slow and sedate, d can be assumed to be a slow-varying disturbance, which can be estimated on-line by using an adaptive law. Let d ^ be the estimated value of the uncertain parameter d; d ˜ = d ^ d is the associated estimated error. Finally, Г denotes the disturbance coefficient matrix. The faults can be modeled as abrupt changes of the nominal control action from u a [26,27],
u a f ( t ) = u a ( t ) + ( I ρ ¯ a ( t ) ) ( f a ( t ) u a ( t ) )
that is,
u a f ( t ) = u h + u f = ρ ¯ a u a ( t ) u h + ( I ρ ¯ a ( t ) ) f a ( t ) u f
where u a = [ u a , 1 , u a , 2 , , u a , m ] R m is the actual control output; the diagonal matrix ρ ¯ a = diag { ρ a , 1 , ρ a , 2 , , ρ a , m } denotes the operational effectiveness of the actuators; I is the m×m identity matrix; the fault vector f a = [ f a , 1 , f a , 2 , , f a , m ] R m denotes the control action from the failed or un-manipulated actuators; u h = ρ ¯ a u a ( t ) denotes the remaining control of the health actuators; u f = ( I ρ ¯ a ( t ) ) f a ( t ) denotes the fault of the inputs.
Assumption 1.
In Equation (34), for the multiplicative and additive actuator faults, ρ a , i and f a , i are all bounded. In addition, ρ ˙ a , i and f ˙ a , i exist and are bounded, i = 1, 2, …, m.
Substituting Equation (35) into Equation (33) yields:
X ˙ ( t ) = f ( X ) + g ( X ) ρ ¯ a sat ( u a ( t ) ) + g ( X ) ( I ρ ¯ a ( t ) ) f a ( t ) + Γ d
Remark 1.
There are two underactuated cases for the flexible airship. Case 1. The airship without a lateral tilt angle of the propellers (i.e., the thrust direction is fixed) is underactuated in the y-direction (that is, the lateral control force, Tsy = 0); thus, the sway velocity, v, cannot be directly controlled. If the wind is in the presence in this case, then the airship can align against the wind through the yaw motion, reducing the lateral forces requirement to a low and acceptable value; thus, the lateral force input can vanish in stationary conditions. Case 2. The airship works in Case 1 without ailerons or differential actuators (i.e., δeL = δeR, δrU = δrB, and the roll control moment, MTx ≈ 0). Sway velocity, v, and bank angle, φ, cannot be directly controlled. In this case, the disturbance of the roll moment resulting from the wind can be attenuated by the airship roll damp; thus, the roll moment input can vanish in stationary conditions.

2.2. Fault Tolerant Trajectory Tracking Problem

The fault-tolerant trajectory tracking problem is a challenge for the flexible stratospheric airship because it is a hard task to achieve precise trajectory tracking under the actuator faults; the main reason is that the interaction between fluid and structure should be considered. Presently, consider the system of Equations (10) and (36)—the control task is to design a fault-tolerant trajectory-tracking controller, such that the closed-loop system is globally asymptotically stable and the output trajectory, η, is steered towards a given reference trajectory, ηr, with lim t e ( t ) 2 < ε 0 , even in a specified model parameter uncertainty, unknown wind disturbances, and control surface faults; where the tracking error is e ( t ) = η ( t ) η r ( t ) , ε0 is a prescribed constant, and 2-norm is e ( t ) 2 = e T e .

3. BNTSM-Based Trajectory Tracking Design

In this section, the robust BNTSM controller is designed to solve the problem in Section 2.2. Firstly, a command filter is used to meet the magnitude and rate constraints of the input signal. The command filter can be obtained by using a second-order lower pass filter, as follows:
G r ( s ) = X r ( s ) X r 0 ( s ) = ω n 2 s 2 + 2 ζ n ω n s + ω n 2
where ζ n denotes the filter-damping ratio and ωn denotes the undamped nature frequency of the filter. To reduce the filter error e = x r 0 x r , the filter frequency, ωn, is the bandwidth of X r 0 ( s ) , which is generally less than that of Gr(s), and then ωn can be selected.
Secondly, for the system (4), the wind field is assumed as a constant, due to the stable meteorological condition in the stratosphere [28]. Some people regarded wind speed as one kind of external disturbance, but it is hard to separate wind and model uncertainty and they are not easily handled. Here, a wind observer is introduced to improve tracking performances. A wind observer state is defined, [ ξ ^ T υ ^ w T ] T , where ξ ^ and υ ^ w denote estimates of ξ and υ w . The wind can be sensed by the wind speed sensors in the fixed frame, and the observer dynamics are designed, as follows [4]:
[ ξ ^ ˙ υ ^ ˙ w ] = [ R ( η ) υ a 0 ] + [ L ξ I 3 L w 0 ] [ ξ ξ ^ υ ^ w ]
The estimation error of the position and wind speed are defined as e 4 = [ ξ ξ ^ υ w υ ^ w ] , and the estimation error e4 can be obtained as
e ˙ 4 = [ ξ ξ ^ υ w υ ^ w ] = [ L ξ I 3 L w 0 ] e 4 = A ˜ e e 4
where L ξ , L w are gain matrices such that A ˜ e = [ L ξ I 3 L w 0 ] be Hurwitz; I3 is the unit matrix with a dimension of 3 × 3. Thus, there exists a positive definite symmetrical matrix, P e , such that
d d t ( e 4 T P e e 4 ) = e 4 T Q e e 4
where Q e = diag ( Q ξ , Q υ w ) , Q ξ , and Q υ w are symmetric positive definite matrices and the weight matrix, P e , meets the following Lyapunov equation,
A ˜ e T P e + P e A ˜ e = Q e
Thirdly, consider the kinematics in Equation (4) and dynamics in Equation (36). The control objective is to track a reference signal, x 1 r = η r , with the derivative x ˙ 1 r under unknown disturbances and faults. The tracking error vectors of the attitudes are defined as
z 1 = x 1 x 1 r
z 2 = x 2 x 2 r
where x 1 = η x 2 = x ˙ 1 = η ˙ , x 1 r is the reference or desired attitude angle, and x 2 r is the output of a command filter.
Presently, the BNTSM controller is designed, and the kinematics Equation (4) and dynamics Equation (33) are considered. To obtain a fast and transient response and finite time convergence without a singular problem, a non-singular terminal sliding-mode (NTSM) surface or manifold is designed by using the fractional order derivative, as follows [29,30:
s = z ˜ 1 + α γ z ˜ 1 λ + 1 β z ˜ 2 p / q = z ˜ 1 + α γ | z ˜ 1 | λ sgn ( z ˜ 1 ) + 1 β | z ˜ 2 | p / q sgn ( z ˜ 2 )
where z ˜ 1 and z ˜ 2 are defined in (A6), the SMC parameter λ > p/q, α γ > 0 , s = [ s 1 s 2 s 3 ] T R 3 denotes the sliding-mode surface, i.e., s1 = s(φ), s2 = s(θ), s3 = s(ψ) and 1 < p/q < 2, β > 0; sgn(·) is the sign function. In addition, the derivative of the sliding-mode surface, s, is
s ˙ = z ˜ ˙ 1 + α γ λ | z ˜ 1 | λ 1 z ˜ ˙ 1 + 1 β p q | z ˜ 2 | ( p / q 1 ) z ˜ ˙ 2
Using an extended CLF, a control law will be designed to drive the virtual angular rate error, z ˜ 2 , to zero, while ensuring that the current attitude vector, η, will converge to the desired value of ηr. The attitude controller is derived in two steps, and the detailed design is presented in Appendix A.
Presently, we consider the input saturation problem, in order to analyze the effect of the actuator saturation on the closed-loop dynamics; the actuator output is
sa t ( u a ) = { u a f 0 + u ^ f   for   | u a f 0 + u ^ f   | < U lim U lim   for   u a f 0 + u ^ f > U lim U lim   for   u a f 0 + u ^ f < U lim
where u a f 0 is the input of the BNTSM controller as Equation (A13), and u ^ f is the fault estimator as Equation (A19).
Furthermore, defining ∆u as the difference between the desired control input u and the actuator output, i.e., ∆u = ua − sat(ua), an anti-windup compensator is designed to compensate for the actuator saturation, as follows:
v = K s ( u a sat ( u a ) )
where K s is the control gain of the saturation compensator, and then the total control is
u a = u a f 0 + u ^ f + v
Remark 2.
For the fault estimator (A19) (in Appendix A), if z ˜ 2 converges to zero, then u ^ ˙ f will converge to zero, and the fault becomes constant; a disturbance observer method [31] can be used to estimate u ^ f = c o n s t .
Note that the operational effectiveness, ρ ¯ a , (or the multiplicative fault) of the actuators in Equation (34) can also be estimated by using a similar method as Equation (A19) [32]; thus, there is no more detail.
The BNTSM controller for the trajectory-tracking control is designed as Figure 2, which includes the flexible airship dynamics and a BNTSM control module. The BNTSM controller is designed based on the above steps. The disturbance estimator and the fault estimator are governed by Equations (A18) and (A19). A wind observer is also used according to Equation (38). The control output is filtered by the command filter (37).

4. Simulation and Analysis

The considered model is the flexible Skyship-500 airship, and the structure parameters are listed in Table 1; the aerodynamic coefficients are as in [21].
The initial position in Fg of the airship is ξ 0 = [0, 0, 0 m]T, the initial velocity in Fb is υ0 = [8.42 m/s, 0, 0]T, the initial attitude is η0 = [0, 0, 0]T, and the initial angular velocity is ω0 = [0, 0, 0]T. The position constraints for the elevator and rudder are [−30°, 30°]. Three scenarios of bounded wind disturbances, control surface faults, and the variable stiffness of the flexible airship envelope are simulated to illustrate the BNTSM controller performances. The software of MATLAB and Simulink are employed to solve the problems.
Scenario 1.
Attitude tracking control under unknown winds.
This scenario will demonstrate the attitude-tracking performance by using the proposed BNTSM controller and the wind observer. The wind vector is initially set as
υ w = { [ 7 , 0 , 0 ] T ( m / s ) , t 20 s [ 0 , 5 , 0 ] T ( m / s ) ,   20 s < t 50 s [ 0 , 0 , 3 ] T ( m / s ) ,   50 s < t 80 s [ 3 , 1 , 2 ] T ( m / s ) , t > 80 s
Wind parameters in (38) and (39) are chosen as L ξ = d i a g ( 1 , 1 , 1 ) , L w = diag ( 2.5 , 2.5 , 2.5 ) , and Q ξ = I 3 , Q υ w = I 3 . To illustrate the proposed performances, a backstepping integral sliding-mode control (BISMC) [3,4] and a PID controller are employed for comparison. The sliding surface s of the BISMC is defined as follows: s = λ 1 z 1 + λ 2 0 t z 1 d ( t ) + z 2 . In addition, the PID controller structure is as in Figure 3.
The controller parameters are designed as shown in Table 2.
The parameters are selected to satisfy the requirements in Section 3 after several design iterations. The damping ratio and undamped nature frequency of the command filter are selected as ζn = 0.9, ωn = 20 rad/s. The mode number, N, selects N = 2, Pe = 0.1, Γ = 1 . The gain of the actuator saturation compensator is Ks = 1, and control gain for thruster kT = 90,000. A doublet command is predefined as the desired attitude.
For an easier comparison, the averaged tracking error is defined as
E = 1 N k = 1 N ( e ( k ) 2 )
where N is the number of simulation steps, and the elapsed time (ET) is the time that MATLAB has used to complete the 10-s simulation time. This may not reflect the true computational burden of the controllers, but it can be used to provide a general idea about the comparison in computational time among the controllers. The simulation results are shown, as follows.
It can be observed from Figure 4 and Table 3 that the reference signals of pitch and yaw angles are precisely tracked, but the averaged tracking errors of the BNTSM controller are smallest, and then those by the BISMC design are smaller, and the PID controller provides worse tracking performances comparing with the BNTSM control and the BISMC. In view of the computational burden, the BNTSM controller is time consuming, while the cost of the PID controller is the cheapest. The roll motion response is free oscillating and has not been controlled due to the absence of the direct roll moment under the underactuated Case 2.
It also can be observed from Figure 5 and Table 3 and Table 4 that the responses have smaller overshoots and shorter transient times by using the BNTSM controller with wind observer compensation versus those without wind compensation. The wind speeds are precisely estimated by the proposed wind observer, as shown in Figure 6. Figure 7 demonstrates the structure mode shapes in the Oxy and Oxz planes; that is, according to (3), the 1st and 2nd bending mode in the Oxy plane are Φ 1 = [ 0 ϕ 1 0 ] T and Φ 2 = [ 0 ϕ 2 0 ] T ; the 1st and 2nd bending mode in the Oxz plane are Φ 3 = [ 0 0 ϕ 1 ] T and Φ 4 = [ 0 0 ϕ 2 ] T ; the associated nature frequency of the 1st and 2nd bending modes can be obtained by the modal rigidness and mass, i.e., ω n 1 = 22.8 rad / s and ω n 2 = 64.2 rad / s , respectively, and it can be observed from Figure 7 that the rigidness is bigger and the elastic displacement is smaller in the middle airship than those on the two sides of the airship.
Scenario 2.
Trajectory-tracking control under actuator faults with unknown winds.
This scenario will demonstrate fault-tolerant tracking performance by using the proposed BNTSM controller and the fault observer. The faults of actuator bias forces with wind are studied here. The faults are set as follows: Fault 1, the bias from the nominal value is 5 deg and the loss of effectiveness is 25% for the elevator in t > 30 s; Fault 2, the bias from the nominal value is 5 deg and the loss of effectiveness is 50% for the rudder in t > 60 s. The faults are simulated by setting bias values of each actuator at the triggering time. The wind disturbances are set as Equation (49) in Scenario 1. A doublet attitude is predefined as the desired attitudes to verify the tracking performance of the BNTSM control. The controller parameters are the same as Table 2, except for
{ c 2 = diag ( 0.2 , 20 , 0.2 ) , β = diag ( 0.1 ,   0.002 ,   0.2 ) ,   for   BNTSM c 1 = c 2 = diag ( 2 , 2 , 2 ) ,                                  for   BISMC ,  
In addition, the fault observer parameters are designed as
γ f   = 3 ,   P f = 3,000,000 × diag ( 0 , 1 , 1 ) .
By using the above design, the simulation results are demonstrated as follows.
Figure 8 and Figure 9 shows the pitch step input of 0.2, the rad shows at 5 s, and the tracking response converges the desired value at 10 s (within 10 s). Because the doublet command is a changing signal, the steady value of the tracking response also changes; thus, the overall tracking process is dynamically varying, but the response will quickly converge when the command is fixed. It can be observed from Figure 8 and Table 5 that the reference inputs of pitch and yaw angles are well tracked, even with control surface faults and unknown winds. Comparing with the PID controller and the BISMC, the BNTSM controller has smaller averaged tracking errors. In the view of elapsed time, the PID controller is the least, with the BISMC next, while the BNTSM control will cost more time. The roll motion is not controlled due to the underactuated character of the flexible airship. It also can be observed from Figure 9 and Table 5 and Table 6 that the responses have shorter transient times by using the fault observer compensation and the BNTSM controller than those without fault compensation. This demonstrates the fault tolerant control capability of the BNTSM controller. The control inputs are shown in Figure 10 and Figure 11. Figure 11 shows that the control surface faults are estimated by using the proposed estimator of (A19) (in Appendix A). There are some differences initially because the corresponding tracking errors z ˜ 2 are affected by the variable reference signals. The fault estimations are gradually close to the true values. Figure 12 shows the response of the sliding-mode variable, s, which asymptotically approach zero except for s(φ), due to the underactuated roll motion.
Scenario 3.
Trajectory-tracking control under variable stiffness and unknown faults and winds.
In this scenario, we investigate the significance envelope flexibility has on the effect of rigid-body dynamics of the airship. What effect does the hull stiffness of the airship have on the structure modes? Four hull stiffness variants, plus the baseline stiffness case, are compared. The bending stiffness of the envelope is scaled for the model variants by ±30% and ±50%, respectively. The model was trimmed in steady level flight at 9 m/s and sea level. Winds and faults are also considered as Scenario 2. The controller parameters are the same as Scenario 2 except for
{ c 2 = diag ( 0.2 ,   0.2 ,   0.4 ) , β = diag ( 0.1 ,   0.1 ,   0.2 ) ,   for   BNTSM c 1 = c 2 = diag ( 2 ,   2 ,   2 ) ,                                                 for   BISMC ,
By using the BNTSM control design, the simulation results are shown in Figure 13, Figure 14, Figure 15 and Figure 16.
From Figure 13, it can be observed that the variable stiffness of the airship envelope has a small effect on the rigid-body attitude and the position motions. However, the stiffness variation has a great effect on the structure modes, as shown in Figure 14 and Figure 15, where it can be found that the bending mode responses of qi will become a smaller and smaller amplitude of oscillation with the increasing stiffness. The general coordinate, q1, for the first bending mode, is increased in amplitude by 30% for the −50% stiffness variant, and the general coordinate, q1, reduced in amplitude by 8.3% for the +30% stiffness case. The same case happens for the generalized coordinate velocity responses. This demonstrates that stiffness can suppress structural mode oscillation. Figure 16 shows that the control inputs change little when the stiffness of the flexible airship is variable, where * denotes multiplication sign.

5. Conclusions

This paper proposes a nonlinear attitude tracking control approach for the flexible airship. A nonlinear backstepping nonsingular terminal sliding-mode control method is used to design the controller to track desired attitudes and stabilize structural modes of the flexible airship. Meanwhile, a wind observer is designed to estimate variable wind speeds. An adaptive fault estimator is designed to deal with nominal offset and loss of effectiveness faults of the control surfaces. Saturation compensator is designed to reduce effect of the actuator saturation. The stability analysis demonstrates that the closed-loop attitude tracking error dynamics are globally exponentially stable. Simulation results demonstrate that the BNTSM control designed for three scenarios of the flexible airship can achieve better attitude tracking performances by using the corresponding compensation controller, even though the airship is affected by unknown winds, control surface faults, and variable stiffness of the airship envelope. Therefore, the effectiveness and availability of the BNTSM control design are verified.

6. Future Recommendation

The future work is to expand the control strategy into a position control loop and velocity control loop for the flexible airship.

Author Contributions

Conceptualization, S.L. and J.F.W.; methodology, S.L., J.F.W., S.S. and W.L.; software, S.L. and S.S.; validation, S.L., J.F.W. and S.S.; formal analysis, J.F.W.; investigation, S.L. and W.L.; writing—original draft preparation, S.L.; writing—review and editing, S.L. and J.F.W.; visualization, S.L.; supervision, J.F.W. All authors have read and agreed to the published version of the manuscript.

Funding

This search is supported by the National Natural Science Fund under the grant number 10577012.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Sample data available on request.

Acknowledgments

This work was supported and developed by the Shanghai Jiao Tong University and Cranfield University.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

BISMCBackstepping integral sliding mode control
BNTSMBackstepping nonsingular terminal sliding mode
CLFsControl Lyapunov functions
CVCenter of volume
EIBending stiffness, E means the elastic modulus, I means the area moment of inertia
NTSMNonsingular terminal sliding mode
SMCSliding-mode control
RBFNNRadial basis function neural network

Appendix A

The attitude controller designed by using BNTSM control.
Step 1 Backstepping for the variant of z1.
Using the standard backstepping procedure, select the positive definite function, V1, as
V 1 = 1 2 z 1 T z 1
To make the derivative function V 1 / z 1 a negative definite, a virtual control is defined as
α 1 = G 1 1 ( x ˙ 1 r c 1 z 1 ) ,
where G 1 = J , J is in Equation (4). Then,
V 1 z 1 = z 1 T ( x ˙ 1 x ˙ 1 r ) = z 1 T ( G 1 α 1 x ˙ 1 r ) = z 1 T c 1 z 1 < 0
However, instead of directly applying this to virtual control α1, a new signal, x 2 r 0 , is defined as
x 2 r 0 = α 1 ζ 2 ,
where ζ2 will be designed in Step 2. The command signal, x 2 r 0 , is filtered to produce the reference signal, x 2 r , and its derivative, x ˙ 2 r . It can be implemented to enforce magnitude and rate limits through the command filter (37). By the design of the command filter (37), the signal ( x 2 , r x 2 , r 0 ) is bounded and small.
The effect of filtering on the tracking error z1 is estimated by the following stable linear filter
ζ ˙ 1 = c 1 ζ 1 + G 1 ( x 2 , r x 2 , r 0 )
where c1 > 0. To remove the effect of filtering the stabilizing functions from the tracking error, the compensated tracking error is defined as
z ˜ i = z i ζ i ,   ( i = 1 , 2 )
Re-selecting the first Lyapunov function V1 as a quadratic function of the compensated tracking error is
V 1 = 1 2 z ˜ 1 T z ˜ 1 ,
whose derivative is
V ˙ 1 = z ˜ 1 T z ˜ ˙ 1 = z ˜ 1 T ( x ˙ 1 x ˙ 1 r ζ ˙ 1 ) = z ˜ 1 T ( G 1 x 2 x ˙ 1 r ( c 1 ζ 1 + G 1 ( x 2 r x 2 r 0 ) ) ) ,
Substituting Equations (A2) and (A4) into Equation (A8) yields
V ˙ 1 = z ˜ 1 T ( G 1 ( z 2 + α 1 ζ 2 ) ) + z ˜ 1 T ( c 1 ζ 1 x ˙ 1 r ) = z ˜ 1 T c 1 z ˜ 1 + z ˜ 1 T G 1 z ˜ 2
Step 2 The effect of filtering the command signal on the tracking error, z2, is estimated by the following stable linear filter:
ζ ˙ 2 = c 2 ζ 2 + G 2 ( u a u a 0 ) ,
where G 2 = diag ( 0 , C m y δ e I y + m 55 , C m z δ r I z + m 66 ) q ¯ S L , C m y δ e , and C m z δ r denote the aerodynamic derivative coefficients of elevator and rudder, respectively, and L is the length of the airship. Note that there is no aileron input for the flexible airship, thus C m x δ a = 0 , that is, the first roll moment term, is C m x δ a I x + m 44 = 0 .
Now let us construct the second augmented quadratic CLF by the compensated tracking error, fault estimation error, and parameter disturbance estimation error, as
V 2 = V 1 + 1 2 s T s + 1 2 d ˜ T d ˜ + 1 2 e 4 T P e e 4 + 1 2 u ˜ f T P f 1 u ˜ f ,
where u ˜ f = u ^ f u f is the fault estimation error, u ^ f denotes the estimation of the fault input u f . In addition, Pe and Pf are positive definite weight matrices. Considering Equation (A9) and substituting Equations (35), (45), (A4), (A6) and (A10) into the derivative of V2, yields
V ˙ 2 = z ˜ 1 T z ˜ ˙ 1 + s T s ˙ + d ˜ T d ˜ ˙ e 4 T Q e e 4 + u ˜ f T P f 1 u ˜ ˙ f = z ˜ 1 T c 1 z ˜ 1 + z ˜ 1 T G 1 z ˜ 2 + s T ( ( c 1 z ˜ 1 + G 1 z ˜ 2 ) ( 1 + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) + 1 β p q | z ˜ 2 | ( p / q 1 ) ( f 2 + G 2 u a f + Γ 2 d x ˙ 2 r + c 2 ζ 2 G 2 ( u a f u a f 0 ) ) ) + d ˜ T d ^ ˙ ξ ˜ T Q ξ ξ ˜ υ ˜ w T Q υ w υ ˜ w + u ˜ f T P f 1 u ˜ ˙ f = z ˜ 1 T ( c 1 z ˜ 1 + G 1 z ˜ 2 ) + ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) ( s T G 1 z ˜ 2 c 1 z ˜ 1 T s ) + s T 1 β p q | z ˜ 2 | ( p / q 1 ) ( f 2 + Γ 2 d x ˙ 2 r + c 2 ζ 2 + G 2 u a f 0 ) + d ˜ T d ^ ˙ ξ ˜ T Q ξ ξ ˜ υ ˜ w T Q υ w υ ˜ w + u ˜ f T P f 1 u ˜ ˙ f ,
where ξ ˜ = ξ ξ ^ , υ ˜ w = υ w υ ^ w , u a f 0 = u h 0 + u f 0 . According to Equation (11), f 2 = 1 M 22 ( M I + M G + M A S + M A D ) ; Γ 2 denotes elements of the matrix Γ from the 4th row to the 6th one, and
M 22 = [ m 44 + I x 0 0 0 m 55 + I y 0 0 0 m 66 + I z ]
As d is an unknown constant or a slow-varying disturbance, then d ˜ ˙ = d ^ ˙ d ˙ d ^ ˙ . If d is a fast-varying process, a RBF-neural network-based backstepping control can be used [33]. The stabilizing function is defined as
u a f 0 = 1 G 2 ( β q p ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) G 1 | z ˜ 2 | ( 2 p / q ) f 2 G 2 u ^ f c 2 z ˜ 2 + x ˙ 2 r Γ 2 d ^ ( h s + ς sgn ( s ) ) )
where c2 > 0, h, ς are sliding-mode surface parameters with h > 0, ς > 0. The control input of uh is generated by using the above Equation (A13), and then
V ˙ 2 = z ˜ 1 T c 1 z ˜ 1 + z ˜ 1 T G 1 z ˜ 2 ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 z ˜ 1 T s + s T 1 β p q | z ˜ 2 | ( p / q 1 ) ( Γ 2 ( d d ^ ) + G 2 ( u f 0 u ^ f ) c 2 z ˜ 2 ) h s 2 ς | s | + d ˜ T d ^ ˙ ξ ˜ T Q ξ ξ ˜ υ ˜ w T Q υ w υ ˜ w + u ˜ f T P f 1 u ˜ ˙ f = z ˜ 1 T c 1 z ˜ 1 + z ˜ 1 T G 1 z ˜ 2 ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 z ˜ 1 T s s T 1 β p q | z ˜ 2 | ( p / q 1 ) c 2 z ˜ 2 + s T 1 β p q | z ˜ 2 | ( p / q 1 ) G 2 u ˜ f s T 1 β p q | z ˜ 2 | ( p / q 1 ) Γ 2 d ˜ + d ˜ T d ^ ˙ h s 2 ς | s | ξ ˜ T Q ξ ξ ˜ υ ˜ w T Q υ w υ ˜ w + u ˜ f T P f 1 u ˜ ˙ f
As
z ˜ 1 T c 1 z ˜ 1 + z ˜ 1 T G 1 z ˜ 2 ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 z ˜ 1 T s s T 1 β p q | z ˜ 2 | ( p / q 1 ) c 2 z ˜ 2 = z ˜ 1 T c 1 z ˜ 1 + z ˜ 1 T G 1 z ˜ 2 ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 z ˜ 1 T ( z ˜ 1 + α γ z ˜ 1 λ + 1 β z ˜ 2 p / q ) ( z ˜ 1 + α γ z ˜ 1 λ + 1 β z ˜ 2 p / q ) T 1 β p q | z ˜ 2 | ( p / q 1 ) c 2 z ˜ 2 = z ˜ 1 T c 1 z ˜ 1 + z ˜ 1 T G 1 z ˜ 2 ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) z ˜ 1 T c 1 z ˜ 1 α γ ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) z ˜ 1 T c 1 z ˜ 1 λ 1 β ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 z ˜ 1 T z ˜ 2 p / q 1 β p q diag ( | z ˜ 2 , i | ( p / q 1 ) z ˜ 1 T c 2 z ˜ 2 1 β p q diag ( | z ˜ 2 , i | ) ( p / q 1 ) α γ ( z ˜ 1 λ ) T c 2 z ˜ 2 1 β 2 p q diag ( | z ˜ 2 , i | ) ( p / q 1 ) ( z ˜ 2 p / q ) T c 2 z ˜ 2 = z ˜ 1 T G 1 z ˜ 2 ( 2 I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 z ˜ 1 T z ˜ 1 α γ ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 | z ˜ 1 | λ 1 z ˜ 1 T z ˜ 1 1 β ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 | z ˜ 2 | p / q 1 z ˜ 1 T z ˜ 2 1 β p q diag ( | z ˜ 2 , i | ( p / q 1 ) c 2 z ˜ 1 T z ˜ 2 1 β p q diag ( | z ˜ 2 , i | ) ( p / q 1 ) α γ | z ˜ 1 | λ 1 c 2 z ˜ 1 T z ˜ 2 1 β 2 p q diag ( | z ˜ 2 , i | ) ( p / q 1 ) | z ˜ 2 | p / q 1 c 2 z ˜ 2 T z ˜ 2 = [ ( 2 I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 α γ ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 | z ˜ 1 | λ 1 ] z ˜ 1 T z ˜ 1 + ( 1 β ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 | z ˜ 2 | p / q 1 1 β p q diag ( | z ˜ 2 , i | ( p / q 1 ) c 2 1 β p q diag ( | z ˜ 2 , i | ) ( p / q 1 ) α γ | z ˜ 1 | λ 1 c 2 ) z ˜ 1 T z ˜ 2 + z ˜ 1 T G 1 z ˜ 2 1 β 2 p q diag ( | z ˜ 2 , i | ) ( p / q 1 ) | z ˜ 2 | p / q 1 c 2 z ˜ 2 T z ˜ 2
and define the positive definite matrix
Q = [ Q 11 1 β ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) | z ˜ 2 | p / q 1 ( c 1 + c 2 p q ) G 1 1 β 2 p q diag ( | z ˜ 2 , i | ) ( p / q 1 ) | z ˜ 2 | p / q 1 c 2 ] where   Q 11 = [ ( 2 I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 + α γ ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 | z ˜ 1 | λ 1 ] ,   then   it   yields
z ˜ 1 T c 1 z ˜ 1 + z ˜ 1 T G 1 z ˜ 2 ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 z ˜ 1 T s s T 1 β p q | z ˜ 2 | ( p / q 1 ) c 2 z ˜ 2 = [ z ˜ 1 T z ˜ 2 T ] Q [ z ˜ 1 z ˜ 2 ]
Choose the update law as
d ^ ˙ = ( 1 β p q | z ˜ 2 | ( p / q 1 ) Γ 2 ) T s
In addition, the fault estimator is designed as
u ^ ˙ f = P f ( 1 β p q | z ˜ 2 | ( p / q 1 ) G 2 ) T z ˜ 2 γ f u ^ f
which yields
u ˜ f T P f 1 u ˜ ˙ f = u ˜ f T P f 1 ( u ^ ˙ f u ˙ f ) = u ˜ f T P f 1 P f ( 1 β p q | z ˜ 2 | ( p / q 1 ) G 2 ) T z ˜ 2 u ˜ f T P f 1 γ f u ^ f u ˜ f T P f 1 u ˙ f = u ˜ f T ( 1 β p q | z ˜ 2 | ( p / q 1 ) G 2 ) T z ˜ 2 P f 1 γ f u ˜ f T u ^ f P f 1 u ˜ f T u ˙ f
As | | u f | | u ¯ max < , | | u ˙ f | | u ˙ ¯ max < , u ¯ max and u ˙ ¯ max denote upper bounds of the actuator position and rotation rate, then
P f 1 γ f u ˜ f T u ^ f = P f 1 γ f u ˜ f T ( u ˜ f + u f ) 1 2 ( P f 1 γ f u ˜ f T u ˜ f + P f 1 γ f u ¯ max 2 )
P f 1 u ˜ f T u ˙ f 1 2 ( P f 1 u ˜ f T u ˜ f + P f 1 u ˙ ¯ max 2 )
Substitute (A21); (A22) into (A20) yields
u ˜ f T P f 1 u ˜ ˙ f u ˜ f T ( 1 β p q | z ˜ 2 | ( p / q 1 ) G 2 ) T z ˜ 2 + c ¯
where
c ¯ = 1 2 ( P f 1 γ f u ˜ f T u ˜ f + P f 1 γ f u ¯ max 2 ) + 1 2 ( P f 1 u ˜ f T u ˜ f + P f 1 u ˙ ¯ max 2 )
Let Z 12 = [ z ˜ 1 T z ˜ 2 T ] T , and thus, it finally yields
V ˙ 2 = Z 12 T Q Z 12 h s 2 ς | s | ξ ˜ T Q ξ ξ ˜ υ ˜ w T Q υ w υ ˜ w + c ¯ Z 12 T Q Z 12 h s 2 ς | s | i = 1 3 q ξ i ξ ˜ i 2 i = 1 3 q w i υ ˜ w i 2 + c ¯
where q ξ i and q w i are diagonal elements of matrices of Qξ and Q υ w , respectively. Using (A16), it is obtained that
| Q | = | [ Q 11 1 β ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) | z ˜ 2 | p / q 1 ( c 1 + c 2 p q ) G 1 1 β 2 p q diag ( | z ˜ 2 , i | ) ( p / q 1 ) | z ˜ 2 | p / q 1 c 2 ] | = [ ( 2 I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 + α γ ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 | z ˜ 1 | λ 1 ] ( 1 β 2 p q diag ( | z ˜ 2 , i | ) ( p / q 1 ) | z ˜ 2 | p / q 1 c 2 ) + 1 β ( I n + α γ λ diag ( | z ˜ 1 , i | λ 1 ) ) c 1 | z ˜ 2 | p / q 1 G 1 + 1 β p q diag ( | z ˜ 2 , i | ( p / q 1 ) c 2 ( 1 + α γ | z ˜ 1 | λ 1 ) G 1 = ( 1 β diag ( | z ˜ 2 , i | ( p / q 1 ) ) ( 1 β p q c 1 c 2 ( 2 I n + α γ ( λ + I ) | z ˜ 1 | λ 1 + α γ 2 λ | z ˜ 1 | 2 λ 2 ) ( | z ˜ 2 | p / q 1 ) + ( I n + α γ λ | z ˜ 1 | λ 1 ) ( c 1 + c 2 p q ) G 1 )
If the following condition is satisfied,
( I n + α γ λ | z ˜ 1 | λ 1 ) ( c 1 + c 2 p q ) G 1 > 1 β p q c 1 c 2 ( 2 I n + α γ ( λ + I ) | z ˜ 1 | λ 1 + α γ 2 λ | z ˜ 1 | 2 λ 2 ) ( | z ˜ 2 | p / q 1 )
then Q is a positive definite, and if ς > 0 , then
V ˙ 2 , 1 = z 12 T Q z 12 ς | s | < 0 .
Let ck = min{2Q, 2ς, –γf+1, 2 q ξ i , 2 q w i } > 0 (i = 1, 2, 3), then
V ˙ 2 c k V 2 + c ¯ .
According to LaSalle-Yoshizawa Lemma [34] and referring to [32], the closed-loop system tracking error will exponentially converge, and
V 2 ( V 2 ( 0 ) c ¯ / c k ) e c k t + c ¯ / c k
Hence
1 2 z ˜ i T z ˜ i V 2 , z ˜ i 2 c ¯ c k
that is, the desired tracking error will exponentially converge to the set z ˜ i 2 c ¯ c k , and
x 1 = z ˜ 1 + ζ 1 z ˜ 1 + ζ 1 , x 2 = z ˜ 2 + ζ 2 z ˜ 2 + ζ 2

References

  1. Khoury, G.A.; Gillett, J.D. Airship Technology, 2nd ed.; Cambridge University Press: Cambridge, UK, 2018; pp. 431–519. [Google Scholar]
  2. Alexandra, M.; Azinheira, J.R.; Paiva, E.D.; Bueno, S.S. Airship Robust Path-tracking: A Tutorial on Airship Modeling and Gain-scheduling Control Design. Contr. Eng. Pract. 2016, 50, 22–36. [Google Scholar]
  3. Liu, S.Q.; Sang, Y.G. Underactuated Stratospheric Airship Trajectory Control Using an Adaptive Integral Backstepping Approach. J. Aircr. 2018, 55, 2357–2371. [Google Scholar] [CrossRef]
  4. Liu, S.Q.; Sang, Y.J.; Whidborne, J.F. Adaptive sliding-mode-backstepping trajectory tracking control of underactuated airships. Aerosp. Sci. Technol. 2020, 97, 105610. [Google Scholar] [CrossRef]
  5. Liu, C.; Wei, Y.; Wang, Y. Finite-time Trajectory Tracking Control for Autonomous Airships with Uncertainties and External Disturbances. IET Intell. Trans. Syst. 2020, 14, 440–448. [Google Scholar]
  6. Amir, S.; Liu, Y.; Shah, Z.M.; Wang, L.; Zuo, Z.Y.; Wang,, Q.G. Higher Order Sliding Mode based Lateral Guidance and Control of Finless Airship. Aerosp. Sci. Technol. 2021, 113, 106670. [Google Scholar] [CrossRef]
  7. Yuan, J.; Zhu, M.; Guo, X.; Lou, W.J. Finite-time Trajectory Tracking Control for a Stratospheric Airship with Full-state Constraint and Disturbances. J. Frankl. Inst. 2021, 358, 1499–1528. [Google Scholar] [CrossRef]
  8. Liu, S.Q.; Whidborne, J.F. Neural Network Adaptive Backstepping Fault Tolerant Control for Unmanned Airships with Multi-vectored Thrusters. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2021, 235, 1507–1520. [Google Scholar] [CrossRef]
  9. Hou, W.; Zhou, P.; Wang, Y.; Wang, N.; Duan, D.P. Station-keeping Control of an Underactuated Stratospheric Airship. Int. J. Fuzzy Syst. 2019, 21, 715–732. [Google Scholar]
  10. Krstic, M.; Kanellakopoulo, S.I.; Kokotovic, P. Nonlinear and Adaptive Control Design; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1995; pp. 22–66. [Google Scholar]
  11. Liu, S.Q.; Whidborne, J.F. Observer-based Incremental Backstepping Sliding Mode Fault-tolerant Control for Blended-wing-body Aircrafts. Neurocomputing 2021, 464, 546–561. [Google Scholar] [CrossRef]
  12. Liu, S.; Liu, Y.; Wang, N. Nonlinear Disturbance Observer-based Backstepping Finite-time Sliding Mode Tracking Control of Underwater Vehicles with System Uncertainties and External Disturbances. Nonlinear Dyn. 2017, 88, 465–476. [Google Scholar] [CrossRef]
  13. Yang, Y.N. A Time-specified Nonsingular Terminal Sliding Mode Control Approach for Trajectory Tracking of Robotic Airships. Nonlinear Dyn. 2018, 92, 1359–1367. [Google Scholar] [CrossRef]
  14. Mofid, O.; Mobayen, S.; Zhang, C.W.; Esakki, B. Desired Tracking of Delayed Quadrotor UAV under Model Uncertainty and Wind Disturbance using Adaptive Super-twisting Terminal Sliding Mode Control. ISA Trans. 2021, 1–17, in press. [Google Scholar] [CrossRef]
  15. Alattas, K.A.; Mofid, O.; Alanazi, A.K.; Abo-Dief, H.M.; Bartoszewicz, A.; Bakouri, M.; Mobayen, S. Barrier Function Adaptive Nonsingular Terminal Sliding Mode Control Approach for Quad-Rotor Unmanned Aerial Vehicles. Sensors 2022, 22, 909. [Google Scholar] [CrossRef]
  16. Pouzesh, M.; Mobayen, S. Event-triggered Fractional-order Sliding Mode Control Technique for Stabilization of Disturbed Quadrotor Unmanned Aerial Vehicles. Aerosp. Sci. Technol. 2022, 121, 107337. [Google Scholar] [CrossRef]
  17. Vieira, H.S.; Paiva, E.C.D.; Moriguchi, S.K.; Carvalho, J.R.H. Unified Backstepping Sliding Mode Framework for Airship Control Design. IEEE Trans. Aerosp. Electr. Syst. 2020, 56, 3246–3258. [Google Scholar] [CrossRef] [Green Version]
  18. Liu, S.Q.; Whidborne, J.F.; He, L. Backstepping Sliding-mode Control of Stratospheric Airships using Disturbance-observer. Adv. Space Res. 2021, 67, 1174–1187. [Google Scholar] [CrossRef]
  19. Guo, X.G.; Xu, W.D.; Wang, J.L.; Park, J.H. Distributed Neuroadaptive Fault-tolerant Sliding-mode Control for 2-D Plane Vehicular Platoon Systems with Spacing Constraints and Unknown Direction Faults. Automatica 2021, 129, 109675. [Google Scholar] [CrossRef]
  20. Xiao, C.; Han, D.; Wang, Y.; Zhou, P.F.; Duan, D.P. Fault-tolerant Tracking Control for a Multi-vectored Thrust Ellipsoidal Airship using Adaptive Integral Sliding Mode Approach. Proc. IMechE Part G J. Aerosp. Eng. 2018, 232, 1911–1924. [Google Scholar] [CrossRef]
  21. Li, Y.W.; Nahon, M.; Sharf, I. Dynamics Modeling and Simulation of Flexible Airships. AIAA J. 2009, 47, 592–605. [Google Scholar] [CrossRef] [Green Version]
  22. Li, Y.W.; Nahon, M.; Sharf, I. Airship Dynamics Modeling: A Literature Review. Progr. Aerosp. Sci. 2011, 47, 217–239. [Google Scholar] [CrossRef]
  23. Bennaceur, S.; Azouz, N.; Abichou, A. An Efficient Modelling of Flexible Blimps: Eulerian Approach. In Proceedings of the ASME 2007 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference, IDETC/CIE 2007, Las Vegas, NV, USA, 4–7 September 2007; pp. 1–10. [Google Scholar]
  24. Bennaceur, S.; Azouz, N. Contribution of the Added Masses in the Dynamic Modeling of Flexible Airships. Nonlinear Dyn. 2012, 67, 215–226. [Google Scholar] [CrossRef]
  25. Han, D.; Yan, L.; Yan, G.Z.; Wang, X.L.; Duan, D.P. Integral Command Filtered Backstepping Control of a Flexible UAV. In Proceedings of the International Conference on Electrical Engineering, Control and Robotics, Chengdu, China, 12–14 January 2018; pp. 1–7. [Google Scholar]
  26. Zhang, Y.M.; Jiang, J. Integrated Design of Reconfigurable Fault Tolerant Control Systems. J. Guid. Contr. Dyn. 2001, 24, 133–136. [Google Scholar] [CrossRef] [Green Version]
  27. Tao, G.; Ma, X.; Joshi, S. Adaptive State Feedback and Tracking Control of Systems with Actuator Failures. IEEE Trans. Autom. Contr. 2001, 46, 78–95. [Google Scholar] [CrossRef]
  28. Zheng, Z.W.; Zhu, M.; Shi, D.L.; Wu, Z. Hovering Control for a Stratospheric Airship in Unknown Wind. In Proceedings of the AIAA Guidance, Navigation, and Control Conference, AIAA SciTech Forum, National Harbor, MD, USA, 13–17 January 2014; pp. 1–16. [Google Scholar]
  29. Van, M.; Mavrovouniotis, M.; Ge, S.S. An Adaptive Backstepping Nonsingular Fast Terminal Sliding Mode Control for Robust Fault Tolerant Control of Robot Manipulators. IEEE Trans. Syst. Man Cybern. Part B Cybern. 2019, 49, 1448–1458. [Google Scholar] [CrossRef] [Green Version]
  30. Chen, S.Y.; Lin, F.J. Robust Nonsingular Terminal Sliding Mode Control for Nonlinear Magnetic Bearing System. IEEE Trans. Control Syst. Technol. 2011, 19, 636–643. [Google Scholar] [CrossRef]
  31. Chen, W.H.; Yang, J.; Guo, L.; Li, S.H. Disturbance-Observer-Based Control and Related Methods—an overview. IEEE Trans. Ind. Electron. 2016, 63, 1083–1095. [Google Scholar] [CrossRef] [Green Version]
  32. Jin, X. Adaptive Fault Tolerant Control for a Class of Input and State Constrained MIMO Nonlinear Systems. Int. J. Robust Nonlinear Control 2016, 26, 286–302. [Google Scholar] [CrossRef]
  33. Liu, S.Q.; Whidborne, J.F.; Chumalee, S. Disturbance Observer Enhanced Neural Network LPV Control for a Blended-wing-body large Aircraft. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 2689–2703. [Google Scholar] [CrossRef]
  34. Khalil, H.K. Nonlinear Systems, 3rd ed.; Pearson Education: London, UK, 2002; p. 124. [Google Scholar]
Figure 1. Structure of the flexible airship and the reference frames.
Figure 1. Structure of the flexible airship and the reference frames.
Aerospace 09 00209 g001
Figure 2. Block diagram of the proposed BNTSM controller.
Figure 2. Block diagram of the proposed BNTSM controller.
Aerospace 09 00209 g002
Figure 3. The PID controller configuration.
Figure 3. The PID controller configuration.
Aerospace 09 00209 g003
Figure 4. Euler angles in Scenario 1 without wind compensation.
Figure 4. Euler angles in Scenario 1 without wind compensation.
Aerospace 09 00209 g004
Figure 5. Tracking responses of Euler angles with and without wind compensation.
Figure 5. Tracking responses of Euler angles with and without wind compensation.
Aerospace 09 00209 g005
Figure 6. Wind speed estimation, in Scenario 1 with wind.
Figure 6. Wind speed estimation, in Scenario 1 with wind.
Aerospace 09 00209 g006
Figure 7. Structure mode shapes in the Oxy and Oxz planes.
Figure 7. Structure mode shapes in the Oxy and Oxz planes.
Aerospace 09 00209 g007
Figure 8. Attitude responses with fault compensation.
Figure 8. Attitude responses with fault compensation.
Aerospace 09 00209 g008
Figure 9. Tracking responses of Euler angles with and without fault compensation.
Figure 9. Tracking responses of Euler angles with and without fault compensation.
Aerospace 09 00209 g009
Figure 10. Control inputs with faults.
Figure 10. Control inputs with faults.
Aerospace 09 00209 g010
Figure 11. Control surface faults estimation.
Figure 11. Control surface faults estimation.
Aerospace 09 00209 g011
Figure 12. Response of sliding-mode variable, s.
Figure 12. Response of sliding-mode variable, s.
Aerospace 09 00209 g012
Figure 13. Attitude responses with five stiffness variants, winds and control faults.
Figure 13. Attitude responses with five stiffness variants, winds and control faults.
Aerospace 09 00209 g013
Figure 14. The generalized coordinate responses of structure modes.
Figure 14. The generalized coordinate responses of structure modes.
Aerospace 09 00209 g014
Figure 15. The generalized coordinate velocity responses of structure modes.
Figure 15. The generalized coordinate velocity responses of structure modes.
Aerospace 09 00209 g015
Figure 16. Control inputs with five stiffness variants.
Figure 16. Control inputs with five stiffness variants.
Aerospace 09 00209 g016
Table 1. Structure parameters for the studied airship.
Table 1. Structure parameters for the studied airship.
ParameterValueParameterValue (kg)
Ix1.0834 × 105 (kg·m2)m2.682 × 103
Iy2.8819 × 105 (kg·m2)m115.0575 × 102
Iz2.1722 × 105 (kg·m2)m225.0737 × 103
Ixz0m335.0737 × 103
L50 (m)m446.8879 × 104
(xG, yG, zG)(0, 0, 3.605) (m)m555.8540 × 105
VB5.131 × 103 (m3)m665.8540 × 105
ET0433,440 (N/m)
Table 2. Selected parameters of the controller.
Table 2. Selected parameters of the controller.
ControllerParametersValue
PIDkp,zj, ki,zj (j = 1, 2)2, 0.04
BNTSMc1, c2diag(2, 2, 2), diag(0.2, 0.2, 0.2)
h, ς, φs0.01 * diag(1,1,1), 0.01 * diag(1, 1.2, 1), 0.4
p, q, β5, 3, 0.1 * diag(1, 1, 1)
αr, γd, λ5, 1, 1. 8 * diag(1, 1, 1)
BISMCc1, c20.1 * diag(1, 1, 1), 0.01 * diag(1, 1, 1)
h, ς0.1 * diag(1, 1, 1), diag(1, 1.2, 1)
λ1,λ20.5 * diag(1, 1, 1),0.05 * diag(1, 1, 1)
φs, ki,smc0.4, 0.1
Table 3. Tracking errors and ET of the airship under the input of the controller without wind compensation.
Table 3. Tracking errors and ET of the airship under the input of the controller without wind compensation.
Controller Error Eφ (rad)Eθ (rad)Eψ (rad)ET (s)
PID0.05320.08430.0893 × 10−32153.6
BNTSM0.05250.08340.0887 × 10−32367.2
BISMC0.05280.08390.0890 × 10−32226.2
Table 4. Tracking errors and ET of the airship under the input of the controller with wind compensation.
Table 4. Tracking errors and ET of the airship under the input of the controller with wind compensation.
Controller Error Eφ (rad)Eθ (rad)Eψ (rad)ET (s)
PID0.05770.07860.0817 × 10−31993.2
BNTSM0.05740.07740.0805 × 10−32050.4
BISMC0.05750.07800.0811 × 10−32048.9
Table 5. Tracking errors and ET of the airship under the input of the controller without fault compensation.
Table 5. Tracking errors and ET of the airship under the input of the controller without fault compensation.
Controller Error Eφ (rad)Eθ (rad)Eψ (rad)ET (s)
PID0.05100.08100.0788 × 10−32171.4
BNTSM0.05110.08030.0775 × 10−32354.3
BISMC0.05100.08050.0780 × 10−32193.0
Table 6. Tracking errors and ET of the airship under the input of the controller with fault compensation.
Table 6. Tracking errors and ET of the airship under the input of the controller with fault compensation.
Controller Error Eφ (rad)Eθ (rad)Eψ (rad)ET (s)
PID0.04880.08040.0829 × 10−31821.0
BNTSM0.04900.07970.0817 × 10−32114.1
BISMC0.04890.07980.0819 × 10−32016.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, S.; Whidborne, J.F.; Song, S.; Lyu, W. Adaptive Backstepping Nonsingular Terminal Sliding-Mode Attitude Control of Flexible Airships with Actuator Faults. Aerospace 2022, 9, 209. https://doi.org/10.3390/aerospace9040209

AMA Style

Liu S, Whidborne JF, Song S, Lyu W. Adaptive Backstepping Nonsingular Terminal Sliding-Mode Attitude Control of Flexible Airships with Actuator Faults. Aerospace. 2022; 9(4):209. https://doi.org/10.3390/aerospace9040209

Chicago/Turabian Style

Liu, Shiqian, James F. Whidborne, Sipeng Song, and Weizhi Lyu. 2022. "Adaptive Backstepping Nonsingular Terminal Sliding-Mode Attitude Control of Flexible Airships with Actuator Faults" Aerospace 9, no. 4: 209. https://doi.org/10.3390/aerospace9040209

APA Style

Liu, S., Whidborne, J. F., Song, S., & Lyu, W. (2022). Adaptive Backstepping Nonsingular Terminal Sliding-Mode Attitude Control of Flexible Airships with Actuator Faults. Aerospace, 9(4), 209. https://doi.org/10.3390/aerospace9040209

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