Mathematical Problems in Engineering

Volume 2015 (2015), Article ID 237453, 16 pages

http://dx.doi.org/10.1155/2015/237453

## Nonlinear Constrained Adaptive Backstepping Tracking Control for a Hypersonic Vehicle with Uncertainty

^{1}State Key Laboratory of Metastable Materials Science and Technology, School of Mechanics, Yanshan University, Qinhuangdao, Hebei 066004, China^{2}School of Science, Yanshan University, Qinhuangdao, Hebei 066004, China^{3}School of Electric Engineering and Automation, Tianjin University, Tianjin 300072, China

Received 21 May 2015; Accepted 14 September 2015

Academic Editor: Yang Tang

Copyright © 2015 Qin Zou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The control problem of a flexible hypersonic vehicle is presented, where input saturation and aerodynamic uncertainty are considered. A control-oriented model including aerodynamic uncertainty is derived for simple controller design due to the nonlinearity and complexity of hypersonic vehicle model. Then it is separated into velocity subsystem and altitude subsystem. On the basis of the integration of robust adaptive control and backstepping technique, respective controller is designed for each subsystem, where an auxiliary signal provided by an additional dynamic system is used to compensate for the control saturation effect. Then to deal with the “explosion of terms” problem inherent in backstepping control, a novel first-order filter is proposed. Simulation results are included to demonstrate the effectiveness of the adaptive backstepping control scheme.

#### 1. Introduction

Air-breathing hypersonic vehicles (AHVs) are characterized by their unique design, incorporating a supersonic combustion ramjet engine located beneath the fuselage. This esoteric configuration results in strong coupling between the thrust and pitch dynamics of the vehicle, which in combination with flexible effects and static instability make the vehicle a challenging application for control [1]. In addition, there are sensitivity changes in the flight conditions, uncertain aerodynamic characteristics of the vehicle, and highly nonlinear nature of hypersonic vehicle dynamics. Thus the problem of control design is one of the key techniques for the application of flexible air-breathing hypersonic vehicle (FAHV), and the control system is required to have robustness to uncertainty. Since there is difficulty in accurately measuring and estimating vehicle’s aerodynamic characteristics, only the longitudinal analytical model of FAHV proposed in [2] has been used for controller design. For the enormous complexity of the nonlinear dynamics of FAHV, linear control theory has been widely employed for flight control design based on a linearized model [3–7]. These controllers are designed based on a linearized model which is obtained at specified trim condition or obtained by feedback linearization technique. With the development of nonlinear control theory, nonlinear control schemes are used to design the controller of FAHV [8–14].

Though the research mentioned above achieved satisfactory control performance, it has not considered input saturation. And it usually appears in many practical systems because the amplitudes of control inputs of almost all practical control systems are limited. The closed-loop system performance may be degraded severely or even lose stability if the input constraint is ignored. Some control methods are applied to handle input constraints [15–17]. For the flight control system, under the occurrence of input saturation, it is in open loop state. If the output of actuator does not return to linear work space, the hypersonic vehicle may lose stability or even disintegrates. Thus it is necessary to design the high reliability control system with input constraints. Many control approaches have been presented to design the controller for the hypersonic vehicle with the consideration of input constraints. Anti-windup control was developed to handle input constraint of hypersonic vehicles while the uncertainty was not considered [18]. Model predictive control has been used popularly because of its inherent capability to implement input constraint directly at the level of control design [19]. However, it depends on the real-time receding horizon optimization, and the main barrier of it applied to hypersonic vehicle is online optimization and the determination of time-domain step size [20]. approach was proposed for a linearized FAHV model in the presence of uncertain parameters and input constraints, where the linearized model was obtained by the feedback linearization approach [21]. It should be pointed out that high-order derivatives of outputs need to be computed. By using the differential geometry principle and the total energy theory, advanced flight control laws were designed for hypersonic vehicle in the presence of actuator limitations [22]. Three adaptive fault control schemes were proposed for AHV in considering external disturbances, actuator faults, and input saturation [23]. The latter two control approaches did not need to know the upper bound of the external disturbances and the real minimum value of actuator efficiency factor in advance. An adaptive backstepping attitude controller was proposed for reentry RLV with input constraint and external disturbance in [24]. And in [25], an adaptive dynamic surface controller was proposed for a generic hypersonic flight vehicle with consideration of magnitude, rate, and bandwidth constraints on actuator signals. Then in [26], a novel integral term was introduced during dynamic surface control (DSC) scheme design procedure to improve the tracking performance of designed controller and avoid a large initial control signal. Moreover, a robust adaptive dynamic surface controller was investigated for a hypersonic vehicle in the presence of parametric model uncertainty and input saturation, where a compensation design was employed when the input saturations occurred [27]. An adaptive DSC scheme based on radial basis function neural network (RBFNN) was presented for a hypersonic vehicle under the magnitude, rate, and bandwidth constraints on actuator in [28]. Furthermore, nonlinear disturbance observer and RBFNN based sliding mode control were designed for a near space vehicle in [29], where RBFNN was constructed as a compensator to avoid the saturation nonlinearity of rudders.

