Adaptive Feed-Forward Control for Gust Load Alleviation on a Flying-Wing Model Using Multiple Control Surfaces

: Based on measured gust information, a multi-input multi-output (MIMO) adaptive feed-forward control scheme for gust load alleviation (GLA) on a semi-span ﬂying-wing aircraft using multiple control surfaces is proposed. In order to remedy weight drift and biased estimation problems that are commonly encountered in adaptive control, the circular leaky LMS (CLLMS) algorithm is employed, which utilizes gust measurement information, ﬁltered reference signals, and error signals to update controller parameters online. The results demonstrate that good load reductions are achieved in both continuous and discrete gust environments. For instance, the designed GLA control system leads to an 80.72% reduction in the root-mean-square (RMS) values of wing-root bending moment in the Dryden gust environment and a 77.59% reduction of its maximum value in the 1-cos discrete gust condition. Based on the limited power of the actuator and the limited authority for control surface deﬂections when integrating GLA into the ﬂight control system, a weight-updating algorithm with deﬂection angle and rate constraints on control surfaces is proposed. The simulation results show that the strict constraints on control surface deﬂections will degrade the GLA performance. Finally, the inﬂuence of the partial jamming fault of actuators on GLA performance is studied. It is found that good GLA performance can be preserved despite the degraded performance during the initial stage of the actuator jamming fault. This is due to the robustness brought about by multiple control surfaces and the adaptability of the control algorithm.


Introduction
In aviation parlance, a gust refers to any sudden change in the magnitude or direction of air flow motion [1].From the perspective of aerodynamics, gust can produce a perturbation in the free-stream velocity experienced by the lifting surfaces on a vehicle, which will result in an additional perturbation in the total angle of attack of the lifting surfaces and dynamic pressure.The change in lift caused by variation in the angle of attack can excite structural vibrations, producing incremental internal loads in the aircraft structure, thereby reducing ride comfort.For example, compared to the conventional wing-tube configuration, the blended wing body (BWB) aircraft provides a substantial enhancement in fuel efficiency, but it exhibits pronounced susceptibility to gusts.In this case, gust load alleviation (GLA) technology must be used to minimize the impact of such instantaneous loads as much as possible during the aircraft life cycle [2,3].
The active load alleviation technique has been established as an effective approach for mitigating gust loads in a wider frequency band, in which a set of control surfaces (leading edge or trailing edge control surfaces) with control authority can be used to regulate local lift [4].This technique has been successfully implemented since the 1970s with the LockheedL-1011 and recently in the Airbus A380 and the Boeing 787 [5].Utilizing the active GLA technique, the B-2 Stealth Bomber achieves up to a 50% reduction in gust loads by combining the elevators and the other control surfaces [6].
The feedback control, as a straightforward method, has been extensively utilized for GLA, in which the controllers use the feedback signals to form the deflection commands of control surfaces to produce the incremental aerodynamic forces to suppress the gustinduced responses.Typical control design uses the empirical tuning-based PID method [7] and the state-space and linear model-based LQG method [8].As a new feedback-based control strategy, Yue and Zhao [9] first employed the equivalent input disturbance (EID) method to alleviate the gust loads of a folding wing model.The advantage of feedback control lies in the robustness towards the uncertainties in the controlled object.However, its limitation lies in the delayed control actions since the disturbance effects must propagate through the entire system before exhibiting their effects, which can be compensated by the feed-forward control strategy based on a device to measure the disturbance in advance.Alam et al. [10] investigated a mixed feed-forward/feedback approach for a flexible aircraft, which indicates a significant improvement in gust alleviation.
The hardware devices capable of measuring turbulence velocity ahead of an aircraft can provide a necessary prerequisite for a feed-forward control design of GLA.Currently, there are two types of onboard sensors available for measuring gust disturbances: LIDAR (Light Detection And Ranging) and alpha sensors [11,12].LIDAR is capable of detecting gust velocity fields approximately 50~150 m ahead of an aircraft, providing sufficient lead time for the control system to make appropriate adjustments.On the other hand, alpha sensors are used to measure the angle of attack at the aircraft's nose, which enables the determination of the gust angle of attack after combining with the current flight parameters.
Feed-forward control can provide an improved vibration reduction when suitable reference signals are available [13,14].With the gust disturbance information available, feedforward control frequently exhibits a greater ability to minimize the impact of measured disturbances on the output compared to feedback control alone [15].The drawback of a feedback-based controller is caused by the time delay between the sensed information; a feed-forward control strategy may overcome this disadvantage when a priori knowledge of the disturbance is available [16,17].Wildschek [18][19][20] developed a feed-forward vibration control system, which has been successfully tested in flight.Usually, the adaptive feedforward controller is implemented by a finite impulse response (FIR) filter, in which the filter coefficients are updated online using the least mean squares (LMS) algorithm [21,22].However, due to errors in quantization and truncation, the LMS algorithm tends to suffer from issues of weight coefficient drift and overflow during long-term operation.The leaky LMS algorithm can remedy this problem but it will bring a biased estimation of filter coefficients [23].Aiming at the weight drift and bias problems, Zhao et al. [24] used the circular-leaky LMS algorithm (CLLMS) to update the weights of the adaptive feed-forward controller.The study showed that gust loads were significantly reduced in both stationary and nonstationary gusty environments.
To date, many studies have dealt with GLA control utilizing multiple control surfaces [25][26][27].Fournier et al. [2] studied the multiple-input multiple-output robustness of the controller designed by H ∞ and µ synthesis based on the available measurements.In fact, multiple control surfaces can significantly enhance the GLA performance.When any failure of a control surface appears in flight, the corresponding closed-loop system may still have good GLA performance, because multiple control surfaces provide a certain degree of robustness for a faulty system.In view of the above, this paper focuses on the GLA control of a flying-wing BWB aircraft equipped with multiple control surfaces along the wing span.To adapt to the layouts of multiple control surfaces and multi-sensors, a multi-input multi-output (MIMO) adaptive feed-forward control scheme is proposed to alleviate the gust loads.Assuming that an alpha sensor located at the nose of the aircraft can be used to measure the gust-induced angle-of-attack signal, which is used to generate a set of reference signals for feed-forward controllers, the error signals used for updating the adaptive controllers are the outputs of the acceleration sensors on the wing.To ensure that the GLA system operates in a stable environment, the CLLMS algorithm is employed to update the weight coefficients in each FIR filter online.When integrating GLA into the flight control system, it must be remembered that no active flight control system is generally allowed to make use of the full-range deflection of control surfaces.That is to say, in practice, GLA has limited control authority.Therefore, constraints on control surface deflections should be considered.To this end, the output constraint algorithm is developed to account for the allowed limits of the control surface deflections.Finally, in order to verify the robustness of the designed GLA control system to the specified failure mode of the control surface, the GLA performance under a control surface jamming failure is examined.All simulation results demonstrated that the proposed MIMO adaptive feed-forward GLA using multiple control surfaces is effective for a flying-wing BWB aircraft.

