Actuator Fault Diagnosis in a Boeing 747 Model via Adaptive Modified Two-Stage Kalman Filter

1 Department of Control Engineering, Istanbul Technical University, Maslak, 34469 Istanbul, Turkey 2Department of Mechanical & Industrial Engineering, Concordia University, 1455 de Maisonneuve Boulevard W., Montreal, QC, Canada H3G 1M8 3Department of Electrical and Computer Engineering, Binghamton University (SUNY), Binghamton, NY 13902-6000, USA 4Gulfstream Aerospace Corporation, Savannah, GA 31402-2206, USA


Introduction
The Integrated Resilient Aircraft Control (IRAC) Project by NASA has conducted research to advance the state of aircraft flight control to provide on-board control resilience for ensuring safe flight in the presence of unforeseen, adverse conditions.The focus of the IRAC project research is to develop a set of validated multidisciplinary integrated aircraft control design tools and techniques for enabling safe flight in the presence of adverse conditions such as structural damage, control surface failures, icing, or aerodynamic upsets.
Aircraft have multiple control surfaces for high performance as well as redundancy.When one or more control surfaces have problems during the operation, a Fault Detection and Identification (FDI) scheme should be carried out to continue a safe flight.Figure 1 shows some basic control surfaces in aircraft.FDI methods based on analytic redundancy have been actively studied using various estimation theories.There are a number of model-based FDI methods that utilize observers, parity equations, or parameter estimation techniques.
This paper reports some results on the simultaneous estimation of the states and parameterized faults injected into the actuators of a nonlinear model of the longitudinal motion of a Boeing 747 [1], modifying the linear adaptive two-stage Kalman filter in [2].Fault parameterization in the design model of the linear estimator aims to capture the degree at which a control surface is stuck and the extent at which the control effectiveness is lost.The estimates serve to support real-time assessment of postfailure flight envelope, in which a controlled aircraft can maintain stability and safe flight [3].
Among very rare studies, Vachtsevanos et al. [4] considered the control role portrayed by mode transitioning, envelope protection, real-time adaptation, and fault detection/control reconfiguration algorithms which are intended to safeguard unmanned aerial vehicle's integrity in the event of component failures, extreme operating conditions, or external disturbances.A fault-tolerant control system architecture was designed to combat faults in the flight control actuators.Flight test results demonstrated that the fault-tolerant architecture can accommodate stuck actuator malfunctions.
Chen and Jiang [5] investigated a linear fault-tolerant control system design technique by using an iterative learning observer when stuck actuator faults had occurred.This observer estimates both state information and stuck actuator values.The technique can accommodate multiple stuck actuators provided that enough redundant actuators exist.Boskovic and Mehra [6] presented an FDI and Adaptive Reconfigurable Control (ARC) scheme for accommodation of control effector failures such as lock-in-place and hardover.The overall system consisting of online FDI observers for all control effectors guarantees asymptotic tracking and ensures convergence of the failure-related parameter estimate to its true value.The scheme can rapidly and accurately identify different control effector failures like float, hard-over, lock-in-place, and loss of effectiveness.They used adaptive interacting multiple observers to estimate the actuator faults and tested the scheme on a linearized model of the Boeing's tailless advanced fighter aircraft.
Jiang and Chowdhury [7] developed a real-time fault estimation module to estimate the actuator effectiveness.They transformed the system into two separate subsystems.One subsystem is not affected by actuator faults.The states were estimated by a reduced order Kalman filter.They demonstrated the simulation results on a linear helicopter model.
Yang and Wang [8] studied an iterative Linear Matrix Inequality-(LMI-) based stuck actuator fault detection method without capability of isolating and identifying the faults.The considered systems are modeled via multiple modes, namely, nominal case and faulty case.The actuator stuck faults are considered.A new LMI-based fault detection method for the control systems is obtained and applied to the fault detection for flight control systems.
Wang and Lum [9] developed an adaptive unknown input observer approach to detect, isolate, and identify the loss of effectiveness and stuck faults in the actuators.The approach is useful for linear systems, but the attempt to extension to nonlinear systems has not been successful [10].
The two-stage Kalman filter of Keller and Darouach [11] was applied to estimate simultaneously the state and the control effectiveness of a linear aircraft model in [2], which in addition introduced covariance-dependent forgetting factors into the filtering algorithm to accelerate convergence of the fault parameter estimates [12].Thus, the filter acquired the name adaptive two-stage Kalman filter (TSKF).The state and control effectiveness estimates with the adaptive TSKF approach were utilized to help achieve fault-tolerant control of impaired linear aircraft model in a number of papers by the researchers such as [13,14].Amoozgar et al. [15] addressed the problem of Fault Detection and Diagnosis (FDD) of a quadrotor helicopter system in the presence of actuator faults.A TSKF is used to simultaneously estimate and isolate possible faults in each actuator.The faults are modelled as losses in control effectiveness of rotors.The developed FDD algorithm is evaluated through experimental application to an unmanned quadrotor helicopter test bed.The obtained results show the effectiveness of the proposed FDD method.
Hajiyev [16] presented an active fault-tolerant control technique against actuator/surface failures.A TSKF was designed in order to identify the control distribution matrix elements that correspond to the faulty actuator/surface; thus, the control reconfiguration was carried out using the identified control distribution matrix.The actuator/surface fault identification problem was solved through two jointly operating Kalman filters: the first one was for the estimation of the control distribution matrix elements that correspond to the faulty actuator/surface and the second one was for the estimation of the state variables of the aircraft model.In simulations, the nonlinear flight dynamics of an AFTI/F-16 fighter model was considered, and the performance of the proposed actuator/surface failure identification and reconfigurable control schemes were examined for this model.
Jiang et al. [17] presented a framework for actuator fault accommodation using a multiple model-based direct adaptive control approach that does not include explicit fault detection and isolation.The approach does not require exact information about the controlled system and persistent input excitations.Robustness is increased, and stable adaptation of unknown faults is provided.Asymptotic model following conditions and adaptive rules are derived and system stability is guaranteed, while appropriate switching of the multiple models ensures asymptotic tracking for system outputs.The efficiency of the proposed method is illustrated by simulations on an aircraft model.
This paper generalizes the design model for the linear adaptive TSKF of Wu et al. [2] which, from now on, will be referred to as a modified TSKF (MTSKF) to include stuck faults in control surfaces.The approach is applied to state and parameter estimation for a B747 nonlinear aircraft model.In order to distinguish the stuck faults from the loss of control effectiveness, persistent excitations must be present.Our findings indicate that the adaptive MTSKF can correctly estimate the aircraft states, as well as the parameters representing less severe faults near a trim point.The identification of the aerodynamic coefficients is not studied in this paper.
The paper is organized as follows.In Section 2, the linear adaptive MTSKF algorithm is described, and the linear TSKF algorithm of Wu et al. [2] is presented.Simulation results are given in Section 3. In conclusion, we summarize and give comments about the results.

