Next Article in Journal
An Enhanced and Secure Trust-Aware Improved GSO for Encrypted Data Sharing in the Internet of Things
Next Article in Special Issue
Experimental and Numerical Study on Stress Distribution Characteristics of Traveling Wave Resonance of High-Speed Bevel Gear in Aero-Engine
Previous Article in Journal
Synthesis of the Energy-Saving Dry Dual Clutch Control Mechanism
Previous Article in Special Issue
Labeling Expert: A New Multi-Network Anomaly Detection Architecture Based on LNN-RLSTM
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Fault-Tolerant Control of Hypersonic Vehicles with Unknown Model Inertia Matrix and System Induced by Centroid Shift

1
School of Science, Jiangsu University of Science and Technology, Zhenjiang 212000, China
2
Shanghai Key Laboratory of Aerospace Intelligent Control Technology, Shanghai Aerospace Control Technology Institute, Shanghai 201108, China
*
Author to whom correspondence should be addressed.
Submission received: 15 November 2022 / Revised: 30 December 2022 / Accepted: 31 December 2022 / Published: 6 January 2023

Abstract

:
In this paper, the fault-tolerant control problem of hypersonic vehicle (HSV) in the presence of unexpected centroid migration, actuator failure and external interference is studied in depth. First, the proposed dynamics for HSV with the aforementioned unexpected factors are modeled to demonstrate the peculiar nature of the subject under study. The adverse effects of accidental centroid migration are mainly reflected in the following aspects: (1) the change of inertia matrix of the system, (2) the uncertainty of the system and (3) the eccentric moment, which are coupled and unknown. Subsequently, to account for the effect of unexpected centroid shifts, a sliding-mode observer and an adaptive estimator are designed to obtain unknowns useful for subsequent FTC controller designs. Later, we derived an innovative adaptive FTC scheme by employing the observer in conjunction with a specific adaptive controller consisting of a sixth-order dynamic compensator, which can guarantee the achievement of the control objective without resorting to the exact knowledge of the inertial matrix. Moreover, the analysis of boundedness with respect to the entire signal in this closed system is performed by means of the Lyapunov stability theory. Ultimately, simulation results show that the proposed FTC strategy is efficient and powerful.

1. Introduction

As a deeply coupled nonlinear aircraft, hypersonic vehicles have the advantages of fast speed, strong maneuverability and good penetration. It has important applications in both military and civil fields, attracting many scholars to study its modeling and controller design [1,2]. The super high speed flight and corresponding complex flight airspace make HSV prone to failure, such as actuator fault, sensor fault, etc. Thus, fault-tolerant control (FTC) systems are central to ensuring HSV safety and achieving its goals.
In a fault event on an aircraft, a significant portion of the aircraft’s aerodynamic lift surface may become detached when it is damaged [3], resulting in an unexpected shift of the aircraft’s center of mass. It is well known that the position of the center of gravity of an aircraft is crucial to the calculation of its control moment [4]. Due to changes in the system, inertia matrix, eccentricity torque and system uncertainty, any unexpected center of gravity movement may put the aircraft in a non-nominal flight state, thus challenging the existing control system to maintain its control performance [3]. A similar situation may arise from fuel consumption, moving mass control [5] as well as the liquid fuel sloshing [6], etc. To ensure the safe flight of the aircraft, adaptive FTC schemes have been developed even in the presence of unexpected centroid shifts, actuator faults and external disturbances.
For variants of aircraft centroids, some previous works deliver the following. An aircraft with unexpected centroid shift is modeled in [7], which is derived from the aircraft body coordinate system, which is not convenient to analyze the influence of centroid migration on HSV flow Angle (Angle of attack and Angle of slip). In addition, a hybrid direct–indirect adaptive control method based on neural network is proposed to deal with the damage-induced centroid migration [7,8]. However, irreversible phenomena may appear in the estimated inertia matrix, leading to singularities in the existing FTC controller resting on the inverse of an inertia matrix. In addition, a linearized aircraft model with centroid shift due to damage is proposed in [9]. Nevertheless, the linearization of the multi-model reference adaptive control makes the system free from the coupling problem caused by centroid migration, leading to a degradation of the control performance. Based on the above studies, the influence of the centroid shift on the HSV motion is mainly reflected in three aspects: (1) the change of inertia matrix of the system, (2) the uncertainty of the system and (3) the eccentric moment, especially in the situation where the inertia matrix appears irreversible. In order to handle the challenge of inertia matrix, RBFNN is adopted in [8] to estimate the inverse of the inertia matrix in order to ensure the normal operation of the controller. To solve this problem, Taylor expansion of the inverse matrix of the inertial matrix is used in [10], but it comes with corresponding Taylor expansion errors and a large number of online computational parameters of the neural network. In addition, Ref. [11] addresses the variation of an inertial matrix by assuming that the minimum and maximum bounds of its eigenvalues are known; however, the mentioned bounds are difficult to measure. With the development of data driven methods [12], a model-free control scheme is developed to handle this challenge, but the training data of the inertial matrix of HSV with centroid-shift is difficult to obtain. Additionally, an observer-based FTC strategy [13], fuzzy approximation [14] and iterative learning [15] are developed to handle this challenge. In summary, how to identify an HSV inertial matrix with uncertainty online and ensure its function in the FTC controller remains elusive.
In the past few decades, the fault-tolerant control of actuators, sensors, control surfaces or other failure systems has been extensively studied, and a large number of FTC control strategies have appeared in the literature. For example, predictive control [1], backstepping control [16], adaptive control [3], fuzzy control [17], control [18] based on neural network, etc. In [16], a second-order system of actuators is established to describe actuator failure (LOE) and stuck faults, which provides a basis for subsequent fault diagnosis and identification (FDI) design.Then, fault-tolerant control strategies resting on FDI and backstepping control are developed to achieve the control objective; however, such accurate actuator models are difficult to establish. In [17], the Takagi–Sugeno (T–S) fuzzy model is applied to describe the attitude dynamics of a near-space vehicle in the case of actuator failure. Then, the adaptive fault-tolerant tracking control is implemented by using the estimated actuator faults, but there is no suitable fuzzy rule to describe the faulty aircraft. In addition, an adaptive output feedback fault-tolerant controller is designed in [19] to handle parameter uncertainties, actuator failures and external disturbances in HSV attitude systems. Therefore, such a large number of adaptive learning parameters can create a heavy online computational load. In addition, based on the strong robustness of sliding mode control and its modified control strategies are more applied in fault-tolerant HSV [20,21]. A new adaptive nonlinear generalized predictive controller is proposed for the longitudinal dynamics of hypersonic vehicles with unknown parameters and constraints [1], showing a good control effect. Nonlinear generalized predictive control [22] shows the good dynamic performance and hugely decreases the computation burden compared with the usual model predictive control. Nonetheless, since inertial matrix variations are not taken into account, they may end up delivering unacceptable performance for HSV attitude control in the face of an unexpected centroid shift. In [23,24], a six-order dynamic compensation control law without considering the inertia of HSV is proposed, which provides a modern treatment method for the change of inertia matrix. However, systematic and eccentric moment uncertainties are not taken into account.
In conclusion, it is significant and important to study the FTC scheme of a hypersonic vehicle affected by unexpected centrocentric migration, actuator failure and system input saturation. The main contributions are as follows:
  • The topic of pose system modeling of rigid HSV with unexpected centroid shifts and actuator faults is explored in detail. In contrast to the common practice in the FTC literature for actuator failures, the challenge of the proposed FTC design manifests itself in the simultaneous variation of the system inertia matrix, the moment of eccentricity and the strong system uncertainty.
  • Due to the above-mentioned coupling and unknown adverse effects, a higher-order sliding mode observer in slow-loop and the adaptive estimators are designed by the nature of eccentric moments and system uncertainties and the actuator faults instead of merging them into one item to improve the control performance of the system.
  • The algorithm uses an adaptive observer and estimator, combined with adaptive feedback control, and can achieve attitude tracking without knowing the inertia matrix J, despite the existence of system inertia matrix, eccentric moment and actuator failures. In addition, the developed FTC scheme has fixed control structure and few design parameters, which is an effective method to solve similar problems.
In this paper, we focus on fault-tolerant control of HSV in the presence of unexpected centroids, actuator faults and external disturbances. The paper is organized in the following subsections: Section 2 is devoted to the modeling of HSV with unexpected centroid shifts and the analysis of their properties, nd some details with respect to the process of modeling are delivered in Appendix A and Appendix B. Section 3 introduces an adaptive fault-tolerant control strategy, which does not require an accurate understanding of the inertial matrix, followed by Section 4, where numerical simulations are performed to assess the effectiveness of the designed FTC algorithm. Finally, the conclusions are provided in Section 5.

2. Problem Formulation

2.1. Nonlinear Attitude Model of HSV with Unexpected Centroid Shift

In order to study the effect of an unexpected centroid shift on the motion of the HSV’s attitude system, a detailed study with respect to HSV attitude dynamics with unexpected centroid shift is presented, where the centroid shift refers to the phenomenon that HSV’s original centroid ( O b ) moves to the current position ( O b ), as shown in Figure 1. Based on previous studies [7,9,25], the effect of the unexpected centroid shift can be transmitted in the following two affine subsystems: (1) attitude angle system and (2) attitude angular rate system.

2.1.1. Attitude Angles Affine Nonlinear Equations of HSV with Centroid Shift