Aeroservoelastic Modeling of the Flying Wing 2.1. Structural and Aerodynamic Models
Similar to the UAV demonstrator designed by the University of Minnesota [28], a flying-wing layout aircraft is used for GLA control.The model has a root chord of 878 mm, a wingtip chord of 240 mm, and a half-span length of 1578 mm.For the convenience of 3D printing, the entire wing is made of ULTEM 9085 filament.Only the semi-span wind tunnel model is considered in this paper, and the corresponding NASTRAN finite element model (FEM) is shown in Figure 1a.The first six elastic mode shapes and natural frequencies of the wing clamped at the wing root are shown in Figure 2 and Table 1, respectively.the GLA system operates in a stable environment, the CLLMS algorithm is employed to update the weight coefficients in each FIR filter online.When integrating GLA into the flight control system, it must be remembered that no active flight control system is generally allowed to make use of the full-range deflection of control surfaces.That is to say, in practice, GLA has limited control authority.Therefore, constraints on control surface deflections should be considered.To this end, the output constraint algorithm is developed to account for the allowed limits of the control surface deflections.Finally, in order to verify the robustness of the designed GLA control system to the specified failure mode of the control surface, the GLA performance under a control surface jamming failure is examined.All simulation results demonstrated that the proposed MIMO adaptive feed-forward GLA using multiple control surfaces is effective for a flying-wing BWB aircraft.

Structural and Aerodynamic Models
Similar to the UAV demonstrator designed by the University of Minnesota [28], a flying-wing layout aircraft is used for GLA control.The model has a root chord of 878 mm, a wingtip chord of 240 mm, and a half-span length of 1578 mm.For the convenience of 3D printing, the entire wing is made of ULTEM 9085 filament.Only the semi-span wind tunnel model is considered in this paper, and the corresponding NASTRAN finite element model (FEM) is shown in Figure 1a.The first six elastic mode shapes and natural frequencies of the wing clamped at the wing root are shown in Figure 2 and Table 1, respectively.In this study, the unsteady aerodynamic forces are computed using the Doub tice Method (DLM), which is a subsonic model for harmonically oscillating lifting s [29].Figure 1b shows the aerodynamic model of the flying-wing aircraft, in wh lifting surface is represented by a number of aerodynamic boxes of trapezoidal sha

State-Space Open-Loop Aeroservoelastic Equations
In the DLM aerodynamic model, the non-penetrating boundary condition is u establish the relationship between the dimensionless downwash velocity and the ha ically oscillating lifting surface at the control points on the aerodynamic boxes.Th pressure coefficients on aerodynamic boxes can be calculated via the aerodynami ence coefficient (AIC) matrix and dimensionless downwash velocity.Finally, a spli trix is constructed to establish the data transfer relationship between the structu placements at the finite element nodes and the displacements at the control points aerodynamic boxes.Using the equivalent principle of virtual work and introduc In this study, the unsteady aerodynamic forces are computed using the Doublet Lattice Method (DLM), which is a subsonic model for harmonically oscillating lifting surfaces [29].Figure 1b shows the aerodynamic model of the flying-wing aircraft, in which the lifting surface is represented by a number of aerodynamic boxes of trapezoidal shape.

