Research Article  Open Access
M. Bahrami, B. Ebrahimi, G. R. Ansarifar, "Sliding Mode Observer and Control Design with Adaptive Parameter Estimation for a Supersonic Flight Vehicle", International Journal of Aerospace Engineering, vol. 2010, Article ID 474537, 9 pages, 2010. https://doi.org/10.1155/2010/474537
Sliding Mode Observer and Control Design with Adaptive Parameter Estimation for a Supersonic Flight Vehicle
Abstract
Design and synthesis of a nonlinear generic supersonic flight vehicle longitudinal dynamics control for angleofattack, AOA, output tracking in the atmospheric flight is presented based on sliding mode control. A sliding mode observer is invoked to estimate AOA which is difficult to measure in practice. Large parameter uncertainties accommodation envisaged by designing adaptive mechanisms for both the control and observer and high chattering authority due to large deviations of aerodynamic coefficients arising from windtunnel measurements are inhibited. The employed method enables the sliding mode control design to exhibit the desired dynamic properties during the entire outputtracking process. Simulations results are presented to demonstrate the performance, robustness, and stability.
1. Introduction
Tailcontrolled skidtoturn flight vehicles are commonly agile systems that should be equipped with control systems whose tasks are to provide stability, disturbance attenuation, and reference signal tracking, while their aerodynamic coefficients vary over a wide dynamic range due to large Machaltitude fluctuations and due to aerodynamic coefficient uncertainties resulting from inaccurate windtunnel measurements.
It is common practice, when designing a control system for a nonlinear flight vehicle, to represent the flight envelope by a grid of Machaltitude operating points and then to perform a linearization of the nonlinear state equations at trim points of the gridded flight envelope. The plant in fact becomes a differential inclusion under continuously varying flight conditions. There are many possible ways of dealing with the control of such linear and time varying plants. The classical approach is to design a controller for a certain point and then to schedule the controller’s gain to place eigenvalues deep in the left half plane and the near realaxis according to a measured or derived parameters that represent flight conditions, such as angle of attack or Mach number [1]. In another method, methods are invoked to design a collection of controllers, where, for each operating point in the flight envelope grid, a controller with a fixed structure results [2]. The ensuing set of controllers is then transformed to a single gain scheduled controller by obtaining a leastsquare fit of its parameters with respect to angle of attack, or Mach number, and so forth. For highly agile air vehicles, these techniques would result in an extensive number of controller design points to be able to cope with the drastically changing dynamic behavior throughout the flight envelop.
All of the aforementioned methods are linear design techniques that require either exact knowledge of system parameters or, alternatively, assumption of some uncertainty model such as norm boundedness, thus allowing for robust controller design. On the other hand, sliding mode control may, in principle, be implemented for dynamic systems having only a qualitative description and a number of inequality restrictions. As will be shown in sequel, this approach will provide a simple and highperformance controller on nonlinear system using straightforward design procedure.
Sliding modes are the primary form of operational variable structure systems. A sliding mode is a motion on a discontinuity set of a dynamic system and is characterized by a suite of feedback control law and a decision rule known as switching function. Such modes are used to maintain the given constraints with utmost accuracy. The sliding mode controller is an attractive robust control algorithm because of its inherent insensitivity and robustness to plant uncertainties and external disturbances.
The nonlinear pitch dynamics of a hypothetical tailcontrolled flight vehicle is derived in Section 2. In Section 3 sliding manifold is employed to control pitchaxis dynamics of the air vehicle. Adaptive sliding mode control is designed in Section 4. Sliding mode observer is designed and utilized in control scheme in Section 5 on the basis of pitch rate measurement for the nonlinear system to estimate angleofattack which is difficult to measure in practice. Integration of sliding mode control with adaptation mechanism based on the designed observer is presented in Section 6.
2. Problem Formulation
The missile model employed in this analysis is based on the hypothetical tailcontrolled air vehicle which has been used as a benchmark in a number of recent studies on nonlinear design techniques [1–3]. The control objective is to force the air vehicle to track a desired motion path generated by the guidancenavigation system as the reference angleof attack. The problem is first formulated and performance objectives are specified.
2.1. Dynamic Model of the Airframe
The model assumes constant mass, no roll rate, zero roll angle, no side slip, and no yaw rate. Under these assumptions, the longitudinal nonlinear equation of motion is reduced to two forces and one moment. Using body axis coordinates, Figure 1, these equations are
The dynamic pressure is defined as
Note that velocity and air density are not assumed to be constant or slowvarying, but a standard atmospheric model is assumed in simulation based on previously reported data [4].
Aerodynamic polynomials resulting from windtunnel measurements are given as [1] The numerical values of the various constant parameters in the dynamic equations are given in Table 1. Note that all coefficients are dimensionless and all angles are in radians.
Consequently, the differential equations describing the body motion are The tailfin actuator dynamics describing the tail deflection is where actuator fin deflection is limited to Table 2 gives the flight vehicle constants.