The motivation of the research is to develop practical nonlinear robust control method for a FAHV model with aerodynamic uncertainty. The main contributions are summarized as follows.

First of all, a nonlinear control-oriented model derived from the curved-fitted model of FAHV without obtaining the linearized model at a trim condition or computing the high-order time derivatives of outputs. Based on the analysis of control-oriented model, it is reasonable to decompose it into two low-order subsystems: velocity subsystem and altitude subsystem. Then the available inputs are designed for subsystems with the reduction of computational burden.

In the second place, auxiliary signals are introduced to cope with input constraints which are provided by auxiliary systems, and the auxiliary signals are employed during the controller design and stability analysis procedure. Although input constraint is also handled by the additional system in [30], in this paper it does not need to construct dynamic robust term during the controller and additional system design procedure. So it simplifies the controller design and stability analysis.

Last but not least, the upper bound of uncertainties is not required to be known in advance. Adaptive law is designed to estimate the upper bound, and the robustness is ensured at the same time. The difference from our previous work [30] is that adaptive technique is employed to approximate compounded uncertainty. The estimation ability of the adaptive law can be clearly shown via theory analysis, and the estimation accuracy can be improved by choosing parameters. From theoretical and simulation aspects, the parameters of adaptive laws are determined more simply than that of RBFNN used in [30]. The “explosion of terms” problem is avoided by developing the novel first-order filter, and its advantage over the traditional first-order filter in DSC method is testified by simulation result.

#### 2. Problem Formulation

##### 2.1. Curved-Fitted Model of FAHV

The nonlinear equations of motion of FAHV used in this study are mentioned in [31]. The longitudinal dynamic equations of a FAHV, which describe velocity, altitude, flight path angle (FPA), angle of attack (AOA), pitch rate, and flexible modes are given as follows:

This model is called a curve-fitted model (CFM) for FAHV and utilized for simulation only. In (1)–(6), the thrust , drag , lift , pitching moment ,and three generalized forces , , are complex algebraic functions of both system states and inputs that must be simplified to render the model analytically tractable. The flexible states , , are related to the deflections of the fore-body turn angle and aft-body vertex angle , denoted by and , respectively. The approximations of the forces and moments are the same as those provided in [32], which can be expressed asthe coefficients obtained from fitting the curves are given in the following:Here, the function arguments are removed in brief, is the free stream Mach number, and is the dynamic pressure; they are defined as and . Herein, is the altitude dependent air-density, and is the speed of sound at a given altitude and temperature.

##### 2.2. Control-Oriented Model and Control Objective

As discussed in [5], the simplified model is developed for controller design. It is called control-oriented model (COM) and derived from the CFM by neglecting the flexible dynamics since the measurements of the flexible states are not assumed to be available for feedback [13]. That is to say, the flexible dynamics are not considered directly at the control design level but are taken as perturbations on the COM, and their effects are evaluated in simulation. Besides, the input coupling terms are represented as perturbations. Namely, the COM is comprised of five rigid body dynamic equations (1)–(5).