State-Space Open-Loop Aeroservoelastic Equations
In the DLM aerodynamic model, the non-penetrating boundary condition is used to establish the relationship between the dimensionless downwash velocity and the harmonically oscillating lifting surface at the control points on the aerodynamic boxes.Then, the pressure coefficients on aerodynamic boxes can be calculated via the aerodynamic influence coefficient (AIC) matrix and dimensionless downwash velocity.Finally, a spline matrix is constructed to establish the data transfer relationship between the structural displacements at the finite element nodes and the displacements at the control points on the aerodynamic boxes.Using the equivalent principle of virtual work and introducing the modal coordinate transformation, one obtains the following aeroelastic equation expressed in the modal space: in which M hh , D hh , and K hh are modal mass, damping, and stiffness matrices, respectively.M hc is the mass coupling matrix.Q hh , Q hc , and Q hg are matrices of generalized aerodynamics force (GAF) due to structural motion, deflections of control surfaces, and gust disturbance, respectively.
In order to cast the aeroelastic equation of motion in a time-domain constant coefficient equation, the minimum state (MS) [30] approximation is used.It assumes that the frequencydomain aerodynamic matrices can be approximated as follows: where where the state vector x ae includes modal displacement ξ, modal velocity .ξ, and the aerodynamic states.u ae is the actual control surface deflection vector.The gust disturbance vector w g contains vertical gust velocity w g and its time derivative .w g , namely The dynamic model of the actuator driving the ith control surface is specified by a third-order transfer function, having the form where u ci is the servo-commanded control surface deflection and δ i is the actual deflection.Note that the DC (s = 0) gain of the transfer function is one and is, thus, in steady state u ac,i = δ i .The assembled actuator state space equation is, thus, cast in the following form: In order to cast the aeroelastic equation of motion in a time-domain constant coefficient equation, the minimum state (MS) [30] approximation is used.It assumes that the frequency-domain aerodynamic matrices can be approximated as follows: where [ ] where the state vector ae x includes modal displacement ξ , modal velocity ξ  , and the aerodynamic states.ae u is the actual control surface deflection vector.The gust disturb- ance vector g w  contains vertical gust velocity g w and its time derivative g w  , namely The dynamic model of the actuator driving the th i control surface is specified by a third-order transfer function, having the form where ci u is the servo-commanded control surface deflection and i δ is the actual deflec- tion.Note that the DC ( 0 s = ) gain of the transfer function is one and is, thus, in steady state The assembled actuator state space equation is, thus, cast in the following form: or The augmentation of ae x in Equation (3) to include the actuator states ac x yields the following open-loop aeroelastic equation: where The output equation is The augmentation of x ae in Equation (3) to include the actuator states x ac yields the following open-loop aeroelastic equation: where The output equation is The output vector y p can be sensor readings, such as displacement, velocity, and acceleration responses at the specified finite element nodes, as well as the virtual load outputs, such as wing-root shear force, bending, and torsional moments.Equations (7) and (9) constitute the final open-loop aeroelastic equations of motion excited by control surface deflections and atmospheric gust.The purpose of GLA is to seek a control law that drives the control surfaces to deflect appropriately, thereby generating additional aerodynamic forces to minimize the outputs of the system.

Control Surface and Sensor Layouts
The design of a MIMO feed-forward control system requires an appropriate reference signal that can provide the necessary disturbance information for controllers to calculate the control commands.This study assumes that there is an alpha probe at the aircraft nose available to measure the vertical gust velocity signal [32], as shown in Figure 3. Afterwards, the effective angle induced by gust can be computed from alpha probe measurements, which is used as the reference signal for feed-forward control.
acceleration responses at the specified finite element nodes, as well as the virtual load outputs, such as wing-root shear force, bending, and torsional moments.Equations (7) and (9) constitute the final open-loop aeroelastic equations of motion excited by control surface deflections and atmospheric gust.The purpose of GLA is to seek a control law that drives the control surfaces to deflect appropriately, thereby generating additional aerodynamic forces to minimize the outputs of the system.

Control Surface and Sensor Layouts
The design of a MIMO feed-forward control system requires an appropriate reference signal that can provide the necessary disturbance information for controllers to calculate the control commands.This study assumes that there is an alpha probe at the aircraft nose available to measure the vertical gust velocity signal [32], as shown in Figure 3. Afterwards, the effective angle induced by gust can be computed from alpha probe measurements, which is used as the reference signal for feed-forward control.In this study, a total of four control surfaces are used, and they are depicted in Figure 3, whose cut-off frequency is set as fifteen hertz.Four acceleration sensors located at nodes 647, 5655, 11098, and 276 are used to monitor the responses of the system.The output signals also include the incremental internal loads relative to the static equilibrium of the structure, which are not measured in the GLA control system and can be regarded as the virtual outputs of the system.In addition, acceleration signals at nodes 647 and 11098 are also employed to form the error signals in the subsequent adaptive feed-forward control design.

MIMO Adaptive Feed-Forward Control Scheme
Figure 3 illustrates the framework of the designed MIMO adaptive feed-forward control for GLA.Due to the high proportion of control surfaces in the span, two acceleration sensors at nodes 647 and 11098 are employed as error sensors.From the open-loop aeroelastic equation (see Equation ( 7)), we can see that both gust and control surface deflections can affect the outputs of the system.Therefore, the inputs of the system affect the outputs through two paths.The vertical gust velocity g w causes the acceleration responses i d , In this study, a total of four control surfaces are used, and they are depicted in Figure 3, whose cut-off frequency is set as fifteen hertz.Four acceleration sensors located at nodes 647, 5655, 11098, and 276 are used to monitor the responses of the system.The output signals also include the incremental internal loads relative to the static equilibrium of the structure, which are not measured in the GLA control system and can be regarded as the virtual outputs of the system.In addition, acceleration signals at nodes 647 and 11098 are also employed to form the error signals in the subsequent adaptive feed-forward control design.