Adaptive Modified Two-Stage Kalman Filtering (MTSKF) Algorithm
Based where   ∈   ,   ∈   , and  +1 ∈   are the state, control input, and output variables, respectively.  and   are bias vectors of dimension , representing faults entering actuators.
The noise sequences   ,   , and V are assumed to be zero mean uncorrelated white Gaussian noise sequences with where   > 0,   > 0,  > 0, and   is the Kronecker delta.The initial states (0) and (0) are assumed to be uncorrelated with the white noise processes   ,   , and V and have covariance P 0 and, covariance P 0 , respectively.The components of bias vector   describe percent reduction in control effectiveness when the terms     +     are considered together, where Estimator design model ( 1) is inherited from (3) of Wu et al. [2] with a set of new bias components added in this paper to denote the degrees at which control surfaces are stuck.A stuck fault is modelled by the combination of the following three terms (1): We now summarize the cases the fault parameter values represented.Specifically, no-fault case is represented by 100   % loss of control effectiveness in the th actuator is represented by control surface stuck at magnitude    degrees is represented by The linear adaptive TSKF of Wu et al. [2] can be directly applied to estimate both   and   via generalizing     to where which model the stuck fault as well, in addition to the loss of control effectiveness fault originally considered in Wu et al. [2].Since   =     , therefore, the entries of   = diag( 1  , . . .,    ) must vary and evolve in time independently (persistently excited) to allow the estimator to distinguish between   and   and among their components.
To make the two-stage filter algorithm more responsive to the changes of bias parameter, like the fault of parameters, following the idea in Parkum et al. [12], a bias covariance matrix related forgetting factor technique can be used in (16).Consider The forgetting factor    can be chosen as a decreasing function of the amount of information received in the direction    .Since the eigenvalue   / of   / is a measure of the uncertainty in the direction of    , a choice of forgetting factor    can be chosen as The linear TSKF algorithm of Wu et al. [2] where the filter residual vector and its covariance are given as

