A Novel Hybrid Robust Control Design Method for F-16 Aircraft Longitudinal Dynamics

is paper presents a hybrid robust control designmethod for a third-order lower-triangular model of nonlinear dynamic systems in the presence of disturbance. In this paper, a novel control design is presented systematically to synthesize a robust nonlinear feedback controller, called backstepping sliding mode control (BSMC), for the proposed system by a combined approach of backstepping design and sliding mode control. In this approach, a family of the “sliding surface” is introduced in state transformations. en, a smooth switching function of the sliding surface is introduced and enforced to include in virtual feedbacks and a real control law from the control selection phrases of the backstepping design loop. e achieved control method proves a well-tracking command with asymptotic stability, provides a robustness in the presence of uncertainties, and eliminates completely a chattering phenomenon. e application of flight-path angle control corresponding to the longitudinal dynamics of a high-performance F-16 aircraft simulationmodel is implemented. Under some assumptions, full nonlinear longitudinal dynamics is reformed into a lower-triangular system for a direct application to formulate a control law. A closed-loop system is achieved for in-flight simulation with different flight profiles for a comparison of the existing methods. Also, an external disturbance on different loading/unloading conditions in flight is applied to verify and validate robustness of the proposed control method.


Introduction
A standard design method of the flight control system for nonlinear aircraft dynamic systems is based on gain scheduling from linear control system designs. In this method, a linearized model of the nonlinear flight dynamic system at key operating points within the flight envelope is achieved. From those selected points, linear controllers are designed and then combined continuously as the aircraft flies from one operating point to another. In [1], it is indicated that actual system performance and stability can be significantly different from design results due to linearized nonlinearities. erefore, a direct application of a nonlinear flight system should be addressed for a new generation of flight control design [2][3][4]. However, these works neglected the effects of uncertainties that can lead to instability and poor performance of an aircraft system due to complex variations in aerodynamic models. Recently, robust control design methods for nonlinear aircraft flight dynamics have been addressed by many researchers [5][6][7]. Twenty-eight uncertain parameters are used to model aircraft motion in nonlinear longitudinal dynamics [5]. en, a genetic algorithm is used to search a feasible design coefficient space that is evaluated by Monte Carlo for stability and performance robustness. A better approach in robust control methods [6,7] is proposed to improve the handling qualities across the specified flight envelope without the use of gain scheduling. In this approach, the control parameters of dynamic inversion are designed to minimize the probability of violating design specifications and so provide a better design with good robustness in stability and performance subject to modeling uncertainties. However, mentioned control design methods are mostly focused on design of control parameters that are not fully able to adapt internal dynamics in order to ensure the stability and performance robustness.
us, a new type of the robust controller, in which a robust property for a nonlinear dynamic system must be included, is necessary to improve stability and performance of aircraft longitudinal dynamics. In recent years, backstepping or recursive control design for a strict-feedback form of nonlinear dynamic systems is introduced firstly by authors [1] and is applied extensively for flight control of the nonlinear flight dynamic system [8][9][10][11][12][13]. In approaches [8][9][10][11][12], a full nonlinear dynamic system is transformed directly or indirectly into a strict-feedback form or a lower-triangular model for applicability via assumptions or state transformations.
e achieved model will be divided into subsystems for control objectives. By considering each subsystem in relation with interconnected subsystems, a virtual feedback control is selected from satisfying the requirement of Lyapunov's theorem. After a series of control design for subsystems, a nonlinear feedback controller is obtained. By numerical simulation and experimental study, the backstepping control method can prove a very good performance and stability if system parameter modeling is accurate. On the other hand, the method might lead to a poor quality across the full flight envelope due to uncertain system parameters.
us, a better approach [13] is proposed to improve the performance in this situation. In this design, an integrator with a tracking error is introduced and added into the control law via feedback selection. e achieved control system can improve robust performance and stability in the presence of disturbance, but time response of the flightpath angle or the attack angle is slower than that of the traditional methods. is behavior or robustness should be improved in designing the control system for a highperformance aircraft system.
Sliding mode control (SMC) was proposed firstly and applied by Slotine [14]. is method provides a robust control technique to deal with unmodeling dynamics or uncertain system parameters for nonlinear dynamic systems and has been widely used and extended for flight control design [15][16][17][18][19][20][21]. A linearized model for longitudinal dynamics is achieved for applying the sliding mode control design [17]. e achieved control method is simulated on the nonlinear flight longitudinal model of the F-16 aircraft. An extra controller of the integral sliding mode control allocation scheme [18], which is active if faults occurs, is built around a traditional controller. Although numerical simulation indicates that the proposed control shows good results for nominal and fault/failure conditions, a guideline on how it works or a stability proof is a concern in the paper. In order to solve for fuzzy stochastic systems subjected to matched/ mismatched uncertainties, Wang et al. [20] propose a new method to remove these assumptions and present a new integral sliding mode control (ISMC) method. As discussed above, it is clear that each work has some advantage and some limitation, but mostly based on the traditional control design method. So, the question arises if there exists a systematic approach to formulate a robust nonlinear feedback controller that can inherit advances of the backstepping control and robust aspects of a sliding mode control design. is paper will provide a novel solution to deal with the proposed challenge by a combined design of the backstepping method and sliding mode control. In this proposal, a family of the "sliding surface" is introduced in state transformations. en, the sliding surface and a smooth switching function are enforced to include in virtual feedbacks and a real control law from control selection phrases. Due to the inheritance of the properties of the backstepping control and sliding mode control, the achieved control method proves a well-tracking command with asymptotic stability, provides a robustness in the presence of uncertainties, and eliminates completely a chattering phenomenon. e application of flight-path angle control of a high-performance F-16 aircraft model is implemented to show advance of the proposed control method.
is paper is organized as follows: a hybrid controller formulation for a general class of a third-order lower-triangular model of a nonlinear dynamic system is presented in Section 2. Study of flight-path angle control of longitudinal dynamics of the F-16 aircraft mode is investigated to show applicability in Section 3. A numerical simulation of the full nonlinear F-16 aircraft model is implemented and compared to the integrator-backstepping control method; finally, conclusion and discussions are presented in the end of the paper.