MIMO Adaptive Feed-Forward Control Scheme
Figure 3 illustrates the framework of the designed MIMO adaptive feed-forward control for GLA.Due to the high proportion of control surfaces in the span, two acceleration sensors at nodes 647 and 11098 are employed as error sensors.From the open-loop aeroelastic equation (see Equation ( 7)), we can see that both gust and control surface deflections can affect the outputs of the system.Therefore, the inputs of the system affect the outputs through two paths.The vertical gust velocity w g causes the acceleration responses d i , i = 1, 2 at nodes 647 and 11098 via the primary path (PP) transfer function G d1 (s) and G d2 (s), respectively.From gust velocity measurements, the effective angle of attack α w induced by gust can be calculated using w g /U.Afterwards, the obtained reference signal α w is fed forward to the adaptive controllers H ci (s) to generate the servo-commanded control surface deflections u ci , i = 1 ∼ 4. On the other hand, each control surface deflection u ci causes the outputs at nodes 647 and 11098 via the secondary path (SP) transfer function G c1i (s) and G c2i (s), respectively.Thus, the total signals y i , i = 1, 2 are obtained by summing these accelerations due to these control surface deflections.The error signals at nodes 647 and 11098, denoted as e 1 (t) and e 2 (t), are formed by summing d i and y i , i = 1, 2, respectively.In Figure 4, Ĝc1i (s) and Ĝc2i (s), as the estimated transfer functions of G c1i (s) and G c2i (s), respectively, are used to generate the filtered reference signals rc1i and rc2i .The filtered reference signals together with the error signals are used to update the parameters of adaptive controllers.The objective of GLA control is to adjust the controller parameters to minimize the error signals as much as possible.
tained by summing these accelerations due to these control surface deflections.The error signals at nodes 647 and 11098, denoted as 1 ( ) e t and 2 ( ) e t , are formed by summing i d , respectively.In Figure 4,

FIR Filter
The adaptive feed-forward controllers ( ) ci H s , 1 4 i =  , are implemented by a set of digital filters.In this study, a finite impulse response (FIR) filter is used to construct the adaptive control law due to its prominent advantage in stability.
The transfer function of an FIR filter can be written as , where ,0 , , , , are the weight coefficients of filter j , N is the order of the filter, and c N is the number of the controllers.

FIR Filter
The adaptive feed-forward controllers H ci (s), i = 1 ∼ 4, are implemented by a set of digital filters.In this study, a finite impulse response (FIR) filter is used to construct the adaptive control law due to its prominent advantage in stability.
The transfer function of an FIR filter can be written as where h j,0 , h j,1 , • • • , h j,N−1 are the weight coefficients of filter j, N is the order of the filter, and N c is the number of the controllers.The feed-forward control command is generated by filtering a sampled version of the sensed gust-induced angle-of-attack signal α w (t) via an FIR filter.Thus, the commanded deflection of a control surface u cj (n) at discrete time step n can be written as where the weight coefficient vector of the filter is and the gust-induced angle-of-attack vector α(n) is defined as

Adaptive Updating Law of the Weight Coefficients
As is well known, the LMS algorithm is widely used in the active control of vibration and noise due to its easy implementation and low computational cost.However, the LMS algorithm will encounter stability issues during long-term operation in a quantized digital computing environment.Specifically, the weight coefficients of the filter are prone to slowly drift, which will eventually lead to overflow.This drawback can be remedied by the leaky least mean squares (LLMS) algorithm, in which a leakage factor is introduced to suppress the weight-drift problem [33].
The cost function J(n) in the LLMS algorithm can be written as where γ ∈ (0, 1) is the introduced leakage factor γ that leaks the excessive energy associated with weight drift.N e is the number of the error signals (N e = 2 in this study).The error signals e i are defined as where The element rcij (n) is the output of the system Ĝcij (s) at discrete time step n.Substituting Equation (15) into Equation ( 14), the cost function is updated as follows Differentiating this with respect to the jth weight coefficient vector h cj gives The gradient descent algorithm yields the following weight updating law or where µ is the convergence coefficient.For simplicity, coefficient 2 is omitted in Equation (20).However, the obtained weights will produce a biased estimate of the Wiener solution in the LLMS algorithm.As an improved version of the LLMS algorithm, the CLLMS algorithm used in this study can ensure an unbiased estimate of the filter coefficients [34].In the CLLMS algorithm, the updating rules for the weight coefficients of the filter can be written as where with k = mod(n, N), satisfying h j,k (n) ≥ C 1 .The leakage factor γ c is calculated by where

Flutter Analysis
In order to ensure that the subsequent GLA simulations can be performed within the safe flight envelope, flutter analysis of the flying wing is carried out by tracing the eigenvalues of the system matrix A p in Equation ( 7) at different flow speeds.Figure 5a shows a root locus diagram of the system matrix A p , in which each arrow indicates the direction of an increase in flow speed.It can be seen that, as the flow speed increases, the second-order mode passes through the imaginary axis and becomes unstable first.The eigenvalue with zero real part on the mode 2 branch is a critical point, at which the flutter speed is found to be 75.41m/s with a flutter frequency of 10.1 Hz. Figure 5b illustrates that the flutter shape is primarily a superposition of the second vertical bending and the first torsional modes.The V-g and V-f diagrams obtained from frequency-domain flutter analysis are also shown in Figure 5c,d.

Simulation Environment for GLA Control
After flutter analysis is completed, GLA control simulations can be carried out.Figure 6 illustrates a simulation block diagram of the adaptive feed-forward GLA control system.The feed-forward controller utilizes three different types of signals, namely acceleration error signals e i (n), filtered reference signals rcij (n), and the gust-induced angle-of-attack vector α(n), to generate deflection commands of control surfaces.The adaptive controller is implemented using a level-2 S-function in the MATLAB R2018a Simulink platform.First, the time-series vector α(n) is constructed by the sensed gust-induced angle-of-attack signal α w (n).Afterwards, the weight coefficients of each FIR filter are updated.Finally, the commanded control surface deflections at discrete time step n are generated and used as the outputs of the level-2 S-function.
of an increase in flow speed.It can be seen that, as the flow speed increases, the secondorder mode passes through the imaginary axis and becomes unstable first.The eigenvalue with zero real part on the mode 2 branch is a critical point, at which the flutter speed is found to be 75.41m/s with a flutter frequency of 10.1 Hz. Figure 5b illustrates that the flutter shape is primarily a superposition of the second vertical bending and the first torsional modes.The V-g and V-f diagrams obtained from frequency-domain flutter analysis are also shown in Figure 5c,d.

