Fault-Tolerant Attitude Control Incorporating Reconﬁguration Control Allocation for Supersonic Tailless Aircraft

: This paper presents a fault-tolerant attitude control scheme, incorporating reconﬁguration control allocation for supersonic tailless aircraft subject to nonlinear characteristics, actuator constraint, uncertainty, and actuator faults. The main idea is to propose an incremental reconﬁguration closed-loop control allocation scheme, coupled with a basic backstepping attitude controller, to achieve attitude control. Based on the virtual control input generated by the basic backstepping attitude controller, ﬁrstly, the incremental nonlinear control allocation method is adopted to deal with the nonlinear characteristics and actuator constraint. Secondly, a distribution error feedback loop is constructed in the incremental nonlinear control allocation method to enhance the robustness against the uncertainty of the control effectiveness matrix. Thirdly, the control effectiveness matrix is reconstructed by different kinds of fault information to deal with actuator faults, and the proper combination of actuator deﬂections is generated to achieve accurate command tracking. The stability of the proposed scheme is guaranteed by the Jury stability criterion and the Lyapunov stability analysis. Finally, in comparison with the three existing approaches, the simulation results of two cases are provided to show the effectiveness of the proposed scheme.


Introduction
Supersonic tailless aircraft are over-actuated systems equipped with redundant actuators [1].Distinguishing them from conventional fixed-wing aircraft, the innovative tailless layout optimizes flight safety and aerodynamic performance by increasing the number of actuators.However, due to the lack of vertical and horizontal tails, the relationship between actuators and flight control becomes less obvious, leading to less stability and greater coupling between the lateral and longitudinal dynamics [2].Additionally, actuator faults are more likely to happen with the increase in actuator amount.Therefore, the control design of supersonic tailless aircraft is a significant but extremely challenging problem.
Recently, the flight control of supersonic tailless aircraft has drawn great attention among a growing number of researchers, and many flight controllers are designed to solve these problems [2][3][4][5][6][7].In [3], a robust adaptive controller was designed based on online reinforcement learning control to guarantee the stability of the supersonic tailless aircraft.In [4], a sample entropy-based prescribed performance controller (SPPC) was proposed for the longitudinal control of a supersonic tailless aircraft subject to model uncertainty and nonlinearity.In [5], an incremental backstepping sliding-mode controller was designed for unknown disturbances and model uncertainties.All these methods have had a significant effect on improving the flight stability of supersonic tailless aircraft.
Moreover, fault-tolerant control (FTC) is applied to improve the safety and reliability of supersonic tailless aircraft against the actuator fault.Generally speaking, the FTC can be categorized into two types: passive FTC and active FTC.The passive FTC designs the controller based on the prior information of the actuator fault to make the system stable [2,[8][9][10].In [2] a fixed-time IT2 fuzzy fault-tolerant control scheme was proposed to deal with the uncertainty of the supersonic tailless aircraft including time-varying actuator failures.However, the controllers of passive FTC are usually designed with fixed structures and parameters, which decreases the flexibility of controllers.In active FTC, the controller is reconstructed based on the fault information to react to the actuator fault [11][12][13][14][15][16][17][18].Zhang et al. [18] proposed an incremental adaptive fault observer to estimate the failure information of tailless aircraft actuators.These control schemes all have excellent performance to deal with the actuator fault.
However, these methods result in control laws specifying the force and moments, rather than the deflection of actuators.In the majority of the approaches, the relationship between the actuator deflection and their generated moment is required to be obvious.Additionally, the controller design becomes complicated due to the inevitable adjustment of the basic control law.Hence, to distribute the virtual command, such as moment and force, to the available actuators, and reduce the difficulty of controller design, control allocation (CA) methods are proposed.
The CA method has been elaborated on and summarized in numerous studies [19][20][21][22].The principle of CA for aircraft can be divided into three categories: direct allocation [23,24], the optimization-based method [25,26], and dynamic control allocation [27][28][29].In direct allocation (DA), the boundary of the moment command is determined by the attainable moments set (AMS), and the optimal actuator command is obtained to reach the desired moment in a certain direction under the constraint of AMS.However, the computation is too complex to fit real-time applications for a large number of actuator situations.In the optimization-based method, CA is solved by minimizing the desired value function.In an unconstraint situation, the CA is described as a 2-norm minimization problem by the pseudo-inverse method.When considering the constraints, the CA problem can be solved by the redistributed weighted pseudo-inverse method or converted to a quadratic program.In the quadratic method, a secondary objective is considered in moment distribution, which can be used to achieve multi-objective optimization.Moreover, with the development of optimization algorithms, the advanced optimization algorithms, such as the butterfly optimization algorithm [30] and deep neural network [31], can be applied to the optimization-based CA method.In dynamic control allocation, the cases where the actuator dynamic is not much faster than the system dynamic is considered, which means that the actuator dynamic cannot be neglected.Dynamic control allocation was first purposed in [27] and further developed in two directions: quadratic programming [28], and model predictive control allocation [29].For an intact system, the redundant actuator is traded for efficiency maximization by the above control allocation method.However, the optimization-based method and dynamic control allocation mentioned above are the open-loop control allocation, relying on the accuracy of the aircraft model.With the actuator fault and uncertainty of the control effectiveness matrix, the accuracy of the aircraft model cannot be guaranteed.
For actuator fault and uncertainty, the control allocation method was developed in two directions: the closed-loop control allocation (CCA) method [15,32,33] and the reconfiguration control allocation [34][35][36].Similar to passive FTC, the CCA applies state feedback to enhance system robustness to compensate for the actuator fault and the uncertainty.However, the CCA has a fixed structure, which limits the fault-tolerant capability.Referring to the active FTC, the reconfiguration control allocation method is proposed.Different from the CCA method, the reconfiguration control allocation redistributes the virtual control command to the available actuator by reconstructing the control effectiveness matrix with the information of the fault actuator, taking full advantage of the available actuators.In [36], a reconfiguration control allocation (RCA) scheme was constructed to compensate for the actuator fault.All these methods are effective to compensate and actuator fault and uncertainty.However, these methods fail to consider the nonlinear relationship between the virtual control command and the deflection of actuators.
For the nonlinear characteristic, many nonlinear CA methods [37][38][39][40] have been proposed, including nonlinear direct allocation [37], affine control allocation with intercept correction [38,39], and nonlinear programming [40].Although all these methods can solve the nonlinear characteristic, they are not applicable to real-time systems because of the computational complexity.To overcome this limitation, the incremental nonlinear control allocation (INCA) method [41] was proposed.INCA transforms the CA problem into an incremental scheme, applying the linear CA method to deal with the nonlinear characteristics to reduce computational complexity.However, these methods fail to consider the actuator faults.
These methods have limited compensation for faults in supersonic tailless aircraft.Hence, the fault-tolerant CA issue for the over-actuated supersonic tailless aircraft with nonlinear characteristics is still an open problem, which has not been solved or at least well solved at present.
Hence, in this paper, a new fault-tolerant CA method scheme, an incremental reconfiguration closed-loop control allocation (IRCCA), is proposed.Combined with a basic backstepping attitude controller, the attitude control of supersonic tailless aircraft with actuator faults is achieved.Firstly, a basic backstepping controller is designed to generate virtual control input according to the attitude command.Secondly, contrary to the uncertainty of the control effectiveness matrix, a distribution error feedback loop is constructed in the INCA method to enhance the robustness.Thirdly, the control effectiveness matrix is reconstructed by different kinds of fault information to achieve an accurate command distribution.Finally, the stability of IRCCA is guaranteed by the Jury stability criterion.By using the IRCCA, the uncertainty and the actuator fault are all compensated in the CA method, which reduces the difficulty of the controller design.This paper is organized as follows: Section 2 introduces a basic backstepping attitude controller design and the faulttolerant control allocation problem.Section 3 introduces the IRCCA.Section 4 presents the stability analysis of the IRCCA.A high-fidelity aerodynamic model of the Innovative Control Effectors (ICE) aircraft is used in Section 5 to demonstrate the performance of the IRCCA method compared with INCA and CCA in the simulation environment.Finally, a summary is given in Section 6.

Backstepping Attitude Controller Design
The supersonic tailless aircraft attitude model is described in the input affine form: where x 1 ∈ R 3×1 denotes the roll, pitch, and yaw attitude angle, x 2 ∈ R 3×1 are the roll, pitch, and yaw angular rates, and ν ∈ R l is the virtual control input, such as moment coefficients.
A basic backstepping controller is designed as the following steps [42,43]: Step 1: The x 2 variable is regarded as a control input of the first relation in Equation (1).By introducing the error signal ξ 1 as: Taking the time derivative of both sides of Equation ( 2) and combining it with the first relation in Equation (1): .
A Lyapunov function is chosen as: Obviously, V 1 > 0 for all x 1 = x 1d .By taking the time derivative of Equation ( 4) and combing with Equation (3): .
According to Equation ( 5), a virtual control law for x 2 can be chosen as: where ω 1 > 0.
Substituting Equation (6) into Equation ( 5), and it is easy to prove that the tracking error converges to zero.
Step 2: Introducing the error signal ξ 2 as: Taking the time derivative of both sides of Equation ( 7) and combining it with the second relation in Equation (1): .
A Lyapunov function is chosen as: Obviously, V 2 > 0 for all x 2 = x 2d .Taking the time derivative of Equation ( 9) and combing it with Equation (8) .
It is easy to prove that the tracking error converges to zero.

Fault-Tolerant Control Allocation Problem
The relationship between control input and the virtual control input of supersonic tailless aircraft is described in the input affine form [36,44]: where the state influencing the virtual input, x ∈ R n , the control input vector, δ ∈ R m , and the virtual controls, ν ∈ R l , are assumed to be a nonlinear function of the aircraft state and control input G(•).
Assumption 1.The G(x, δ) can be expressed as where B ∈ R l×m is the control effectiveness matrix.
Considering the uncertainty of the control effectiveness matrix: where ∆B is the uncertainty of control effectiveness matrix, B r ∈ R l×m is the actual control effectiveness matrix, and The actuators of the aircraft can be influenced by many types of faults in flight missions.There are three typical faults which are the loss of effectiveness, lock-in-place fault, and loose fault.All of them are considered here.Assumption 2. Two actuator faults do not happen simultaneously on the same actuator.
The fault-tolerant control allocation problem can be described as: gain the virtual control input ν by the control law, then contrary to the actuator fault, design the control allocation scheme to determine the control input δ to reach the desired virtual control input ν.

Incremental Reconfiguration Closed-Loop Control Allocation Scheme Design
Consider the nonlinear control allocation problem without the fault and uncertainty expressed in Equation (3).With the control system observed to work at the discrete time, the system is described in discrete form, and the kth time step is denoted by k.The Equation ( 12) can be locally linearized at every time step according to first-order Taylor expansion: With the time scale separation principle [45], a time scale separation between the state on which the control action has direct effect and the controlled states is assumed.This means that for fast actuators and a small-time increment, a change of control input causes the state that it affects directly to change much faster than the controlled one.As a consequence, the difference between the state value is negligible.For the supersonic tailless aircraft, it is stated that a change in control input has a change in moment as effect.The change in moment is directly affecting the angular accelerations.The angular rates only change by integrating the angular accelerations, hence by integrating the control surface deflection component.It means that the control surface deflection affects the angular rate directly, which makes the (x(k +1) − x(k)) component, the change in angular rates, negligible at high sampling rates, which turned out to be true in [46] at sampling frequencies of 100 Hz.Hence, Equation ( 16) can be simplified to where ∆δ(k) = δ(k + 1) − δ(k).To simplify the notation, take the following substitution: ∂δ Define the input of incremental control allocation: The incremental nonlinear control allocation (INCA) is formulated as: given the current state x, the current control input δ, and virtual control command ν c , determine an incremental change in the actual control input vector ∆δ.
where B d can be computed by the data gained in the database of the over-actuated system, and (∆δ, ∆δ) are the constraints of ∆δ calculated by the (δ, δ) and where ∆t is the simulation step size.The problem can be transformed into solving the following constrained quadratic optimization problems: where W v ∈ R n×n and W δ ∈ R l×l are nonsingular weight matrices.Many mature linear control allocation algorithms [25][26][27][28][29] can solve Equation (17).
The common way is the pseudo-inverse method: where Based on this method, many optimization-based methods are proposed to match the constraints.
According to the null-space concept, the solution of this problem can be expressed as ∆δ = B + d v c + ∆δ 0 , where B d • ∆δ 0 = 0. ∆δ 0 is the adjustment to the pseudo-inverse solution B + d v c according to the constraints.The control input δ(k + 1) is computed as: And the actual output v r (k) can be derived as: Considering the uncertainty of the control effectiveness matrix in Equation ( 14), the B is transformed to be B r .Hence, the actual output of the IRCCA is derived as Equation (19).There is a discrepancy between B r and B. Despite the mismatch situation, a state feedback loop was designed.The accuracy of the control command allocation is implemented by the feedback of allocation error v, created by the uncertainty.Due to the incremental input of CA, the output of the linear CA method is the incremental control input, and the allocation error v cannot be directly feedbacked by the actual output v r (k).In this paper, the difference between the actual output of two adjacent steps (v r (k) − v r (k − 1)) and the incremental virtual input ν c is used to calculate the allocation error v in Equation (20).Applying the allocation error feedback, the ν c (k) is adjusted by the actual output of the previous step v r (k − 1) and the allocation error v(k − 1) in Equation (21).
Considering the actuator fault, the actuators are reconstructed to isolate faults and eliminate the error caused by the actuator fault.Based on the fault information, the control effectiveness matrix is reconfigured.Additionally, the input and output of the CA algorithm also change simultaneously.
(1) jth actuator with lock-in-place fault Referring to Equation (15), under the lock-in-place fault, the control effectiveness matrix , where E 1 contains the lock-in-place fault information.The actual output of IRCCA is derived as Equation (22).The additional moment generated by the fault actuator due to the lock-in-place fault reduces the accuracy of control allocation.Therefore, the virtual input is demanded to subtract the additional part.Referring to Equation ( 15), the additional part is derived as b j δ f j .However, b j cannot be gained in the actual situation due to B r being unknown.We replace B r by B d , and b j can be gained by Then the v(k) and v c (k) are rewritten in Equations ( 23) and (24).
(2) jth actuator with loose fault Referring to Equation (15), under the loose fault, the control effectiveness matrix B d is reconfigured to be B d (k) , where E 2 contains the loose fault information.Additionally, with the loose fault, the fault actuator deflection is 0 deg.The control input δ(k − 1) is transformed into (I − E 2 )δ(k − 1).The actual output of IRCCA is derived as Equation (25).The actuator with a loose fault cannot generate additional moments.Therefore, the v c (k) and v(k) retain the form of Equations ( 20) and (21).
(3) jth actuator with a loss of effectiveness fault Referring to Equation (15), under the loss of effectiveness fault, the control effectiveness matrix B d is reconfigured to be , where E 3 contains the loss of effectiveness fault information.The actual output of IRCCA is derived as Equation ( 26).The actuator with the loss of effectiveness fault cannot generate additional moments.Therefore, the v c (k) and v(k) retain the form of Equations ( 20) and (21).
(4) multiple fault scenarios Considering the multiple fault scenarios, the control effectiveness matrix Additionally, with the loose fault, the fault actuator deflection is 0 deg.The control input δ(k − 1) is transformed into (I − E 2 )δ(k − 1).The actual output of IRCCA is derived as Equation ( 27).The actuator with the loss of effectiveness fault and loose fault cannot generate additional moments.Therefore, the v c (k) and v(k) retain the form of Equations ( 23) and (24).
the whole fault-tolerant attitude control incorporating reconfiguration control allocation scheme is shown in Figure 1.
effectiveness fault and loose fault cannot generate additional moments.Therefore, the c ( ) retain the form of Equations ( 23) and (24).(27) the whole fault-tolerant attitude control incorporating reconfiguration control allocation scheme is shown in Figure 1.
Linear CA Actuator Model

Stability in the Absence of Actuator Faults
In this section, we consider the stability of IRCCA without the control effectiveness matrix uncertainty and actuator fault.

and the characteristic roots as being
( ) . For the IRCCA without the control effectiveness matrix uncertainty and actuator fault, the control allocation is stable if Proof of Theorem 1.To guarantee that the over-actuated attitude system is controllable, the full rank of B and d B are required [15].Hence, suppose can be rewritten by Equation (18) as: Based on Equations ( 20), ( 21) and ( 28), c ( ) k v can be also derived as: [ ( 1 Combining Equations ( 28) and ( 29), we have: In the actual system, the control effectiveness matrix is changing slowly, so the variation of B and d B can be neglected during a time step.Despite the small step size, the variation of δ and 0 Δ δ can be also neglected.In this section, we consider the stability of IRCCA without the control effectiveness matrix uncertainty and actuator fault.
, and the characteristic roots as being For the IRCCA without the control effectiveness matrix uncertainty and actuator fault, the control allocation is stable if λ i > 1.5 is satisfied.
Proof of Theorem 1.To guarantee that the over-actuated attitude system is controllable, the full rank of B and B d are required [15].Hence, suppose B(k)B + d (k) is a non-singular array; v c (k) can be rewritten by Equation (18) as: Based on Equations ( 20), ( 21) and (28), v c (k) can be also derived as: Combining Equations ( 28) and ( 29), we have: In the actual system, the control effectiveness matrix is changing slowly, so the variation of B and B d can be neglected during a time step.Despite the small step size, the variation of δ and ∆δ 0 can be also neglected.
Equation ( 30) can be simplified as: The input-output relation is obtained according to the z transform of Equation ( 32): According to the Jury stability criterion [47], the stability of the control allocation system is determined by the location of the closed-loop poles or the roots of the characteristic equation in the z plane.The characteristic equation is: The system is stable if any of the closed-loop characteristic roots lie inside the unit circle.According to Equation (31), the variation of Q(k) can be ignored, and Equation (33) can be transformed as: There is an invertible matrix P which can transform Q into an upper triangular matrix.
This completes the Proof of Theorem 1.

Stability in the Presence of Control Effectiveness Matrix Uncertainty
In this part, we consider the stability of IRCCA with the control effectiveness matrix uncertainty.
, and the characteristic roots of Q 1 are For the IRCCA with control effectiveness matrix uncertainty, the control allocation is stable if λ 1 i > 1.5 is satisfied.

Proof of Corollary 1.
Referring to Proof of Theorem 1, the v(k) is formulated as Equation ( 36) based on Equations ( 19)-( 21) with the control effectiveness matrix uncertainty: The input-output relation is obtained according to the z transform of Equation ( 36), and the stability problem can be transformed into guaranteeing all roots of lie inside the unit circle.Then, reach the solution: This completes the Proof of Corollary 1.

Stability in the Presence of Actuator Fault
In this part, we consider the stability for IRCCA with the control effectiveness matrix uncertainty and actuator faults.
For the IRCCA with the control effectiveness matrix uncertainty and actuator faults, the control allocation is stable if Proof of Corollary 2. Based on Equation ( 27), the v c (k) is rewritten in (37): Based on Equations ( 23), ( 24) and (37), and according to (31): Combing Equations ( 37) and (38), and according to (31): The input-output relation is obtained according to the z transform of Equation (39), and the stability problem can be transformed into guaranteeing all roots of lie inside the unit circle.Then, reach the solution: This completes the Proof of Corollary 2.

Stability Analysis for Fault-Tolerant Attitude Control System
Referring to the [48], the stability of the fault-tolerant attitude control system is analyzed subject to quantization.
The control system subject to quantization can be expressed as: .
where q τ (v) is the quantified virtual control input.
Defining the quantization error ξ q = q τ (v) − v, and the control system can be expressed as: .
And the dynamic expression of tracking error is: .
V ξ ≤ 0, where ρ( ξ q ) = 1 ω 3 H • ξ q , ω 3 equals to the smaller one of ω 1 and ω 2 , which means that the control system is input-to-state stable if ξ q is treated as the input of the system (40).
This completes the Proof of Theorem 2.

Simulation Results and Discussion
In this section, we consider the ICE supersonic tailless aircraft model (see [1] for more details) as the simulation object.The simulation model is described in Equations ( 42) and (43) referring to [2].
where µ, α, β are roll angle, angle of attack and sideslip angle; p, q, r are the roll, pitch, and yaw angular rates; χ, γ are the heading angle and flight-path angle; I xyz is the inertia matrix and expressed as Equation ( 44); l, m, n is the aerodynamic moment around the body axes; T BV is the transformation matrix from the body coordinate system to the velocity coordinate system and expressed as Equation (45); ρ is the air density, S is the aircraft wing area, V is the aircraft velocity, C l , C m , C n are the dimensionless rolling, pitching, yawing moment coefficients; b is the wing span; and c is the mean aerodynamic chord.The wing planform and inertia parameters are shown in Table 1.The operation ranges of the angular rates and attitude angles are illustrated in Table 2.These parameters are all derived from [1].Different from the [2], the specific aerodynamic moment coefficients of the ICE aircraft are described as a nonlinear form in Equation ( 46) [1].All actuators except multi-axis thrust vectoring (MTV) are used as the main research object in our study.The control inputs δ consist of the deflection of inboard leading-edge flap (ILEF) δ lilef , δ rilef , the outboard leadingedge flap (OLEF) δ lolef , δ rolef , the all-moving wing tips (AMT) δ lamt , δ ramt , the elevons δ le , δ re , the spoiler-slot deflectors (SSD) δ lssd , δ rssd , and pitch flaps (PF), δ pf .The position and rate limit of the control inputs are shown in Table 3.
where C la (x), C ma (x), C na (x) are the remaining part of the roll, pitch, yaw moment coefficients only influenced by the state x, and C lδ (x, δ), C mδ (x, δ), C nδ (x, δ) are the roll, pitch, yaw moment coefficients generated by the actuators' deflection under the state x.A highfidelity aerodynamic database of the ICE aircraft mainly obtained from wind tunnel tests was released by Lockheed Martin for academic use.This aerodynamic model is described in detail in [1].Additionally, the elastic influence is considered to be the uncertainty of the control effectiveness matrix, which is described in [1].As article [1] points out, the control effectiveness of the actuator is strongly nonlinear and coupled.Additionally, during the supersonic dash process, the influence of elastic cannot be ignored.
According to Equations ( 42) and ( 43), the backstepping control law can be set as: The attitude controller parameters are set as follows, referring to [47]: In the following stabilization control simulations, the proposed method is compared with the INCA in [41] and the RCA in [36].The RCA is also used with the fault information.Compared with the INCA, the effectiveness of reconfiguration design is demonstrated.Compared with the RCA, the consideration of the nonlinear characteristics is demonstrated.For INCA and IRCCA, the constant matrices in Equation ( 7) are set as: (200000,200000,200000) and W δ = diag(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), referring to [49].The constant matrices for RCA in [40] are set as: (200000,200000,200000).The initial attitude angles are chosen as µ = 6deg, a = 9deg, β = 5deg, and angular rates are set as p, q, r = 0.In addition, during the simulation process, the set χ, γ, .χ, .γ = 0, and the altitude and velocity are 500 ft and 1240 ft/s.The active set algorithm-based quadratic programming [26] and the backstepping controller are uniformly used to ensure consistency.The simulations are run in real-time with an ode4 solver with a fixed-step size of 0.01.The simulations were performed in a 64-bit computer with AMD Ryzen 7 5800X 8-Core Processor @ 4.20 GHz and 16 GB RAM.
The actual control effectiveness matrix is variable with the state and the deflection of control input.However, for the adaption of the RCA method, the actual control effectiveness matrix of the ICE model is transformed into the linearization form based on the aerodynamic database of the ICE aircraft and given by: 2.1 0.07 −0.32 −0.To verify the effectiveness and superiority of the proposed IRCCA method, the simulation process is divided into two parts: no actuator fault case and actuator fault case.In the first part, the performances of the three methods are tested without the fault case to verify their plausibility of the B r and the performance of the three methods.In the second part, with the actuator fault being set as Table 4, the effectiveness and superiority of the proposed method are verified by comparisons with the existing algorithm.

Fault Matrix
Parameters of Actuator Faults With the fault information, the B r used in RCA is reconfigured as: Firstly, the performances of methods are tested in the normal case.The attitude angle, the attitude angle rate, and the moment coefficient gained by the three methods are shown in Figures 2-4. Figure 2 shows that there is a bigger overshot and a longer settling time for the RCA method.It is because the RCA uses the fixed effectiveness matrix to allocate the virtual command, which ignores the nonlinear relationship between the actuator deflection and their generated moment.Moreover, in contrast to INCA method, there is a smaller overshot for the IRCCA method.It is because with the distribution error feedback loop, the uncertainty of the control effectiveness matrix is compensated, and the distribution of the virtual command can be more accurate than INCA. Figure 3 shows that all attitude angle rates met the limitation and the changes in amplitude in IRCCA and INCA are smaller than that in RCA. Figure 4 shows that the change in the moment coefficient in RCA always lags behind INCA and IRCCA.With consideration of nonlinear characteristics, INCA and IRCCA can generate an accurate moment coefficient according to the virtual command.Conversely, RCA can only compensate for the error by the outer control loop.In view of these results, as it can be observed all three methods can make the system stable and the IRCCA achieves the best performance to enforce the attitude system to be stable.Table 5 shows the actual simulation time with the three methods.It can be seen that the three methods meet the real-time requirement.Moreover, the simulation time with RCA is less than the other two methods.It is because RCA saves computation time by applying the fixed control effectiveness matrix, whereas the INCA and IRCCA need to compute and gain the control effectiveness matrix in real-time.

Lock-in-Place Fault
The performances of methods are tested with lock-in-place fault.The left elevon is stuck at 15 deg.The attitude angle, the attitude angle rate, and the moment coefficient gained by the three methods are shown in Figures 5-7.As can be observed the INCA cannot maintain the system stability without the fault information.In contrast, IRCCA and RCA can still maintain the stability of the system.Moreover, there is a smaller overshot for IRCCA than RCA. Figure 8 shows the expected control output δ le gained by the three methods.It can be seen that IRCCA and INCA can identify the accuracy of the fault information, whereas the δ le gained by the RCA cannot be stable.In view of these results, one can assume that the proposed IRCCA scheme makes a better performance to deal with the lock-in-place fault than the other two methods.

Lock-in-Place Fault
The performances of methods are tested with lock-in-place fault.The left ele stuck at 15 deg.The attitude angle, the attitude angle rate, and the moment coe gained by the three methods are shown in Figures 5-7.As can be observed the cannot maintain the system stability without the fault information.In contrast, I and RCA can still maintain the stability of the system.Moreover, there is a smalle shot for IRCCA than RCA. Figure 8 shows the expected control output le δ gained three methods.It can be seen that IRCCA and INCA can identify the accuracy of th information, whereas the le δ gained by the RCA cannot be stable.In view of th sults, one can assume that the proposed IRCCA scheme makes a better performa deal with the lock-in-place fault than the other two methods.

Loose Fault
The performances of the methods are tested with the loose fault.It is assume the right elevon has a loose fault.The attitude angle, the attitude angle rate, and th ment coefficient gained by the three methods are shown in Figures 9-11.As can served the output of INCA presents a chattering phenomenon.In contrast, IRCC RCA can still maintain the stability of the system.Moreover, there is a smaller ov for IRCCA than RCA. Figure 12 shows the expected control output re δ gained three methods.It can be seen that IRCCA and INCA can identify the accuracy of th information, whereas the re δ gained by the RCA cannot be stable.In view of th sults, one can obtain that the proposed IRCCA scheme makes a better performance with the loose fault than the other two methods.

Loose Fault
The performances of the methods are tested with the loose fault.It is assumed that the right elevon has a loose fault.The attitude angle, the attitude angle rate, and the moment coefficient gained by the three methods are shown in Figures 9-11.As can be observed the output of INCA presents a chattering phenomenon.In contrast, IRCCA and RCA can still maintain the stability of the system.Moreover, there is a smaller overshot for IRCCA than RCA. Figure 12 shows the expected control output δ re gained by the three methods.It can be seen that IRCCA and INCA can identify the accuracy of the fault information, whereas the δ re gained by the RCA cannot be stable.In view of these results, one can obtain that the proposed IRCCA scheme makes a better performance to deal with the loose fault than the other two methods.

Loss of Effectiveness Fault
The performances of the methods are tested with the loss of effectiveness fau assumed that the right SSD has a 50% loss of effectiveness fault.The attitude ang attitude angle rate, and the moment coefficient gained by the three methods are sho

Loss of Effectiveness Fault
The performances of the methods are tested with the loss of effectiveness fault.It is assumed that the right SSD has a 50% loss of effectiveness fault.The attitude angle, the attitude angle rate, and the moment coefficient gained by the three methods are shown in Figures 13-15.As can be observed the three methods can maintain the stability of the system.The tracking performance of INCA is almost the same as IRCCA.It is because the loss of effectiveness fault just changes the control effectiveness matrix, and does not instantaneously change the actuator deflections.Due to the high robustness against the change of control effectiveness matrix, the INCA can also attain a superior performance to deal with the loss of effectiveness fault.Figure 16 shows the expected control output δ rssd gained by the three methods.It can be seen that the expected control output δ rssd gained by RCA is different from INCA and IRCCA.This different distribution result is caused by the different CA methods.In view of these results, one can assume that the INCA and the proposed IRCCA make better performance to deal with the loss of effectiveness fault than the RCA methods.

Multiple Faults
Finally, the fault case is considered under a complex and harsh situation: the von is stuck at 15 deg, the right elevon has loose fault, and the right SSD has a 50% effectiveness fault.
The attitude angle, the attitude angle velocity, and the moment coefficient ga the four methods are shown in Figures 17-19.As can be observed the IRCCA achie best performance to enforce the attitude system to be stable, whereas the INCA maintain stability under the severe situation.The effectiveness of the reconfigurat sign can be seen by the comparison.Moreover, there is a smaller overshot in pitc mand for IRCCA contrasting with RCA, which also demonstrates the superiority proposed method.It can be seen in Figure 20

Multiple Faults
Finally, the fault case is considered under a complex and harsh situation: the left elevon is stuck at 15 deg, the right elevon has loose fault, and the right SSD has a 50% loss of effectiveness fault.
The attitude angle, the attitude angle velocity, and the moment coefficient gained by the four methods are shown in Figures 17-19.As can be observed the IRCCA achieves the best performance to enforce the attitude system to be stable, whereas the INCA cannot maintain stability under the severe situation.The effectiveness of the reconfiguration design can be seen by the comparison.Moreover, there is a smaller overshot in pitch command for IRCCA contrasting with RCA, which also demonstrates the superiority of the proposed method.It can be seen in Figure 20 that the stability condition of the control system ξ ≥ ρ( ξ q ) is satisfied throughout the whole process.As what is expected, the IRCCA successfully compensates the actuator fault with the fault information.

Conclusions
In this paper, an incremental reconfiguration closed-loop control allocation is proposed for supersonic tailless aircraft attitude control with actuator faults.With a closedloop distribution error feedback design in the INCA, the nonlinear control allocation problem is transformed into an incremental form and the uncertainty is compensated.Utilizing the fault information, the control allocation matrix is reconfigured to ensure that the

Conclusions
In this paper, an incremental reconfiguration closed-loop control allocation is proposed for supersonic tailless aircraft attitude control with actuator faults.With a closed-loop distribution error feedback design in the INCA, the nonlinear control allocation problem is transformed into an incremental form and the uncertainty is compensated.Utilizing the fault information, the control allocation matrix is reconfigured to ensure that the control system can deal with the actuator fault.With a basic backstepping attitude controller, the attitude control of supersonic tailless aircraft is achieved.Numerical simulation results illustrate that the proposed scheme can guarantee that the closed-loop system is stable, achieves accurate command tracking in the presence of actuator faults, and has better performance than the existing methods.Hence, the IRCCA achieves a more convenient design and better performance than the traditional FTC.However, the control input in this paper does not consider the multi-axis thrust vectoring.This paper mainly focuses on the nonlinearity and uncertainty of the control effectiveness matrix and the actuator fault.Severe external disturbance is not considered.Therefore, the robust control containing the multi-axis thrust vectoring towards severe external disturbance will be the focus of our future research.

Figure 2 .
Figure 2. Time response of the attitude angle with no actuator fault.

Figure 3 .
Figure 3.Time response of the attitude angle rate with no actuator fault.

Figure 2 .Figure 2 .
Figure 2. Time response of the attitude angle with no actuator fault.

Figure 3 .
Figure 3.Time response of the attitude angle rate with no actuator fault.Figure 3. Time response of the attitude angle rate with no actuator fault.

Figure 3 .
Figure 3.Time response of the attitude angle rate with no actuator fault.Figure 3. Time response of the attitude angle rate with no actuator fault.

Figure 4 .
Figure 4. Time response of the moment coefficient with no actuator fault.

Figure 4 . 1 Figure 5 .
Figure 4. Time response of the moment coefficient with no actuator fault.Aerospace 2023, 10, x FOR PEER REVIEW 1

Figure 5 .
Figure 5.Time response of the attitude angle with lock-in-place fault.Figure 5. Time response of the attitude angle with lock-in-place fault.

Figure 5 .
Figure 5.Time response of the attitude angle with lock-in-place fault.

Figure 6 .
Figure 6.Time response of the attitude angle rate with lock-in-place fault.Figure 6.Time response of the attitude angle rate with lock-in-place fault.

Figure 6 .Figure 7 .
Figure 6.Time response of the attitude angle rate with lock-in-place fault.Figure 6.Time response of the attitude angle rate with lock-in-place fault.Aerospace 2023, 10, x FOR PEER REVIEW

Figure 7 .
Figure 7. Time response of the moment coefficient with lock-in-place fault.Figure 7. Time response of the moment coefficient with lock-in-place fault.

Figure 7 .
Figure 7. Time response of the moment coefficient with lock-in-place fault.

Figure 8 .
Figure 8.Time response of the fault actuator le δ deflection.

Figure 8 .
Figure 8.Time response of the fault actuator δ le deflection.

Aerospace 2023 ,
10, x FOR PEER REVIEW

Figure 9 .
Figure 9.Time response of the attitude angle with loose fault.Figure 9. Time response of the attitude angle with loose fault.

Figure 9 .
Figure 9.Time response of the attitude angle with loose fault.Figure 9. Time response of the attitude angle with loose fault.

Figure 9 .
Figure 9.Time response of the attitude angle with loose fault.

Figure 10 .
Figure 10.Time response of the attitude angle rate with loose fault.Figure 10.Time response of the attitude angle rate with loose fault.

Figure 10 . 2 Figure 11 .
Figure 10.Time response of the attitude angle rate with loose fault.Figure 10.Time response of the attitude angle rate with loose fault.Aerospace 2023, 10, x FOR PEER REVIEW 2

Figure 11 .
Figure 11.Time response of the moment coefficient with loose fault.Figure 11.Time response of the moment coefficient with loose fault.

Figure 11 .
Figure 11.Time response of the moment coefficient with loose fault.

Figure 12 .
Figure 12.Time response of the fault actuator re δ deflection.
Figures 13-15.As can be observed the three methods can maintain the stability of th tem.The tracking performance of INCA is almost the same as IRCCA.It is becau loss of effectiveness fault just changes the control effectiveness matrix, and does stantaneously change the actuator deflections.Due to the high robustness again change of control effectiveness matrix, the INCA can also attain a superior perfor to deal with the loss of effectiveness fault.Figure 16 shows the expected control rssd δ gained by the three methods.It can be seen that the expected control outpu gained by RCA is different from INCA and IRCCA.This different distribution re caused by the different CA methods.In view of these results, one can assume th INCA and the proposed IRCCA make better performance to deal with the loss of tiveness fault than the RCA methods.

Figure 12 .
Figure 12.Time response of the fault actuator δ re deflection.

Aerospace 2023 ,
Figure 13.Time response of the attitude angle with a loss of effectiveness fault.Figure 13.Time response of the attitude angle with a loss of effectiveness fault.

Figure 13 .
Figure 13.Time response of the attitude angle with a loss of effectiveness fault.Figure 13.Time response of the attitude angle with a loss of effectiveness fault.

Figure 13 .
Figure 13.Time response of the attitude angle with a loss of effectiveness fault.

Figure 14 .
Figure 14.Time response of the attitude angle rate with a loss of effectiveness fault.Figure 14.Time response of the attitude angle rate with a loss of effectiveness fault.

Figure 14 .Figure 15 .
Figure 14.Time response of the attitude angle rate with a loss of effectiveness fault.Figure 14.Time response of the attitude angle rate with a loss of effectiveness fault.Aerospace 2023, 10, x FOR PEER REVIEW

Figure 15 .
Figure 15.Time response of the moment coefficient with a loss of effectiveness fault.Figure 15.Time response of the moment coefficient with a loss of effectiveness fault.

Figure 15 .
Figure 15.Time response of the moment coefficient with a loss of effectiveness fault.

Figure 16 .
Figure 16.Time response of the fault actuator rssd δ that the stability condition of the system || || (|| ||) the whole process.As what is ex the IRCCA successfully compensates the actuator fault with the fault information.

Figure 16 .
Figure 16.Time response of the fault actuator δ rssd deflection.

Aerospace 2023 ,Figure 17 .
Figure 17.Time response of the attitude angle with multiple faults.Figure 17.Time response of the attitude angle with multiple faults.

Figure 17 .
Figure 17.Time response of the attitude angle with multiple faults.Figure 17.Time response of the attitude angle with multiple faults.

Figure 17 .
Figure 17.Time response of the attitude angle with multiple faults.

Figure 18 .
Figure 18.Time response of the attitude angle rate with multiple faults.Figure 18.Time response of the attitude angle rate with multiple faults.

Figure 18 .Figure 19 .
Figure 18.Time response of the attitude angle rate with multiple faults.Figure 18.Time response of the attitude angle rate with multiple faults.Aerospace 2023, 10, x FOR PEER REVIEW

Figure 19 .
Figure 19.Time response of the moment coefficient with multiple faults.Figure 19.Time response of the moment coefficient with multiple

Figure 19 .
Figure 19.Time response of the moment coefficient with multiple faults.

Table 1 .
Wing planform and inertia parameters for the model.

Table 2 .
The operation ranges of the angular rates and attitude angles.

Table 3 .
Dynamic characteristics of the ICE control actuators.

Table 4 .
The parameters of actuator faults.

Table 5 .
The simulation time with no actuator fault.

Table 5 .
The simulation time with no actuator fault.