A Hybrid Controller Formulation for a Single-Input Single-Output (SISO) Third-Order Lower-Triangular Model of Nonlinear Dynamic Systems
In this section, a problem formulation is considered firstly. Secondly, a new control design method is presented for formulating a backstepping sliding mode control (BSMC) law for the (SISO) third-order lower-triangular form of the nonlinear dynamic system. Some important remarks are given in the end of the section.

Problem Formulation.
Consider an SISO third-order lower-triangular model [1] of nonlinear dynamics systems as where the unknown valued functions f i (i � 1, 2, 3) and g i (i � 1, 2, 3) are approximated as known functions f i (i � 1, 2, 3) and g i (i � 1, 2, 3), respectively, in this paper. is approximation is based on working condition and experience for a specific system. x i (i � 1, 2, 3) ∈ R are the state variables, and the functions f i (i � 1, 2, 3) ∈ R are known real-valued, continuous, and differentiable. e functions g i (i � 1, 2, 3) ∈ R are known real-valued, continuous, differentiable, and invertible. e variable u ∈ R is the control input. e objective of this research is to synthesize a robust control law for systems (1)-(3) such that the outputs y � h(x) � x 1 track well the desired value or command, x d , asymptotically, where x d is a triple-differentiable function in time. Also, a global asymptotic stability is achieved with no or acceptably small overshoot in the presence of the model parameter errors and external disturbances.