Simulation Environment for GLA Control
After flutter analysis is completed, GLA control simulations can be carried out.  .The time history of gust velocity for the Dryden gust model is generated using a shaped filter excited by white noise.In Simulink simulations, the sampling period is taken as 0.001 s.The order of the FIR filter for each adaptive controller is set as N = 42.When the GLA control system is activated, the filter coefficients are updated online by the CLLMS algorithm, in which the convergence Since GLA control always operates below the flutter speed of the system, all numerical simulations are performed at a cruise speed of 27.2 m/s and sea level.The continuous gust scale in the Dryden model introduced in Section 4.3 is set to L w = 53.3m, and the RMS value of the gust velocity is taken as σ w = 0.5 m/s.The time history of gust velocity for the Dryden gust model is generated using a shaped filter excited by white noise.In Simulink simulations, the sampling period is taken as 0.001 s.The order of the FIR filter for each adaptive controller is set as N = 42.When the GLA control system is activated, the filter coefficients are updated online by the CLLMS algorithm, in which the convergence coefficient is µ = 0.005, leakage factor is γ = 0.5, and constants C 1 and C 2 are taken as 0.5 and 0.7, respectively.

GLA Control in a Continuous Dryden Gust Field
The widely used standard gust models can be divided into two categories: continuous gusts and discrete gusts.The former employs a statistical approach to describe the gust disturbance, including two standard forms: Dryden gust and von Karman gust.In this study, the Dryden model is used because it can be represented as a rational function form and it is easier to obtain a realization in the time domain.The power spectral density (PSD) of vertical gust velocity in the Dryden model can be written as [35] where τ w = L w /U, and L w are the gust scale.σ w represents the root-mean-square (RMS) value of the gust velocity.ω is the angular frequency.
In order to perform time-domain simulations, Dryden gust can be expressed in the following state-space form where .
x g (t) x g (t) .
w g (t) , where a is the cut-off frequency of the low-pass filter.r(t) is the Gaussian white noise input.Dryden gust distribution experienced by the wing is shown in Figure 7.
Aerospace 2023, 10, x FOR PEER REVIEW 12 of 24 where where a is the cut-off frequency of the low-pass filter.( ) r t is the Gaussian white noise input.Dryden gust distribution experienced by the wing is shown in Figure 7.The simulation parameters in the Dryden gust condition are shown in Section 4.2. Figure 8 shows the convergence behaviors of the selected weight coefficients of an FIR controller during the training process.The initial values of the FIR filter coefficients are all set to zero.Obviously, a converged controller is obtained at t = 870 s.The converged coefficients still exhibit oscillations around their average values due to the stochastic nature of the gust disturbance and instantaneous approximation of the gradient of cost function during the weight updating.The simulation parameters in the Dryden gust condition are shown in Section 4.2. Figure 8 shows the convergence behaviors of the selected weight coefficients of an FIR controller during the training process.The initial values of the FIR filter coefficients are all set to zero.Obviously, a converged controller is obtained at t = 870 s.The converged coefficients still exhibit oscillations around their average values due to the stochastic nature of the gust disturbance and instantaneous approximation of the gradient of cost function during the weight updating.In simulations, Dryden gust velocity and its first derivative with respect to time are computed using Equations ( 26) and (27).Subsequently, the trained CLLMS controller is employed to verify the effectiveness of GLA control in the Dryden gust condition.Figure 9 illustrates the acceleration responses at the monitoring nodes when the flying wing equipped with a GLA system travels through the Dryden gust field.Obviously, the GLA control can effectively suppress vibrations at the monitoring nodes on the wing.As shown in Figure 10, after activating the GLA control, the wing-root loads are significantly reduced.The actual deflections of control surfaces, namely the outputs of actuators, are shown in Figure 11.The statistical results in Table 2 indicate that the error signals, RMS values of the accelerations at monitoring nodes, and wing-root internal loads are significantly reduced.For example, an 80.72% reduction in the RMS value of the wing-root bending moment is achieved.This is of great significance for the weight reduction design of the aircraft.The GLA results demonstrate that the designed CLLMS-based controllers are very effective.The statistical results in Table 2 indicate that the error signals, RMS values of the accelerations at monitoring nodes, and wing-root internal loads are significantly reduced.For example, an 80.72% reduction in the RMS value of the wing-root bending moment is achieved.This is of great significance for the weight reduction design of the aircraft.The GLA results demonstrate that the designed CLLMS-based controllers are very effective.It is noted that a set of time-history curves cannot easily be used to simultaneously display the correlations of the GLA performances between two physical quantities.This issue can be remedied by plotting a convex hull diagram about two physical variables [36].First, the gust-induced response data for a set of physical quantities are computed by open-and closed-loop simulations, respectively.Then, the discrete points corresponding to two selected physical quantities are plotted on a two-dimensional plane.Finally, a convex hull is formed by the minimum convex polygon, which encompasses all these discrete points.The correlations of GLA performance between two selected physical quantities are clearly demonstrated in Figure 12.The statistical results in Table 2 indicate that the error signals, RMS values of celerations at monitoring nodes, and wing-root internal loads are significantly re For example, an 80.72% reduction in the RMS value of the wing-root bending mo achieved.This is of great significance for the weight reduction design of the aircr GLA results demonstrate that the designed CLLMS-based controllers are very effe

