Next Article in Journal
Auto Sweptback Wing Based on Low Scattering Demand for an Unmanned Aerial Vehicle in Phase Flight
Next Article in Special Issue
Inverse Optimal Zero Effort Miss Guidance Based on Disturbance Observer
Previous Article in Journal
Playful Probing: Towards Understanding the Interaction with Machine Learning in the Design of Maintenance Planning Tools
Previous Article in Special Issue
Longitudinal Aerodynamic Parameter Identification for Blended-Wing-Body Aircraft Based on a Wind Tunnel Virtual Flight Test
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sliding Mode Backstepping Control for the Ascent Phase of Near-Space Hypersonic Vehicle Based on a Novel Triple Power Reaching Law

1
College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 211106, China
2
College of Electronic Information and Automation Engineering, Guilin University of Aerospace Technology, Guilin 541004, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(12), 755; https://doi.org/10.3390/aerospace9120755
Submission received: 14 October 2022 / Revised: 21 November 2022 / Accepted: 24 November 2022 / Published: 26 November 2022
(This article belongs to the Special Issue Flight Control)

Abstract

:
This paper presents a novel sliding mode backstepping control scheme for the ascent phase of a near-space hypersonic vehicle (NSHV) base on a triple power reaching law (TPRL). A new model transformation is proposed for NSHV with uncertain parameters subject to uncertainties during ascent phase. To shorten the reaching time and reduce the chattering in sliding mode scheme, TPRL is proposed. Then, based on TPRL, a sliding mode backstepping control scheme is proposed, which is combined with new adaptive laws to further reduce the adverse impact of uncertainties. Simulation results demonstrate that TPRL is effective, and the proposed controller for the ascent phase of NSHV is robust with respect to uncertainties.

1. Introduction

A near-space hypersonic vehicle (NSHV) refers to an aircraft operating in near-space, the flight speed of which can exceed Mach-5 [1,2,3]. Great concerns have been raised regarding the development of NSHV, due to its significant military value and potential economic value [4,5,6]. However, NSHV is of high couplings, strong nonlinearities, and complex uncertainties. Especially in ascent phase, the altitude and speed of NSHV change dramatically, which results in a large consumption of fuel, and the aerodynamic characteristics of NSHV vary greatly [7,8,9]. Therefore, it is a great challenge to design a high-performance flight controller for NSHV in the ascent phase.
The controller design for the ascent phase of NSHV encounters challenges originating from model uncertainties and nonlinearities. Therefore, a controller with robustness and stability is especially important for the ascent phase of NSHV. In recent years, there have been many methods used to study NSHV in the literature, such as H∞ [10], predictive flight control [11,12], fuzzy control [13,14], sliding mode control [15,16], backstepping control [17,18], neural network control [19,20], and adaptive control [21,22].
The sliding mode control provides a systematic approach to modeling imprecision, of which the main advantage is that it is robust subject to model uncertainties and disturbances [23,24]. Yin et al. proposed a high-order sliding mode control scheme for hypersonic vehicles, which combine the advantages of the disturbance observer and gain adaptation to estimate the uncertainties and external disturbance [25]. Guo et al., in 2018, proposed adaptive twisting sliding mode control for hypersonic re-entry vehicles, which adopted an adaptive finite-time observer to estimate the unknown state, guaranteed the stability of system and reduced the chattering phenomenon [26]. In the literature [27], Guo et al. estimated the uncertainty through the observer and proposed a fixed time sliding mode scheme, which reduced the interference of the uncertainty and obtained good tracking effect. However, the impact of external disturbances on the vehicle needs to be further explored. Basin et al. proposed a higher order sliding mode observer using fixed-settling time differentiators, which guarantee the robustness and high-accuracy output tracking in the presence of external disturbances and missile model uncertainties [28]. In the literature [29], Wang et al. proposed a terminal sliding mode controller based on sliding mode disturbance observer for a hypersonic flight vehicles. This approach obtains expected control performance in the presence of uncertainties and external disturbances.
The principle of sliding mode control is that the state of the system approaches the sliding surface under reaching law, so the stability and convergence time of the system are directly affected by the reaching law [30,31]. Many approaches have been proposed to eliminate or attenuate chattering and accelerate convergence time. Zhang et al. developed exponential reaching law in sliding mode control to enable the hypersonic vehicle to track the command input and reach a steady state along the designed sliding surface in the case of mismatched uncertainty [32]. Mozayan et al. proposed a sliding mode controller using a new enhanced exponential reaching law, which can mitigate chattering and approach the sliding surface quickly [33]. In the literature [34], fast terminal sliding-mode control was proposed to ensure that the system states converge within a finite time. An improved double power reaching law was proposed in the literature [35], which can shorten the time to approach the equilibrium point when the system is away from the sliding surface.
As mentioned above, the technique of sliding mode control has good properties, such as its robustness, order reduction, insensitivity to parameter variations, finite-time convergence, disturbance rejection, and good dynamic behavior. However, pure sliding mode control presents drawbacks that include large control authority requirements and chattering [24,33,36]. The performance of pure sliding mode control can be improved by coupling with adaptive control scheme to estimate the parameters of uncertainties and disturbances. In the literature [28], Basin et al. designed a non-recursive higher order sliding mode observer, which was combined with an adaptive algorithm to estimate the control gain to reduce control chattering. Jiao et al. proposed a new adaptive sliding mode control scheme which was combined with type-2 fuzzy approach to make the controller show good robustness in the presence of uncertain parameters. However, this control method was only applied to the aircraft in the cruise phase [37]. Ding et al. in 2020 proposed an improved continuous sliding mode control scheme to ensure the convergence precision, in which nonsingular fast fixed-time is adopted and dual-layer adaptive continuous twisting reaching law is developed [38]. Sagliano et al. proposed an adaptive disturbance-based sliding-mode controller for hypersonic-entry vehicles, which was able to accurately reconstruct disturbances and has good tracking performance [39].
The backstepping approach demonstrates great performance in complicated nonlinear systems control, such as NSHV in the presence of uncertainties and disturbances. In the current years, in order to improve the performance of flight controller, many advanced control methods based on backstepping have been proposed. Zhang et al. proposed an anti-disturbance backstepping control approach for hypersonic vehicle in presence of the large uncertainties and the external disturbances, which obtain the superior tracking performance [40]. In the literature [41], adaptive backstepping controller and dynamic inverse controller were developed for hypersonic vehicle subject to input constraints and aerodynamic uncertainties. Based on backstepping control, Wu et al. developed a nonlinear disturbance observer for a flexible air-breathing hypersonic vehicle [42]. Yu et al. proposed a novel adaptive backstepping control method for hypersonic vehicle with mismatched uncertainties to guarantee asymptotical command tracking [43]. Lu proposed a disturbance observer-based backstepping control scheme for hypersonic flight vehicles, where a nonlinear disturbance observer is conducted to estimate the model uncertainties [44]. However, the literature mentioned above has not focused on the design of hypersonic vehicles in the ascent phase.
In this paper, our scope is to propose a backstepping-based adaptive sliding mode controller adopting triple power reaching law (TPRL), which is able to meet the requirements in tracking control for the ascent phase of NSHV. The contributions of this paper are as follows:
(1)
The proposed TPRL shortens the reaching time and reduces the chattering of sliding mode approach. To shorten the reaching time of system converging to the equilibrium point, the reaching process can be divided by TPRL into three stages, and the reaching speed of each stage can be set, respectively. Then, TPRL can accelerate the speed of convergence effectively. Meanwhile, TPRL has the capability of mitigating the chattering of the sliding mode approach with respect to the traditional reaching law.
(2)
A novel tracking controller for the ascent phase of NSHV is proposed. At present, due to the limitation of data, little research has focused on the tracking control for the ascent phase of hypersonic vehicle in the presence of uncertainties. However, as we know, the ascent mission of NSHV could significantly affect the whole flight performance of the NSHV. In order to improve the tracking accuracy and robustness in the ascent mission of NSHV, a novel backstepping-based adaptive sliding mode controller is developed. Firstly, a transformation model is proposed for NSHV in face of uncertainties. The strict feedback form with uncertain parameters is adopted in the transformed mode, which is more suitable to accommodate uncertainties on NSHV than traditional model transformation. Furthermore, based on the proposed model transformation, a backstepping control scheme is proposed, which is combined with sliding mode control using TPRL. To further attenuate the adverse impact of uncertainties, new adaptive laws are designed. The flight controller proposed improves the tracking performance for NSHV in ascent phase.
This paper is organized as follows. The linearized model of the NSHV is established and a model transformation is proposed in Section 2. In Section 3, a novel TPRL is proposed, and a new adaptive slide mode controller based on backstepping is designed. In Section 4, based on the proposed method, a controller of the NSHV is developed. Numerical simulation results are reported in Section 5. Finally, the summary of key features of the proposed scheme are given in Section 6.