Simulation Results
Model.The adaptive MTSKF is applied to the longitudinal model of nonlinear B747 aircraft linearized around a trim point.The overall structure of the aircraft subject to actuator failures is given in Figure 2.
The model used in this study is a Boeing 747 aircraft, which is an intercontinental wide-body transport with four fan jet engines designed to operate from international airports.The nonlinear model of the aircraft given in [18], including the set of aerodynamic coefficients, is linearized around a trim point and then used in the simulations to assess the actuator fault isolation and identification performance of the adaptive MTSKF.The state space model of the linearized longitudinal dynamics of the B747 is as follows: where is angle of attack (rad),  is pitch rate (rad/sec),  is true air speed (m/sec),  is pitch angle (rad), ℎ is altitude (m),   is elevator deflection (deg),   is stabilizer deflection (rad), and   is engine thrust (N).Note that the states and the inputs are perturbations around the trim values.
Throughout the simulations, the symbols seen in the figures are as follows: the state vector is The data corresponding to the states and control inputs were collected from a Simulink model of the closed-loop nonlinear B747 model provided by Shin et al. [19].The data were used in the adaptive MTSKF algorithm coded in Matlab for a sampling time of 0.05 seconds.The same results were also obtained by online identification settings in Simulink.The actuator faults are assumed to occur symmetrically, for example, right and left stabilizers, or right and left elevators are considered to move together.
Several experiments were performed as follows.
Experiment 1 (stuck fault in stabilizer at 0.05 (rad) at 20 seconds).Under the fault condition, simulation results for state and stuck fault estimation are shown in Figures 3 and  4, respectively.
In Figure 3, it can be seen that the states are estimated correctly before and after occurrence of the stuck fault in the stabilizer.
From Figure 4, it can be concluded that a stuck actuator fault with magnitude 0.05 (rad) has occurred in the stabilizer because ζ2 = γ2 = −1 and ζ5 = β2 = 0.05.When γ becomes −1, we can obtain the estimated stuck magnitude in the corresponding control channel.
The estimation error  6 − ζ6 seems large, but, in fact, it is relatively small since the actual value for engine thrust is 41631 (N).Therefore, the estimation error in the engine thrust under the considered stuck magnitude is acceptable.Experiment 2 (50% partial loss in stabilizer at 20 seconds).Under the fault condition, simulation results for state and stuck fault estimation are shown in Figures 5 and 6, respectively.
In Figure 5, it can be seen that the states are estimated correctly before and after occurrence of the partial loss in the stabilizer.The operating point has not changed much.
Experiment 3 (stuck fault in elevator at 0.0001 (rad) at 20 seconds).Due to the significant effect of elevator on the dynamics of the B747, the operating point has considerably changed due to the stuck fault in elevator.The adaptive MTSKF could not estimate the accurate stuck magnitude due to the deviation of the linearized model used in the adaptive MTSKF although the fault magnitude is small.However, satisfactory state estimates are obtained as shown in Figure 7.It is obvious that the trim point has changed after 20 seconds.The elevator stuck at trim has also been tested but the operating point has considerably changed again and the plant has become unstable.The filter is not successful in identifying stuck elevator faults.Therefore, the controller should be redesigned to take care of instability.
The operating point has changed severely.Since actual and estimated state variables are too close, two curves are seen as one curve.Although state estimation is highly accurate, the parameter estimation is not accurate.Therefore, the simulation results are not presented.To overcome this problem, a general two-stage extended Kalman filter (GTSEKF), which extends the linear general two-stage Kalman filter to nonlinear systems, may be proposed.As, in the linear one, the GTSEKF is derived mainly by applying the nonlinear twostage transformation to the well-known extended Kalman filter (EKF).A nonlinear filter for estimating constant parameters in dynamic systems is presented in [20].
International Journal of Aerospace Engineering It has been observed from the simulations that actuator dynamics play very important role in the identification of the parameters.The excitation and control signals need to be prefiltered to be able to obtain accurate estimation results under the assumption that the adaptive MTSKF does not have actuator model.It has also been observed that larger   makes the convergence faster.

Conclusions
In this paper, an adaptive MTSKF was designed and used to estimate the reduction of control effectiveness and the magnitude of stuck faults for a nonlinear B747 model.Some simulation results on fault parameter estimation performed on a nonlinear model of the longitudinal motion of the aircraft were presented.The linear adaptive MTSKF is successful in identifying the magnitude of the stuck fault and the percentage of the partial loss in the stabilizer.Control inputs have been excited by a pseudorandom noise such that successful estimation of stuck magnitudes and loss of effectiveness can be achieved for a limited linear region.The study showed that the linear adaptive MTSKF is successful in identifying low degree stuck fault and partial loss in the stabilizer.However, the plant nonlinearity must be included in the design model of the filter to allow the parameters represent more severe faults to be correctly estimated.While this study focuses on Boeing 747 aircraft, the estimation procedure may also have applications for other aircraft types and to the military and space sectors.

Figure 3 : 6 InternationalFigure 4 :
Figure 3: Actual state variables and estimations when a stuck fault in stabilizer occurred at 20 seconds.Since actual and estimated state variables are too close, two curves are seen as one curve.

Figure 5 :
Figure 5: Actual state variables and estimations when 50% partial loss in stabilizer occurred at 20 seconds.Since actual and estimated state variables are too close, two curves are seen as one curve.

NomenclatureFDI:
Fault Detection and Identification ARC: Adaptive Reconfigurable Control LMI: Linear Matrix Inequality TSKF: Two-stage Kalman filter FDD: Fault Detection and Diagnosis
the following discrete time model is used as the design model of the adaptive MTSKF  +1 =     +     +     +     + on the linearized model of the open-loop aircraft around a trim point, and a parameterization of two types of International Journal of Aerospace Engineering 3 actuator faults, is as follows.Bias-free state estimator isx+1/ =   x/ +     + [  −  +1/ ] ζ/ , Modified Two-Stage Kalman filter   ,   , and  +1 : The state, control input, and output variables   and   : Bias vectors representing faults entering actuators   ,   : MTSKF: