Reliable Robust Flight Tracking Control: An LMI Approach
Reliable Robust Flight Tracking Control: An LMI Approach
Reliable Robust Flight Tracking Control: An LMI Approach
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 10, NO. 1, JANUARY 2002
AbstractThis paper studies the reliable robust tracking controller design problem against actuator faults and control surface
impairment for aircraft. First, models of actuator faults and
control surface impairment are presented. Then a reliable robust
tracking controller design method is developed. This method
tracking
is based on the mixed linear quadratic (LQ)/
performance indexes and multiobjective optimization in terms of
linear matrix inequalities. Flight control examples are given, and
both linear and nonlinear simulations are given.
Index TermsActuator faults, control surface impairment,
flight control, linear matrix inequality (LMI), reliable control,
tracking.
I. INTRODUCTION
N THE design of the advanced tactical fighter (ATF), reliability, maintainability, and survivability are three basic
design requirements. A fighter aircraft with reliability should
be able to maintain a satisfactory closed-loop system performance in the presence of sensor/actuator faults or control surface impairment in the battle environment. Aircraft has certain
inherent redundancy built in. For example, conventional aircraft
has mainly three sets of control surfaces: the elevator (horizontal
stabilator) for pitch control; the differential aileron for roll control; and the rudder for yaw (directional) control. For ATF, there
are usually also canard and flaps for pitch and yaw control. If the
movement of the left and right elevators and ailerons are made
independent of each other, the inherent redundancy of aircraft
control surfaces can be made use of to tolerate the loss of certain
control surfaces. To a certain extent, the left and right elevators
provide roll (roll and pitch) control if they deflect the same (different) amount toward different direction, and the left and right
ailerons also provide pitch (pitch and roll) control if they deflect
the same (different) amount toward the same direction.
Currently there are several fault-tolerant/reliable flight
control system design conceptions, such as self-repairing flight
control systems (SRFCS) [1], [17], reliable flight control systems [18], [26], and so on. SRFCS design conception involves
such procedures as real-time fault detection, isolation (FDI),
Manuscript received May 29, 2001. Manuscript received in final form October 29, 2001. Recommended by Associate Editor A. Ray. This work was supported in part by the Academic Research Fund of the Ministry of Education,
Singapore, under Grant MID-ARC 3/97 and DSO National Laboratories, Singapore, under Grand DSOCL01144.
F. Liao and J. L. Wang are with the School of Electrical and Electronic Engineering, Nanyang Technological University, Singapore 639798 (e-mail: ejlwang@ntu.edu.sg).
G.-H. Yang is with the Temasek Laboratories, National University of Singapore, Singapore 119260.
Publisher Item Identifier S 1063-6536(02)00339-1.
77
(5)
vertices, of which, one vertex corresponds
Here, there are
to the normal case (no fault), and the remaining vertices correspond to the fault cases. Uncertain system with polytopic uncertainty has been studied by many researchers [9], [11], [12]. In
this paper, we use it to model control surface impairment faults.
are known
The matrices , , , and
constant matrices of appropriate dimensions, which represent
the vertices of possible control surface impairment (including
the nominal case of no control surface impairment). Without
corresponds to the nomloss of generality, suppose that
,
,
,
.
inal case, i.e.,
Note that there is a certain degree of over-design in adopting
this polytopic uncertainty model. However, this ensures that the
designed performance can be achieved for intermediate control
surface impairment faults (as will be demonstrated by the nonlinear simulation results in Section V-B).
Hence, the aircraft dynamics with both actuator faults (2) and
control surface impairment (4) is described by
(2)
(6)
or
(3)
for
, the fault model
Obviously when
(2) corresponds to the case of the th actuator outage. When
, it corresponds to the case of no fault in the th actuator.
, namely,
Without loss of generality, we assume that
corresponds to the normal control input vector
.
The usual control surface fault is the control surface impairment. Unlike actuator faults of the form (2), aircraft control surface impairment will change the aerodynamic characteristics of
the aircraft (i.e., , , , and matrices). Control surface impairment can be characterized by the percentage loss of the total
control surface area. Models of the aircraft can be obtained for
a number of percentage losses of certain control surfaces. By
using these models, aircraft models corresponding to intermediate percentage losses of such control surfaces can be obtained
approximately by linear interpolation. To ensure that the design
method can handle the inherent modeling errors in these interpolated models, we adopt the following polytopic uncertainties
in the matrices , , , and :
(4)
78
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 10, NO. 1, JANUARY 2002
Fig. 1. Block diagram of fight control system with independent control surfaces.
It is easy to see from (4) and (10) that this augmented system
can be expressed in the polytopic form as follows:
(11)
where
(8)
(12)
Fig. 1 shows the block diagram of the flight control system. De,
fine the augmented state vector
and disturbance vector
,
. The augmented
system (8) can be rewritten as
(9)
where
, and
(13)
(10)
and
are symmetric positive
where
is symmetric positive definite.
semidefinite and
of the
Our objective is to minimize the performance index
nominal system (1), while simultaneously guarantees that the
) meet certain upper
performance indexes
.
bounds for all probable actuator faults and for all
79
(21)
are satisfied. By Schur complement, (21) is equivalent to
(14)
where
system is given by, for
(15)
Adopting the multiconvexity concept in [2], the following
Lemma 1 is easily derived.
Lemma 1: Assume that the hypotheses [H1] and [H2] hold.
(
;
For every given upper bound
), if there exist a matrix
and
with
, satissymmetric matrices
and
fying, for all
(22)
,
Since
(22) we have that
and
, from
(23)
affinely quadratically stabilizes
Hence controller
the augmented system (9). Furthermore, substituting
into (13) yields, for all
(16)
(17)
where
(18)
then the closed-loop system (15) achieves the upper bounds of
performance indexes
(19)
where
is as in (5) and
(20)
80
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 10, NO. 1, JANUARY 2002
(24)
(
;
subject to (16), (17), and
) as well as
(
;
;
) where
are performance bounds
in the fault cases.
However, the matrix inequalities (16) and (17) in Lemma 1
are not jointly convex. In order to get around this problem and
reduce the conservativeness, an iterative method is adopted here.
Theorem 1: Consider the closed-loop system (15). Assume
that the hypotheses [H1] and [H2] hold. For given upper
,
and a
bounds
given sufficiently small positive constant as well as a given
and given symmetric positive-definite
matrix
, suppose that there exist a
matrices
and symmetric matrices
feedback gain matrix
satisfying, for all
and
,
where ,
,
and
are as in (18). Then, the
controller (14) robustly stabilizes the closed-loop system of the
augmented system (9). Furthermore, the upper bounds of perare given by
formance indexes
(27)
, and
and
as in (20).
Proof: By the Schur complement formula, (25) and (26)
are equivalent to
with
(28)
and
(29)
(25)
(26)
with
(30)
(31)
Let
2) Minimize
and
81
) subject to
and
(32)
(33)
m/s
.
then we obtain
and performance
3) Choose a small error tolerance
(
;
bounds in the fault cases
;
). Let
(
;
) and
.
), minimize
4) At the th iteration (
subject to the LMI constraints
and
for
and
, and
for
and
,
. Here
and
are defined by (25) and (26), respectively.
, set
and stop.
5) If
and
for
Otherwise, let
and
, and go back to
Step 4.
,
,
It should be noted that the existence of solutions
and
in (25) and (26) is equivalent to the existence of
and
in (16) and (17). Therefore, Theorem 1
solutions
does not introduce any additional conservativeness other than
and
, (25) and
those in Lemma 1. However, for given
, which can be solved by using the
(26) are LMIs for and
LMI Toolbox in MATLAB environment. Hence, in this sense,
Theorem 1 converts a nonconvex problem into a convex problem
tends to
and its conservativeness can be minimized when
(
;
) and
tends to .
(
;
In Algorithm 1, the initial values of
) in (25) and (26) are obtained by using the
common LMI solution.
be defined in Step 4
Theorem 2: Let
is converof Algorithm 1. Then the sequence
gent.
,
is a feasible soProof: For a given integer
) in step 4, and
is
lution for the th optimization (
.
the optimal solution. So
must be convergent.
Hence, the sequence
Theorem 2 shows that Algorithm 1 is convergent. It should
also be noted that Algorithm 1 is guaranteed to achieve a local
optimum only. There is no guarantee for achieving the global
optimum.
V. APPLICATION EXAMPLES
In modern tactical fighter aircraft design, the control system
must be able to perform challenging maneuvers such as high
flight under a high
angle of attack flight and high roll rate
low speed condition. Under these conditions, there are significant kinematics and inertial couplings, so the variables stability
is the total
where is the engine throttle, is the altitude,
is the reference genter of gravity location, and
airspeed,
and
are, respectively, the horizontal stabilizer, the
,
and
aileron, and the rudder. Let
where
82
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 10, NO. 1, JANUARY 2002
TABLE I
PERFORMANCE INDICES COMPARISON FOR FAULTS 18 (! ! ) WITH
1 = 0:01
83
Fig. 2. Response curves of the nominal and fault 18 cases with Standard tracker
,
represents the vertical gust
(36)
, and , , , and
(
) being as
with
given in the Appendix.
By using the iterative algorithm, a reliable tracking controller
is designed to tolerate the above surface impairment faults in the
K.
84
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 10, NO. 1, JANUARY 2002
Fig. 3.
Response curves of the nominal and fault 18 cases with Reliable tracker
TABLE II
PERFORMANCE COMPARISON THE HORIZONTAL STABILATOR IMPAIRMENT
WITH
:
1 = 0 01
command for , , and are all set as one. It can be seen that
the response in Fig. 4 using the standard controller becomes
very oscillatory and is almost unstable when a 75% loss of the
horizontal stabilator occurs. But for our reliable controller, response remains well damped.
To verify the robustness of our reliable tracking controller,
simulation results using the original nonlinear aircraft model
and data are carried out, and the results are given in Figs. 6 and
7. In these simulations, a 25% loss of the horizontal stabilator
is introduced at 12 s, followed by a 50% loss at 16 s and a 75%
Fig. 4. Response curves of the nominal and horizontal stabilator impairment cases with the standard tracker.
Fig. 5. Response curves of the nominal and horizontal stabilator impairment cases with the reliable tracker.
Fig. 6. Simulation result for the standard tracker using nonlinear aircraft model.
85
86
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 10, NO. 1, JANUARY 2002
Fig. 7. Simulation result for the reliable tracker using nonlinear aircraft model.
APPENDIX
System Matrices in Example 1
m
m/s
87
m
m/s
88
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 10, NO. 1, JANUARY 2002
REFERENCES
[1] A. K. Caglayan, S. M. Allen, and K. Wehmuller, Evaluation of a second
generation reconfiguration strategy for aircraft flight control systems
subjected to actuator failure/surface damage, in Proc. IEEE 1988 Nat.
Aerospace Electron. Conf. (NAECON 1988), vol. 2, 1988, pp. 520529.
[2] P. Gahinet, P. Apkarian, and M. Chilali, Affine parameter-dependent
Lyapunov functions and real parametric uncertainty, IEEE Trans. Automat. Contr., vol. 41, pp. 436442, Mar. 1996.
[3] P. Gahinet, A. Nemirovskii, A. J. Laub, and M. Chilali, LMI Control
Toolbox. Natick, MA: The MathWorks, Inc., 1994.
[4] I. M. Horowitz, Synthesis Feedback Systems. New York: Academic,
1963.
[5] H. Kwakernaak, A polynomial approach to H optimization of control
systems, in NATO Workshop on Modeling, Robustness and Sensitivity
Reduction in Control Systems, Groningen, Dec. 1986.
[6] J. Lam and Y. Y. Cao, Simultaneous linear-quadratic optimal control
design via static output feedback, Int. J. Robust Nonlinear Contr., vol.
9, pp. 551558, 1999.
[7] J. W. Lee, J. H. Lee, and W. H. Kwon, Robust tracking control for uncertain linear systems without matching conditions, in Proc. 33rd IEEE
Conf. Decision Contr., Lake Buena Vista, FL, 1994, pp. 41814186.
[8] F. Liao, J. L. Wang, and G.H. Yang, LMI-based reliable robust preview
tracking control against actuator faults, in Proc. Amer. Contr. Conf.,
Arlington, VA, 2001, pp. 10471052.
[9] T. Mori and H. Kokame, A parameter-dependent Lyapunov function
for a polytope of matrices, IEEE Trans. Automat. Contr., vol. 45, pp.
15161519, Aug. 2000.
[10] L. T. Nguyen and M. E. Ogburn et al., Simulator study of stall/poststall characteristics of a fighter airplane with relaxed longitudinal static
stability, NASA Tech. Paper 1538, 1979.
[11] T. Pancake, M. Corless, and M. Brockman, Analysis and control of
polytopic uncertain/nonlinear systems in the presence of bounded disturbance inputs, in Proc. Amer. Contr. Conf., Chicago, IL, June 2000,
pp. 159163.
[12] U. Shaked, Improved LMI representations for the analysis and the design of continuous-time systems with polytopic type uncertainty, IEEE
Trans. Automat. Contr., vol. 46, pp. 652656, Apr. 2001.
Fang Liao received the B.E. degree in automatic control and navigation in 1992
and the M.E. degree in flight control, guidance, and simulation in 1995, both
from Beijing University of Aeronautics and Astronautics (BUAA). Currently,
she is pursuing the Ph.D. degree in the School of Electrical and Electronic Engineering at Nanyang Technological University.
From 1995 to 1999, she was an Engineer with the Research Institute of Unmanned Air Vehicles at BUAA. Her research interest includes robust and reliable control and tracking as well as genetic algorithms and their applications to
aircraft flight control system design.
Jian Liang Wang (M91SM98) received the B.E. degree in electrical engineering from Beijing Institute of Technology, China, in 1982. He received the
M.S.E. and Ph.D. degrees in electrical engineering from The Johns Hopkins
University, Baltimore, MD, in 1985 and 1988, respectively.
From 1988 to 1990, he was a Lecturer with the Department of Automatic
Control at Beijing University of Aeronautics and Astronautics, China. Since
1990, he has been with the School of Electrical and Electronic Engineering at
Nanyang Technological University, Singapore, where he is currently an Associate Professor. His current research interest includes robust and reliable control
and filtering, nonlinear control, two-time-scale systems, and their application to
flight control system design.
Dr. Wang was the Technical Program Chairman for International Conference
on Control, Automation, Robotics and Vision, Singapore, December 2000.
Guang-Hong Yang was born in Jilin, China, on September 2, 1963. He received the B.S. and M.S. degrees in mathematics from Northeast University of
Technology, China in 1983 and 1986, respectively, and the Ph.D. degree in control engineering from Northeastern University (formerly Northeast University
of Technology) in 1994.
From 1986 to 1995, he was with the Northeastern University. From 1996
to July 2001, he was a Postdoctoral/Research Fellow with the School of Electric and Electronic Engineering, Nanyang Technological University, Singapore.
From August 2001, he joint the Tamasek Laboratories at National University of
Singapore as a Research Scientist. His research interests include decentralized
control, symmetric systems, fault-tolerant control, nonlinear control, and nonfragile controller design with applications to aircraft flight control.
89