2. Model of NSHV

2.1. Longitudinal Model of NSHV

The three-view drawing of NSHV studied in this paper is shown in Figure 1.
Using the longitudinal force and moment equilibrium of hypersonic vehicle, the longitudinal model of NSHV is obtained as follows [45,46]:
    V ˙ = T cos α D m μ r 2 sin γ γ ˙ = L + T sin α m V μ V 2 r V r 2 cos γ α ˙ = q γ ˙ q ˙ = M y y I y y h ˙ = V sin γ m ˙ = T g I s p
where V , h , m , γ , α , and q are velocity, altitude, mass, flight path angle, attack angle, and pitch angle rate, respectively. The μ , g , and r are the Earth’s gravity constant, the acceleration of gravity and radial distance from the Earth’s center. The M y y , I y y   L , D , T , and I s p are pitch moment, rotation inertia, lift, drag, and thrust and fuel ratio, respectively. The M y y , L , D , and T can be expressed as follows:
    M y y = 1 2 ρ V 2 s c ¯ C M α + C M q + C M δ e L = 1 2 ρ V 2 s C L D = 1 2 ρ V 2 s C D T = 1 2 ρ V 2 s C T
where, ρ , s , and c ¯ are air density, wing area, and mean aerodynamic chord, respectively; C M α , C M q , and C M δ e are moment coefficient due to attack angle, moment coefficient due to pitch rate, and moment coefficient due to elevator deflection, respectively. C L , C D , and C T are lift coefficient, drag coefficient, and engine thrust coefficient, respectively.
In the ascent phase, the thrust coefficient and fuel ratio are set as follows:
    C T = 0.12886 β                                                             ( β < 1 ) 0.112 + 0.0168 β                                   β 1 I s p = 4700 h 100                                                     ( M a < 4 ) 0.112 + 0.0168 β                                   M a 4
where β and M a are the state of engine and Mach, respectively. The engine dynamics are modelled as follows:
β ¨ = 2 ξ ω β ˙ ω 2 β + ω 2 β c
where β ,  ω , ξ and β c are the state of engine, natural frequency, damp ratio, and throttle setting, respectively.
The linearized model of NSHV is developed by repeatedly differentiating V three times and h four times. The output dynamics variables of V and h can be expressed in a form in which control input β c and δ e appear explicitly as follows [32]:
    V = f V + b 11 β c + b 12 δ e h 4 = f h + b 21 β c + b 22 δ e
where the detailed expression of the f V , f h , b 11 , b 12 , b 21 , and b 22 can be found in [32].
In ascent phase, uncertainties are modelled as additive variance Δ to the nominal value, which is expressed as follows:
i = i 0 + Δ i
where i = m ,   I y y ,   ρ ,   s ,   c ¯ ,   c e ,   C L ,   C D ,   C T ,   C M α ,   C M q , C M δ e , i is the real value, i 0 is the nominal value, and Δ i represents the value of uncertainty.

2.2. Novel Model Transformation

According to Equation (5), traditional model transformation can be expressed as follows:
    ξ ˙ 1 = ξ 2 ξ ˙ 2 = ξ 3 ξ ˙ 3 = f V + b 11 β c + b 12 δ e η ˙ 1 = η 2 η ˙ 2 = η 3 η ˙ 3 = η 4 η ˙ 4 = f h + b 21 β c + b 22 δ e