Formulation of Backstepping Sliding Mode Control Law
ere exists coordinate transformations (4), (5), and (6) and virtual feedbacks α 1 and α 2 in (5) and (6), respectively, and real control u in (7): Such that systems (1)-(3) can be transformed into a stable linear decoupling system in terms of a new set of state variables e 1 , e 2 , and e 3 as where p 1 , p 2 , p 3 , and λ are the strict positive values and tanh(·) is a hyperbolic tangent function and defined as tanh(σ) � (e σ − e − σ )/(e σ + e − σ ). e following is a proof for eorem 1. e proof of eorem 1 consists of three steps. Coordinate transformations are introduced firstly. en, virtual feedbacks and real control can be selected in terms of recursive design approach [11]. By doing so, systems (1)-(3) can be reformed into systems (8)-(10).
Step 1. Taking derivative of equation (4) in association with In Lyapunov's approach, a virtual control law can be selected such that the tracking error e 1 (t) in equation (11) converges asymptotically to zero with small or no overshoot. e virtual controller can be chosen from satisfying the condition in which time differentiation of a control Lyapunov function (CLF) [1], L 1 (e 1 (t)), is less than or equal to a negative definitive function, w(e 1 (t)). In fact, a CLF is selected as Combining the result from time differentiation of the CLF, L(e 1 (t)), and equation (11) achieves en, the virtual control α 1 can be obtained from enforcing _ L 1 (e 1 (t)) to be a negative definitive function, or Substituting equation (14) into equation (13) results in is means that an achieved virtual control x d , p 1 ) assures to make tracking error tend to zero.
Step 2. By choosing feedback (14) and state transformations (4) and (5), systems (1)-(2) are resynthesized in (e 1 , e 2 ) coordinates as follows: Noting that the functions f 1 and α 1 are only functions of states and command. us, the term _ α 1 can be represented explicitly in variables of x 1 , x 2 , x d , _ x d , € x d , p 1 , and λ, and _ α 1 is calculated in the following equation: As a similar approach, a virtual control α 2 : � x 3 and CLF L 2 (e 1 (t), e 2 (t)) are selected to ensure that systems (16)- (17) are stabilized at origin, i.e., By taking derivative of equation (19) in consideration with equations (15)- (17), one gets en, a virtual control α 2 can be obtained from enforcing _ L 2 (e 1 (t), e 2 (t)) to be a negative definitive function, or Mathematical Problems in Engineering Substituting equation (21) into equation (20) results in Note that w 2 (e 1 (t), e 2 (t)) � p 1 e 1 (t) tanh(λe 1 (t))+ p 2 e 2 (t)tanh(λe 2 (t)) is a positive definitive function. us, _ L 2 (e 1 (t), e 2 (t)) is a negative definitive function. is means that an achieved virtual control α 2 assures to make tracking error, e 2 (t), tend to zero or the first subsystem and the second subsystem [1] are interconnected with the proposed method.
With the achieved control law in (28), substituting this control law into systems (1)-(3) results in a closed-loop system that is used to implement numerical simulation with different commands and initial conditions.

Remark 1.
With positive values of p 1 , p 2 , p 3 , and λ, the solutions e 1 (t), e 2 (t), and e 3 (t) of systems (30)-(32) are stabilized at origin (0, 0, 0) at a certain time. is means that the BSMC will enforce the outputs y � h(x) � x 1 to have a well-tracking command x d .

Remark 2.
It is observed that the tracking errors e 1 (t), e 2 (t), and e 3 (t) play a similar role as "sliding surfaces" in sliding mode control and the smooth switch function tanh(λe 2 ) ≈ tanh(λ[x 2 + p 1 (x 1 − x d ) + · · ·]) plays a "similar" role as switch function sign(·) of the sliding surface in traditional SMC. is term will ensure the robustness of the proposed method as system parameter errors or external disturbance is present in the system.

Remark 4.
Noting that the hyperbolic tangent function tanh(λσ) with a large value λ will reach to the discontinuous sign function sign(·). us, a very large value of λ might lead to a chattering problem. A suitable parameter of λ must be traded off for a good performance. (8)-(10), it is clear that the convergence speed of the proposed control algorithm will follow the exponential rule with gains p 1 , p 2 , and p 3 . e greater the gains are, the faster the convergence rate is. Besides, the convergence speed of the proposed control algorithm also depends on the function tanh(e). If you would like to gain a faster response, you can increase λ, but a trade-off between λ and chattering issue must be considered for good performance.

Flight-Path Angle Control of F-16 Aircraft Model
In this section, the lower-triangular model of nonlinear longitudinal dynamics is introduced and explained for design suitability. An explicit BSMC law is presented by using the proposed theory in Section 2. en, a numerical demonstration of flight-path angle control of the F-16 aircraft simulation model is implemented to show the advances, robustness, and applicability of the proposed method to flight dynamic systems as compared to the existing methods.

Problem Formulation.
Consider longitudinal dynamics of the aircraft flight system as shown in Figure 1.
Assumptions on aerodynamic forces [10] and use of the relationship α � θ − c are applied to transform an original longitudinal motion into the lower-triangular systems (1)- (3): where e performance specifications of the system should achieve wellbehaved command tracking with zero or small acceptable overshoot.

