MA Xinjian, LIU Shiqian, and CHENG Huihui
School of Aeronautics and Astronautics, Shanghai Jiao Tong University, Shanghai 200240, China
Abstract: For the problem of sensor faults and actuator faults in aircraft attitude control, this paper proposes a fault tolerant control (FTC)scheme based on extended state observer (ESO)and nonlinear dynamic inversion (NDI).First, two ESOs are designed to estimate sensor faults and actuator faults respectively.Second, the angular rate signal is reconstructed according to the estimation of sensor faults.Third, in angular rate loop, NDI is designed based on reconstruction of angular rate signals and estimation of actuator faults.The FTC scheme proposed in this paper is testified through numerical simulations.The results show that it is feasible and has good fault tolerant ability.
Keywords: fault tolerant control (FTC), signal reconstruction,extended state observer (ESO), nonlinear dynamic inversion(NDI).
To ensure flight safety of civil aircrafts, fault tolerant control (FTC)are designed to handle faults of important sensors and actuators.Typical researches on actuator FTC are conducted [1-9].Alwi et al.used a sliding mode observer (SMO)to estimate actuator fault [1,2], and Quan applied this method on flight control [3], but this method was only suitable for systems with linear observable terms.Lu et al.proposed an incremental nonlinear dynamic inversion (NDI)control to realize fault tolerant trajectory tracking [4], and Cheng et al.combined incremental NDI with neural network to deal with actuator fault and pilot mishandling simultaneously [5].However,this method had heavy work load on parameter adjustment.Lee et al.used an unknown input observer to detect actuator fault [6], but this method was also based on the linearized system.Yu et al.combined fuzzy neural networks and sliding mode controller to realize the FTC [7], but this method involved complex calculation.Wang et al.used stable inversion on fault tolerant auto landing of civil aircrafts [8], but it did not solve the problem of nonlinearity.Liu et al.proposed a sliding mode backstepping to do trajectory FTC [9] based on a nonlinear aircraft.In order to deal with the nonlinearity, recently, researchers have been working on actuator FTC by using extended state observer (ESO)[10-17].On behalf of these researches, Yuan et al.combined ESO with NDI to improve the control robustness [16], and Falconi et al.further realized FTC through this method [17].Theoretical analysis of ESO were given in [18,19].These ESO-based FTC methods were designed to solve the problem of nonlinearity.
However, few of these researches mentioned above involved the sensor faults.Hence, in order to solve the problems of sensor faults, Edwards et al.proposed an SMO-based fault detection method for sensor faults [20],and Xia et al.proposed an observer-based signal reconstruction method [21], but they were only suitable for the linear system.Johansen et al.proposed an airflow angle estimation method which can be used when airflow angle sensor have faults [22], but its accuracy still needs to be improved.Li et al.used Kalman filter to do the signal estimation under white noise [23], but this method has difficulty dealing with bias sensor fault.Chen et al.proposed a nonlinear adaptive observer-based sensor FTC scheme [24].Yang et al.proposed a fault diagnosis method for angle-of-attack based on the transfer function of linear model [25], but this method has difficulty dealing with nonlinear systems.Lu et al.proposed a sensor fault detection method in [26,27] and developed an FTC scheme for simultaneous fault of sensors and actuators by decoupling design [28].However, Lu’s method needs sophisticated parameter adjustment.Generally, most of these methods either have difficulty dealing with nonlinear systems or need complex calculation.Therefore, how to solve simultaneous sensor and actuator faults in aneasy way has become a heated topic.
The control scheme of [16] and [17] based on NDI and ESO has a good performance on actuator FTC.However,neither of them has taken sensor faults into consideration.Thus, motivated by the work of [16] and [17], this paper further improves the control law design so that it can deal with actuator faults and sensor faults simultaneously.First, two ESOs are designed to estimate sensor faults and actuator faults respectively.Second, the angular rate signal is reconstructed according to the estimation of sensor faults.Third, the NDI in angular rate loop is designed based on the reconstruction of angular rate signals and the estimation of actuator faults.Finally, another NDI is designed to realize the attitude tracking.The proposed FTC scheme is verified through simulations.The results show that it is feasible and has better fault tolerant ability than conventional ESO-NDI design.
Main innovation and contribution of this paper can be summarized as follows:
(i)To solve the problem of sensor faults, this paper proposes a novel ESO-based sensor fault estimation method and realizes the reconstruction of aircraft angular rate signal using a nonlinear aircraft model.
(ii)Conventional ESO-NDI design has difficulty dealing with simultaneous sensor and actuator faults.To overcome this weakness, the FTC is modified, and the fault tolerant attitude control of an aircraft is realized with simultaneous sensor and actuator faults.
This paper is organized as follows: Section 2 gives the aircraft model; Section 3 gives the design of ESO and NDI; Section 4 gives some simulations on a Cessna172 aircraft model; Section 5 presents conclusions.
The aircraft attitude control scheme includes an attitude angle control loop and an angular rate control loop, as shown in Fig.1, whereωandηdenote the angular rate and attitude vector respectively, the associated commands areωcandηc, and measurements areωmandηmrespectively.ucanduare command and actual inputs of actuators.
Fig.1 Attitude control scheme
According to [29], the angular rate dynamics of aircraft is presented as follows:
wherep,q, andrrepresent the roll, pitch, and yaw rate of the aircraft respectively.L,M, andNrepresent the roll,pitch, and yaw moments of the aircraft.Ixx,Iyy,Izz, andIxzare the moments of inertia of the aircraft.Denoteas the state vector of angular rate loop.
Assume that the moments on the aircraft are only generated by aerodynamic force, which means the moments generated by thrust are ignored.Hence, the aerodynamic moments can be calculated as follows:
whereRbwis the coordinate transform matrix from wind frame to body frame.The transform matrix can be found in [29].is the dynamic pressure,Sis wing reference area,bis wing span,is the mean aerodynamic chord.CL,CM, andCNare aerodynamic moment coefficients,CD,CY, andClare aerodynamic coefficients of drag, lateral and lift force in wind frame respectively.And Δlis the distance vector from pressure center to gravity center that can be described as follows:
whereRwbis the transform matrix from body frame to wind frame, andxcgandxcpare the position of gravity center and pressure center alongxaxis in body frame.
The aerodynamic moment coefficientCL,CM, andCNare coefficients that directly affect angular motion, which can be calculated as follows:
whereCL0,CM0, andCN0are the constant parts of aerodynamic moment coefficient.CLβ,CMα, andCNβare aerodynamic moment coefficients generated by airflow angles.Vtis the true airspeed.CLp,CNp,CMq,CLr, andCNrare aerodynamic moment derivatives.CLδa,CNδa,CNδe,CLδr, andCNδrare control surface related derivatives, δa, δe, andδrare the deflection of aileron, elevator, and rudder, respectively.Denoteas the input to aircraft.Substituting (2)and (4)into (1)yields
whereF1(ω)andG1(ω)are nonlinear functions that can be described as follows:
F01,MI,Mω, andMuare described as follows:
Since sensor and actuator faults may occur in the angular rate control loop, the angular rate loop can be described as follows:
wherefais actuator fault, which means the difference between actual deflection and desired deflection of aileron, elevator and rudder.fsis sensor fault, which means the difference between angular rate sensor measurements and true value.ucis the command of control input,uis the eventual input (u=uc+fa), and ωmis measurement value of angular rate.
The attitude control loop uses Euler angles to describe the attitude.Denoteas the state vector for the attitude loop, where φ, θ, and ψ are the roll, pitch,and yaw angle respectively.Attitude dynamics can be described as (13)according to [29]
where
Assume that there are no sensor faults in attitude variables.Hence, the attitude loop can be described as
where ηmis the measurement of the attitude variables.
The FTC design includes four parts.ESOsis for sensor fault estimation, while ESOais for actuator fault estimation.NDIωis the controller of angular rate loop, while NDIηis the controller of attitude loop, as shown in Fig.2.
Fig.2 FTC design based on ESO and NDI
ESOsis used to estimate the sensor faults and reconstruct the angular rate signals.Substitute (12)into (15), and the attitude dynamics can be rewritten as follows:
Denote (16)as System A.Decompose System A into three channels:
whereG2i(η)are theith row ofG2(η)(i= 1, 2, 3).
Denote:
wheredφ,dθ,dψdenote the influence of faults in channels of roll, pitch, and yaw respectively.
Design an equivalent system to System A as (19), and denote it as System B
According to (18), sensor faults can be estimated as
Finally, angular rate signal can be reconstructed as
ESOais designed to estimate actuator faults.Denote (12)as System C, and decompose System C into three channels:
whereF1i(ω)andG1i(ω)are theith row ofF1(ω)andG1(ω)(i= 1, 2, 3).Denote:
wheredp,dq,drare the influence of actuator faults in roll,pitch, and yaw rate channel respectively.Then System D can be designed as an equivalent system to System C:
NDIωis for angular rate control.The design of NDIωis based on the signal reconstruction and fault estimation.It can be described in the following steps:
Step 1Calculate tracking error of angular rates.
Step 2Calculate desired derivatives of angular rates according to proportion-integral-derivative (PID)rules.
wherekp1,ki1, andkd1are PID gains.
Step 3Let the desired derivatives equal to the true derivatives.
Step 4Substitute (12)into (30).
Step 5Transform (31), and replace ω andfawithandrespectively.
NDIηis for attitude control.The design of NDIηcan be described as following steps:
Step 1Calculate tracking error of attitude angles.
Step 2Calculate desired derivatives of attitude angles according to PID rules.
wherekp2,ki2, andkd2are PID gains.
Step 3Let desired derivatives equal to true derivatives.
Step 4Substitute (13)into (35).
Step 5Transform (36), and replace ω with ωc.
The simulations are carried out on a Cessna172 aircraft,cruising at the height of 500 m and at the speed of 50 m/s.The aircraft model is established based on (12)and (13),with reference to [29] and [30].
In Scenario A, the ESO-NDI-based FTC scheme is compared with normal NDI to verify the feasibility and fault tolerant ability of ESO-NDI-based FTC proposed in this paper.The faults are set as (38)to (40).
The control parameters are set as follows:
Angular rate responses under simultaneous sensor and actuator faults are shown in the following:
Fig.3 gives the comparison of angular rate tracking results based on the FTC proposed in this paper and conventional ESO-NDI design in [16] or [17].At the beginning the tracking results of proposed FTC have error in pitch rate channel due to the initial value error of ESOs,but the error diminishes in two seconds.At 5 s, actuator faults cause tracking error in pitch rate channel, and both methods have eliminated the influence of actuator fault.After 10 s, the sensor fault occurred.It can be seen that conventional ESO-NDI design failed to follow the command, while the proposed FTC eliminate the influence of sensor faults and follows the command well.The comparison shows that conventional ESO-NDI design cannot deal with sensor faults, while the proposed FTC scheme can deal the simultaneous sensor and actuator faults,which means it has better fault tolerant ability.
Fig.3 Angular rate tracking responses
Fig.4 gives the angular rate signal reconstruction results.At the beginning, the initial value error caused some error, but it soon diminished.At 10 s and 20 s, the measurement suddenly changed and caused reconstruction error, but the reconstruction soon get close to the true value, which shows the effectiveness of the signal reconstruction.
Fig.4 Angular rate signal reconstruction results
In Scenario B, the attitude tracking response is given.The faults are set as (44)to (48).
The parameters of ESOa, ESOs, and NDIωare set as(41)to (43).The parameters of NDIηis given as follows:
During 0-25 s in angular rate loop, conventional ESONDI design is used and after 25 s the proposed FTC is used.The simulation response are shown in the following:
Fig.5 Attitude tracking response
Fig.6 Angular rate tracking response
It can be seen that sudden actuator faults at 5 s cause little tracking errors, but soon are eliminated.During 10-25 s, sensor faults cause tracking error in angular rate loop and attitude loop, and huge attitude deviation is caused, which again shows conventional ESO-NDI design has difficulty in dealing with sensor faults.After 25 s, the proposed FTC was used.As a result, angular rates and Euler angles follows the command successfully,which proves the proposed FTC has the ability to deal with simultaneous sensor and actuator faults.
In this paper, an FTC scheme based on ESO and NDI for civil aircrafts with sensor and actuator faults is proposed.First, an ESO-based signal reconstruction is proposed to solve the sensor faults.Second, another ESO is designed to estimate actuator faults.The angular rate controller and attitude controller are designed by using NDI.In angular rate loop, the NDI is designed based on the signal reconstruction and the estimation of actuator faults.Simulation shows that this FTC scheme is feasible and has better fault tolerant ability than conventional ESO-NDI design in dealing with sensor faults.The future work will be focused on faults of different kinds of sensors and how to deal with sensor noises.
Journal of Systems Engineering and Electronics2022年1期