According to [8,25], the attitude angle subsystem of HSV with unexpected centroid shift can be developed as:
α ˙ = q ( p cos α + r sin α ) tan β + ( c 31 Δ x + c 32 Δ y + c 33 Δ z ) ( c 31 Δ x + c 32 Δ y + c 33 Δ z ) ( m V cos β ) ( m V cos β )             + ( L + m g cos γ cos μ / ( m V cos β )
β ˙ = p sin α r cos α + ( Y + m g cos γ sin μ ) ( Y + m g cos γ sin μ ) m V m V + ( c 21 Δ x + c 22 Δ y + c 23 Δ z ) ( c 21 Δ x + c 22 Δ y + c 23 Δ z ) m V m V
μ ˙ = sec β ( p cos α + r sin α ) tan β ( L + m g cos γ cos μ ) tan β ( L + m g cos γ cos μ ) m V m V + tan γ ( σ 1 σ 2 + σ ˙ 3 ) ρ ¯ x p + ( σ 1 2 + σ 3 2 ) ρ ¯ y p ( σ 2 σ 3 σ ˙ 1 ) ρ ¯ z p tan γ ( σ 1 σ 2 + σ ˙ 3 ) ρ ¯ x p + ( σ 1 2 + σ 3 2 ) ρ ¯ y p ( σ 2 σ 3 σ ˙ 1 ) ρ ¯ z p V V tan β ( c 31 Δ x + c 32 Δ y + c 33 Δ z ) tan β ( c 31 Δ x + c 32 Δ y + c 33 Δ z ) m V m V + tan γ ( Y cos μ + L sin μ ) tan γ ( Y cos μ + L sin μ ) m V m V
where V is the velocity vector of the system and γ represents the flight path angle. α , β and μ are the attack angle, sideslip angle and roll angle, respectively, and ω = [ p , q , r ] T stands for the attitude angular rate vector that consists of rolling rate (p), pitching rate (q) and yawing rate (r). L and Y are the lift and side force of HSV. As for the influences of unexpected centroid shift, ρ ¯ = [ Δ x , Δ y , Δ z ] T denotes the variation of the center-of-mass location ( O b ) with respect to the original centroid ( O b , as shown in Figure 1). Moreover, C i , j , i , j = 1 , 2 , 3 are intermediate variables, shown in (A1), and σ 1 , σ 2 and σ 3 are the angular velocities in view of flight-path coordinate frame with respect to the earth-surface inertial coordinate frame, provided in (A2). ρ ¯ x p , ρ ¯ y p and ρ ¯ z p are the functions of Δ x , Δ y and Δ z , standing for the influence of centroid shift on the system of HSV, seen in (A3). Finally, m is the mass of HSV.
In what follows, in order to facilitate the subsequent theoretical analysis, based on (1)–(3), the affine attitude angle subsystem of HSV is obtained as
Ω ˙ = f s + Δ f s + g s ( I F ) ω c + d s
where Ω = [ α , β , μ ] T is defined to describe the attitude angle vector. f s = [ f α f β f μ ] T denotes the system matrix, and Δ f s = [ Δ f α Δ f β Δ f μ ] T represents the uncertainty of f s . In addition, ω c is the tracking index vector of attitude angular rate based on attitude angle tracking error. F stands for the fault factor of the slow-loop controller and the details are seen in (5), which represents a kind of signal gain fault [16]. D s = Δ f s + g s F ω c + d s denotes the influence of centroid shift, external disturbance and virtual controller gain fault of slow-loop of HSV, satisfying that D s is bounded. The details are described as follows:
f α = ( L + m g cos γ cos μ ) ( L + m g cos γ cos μ ) ( m V cos β ( m V cos β ) , f β = ( Y + m g cos γ sin μ ) ( Y + m g cos γ sin μ ) m V m V f μ = ( Y cos μ + L sin μ ) tan γ tan γ m V p m V p ( L + m g cos γ cos μ ) tan β ( L + m g cos γ cos μ ) tan β m V p m V p Δ f α = ( c 31 Δ x + c 32 Δ y + c 33 Δ z ) ( c 31 Δ x + c 32 Δ y + c 33 Δ z ) ( m V cos β ) ( m V cos β ) , Δ f β = ( c 21 Δ x + c 22 Δ y + c 23 Δ z ) ( c 21 Δ x + c 22 Δ y + c 23 Δ z ) m V m V Δ f μ = tan γ ( σ 1 σ 2 + σ ˙ 3 ) ρ ¯ x p + ( σ 1 2 + σ 3 2 ) ρ ¯ y p ( σ 2 σ 3 σ ˙ 1 ) ρ ¯ z p tan γ ( σ 1 σ 2 + σ ˙ 3 ) ρ ¯ x p + ( σ 1 2 + σ 3 2 ) ρ ¯ y p ( σ 2 σ 3 σ ˙ 1 ) ρ ¯ z p V p V p ( c 31 Δ x + c 32 Δ y + c 33 Δ z ) tan β ( c 31 Δ x + c 32 Δ y + c 33 Δ z ) tan β m V p m V p g s = tan β cos α 1 tan β sin α sin α 0 cos α sec β cos α 0 sec β sin α , F i t T = F i ( t T ) = 0 , F i ( t T ) = 1 , t T t T
T denotes the time of the fault, which occurs at some unknown time. F i is the signal gain fault of the slow-loop controller, satisfying 0 < F i < 1 . When F i = 0 , the output of ω c is in normal condition; while at F i = 1 , the output of ω c is entirely ineffective.

2.1.2. Attitude Angular Rates Affine Nonlinear Equations of HSV with Centroid Shift

In standard aircraft modeling, the mass distribution of an aircraft is assumed to be symmetric with respect to the x–z plane, such that the cross products of the inertial forces J x y and J y z are zero [26]. However, with the unexpected centroid shift ( O b moving to the original location O b , as shown in Figure 1), J x y and J y z become nonzero leading to the following moment equations for HSV [7,9,25]:
J x x p ˙ J x y q ˙ J x z r ˙ + J x y p r J x z p q + ( J z z J y y ) q r + ( r 2 q 2 ) J y z + m ( ( p V A q U A + W ˙ A ) Δ y + ( p W A r U A V ˙ A ) Δ z ) = L + m g cos θ cos ϕ Δ y m g cos θ sin ϕ Δ z J x y p ˙ + J y y q ˙ J y z r ˙ + J y z p q J x y q r + ( J x x J z z ) p r + ( p 2 r 2 ) J x z + m ( ( q U A p V A W ˙ A ) Δ x + ( q W A r V A + U ˙ A ) Δ z ) = M m g cos θ cos ϕ Δ x m g sin θ Δ z J x z p ˙ J y z q ˙ + J z z r ˙ + J x z q r J y z p r + ( J y y J x x ) p q + ( q 2 p 2 ) J x y + m ( ( r U A p W A + V ˙ A ) Δ x + ( r V A q W A U ˙ A ) Δ y ) = N + m g cos θ sin ϕ Δ x + m g sin θ Δ y
where J x x , J y y and J z z represent the moments of inertia. J x y , J x z and J y z are the products of inertia of HSV. In addition, [ U A V A W A ] T and [ U ˙ A V ˙ A W ˙ A ] T are the states of HSV with respect to the aircraft body frame, as shown in (A4). For the convenience of theoretical analysis, v A = [ U A V A W A ] T is introduced as the velocity vector of HSV.
Then, through a close inspection of (6), the affine attitude angular rate subsystem of HSV can be reshaped as
ω ˙ = J 1 ω × J ω J 1 ( m ρ ¯ × v ˙ A + m ω × ( ρ ¯ × v A ) + m v A × ( ω × ρ ¯ ) ρ ¯ × W ) + J 1 ( τ + Δ τ ) + J 1 d f = J 1 ω × J ω + J 1 τ J 1 f ( ρ ¯ ) + J 1 d f + J 1 Δ τ
where J = J * + Δ J represents the inertial matrix of HSV, where J * stands for the inertial matrix without centroid shift and Δ J is the uncertainties of inertial matrix due to centroid shift, and the details are explained later. d f is the external disturbance. Furthermore, W = [ sin θ cos θ sin ϕ cos θ cos ϕ ] T is the gravitational force vector.

2.2. Analysis of the Influence of the Centroid Shift on HSV Motion

Before proceeding further, a relative analysis is presented to explain the challenges of this paper. According to previous research about the influence of centroid shift [3,25] and a close inspection of (4) and (7), it can be concluded in the following three aspects: (1) variation of inertial matrix; (2) eccentric moment; and (3) system uncertainty, and the details are delivered as follows:

2.2.1. Variation of Inertial Matrix

As for the Δ J R 3 × 3 , it represents an uncertain portion of J due to the centroid shift (from O b to O b , as shown in Figure 1). Details are given as follows:
J = J * + Δ J = J x x 0 J x z 0 J y y 0 J x z 0 J z z + Δ J x x J x y 0 J x y Δ J y y J y z 0 J y z Δ J z z
where the details of Δ J are shown as (A5). It is worth pointing out that changes in the inertia matrix of the system ( J * + Δ J ), due to the potential irreversibility of ( J * + Δ J ), will pose challenges to the controller that depends on the inverse matrix of the inertia matrix of the system. Such as in sliding mode control [16], robust adaptive control [27], and so on, resulting in the singularity of the controller or giving rise to the underactuated characteristic of the system. Existing controllers struggle to function properly.

2.2.2. Eccentric Moment

The eccentricity Δ τ due to the accidental centroid shift is derived as follows
Δ τ = 0 Z a Y a Z a 0 X a Y a X a 0 Δ x Δ y Δ z = Θ × ρ ¯
where X a , Y a and Z a are the three components of the system control input vector τ along the aircraft body frame O b x y z , as illustrated in Figure 1. The skew symmetric matrix Θ × is composed of X a , Y a and Z a and the relative details are displayed as [28]
X a = D m g sin γ , Y a = Y + m g cos γ sin μ , Z a = L + m g cos γ cos μ
where D, Y and L are the drag, side force and lift, respectively, which are shown as (A8) and (A9).

2.2.3. The Uncertainty of the System

In contrast to the usual decoupling of system uncertainty from the normal case, by multiplying ( J * + Δ J ) 1 to both sides of (7), it is obvious that the system uncertainty caused by Δ J is difficult to decouple from ( J * + Δ J ) 1 ω × ( J * + Δ J ) ω .
In addition, the remaining part of uncertainty is remarked as [7,8]
f ( ρ ¯ ) = m ρ ¯ × v ˙ A + m ω × ( ρ ¯ × v A ) + m v A × ( ω × ρ ¯ ) ρ ¯ × W
where ρ ¯ = [ Δ x , Δ y , Δ z ] T represents the offset vector of the center-of-mass. Furthermore, when [ Δ x , Δ y , Δ z ] T 0 , J x y and J y z are non-zero, which means that ( J * + Δ J ) evolves from the decoupling matrix J * to the coupling matrix J * , resulting in a strong coupling between the longitudinal and transverse dynamics of the HSV.

2.3. Actuator Faults

In this work, to describe actuator faults, the control inputs of the system (7) are further modified as [29]
τ i = λ i ( 1 ρ i ) u τ i + ρ i u ¯ τ i
where u τ i indicates the designed controller.
ρ i t t f = ρ i ( t t f ) = 0 , ρ i ( t t f ) = 1 , t t f t t f
Therefore, when t > t f , multiple critical and common fault modes can appear as
  • Stuck or lock in place (LIP) fault: ρ i = 1 , τ i = u ¯ τ i ( t f ) is a constant, and τ i denotes the disturbance moment;
  • Loss of effectiveness (LOE) fault: ρ i = 0 , τ i = λ i u τ i ( t ) , 0 < λ i < 1 , with λ i being a positive constant.

2.4. Control Objective

The purpose of this paper is to propose an adaptive compensation-fault-tolerant control scheme that enables the output signal Ω to follow the reference signal Ω c without resorting to the exact information of inertia matrix and despite the occurrence of eccentric moments, uncertain system parameters, actuator failures and external disturbances. In addition, signals in a closed-loop system are bounded. To do this, the following assumptions are required.
Assumption 1.
Before and after the occurrence of the unexpected centroid shift, the aerodynamic configuration of the HSV remains the same, that is, there is no variation in the aerodynamic parameters.
Assumption 2.
The inverse of the inertial matrix ( J * + Δ J ) exists.
Assumption 3.
The complex disturbance D s satisfies D s κ s , where κ s is an unknown positive real constant. Besides, the uncertainty of quick-loop system f ( ρ ¯ ) meets f ( ρ ¯ ) κ f , where κ f is an unknown constant, satisfying κ f > 0 . The eccentric disturbance moment of HSV satisfies Δ τ κ τ , where κ τ is the unknown positive real constant. In addition, the disturbance of quick-loop d f is also bounded.
Assumption 4.
The desired tracking command signals Ω r , ω c and their time-derivative Ω ˙ r , ω ˙ c is continuous and bounded.
Remark 1.
Under Assumption 2, the inertia matrix ( J * + Δ J ) is an invertible matrix, which plays a crucial role in controller; however, the inverse function of ( J * + Δ J ^ ) may not exist because Δ J ^ is estimated by the estimator. For this challenge, a special negative feedback control is developed in this work without resorting to the exact knowledge of ( J * + Δ J ^ ) 1 . In addition, based on an auxiliary integrator [30] or necessary gain [31], they can be another way to deal with this problem.
Before the fault-tolerant controller design, a necessary lemma is delivered as follows.
Lemma 1
([32]). The higher-order sliding mode differentiator is defined as
z ˙ 0 = υ 0 = κ 0 z 0 f ( t ) ( n n ( n + 1 ) ( n + 1 ) ) sgn ( z 0 f ( t ) ) + z 1 z ˙ i = υ i = κ i z i υ i 1 ( n i n i ( n i + 1 ) ( n i + 1 ) ) sgn ( z i υ i 1 ) + z i + 1 z ˙ n 1 = υ n 1 = κ n 1 z n 1 υ n 2 1 12 2 sgn ( z n 1 υ n 2 ) + z n z ˙ n = κ n sgn ( z n υ n 1 )
where z i and υ i are the states of higher-order sliding mode differentiator (13) and κ i is the design parameter. Assisted by the proper parameters, z 0 = f ( t ) and z i = υ i = f ( i ) ( t ) , i = 1 , , n are true with no input noise after the finite time t σ of the transient process. In addition, the corresponding solution of this dynamic system is Lyapunov stable, i.e., namely finite time stability.

3. Fault-Tolerant Control Design for the Attitude System of HSV

As shown in Figure 2, the fault tolerant control framework with unexpected centroid migration is divided into two prats: slow-Loop control strategy and slow-Loop control strategy. In slow-loop of HSV, the control strategy is consisted of nonlinear generalized predictive control scheme and a high-order sliding mode observer. In fast-loop of HSV, an adaptive feedback fault-tolerant controller that does not need the exact knowledge of ( J * + Δ J ) 1 is proposed. The details are shown as follows.

3.1. Fault-Tolerant Control Design for Slow-Loop of HSV

To achieve control target, the tracking error vector ( e s ) of quiet-loop is defined as
e s = Ω Ω r = [ e α , e β , e μ ] T = [ α α r , β β r , μ μ r ] T
where Ω r = [ α r , β r , μ r ] T represents the expected index.

3.1.1. Baseline Controller under Nominal Condition

In this work, the nonlinear generalized predictive control, one kind of optimal control is utilized to achieve the control target of the low-loop of HSV due to its excellent dynamic performance [22]. Initially, the receding horizon performance metric was introduced as
J = 0 T e ^ s T t + e ^ s t + d 0 T e ^ T t + τ e ^ s t + τ d τ 2 2
where T is the predictive period and e ^ s represents the predictive tracking error that is defined as
e ^ s t + = y ^ t + y ^ r t +
where y ^ t + is the predictive output, and y ^ r t + stands for the reference output.
  • Obtaining future behaviors. To this aim, by resorting to the Taylor series expansion, the future behavior of y ^ t + can be predicted by
    y ^ t + = Γ Y ¯ t
    where Γ = 1 , ¯ , , ¯ ρ + r ρ + r ! R m × m ρ + r + 1 , = d i a g { , , } R m × m , Y ¯ = [ y ^ , y ^ ˙ , , y ^ [ ρ + r ] ] T , as shown in (A6), where ρ is the relative degree corresponding to the system output, and r is the control order [22]. Then, the index signal y r is also predicted by
    y ^ r t + = Γ Y ¯ r t
    where Y ¯ r = [ y ^ r , y ^ ˙ r , , y ^ r [ ρ + r ] ] T , seen in (A6). From this, the tracking error associated with the predicted behavior is derived as
    e ^ ( t + ) = Γ ( Y ¯ Y ¯ r )
    Moreover, by some algebraic manipulations and rearrangements of the resulting formulas, the prediction error in the performance index (15) is rewritten as
    J = 1 2 ( Y ¯ ( t ) Y ¯ r ( t ) ) T 0 T Γ T ( ) Γ ( ) d ( Y ¯ ( t ) Y ¯ r ( t ) )
    For analysis purpose, Γ ¯ is defined as Γ ¯ = 0 T Γ T ( ) Γ ( ) d and its matrix blocks are split by
    Γ ¯ = Γ ¯ ρ ρ Γ ¯ ρ r Γ ¯ ρ r Γ ¯ r r
  • Minimizing performance index. In order to minimize the performance index (20) through the proper u ^ ( z ) , the optimal u ^ ( z ) is solved by the following partial differential equation.
    J u ^ z u z = u z * = 0
    According to (20)–(22), it is further calculated as
    ( H ˜ ( u ^ z ) u ^ z ) T [ Γ ¯ ρ r T Γ ¯ r r ] M + ( H ˜ ( u ^ z ) u ^ z ) T Γ ¯ r r H ˜ ( u ^ z ) = 0
    where M = [ h ( x ) y r ( t ) , L f s h ( x ) y ˙ r ( t ) , , L f s ρ + r h ( x ) y r ρ + r ( t ) ] T R m ( ρ + r + 1 ) , and the details of H ˜ ( u ^ z ) u ^ z are seen in (A7).
    Based on the assumption that L g s L f s ρ 1 h ( x ) and H ˜ ( u ^ z ) u ^ z are reversible, it is obtained that
    H ˜ ( u ^ z ) = [ Γ ¯ r r 1 Γ ρ r T I m ( r + 1 ) × m ( r + 1 ) ] M
    By invoking (A6), the front m rows of (24) can be described as follows
    L g s L f s ρ 1 h ( x ) u ^ z ( t ) + K M ρ + L f s ρ h ( x ) y r [ ρ ] + L f s ρ 1 h ( x ) D s = 0
    where K represents the front m rows of Γ r r 1 Γ ρ r T .
    M ρ = [ ( h ( x ) y r ( t ) ) T , ( L f s h ( x ) y ˙ r ( t ) ) T , , ( L f s ρ 1 h ( x ) y r ρ 1 ( t ) ) T ] T R m ρ .
    Eventually, the optimal control law u ^ z ( t ) can be obtained as [33]
    u ^ z ( t ) = ( L g s L f s ρ 1 h ( x ) ) 1 ( L f s ρ h ( x ) + K M ρ y r [ ρ ] + L f s ρ 1 D s )
Remark 2.
The actual control input u z ( t ) , is obtained by the initial value of the optimal control input u ^ z ( t + ) , 0 T . Namely, u z ( t + ) = u ^ z ( t + ) when = 0 [22].
Thereafter, a closed inspection of (4) reveals that relative degree ρ = 1 , and thus, by keeping (4) and (26) in mind and setting control order r = 0 , the final baseline controller under nominal condition can be readily designed as
u ( t ) = ( g s ) 1 ( f s + K e s Ω ˙ r + D s )
where the unknown complex disturbance D s of (4) can give rise to the huge challenges for the control performance of the designed controller (27). To deal with such problems, fault-tolerant controllers with the help of observers are utilized in this work, the details of which are communicated in the following sections.

3.1.2. Fault-Tolerant Controller Design under Unexpected Centroid Shift and System Fault

For the adverse effects caused by complex disturbance D s , as shown in (4), a top-order sliding mode observer is designed in this work to obtain unknown knowledge serving for the basic controller to realize the fault-tolerant control targets.
  • High-order sliding mode observer design. To estimate the complex perturbation D s , a top-order sliding mode observer is designed as follows [21]
    σ 0 = Ω z , z ˙ = f s + g s ω + D ^ s , σ 1 = σ 0 + θ σ ˙ 0 + ζ σ ˙ 0 p q D ^ ˙ s = ( θ + p q ζ d i a g ( σ ˙ 0 p q 1 ) ) 1 ( σ ˙ 0 + l 1 σ 1 + l 2 sgn ( σ 1 ) ) + ϕ ^ sgn ( σ 1 ) ϕ ^ ˙ = ϑ ( i = 1 3 ( θ i + p q ζ i d i a g ( σ ˙ 0 i p q 1 ) ) σ 1 i ) , i = 1 , 2 , 3
    where σ 0 and σ 1 are the designed sliding mode surfaces based on the terminal sliding mode trick, and z denote the state of the designed observer. θ and ζ are the designed system matrix, satisfying θ = d i a g { θ 1 , θ 2 , θ 3 } , θ i > 0 and ζ = d i a g { ζ 1 , ζ 2 , ζ 3 } , ζ i > 0 . In addition, p and q are odd positive integers, satisfying 0 < p < 2 q . Furthermore, l 1 and l 2 are also the designed positive constants to regulate the estimation effect of D s . Finally, ϕ ^ is the estimation of the adaptive parameter ϕ .
    Before analyzing the stability of the disturbance estimate error σ 0 , combined with (4) and (28), the derivative of σ 0 is deduced as
    σ 0 ˙ = Ω ˙ z ˙ = D s D ^ s , σ ¨ 0 = D ˙ s D ^ ˙ s
    Furthermore, by keeping the exponential reaching law σ ˙ 1 = l 1 σ 1 l 2 sgn ( σ 1 ) and (29) in mind, σ ˙ 1 can be obtained as
    σ ˙ 1 = l 1 σ 1 l 2 sgn ( σ 1 ) + ( θ + p q ζ d i a g ( σ ˙ 0 p q 1 ) ) ( D s ˙ ϕ ^ sgn ( σ 1 ) )
    In what follows, the stability of the designed observer is analyzed through the following Lyapunov function
    V σ = σ 1 T σ 1 σ 1 T σ 1 2 2 + ϕ ˜ 2 ϕ ˜ 2 2 ϑ 2 ϑ
    where ϑ is a designed parameter, satisfying ϑ > 0 , and we define that ϕ ˜ = ϕ ϕ ^ .
    Evaluating the time derivative of V σ along (28) and (30) leads to
    V ˙ σ = σ 1 T σ ˙ 1 ϕ ˜ ϕ ^ ˙ ϕ ˜ ^ ˙ ϕ ϑ ϑ l 1 σ 1 T σ 1 l 2 i = 1 3 σ 1 i T + i = 1 3 ( θ i + p q ζ i d i a g ( σ ˙ 0 i p q 1 ) ) ( ϕ i σ 1 i ϕ ^ i σ 1 i ) ϕ ˜ ϕ ^ ˙ ϕ ˜ ^ ˙ ϕ ϑ ϑ
    where i = 1 , 2 , 3 . Furthermore, invoking ϕ ^ ˙ = ϑ ( i = 1 3 ( θ i + p q ζ i d i a g ( σ ˙ 0 i p q 1 ) ) σ 1 i ) , (32) can be further deduced as
    V ˙ 1 l 1 σ 1 T σ 1 l 2 i = 1 3 σ 1 i T
    The differential term σ ˙ 0 in the designed higher-order sliding-mode observer is not directly obtained by using the derivative method. Moreover, according to Lemma 1, applying a higher-order sliding mode differentiator yields the corresponding derivative. Then, the states z converge to σ 0 in finite time t σ .
    • According to (31), if σ 1 i 0 , i = 1 , 2 , 3 at time t σ , the states will arrive at the sliding mode σ 1 i = 0 in finite time t 1 i , namely
      t 1 i = 1 l 1 [ ln ( σ 1 i ( t σ ) + l 2 l 1 ) ln ( l 2 l 1 ) ]
      Examining (34), it is straightforward to show that the sliding mode σ 1 = 0 can be reached within the finite time t a 1 = max ( t 11 , t 12 , t 13 ) .
    • When the sliding mode σ 1 = 0 is arrived, the convergence time of the system is up to the terminal sliding mode equation
      σ 1 = σ 0 + θ σ ˙ 0 + ζ σ ˙ 0 p p q q = 0
      where σ 0 i = 0 denotes a terminal attractor of the system of (35). According to [21], the convergence time t 0 i from t σ + t a 1 when σ 0 i ( t σ i + t a 1 i ) 0 to t σ + t a 1 + t 0 i when σ 0 i ( t σ + t a 1 + t 0 i ) = 0 satisfies
      t 0 i ( ζ i ) q q p p p p q σ 0 i ( t σ + t a 1 ) 1 q p
      In a similar way, from the time point t σ + t a 1 when σ 0 ( t σ i + t a 1 i ) 0 move to the time point t σ + t a 1 + t a 0 when σ 0 ( t σ i + t a 1 i ) = 0 , the convergence time t a 0 is obtained as t a 0 = max ( t 01 , , t 03 ) .
    From this, we can obtain that the estimated value D ^ s be approximated by an external perturbation D s in finite time t σ = max ( t a 0 , t a 1 ) . Combined with Assumption 3, we can conclude that the disturbance estimation error is bounded and can be defined as D s D ^ s < λ , where λ is a positive unknown constant.
  • Fault-tolerant controller design of slow-loop. Recall that Assumption 3, together with the estimate of complex disturbance D ^ s from observer (28) and the baseline controller (27), the final control law and adaptive law are synthesized as
    u ( t ) = ( g s ) 1 ( f s + K e s Ω ˙ c + D ^ s + λ ^ ) λ ^ ˙ = τ s e s T P s
    where ( g s ) 1 D ^ s is applied to suppress uncertainties due to unexpected centroid shifts, external perturbations, etc. As for the adaptive term λ , it is applied to handle the estimation error of D s , where λ ^ represents the estimation of λ .
To analyze the stability of the observer and the resulting closed-loop system, the following theorem summarizes the main results of this subsection.
Theorem 1.
Consider the slow-loop model of HSV describer by (4), with unexpected centroid shift, system fault and disturbance, satisfying Assumptions 1–4, under the designed controller (37) and the observer (28), the slow-loop system of HSV can track expected signals and the whole signals in closed-loop system are bounded.
Proof. 
Consider the following Lyapunov function
V s = 1 2 e s T P s e s + 1 2 τ s λ ˜ T λ ˜ + V σ
where P s is a positive definite symmetric matrix, and V σ = σ 1 T σ 1 σ 1 T σ 1 2 2 + ϕ ˜ 2 ϕ ˜ 2 2 ϑ 2 ϑ is applied to analyze the estimate error of D s .
Take the derivative of (38) and we have
V ˙ s = 1 2 e ˙ s T P s e s + 1 2 e s T P s e ˙ s + 1 τ s λ ˜ T λ ˜ ˙ + V ˙ σ = 1 2 e s T ( K T P s P s K ) e s + e s T P s ( D s D ^ s ) e s T P s λ ^ + 1 τ s λ ˜ T λ ^ ˙ + V ˙ σ 1 2 e s T Q e s + e s T P s ( λ λ ^ ) + 1 τ s λ ˜ T λ ^ ˙ + V ˙ σ 1 2 e s T Q e s l 1 σ 1 T σ 1 l 2 i = 1 3 σ 1 i T
Based on the convergence analysis of V σ , as shown in (31)–(36), combined with (39), we have that
V ˙ s a j i = 1 3 e s i 2
Therefore, due to the positive definition of V s , V ˙ s is negative semi-definite, which means V s is an increasing function of time; that is, the tracking error e s , λ ˜ , σ 1 and ϕ ˜ are all bounded. Next, the parameter projection algorithm is applied to modify the adaptive updating laws in both observer and controller to guarantee the boundedness of the designed adaptive parameters D ^ s , λ ^ , Φ ^ . Furthermore, due to e ˙ s = K e s + D s D ^ s λ ^ , combined with Assumptions 3 and 4, we can infer that e ˙ s is also bounded. Subsequently, based on the Barbalat lemma and the bounded e ˙ s and definition of V s , it can be obtained that lim t e s = 0 to imply that the slow-loop system of HSV can track expected signals. In this way, Theorem 1 is proved and the argument projection is delivered as follows. □

3.2. Fault-Tolerant Control Design for Fast Loop of HSV

In this part, for the challenges caused by unexpected centroid shift, as shown in (8)–(12), the fault-tolerant controller of the quick-loop is deduced from the three aspects: (1) uncertainties of inertial matrix (J); (2) system uncertainties and eccentric moments; (3) system input saturation.

3.2.1. Fault-Tolerant Controller without Exact Knowledge of Inertial Matrix under Nominal Condition

As for the crucial role of the inertial matrix J on the controller design, once it is unexpectedly moved, as variation ( Δ J ) can present great challenges to controllers that rely on the inverse matrix of the system’s inertial matrix due to the potential irreversibility of ( J * + Δ J ), as shown in (8), such as sliding mode control [13], robust adaptive control [27], etc. In order to cope with this question, an adaptive feedback fault-tolerant controller that does not need the exact knowledge of ( J * + Δ J ) 1 is developed in this work to handle this question.
Before the fault-tolerant controller design, the angular velocity error ω ˜ is introduced as
ω ˜ = ω ω c
where ω c is the desired angular velocity. Furthermore, combined with Assumption 4 and (7), taking the derivative of ω ˜ yields [23]
ω ˜ ˙ = J 1 ( ω ˜ + ω c ) × J ( ω ˜ + ω c ) + J 1 τ ω ˙ c J 1 f ( ρ ¯ ) + J 1 d f + J 1 Δ τ
Considering (41) and (42), by assuming that system uncertainty f ( ρ ¯ ) , external disturbance d f and eccentric moment Δ τ are all known, the nominal controller τ without the knowledge of J is developed to achieve that ω ˜ is bounded and the other signals are bounded. The details are conveyed in the following three steps.
  • Step 1, a linear operator L : R 3 R 3 × 6 acting on a = [ a 1 , a 2 , a 3 ] T is introduced by
    L ( a ) = a 1 0 0 0 a 3 a 2 0 a 2 0 a 3 0 a 1 0 0 a 3 a 2 a 1 0
    Then, based on (8), an intermediate vector Ξ = [ J x x , J y y , J z z , J y z , J x z , J x y ] T is introduced to obtain the following transformation.
    J a = L ( a ) Ξ
  • Step 2, based on (44), the derivative of ω ˜ seen in (42) can be rewritten as
    ω ˜ ˙ = J 1 F ( ω ˜ , ω c , ω ˙ c ) Ξ + J 1 τ J 1 f ( ρ ¯ ) + J 1 d f + J 1 Δ τ
    where F : R 3 × R 3 × R 3 R 3 × 6 is defined by
    F ( ω ˜ , ω c , ω ˙ c ) = ( ω ˜ + ω c ) × L ( ω ˜ + ω c ) + L ( ω ˙ c )
  • Step 3, combined with Assumption 2, an adaptive feedback control law and its adaptive law are synthesized as
    τ = K ω ˜ + F ( ω ˜ , ω c , ω ˙ c ) Ξ ^ + f ( ρ ¯ ) d f Δ τ , Ξ ^ ˙ = Q ( F T ( ω ˜ , ω c , ω ˙ c ) ω ˜ k Ξ ^ Ξ ^ )
    where Q R 6 × 6 is a designed as positive definite symmetric matric, and k Ξ ^ represents the designed parameter, satisfying k Ξ ^ > 0 .

3.2.2. Adaptive Fault-Tolerant Controller without Exact Knowledge of Inertial Matrix

Due to the unknown effects of system uncertainty f ( ρ ¯ ) , external disturbance d f and eccentric moment Δ τ deriving from unexpected centroid shifts, actuator failures and external disturbances have a huge impact on the control performance of the designed controller (47). In order to overcome this obstacle, based on the characteristics of system uncertainty (11), eccentric moment (9) and LOE fault of actuator (12), the following adaptive laws are designed as
  • System uncertainty. By means of the Jacobi identity of the cross product, a × ( b × c ) + b × ( c × a ) + c × ( a × b ) = 0 , where a , b are the vectors. In view of (11), we have
    f ( ρ ¯ ) = m ρ ¯ × v ˙ A + m ω × ( ρ ¯ × v A ) + m v A × ( ω × ρ ¯ ) ρ ¯ × W = m ρ ¯ × v ˙ A m ρ ¯ × ( v A × ω ) ρ ¯ × W = m ρ ¯ × Φ
    where Φ = v ˙ A v A × ω W m , remarked as Φ = [ ϕ 1 ϕ 2 ϕ 3 ] T . From this it follows that
    ρ ¯ × Φ = ϕ 2 Δ z ϕ 3 Δ y ϕ 3 Δ x ϕ 1 Δ z ϕ 1 Δ y ϕ 2 Δ x = 0 ϕ 3 ϕ 2 ϕ 3 0 ϕ 1 ϕ 2 ϕ 1 0 Δ x Δ y Δ z = Φ × ρ ¯
    Subsequently, its adaptation law was developed as
    θ ^ ˙ 1 = Ψ 1 ( Φ × T ω ˜ k θ 1 ^ θ ^ 1 )
    where θ ^ 1 is a designed adaptive variable to estimate the unexpected effects of system uncertainties caused by centroid shift. Further, Ψ 1 R 3 × 3 represents a designed positive definite symmetric parameter matrix, and k θ 1 ^ stands for a constant coefficient, satisfying k θ 1 ^ > 0 .
  • Eccentric moment. For clarity, the effect of the eccentricity moment is rewritten as follows.
    Δ τ = 0 Z a Y a Z a 0 X a Y a X a 0 Δ x Δ y Δ z = Θ × ρ ¯
    and the corresponding adaptation law is developed as
    θ ^ ˙ 2 = Ψ 2 ( Θ × T ω ˜ k θ 2 ^ θ ^ 2 )
    where θ ^ 2 is also an adaptive state to obtain the effect of the eccentricity moment. Ψ 2 R 3 × 3 represents a design with positive definite symmetric coefficient matrix, and k θ 2 ^ is also a positive parameter.
  • LOE fault of actuator and external disturbance. For actuator faults and external perturbations d s , the complex perturbation D f is defined as
    D f = d f λ u τ
    where λ stands for the failure factor of the actuator fault, as shown in (12). Furthermore, combined with Assumption 3, the bounds of D f are introduced as κ d , satisfying D f κ d , where κ d > 0 . Furthermore, an adaptive estimator is introduced as
    κ ^ ˙ d = P ( ω ˜ T tanh ( ω ˜ ω ˜ ς 1 ς 1 ) k κ d κ ^ d )
    where κ ^ d is also utilized to estimate the influence of disturbance and LOE fault of actuator. k κ d and ς 1 are both positive constants.
In conclusion, combined with (50), (52) and (54), the nominal fault-tolerant controller of quick-loop (47) can be further developed as
Ξ ^ ˙ = Q ( F T ( ω ˜ , ν , ν ˙ ) ω ˜ k Ξ ^ Ξ ^ ) , θ 1 ^ ˙ = Ψ 1 ( Φ × T ω ˜ k θ 1 ^ θ 1 ^ ) , θ 2 ^ ˙ = Ψ 2 ( Θ × T ω ˜ k θ 2 ^ θ 2 ^ ) , κ ˙ d = P ( ω ˜ T tanh ( ω ˜ ω ˜ ς 1 ς 1 ) k κ d κ ^ d ) u τ = K ω ˜ + F ( ω ˜ , ν , ν ˙ ) Ξ ^ Φ × θ 1 ^ κ d ^ tanh ( ω ˜ ω ˜ ς 1 ς 1 ) Θ × θ 2 ^
Eventually, based on the Lyapunov stability theorem, a theoretical analysis is performed to summarize the main results of this section and guarantee the stability of the closed-loop control system.
Theorem 2.
Considering the quick-loop dynamics of a hypersonic vehicle (7) suffered from the variation of inertial matrix (8), eccentric moment (9), system uncertainty (11), actuator failures (12) and external disturbances, d f , adaptive feedback controls (55) without resorting to the exact knowledge of inertial matrix J by means of the adaptive laws (50), (52) and (54) developed to ensure that the quick-loop system of HSV can track expected signals, converging to a tiny compact set of the origin, and all the closed signals are bounded.
Proof. 
First, let us choose a positive definite Lyapunov function candidate as
V f = 1 2 ω ˜ T J ω ˜ + 1 2 Ξ ˜ T Q 1 Ξ ˜ + 1 2 θ 1 T Ψ 1 1 θ ˜ 1 + 1 2 θ 2 T Ψ 2 1 θ ˜ 2 + 1 2 κ ˜ d T P 1 κ ˜ d
where the estimate errors of Ξ ˜ , θ 1 ˜ , θ 2 ˜ and κ d ˜ are defined as Ξ ˜ = Ξ Ξ ^ , θ 1 ˜ = θ 1 θ 1 ^ , θ 2 ˜ = θ 2 θ 2 ^ , κ d ˜ = κ d κ d ^ .
Then, invoking Assumption 2 and (55), differentiating (56) gives
V ˙ f = ω ˜ T J ( J 1 ( ω ˜ + ν ) × J ( ω ˜ + ν ) ν ˙ + J 1 u τ + J 1 Φ × ρ ¯ + J 1 D f + J 1 Θ × ρ ¯ ) + Ξ ˜ T Q 1 Ξ ˜ ˙ + θ 1 ˜ T Ψ 1 1 θ 1 ˜ ˙ + θ 2 ˜ T Ψ 2 1 θ 2 ˜ ˙ + κ d ˜ T P 1 κ d ˜ ˙ = ω ˜ T J ( J 1 F ( ω ˜ , ν , ν ˙ ) Ξ + J 1 Φ × ρ ¯ + J 1 Θ × ρ ¯ + J 1 D f + J 1 ( K ω ˜ κ d ^ tanh ( ω ˜ ω ˜ ς 1 ς 1 ) + F ( ω ˜ , ν , ν ˙ ) Ξ ^ Φ × θ 1 ^ Θ × θ 2 ^ ) + Ξ ˜ T Q 1 Ξ ˜ ˙ + θ 1 ˜ T Ψ 1 1 θ 1 ˜ ˙ + θ 2 ˜ T Ψ 2 1 θ 2 ˜ ˙ + κ d ˜ T P 1 κ d ˜ ˙
where ω ˜ T ( D f κ d ^ tanh ( ω ˜ ς 1 ) ) 3 κ d ν ς 1 + κ d ˜ ω ˜ T tanh ( ω ˜ ς 1 ) , ν = 0.2785 [34].
To continue, submitting the adaptive laws (50), (52) and (54) into (57) can further give
V ˙ f ω ˜ T K ω ˜ ω ˜ T F ( ω ˜ , ν , ν ˙ ) Ξ ˜ + ω ˜ T Φ × ρ ¯ ˜ + ω ˜ T Θ × ρ ¯ ˜ + ω ˜ T κ ˜ d Ξ ˜ T Q 1 Ξ ^ ˙ θ 1 ˜ T Ψ 1 1 θ 1 ^ ˙ θ 2 ˜ T Ψ 2 1 θ 2 ^ ˙ κ d ˜ T P 1 κ d ^ ˙ ω ˜ T K ω ˜ + Ξ ˜ T k Ξ ^ Ξ ^ + θ ˜ 1 T k θ ^ 1 θ ^ 1 + θ ˜ 2 T k θ ^ 2 θ ^ 2 + κ ˜ d T k κ ^ d κ ^ d + 3 κ d ν ς 1
Subsequently, by taking the following inequalities
k Ξ ^ Ξ ˜ T Ξ ^ k Ξ ^ 2 Ξ ˜ F 2 + k Ξ ^ 2 Ξ F 2 , k θ ^ 1 θ ˜ 1 T θ ^ 1 k θ ^ 1 2 θ ˜ 1 F 2 + k θ ^ 1 2 θ 1 F 2 k θ ^ 2 θ ˜ 2 T θ ^ 2 k θ ^ 2 2 θ ˜ 2 F 2 + k θ ^ 2 2 θ 2 F 2 , k κ ^ d κ ˜ d T κ ^ d k κ ^ d 2 κ ˜ d F 2 + k κ ^ d 2 κ d F 2
into consideration, (58) can be further deduced as
V ˙ f ω ˜ T K ω ˜ k Ξ ^ 2 Ξ ˜ F 2 i = 1 2 k θ ^ i 2 θ ˜ i F 2 k κ ^ d 2 κ ˜ d F 2 + k Ξ ^ 2 Ξ F 2 + i = 1 2 k θ ^ i 2 θ i F 2 + k κ ^ d 2 κ d F 2 + 3 κ d ν ς 1
Furthermore, two indirect designed parameters are adopted as
H = { K min λ max ( J ) , k Ξ ^ min λ max ( Q 1 ) , k θ ^ 1 min λ max ( Ψ 1 1 ) , k θ ^ 2 min λ max ( Ψ 2 1 ) , k κ ^ d min λ max ( p 1 ) } M = k Ξ ^ 2 Ξ F 2 + i = 1 2 k θ ^ i 2 θ i F 2 + k κ ^ d 2 κ d F 2 + 3 κ d ν ς 1
Finally, (60) is reshaped as
V ˙ f ( t ) H V ( t ) + M
By multiplying e H t into both sides of (62) and integrating at the interval [ t 0 , t ) , we have
ω ˜ T ( t ) J ω ˜ ( t ) V f ( t ) V f ( t 0 ) e H ( t t 0 ) + M H ( 1 e H ( t t 0 ) )
where t 0 is the initial actuation time and t > t 0 .
Furthermore, by defining that V f ( t 0 ) e H ( t t 0 ) + M H ( 1 e H ( t t 0 ) ) = N ( t ) , combined with the definition of V f in (56), the boundedness of attitude angle rate error can be obtained as
λ min ( J ) ω ˜ ( t ) 2 2 ω ˜ T ( t ) J ω ˜ ( t ) N ( t )
where λ min ( J ) is the minimum eigenvalue of inertia matrix J. Moreover, note that N ( t ) M M H H as t . As a result, ω ˜ i ( t ) N ( t ) N ( t ) λ min ( J ) λ min ( J ) , i = 1 , 2 , 3 .
In addition, to guarantee the boundedness of the system parameters, the adaptive update law is modified by the parameter projection algorithm. The compact sets Π Ξ ^ , Π θ 1 , Π θ 2 and Π d ^ f are defined to restrain the above adaptive laws and the details are shown as follows [19]
Π Ξ = { Ξ ^ : Ξ ̲ Ξ ^ Ξ ¯ } , Π θ 1 = { θ 1 ^ : θ 1 ̲ θ 1 ^ θ 1 ¯ } , Π θ 2 = { θ 2 ^ : θ 2 ̲ θ 2 ^ θ 2 ¯ } , Π κ d = { κ d ^ : κ d ̲ κ d ^ κ d ¯ }
where Ξ ̲ , θ 1 ̲ , θ 2 ̲ , κ d ̲ and Ξ ¯ , θ 1 ¯ , θ 2 ¯ , κ d ¯ are the lower and upper bounds of adaptive laws. Then, the adaptive laws of Ξ ^ ˙ , θ 1 ^ ˙ , θ 2 ^ ˙ and κ d ^ ˙ are modified as follows
Ξ ^ ˙ = Proj ( Ξ ^ , Q ( F T ( ω ˜ , ν , ν ˙ ) ω ˜ k Ξ ^ Ξ ^ ) ) , θ 1 ^ ˙ = Proj ( θ 1 ^ , Ψ 1 ( Φ × T ω ˜ k θ 1 ^ θ 1 ^ ) , θ 2 ^ ˙ = Proj ( θ 2 ^ , Ψ 2 ( Θ × T ω ˜ k θ 2 ^ θ 2 ^ ) ) , κ d ^ ˙ = Proj ( κ d ^ , P ( ω ˜ T tanh ( ω ˜ ω ˜ ς 1 ς 1 ) k κ d κ ^ d ) )
where
Pr oj ( Ξ ^ , Ξ ^ ˙ ) = 0 0 Q ( F T ( ω ˜ , ν , ν ˙ ) ω ˜ k Ξ ^ Ξ ^ ) i f i f Ξ ^ = Ξ ¯ , Ξ ^ ˙ > 0 Ξ ^ = Ξ ̲ , Ξ ^ ˙ < 0 Pr oj ( θ i ^ , θ i ^ ˙ ) = 0 0 Ψ i ( Φ × T ω ˜ k θ i ^ θ i ^ ) i f i f θ i ^ = θ i ¯ , θ i ^ ˙ > 0 θ i ^ = θ i ̲ , θ i ^ ˙ < 0 Pr oj ( κ d ^ , κ d ^ ˙ ) = 0 0 P ( ω ˜ T tanh ( ω ˜ ω ˜ ς 1 ς 1 ) k κ d κ ^ d ) i f i f κ d ^ = κ d ¯ , κ d ^ ˙ > 0 κ d ^ = κ d ̲ , κ d ^ ˙ < 0
Finally, since by assumption ω c and ω ˙ c are bounded, it can be obtained that F ( ω ˜ , ω c , ω ˙ c ) is bounded as well. Since the boundedness of Ξ ˜ , θ 1 ˜ , θ 2 ˜ , κ d ˜ is already guaranteed by the use of the parametric projection operator, the stability of the closed-loop system is ensured and all the signals of the HSV fast-loop are bounded. □

4. Simulation

Considering that the HSV is carrying out a hypersonic re-entry flight with the speed of 3   km / s , there is an altitude of 30 km and a mass of m = 136,820 kg. The initial altitude vector is Ω = [ 0.017 , 0.026 , 0.017 ] rad and the angular rate is ω = [ 0 , 0 , 0 ] T rad / s . The reference signals are set as α c = 0.087 rad during 1 s~8 s, α c = 0.02 rad during 8 s∼20 s and α c = 0.05 rad during t > 20 s, besides, β is always 0 rad and μ is keeping 0.035   rad . In addition, the loss of efficiency fault of the actuator was set to λ = d i a g { 0.2 , 0.17 , 0.2 } as the 5th second.
The external perturbations of low-loop and quick-loop systems of HSV are described as follows d s 1 = 0.005 · sin t + 1 rad / s , d s 2 = 0.003 · cos t + 1 rad / s , d s 3 = 0.002 · sin t + 1 rad / s , d f 1 = 2 × 10 5 sin 4 t + 0.2 N · m , d f 2 = 4 × 10 5 sin 11 t 0.6 N · m , d f 3 = 4 × 10 5 sin 5 t + 0.2 N · m . The design parameters of the low-loop controller of HSV are as follows K s = diag { k 1 , k 2 , k 3 } , k i = 3 / ( 2 T ) , i = 1 , 2 , 3 , T = 2.4 . P s = diag { 5 , 5 , 5 } , θ = 25 diag { 3 , 3 , 3 } , ζ = 16 diag { 2 , 2 , 2 } , l 1 = 3.5 , l 2 = 2 , p = 5 , q = 3 . As for the fast-loop controller, K f = diag { 5 , 6.5 , 6.5 } , Q = 50 diag { 2.5 , 2.5 , 2.5 , 3 , 3 , 6 } , ψ 1 = diag { 2.5 , 2.5 , 2.5 } , ψ 2 = diag { 3.2 , 3.2 , 3.2 } , P = diag { 4.5 , 4.5 , 4.5 } .
To study the effect of the centroid shift on the HSV motion and validate the effectiveness of the designed fault-tolerant controller, simulations were performed. For analyzing the effects of unexpected centroid shift on the longitudinal movement of HSV, compared with the tracking effects of the FTC strategy consisted of sliding mode control combined with nonlinear disturbance observer and the FTC controller based on RBFNN [10], the simulation results demonstrate the effectiveness of the designed FTC strategy.
As shown in Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11, ( · ) n m l is the expected index of HSV’s slow and fast loops. ( · ) o b represents the fault-tolerant control effect of sliding mode controller (SMC) combined with the nonlinear disturbance observer (NDO) [35] by moving the uncertainties of J * on the right hands of attitude dynamics, and ( · ) n n represents the simulation results of a control scheme that is composed of nonlinear generalized predictive control (NGPC) and radial basis function neural network (RBFNN), as shown in [10]. In addition, ( · ) a d denotes the control effect of the controller designed in this paper. On the basis of Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8, in case of actuator failure only at the fifth second, attitude angles ( α , β , μ ) with ( · ) o b and ( · ) a d tracking effect, ( · ) n n has faster convergence speed. The main reason may be the heavy line load calculation of RBFNN and the horizontal regression principle of NGPC. Before the centroid shift occurs, it is straightforward to verify that the above controller presents an acceptable FTC tracking effect based on the arguments above.
When the unforeseeable centroid displacement occurs in the 10 t h second, the eccentric movement and variation of inertial matrix combined with the LOE fault of the actuator give rise to the fluctuation of the control surface of ( · ) o b compared with that of ( · ) n n and ( · ) a d , as shown in Figure 9, Figure 10 and Figure 11. However, this does not affect the steady flight of the HSV considerably for the time being, as shown by the trajectories in Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8 during the 10th and 20th seconds. Once the attitude angle index of HSV changes in the 20 t h second, found in Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8, the attitude angles and attitude angular rates both exhibit the phenomenon of rapid fluctuation in ( · ) o b compared with that of ( · ) a d . Besides, the tracking effects of ( · ) n n are of high tracking accuracy, but the relative low convergence speed, for instance, the tracking effect of α n n of Figure 3. It can be concluded that when the centroid shift and actuator fault both occur, the usual control strategy consisted of sliding mode controller and NDO can not work well under the variation of the inertial matrix J caused by unexpected centroid shift, which needs a more stable system control structure (exact knowledge of J and J 1 ), but the FTC strategies of ( · ) a d and ( · ) n n do not need.
Figure 9. Control surface deflection angles of δ e .
Figure 9. Control surface deflection angles of δ e .
Applsci 13 00830 g009
Figure 10. Control surface deflection angles of δ a .
Figure 10. Control surface deflection angles of δ a .
Applsci 13 00830 g010
Figure 11. Control surface deflection angles of δ r .
Figure 11. Control surface deflection angles of δ r .
Applsci 13 00830 g011
Furthermore, due to unknown effects being integrated into one item to facilitate the observation of NDO, the larger deflection angles of control surfaces are needed to cope with the estimated effects through NDO, as shown in the response curves of control surfaces ( · ) o b of Figure 9, Figure 10 and Figure 11 after the 20th second. Moreover, given the heavy online computational load of neural networks, the application of FTC controller resting on RBFNN is utilized to handle effects ( Δ J , Δ τ and f ( ρ ¯ ) ) one-by-one obtain excellent tracking effects, but this FTC scheme shows a low convergence speed, shown in the tracking effects of ( · ) n n in Figure 9, Figure 10 and Figure 11. To remedy this, the FTC algorithm of this work is developed by employing adaptive observers and estimators in conjunction with adaptive feedback control, which has a known and fixed control structure and does not need the exact knowledge of inertial matrix (J and J 1 ). Finally, the simulation results show that the proposed FTC strategy in this work can suppress unexpected effects and maintain excellent tracking performance.

5. Conclusions

In this paper, we develop a useful solution to the HSV problem with unexpected centroid shifts, actuator faults and external perturbations. Resorting to a designed sixth-order dynamic compensator, the proposed adaptive FTC scheme can achieve the attitude tracking objective without the exact knowledge of the inertia matrix J, despite the presence of variations in the system inertia matrix, eccentricity and actuator faults. Subsequently, based on the properties of eccentric moments, systematic uncertainties and external perturbations, the performance of the proposed FTC controller can be improved by adapting the corresponding adaptive estimator to support the designed FTC controller using adaptive techniques.Then, the stability of this designed FTC strategy is guaranteed by the Lyapunov stability theorem. Finally, the simulation results show that the proposed FTC scheme is effective for the reconfigurable HSV flight control design under the condition of the inertia matrix and eccentricity.
In future work, we will focus on the FTC of HSV with structural failures caused by variable centroid displacement (the structure and size of the elements of the inertial matrix J are uncertain). Furthermore, the problem of deflection saturation of the control surface should also be taken into consideration [24].

Author Contributions

Writing—original draft preparation, H.Y.; writing—review and editing, Y.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (61773012), the National natural sciences fund youth fund project (61807017), the research start-up fund of Jiangsu University of science and technology (1052932101).

Institutional Review Board Statement

This study did not involve humans or animals.

Informed Consent Statement

This study did not involve humans or animals.

Data Availability Statement

Some or all data, models, or code that support the findings of this study are available from the corresponding author upon reasonable request. The database is identified by References [8,10].

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A

c 21 = ( q 2 + r 2 ) cos α sin β + ( p q r ˙ ) cos β ( p r + q ˙ ) sin α sin β c 22 = ( p 2 + r 2 ) cos β ( q p + r ˙ ) cos α sin β ( p ˙ q r ) sin α sin β c 23 = ( p ˙ q r ) cos β ( q ˙ p r ) cos α sin β ( p 2 + q 2 ) sin α sin β c 31 = ( p r + q ˙ ) cos α ( q 2 + r 2 ) sin α c 32 = ( p ˙ q r ) cos α ( q p + r ˙ ) sin α c 33 = ( p 2 + q 2 ) cos α ( q ˙ r p ) sin α
σ 1 σ 2 σ 3 = 0 d γ d t 0 + cos γ 0 sin γ 0 1 0 sin γ 0 cos γ 0 0 d χ d t = d χ d t sin γ d γ d t d χ d t cos γ
ρ ¯ x p = Δ x cos α cos β + Δ y sin β + Δ z sin α cos β ρ ¯ y p = Δ x ( cos α sin β cos μ + sin α sin μ ) + Δ y cos β cos μ + Δ z ( cos μ sin α sin β cos α sin μ ) ρ ¯ z p = Δ x ( cos μ sin α cos α sin β sin μ ) + Δ y cos β sin μ + Δ z ( cos α cos μ sin α sin β sin μ )
U A V A W A = V cos α cos β V sin β V sin α cos β , U ˙ A V ˙ A W ˙ A = S α , β T V ˙ V β ˙ V cos β α ˙
S α , β T = cos α cos β cos α sin β sin α sin β cos β 0 sin α cos β sin α sin β cos α .
J x x = J x x * + m ( Δ y 2 + Δ z 2 ) , J y y = J y y * + m ( Δ x 2 + Δ z 2 ) J z z = J z z * + m ( Δ x 2 + Δ y 2 ) , J x y = m Δ x Δ y J x z = m Δ x Δ z , J y z = m Δ y Δ z
Y ¯ ( t ) = y ^ 0 y ^ [ 1 ] y ^ [ ρ + r ] = h ( x ) L f s h ( x ) L f s ρ + r h ( x ) + 0 m × 1 0 m × 1 H ˜ ( u ^ z )
where
H ˜ ( u ^ z ) u ^ z = L g s L f s ρ 1 h ( x ) 0 m × m 0 m × m 0 m × m * m × m L g s L f s ρ 1 h ( x ) 0 m × m 0 m × m * m × m * m × m * m × m * m × m L g s L f s ρ 1 h ( x )

Appendix B

L , D and Y are described as follows [33]
L = q ^ S ( Δ C L + C L * ) , D = q ^ S ( Δ C D + C D * ) , Y = q ^ S ( Δ C Y + C Y * )
where q ^ = ρ V 2 ρ V 2 2 2 . C L * = C L , α + C L , δ e δ e + C L , δ a δ a , C D = C D , α + C D , δ e δ e + C D , δ a δ a + C D , δ r δ r , and C Y = C Y , β β + C Y , δ e δ e + C Y , δ a δ a + C Y , δ r δ r . δ e , δ a and δ r are the control deflections of the aileron, elevator and rudder, respectively.
The rolling, pitching and yawing moments are described as [28]
l A = q ^ b S ( C l * + Δ C l ) + Δ l A , m A = q ^ c S ( Δ C m + C m * ) X c g Z + Δ m A n A = q ^ b S ( Δ C n + C n * ) X c g Y + Δ n A
where C l * = C l , β β + C l , δ e δ e + C l , δ a δ a + C l , δ r δ r + C l , p p · b 2 V + C l , r r · b 2 V , C m = C m , a + C m , δ e δ e + C m , δ a δ a + C m , δ r δ r + C m , q q · c 2 V , C n = C n , β β + C n , δ e δ e + C n , δ a δ a + C n , δ r δ r + C n , p p · b 2 V + C n , r r · b 2 V . Z = D sin α L cos α , X c g = 1.6469 e 10 * M 2 5.5715 e 05 * M + 7.3672 .

References

  1. Tang, T.; Qi, R.Y.; Jiang, B. Adaptive nonlinear generalized predictive control for hypersonic vehicle with unknown parameters and control constraints. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2017, 233, 510–532. [Google Scholar] [CrossRef]
  2. Yu, X.; Li, P.; Zhang, Y.M. The design of fixed-Time observer and finite-time fault-tolerant control for hypersonic gliding vehicles. IEEE Trans. Ind. Electron. 2018, 65, 4135–4144. [Google Scholar] [CrossRef] [Green Version]
  3. Nguyen, N.; Krishnakumar, K.; Kaneshige, J.; Nespeca, P. Dynamics and adaptive control for stability recovery of damaged asymmetric aircraft. In Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, Keystone, CO, USA, 21 August 2006; p. 6049. [Google Scholar]
  4. Li, J.; Gao, C.; Jing, W.; Wei, P. Dynamic analysis and control of novel moving mass flight vehicle. Acta Astronaut. 2017, 131, 36–44. [Google Scholar] [CrossRef]
  5. Li, J.Q.; Gao, C.S.; Li, C.Y. A survey on moving mass control technology. Aerosp. Sci. Technol. 2018, 82, 594–606. [Google Scholar] [CrossRef]
  6. Gasbarri, P.; Sabatini, M.; Pisculli, A. Dynamic modelling and stability parametric analysis of a flexible spacecraft with fuel slosh. Acta Astronaut. 2016, 127, 141–159. [Google Scholar] [CrossRef]
  7. Bacon, B.J.; Gregory, I.M. General equations of motion for a damaged asymmetric aircraft. AIAA Paper 2007-6306. In Proceedings of the 2007 AIAA Guidance, Navigation, and Control Conference, Hilton Head, SC, USA, 20 August 2007. [Google Scholar]
  8. Nguyen, N.; Krishnakumar, K.; Kaneshige, J.; Nespeca, P. Flight dynamics modeling and hybrid adaptive control of damaged asymmetric aircraft. J. Guid. Control. Dyn. 2008, 31, 171–186. [Google Scholar]
  9. Liu, Y.; Tao, G.; Joshi, S.M. Modeling and model reference adaptive control of aircraft with asymmetric damage. J. Guid. Control. Dyn. 2012, 33, 1500–1517. [Google Scholar] [CrossRef]
  10. Meng, Y.Z.; Jiang, B.; Qi, R.Y. Adaptive non-singular fault-tolerant control for hypersonic vehicle with unexpected centroid shift. IET Control Theory Appl. 2019, 13, 1773–1785. [Google Scholar] [CrossRef]
  11. Lee, T. Robust adaptive attitude tracking on SO(3) with an application to a quadrotor UAV. J. Guid. Control. Dyn. 2013, 21, 1924–1930. [Google Scholar]
  12. Yin, S.; Ding, S.X.; Xie, X.; Luo, H. A review on basic data-driven approaches for industrial process monitoring. IEEE Trans. Ind. Electron. 2014, 61, 6418–6428. [Google Scholar] [CrossRef]
  13. Xiao, B.; Yin, S.; Hao, H. Reconfigurable tolerant control of uncertain mechanical systems with actuator faults: A sliding mode observer-based approach. IEEE Trans. Control Syst. Technol. 2018, 26, 1249–1258. [Google Scholar] [CrossRef]
  14. Bounemeur, A.; Chemachema, M.; Essounbouli, N. Indirect adaptive fuzzy fault-tolerant tracking control for MIMO nonlinear systems with actuator and sensor failures. ISA Trans. 2018, 79, 45–61. [Google Scholar] [CrossRef] [PubMed]
  15. Hu, Q.; Niu, G.; Wang, C. Spacecraft attitude fault-tolerant control based on iterative learning observer and control allocation. Aerosp. Sci. Technol. 2018, 75, 245–253. [Google Scholar] [CrossRef]
  16. Xu, D.Z.; Jiang, B.; Shi, P. Robust NSV fault-tolerant control system design against actuator faults and control surface damage under actuator dynamics. IEEE Trans. Ind. Electron. 2015, 62, 5919–5928. [Google Scholar] [CrossRef]
  17. Jiang, B.; Gao, Z.; Shi, P.; Xu, Y. Adaptive fault-tolerant tracking control of near-space vehicle using takagi-sugeno fuzzy models. IEEE Trans. Fuzzy Syst. 2010, 18, 1000–1007. [Google Scholar] [CrossRef]
  18. Yin, S.; Yang, H.; Gao, H.; Qiu, J.; Kaynak, O. An adaptive NN-based approach for fault-tolerant control of nonlinear time-varying delay systems with unmodeled dynamics. IEEE Trans. Neural Netw. Learn. Syst. 2017, 28, 1902–1913. [Google Scholar] [CrossRef]
  19. He, J.; Qi, R.; Jiang, B.; Qian, J. Adaptive output feedback fault-tolerant control design for hypersonic flight vehicles. J. Frankl. Inst. 2015, 352, 1811–1835. [Google Scholar] [CrossRef]
  20. Zhang, A.; Hu, Q.L.; Friswell, M.I. Finite-time fault tolerant attitude control for over-activated spacecraft subject to actuator misalignment and faults. Iet Control Theory Appl. 2013, 7, 2007–2020. [Google Scholar] [CrossRef] [Green Version]
  21. Chen, M.; Yu, J. Disturbance observer-based adaptive sliding mode control for near-space vehicles. Nonlinear Dyn. 2015, 82, 1671–1682. [Google Scholar] [CrossRef]
  22. Chen, W.H.; Ballance, D.J.; Gawthrop, P.J. Brief optimal control of nonlinear systems: A predictive control approach. Automatica 2003, 39, 633–641. [Google Scholar] [CrossRef] [Green Version]
  23. Chaturvedi, N.A.; Bernstein, D.S.; Ahmed, J.; Bacconi, F.; McClamroch, N.H. Globally convergent adaptive tracking of angular velocity and inertia identification for a 3-DOF rigid body. IEEE Trans. Control. Syst. Technol. 2006, 14, 841–853. [Google Scholar] [CrossRef] [Green Version]
  24. Hu, Q.; Shao, X.; Guo, L. Adaptive fault-tolerant attitude tracking control of spacecraft with prescribed performance. IEEE/ASME Trans. Mechatron. 2018, 23, 331–341. [Google Scholar] [CrossRef]
  25. Meng, Y.Z.; Jiang, B.; Qi, R.Y. Modeling and control of hypersonic vehicle dynamic under centroid shift. Adv. Mech. Eng. 2018, 10, 1–21. [Google Scholar] [CrossRef]
  26. Zong, Q.; Zeng, F.L.; Zhang, B.X.; You, M. Modeling and Model Verification of Hypersonic Vehicle; Science China Press: Beijing, China, 2016. [Google Scholar]
  27. Yang, H.; Jiang, B.; Liu, H.H.; Yang, H.; Zhang, Q. Attitude synchronization for multiple 3-DOF helicopters with actuator faults. IEEE/ASME Trans. Mechatron. 2019, 24, 597–608. [Google Scholar] [CrossRef]
  28. Shaughnessy, J.D.; Pinckney, S.Z.; McMinn, J.D. Hypersonic Vehicle Simulation Model: Winged-Cone Configuration; NASA TM-102610; NASA Langley Reserch Center: Hamption, VA, USA, 1990; pp. 1–140.
  29. Xu, B.; Qi, R.; Jiang, B. Adaptive fault-tolerant control for HSV with unknown control direction. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 2743–2758. [Google Scholar] [CrossRef]
  30. Liu, Y.H.; Huang, L.; Xiao, D. Adaptive dynamic surface control for uncertain nonaffine nonlinear systems. Int. J. Robust Nonlinear Control 2017, 27, 535–546. [Google Scholar] [CrossRef]
  31. Xian, B.; Zhang, Y. Continuous asymptotically tracking control for a class of nonaffine-in-input system with non-vanishing disturbance. IEEE Trans. Autom. Control 2017, 62, 6019–6025. [Google Scholar] [CrossRef]
  32. Levant, A. Higher-order sliding modes, differentiation and output-feedback control. Int. J. Control 2003, 76, 924–941. [Google Scholar] [CrossRef]
  33. Du, Y.L. Study of Nonlinear Adaptive Attitude and Trajectory Control for Near Space Vehicles; Nanjing University of Aeronautics and Astronautics: Nanjing, China, 2010. [Google Scholar]
  34. Lai, G.; Liu, Z.; Zhang, Y.; Chen, C.P. Adaptive position/attitude tracking control of aerial robot with unknown inertial matrix based on a new robust neural identifier. IEEE Trans. Neural Netw. Learn. Syst. 2016, 27, 18–31. [Google Scholar] [CrossRef]
  35. Lee, D. Nonlinear disturbance observer-based robust control for spacecraft formation flying. Aerosp. Sci. Technol. 2018, 76, 82–90. [Google Scholar] [CrossRef]
Figure 1. Three-dimensional view of the displacement of the center of mass, O b O b .
Figure 1. Three-dimensional view of the displacement of the center of mass, O b O b .
Applsci 13 00830 g001
Figure 2. Frame of fault-tolerant controller for HSV with centroid shift.
Figure 2. Frame of fault-tolerant controller for HSV with centroid shift.
Applsci 13 00830 g002
Figure 3. Tracking effects of attack angle α .
Figure 3. Tracking effects of attack angle α .
Applsci 13 00830 g003
Figure 4. Tracking effects of bank angle μ .
Figure 4. Tracking effects of bank angle μ .
Applsci 13 00830 g004
Figure 5. Tracking effects of sideslip angle β .
Figure 5. Tracking effects of sideslip angle β .
Applsci 13 00830 g005
Figure 6. Tracking effects of roll rate p.
Figure 6. Tracking effects of roll rate p.
Applsci 13 00830 g006
Figure 7. Tracking effects of pitch rate q.
Figure 7. Tracking effects of pitch rate q.
Applsci 13 00830 g007
Figure 8. Tracking effects of yaw rate r.
Figure 8. Tracking effects of yaw rate r.
Applsci 13 00830 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ye, H.; Meng, Y. Adaptive Fault-Tolerant Control of Hypersonic Vehicles with Unknown Model Inertia Matrix and System Induced by Centroid Shift. Appl. Sci. 2023, 13, 830. https://0-doi-org.brum.beds.ac.uk/10.3390/app13020830

AMA Style

Ye H, Meng Y. Adaptive Fault-Tolerant Control of Hypersonic Vehicles with Unknown Model Inertia Matrix and System Induced by Centroid Shift. Applied Sciences. 2023; 13(2):830. https://0-doi-org.brum.beds.ac.uk/10.3390/app13020830

Chicago/Turabian Style

Ye, Hui, and Yizhen Meng. 2023. "Adaptive Fault-Tolerant Control of Hypersonic Vehicles with Unknown Model Inertia Matrix and System Induced by Centroid Shift" Applied Sciences 13, no. 2: 830. https://0-doi-org.brum.beds.ac.uk/10.3390/app13020830

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