94.14% --
It is noted that a set of time-history curves cannot easily be used to simultan display the correlations of the GLA performances between two physical quantiti issue can be remedied by plotting a convex hull diagram about two physical variab First, the gust-induced response data for a set of physical quantities are compu open-and closed-loop simulations, respectively.Then, the discrete points corresp to two selected physical quantities are plotted on a two-dimensional plane.Finally vex hull is formed by the minimum convex polygon, which encompasses all these d points.The correlations of GLA performance between two selected physical quanti clearly demonstrated in Figure 12.It is noted that a set of time-history curves cannot easily be used to simultaneously display the correlations of the GLA performances between two physical quantities.This issue can be remedied by plotting a convex hull diagram about two physical variables [36].First, the gust-induced response data for a set of physical quantities are computed by openand closed-loop simulations, respectively.Then, the discrete points corresponding to two selected physical quantities are plotted on a two-dimensional plane.Finally, a convex hull is formed by the minimum convex polygon, which encompasses all these discrete points.The correlations of GLA performance between two selected physical quantities are clearly demonstrated in Figure 12.The statistical results in Table 2 indicate that the error signals, RMS values of the accelerations at monitoring nodes, and wing-root internal loads are significantly reduced.For example, an 80.72% reduction in the RMS value of the wing-root bending moment is achieved.This is of great significance for the weight reduction design of the aircraft.The GLA results demonstrate that the designed CLLMS-based controllers are very effective.It is noted that a set of time-history curves cannot easily be used to simultaneously display the correlations of the GLA performances between two physical quantities.This issue can be remedied by plotting a convex hull diagram about two physical variables [36].First, the gust-induced response data for a set of physical quantities are computed by open-and closed-loop simulations, respectively.Then, the discrete points corresponding to two selected physical quantities are plotted on a two-dimensional plane.Finally, a convex hull is formed by the minimum convex polygon, which encompasses all these discrete points.The correlations of GLA performance between two selected physical quantities are clearly demonstrated in Figure 12.

GLA Control in a Discrete Gust Field
When an aircraft travels through the 1-cos gust field with a cruise speed U, the vertical disturbance velocity sensed by the gust reference point on the aircraft is given by where w g is the amplitude of gust velocity, L g is the gust scale, and f g = U/L g is the gust excitation frequency.
The aim of this subsection is to verify whether the converged CLLMS-based controllers trained in the continuous gust environment are also applicable in the case of a discrete gust condition.In simulations, the maximum gust velocity w g and the gust excitation frequency f g are taken as 1 m/s and 3 Hz, respectively.The alleviation results for acceleration and wing-root internal load responses are shown in Figure 13 and Figure 14, respectively.It can be observed that the responses of the flying-wing system are significantly reduced after applying the GLA control.It can be seen from Figure 15 that among all control surfaces, aileron 4 has the largest deflection amplitude because the outermost control surface has the highest effectiveness in suppressing bending vibration response.
where g w is the amplitude of gust velocity, g L is the gust scale, and is the gust excitation frequency.The aim of this subsection is to verify whether the converged CLLMS-based controllers trained in the continuous gust environment are also applicable in the case of a discrete gust condition.In simulations, the maximum gust velocity g w and the gust excitation frequency g f are taken as 1 m/s and 3 Hz, respectively.The alleviation results for accel- eration and wing-root internal load responses are shown in Figures 13 and 14, respectively.It can be observed that the responses of the flying-wing system are significantly reduced after applying the GLA control.It can be seen from Figure 15 that among all control surfaces, aileron 4 has the largest deflection amplitude because the outermost control surface has the highest effectiveness in suppressing bending vibration response.Table 3 gives a summary of the GLA performance in the 1-cos gust condition.It is, thus, clear that the CLLMS-based controllers trained in the continuous gust condition still maintain good performance in the discrete gust environment.Table 3. Performance of the CLLMS-based adaptive feed-forward control system (1-cos gust).

Constraints on Control Surface Deflections
Note that the effort constraints are not involved in the above control algorithm.In practical applications, the total power of the control signal is limited because the actuator has limited driving capability.In addition, when incorporating GLA control into a flight control system, it is necessary to limit the deflection of the control surface.To this end, assume that the deflection angle Let T be the sampling interval.Using the definition of time derivative, Equation (31) can be transformed into the following form: Table 3 gives a summary of the GLA performance in the 1-cos gust condition.It is, thus, clear that the CLLMS-based controllers trained in the continuous gust condition still maintain good performance in the discrete gust environment.

Constraints on Control Surface Deflections
Note that the effort constraints are not involved in the above control algorithm.In practical applications, the total power of the control signal is limited because the actuator has limited driving capability.In addition, when incorporating GLA control into a flight control system, it is necessary to limit the deflection of the control surface.To this end, assume that the deflection angle δ c (t) and the deflection rate .δ c (t) of a control surface are bounded by: where δ cmin and δ cmax represent the minimum and maximum allowed deflections, respectively.ρ cmin and ρ cmax represent the minimum and maximum allowed deflection rates, respectively.Let T be the sampling interval.Using the definition of time derivative, Equation (31) can be transformed into the following form: where In this paper, an algorithm for bounding the controller output is given by the following coefficient update and control output equations: where u c (n) is the controller output, namely the command control surface deflection.
To verify the effectiveness of the above weight-updating rules, first, the deflection angle for each aileron is bounded by ±3.0 deg, and the deflection rate is bounded by ±50 deg/s.From the solid red line in Figures 16 and 17, the algorithm does limit the deflection angle of aileron 1.The algorithm did not affect the deflection rate because it did not exceed the amplitude of the constraint.Constraints on deflection motion will result in the degradation of GLA performance, as illustrated in Figure 18.Table 4 lists the percentage of RMS reduction under the constraint of the amplitude of defection angle.