As shown in the equations of motion, the timevarying aerodynamic parameters contribute heavily to the variation of dynamic forces and moments exerted on the vehicle airframe. It is assumed here as is often customary in air vehicle autopilot design that the vehicle velocity (or Mach number) is constant and the nonlinear state equation associated with is dropped from the design model [1]. In this case, the vehicle velocity becomes an independent (external) parameter upon which the state dynamics depend.
By rewriting (4) and (5), it is easy to show that where , . Numerical values of and are 0.0206 and 1.2326 , respectively, at an altitude of 6100 m.
2.2. Performance Objective
The control objective is to force the air vehicle to track a desired motion path. To this end, the guidancenavigation system generates the reference angleofattack, denoted by This reference signal is compared with the actual signal measured in order to reveal a tracking error. The control problem consists of generating a tail deflection that corresponds to the required maneuver. Considering the air vehicle flying at an altitude of 6100 m (20,000 ft), the closedloop system should maintain stability over the operating range.
3. Control Design
The longitudinal model of the flight vehicle described by (10) is a special case of a general nonlinear system of the form where f, g, and h are sufficiently smooth functions.
Remark 1. Relative degree, for nonlinear system (11), is defined as follows [5]:
Utilizing (12), the relative degree of the system (10) can be obtained as Hence, the switching surface can be designed as where and λ is strictly positive constant defining the bandwidth of the error dynamics. Differentiating s, we have In (14) derivatives of the desired angleofattack, generated by guidance subsystem can be obtained by conventional numerical approaches and there would be no measurement noise which is common for navigation sensors. The sliding mode autopilot design is then choosing the control input such that the following attractive equation is satisfied: where is strictly positive constant that determines the desired reaching time to the sliding surface. The attractive equation (15) is also called sliding condition and implies that the distance to the sliding surface decreases along all system trajectories. Furthermore, the sliding condition makes the sliding surfaces an invariant set; that is, once a system trajectory reaches the surface, it will remain on it for the rest of time. In addition, for any initial condition, the sliding surface will be reached in a finite time. When (14) is used, the controller that satisfies the sliding conditions (15) can be chosen as This control can be viewed as consisting of two parts: the term, called equivalent control, which guaranties for the nominal model, and the other term, incorporated to deal with uncertainties. However, because of the discontinuity across the sliding surfaces, the preceding control law may result in control chattering. As a practical matter, chattering is undesirable phenomena because it involves very high control action and may excite highfrequency dynamics neglected in the modeling. The discontinuity in the control law can be dealt with by defining a thin boundary layer of width around the sliding surface, that is, replacing sign(s) with continuous saturation functions where otherwise. Thus, our sliding controller is modified as follows: Figures 2 and 3 show the response of the nominal model, without parameter uncertainty. In the simulation studies for the nominal model, is taken as 10 and is chosen as 0.05. Figure 4 shows the simulation results for the same maneuver albeit in the four various possible parameter uncertainties. The results show that robustness and stability have indeed been achieved. The price paid to achieve robustness in this case is in the form of large gain ρ and control chattering. However, the control gain cannot be chosen arbitrary large due to practical considerations such as reaching control surface deflection limits. Furthermore, control chattering is undesirable in practice because it involves high control activity and may excite highfrequency unmodeled dynamics. In simulation, the width of the boundary layer has expanded to reduce chattering. To eliminate chattering, the width of the boundary layer has to be further expanded, which will result in large steadystate errors. To circumvent this drawback an adaptation mechanism will be utilized in sequel.
4. Adaptive Sliding Mode Control Design
The sliding mode controller exhibits good tracking in the presence of parametric uncertainty at the expense of high gains and control chattering. The performance of the sliding mode controller is dependent on the size of the parametric uncertainties involved. An undesirable characteristic of this controller is the quick erosion of performance to gain robustness. This shortcoming of the sliding controller has urged combining the given control law with online parameter adaptation. In this section, an adaptive sliding mode controller will be developed that combines an online parameter estimator with the sliding mode controller of the preceding section. The adaptive laws for updating parameters are generated commonly using the Lyapunov approach. The dynamics of the system is more sensitive to the uncertainty parameters in the aerodynamic coefficients achieved by windtunnel measurements. C_{N} embodies parameters , and embodies parameters , and 1,2,3. To deal with the parametric uncertainty, we will combine the sliding mode controller with an online parameter estimator forming an adaptive sliding mode controller: where
In order to proceed with adaptive parameter estimation for aerodynamic coefficients we will follow as where The switching dynamics based on the control input (18) can be written as where, , and is the resulting regression vector: The adaptive laws can be derived using the Lyapunov synthesis. Consider the Lyapunovlike function: where , and is a definite positive scalar. With the assumption that the variation of is small, (i.e., and , the following can be deduced: Performance of parameter estimations based on the achieved dynamics in (25) has been demonstrated in Figures 5 and 6. The flight control behavior under the parameters variation has been shown in Figures 7, 8, and 9.
5. Sliding Mode ObserverControl Design
The sliding mode controller developed in the preceding sections assumes that full states are available for measurement. In practice, however, only the states that correspond to pitch rate are expected to be measured. The state variables corresponding to the angleofattack may be difficult to measure in an actual supersonic flight vehicle. Accurate measurements are costly and difficult if at all possible to obtain practically. In this section, a sliding mode observer is designed that can provide estimation of the angleofattack based on the measurement of pitch rate. When the techniques developed in [6–8] are applied, a sliding mode observer is designed whose structure is given by where and .
and are positive constants selected as 20 and 0.6. and are sliding gains and chosen as 0.0001 and 0.001, respectively. The sliding surface of the observer is defined by Using the same procedure as in the preceding section the dynamics of the observer can be obtained as The control input then will be conducted as where and Performance of the observer for estimation of the angleofattack has been illustrated in Figure 9. By employing observed angleofattack from designed sliding mode observer and utilizing the tracking error as , output tracking of the system, control input, climb rate, and horizontal speeds can be achieved for the nominal system as shown through Figures 10, 11, 12, to Figure 13, respectively, where the initial condition for the observer has been considered as . Figure 14 shows the simulation results for four various possible parameter variations in the system that influences performance of the observer. The results show that robustness and stability have been achieved.
6. Adaptive Sliding Mode ObserverControl Design
In this section adaptive parameter estimation will be conducted on the sliding mode observer designed in the preceding section. This will inhibit performance decrease against parameter variation introduced by C_{N} and C_{M}. Now (20) can be rewritten as where Consider Lyapunov synthesis as With the assumption that the variation of is small (i.e., and ), the following can be deduced Performance of parameter estimations based on the achieved dynamics in (33) has been demonstrated in Figures 15 and 16. The flight control behavior under the parameters variation has been shown in Figures 17 and 18. Performance of the adaptive sliding mode controlobserver demonstrated both for the nominal and perturbed system accordingly. Output tracking errors for the perturbed system following the designed approaches are shown in Figure 19 with robustness against parameters uncertainties.
7. Conclusion
A robust control is designed for a hypothetical nonlinear supersonic flight vehicle using sliding mode and adaptive control theories based on the Lyapunov theory. First a robust sliding mode control is designed based on the assumption that all system states are available. An adaptation mechanism is then integrated with the robust control to circumvent large parameters uncertainties and decrease both chattering phenomena and boundary layer thickness expansion which exhibits large steady output tracking error. In order to deal with shortcomings arising from system states unavailability a sliding mode observer is developed to estimate angleofattack which is difficult to measure in practice. The sliding mode observer finally is integrated with an adaptation scheme that makes system to be adaptive against large parameter variations.
Nomenclature
:  Axial force dimensionless aerodynamic 
:  Normal force dimensionless aerodynamic polynomial 
:  Pitching moment dimensionless aerodynamic polynomial 
:  Reference diameter, m 
:  Axial force, N 
:  Normal force, N 
:  Gravity, m/s^{2} 
:  Altitude, m 
:  Moment of inertia, kg m^{2} 
:  Vehicle constant parameters 
:  Vehicle mass, kg 
:  Mach number 
:  Pitching moment, 
:  Static pressure, N/m^{2} 
:  Pitch rate, rad/s 
:  Dynamic pressure, kg/m s^{2} 
:  Reference surface area, m^{2} 
:  Speed of sound, m/s 
:  Speed of vehicle, m/s 
:  Angleofattack, rad 
:  Angle of attack component of normal force dimensionless aerodynamic polynomial coefficients 
:  Angle of attack component of pitching moment dimensionless aerodynamic polynomial coefficients 
:  Actual tail fin deflection, rad 
:  Commanded tail fin deflection, rad 
:  Pitch angle, rad 
:  Tail actuator time constant, sec. 
References
 A. R. Mehrabian and J. Roshanian, “Skidtoturn missile autopilot design using scheduled eigenstructure assignment technique,” Proceedings of the Institution of Mechanical Engineers G, vol. 220, no. 3, pp. 225–239, 2006. View at: Publisher Site  Google Scholar
 R. A. Nichols, R. T. Reichert, and W. J. Rugh, “Gain scheduling for Hinfinity controllers: a flight control example,” IEEE Transactions on Control Systems Technology, vol. 1, no. 2, pp. 69–78, 1993. View at: Publisher Site  Google Scholar
 J. J. Zhu and M. C. Mickle, “Missile autopilot design using a new linear timevarying control technique,” Journal of Guidance, Control, and Dynamics, vol. 20, no. 1, pp. 150–157, 1997. View at: Google Scholar
 US Standard Atmosphere, US government printing office, Washington, DC, USA, 1976.
 A. Isidori, Nonlinear Control Systems, Springer, New York, NY, USA, 1995.
 J.J. E. Slotine, J. K. Hedrick, and E. A. Misawa, “On sliding observers for nonlinear systems,” ASME Journal of Dynamics Systems and Control, vol. 109, pp. 245–252, 1987. View at: Google Scholar
 I. Haskara, U. Ozguner, and V. Utkin, “On variable structure observers,” in Proceedings of the IEEE International Workshop on Variable Structure Systems (VSS '96), pp. 193–198, December 1996. View at: Google Scholar
 V. I. Utkin, J. Guldner, and J. Shi, Sliding Modes in Electromechanical Systems, Taylor & Francis, London, UK, 1998.
Copyright
Copyright © 2010 M. Bahrami 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.