Derivation of Backstepping Sliding Mode Control Law.
With state variables (c, θ, q) and control input δ E , systems (33)-(35) are exactly the ones in (1) 1, 2, 3) and g i (i � 1, 2, 3) are satisfying conditions in Section 2.2. us, there is a direct application of eorem 1 in Section 2 or state transformations and feedbacks as where a virtual feedback α 1 in equation (37) is determined as in the following equation: where virtual feedback α 2 in equation (39) is determined as in the following equations: In addition, a BSMC law for systems (33)-(35) can be achieved in the following equation:

Mathematical Problems in Engineering
where α 2 is determined in (40) and _ α 2 is determined in equations (44)-(46). e achieved BSMC law in equation (43) is used for further numerical study in the next section.
By substituting the BSMC in equation (43) into systems (33)-(35), a closed-loop system is achieved for numerical simulation. Figure 2 shows a control structure design for implementing a numerical simulation. A further explanation on this block diagram can be found in [13].

Numerical Simulation of Flight-Path Angle of F-16 Aircraft Model.
For a further verification of advance and robustness of the proposed method, a full nonlinear F-16 aircraft model or a true model is used for numerical solutions. e aerodynamic data of the F-16 aircraft model are provided in [22] in which data are achieved from low-speed static and dynamic wind-tunnel tests at the NASA Langley Research Center. In this paper, the aircraft is assumed to be in level fight at Mach 0.5 and at an altitude of 25000 ft. Based on this assumption, control surfaces, aerodynamic coefficients, and other state variables at trim condition or level flight can be calculated for further simulation. For simplicity, those numbers are not listed in this paper. A closed-loop system with the BSMC law in equation (43) and a control structure in Figure 2 are used for four different flight-path command profiles away from the trim condition as follows.
Profile 1: a flight-path angle command of 10 (deg) is applied for 5 seconds; Profile 2: a multistep command is applied for the flight-path angle, in which c ref is put at 5 (deg) for the first 2.5 seconds, at 10 (deg) for the next 2.5 seconds, and then 5 (deg) in 5 seconds; Profile 3: a similar multistep command as in Profile 2 is applied for flight-path angle, but aircraft speed is varied from 600 ft/s to 500 ft/s during the simulation time; Profile 4: a doublet-step command is applied for flight-path angle for 10 seconds with change of aircraft center in which x cm � 0.3 is used in control design, but x cm � 0.35 will be used for this case for testing robustness. In the simulation, optimal positive gains of the BSMC in (43) p 1 � 2.134; p 2 � 2.213; p 3 � 2.432 are determined from the design model by using the modified genetic algorithm (MGA) [23]. e value λ � 3 is selected by the try and error technique to have a good performance. Numerical results are simulated both within and without the presence of disturbances. Also, numerical results from the BSMC-based method will be compared to the integratorbackstepping control (IBSC) method [13].

Profile 1: Single-Step Command of Flight-Path Angle.
Solid line and dash line in Figure 3 show the time response of flight-path angle of the IBSC-based control and BSMCbased control methods, respectively. Both methods can provide a well-tracking command. Results show that the settling time t s � 2 seconds of BSMC is shorter than that of the IBSC of t s � 3.5 seconds. is has been predicted in [13] due to an introduction of an integrator of the tracking error.
is is a limitation of the IBSC-based control system for a high-performance aircraft and is improved by the proposed BSMC-based method. Another advance of the proposed method is that control parameters are reduced from 4 (c 0 , c 1 , c 2 , c 3 )[13] to 3 (p 1 , p 2 , p 3 ).
is will make MGA more robust as searching for optimal control parameters. Solid line and dash line in Figure 4 show the time response of the elevator surface of the IBSC-based control and BSMCbased control methods, respectively. Results indicate that the control signal from IBSC is saturated at an initial period, and this is a reason why there is a decrease in slope for rising time of time response after that. In a similar input signal of the BSMC, a saturation occurs at t � 1.5 seconds and leads to a decrease in magnitude of flight-path angle after that. However, saturation occurs in a short time, and thus the command tracking achieves quickly in this case. Also, results in Figures 3 and 4 indicate there is no chattering problem in the proposed control method.
is is also an advance as compared to a traditional sliding mode control method [14]. Figure 5 shows time response of flight-path angle with multistep command. e dot-dash line is the flight-path angle command, the dash line is associated with the IBSC-based control method, and the solid line is associated with the BSMC-based control methods. Results indicate the IBSC method results in a poor and less robust performance for a short period of time, but the BSMC can provide a much better performance than that of the IBSC. Solid line and dash line in Figure 6 show the time response of the elevator surface of the IBSC-based control and BSMC-based control methods, respectively. In order to meet demand of a given performance specification, the control signal of the IBSC tends to generate a high magnitude in the negative side, but the behavior of the BSMC is relatively different. We can see that BSMC-based response of the elevator surface is more effective and suitable to the practical observation. is is predicted in Section 1, in which a suitable control gain might not able to adapt with the change in internal dynamics of the system.