2.
( ) u n is the controller output, namely the command control surface deflection.
To verify the effectiveness of the above weight-updating rules, first, the deflection angle for each aileron is bounded by 3.0 deg ± , and the deflection rate is bounded by 50 deg/ s ± . From the solid red line in Figures 16 and 17, the algorithm does limit the deflection angle of aileron 1.The algorithm did not affect the deflection rate because it did not exceed the amplitude of the constraint.Constraints on deflection motion will result in the degradation of GLA performance, as illustrated in Figure 18.Table 4 lists the percentage of RMS reduction under the constraint of the amplitude of defection angle.
To verify the effectiveness of the above weight-updating rules, first, the deflection angle for each aileron is bounded by 3.0 deg ± , and the deflection rate is bounded by 50 deg/ s ± . From the solid red line in Figures 16 and 17, the algorithm does limit the deflection angle of aileron 1.The algorithm did not affect the deflection rate because it did not exceed the amplitude of the constraint.Constraints on deflection motion will result in the degradation of GLA performance, as illustrated in Figure 18.Table 4 lists the percentage of RMS reduction under the constraint of the amplitude of defection angle.The statistical results in Table 2 indicate that the error signals, RMS values of the accelerations at monitoring nodes, and wing-root internal loads are significantly reduced.For example, an 80.72% reduction in the RMS value of the wing-root bending moment is achieved.This is of great significance for the weight reduction design of the aircraft.The GLA results demonstrate that the designed CLLMS-based controllers are very effective.It is noted that a set of time-history curves cannot easily be used to simultaneously display the correlations of the GLA performances between two physical quantities.This issue can be remedied by plotting a convex hull diagram about two physical variables [36].First, the gust-induced response data for a set of physical quantities are computed by open-and closed-loop simulations, respectively.Then, the discrete points corresponding to two selected physical quantities are plotted on a two-dimensional plane.Finally, a convex hull is formed by the minimum convex polygon, which encompasses all these discrete points.The correlations of GLA performance between two selected physical quantities are clearly demonstrated in Figure 12.The statistical results in Table 2 indicate that the error signals, RMS values of celerations at monitoring nodes, and wing-root internal loads are significantly re For example, an 80.72% reduction in the RMS value of the wing-root bending mo achieved.This is of great significance for the weight reduction design of the aircr GLA results demonstrate that the designed CLLMS-based controllers are very effe It is noted that a set of time-history curves cannot easily be used to simultan display the correlations of the GLA performances between two physical quantiti issue can be remedied by plotting a convex hull diagram about two physical variab First, the gust-induced response data for a set of physical quantities are compu open-and closed-loop simulations, respectively.Then, the discrete points corresp to two selected physical quantities are plotted on a two-dimensional plane.Finally vex hull is formed by the minimum convex polygon, which encompasses all these d points.The correlations of GLA performance between two selected physical quanti clearly demonstrated in Figure 12.Furthermore, in order to study the GLA performance when the deflation rate reaches saturation, the deflection angle for each aileron is bounded by ±3.0 deg, and the deflection rate is bounded by ±20 deg/s, which means the fast rate of the deflection is constrained.As illustrated in Figures 19 and 20, the algorithm takes account of both amplitude and rate constraints.Figure 21 shows that the performance of the load reduction deteriorates compared with the unsaturation situation of the deflection rate.Table 5 lists the percentages of RMS reduction under the constraints of the amplitude and rate of defection.The reduction is only 28.27% and 31.43% in wing-root bending moment and torsional moment, respectively.It can be seen that the deflection motion constraints have a significant impact on GLA performance.19 and 20, the algorithm takes account of both amplitude and rate constraints.Figure 21 shows that the performance of the load reduction deteriorates compared with the unsaturation situation of the deflection rate.Table 5 lists the percentages of RMS reduction under the constraints of the amplitude and rate of defection.The reduction is only 28.27% and 31.43% in wing-root bending moment and torsional moment, respectively.It can be seen that the deflection motion constraints have a significant impact on GLA performance.The statistical results in Table 2 indicate that the error signals, RMS values of the accelerations at monitoring nodes, and wing-root internal loads are significantly reduced.For example, an 80.72% reduction in the RMS value of the wing-root bending moment is achieved.This is of great significance for the weight reduction design of the aircraft.The GLA results demonstrate that the designed CLLMS-based controllers are very effective.It is noted that a set of time-history curves cannot easily be used to simultaneously display the correlations of the GLA performances between two physical quantities.This issue can be remedied by plotting a convex hull diagram about two physical variables [36].First, the gust-induced response data for a set of physical quantities are computed by open-and closed-loop simulations, respectively.Then, the discrete points corresponding to two selected physical quantities are plotted on a two-dimensional plane.Finally, a convex hull is formed by the minimum convex polygon, which encompasses all these discrete points.The correlations of GLA performance between two selected physical quantities are clearly demonstrated in Figure 12.The statistical results in Table 2 indicate that the error signals, RMS values of celerations at monitoring nodes, and wing-root internal loads are significantly re For example, an 80.72% reduction in the RMS value of the wing-root bending mo achieved.This is of great significance for the weight reduction design of the aircr GLA results demonstrate that the designed CLLMS-based controllers are very effe It is noted that a set of time-history curves cannot easily be used to simultan display the correlations of the GLA performances between two physical quantiti issue can be remedied by plotting a convex hull diagram about two physical variab First, the gust-induced response data for a set of physical quantities are compu open-and closed-loop simulations, respectively.Then, the discrete points corresp to two selected physical quantities are plotted on a two-dimensional plane.Finally vex hull is formed by the minimum convex polygon, which encompasses all these d points.The correlations of GLA performance between two selected physical quanti clearly demonstrated in Figure 12.The actuator has several failure modes, in which the control surface jamming is the most common failure mode.This subsection focuses on the influence of the jamming failure of the outmost control surface (aileron 1) on the GLA performance.
We assume that at the 25th second after GLA control is activated, aileron 1 is jammed at the 2 deg position; see the blue solid line in Figure 22.It can be seen that when aileron 1 fails suddenly, the deflection amplitudes of the other control surfaces gradually increase to compensate for the adverse effects resulting from aileron 1 failure.This indicates that the proposed MIMO adaptive feed-forward control has the self-regulating function for control energy in the case of jamming failure on the control surface.Figure 23 shows the displacement response at node 647 in the case of aileron 1 jamming failure.Before 25 s, the wing-tip displacement is significantly suppressed.Between 25 s and 40 s, the controller begins to adapt to this new fault environment.In this adaptation stage, there is a partial degradation of GLA performance.After 40 s, good GLA performance is restored.Since aileron 1 is jammed at the 2-degree position, a constant static displacement is found, As shown in Figure 23.It can be seen that the utilization of multiple control surfaces for GLA is beneficial for maintaining a stable GLA performance when actuator failure occurs.This is attributed to the robustness brought about by multiple control surfaces and the adaptability of the control algorithm.

