www.ietdl.org
is still in progress, much work should be done. For NFTC,
the faults are treated as additive actuator faults [15, 16]at
first, and the non-linear regulation theory is used to solve the
NFTC problem; a Lyapunov reconstruction technique, which
based on the a priori knowledge of a stabilising feedback for
the nominal safe model and the knowledge of the associ-
ated Lyapunov function, is used to solve NFTC problem in
[17]. In [18], a model predictive control (MPC)-based online
reference reshaping and controller reconstruction method
is presented. The proposed method has two main stages.
Firstly, the reference command of the faulty system is
reshaped online with respect to system faults; secondly,
based on a non-linear MPC strategy, the control of the plant
is reallocated according to the new reference command.
Although the mentioned method needs real-time calculation,
it provides an efficient way to the NFTC.
Motivated by the above discussions, in this paper, an
NFTC strategy for AHVs will be proposed. After presenting
the non-linear dynamics of AHVs, the fault model and the
control objective of the paper are discussed. Since the fault
model considered in this paper is a general form, and the
dynamics of AHVs are really complex, we propose an MPC-
based FTC strategy. The reference command is reshaped
with respect to the faults, firstly. Then an optimal prob-
lem is obtained by MPC, and through solving the optimal
problem online, the input of the plant can be obtained in
real time. A simulation analysis is provided to confirm the
effectiveness of the proposed control design approach. From
the simulation results we can see that the proposed control
strategy can guarantee a good tracking performance in the
existence of faults. The main contributions of the paper can
be summarised as follows:
(1) A non-linear FTC method is proposed for AHVs. Via
reshaping the reference command online and reconstructing
the input of the plant real timely, the FTC of AHVs can be
solved.
(2) The proposed method can deal with not only the FTC
problem, but also the input saturation of non-linear system.
(3) Though need real-time calculation, the proposed method
is easily to be carried out in practice.
The rest of this paper is organised as follows. The non-
linear model of AHVs, the fault model and the control
objective of this paper are presented in Section 2. Section 3
gives the main results of the paper, the reference reshaping
method and the MPC-based FTC strategy. In Section 4, a
numeric simulation is given and we conclude this paper in
Section 5.
Notation: The notations used throughout the paper are fairly
standard. The superscript ‘T’ denotes matrix transposition;
R
n
denotes the n-dimensional Euclidean space; and diag{...}
stands for a block-diagonal matrix. Matrices are assumed to
be compatible for algebraic operations if their dimensions
are not explicitly stated.
2 Problem formulation
2.1 Non-linear model of AHVs
The rigid-body equation of motion for AHVs considered in
this paper is developed by NASA Langley Research Center
[19]. The non-linear equations are described as follows
⎧
⎪
⎪
⎪
⎪
⎪
⎨
⎪
⎪
⎪
⎪
⎪
⎩
˙
V = (T cos α − D)/m − μ sin γ/r
2
˙γ = (L + T sin α)/mV − (μ − V
2
r) cos γ/Vr
2
˙
h = V sin γ
˙α = q −˙γ
˙q = M
yy
/I
yy
(1)
and the engine dynamics can be written as a second-order
system
¨
β =−2ξω
n
˙
β − ω
n
2
β + ω
n
2
β
c
(2)
where h and V represent the flight altitude and velocity of
AHVs, respectively; α is the angle of attack of the vehicle,
γ is the flight path angle and q represents the pitch rate. T ,
L, D and M
yy
are the thrust, lift, drag and pitching moment,
respectively. I
yy
is the moment of inertia. Equation (2) rep-
resents the dynamics of the actuator, where ξ is the damping
ratio of the actuator dynamics, ω
n
is the natural frequency
and β
c
is the throttle setting.
The expressions of L, D, T and M
yy
are modelled as
⎧
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎨
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎩
L =
1
2
ρV
2
SC
L
D =
1
2
ρV
2
SC
D
T =
1
2
ρV
2
SC
T
M
yy
=
1
2
ρV
2
S ¯c[C
M
(α) + C
M
(δ
e
) + C
M
(q)]
r = h + R
E
(3)
where
⎧
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎨
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎩
C
L
= 0.6203α
C
D
= 0.6450α
2
+ 0.0043378α + 0.003772
C
T
=
0.02576β, β<1
0.0224 + 0.00336β, β>1
C
M
(α) =−0.035α
2
+ 0.036617α + 5.3216 × 10
−6
C
M
(δ
e
) = ce(δ
e
− α)
C
M
(q) = (¯c/2V )q(−6.796α
2
+ 0.3015α − 0.2289)
(4)
¯c represents mean aerodynamic chord and δ
e
means elevator
deflection of AHVs. This non-linear model is composed of
five rigid-body state variables x =[V , γ , α, β, h]
T
, the output
to be controlled is selected as the velocity V and the altitude
h, then y =[V , h]
T
. The control input u =[β
c
, δ
e
]
T
does not
appear explicitly in the equations. For a simple description,
the following equations are used to represent the non-linear
dynamics of AHVs in this paper
˙x = f (x, , u)
y = h(x)
where is a family of parameters included in (2)–(4).
2.2 FTC objective
Faults are inevitable, and may influence stability of AHVs.
Faults may locate in the actuators, the components or the
elevators of AHVs and may be caused by partial damage of
1148 IET Control Theory Appl., 2014, Vol. 8, Iss. 13, pp. 1147–1153
© The Institution of Engineering and Technology 2014 doi: 10.1049/iet-cta.2013.0986