It is noted that in this model a canard is used to cancel the lift-elevator coupling to cancel nonminimum phase. The canard deflection is a function of elevator deflection so that , where , is an interconnect gain. From (8), it is obvious that if is precisely specified, the canard would exactly cancel the lift due to the elevator deflection. However, in practice, an ideal interconnect gain is hard to achieve and thus it is impossible to exactly cancel the lift-elevator coupling [33]; it is assumed that is uncertain. Herein, the uncertainty of aerodynamic parameters except , and is considered. For example, the uncertainty of aerodynamic parameter is denoted as . Based on it, we can obtain the following expressions of aerodynamic uncertainty:wherewith

For the hypersonic vehicles control system, it is unavoidable that the actuator output is limited, especially the magnitude constraint of actuator input. The input constraints studied herein include the constraints on fuel equivalence ratio and elevator deflection. Since the propulsion system of hypersonic vehicle is required to maintain the conditions that sustain scramjet operation, the limit on fuel equivalence ratio is naturally induced. The thermal choking will occur if the constraint is violated, which could lead engine to be unstart which could jeopardize the mission, the vehicle, and its contents. The constraint on elevator deflection is mainly caused by the limits on control surface displacement. Because the actuator outputs are constrained, the input constraint of fuel equivalence ratio and elevator deflection are denoted as and , respectively.

The expression of input saturation is as follows:where is the desired control input to be designed in the following section and and are the minimum value and maximum value of fuel equivalence ratio, respectively.

The expression of input saturation is as follows:where is the desired control input to be designed in the following section and and are the minimum value and maximum value of elevator deflection, respectively.

There are four inputs in (1)–(5) and they are the diffuser-area-ratio (it is fitted as in this study), canard deflection, fuel equivalence ratio, and elevator deflection. The outputs to be controlled are selected as velocity and altitude. It is assumed that the states of the rigid body system are available and the controller design only utilizes the feedback from the rigid body states. Because the measurements of the flexible states are not assumed to be available for feedback [13], the flexible states are treated as disturbances. The control objective is to design fuel equivalence ratio and elevator deflection to make velocity and altitude track their command trajectories with aerodynamic uncertainty and input saturation.

#### 3. Controller Design

It can be obtained from the aircraft model (1)–(5) and the aerodynamic formulations (7)-(8) that thrust affects the velocity and the fuel equivalence ratio affected by the thrust, so the velocity is mainly affected by fuel equivalence ratio. Moreover, it is reasonable to implement a separate control design, since elevator deflection has a dominant contribution towards the altitude change. The COM is decomposed into two subsystems firstly, and they are the velocity subsystem and the altitude subsystem. Then every subsystem is controlled separately by the available input. Dynamic inversion and robust adaptive control are synthesized to design control input (fuel equivalence ratio) for the first subsystem. Backstepping control and robust adaptive control are combined to design control input (elevator deflection) after the second subsystem is transformed into strict-feedback form. With the consideration of input saturation, the auxiliary signals which are provided by the auxiliary system are applied to cope with them. The auxiliary signals are employed during the controller design and stability analysis procedure. The detailed design procedures are given in the following two subsections.

##### 3.1. Robust Adaptive Control for Velocity Subsystem

Dynamic inversion can achieve the decoupling between the input and output of a system. It can be used for different types of aircraft and adapted to a model change from simulation and flight tests, but it lacks robustness to uncertainty. The adaptive control can overcome the uncertainty, disturbance, and unmodeled dynamics. Therefore, in this subsection, the robust adaptive control is incorporated into the dynamic inversion to design the fuel equivalence ratio.

The dynamic of velocity can be written aswhere

It is obvious that the desired control input may be larger than the actual control energy provided. Thus there is a difference between the desired control input and the actual control input, and it is described as

According to physical backgrounds of FAHV, it is reasonable to make the following assumption.

*Assumption 1. *For the uncertain term , there is a constant such that .

The tracking error of velocity is defined aswhere is the reference command of velocity. The time derivative of (17) is