Conclusions
To investigate the GLA problem for a flying-wing model with multiple control surfaces, this paper proposed a MIMO adaptive control scheme based on the CLLMS algo-

Conclusions
To investigate the GLA problem for a flying-wing model with multiple control surfaces, this paper proposed a MIMO adaptive control scheme based on the CLLMS algorithm.The simulation results demonstrate that the developed GLA control system can

Figure 1 .
Figure 1.Structural and aerodynamic models of the flying wing.(a) FEM model; (b) aerodynamic model.

Figure 3 .
Figure 3. Locations of the alpha probe and monitoring FEM nodes.

Figure 3 .
Figure 3. Locations of the alpha probe and monitoring FEM nodes.
respectively, are used to generate the filtered reference signals 1 ˆc i r and 2 ˆc i r .The filtered reference signals together with the error signals are used to update the parameters of adaptive controllers.The objective of GLA control is to adjust the controller parameters to minimize the error signals as much as possible.

Figure 4 .
Figure 4. Block diagram of the designed MIMO adaptive feed-forward control system.

Figure 4 .
Figure 4. Block diagram of the designed MIMO adaptive feed-forward control system.

Figure 6 .
Figure 6.Block diagram for MIMO adaptive feed-forward GLA control.Since GLA control always operates below the flutter speed of the system, all numerical simulations are performed at a cruise speed of 27.2 m/s and sea level.The continuous gust scale in the Dryden model introduced in Section 4.3 is set to 53.3 m w L = , and the RMS value of the gust velocity is taken as 0.5 m/s w σ =

Figure 6 .
Figure 6.Block diagram for MIMO adaptive feed-forward GLA control.

Figure 7 .
Figure 7. Simulation results of the Dryden gust.

Figure 7 .
Figure 7. Simulation results of the Dryden gust.

Figure 12 .
Figure 12.Convex hulls about two physical quantities.(a) Wingtip acceleration and bending moment; (b) accelerations at node 647 and node 11098; (c) bending and torsional moments; (d) bending moment and shear force.
and maximum allowed deflection rates, respectively.

Figure 18 .
Figure 18.Wing-root bending moment with constraint on amplitude of defection angle.

Figure 18 .
Figure 18.Wing-root bending moment with constraint on amplitude of defection angle.

Figure 21 .
Figure 21.Wing-root bending moment with constraints on amplitude and rate of defection.

Figure 21 .
Figure 21.Wing-root bending moment with constraints on amplitude and rate of defection.

Figure 21 .
Figure 21.Wing-root bending moment with constraints on amplitude and rate of defection.

6 .
Case of Partial Failure of Actuators

Aerospace 2023 ,
10, x FOR PEER REVIEW 21 of 24is attributed to the robustness brought about by multiple control surfaces and the adaptability of the control algorithm.

Table 1 .
The first six natural frequencies of the clamped wing.

Table 1 .
The first six natural frequencies of the clamped wing.
w is the aerodynamic lag root matrix, and the other matrices are obtained by solving the least squares (LS) problem.Using the MS approximation and introducing the additional aerodynamic states, the state-space open-loop aeroelastic equation can be written as [31] ae = A ae x ae + B ae u ae + B aw w g .x and C 2 are the positive constants satisfying C 1 < C 2 .Note that the CLLMS algorithm requires four constants: µ, γ, C 1 , and C 2 [34].

Table 2 .
Performance of the CLLMS-based adaptive feed-forward control system (Dryden gust).

Table 2 .
Performance of the CLLMS-based adaptive feed-forward control system (Dryden gust).

Table 2 .
Performance of the CLLMS-based adaptive feed-forward control system (Dryden g y M , Nm

Table 2 .
Performance of the CLLMS-based adaptive feed-forward control system (Dryden gust).

Table 4 .
Performance of the adaptive feed-forward control system with constraint on amplitude of defection ( min z F , N 55.08% (66.12%) x M , Nm 58.77% (80.72%)

Table 2 .
Performance of the CLLMS-based adaptive feed-forward control system (Dryden gust).

Table 2 .
Performance of the CLLMS-based adaptive feed-forward control system (Dryden g

Table 4 .
Performance of the adaptive feed-forward control system with constraint on amplitude of defection ( min

Table 5 .
Performance of the adaptive feed-forward control system with constraints on amplitude and rate of defection ( min

Table 5 .
Performance of the adaptive feed-forward control system with constraints on amplitude and rate of defection ( min

Table 2 .
Performance of the CLLMS-based adaptive feed-forward control system (Dryden gust).

Table 2 .
Performance of the CLLMS-based adaptive feed-forward control system (Dryden g