where ξ = ξ 1 , ξ 2 , ξ 3 T = V , V ˙ , V ¨ T and η = η 1 , η 2 , η 3 , η 4 T = h , h ˙ , h ¨ ,   h T . It is obvious that the Model (7) is not suitable for the system with uncertainties.
To deal with the drawback of the traditional model (7), we define new states as x 1 = 0 t ( V τ d τ h T , x 2 = V h ˙ T , x 3 = V ˙   h ¨ T , and x 4 = V ¨   h T . Then, the model of NSHV in Equation (5) is written as the following novel model transformation.
    x ˙ 1 = x 2 + φ 1 T x 1 θ 1 x ˙ 2 = x 3 + φ 2 T x 2 θ 2 x ˙ 3 = x 4 + φ 3 T x 3 θ 3 x ˙ 4 = f x , t + G x , t u + φ 4 T x 4 θ 4
where
    f x , t = f V f h T G x , t = b 11 b 12 b 21 b 22 u = β c δ e T
where φ i T x i θ i ,   i = 1 , 2 , 3 , 4 represent the uncertainties of each subsystem; φ i T x i = d i a g φ i , 1 , φ i , 2 are known functions, which are assumed to be sufficient smooth; and θ i = θ i , 1 , θ i , 2 T are unknown constant parameters.

3. Sliding Mode Backstepping Controller Design

In this section, a triple power reaching law is proposed and a sliding mode backstepping tracking control scheme combined with adaptive laws is developed. The system structure of the proposed controller is shown in Figure 2.

3.1. A Novel Triple Power Reaching Law for Sliding Mode

In order to speed up the convergence and reduce the chattering of the sliding mode controller, a novel triple power reaching law (TPRL) is proposed and given by:
s ˙ = l 1 s α s g n s l 2 s β s g n s l 3 s γ s g n s
where s is the sliding surface; s g n s is the sign function; l 1 > 0 , l 2 > 0 , l 3 > 0 , α > 1 , 0 < β < 1 ; and γ is a variable parameter and is given by:
γ =     δ ,         s > ω 1 ,         s ω
where δ and ω are constants and δ > α , ω > 1 .
Theorem 1. 
For the reaching law given as Equation (10), the state of system s converges to the equilibrium point in fixed time and the convergence time T satisfies the following conditions: T < t 1 + t 2 + t 3 , where
    t 1 = ω 1 γ s 0 1 γ γ 1 l 3 t 2 = 1 α 1 α α 1 l 1 t 3 = 1 1 β l 2
Proof of Theorem 1. 
(1) Analysis of accessibility and stability
According to Equation (10), we obtain:
s s ˙ = s l 1 s α s g n s l 2 s β s g n s l 3 s γ s g n s = l 1 s α + 1 l 2 s β + 1 l 3 s γ + 1 0
Only if s = 0 can reach s s ˙ = 0 .
Therefore, the system state s can reach the equilibrium point s = 0 under the reaching law Equation (10).
(2) Fixed-time convergence
It is supposed that the initial state of the system is s 0 > ω > 1 . The process of system convergence may be divided into 3 stages as follows: s 0 s ( t 1 ) = ω s ( t 2 ) = 1 s ( t 3 ) = 0 . For the sake of analysis, we assume that l 1 = l 2 = l 3 .
(i)
s 0 s ( t 1 ) = ω
In this stage, because γ = ω , ω > α > β , then s γ > s α > s β , the third term l 3 s γ s g n s plays a major role in the reaching law. In other words, the reaching speed at this stage is mainly affected by the third term; then, the reaching law Equation (10) can be expressed as:
s ˙ = l 3 s γ s g n s
Taking the integrate of Equation (14), we have:
s 1 γ = 1 γ l 3 t + s 0 1 γ
Then, the max reaching time in this stage can be expressed as follows:
t 1 = ω 1 γ s 0 1 γ γ 1 l 3
Thus, the time required for the system to go from s 0 to s ( t 1 ) is less than t 1 .
(ii)
s ( t 1 ) = ω s ( t 2 ) = 1
In this stage, γ = 1 , α > γ > β > 0   and s α > s γ > s β , and the first term l 1 s α s g n s plays a major role in the reaching law; then, the reaching law Equation (10) can be expressed as:
s ˙ = l 1 s α s g n s
Integrating Equation (17), we obtain:
s 1 α = 1 γ l 3 t + α 1 α
Then, the max reaching time in this stage can be expressed as follows:
t 2 = 1 α 1 α α 1 l 1
Therefore, the time required for the system to go from s ( t 1 ) to s ( t 2 ) is less than t 2 .
(iii)
s ( t 2 ) = 1 s ( t 3 ) = 0
In this stage, because s < 0 , s β > s γ > s α , the second term l 2 s β s g n s plays a major role in the reaching law, and the reaching law Equation (10) can be expressed as:
s ˙ = l 2 s β s g n s
Integrating Equation (20), we have:
s 1 β = 1 β l 2 t + 1
Then, the max reaching time in this stage can be expressed as follows:
t 3 = 1 1 β l 2
Therefore, the time required for the system to go from s ( t 2 ) to s ( t 3 ) is less than t 3 .
Above all, for the reaching law of sliding mode control as Equation (10), and the reaching time T satisfies the following formula:
T < t 1 + t 2 + t 3
The proof is completed. □
According to the analysis above, for the reaching law given by Equation (10), when the initial state of the system is s 0 > ω , the convergence process of sliding surface s can be divided into 3 stages, s 0 s ( t 1 ) = ω s ( t 2 ) = 1 s ( t 3 ) = 0 , and TPRL can speed up the reaching time of each stage. The advantages of TPRL can be conclude as follows:
(1) Compared with traditional reaching law, the TPRL proposed in this paper can shorten the reaching time of system converging to the equilibrium point. Firstly, TPRL can reduce the reaching time when the initial conditions of the system parameters are far from the equilibrium point. According to Equation (10), when the sliding surface   s > ω , γ = δ > α . Compared with the traditional double power reaching law, which has the form s ˙ = l 1 s α s g n s l 2 s β s g n s , the third term l 3 s γ s g n s of TPRL is the largest, which plays a major role in the reaching law and speed up the process of convergence. Secondly, TPRL can reduce the reaching time when the condition of the system is close to the equilibrium point. Compared with the traditional double power reaching law, when   s ω , γ = 1 , TPRL can reach the equilibrium point in shorter time, since the third term of TPRL is set as l 3 s s g n s in this stage, which can increase the value of the reaching law, and thus, shorten the reaching time of the system.
(2) TPRL can attenuate the chattering of the sliding mode approach more effectively than traditional reaching law. In order to reduce chattering, this paper tries to present a new reaching law TPRL, which is continuous when s 0 , lim s 0 ( l 1 s α s g n s l 2 s β s g n s l 3 s γ s g n s ) = 0 ,   s 0 + 0 ,   s 0 , where s is the sliding surface. For the traditional exponential reaching law given as s ˙ = l 1 s g n s l 2 s , the discontinuous functions s g n s may cause some undesired effects, such as the chattering phenomenon due to lim s 0 [ l 1 s g n s l 2 s ] = l 1 ,   s 0 + l 1 ,   s 0 .

3.2. Sliding Mode Backstepping Controller Design

This subsection presents a sliding mode backstepping control scheme for a class of high-order nonlinear system in the presence of uncertainties. The aim is to design a controller to guarantee that the output of the system will track the command signal.
A high order nonlinear system with uncertainties is expressed as follows:
    x ˙ 1 = x 2 + φ 1 T x 1 θ 1 x ˙ 2 = x 3 + φ 2 T x 2 θ 2                           x ˙ i = x i + 1 + φ i T x i θ i                           x ˙ n = f x , t + G x , t u + φ n T x n θ n y = x 1
where x i = x i , 1 , x i , 2 , , x i , m T m , i = 1 , 2 n are states; f x , t and G x , t are given nonlinear functions; u m is control input; y m is the system output; φ i T x i = d i a g φ i , 1 , φ i , 2 , , φ i , m , i = 1 , 2 n are known functions, which are assumed to be sufficient smooth; and θ i = θ i , 1 , θ i , 2 , , θ i , m T m , i = 1 , 2 n are unknown constant parameters, where θ i is defined as follows:
θ i = θ ^ i + θ ˜ i
where θ ^ i is the estimation of θ i , θ ˜ i is estimation error.
The design steps for controller are as follows.
Defining a command signal for the first subsystem as μ 1 = x 1 d , the tracking error for the first subsystem is defined as follows:
z 1 = x 1 μ 1
Taking the derivative of Equation (26), we obtain:
z ˙ 1 = x ˙ 1 μ ˙ 1 = x 2 + φ 1 T x 1 θ 1 μ ˙ 1
To stabilize the first subsystem, the virtual control law is defined as follows:
μ 2 = k 1 z 1 + μ ˙ 1 φ 1 T x 1 θ ^ 1
where k 1 is a positive definite matrix, k i , i = 2 , 3 n , in later steps are also positive definite matrixes. The tracking error for the second subsystem is defined as follows:
z 2 = x 2 μ 2
Substituting Equations (28) and (29) into Equation (27), we get:
z ˙ 1 = k 1 z 1 + z 2 φ 1 T x 1 θ ^ 1 + φ 1 T x 1 θ 1 = k 1 z 1 + z 2 + φ 1 T x 1 θ ˜ 1
A Lyapunov function candidate for the first subsystem is chosen as follows:
V 1 = 1 2 z 1 T z 1 + 1 2 θ ˜ 1 T Γ 1 1 θ ˜ 1
where Γ 1 is a symmetric positive definite matrix, Γ i , i = 2 , 3 n , in later steps are also symmetric positive definite matrixes. Taking the derivative of Equation (31) along Equation (30), we have:
V ˙ 1 = z 1 T k 1 z 1 + z 2 + φ 1 T x 1 θ ˜ 1 + θ ˜ 1 T Γ 1 1 θ ˙ 1 θ ^ 1 ˙ = k 1 z 1 T z 1 + z 1 T z 2 + z 1 T φ 1 T x 1 θ ˜ 1 + θ ˜ 1 T Γ 1 1 θ ˙ 1 θ ^ 1 ˙ = k 1 z 1 T z 1 + z 1 T z 2 + θ ˜ 1 T φ 1 x 1 z 1 Γ 1 1 θ ^ 1 ˙
A Lyapunov function candidate for the 2nd subsystem is chosen as follows:
V 2 = V 1 + 1 2 z 2 T z 2 + 1 2 θ ˜ 2 T Γ 2 1 θ ˜ 2
Taking the derivative of Equation (33), we obtain:
V ˙ 2 = V ˙ 1 + z 2 T k 2 z 2 z 1 + z 3 + φ 2 T x 2 θ ˜ 2 + θ ˜ 2 T Γ 2 1 θ ˙ 2 θ ^ 2 ˙ = V ˙ 1 + z 2 T k 2 z 2 z 1 + z 3 + z 2 T φ 2 T x 2 θ ˜ 2 + θ ˜ 2 T Γ 2 1 θ ˙ 2 θ ^ 2 ˙ = k 1 z 1 T z 1 k 2 z 2 T z 2 + z 2 T z 3 + θ ˜ 1 T φ 1 x 1 z 1 Γ 1 1 θ ^ 1 ˙ + θ ˜ 2 T φ 2 x 2 z 2 Γ 2 1 θ ^ 2 ˙ = i = 1 2 [ k i z i T z i + θ ˜ i T φ i x i z i Γ i 1 θ ^ i ˙ ] + z 2 T z 3
where z 3 = x 3 μ 3 is the tracking error for the 2nd subsystem and μ 3 = k 2 z 2 + μ ˙ 2 z 1 φ 2 T x 2 θ ^ 2 is the virtual control law.
Similarly, a Lyapunov function candidate for the n 1 t h subsystem is chosen as:
V n 1 = V n 2 + 1 2 z n 1 T z n 1 + 1 2 θ ˜ n 1 T Γ n 1 1 θ ˜ n 1
Taking the derivative of Equation (35), we obtain:
V ˙ n 1 = V ˙ n 2 + z n 1 T k n 1 z n 1 z n 2 + z n 1 + θ ˜ n 1 T ( z n 1 Γ n 1 1 θ ^ ˙ n 1 ) = k 1 z 1 T z 1 k 2 z 2 T z 2 k n 1 z n 1 T z n 1 + z n 1 T z n + θ ˜ 1 T φ 1 x 1 z 1 Γ 1 1 θ ^ 1 ˙ + θ ˜ 2 T φ 2 x 2 z 2 Γ 2 1 θ ^ 2 ˙ + + θ ˜ n 1 T ( φ n 1 x n 1 z n 1 Γ n 1 1 θ ^ ˙ n 1 ) = i = 1 n 1 [ k i z i T z i + θ ˜ i T ( φ i x i z i Γ i 1 θ ^ ˙ i ) ] + z n 1 T z n  
where z n = x n μ n is the tracking error for the nth subsystem and μ n = k n 1 z n 1 + μ ˙ n 1 z n 2 φ n T x n θ ^ n is the virtual control law.
Defining a sliding function for system as follows:
S = c 1 z 1 + c 2 z 2 + c i z i + + z n = i = 1 n 1 c i z i + z n
where S m , c i , i = 1 ,   2 , n 1 are positive definite matrixes.
Taking the derivative of Equation (37) along Equation (24), we obtain:
S ˙ = i = 1 n 1 c i z ˙ i + z ˙ n = i = 1 n 1 c i z ˙ i + x ˙ n μ ˙ n = i = 1 n 1 c i z ˙ i + f x , t + G x , t u + φ n T x n θ n μ ˙ n
A Lyapunov function candidate for the nth subsystem is chosen as:
V n = V n 1 + 1 2 S T S + 1 2 θ ˜ n T Γ n 1 θ ˜ n = 1 2 i = 1 n 1 k i z i T z i + 1 2 S T S + 1 2 j = 1 n θ ˜ j T Γ j 1 θ ˜ j
Taking the derivative of Equation (39) along Equations (36) and (38), we obtain:
V ˙ n = V ˙ n 1 + S T S ˙ + θ ˜ n T Γ n 1 θ ~ ˙ n = V ˙ n 1 + S T [ i = 1 n 1 c i z ˙ i + f ( x , t ) + G ( x , t ) u + φ n T x n θ n μ ˙ n ] + θ ˜ n T Γ n 1 θ ˙ n θ ˜ ˙ n = V ˙ n 1 + S T [ i = 1 n 1 c i z ˙ i + f ( x , t ) + G ( x , t ) u + φ n T x n θ ^ n μ ˙ n ] + θ ˜ n T φ n x n S Γ n 1 θ ^ ˙ n = i = 1 n 1 [ k i z i T z i + θ ˜ i T φ i x i z i Γ i 1 θ ^ ˙ i ] + z n 1 T z n + S T [ i = 1 n 1 c i z ˙ i + f ( x , t ) + G ( x , t ) u + φ n T x n θ ^ n μ ˙ n ] + θ ˜ n T [ φ n x n S Γ n 1 θ ^ n ]
Substituting Equation (37) into Equation (40), we have:
V ˙ n = i = 1 n 1 [ k i z i T z i + θ ˜ i T ( φ i x i z i Γ i 1 θ ^ ˙ i ) ] + z n 1 T S i = 1 n 1 c i z i + S T [ i = 1 n 1 c i z ˙ i + f x , t + G x , t u + φ n T x n θ ^ n μ ˙ n ] + θ ˜ n T [ φ n x n S Γ n 1 θ ^ n ˙ ] = i = 1 n 1 [ k i z i T z i + θ ˜ i T ( φ i x i z i Γ i 1 θ ^ ˙ i ) ] z n 1 T i = 1 n 1 c i z i + S T [ z n 1 + i = 1 n 1 c i z ˙ i + f x , t + G x , t u + φ n T x n θ ^ n μ ˙ n ] + θ ˜ n T [ φ n x n S Γ n 1 θ ^ n ˙ ] = i = 1 n 1 k i z i T z i z n 1 T i = 1 n 1 c i z i + S T [ z n 1 + i = 1 n 1 c i z ˙ i + f x , t + G x , t u + φ n T x n θ ^ n μ ˙ n ] + i = 1 n 1 θ ˜ i T [ φ i x i z i Γ i 1 θ ^ ˙ i ] + θ ˜ n T [ φ n x n S Γ n 1 θ ^ n ˙ ]
In order to ensure the stability of system, the control law is designed as follows:
u = G x , t 1 ( z n 1 i = 1 n 1 c i z ˙ i f x , t G x , t φ n T x n θ ^ n + μ ˙ n + u s w )
where u s w m is the triple power reaching law proposed as Equation (10) and u s w is designed as:
u s w , i = l i , 1 s i α i s g n s i l i , 2 s i β i s g n s 1 l i , 1 s i γ i s g n s i
where i = 1 ,   2 m . Adaptive laws are designed as follows:
    θ ^ ˙ 1 = Γ 1 φ 1 x 1 z 1 θ ^ ˙ 2 = Γ 2 φ 2 x 2 z 2         θ ^ ˙ n 1 = Γ n 1 φ n 1 x n 1 z n 1 θ ^ ˙ n = Γ n φ n x n S

3.3. Stability Analysis

Theorem 2. 
The closed-loop system consisting of the nonlinear system given by Equation (24) with the controller given by Equation (42), and the adaptive laws given by Equation (44) can be guaranteed to be stable.
Proof of Theorem 2. 
A Lyapunov function candidate is chosen as Equation (39). Substituting Equations (42) and (44) into Equation (41), we obtain:
V ˙ n = j = 1 n 1 k j z j T z j z n 1 T i = 1 n 1 c i z i S T u s w = j = 1 n 1 k j z j T z j z n 1 T i = 1 n 1 c i z i i = 1 m ( l i , 1 s i α i + 1 + l i , 2 s i β i + 1 + l i , 3 s i γ i + 1 ) = z T Q z i = 1 m ( l i , 1 s i α i + 1 + l i , 2 s i β i + 1 + l i , 3 s i γ i + 1 )
in which
Q = k 1 0 0 0 0 k 2 0 0 0 0 k n 2 0 c 1 c 2 c n 2 k n 1 + c n 1
and
z = z 1 z 2 z n 1 T
Because the principal minors in matrix Q are positive, Q is a positive definite matrix. Then, we obtain:
V ˙ n = z T Q z i = 1 m ( l i , 1 s i α i + 1 + l i , 2 s i β i + 1 + l i , 3 s i γ i + 1 ) 0
The proof is completed. □

4. Flight Controller Design for NSHV

In this section, the flight controller for the ascent phase of NSHV is developed through Theorems 1 and 2.

4.1. Controller Design

For the given transformation model of Equation (8), the command signals are velocity V d and altitude h d and the control object is that the controller can asymptotically track the command signal.
According to Equations (42) and (43), the tracking control law is designed as follows:
u = G x , t 1 ( z 3 i = 1 3 c i z ˙ i f x , t φ 3 T x 3 θ ^ 3 + μ ˙ 4 + u s w )
The reaching laws for sliding mode is designed as follows:
u s w = l 1 , 1 s 1 α 1 s g n s 1 l 1 , 2 s 1 β 1 s g n s 1 l 1 , 3 s 1 γ 1 s g n s 1 l 2 , 1 s 2 α 2 s g n s 2 l 2 , 2 s 2 β 2 s g n s 2 l 2 , 3 s 2 γ 2 s g n s 2
where
γ i =     δ i ,         s i > ω i 1 ,           s i ω i , i = 1 , 2
According to Equation (44), adaptive laws are designed as follows:
    θ ^ ˙ 1 = Γ 1 φ 1 x 1 z 1 θ ^ ˙ 2 = Γ 2 φ 2 x 2 z 2 θ ^ ˙ 3 = Γ 3 φ 3 x 3 z 3 θ ^ ˙ 4 = Γ 4 φ 4 x 4 S
where S = c 1 z 1 + c 2 z 2 + c 3 z 3 + z 4 , and the errors are defined as follows:
    z 1 = x 1 μ 1 z 2 = x 2 μ 2 z 3 = x 3 μ 3 z 4 = x 4 μ 4
The virtual control laws are designed as follows:
    μ 1 = 0 t V d τ d τ h d T μ 2 = k 1 z 1 + μ ˙ 1 φ 1 T x 1 θ ^ 1 μ 3 = k 2 z 2 + μ ˙ 2 z 1 φ 2 T x 2 θ ^ 2 μ 4 = k 3 z 3 + μ ˙ 3 z 2 φ 3 T x 3 θ ^ 3

4.2. Stability Analysis

To analyze the stability of system, according to Equation (39), a Lyapunov function candidate for NSHV is chosen as:
V 4 = 1 2 i = 1 3 k i z i T z i + 1 2 j = 1 4 θ ˜ j T Γ j 1 θ ˜ j + 1 2 S T S
According to Theorem 2, two matrixes are defined as follows:
    z = z 1 z 2 z 3 T Q = k 1 0 0 0 k 2 0 c 1 c 2 k 3 + c 3
Taking the derivative of Equation (55), according to Equation (48), we have:
V ˙ 4 = z T Q z i = 1 2 ( l i , 1 s i α i + 1 + l i , 2 s i β i + 1 + l i , 3 s i γ i + 1 ) 0

5. Simulation Results

For the verification of proposed control method in NSHV applications, simulation is performed by MATLAB/Simulink. This section includes two scenarios, and the NSHV is flying in the ascent phase. Scenario 1 verifies the advantage of the novel reaching law proposed in this paper, and the good tracking performance of the proposed method in this paper for NSHV in the ascent phase is illustrated in Scenario 2.

5.1. Scenario 1: Simulation of Reaching Laws

To illustrate the effectiveness of triple power reaching law proposed in this paper, five different reaching laws are adopted for comparison. The five compared reaching law are expressed as follows:
(a) Triple power reaching law (TPRL):
s ˙ = l 1 s α s g n s l 2 s β s g n s l 3 s γ s g n s
(b) Double power reaching law (DPRL):
s ˙ = l 1 s α s g n s l 2 s β s g n s
(c) Exponential reaching law (ERL):
s ˙ = l 1 s g n s l 2 s
(d) Traditional symbolic function reaching law (TRL):
s ˙ = l 1 s g n s
(e) Fast power reaching law (FPRL):
s ˙ = l 1 s β s g n s l 2 s
The traditional sliding mode controller (SMC) is adopted in this scenario to track the velocity of the NSHV, and the sliding function is defined as follows [24]:
s v = d d t + λ 3 0 t e τ d t
The model parameters of the NSHV are presented in Table 1.
In this scenario, the initial parameters of the NSHV are set as V 0 = 4590   m / s , h 0 = 33 , 528   m , γ = 0 ° , α = 2.75 ° , and q = 0 ° / s . the uncertainties are set as 0, and the command signal of step velocity is set as 10 m/s.
The value of parameters for five different reaching laws are given in Table 2.
The simulation results for five different reaching laws are shown in Figure 3, Figure 4 and Figure 5.
As can be seen from Figure 3, in the five reaching laws, the sliding surface under TPRL can converge to the equilibrium point zero in the shortest time. In Figure 4, compared with the other four reaching laws, the rise-time of velocity under TPRL is shortest and the velocity tracking under TPRL reaches a steady state in a shortest time. In Figure 5, the control input under TPRL exhibits no chattering, while that for TRL and ERL exhibit chattering.
According to the simulation results above, the following conclusions can be drawn as follows:
(1) It can be known from Figure 3 and Figure 4 that TPRL proposed in this paper can guarantee system converging to the equilibrium point in shorter time than other reaching laws mentioned above. In other words, considering similar gains, TPRL approach specifies faster reaching speed with respect to other reaching laws.
(2) It can be seen from Figure 5 that TPRL reduces the chattering of sliding mode approach with respect to the traditional reaching laws. When the system state comes close to the equilibrium point, the coefficients of TPRL methods gradually diminish to mitigate the chattering problem, expressed as lim s 0 ( l 1 s α s g n s l 2 s β s g n s l 3 s γ s g n s ) = 0 , s 0 + 0 , s 0 . However, in Scenario 1, when s→0, ERL is expressed as lim s 0 [ l 1 s g n s l 2 s ] = 1 , s 0 + 1 , s 0 , and it is shown as lim s 0 [ l 1 s g n s ] = 1 , s 0 + 1 , s 0 for TRL. Therefore, TPRL can reduce chattering effectively than traditional reaching law ERL and TRL.

5.2. Scenario 2: Simulation of Controller for NSHV

In this scenario, the NSHV is flying in the ascent phase, and the initial conditions of the NSHV are V 0 = 1700   m / s , h 0 = 12 , 000   m , γ = 0 ° , α = 1.3 ° , and q = 0 ° / s and the final velocity and altitude are V f = 3700   m / s and h f = 32 , 000   m , respectivly. It is worth noting that due to limitations in real systems, the elevator deflection δ e should be limited in the range from 20 ° to 20 ° and β c should be limited in the range from 0 to 1. The parameters of the NSHV are given in Table 1.
To verify effectiveness of the method proposed in this paper, the command signals are set as h c t = h 0 + Δ h t and V c t = V 0 + Δ V t , where Δ h t and Δ V t are generated by inputting step signals into the following filters:
    Δ h s h s t e p s = 0.03 4 s + 0.03 4 Δ V s V s t e p s = 0.03 4 s + 0.03 4
where the altitude step signal is set as h step = 20 , 000   m and the velocity step signal is set as V step = 2000   m / s . Uncertainties and disturbances are considered as follows:
(1)
Parameters of uncertainties are given by Equation (13) and set as 20%, which are presented as follows:
Δ i = 0.2 s i n 0.02 π t
where i = m , I y y , ρ , s , c ¯ , c e , C L , C D , C T , C M α , C M q , C M δ e
(2)
Parameters of disturbances for control input are set as follows:
    β c = β c 0 + 0.1 s i n 0.02 π t δ e = δ e 0 + 0.1 s i n 0.02 π t
The controller is developed in Equation (49) along with the adaptive update laws given by Equation (52). The parameters of the controller are presented in Table 3.
To explore the effectiveness of proposed method, it is compared with three approaches: (a) the traditional sliding mode control using traditional double power reaching law (SMC); (b) the backstepping sliding mode control using traditional double power reaching law (BSMC), where BSMC is based on the conventional model given Equation (7); and (c) the terminal sliding mode control (TSMC). The simulation results of the NSHV in ascent phase are shown in Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14.
The tracking performance of velocity and altitude are shown in Figure 6, Figure 7, Figure 8 and Figure 9. As can be seen in Figure 6 and Figure 7, compared with SMC, BSMC, and TSMC, the proposed method can reduce the tracking error of velocity more effectively and converge in a shorter time. By contrast, there is a 7.4 m/s maximum tracking error under SMC, 3.5 m/s under BSMC, and 3.8 m/s under TSMC, respectively. The tracking errors of velocity for SMC, BSMC, and TSMC are with larger fluctuations. Furthermore, Figure 8 and Figure 9 show that the altitude tracking error of the proposed method in this paper is the smallest and the controller developed by the proposed method can track stably. In comparison, the tracking errors of altitude under SMC, BSMC, and TSMC are seen in larger fluctuations.
Figure 10, Figure 11 and Figure 12 show the responses of attack angle, pitch angle rate, and mass, respectively. From Figure 10 and Figure 11, it obvious that the responses of attack angle and pitch angle under the proposed method are smooth and within a reasonable range. The responses of mass in Figure 12 indicate that NSHV consumes a large amount of fuel in the ascent phase and the mass changes smoothly. The control input responses of the elevator deflection and the throttle setting are shown in Figure 13 and Figure 14. In Figure 13, the responses of the throttle setting change smoothly and fluctuate within the margin range. In Figure 14, it can be seen that the response of the elevator deflection under the method proposed is with smaller fluctuations than those for SMC, BSMC, and TSMC and fluctuates within the margin range.
According to the simulation results above, the following points can be summarized:
(1)
In the ascent phase of NSHV, sliding mode backstepping controller based on TPRL guarantees a better tracking performance in situations of uncertainties in comparison with SMC, BSMC, and TSMC. Firstly, the tracking errors of both velocity and altitude under the proposed method are smaller than those for SMC, BSMC, and TSMC. Secondly, the attack angle and pitch angle rate change more smoothly under the proposed method. Thirdly, the control inputs curve for throttle setting and elevator deflection are smoother under proposed method, compared with SMC, BSMC, and TSMC.
(2)
It is evident that the proposed controller can reduce the adverse influence of uncertainties more effectively in comparison with SMC, BSMC, and TSMC. BSMC is based on the model transformation given by Equation (7), which is not suited to dealing with uncertainties, while the proposed method is based on the new model transformation given by Equation (8), which can compensate the uncertainties more effectively. In the ascent phase of NSHV, uncertainties are the significant reason for tracking errors. It can be seen that the tracking errors of velocity and altitude converge more effectively under the proposed method, while those for SMC, BSMC, and TSMC are with larger fluctuations.

6. Conclusions

In this paper, a backstepping-based adaptive sliding mode controller using a novel reaching law TPRL is proposed for NSHV in the ascent phase. Firstly, a new model transformation for NSHV is proposed in the presence of uncertainties. Secondly, in order to reduce reaching time and attenuate chattering for sliding mode control, a new reaching law TPRL is proposed. Then, to ensure the tracking accuracy and robustness of system, based on TPRL, a sliding mode backstepping controller is developed. Furthermore, the adaptive laws are designed to reduce the adverse effect of uncertainties. Finally, the simulation results verify the advantages of TPRL and also demonstrate that the proposed controller provides good tracking performance for NSHV in the ascent phase and shows robustness in case of uncertainties.

Author Contributions

Conceptualization, S.H.; methodology, S.H. and J.J.; software, S.H.; validation, S.H. and O.L.; formal analysis, O.L.; investigation, S.H.; resources, S.H. and J.J.; data curation, S.H.; writing—original draft preparation, S.H.; writing—review and editing, S.H. and O.L. All authors have read and agreed to the published version of the manuscript. S.H. is a PHD candidate from Nanjing University of Aeronautics and Astronautics, and J.J. is supervisor of S.H.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 61966010 and 61673209.

Institutional Review Board Statement

Not applicable for studies not involving humans or animals.

Informed Consent Statement

Not applicable for studies not involving humans or animals.

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, Y.; Yang, X.; Yan, H. Reliable fuzzy tracking control of near-space hypersonic vehicle using aperiodic measurement information. IEEE Trans. Ind. Electron. 2019, 66, 9439–9447. [Google Scholar] [CrossRef]
  2. Li, Z.; Shi, S. L1 adaptive loss fault tolerance control of unmanned hypersonic aircraft with elasticity. Aerospace 2021, 8, 176. [Google Scholar] [CrossRef]
  3. Zhang, J.; Sun, C.; Zhang, R.; Qian, C. Adaptive sliding mode control for re-entry attitude of near space hypersonic vehicle based on backstepping design. IEEE/CAA J. Autom. Sin. 2015, 2, 94–101. [Google Scholar]
  4. Zhao, K.; Song, J.; Ai, S.; Xu, X.; Liu, Y. Active fault-tolerant control for near-space hypersonic vehicles. Aerospace 2022, 9, 237. [Google Scholar] [CrossRef]
  5. Wang, G.; An, H.; Wang, Y.; Xia, H.; Ma, G. Intelligent control of air-breathing hypersonic vehicles subject to path and angle-of-attack constraints. Acta Astronaut. 2022, 198, 606–616. [Google Scholar] [CrossRef]
  6. Ding, Y.; Yue, X.; Chen, G.; Si, J. Review of control and guidance technology on hypersonic vehicle. Chin. J. Aeronaut. 2022, 35, 1–18. [Google Scholar]
  7. Liu, Y.; Dong, C.; Zhang, W.; Wang, Q. Phase plane design based fast altitude tracking control for hypersonic flight vehicle with angle of attack constraint. Chin. J. Aeronaut. 2021, 34, 490–503. [Google Scholar] [CrossRef]
  8. Zhang, D.; Zheng, X.; Xie, Y.; Hu, X. Angular-Accelerometer-Based Flexible-State Estimation and Tracking Controller Design for Hypersonic Flight Vehicle. Aerospace 2022, 9, 206. [Google Scholar] [CrossRef]
  9. Zuo, R.; Li, Y.; Lv, M.; Liu, Z. Realization of trajectory precise tracking for hypersonic flight vehicles with prescribed performances. Aerosp. Sci. Technol. 2021, 111, 106554. [Google Scholar] [CrossRef]
  10. Wu, H.; Feng, S.; Liu, Z.; Guo, L. Disturbance observer based robust mixed H2/H∞ fuzzy tracking control for hypersonic vehicles. Fuzzy Sets Syst. 2017, 306, 118–136. [Google Scholar]
  11. Zhang, B.; Zhou, D. Optimal predictive sliding-mode guidance law for intercepting near-space hypersonic maneuvering target. Chin. J. Aeronaut. 2022, 35, 320–331. [Google Scholar] [CrossRef]
  12. Guo, J.; Peng, Q.; Zhou, J. Disturbance observer-based nonlinear model predictive control for air-breathing hypersonic vehicles. J. Aerosp. Eng. 2019, 32, 1–13. [Google Scholar] [CrossRef]
  13. Bu, X.; Qi, Q. Fuzzy optimal tracking control of hypersonic flight vehicles via single-network adaptive critic design. IEEE Trans. Fuzzy Syst. 2022, 30, 270–278. [Google Scholar] [CrossRef]
  14. Bu, X.; Xiao, Y.; Lei, H. An adaptive critic design-based fuzzy neural controller for hypersonic vehicles: Predefined behavioral nonaffine control. IEEE/ASME Trans. Mechatron. 2020, 24, 1871–1881. [Google Scholar] [CrossRef]
  15. Shou, Y.; Xu, B.; Liang, X.; Yang, D. Aerodynamic/reaction-jet compound control of hypersonic reentry vehicle using sliding mode control and neural learning. Aerosp. Sci. Technol. 2021, 111, 106564. [Google Scholar] [CrossRef]
  16. Sun, J.; Pu, Z.; Yi, J.; Liu, Z. Fixed-time control with uncertainty and measurement noise suppression for hypersonic vehicles via augmented sliding mode observers. IEEE Trans. Ind. Inform. 2020, 16, 1192–1203. [Google Scholar] [CrossRef]
  17. Chen, K.; Zhu, S.; Wei, C.; Xu, T.; Zhang, X. Output constrained adaptive neural control for generic hypersonic vehicles suffering from non-affine aerodynamic characteristics and stochastic disturbances. Aerosp. Sci. Technol. 2021, 111, 106469. [Google Scholar] [CrossRef]
  18. Sun, J.; Yi, J.; Pu, Z.; Liu, Z. Adaptive fuzzy nonsmooth backstepping output-feedback control for hypersonic vehicles with finite-time convergence. IEEE Trans. Fuzzy Syst. 2020, 28, 2320–2334. [Google Scholar] [CrossRef]
  19. Chai, R.; Tsourdos, A.; Savvaris, A.; Chai, S.; Xia, Y.; Chen, C.L.P. Six-DOF spacecraft optimal trajectory planning and real-time attitude control: A deep neural network-based approach. IEEE Trans. Neural Netw. Learn. Syst. 2020, 31, 5005–5013. [Google Scholar] [CrossRef] [PubMed]
  20. Xu, B.; Wang, D.; Zhang, Y.; Shi, Z. DOB-based neural control of flexible hypersonic flight vehicle considering wind effects. IEEE Trans. Ind. Electron. 2017, 64, 8676–8685. [Google Scholar] [CrossRef]
  21. An, H.; Wu, Q.; Wang, G.; Kao, Y.; Wang, C. Adaptive compound control of air-breathing hypersonic vehicles. IEEE Trans. Aerosp. Electron. Syst. 2020, 56, 4519–4532. [Google Scholar] [CrossRef]
  22. Dong, C.; Liu, Y.; Wang, Q. Barrier Lyapunov function based adaptive finite-time control for hypersonic flight vehicles with state constraints. ISA Trans. 2020, 96, 163–176. [Google Scholar] [CrossRef] [PubMed]
  23. Mu, C.; Zong, Q.; Tian, B.; Xu, W. Continuous sliding mode controller with disturbance observer for hypersonic vehicles. IEEE/CAA J. Autom. Sin. 2015, 2, 45–55. [Google Scholar]
  24. Xu, H.; Mirmirani, M.D.; Ioannou, P.A. Adaptive sliding mode control design for a hypersonic flight vehicle. J. Guid. Control. Dyn. 2004, 27, 829–838. [Google Scholar] [CrossRef]
  25. Yin, X.; Wang, B.; Liu, L.; Wang, Y. Disturbance observer-based gain adaptation high-order sliding mode control of hypersonic vehicles. Aerosp. Sci. Technol. 2019, 89, 19–30. [Google Scholar] [CrossRef]
  26. Guo, Z.; Chang, J.; Guo, J.; Zhou, J. Adaptive twisting sliding mode algorithm for hypersonic reentry vehicle attitude control based on finite-time observer. ISA Trans. 2018, 77, 20–29. [Google Scholar] [CrossRef]
  27. Guo, J.; Yang, S.; Guo, Z. Robust tracking for hypersonic vehicles subjected to mismatched uncertainties via fixed-time sliding mode control. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2021, 235, 2145–2153. [Google Scholar]
  28. Basin, M.V.; Yu, P.; Shtessel, Y.B. Hypersonic missile adaptive sliding mode control using finite and fixed-time observers. IEEE Trans. Ind. Electron. 2018, 65, 930–941. [Google Scholar]
  29. Wang, J.; Wu, Y.; Dong, X. Recursive terminal sliding mode control for hypersonic flight vehicle with sliding mode disturbance observer. Nonlinear Dyn. 2015, 81, 1489–1510. [Google Scholar] [CrossRef]
  30. Ma, H.; Wu, J.; Xiong, Z. A novel exponential reaching law of discrete-time sliding-mode control. IEEE Trans. Ind. Electron. 2017, 64, 3840–3850. [Google Scholar]
  31. Wang, Y.; Feng, Y.; Zhang, X.; Liang, J. A new reaching law for antidisturbance sliding-mode control of PMSM speed regulation system. IEEE Trans. Power Electron. 2020, 35, 4117–4126. [Google Scholar] [CrossRef]
  32. Zhang, Y.; Li, R.; Xue, T.; Lei, Z. Exponential sliding mode tracking control via back-stepping approach for a hypersonic vehicle with mismatched uncertainty. J. Frankl. Inst. 2016, 353, 2319–2343. [Google Scholar] [CrossRef]
  33. Mozayan, S.M.; Saad, M.; Vahedi, H.; Fortin-Blanchette, H.; Soltani, M. Sliding mode control of PMSG wind turbine based on enhanced exponential reaching law. IEEE Trans. Ind. Electron. 2016, 63, 6148–6159. [Google Scholar] [CrossRef]
  34. Xu, S.S.; Chen, C.; Wu, Z. Study of nonsingular fast terminal sliding-mode fault-tolerant control. IEEE Trans. Ind. Electron. 2015, 62, 3906–3913. [Google Scholar] [CrossRef]
  35. Pan, J.; Li, W.; Zhang, H. Control algorithms of magnetic suspension systems based on the improved double exponential reaching law of sliding mode control. Int. J. Control. Autom. Syst. 2018, 16, 2878–2887. [Google Scholar] [CrossRef]
  36. Rabiee, H.; Ataei, M.; Ekramian, M. Continuous nonsingular terminal sliding mode control based on adaptive sliding mode disturbance observer for uncertain nonlinear systems. Automatica 2019, 109, 108515. [Google Scholar] [CrossRef]
  37. Jiao, X.; Fidan, B.; Jiang, J.; Kamel, M. Type-2 fuzzy adaptive sliding mode control of hypersonic flight. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2019, 233, 2731–2744. [Google Scholar] [CrossRef]
  38. Ding, Y.; Wang, X.; Bai, Y.; Cui, N. An improved continuous sliding mode controller for flexible air-breathing hypersonic vehicle. Int. J. Robust Nonlinear Control 2020, 30, 5751–5772. [Google Scholar] [CrossRef]
  39. Sagliano, M.; Mooij, E.; Theil, S. Adaptive disturbance-based high-order sliding-mode control for hypersonic-entry vehicles. J. Guid. Control. Dyn. 2017, 40, 521–536. [Google Scholar] [CrossRef] [Green Version]
  40. Zhang, S.; Wang, Q.; Yang, G.; Zhang, M. Anti-disturbance backstepping control for air-breathing hypersonic vehicles based on extended state observer. ISA Trans. 2019, 92, 84–93. [Google Scholar] [CrossRef]
  41. Hu, Q.; Meng, Y.; Wang, C.; Zhang, Y. Adaptive backstepping control for air-breathing hypersonic vehicles with input nonlinearities. Aerosp. Sci. Technol. 2018, 73, 289–299. [Google Scholar] [CrossRef] [Green Version]
  42. Wu, G.; Meng, X. Nonlinear disturbance observer based robust backstepping control for a flexible air-breathing hypersonic vehicle. Aerosp. Sci. Technol. 2016, 54, 174–182. [Google Scholar] [CrossRef]
  43. Yu, C.; Jiang, J.; Zhen, Z.; Bhatia, A.K.; Wang, S. Adaptive backstepping control for air-breathing hypersonic vehicle subject to mismatched uncertainties. Aerosp. Sci. Technol. 2020, 107, 1–19. [Google Scholar] [CrossRef]
  44. Lu, Y. Disturbance observer-based backstepping control for hypersonic flight vehicles without use of measured flight path angle. Chin. J. Aeronuat. 2021, 34, 396–406. [Google Scholar] [CrossRef]
  45. Jiao, X.; Jiang, J. Design of adaptive switching control for hypersonic aircraft. Adv. Mech. Eng. 2015, 7, 1–10. [Google Scholar] [CrossRef] [Green Version]
  46. Li, Y.; Wu, Y.; Qu, X. Chicken Swarm—Based Method for Ascent Trajectory Optimization of Hypersonic Vehicles. J. Aerosp. Eng. 2017, 30, 1–12. [Google Scholar] [CrossRef]
Figure 1. Three views of a near-space hypersonic vehicle.
Figure 1. Three views of a near-space hypersonic vehicle.
Aerospace 09 00755 g001
Figure 2. The system structure of sliding mode backstepping controller using DPRL.
Figure 2. The system structure of sliding mode backstepping controller using DPRL.
Aerospace 09 00755 g002
Figure 3. Condition of s v under reaching laws.
Figure 3. Condition of s v under reaching laws.
Aerospace 09 00755 g003
Figure 4. Velocity tracking under reaching laws.
Figure 4. Velocity tracking under reaching laws.
Aerospace 09 00755 g004
Figure 5. Controller output β c under reaching laws.
Figure 5. Controller output β c under reaching laws.
Aerospace 09 00755 g005
Figure 6. Responses of velocity tracking.
Figure 6. Responses of velocity tracking.
Aerospace 09 00755 g006
Figure 7. Tracking errors of velocity.
Figure 7. Tracking errors of velocity.
Aerospace 09 00755 g007
Figure 8. Responses of altitude tracking.
Figure 8. Responses of altitude tracking.
Aerospace 09 00755 g008
Figure 9. Tracking errors of altitude.
Figure 9. Tracking errors of altitude.
Aerospace 09 00755 g009
Figure 10. Responses of attack angle.
Figure 10. Responses of attack angle.
Aerospace 09 00755 g010
Figure 11. Responses of pitch angle rate.
Figure 11. Responses of pitch angle rate.
Aerospace 09 00755 g011
Figure 12. Responses of mass.
Figure 12. Responses of mass.
Aerospace 09 00755 g012
Figure 13. Responses of throttle setting.
Figure 13. Responses of throttle setting.
Aerospace 09 00755 g013
Figure 14. Responses of elevator deflection.
Figure 14. Responses of elevator deflection.
Aerospace 09 00755 g014
Table 1. Parameters of NSHV.
Table 1. Parameters of NSHV.
ParameterValueUnites
Mass100,200 kg
Reference area389m
Aerodynamic chord30 m 2
Moment of inertia8,466,900 kg · m 2
Table 2. Parameters of five reaching laws.
Table 2. Parameters of five reaching laws.
TPRLDPRLERLTRLFPRL
l 1 = 1 l 1 = 1 l 1 = 1 l 1 = 1 l 1 = 1
l 2 = 1 l 2 = 1 l 2 = 1 l 2 = 1
l 3 = 1 α = 1.3 β = 0.7
α = 1.3 β = 0.7
β = 0.7 δ = 1.8
ω = 2
Table 3. Controller Parameters.
Table 3. Controller Parameters.
ParameterValue
k 1 1 0 0 1
  k 2 1.2 0 0 1.2
  k 3 1 0 0 1
  c 1 1 0 0 1
  c 2 1.5 0 0 1.5
  c 3 1 0 0 1
l 1 , 1 0.5
l 1 , 2 0.3
l 1 , 3 0.3
l 2 , 1 0.5
l 2 , 2 0.3
l 2 , 3 0.3
α 1 1.2
α 2 1.2
β 1 0.7
β 2 0.7
δ 1 1.8
δ 2 1.8
ω 1 2
ω 2 2
Γ 1 1 0 0 1
Γ 2 1 0 0 1
Γ 3 1 0 0 1
Γ 4 1 0 0 1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, S.; Jiang, J.; Li, O. Sliding Mode Backstepping Control for the Ascent Phase of Near-Space Hypersonic Vehicle Based on a Novel Triple Power Reaching Law. Aerospace 2022, 9, 755. https://doi.org/10.3390/aerospace9120755

AMA Style

Huang S, Jiang J, Li O. Sliding Mode Backstepping Control for the Ascent Phase of Near-Space Hypersonic Vehicle Based on a Novel Triple Power Reaching Law. Aerospace. 2022; 9(12):755. https://doi.org/10.3390/aerospace9120755

Chicago/Turabian Style

Huang, Shutong, Ju Jiang, and Ouxun Li. 2022. "Sliding Mode Backstepping Control for the Ascent Phase of Near-Space Hypersonic Vehicle Based on a Novel Triple Power Reaching Law" Aerospace 9, no. 12: 755. https://doi.org/10.3390/aerospace9120755

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