Considering the saturation of fuel equivalence ratio (12), inspired by the work in [34], the following auxiliary system is constructed:where .

The fuel equivalence ratio is designed asThe adaptive law for iswhere , .

Taking the estimation error and tracking error into account, the Lyapunov function is constructed aswhere is the estimation error of .

Using the derivative of respective to time,

From (18)–(21), we haveSincethenas long as , ,

##### 3.2. Robust Adaptive Backstepping Control for Altitude Subsystem

The altitude subsystem includes dynamic equations of altitude, FPA, AOA, and pitch rate. The altitude is controlled through the tracking of FPA reference command , which is derived from altitude reference command .

The tracking error of altitude is defined aswhere is altitude reference command.

The dynamic of (28) is . Because FPA is very small during the cruise phase, , then

The FPA reference command is defined as and the tracking error of FPA is . Then (29) becomesAnd the FPA reference command is chosen as where is the parameter to be designed.

The Lyapunov function is constructed asBased on (30) and (31), the time derivative of (32) satisfiesIn the next step, the time derivative of needs to be computed. But from (31), the term is needed, and there is aerodynamic uncertainty in the formulations , so it is difficult to compute . Thus, the derivative of is estimated by the following first-order filter:where are the filter estimation error, is the filter time constant, and are constants.

*Remark 2. * is approximated by the first-order filter (34), and is used for the controller design in the next step. Obviously, with assumed to be zero, the first-order filter is reduced to a classical integral filter as used in DSC method . With this filter, the measurement noise can be eliminated in the virtual control effort. Compared with the classical first-order filter in the DSC method, the fast transient response of filter can be obtained. Moreover, compared with the filter that includes saturation function component, the hyperbolic tangent function component is superior to the saturation function component; thus the performances of the filter proposed in this paper can be improved.

The dynamics of FPA, AOA, and pitch rate are rewritten as the following strict-feedback formulation:where

The difference between the desired control input and the actual control input is described as

It is noted that the structure of (35)–(37) possesses a strict-feedback form, where the uncertain terms do not satisfy the matched condition, and it makes backstepping control philosophy applicable. Here, robust adaptive control is incorporated into backstepping control to design control input (elevator deflection) and the states AOA and pitch rate are taken as the virtual control inputs. The order of the altitude subsystem is four; it will induce repeated differentiations of virtual control inputs and may cause “explosion of terms” problem. What is more, the time derivatives of virtual control inputs are needed in the next step, but there are nonlinearity and uncertainty in (35)–(37); it is difficult to obtain the time derivatives. And it may cost large computational load even if the time derivatives can be computed. To cope with this situation, the time derivatives of virtual control inputs are estimated by the novel first-order filter.

According to physical backgrounds, it is reasonable to make the following assumption.

*Assumption 3. *For the uncertain terms , , , there exist constants such that , , .

From (35), the dynamic of tracking error of FPA is written asThe virtual control input is designed aswith the adaptive law of where are parameters to be designed.

Considering the estimation error and tracking error, the Lyapunov function is constructed aswhere is estimation error, .

The time derivative of (43) isFrom (34), filter error, and [35], if , the term yieldswhere is a constant that satisfies ; that is, .

Based on (41), (42), (45), and , (44) yields

The derivative of is estimated by the following first-order filter:where are the filter estimation errors, is the filter time constant, and are constants.

Define the error signal of AOA as From (36), the time derivative of isThe virtual control input is designed aswith the adaptive law of ,where are parameters to be designed.

Considering the estimation error and tracking error, the Lyapunov function is constructed aswhere is the estimation of and is the estimation error.

From (50) and (51), the time derivative of (52) isfrom (47), filter error, and [35], if , the term yields Based on (54) and , (53) yields

The derivative of is estimated by the following first-order filter:where are the filter estimation errors, is the filter time constant, and are constants.

The error signal of pitch rate isOn the basis of (37) and (39), the time derivative of (57) isInspired by the work in [34], an auxiliary system is used to handle saturation of elevator deflection (13): where .

The elevator deflection is designed asThe adaptive law for iswhere ,