Profile 3: Multistep Command of Flight-Path Angle with a Varying Speed of Aircraft.
In order to verify the robust performance, a multistep command is applied for flight-path angle with a change in aircraft speed that is varied from 600 ft/sec to 500 ft/sec during the 10 seconds in simulation. In this simulation condition, dash line in Figure 7 indicates the BSMC control method which can provide a good tracking for a short period. On the other hand, solid line in Figure 7 indicates the IBSC control method which can provide a good tracking for a phugoid (or longer) period of the aircraft system. is is due to a considerable effect of the integrator to improve the tracking as time increases.

Profile 4: Doublet Command of Flight-Path Angle with a Variation of Aircraft Mass
Center, x cm . In previous simulations, the mass center of the aircraft x cm � 0.3 is used for designing control parameters and in-flight simulation. In reality, the mass center might be changed due to the unload of cargoes or other equipment.
us, a variation of this parameter can further verify reliability and robustness of the control methods. In this simulation condition, the mass center is assigned to be x cm � 0.35. Figure 8 shows time response of flight-path angle for a doublet command in the presence of mass center change.
e dash line again   indicates that the BSMC results in a well-behaved command tracking and fast response with no or small acceptable overshoot.

Conclusion and Future Work
A robust and fast-response control system plays an important role in designing and developing a high-performance aircraft.
is requires an advanced control design method for nonlinear flight dynamic systems. In this paper, a hybrid robust control design method is proposed for a third-order SISO lower-triangular model of nonlinear dynamic systems in the presence of disturbances. en, the application of flight-path angle control corresponding to the longitudinal dynamics of a high-performance F-16 aircraft simulation model is implemented. e contributions of this paper are as follows: (i) e proposed control method also provides a diversity of nonlinear control design tools to a specific class of nonlinear dynamic systems and has a high potential for other application in engineering systems.     (ii) A novel controller synthesis method, called backstepping sliding mode control (BSMC) law, is presented systematically for a third-order SISO lower-triangular system. e BSMC-based control method results in an advanced control system which can provide an excellent and robust performance with asymptotic stability in the presence of uncertainties and eliminate completely a chattering phenomenon that is a drawback of a traditional SMC. (iii) A successful application is implemented for the nonlinear flight dynamic system. In particular, a BSMC for flight-path angle control of the F-16 aircraft model is derived. en, a control algorithm is provided for in-flight closed-loop simulation.
Numerical results indicate the BSMC can provide a robust performance and stability with different flight profiles and presence of external disturbances.
With the development of high-performance digital computer, robust nonlinear control methods, and computational software, the nonlinear control system will result in a considerable role in improving accuracy, flexibility, and robustness in control systems. e field of the research focuses only on a lower-triangular model of nonlinear dynamic systems. erefore, an extended version of a generalized form of the nonlinear dynamic system can prove a better solution. Also, application is limited in aircraft motion in longitudinal axis and simulation results neglect effectiveness of lateral dynamics of the aircraft. us, investigation of full nonlinear aircraft fight dynamics should be considered to validate further the proposed method. In addition, this proposed control method can be applied to other engineering systems.

Data Availability
Readers can access the data supporting the conclusions in case readers would like to repeat the results, and data can be requested by sending an e-mail to the corresponding author.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.