*Remark 4. *In the traditional sliding mode control, the* sign* function may cause chattering problem, and it may induce that the virtual control input in the backstepping control cannot be tracked accurately by the next subsystem. What is worse, for the hypersonic vehicle system, the chattering may lead to the disintegration of the vehicle rudder. To eliminate the chattering problem, the continuous robust term with norm-type switched function is used during control input design.

Considering the estimation error and tracking error, the Lyapunov function is constructed aswhere is the estimation of and is the estimation error.

Using the derivative of respective to time,If , the term satisfies From (59)–(61), the inequality (63) satisfiesSincethe following inequality holds:

*Remark 5. *This paper contains the following different aspects compared with existing results [26, 27, 30] that investigate the controller design of FAHV. (a) The robustness of the scheme developed herein is shown through aerodynamic uncertainty, whereas the robustness of the designed scheme in [26] is evaluated through different fuel levels. Although in [30] input constraint is handled by the additional system, we do not need to construct dynamic robust term when designing controller and additional system in this paper. It simplifies the controller design and stability analysis. Moreover, the novel first-order filter that is different from the filter used in [30] is developed to handle “explosion of terms” problem. (b) The difference of this paper compared with [27] is as follows. On the one hand, the two subsystems do not needed to be transformed into linear parameterized form. On the other hand, input constraint is handled by constructing additional system in this paper, whereas the application of compensation technique in [27] coped with input constraint.

The Lyapunov function for the altitude subsystem is constructed asUsing the derivative of it respective to time, thenFrom (27), (33), (46), (55), and (67), the following inequality holds:As long as , , , , , , the above inequality yieldsHere and + .

The convergence domain of can be expressed as the following compact set:So far, the stability of the subsystem is proven.

#### 4. Simulations

Simulations are carried out to illustrate the effectiveness of the robust adaptive backstepping control scheme proposed in the previous section. The equations of motion (1)–(6) are used for simulations. The vehicle model parameters and the initial flight condition of the vehicle dynamics are referred to in [36]. The fuel level is assumed to be 50%. Input constraints that are used to test the capability of the developed controller to handle the input constraints are set as , , , and . Parameters adopted for the control inputs, adaptive laws, and auxiliary systems are given as follows: , , , , , , , , . , , , , .

As shown in adaptive laws (21), is monotone increasing. Its overincrease may cause the control input to largely increase. So the adaptive law is revised to suppress the overincreasing of :It is noted that after the error signal entering the stability region, the estimation change rate remains at zero, which means that will not change; then the control input will not be too large.

To test the performance of the designed strategy, the maneuver simulation is carried out, where the maneuver uses separate reference commands of altitude and velocity. The reference command of velocity is chosen as 2500 ft/s. The reference command of altitude is represented as . Here and is generated as follows: As demonstrated in (74), the maneuver begins at 100 s. The vehicle climbs about 5000 ft between 100 s to 550 s and descents about 5000 ft between 650 s and 1100 s. Finally, the altitude remains constant after 1100 s. In addition, reference commands are generated from a second-order prefilter with a natural frequency rad/s and a damping ratio . Two uncertain cases are included: case 1: 10% of uncertainty of the aerodynamic parameters is taken into consideration; case 2: 20% of uncertainty of the aerodynamic parameters is taken into account.

The simulation results are shown in Figures 1–5. Besides, some local time responses are given in the corresponding figures to clearly demonstrate the dynamic process. Altitude and velocity achieve their stable tracking of their respective reference commands as given in Figure 1. It is noted that the FPA reference commands are approximately the same in case 1 and case 2 from Figure 1(c). As shown in Figure 2, the control inputs remain within their constraints. Other states are shown in Figures 3 and 4. In two cases, the pitch rate converges to 0 rad/s, and three flexible states are bounded between the values −0.1 and 0.8 during the whole maneuver. Compared simulation result between the proposed first-order filter of this paper and the traditional first-order filter in DSC method is given in Figure 5. It is clear that the estimation errors of the proposed filters in this paper are all smaller than that of the traditional first-order filters in DSC method.