# GPC-Based Gust Response Alleviation for Aircraft Model Adapting to Various Flow Velocities in the Wind Tunnel

Authors: Yuting Dai; Chao Yang
Journal: Shock and Vibration
 (2015)

---

## Abstract

A unified autoregressive (AR) model is identified, based on the wind tunnel test data of open-loop gust response for an aircraft model. The identified AR model can be adapted to various flow velocities in the wind tunnel test. Due to the lack of discrete gust input measurement, a second-order polynomial function is used to approximate the gust input amplitude by flow velocity. Afterwards, with the identified online aeroelastic model, the modified generalized predictive control (GPC) theory is applied to alleviate wing tip acceleration induced by sinusoidal gust. Finally, the alleviation effects of gust response at different flow velocities are estimated based on the comparison of simulated closed-loop acceleration with experimental open-loop one. The comparison indicates that, after gust response alleviation, the wing tip acceleration can be reduced up to 20% at the tested velocities ranging from 12 m/s to 24 m/s. Demonstratively, the unified control law can be adapted to varying wind tunnel velocities and gust frequencies. It does not need to be altered at different test conditions, which will save the idle time.


---

## Body

## 1. Introduction

Dynamic response of aircrafts induced by gust or turbulence may reduce the ride quality, and it increases the structural load [1]. Researchers found that the gust loads can be successfully reduced when gust load alleviation systems are applied to aircrafts [2, 3]. In fact, gust load alleviation (GLA) active control is an effective tool to reduce the dynamic gust response, with a minor increasing of aircraft’s weight. For example, they are used in aircraft such as the B-52 and C5 [4, 5]. In order to guarantee real aircraft’s safety, gust response alleviation systems are designed and validated using many wind-tunnel tests [6, 7]. The control law design is an important part in the design process of GLA active technique. Most researches about GLA control law design were concentrated on PID method and linear quadratic Gauss (LQG) theory [8]. H∞ optimal control and μ synthesis are also effective robust control methods to account for variations in the mathematical model. However, all those control laws are based on a known theoretical continuous aeroelastic model. Model-based control laws may fail when the aeroelastic model is not accurate enough. Alternatively, a data-based control law design may be suitable for gust load alleviation. The data-based autoregressive model is a good tool to construct a discrete model without the theoretical continuous mathematical model [9]. It is found to be effective in aeroelastic modeling for online flutter prediction [10]. Therefore, the controller design can be conducted on a data-based AR model. This is the significant advantage of generalized predictive control (GPC) method. GPC can tackle not only with the theoretical continuous model but also for discrete AR model. Notably, it is demonstrated to be a useful controller for linear-parameter-varying nonlinear system [11]. Based on the above advantages, the data-based GPC is welcome to aeroelastic active control, for both gust load alleviation and flutter suppression [12, 13]. It is validated to be effective for gust load alleviation in the simulation [12].Wu et al. conducted a gust response wind tunnel test and they designed a PID controller to alleviate the wing tip acceleration. The open-loop response data was measured at a specific velocity under sinusoidal gust. After all the open-loop responses are measured ranging from 12 m/s to 24 m/s, the closed-loop responses are also measured at these conditions by a PID controller. In this test, there are two challenged problems. Because there are some discrepancies of theoretical aeroelastic model and real aircrafts [14], the best PID parameters on theoretical model may not work well on real wind tunnel test model. They have to be trialed several times in the wind tunnel test. The other problem is that the amplitude of discrete sinusoidal gust disturbance was not measured in this test, which makes the comparison between theoretical results and testing ones difficult. Hence, the data-based GPC is designed to the wind gust response alleviation test. The basis for generalized predictive control is the identification of an autoregressive mode. While the gust input used for AR model identification is unknown, it varies with test velocity. Moreover, if we identify one AR model at one test condition, we will need to alter the control law for other different test conditions. It will be a waste of time to switch the control law manually in the wind tunnel test.Hence, a unified GPC controller is developed. It is adapting to all the test conditions, used for gust response alleviation in the wind tunnel test. The sections of this paper are as follows. First, the standard GPC control law is introduced, and then it is modified to adapt to varying wind tunnel test velocities. Finally, a wind tunnel test of an aircraft model is employed to validate the alleviation effect.


## 2. Control Law Design for Gust Response Alleviation

In this section, firstly, the standard GPC design method for gust response alleviation is derived at a specific flow velocity and at a specific gust frequency. Afterwards, the standard GPC method is modified to adapt to varying flow velocities.


### 2.1. Gust Response Alleviation at a Fixed Flow Velocity

GPC is a data-based method for control law design. The design of it begins by identifying an autoregressive (AR) model, based on open-loop input and output data [12]. After the AR model is identified, a controller is acted on it to minimize the prediction of the system response in the future. In this design process, the control law is not represented by a state-space equation form but is described as a sequence of discrete input-output data.When an external excitation exists, a time-invariant multi-input-multi-output AR model is written as [12](1)yk=α1yk-1+α2yk-2+⋯+αpyk-p+β0uk+β1uk-1+⋯+βpuk-p+γ0dk+γ1dk-1+⋯+γpdk-p,where integer p is called the order of the AR model. y(k-p) is the output of the model and d(k-p) is the measured disturbance. Since the gust response in the wind tunnel test is excited by a sinusoidal-moving gust generator, the level of signal-to-noise ratio is large enough to omit noise d(k-p). u(k-p) is the system input at the time of k-p. It contains two parts. One is the deflection of control surfaces uc, the other is the gust velocity on the aircraft model in the wind tunnel test. As mentioned above, the discrete gust is generated by a biplane which moves as a sinusoidal function. Hence, the gust input ug has the following form:(2)ug=Avsin⁡ωt,where Av is the gust amplitude; it is related with flow velocity. ω is the known gust frequency.Substituting (2) to (1) and dividing the input signal u into two parts, uc and ug, then (1) can be modified as(3)yk=α1yk-1+α2yk-2+⋯+αpyk-p+β0′uck+β1′uck-1+⋯+βp′uck-p+γ10Avsin⁡ωtk+γ11Avsin⁡ωtk-1+⋯+γ1pAvsin⁡ωtk-p.When we keep the flow velocity fixed under one test condition, the gust amplitude can be regarded as a constant. Rearranging (3), we can get(4)yk=α1yk-1+α2yk-2+⋯+αpyk-p+β0′uck+β1′uck-1+⋯+βp′uck-p+γ0′udk+γ1′udk-1+⋯+γp′udk-p,where γ0′=Avγ10 and ud(k)=sin⁡ωt(k). By this rearranging, the gust disturbance ud is known to us. Hence, when the responses of wing tip acceleration and deflections of control surfaces uc are measured, the regressive relationship of input and output in the l time steps can be written as(5)Y=Y-V,where Y=[y(0),y(1),…,y(l-1)] and V is formed of the discrete time series of uc, ud, and y. It is written as(6)V=ud0ud1ud2⋯udp⋯udl-1uc0uc1uc2⋯ucp⋯ucl-10v0v1⋯vp-1⋯vl-200v0⋯vp-2⋯vl-3⋱⋯⋯⋮v0⋯vl-p-1,where the matrix element v(k) is the combination of gust disturbance signal, input signal, and output signal at the kth time step. That is, vk=[ud(k),uc(k),y(k)]T. Y- is the vector of observer Markov parameters to be identified. Associating Y- with (4), we can get(7)Y-=γ0′β0′γ1′β1′α1γ2′β2′α2⋯γp′βp′αp′.The solution of Y- is calculated by employing a least square algorithm. That is [12],(8)Y-=yV+=yVTVVT-1.After the parameter Y- of the AR model at the kth time step in (4) is identified according to the past l time steps, the control system is switched on from the k+1 time step. It drives the control surface to move. In the closed-loop system, the deflection of control surface is still denoted as uc. Then in the future hp time steps, the response of the future time step jth can be also represented as a linear combination of three parts. One part is the response of future j time steps and the last p time steps. The second part is the deflection of control surface uc in the future hc time steps and in the last p time steps. The third part is the gust disturbance in the future j time steps and the past p time steps. The regressive relationship of future response is written as(9)yk+j=α1jyk-1+j+α2jyk-2+j+⋯+αpjyk-p+j+β0′uck+j+β0′1uck-1+j+⋯+β0′juck+β1′juck-1+⋯+βp′juck-p+γ0′udk+j+γ0′1udk-1+j+⋯+γ0′judk+γ1′judk-1+⋯+γp′judk-p.In the above equation, the regressive coefficients are represented by combinations of the observer Markov parameters in (7). That is,elements inA: αpj=αp+1j-1+α1j-1αp,elements inT and B: β0′j=β1′j-1+α1j-1β0′, βp′j=βp+1′j-1+α1j-1βp′,elements inDf and Dp: γ0′j=γ1′j-1+α1j-1γ0′, γp′j=γp+1′j-1+α1j-1γp′.The goal for control law design is to alleviate the gust response in the wind tunnel test. Hence, in the framework of GPC strategy, it is required to minimize the predicted wing tip acceleration in the future hp time steps. Hence, the objective of the optimized control command in the next hc time steps is written as [12] (10)min⁡yhpTkRyhpk+uhcTkQuhck,where(11)yhpk=ykyk+1⋮yk+hp-1,uhck=uckuck+1⋮uck+hc-1,dgfk=udkudk+1⋮udk+hp-1.It is a typical optimal control problem. Equation (10) means that we want to find a control command to obtain a weighted minimized response in the future hp time steps, with relatively low control energy. By solving the optimal problem shown in (10), the optimal deflection of control surface at the time step k+1 is written as(12)uck+1=-TTRT+Q†TTRBupk-p+Aypk-p+Dpdgpk-p+Dfdgfkrc,where T, B, A, Dp, and Df can be calculated in advance and be written into the control program for online gust load alleviation. R and Q are the weighting matrices for predicted response and control effort, respectively. In the wind tunnel test, the wing tip acceleration has to be alleviated. And the aileron and elevator are selected as the control surfaces for gust alleviation.


### 2.2. Gust Response Alleviation at Varying Flow Velocities

The basis of standard GPC is to identify an AR model according to the known input-output data. Hence, when GPC is employed to gust response alleviation, it is assumed that the sinusoidal gust input is already known, not only its frequency, but also its amplitude. The gust frequency is easy to get from measuring the angular velocity of the gust generator. When the flow velocity in the wind tunnel is fixed, the value of gust amplitude can be merged to the identified coefficientsγp′, as shown in (4). However, when the flow velocity varies, the corresponding gust input has been changed. In this case, another AR model should be identified. In this case, a different control law should be redesigned for this different test condition. Therefore, we have to switch the control law frequently from one test condition to another. This may waste a lot of time in a real wind tunnel test. The idea in this paper is to identify a unified AR model to apply only one GPC controller to all the test velocities in the wind tunnel.As mentioned above, the gust velocity in the wind tunnel test is not measured. Moreover, it is difficult to deduce the gust velocity from the comparison between theoretical response and experimental one. Hence, the relationship between gust velocities and flow velocities is unknown for us. In order to solve this problem, the relationship between gust input and flow velocity is approximated in this current work.A direct way is to employ a fitting function to represent the relationship between gust input and flow velocity. For the sake of simplicity, a second-order polynomial function is applied to represent it by the linear combinations of flow velocities. That is,(13)ug=A0sin⁡ωit+A1V0sin⁡ωit+A2V02sin⁡ωit,where ωi is the gust frequency at different test conditions. To combine this expression to the standard AR model, we increase the dimensions of gust disturbance from one to three. Therefore, the gust disturbance in (4) is written as a vector:(14)udk=sin⁡ωitk,V0sin⁡ωitk,V02sin⁡ωitkT.By this approximation, we can represent the gust disturbance by varying flow velocities. The predicted control law also needs to be modified to adjust to three gust inputs, shown in (12). Hence, a unified control law can be adapted to different flow velocities, and there is no need to switch the control laws from one test condition to another.


## 2.1. Gust Response Alleviation at a Fixed Flow Velocity

GPC is a data-based method for control law design. The design of it begins by identifying an autoregressive (AR) model, based on open-loop input and output data [12]. After the AR model is identified, a controller is acted on it to minimize the prediction of the system response in the future. In this design process, the control law is not represented by a state-space equation form but is described as a sequence of discrete input-output data.When an external excitation exists, a time-invariant multi-input-multi-output AR model is written as [12](1)yk=α1yk-1+α2yk-2+⋯+αpyk-p+β0uk+β1uk-1+⋯+βpuk-p+γ0dk+γ1dk-1+⋯+γpdk-p,where integer p is called the order of the AR model. y(k-p) is the output of the model and d(k-p) is the measured disturbance. Since the gust response in the wind tunnel test is excited by a sinusoidal-moving gust generator, the level of signal-to-noise ratio is large enough to omit noise d(k-p). u(k-p) is the system input at the time of k-p. It contains two parts. One is the deflection of control surfaces uc, the other is the gust velocity on the aircraft model in the wind tunnel test. As mentioned above, the discrete gust is generated by a biplane which moves as a sinusoidal function. Hence, the gust input ug has the following form:(2)ug=Avsin⁡ωt,where Av is the gust amplitude; it is related with flow velocity. ω is the known gust frequency.Substituting (2) to (1) and dividing the input signal u into two parts, uc and ug, then (1) can be modified as(3)yk=α1yk-1+α2yk-2+⋯+αpyk-p+β0′uck+β1′uck-1+⋯+βp′uck-p+γ10Avsin⁡ωtk+γ11Avsin⁡ωtk-1+⋯+γ1pAvsin⁡ωtk-p.When we keep the flow velocity fixed under one test condition, the gust amplitude can be regarded as a constant. Rearranging (3), we can get(4)yk=α1yk-1+α2yk-2+⋯+αpyk-p+β0′uck+β1′uck-1+⋯+βp′uck-p+γ0′udk+γ1′udk-1+⋯+γp′udk-p,where γ0′=Avγ10 and ud(k)=sin⁡ωt(k). By this rearranging, the gust disturbance ud is known to us. Hence, when the responses of wing tip acceleration and deflections of control surfaces uc are measured, the regressive relationship of input and output in the l time steps can be written as(5)Y=Y-V,where Y=[y(0),y(1),…,y(l-1)] and V is formed of the discrete time series of uc, ud, and y. It is written as(6)V=ud0ud1ud2⋯udp⋯udl-1uc0uc1uc2⋯ucp⋯ucl-10v0v1⋯vp-1⋯vl-200v0⋯vp-2⋯vl-3⋱⋯⋯⋮v0⋯vl-p-1,where the matrix element v(k) is the combination of gust disturbance signal, input signal, and output signal at the kth time step. That is, vk=[ud(k),uc(k),y(k)]T. Y- is the vector of observer Markov parameters to be identified. Associating Y- with (4), we can get(7)Y-=γ0′β0′γ1′β1′α1γ2′β2′α2⋯γp′βp′αp′.The solution of Y- is calculated by employing a least square algorithm. That is [12],(8)Y-=yV+=yVTVVT-1.After the parameter Y- of the AR model at the kth time step in (4) is identified according to the past l time steps, the control system is switched on from the k+1 time step. It drives the control surface to move. In the closed-loop system, the deflection of control surface is still denoted as uc. Then in the future hp time steps, the response of the future time step jth can be also represented as a linear combination of three parts. One part is the response of future j time steps and the last p time steps. The second part is the deflection of control surface uc in the future hc time steps and in the last p time steps. The third part is the gust disturbance in the future j time steps and the past p time steps. The regressive relationship of future response is written as(9)yk+j=α1jyk-1+j+α2jyk-2+j+⋯+αpjyk-p+j+β0′uck+j+β0′1uck-1+j+⋯+β0′juck+β1′juck-1+⋯+βp′juck-p+γ0′udk+j+γ0′1udk-1+j+⋯+γ0′judk+γ1′judk-1+⋯+γp′judk-p.In the above equation, the regressive coefficients are represented by combinations of the observer Markov parameters in (7). That is,elements inA: αpj=αp+1j-1+α1j-1αp,elements inT and B: β0′j=β1′j-1+α1j-1β0′, βp′j=βp+1′j-1+α1j-1βp′,elements inDf and Dp: γ0′j=γ1′j-1+α1j-1γ0′, γp′j=γp+1′j-1+α1j-1γp′.The goal for control law design is to alleviate the gust response in the wind tunnel test. Hence, in the framework of GPC strategy, it is required to minimize the predicted wing tip acceleration in the future hp time steps. Hence, the objective of the optimized control command in the next hc time steps is written as [12] (10)min⁡yhpTkRyhpk+uhcTkQuhck,where(11)yhpk=ykyk+1⋮yk+hp-1,uhck=uckuck+1⋮uck+hc-1,dgfk=udkudk+1⋮udk+hp-1.It is a typical optimal control problem. Equation (10) means that we want to find a control command to obtain a weighted minimized response in the future hp time steps, with relatively low control energy. By solving the optimal problem shown in (10), the optimal deflection of control surface at the time step k+1 is written as(12)uck+1=-TTRT+Q†TTRBupk-p+Aypk-p+Dpdgpk-p+Dfdgfkrc,where T, B, A, Dp, and Df can be calculated in advance and be written into the control program for online gust load alleviation. R and Q are the weighting matrices for predicted response and control effort, respectively. In the wind tunnel test, the wing tip acceleration has to be alleviated. And the aileron and elevator are selected as the control surfaces for gust alleviation.


## 2.2. Gust Response Alleviation at Varying Flow Velocities

The basis of standard GPC is to identify an AR model according to the known input-output data. Hence, when GPC is employed to gust response alleviation, it is assumed that the sinusoidal gust input is already known, not only its frequency, but also its amplitude. The gust frequency is easy to get from measuring the angular velocity of the gust generator. When the flow velocity in the wind tunnel is fixed, the value of gust amplitude can be merged to the identified coefficientsγp′, as shown in (4). However, when the flow velocity varies, the corresponding gust input has been changed. In this case, another AR model should be identified. In this case, a different control law should be redesigned for this different test condition. Therefore, we have to switch the control law frequently from one test condition to another. This may waste a lot of time in a real wind tunnel test. The idea in this paper is to identify a unified AR model to apply only one GPC controller to all the test velocities in the wind tunnel.As mentioned above, the gust velocity in the wind tunnel test is not measured. Moreover, it is difficult to deduce the gust velocity from the comparison between theoretical response and experimental one. Hence, the relationship between gust velocities and flow velocities is unknown for us. In order to solve this problem, the relationship between gust input and flow velocity is approximated in this current work.A direct way is to employ a fitting function to represent the relationship between gust input and flow velocity. For the sake of simplicity, a second-order polynomial function is applied to represent it by the linear combinations of flow velocities. That is,(13)ug=A0sin⁡ωit+A1V0sin⁡ωit+A2V02sin⁡ωit,where ωi is the gust frequency at different test conditions. To combine this expression to the standard AR model, we increase the dimensions of gust disturbance from one to three. Therefore, the gust disturbance in (4) is written as a vector:(14)udk=sin⁡ωitk,V0sin⁡ωitk,V02sin⁡ωitkT.By this approximation, we can represent the gust disturbance by varying flow velocities. The predicted control law also needs to be modified to adjust to three gust inputs, shown in (12). Hence, a unified control law can be adapted to different flow velocities, and there is no need to switch the control laws from one test condition to another.


## 3. Application Example

### 3.1. The Open-Loop Gust Response Wind Tunnel Test

In 2011, the gust response and alleviation test for the half-span aircraft model was conducted in the FD-09 Wind Tunnel [8]. The dimension of the wind tunnel is 3 × 3 m in the test section. The pitch-plunge supporting mechanism was mounted below the wind tunnel floor without disturbance to the airflow. A feedback control law is acting on the actuator of the elevator, to stabilize the pitch and plunge motions. The responses are measured by the NI PXI-4472B NI data acquisition equipment, with the acquisition frequency of 200 Hz. The mounted aircraft model and gust response alleviation flowchart are shown in Figure 1. In this test, the open-loop gust response is measured first. And then a PI control law is acted on the aileron for gust load alleviation.Figure 1
The elastic aircraft model and experimental control flowchart.In this current work, the open-loop experimental gust response is used for AR model identification. Afterwards, the developed GPC control law is designed based on the AR model and it is switched on after several time steps. The simulated closed-loop gust response is compared with the experimental open-loop one.In order to guarantee the aircraft model’s safety, the experimental flow velocity is from 12 m/s to 24 m/s, much far away from the flutter velocity of 31 m/s. The frequency of gust generator changes from 1.0 Hz to 6.0 Hz with a fixed increasing step of 0.5 Hz. In the open-loop gust response test, the gust generator upwards the aircraft model is switched on all the time. After the gust frequency keeps still, the flow velocity increases to a given value with a fixed increasing step of 2 m/s. At a fixed velocity and gust frequency, we measured the accelerations and deflections of control surfaces for 15 seconds. A Fourier transformation is conducted on all the 15-second data to obtain the response amplitude in the frequency domain. The amplitude of the open-loop wing tip acceleration at different velocities and frequencies is shown in Figure2 and Table 1.Table 1
The open-loop wing tip acceleration data due to sinusoidal gust in the wind tunnel test.
Frequency: Hz
Acceleration at 12 m/s : g
Acceleration at 14 m/s : g
Acceleration at 16 m/s : g
Acceleration at 18 m/s : g
Acceleration at 20 m/s : g
Acceleration at 22 m/s : g
Acceleration at 24 m/s : g
1.0
0.10
0.11
0.12
0.17
0.17
0.18
0.18
1.5
0.21
0.35
0.29
0.28
0.29
0.40
0.42
2.0
0.70
0.69
0.79
0.70
0.79
0.83
0.86
2.5
0.55
0.47
0.80
1.09
0.89
1.30
0.91
3.0
0.42
0.61
0.85
1.28
1.34
1.56
1.82
3.5
0.35
0.37
0.54
1.26
1.24
1.28
1.69
4.0
0.39
0.46
0.73
1.23
1.40
1.80
2.16
4.5
0.32
0.44
0.61
1.06
0.91
1.28
1.82
5.0
0.34
0.42
0.70
0.86
1.32
1.52
1.93
5.5
0.27
0.34
0.49
1.36
0.99
1.49
2.16
6.0
0.33
0.41
0.53
1.41
1.23
1.97
2.48Figure 2
The open-loop wing tip acceleration due to sinusoidal gust in the wind tunnel test.Figure3 is the closed-loop acceleration by the PI controller in the wind tunnel test at the velocity of 18 m/s. From the comparison, it can be seen that the controller cannot suppress the wing tip acceleration at a low gust frequency of 1 Hz or at a high gust frequency of 5 Hz. It works even worse at frequencies higher than 5 Hz. The controller alleviates the wing tip acceleration well at other gust frequencies.Figure 3
Comparison of the experimental open-loop gust response and the closed-loop one by PI controller.


### 3.2. Gust Response Alleviation at Fixed Flow Velocity and Gust Frequency

Gust response alleviation by GPC is simulated at a fixed flow velocity and fixed gust frequency. First, the open-loop wing tip acceleration at the flow velocity of 24 m/s and gust frequency of 2 Hz is employed to identify an AR model. Then, a GPC controller at this specified test condition is designed to alleviate the wing tip response. As mentioned above, at a fixed flow velocity, though the gust input is not measured, it can be given as a known form. That is,ud1=5sin⁡(4πt). Some parameters for model identification are set as below: the input-output data length for identification is 400. The order of AR model p is 5. The number of control inputs is 2. The number of response is one and the number of gust input is one. The weighting matrices for controller design are set to be unity ones. Considering the feasibility of controllers in the wind tunnel test, the lengths of prediction steps and control command steps should not to be too long. Different prediction steps are attempted to choose a better one. Three cases for prediction steps are considered. Their corresponding alleviation results are shown in Table 2. From this table, more prediction steps result in better alleviation results. However, the lengths of prediction window and controller window should not to be too large. It may make controller complicated in the real wind tunnel test. Hence, in the following controller design process, the length of prediction steps hp is fixed to be 5, and the length of control steps hc is also fixed to be 5.Table 2
Alleviation result for different length of prediction steps.
Case number
Length of prediction and control steps
Open-loop acceleration: g
Closed-loop acceleration: g
Alleviation factor
Case  1
h
p
=
1, hc=1
0.858
0.823
4.1%
Case  2
h
p
=
3, hc=3
0.858
0.371
56.8%
Case  3
h
p
=
5, hc=5
0.858
0.11
87.2%The gust response alleviation result for Case 3 is shown in Figure4. The response data in the first 2 seconds is used for AR model identification. After 2 seconds, the prediction controller is switched on. From the figure, it can be seen that the output of the AR model is nearly the same as the experimental acceleration. It indicates that the identified model can represent the aeroelastic behavior at this test condition, and the AR model has a high accuracy for generalized prediction control. After the prediction controller is switched on, the acceleration response reduces to a relatively small value of 0.11 g in a short time. The response comparison in the frequency domain is shown in Figure 5, and the deflections of control surfaces are shown in Figure 6. The deflection of aileron has a shock on the switch-on time step and then it oscillates between −0.1 degrees and 0.1 degrees in the later time. The experimental deflection of elevator before 2 seconds exceeds 12 degrees. The reason is as follows: the elevator is used for stabilizing the aircraft model in the wind tunnel test, and the model is up-mounted which is unstable at the equilibrium point. Hence, its deflection is a little large before the time of 2 seconds. In the closed-loop simulation, it is only used for active gust response alleviation. Hence, after the time of 2 seconds, its deflection remains oscillating between −0.1 degrees and 0.1 degrees, similar to aileron.Figure 4
The comparison of open-loop and closed-loop accelerations in the time domain forhp=5, hc=5.Figure 5
The comparison of open-loop and closed-loop accelerations in the frequency domain.Figure 6
Deflections of control surfaces.
(a)
Elevator
(b)
AileronThe open-loop responses are measured at each fixed gust frequency and each fixed flow velocity. Hence, the AR model is identified straightforward at each test condition too. The control input in (12) is updated with corresponding test conditions to match the current identified AR model. Therefore, there are 77 AR models and 77 discrete control laws by the standard GPC method. The data length for model identification is set to be 2000 in this case. A Fourier transformation is conducted on the closed-loop response to get the response amplitude. The closed-loop response amplitude at each test condition is shown in Table 3.Table 3
Gust response alleviation result at fixed flow velocity and gust frequency.
Frequency: Hz
Acceleration at 12 m/s : g
Acceleration at 14 m/s : g
Acceleration at 16 m/s : g
Acceleration at 18 m/s : g
Acceleration at 20 m/s : g
Acceleration at 22 m/s : g
Acceleration at 24 m/s : g
1.0
0.10
0.11
0.11
0.15
0.16
0.15
0.12
1.5
0.21
0.23
0.23
0.25
0.20
0.15
0.27
2.0
0.68
0.70
0.65
0.74
∗
0.74
0.42
0.59
2.5
0.38
0.46
0.61
0.50
0.80
0.88
0.63
3.0
0.33
0.61
0.82
1.15
0.82
0.89
0.57
3.5
0.27
0.35
0.53
0.70
0.40
0.45
0.86
4.0
0.38
0.36
0.71
1.26
∗
1.12
1.52
1.09
4.5
0.21
0.30
0.41
0.70
0.82
0.60
0.57
5.0
0.34
0.37
0.70
0.88
∗
0.62
0.37
0.55
5.5
0.18
0.27
0.48
0.40
0.80
0.78
0.30
6.0
0.31
0.37
0.38
0.44
0.51
0.81
0.91When the closed-loop response is larger than the open-loop one, it is denoted as bold star in Table3. From the comparison of Tables 1 and 3, it can be seen that the standard GPC controller suppresses the acceleration response at most test conditions. However, it is not always effective at all test conditions. At some frequencies of 18 m/s, the closed-loop response is even larger than the open one. When the data length for identification is only 400, the alleviation effect is even worse. From the alleviation comparison, the standard GPC controller does not work much better than the PI controller at some low frequencies. Moreover, it needs to identify a new AR model when either the flow velocity or gust frequency changes. This may waste a lot of time to switch the controllers frequently in real wind tunnel test.


### 3.3. Gust Response Alleviation at Varying Frequencies

In the above section, we need to identify 77 AR models from one test condition to another. In order to reduce the number of identification models, we need to modify the standard GPC design process. A straightforward way is to merge input-output data of several test conditions to an entire data series. Then we can identify one unified AR model with these merged data. The difficulty is how to tackle with the varying gust input at different test conditions. From (2), two factors may affect the gust velocity. One is the amplitude; the other is the gust frequency. The gust amplitude may keep constant at a fixed flow velocity. Hence, in this section, we merged the test data of all frequencies at a fixed-flow-velocity condition. Consequently, the gust input series is written as(15)ud=sin⁡ω1t1⋯sin⁡ω1tLsin⁡ω2t1⋯sin⁡ω11t1⋯sin⁡ω11tL.Together with the merged experimental acceleration response, the above gust input is applied to identify one unified AR model at a specified flow velocity. By this approach, there will be 7 AR models for 7 flow velocities from 12 m/s to 24 m/s. At each flow velocity, each AR model is employed to design a GPC control law. The identification and controller parameters are set as follows:L=2000, p=5, hp=5, hc=5, and R and Q are unity matrices. Afterwards, each control law is acted on the aircraft at corresponding test velocity. The results are shown in Table 4.Table 4
Gust response alleviation result for varying frequencies and fixed velocity.
Frequency: Hz
Acceleration at 12 m/s : g
Acceleration at 14 m/s : g
Acceleration at 16 m/s : g
Acceleration at 18 m/s : g
Acceleration at 20 m/s : g
Acceleration at 22 m/s : g
Acceleration at 24 m/s : g
1.0
0.08
0.20
∗
0.06
0.27
∗
0.20
∗
0.17
0.14
1.5
0.06
0.17
0.05
0.22
0.15
0.13
0.12
2.0
0.10
0.32
0.09
0.42
0.26
0.25
0.24
2.5
0.08
0.25
0.07
0.33
0.20
0.20
0.19
3.0
0.14
0.41
0.15
0.61
0.34
0.39
0.37
3.5
0.10
0.21
0.14
0.47
0.26
0.33
0.31
4.0
0.18
0.11
0.31
0.77
0.46
0.74
0.64
4.5
0.13
0.05
0.20
0.41
0.33
0.80
0.61
5.0
0.21
0.21
0.22
0.11
0.49
1.58
∗
0.82
5.5
0.15
0.21
0.12
0.31
0.38
0.89
0.45
6.0
0.24
0.39
0.16
0.77
0.67
1.14
0.70In this case, the deflection of aileron is less than 2 degrees at the whole test condition. The bold star mark denotes that the closed-loop response is larger than that of the open-loop one. From this table, it can be seen that the controller does not alleviate the wing tip response at some low frequencies. It may be caused by the turbulence at low frequencies. From Figure3, in the real wind tunnel test, the PI controller does not alleviate gust response at these frequencies, either. At the low or high frequencies, we have to develop other excellent control methods. Note that the data length for identification may influence the alleviation effect.


### 3.4. Gust Response Alleviation at Varying Velocities

In the above section, we identify a unified AR model at a specified flow velocity. Hence, we need to use 7 different AR models for GPC design on the total 77 test conditions. This will save the switch time in the real test. In this section, we still want to merge the discrete data into one data segment at several test conditions. Other than identifying AR models for varying frequencies, alternatively, we identify the AR model for varying velocities at the same gust frequency. In this case, the amplitude of gust disturbance varies with flow velocities. Pointed out in Section2.2, it is represented by the polynomials of flow velocity, and the gust frequency is fixed for a specific AR model. At a specific gust frequency ω0, the gust input is written as(16)udk=sin⁡ω0tk,Visin⁡ω0tk,Vi2sin⁡ω0tkT,where Vi is the flow velocity; it ranges from 12 m/s to 24 m/s. A significant change on the standard GPC method is to increase the dimension of gust disturbance from one to three. Together with merged experimental response, the specific identified AR model can be adaptive to varying velocities. There are total of 11 different gust frequencies in the test. Hence, we have to identify 11 different AR models. By this approach, at each gust frequency in the wind tunnel test, we use the same AR model for gust alleviation. Therefore, we do not need to change the parameters of T, B, A, Dp, and Df at the same gust frequency. It may save the experimental time of switching controllers at different test conditions. The length of each data segment for identification is chosen as 2000. Simulation results indicate that this data length will result in better alleviation effect. The other parameters are set the same as above. Table 5 shows the gust response alleviation results at each test condition for varying velocities. The bold star infers that the closed loop response is not alleviated by the controller.Table 5
Gust response alleviation result for fixed frequency and varying velocities.
Frequency: Hz
Acceleration at 12 m/s : g
Acceleration at 14 m/s : g
Acceleration at 16 m/s : g
Acceleration at 18 m/s : g
Acceleration at 20 m/s : g
Acceleration at 22 m/s : g
Acceleration at 24 m/s : g
1.0
0.13
∗
0.07
0.03
0.01
0.03
0.07
0.14
1.5
0.13
0.12
0.09
0.06
0.07
0.14
0.24
2.0
0.56
0.16
0.24
0.47
0.62
0.70
0.73
2.5
0.20
0.11
0.16
0.19
0.18
0.13
0.19
3.0
0.41
0.23
0.35
0.45
0.68
1.21
2.02
∗
3.5
0.26
0.05
0.05
0.05
0.14
0.40
0.78
4.0
0.22
0.11
0.27
0.25
0.14
0.40
0.94
4.5
0.28
0.04
0.16
0.21
0.22
0.40
0.77
5.0
0.03
0.19
0.26
0.29
0.40
0.67
1.10
5.5
0.20
0.19
0.32
0.33
0.21
0.20
0.57
6.0
0.30
0.32
0.28
0.19
0.30
0.68
1.22From this table, it can be seen that this controller suppresses most of the response. However, the closed-loop response is larger than the open-loop one at 2 test conditions. The amplified closed-loop response mainly occurs at very low and very high frequencies, similar to the results by the PI controller. At these test conditions, the gust frequency is far away from its frequency of first bending mode, which is about 3 Hz.


### 3.5. Gust Alleviation at Varying Flow Velocities and Gust Frequencies

In this section, we merged the whole input-output data to one data series to identify only one AR model. There are total of seven different flow velocities. And, at each velocity, there are 11 different frequencies. Hence, a total 77 varying open-loop responses are used to identify a unified AR model. The first 0.6 seconds of each data segment is chosen to compose the total 46.2 seconds open-loop response data. A fifth-order AR model with 3 gust inputs is identified by the linear square algorithm. The comparison of experimental open-loop response and identified acceleration is shown in Figure7. Figure 7(a) shows the response comparison for lower frequency at varying flow velocity. Figure 7(b) shows the response comparison for higher frequency. From the comparison, it can be seen that though the order of AR model is only five, the experimental response can also agree well with the identified one under varying flow velocities. The identified model can describe the aircraft model’s aeroelastic behavior well. Noting that there is some disturbance when the gust frequency is only 1 Hz, the response of identified model in the low frequency may have some discrepancies from the experimental one, indicated in Figure 7(a). Noting that alleviation factor may fluctuate with the length of each data segment, 0.6-second data here is enough to get a not bad alleviation result.Figure 7
The comparison between the experimental response and the output of identified AR model.
(a)
Response at the first 5 seconds
(b)
Response at the last 5 secondsAfter a unified AR model is identified using the 77 pieces of test data, the modified GPC controller is applied to gust response alleviation at all the experimental conditions. The amplitude of closed-loop response is shown in Table6. From this table, it can be seen that the modified GPC controller is effective at all the test conditions. Comparing it with the open-loop response shown in Table 1, the alleviation factor can be up to 20% for all the test conditions. The closed-loop response and aileron’s deflection are shown in Figure 8. It is compared with the experimental one at the velocity of 24 m/s and the gust frequency of 6 Hz. Though there is a shock at the time of controller acting on the aircraft model, the deflection does not reach 0.5 degrees all the time. It needs only a little control effort to alleviate the wing tip acceleration.Table 6
Gust response alleviation result for varying frequencies and varying velocities.
Frequency: Hz
Acceleration at 12 m/s : g
Acceleration at 14 m/s : g
Acceleration at 16 m/s : g
Acceleration at 18 m/s : g
Acceleration at 20 m/s : g
Acceleration at 22 m/s : g
Acceleration at 24 m/s : g
1.0
0.03
0.01
0.01
0.02
0.05
0.09
0.14
1.5
0.02
0.01
0.01
0.02
0.03
0.06
0.09
2.0
0.03
0.02
0.02
0.03
0.05
0.09
0.14
2.5
0.02
0.01
0.01
0.02
0.03
0.06
0.09
3.0
0.04
0.03
0.03
0.03
0.05
0.09
0.15
3.5
0.03
0.02
0.02
0.02
0.03
0.07
0.12
4.0
0.05
0.04
0.04
0.03
0.05
0.12
0.23
4.5
0.04
0.03
0.03
0.02
0.04
0.10
0.20
5.0
0.07
0.05
0.06
0.03
0.07
0.20
0.39
5.5
0.05
0.04
0.05
0.02
0.05
0.16
0.31
6.0
0.07
0.06
0.08
0.03
0.08
0.27
0.52Figure 8
The closed-loop response (a) and the deflection of aileron (b).
(a)
(b)


## 3.1. The Open-Loop Gust Response Wind Tunnel Test

In 2011, the gust response and alleviation test for the half-span aircraft model was conducted in the FD-09 Wind Tunnel [8]. The dimension of the wind tunnel is 3 × 3 m in the test section. The pitch-plunge supporting mechanism was mounted below the wind tunnel floor without disturbance to the airflow. A feedback control law is acting on the actuator of the elevator, to stabilize the pitch and plunge motions. The responses are measured by the NI PXI-4472B NI data acquisition equipment, with the acquisition frequency of 200 Hz. The mounted aircraft model and gust response alleviation flowchart are shown in Figure 1. In this test, the open-loop gust response is measured first. And then a PI control law is acted on the aileron for gust load alleviation.Figure 1
The elastic aircraft model and experimental control flowchart.In this current work, the open-loop experimental gust response is used for AR model identification. Afterwards, the developed GPC control law is designed based on the AR model and it is switched on after several time steps. The simulated closed-loop gust response is compared with the experimental open-loop one.In order to guarantee the aircraft model’s safety, the experimental flow velocity is from 12 m/s to 24 m/s, much far away from the flutter velocity of 31 m/s. The frequency of gust generator changes from 1.0 Hz to 6.0 Hz with a fixed increasing step of 0.5 Hz. In the open-loop gust response test, the gust generator upwards the aircraft model is switched on all the time. After the gust frequency keeps still, the flow velocity increases to a given value with a fixed increasing step of 2 m/s. At a fixed velocity and gust frequency, we measured the accelerations and deflections of control surfaces for 15 seconds. A Fourier transformation is conducted on all the 15-second data to obtain the response amplitude in the frequency domain. The amplitude of the open-loop wing tip acceleration at different velocities and frequencies is shown in Figure2 and Table 1.Table 1
The open-loop wing tip acceleration data due to sinusoidal gust in the wind tunnel test.
Frequency: Hz
Acceleration at 12 m/s : g
Acceleration at 14 m/s : g
Acceleration at 16 m/s : g
Acceleration at 18 m/s : g
Acceleration at 20 m/s : g
Acceleration at 22 m/s : g
Acceleration at 24 m/s : g
1.0
0.10
0.11
0.12
0.17
0.17
0.18
0.18
1.5
0.21
0.35
0.29
0.28
0.29
0.40
0.42
2.0
0.70
0.69
0.79
0.70
0.79
0.83
0.86
2.5
0.55
0.47
0.80
1.09
0.89
1.30
0.91
3.0
0.42
0.61
0.85
1.28
1.34
1.56
1.82
3.5
0.35
0.37
0.54
1.26
1.24
1.28
1.69
4.0
0.39
0.46
0.73
1.23
1.40
1.80
2.16
4.5
0.32
0.44
0.61
1.06
0.91
1.28
1.82
5.0
0.34
0.42
0.70
0.86
1.32
1.52
1.93
5.5
0.27
0.34
0.49
1.36
0.99
1.49
2.16
6.0
0.33
0.41
0.53
1.41
1.23
1.97
2.48Figure 2
The open-loop wing tip acceleration due to sinusoidal gust in the wind tunnel test.Figure3 is the closed-loop acceleration by the PI controller in the wind tunnel test at the velocity of 18 m/s. From the comparison, it can be seen that the controller cannot suppress the wing tip acceleration at a low gust frequency of 1 Hz or at a high gust frequency of 5 Hz. It works even worse at frequencies higher than 5 Hz. The controller alleviates the wing tip acceleration well at other gust frequencies.Figure 3
Comparison of the experimental open-loop gust response and the closed-loop one by PI controller.


## 3.2. Gust Response Alleviation at Fixed Flow Velocity and Gust Frequency

Gust response alleviation by GPC is simulated at a fixed flow velocity and fixed gust frequency. First, the open-loop wing tip acceleration at the flow velocity of 24 m/s and gust frequency of 2 Hz is employed to identify an AR model. Then, a GPC controller at this specified test condition is designed to alleviate the wing tip response. As mentioned above, at a fixed flow velocity, though the gust input is not measured, it can be given as a known form. That is,ud1=5sin⁡(4πt). Some parameters for model identification are set as below: the input-output data length for identification is 400. The order of AR model p is 5. The number of control inputs is 2. The number of response is one and the number of gust input is one. The weighting matrices for controller design are set to be unity ones. Considering the feasibility of controllers in the wind tunnel test, the lengths of prediction steps and control command steps should not to be too long. Different prediction steps are attempted to choose a better one. Three cases for prediction steps are considered. Their corresponding alleviation results are shown in Table 2. From this table, more prediction steps result in better alleviation results. However, the lengths of prediction window and controller window should not to be too large. It may make controller complicated in the real wind tunnel test. Hence, in the following controller design process, the length of prediction steps hp is fixed to be 5, and the length of control steps hc is also fixed to be 5.Table 2
Alleviation result for different length of prediction steps.
Case number
Length of prediction and control steps
Open-loop acceleration: g
Closed-loop acceleration: g
Alleviation factor
Case  1
h
p
=
1, hc=1
0.858
0.823
4.1%
Case  2
h
p
=
3, hc=3
0.858
0.371
56.8%
Case  3
h
p
=
5, hc=5
0.858
0.11
87.2%The gust response alleviation result for Case 3 is shown in Figure4. The response data in the first 2 seconds is used for AR model identification. After 2 seconds, the prediction controller is switched on. From the figure, it can be seen that the output of the AR model is nearly the same as the experimental acceleration. It indicates that the identified model can represent the aeroelastic behavior at this test condition, and the AR model has a high accuracy for generalized prediction control. After the prediction controller is switched on, the acceleration response reduces to a relatively small value of 0.11 g in a short time. The response comparison in the frequency domain is shown in Figure 5, and the deflections of control surfaces are shown in Figure 6. The deflection of aileron has a shock on the switch-on time step and then it oscillates between −0.1 degrees and 0.1 degrees in the later time. The experimental deflection of elevator before 2 seconds exceeds 12 degrees. The reason is as follows: the elevator is used for stabilizing the aircraft model in the wind tunnel test, and the model is up-mounted which is unstable at the equilibrium point. Hence, its deflection is a little large before the time of 2 seconds. In the closed-loop simulation, it is only used for active gust response alleviation. Hence, after the time of 2 seconds, its deflection remains oscillating between −0.1 degrees and 0.1 degrees, similar to aileron.Figure 4
The comparison of open-loop and closed-loop accelerations in the time domain forhp=5, hc=5.Figure 5
The comparison of open-loop and closed-loop accelerations in the frequency domain.Figure 6
Deflections of control surfaces.
(a)
Elevator
(b)
AileronThe open-loop responses are measured at each fixed gust frequency and each fixed flow velocity. Hence, the AR model is identified straightforward at each test condition too. The control input in (12) is updated with corresponding test conditions to match the current identified AR model. Therefore, there are 77 AR models and 77 discrete control laws by the standard GPC method. The data length for model identification is set to be 2000 in this case. A Fourier transformation is conducted on the closed-loop response to get the response amplitude. The closed-loop response amplitude at each test condition is shown in Table 3.Table 3
Gust response alleviation result at fixed flow velocity and gust frequency.
Frequency: Hz
Acceleration at 12 m/s : g
Acceleration at 14 m/s : g
Acceleration at 16 m/s : g
Acceleration at 18 m/s : g
Acceleration at 20 m/s : g
Acceleration at 22 m/s : g
Acceleration at 24 m/s : g
1.0
0.10
0.11
0.11
0.15
0.16
0.15
0.12
1.5
0.21
0.23
0.23
0.25
0.20
0.15
0.27
2.0
0.68
0.70
0.65
0.74
∗
0.74
0.42
0.59
2.5
0.38
0.46
0.61
0.50
0.80
0.88
0.63
3.0
0.33
0.61
0.82
1.15
0.82
0.89
0.57
3.5
0.27
0.35
0.53
0.70
0.40
0.45
0.86
4.0
0.38
0.36
0.71
1.26
∗
1.12
1.52
1.09
4.5
0.21
0.30
0.41
0.70
0.82
0.60
0.57
5.0
0.34
0.37
0.70
0.88
∗
0.62
0.37
0.55
5.5
0.18
0.27
0.48
0.40
0.80
0.78
0.30
6.0
0.31
0.37
0.38
0.44
0.51
0.81
0.91When the closed-loop response is larger than the open-loop one, it is denoted as bold star in Table3. From the comparison of Tables 1 and 3, it can be seen that the standard GPC controller suppresses the acceleration response at most test conditions. However, it is not always effective at all test conditions. At some frequencies of 18 m/s, the closed-loop response is even larger than the open one. When the data length for identification is only 400, the alleviation effect is even worse. From the alleviation comparison, the standard GPC controller does not work much better than the PI controller at some low frequencies. Moreover, it needs to identify a new AR model when either the flow velocity or gust frequency changes. This may waste a lot of time to switch the controllers frequently in real wind tunnel test.


## 3.3. Gust Response Alleviation at Varying Frequencies

In the above section, we need to identify 77 AR models from one test condition to another. In order to reduce the number of identification models, we need to modify the standard GPC design process. A straightforward way is to merge input-output data of several test conditions to an entire data series. Then we can identify one unified AR model with these merged data. The difficulty is how to tackle with the varying gust input at different test conditions. From (2), two factors may affect the gust velocity. One is the amplitude; the other is the gust frequency. The gust amplitude may keep constant at a fixed flow velocity. Hence, in this section, we merged the test data of all frequencies at a fixed-flow-velocity condition. Consequently, the gust input series is written as(15)ud=sin⁡ω1t1⋯sin⁡ω1tLsin⁡ω2t1⋯sin⁡ω11t1⋯sin⁡ω11tL.Together with the merged experimental acceleration response, the above gust input is applied to identify one unified AR model at a specified flow velocity. By this approach, there will be 7 AR models for 7 flow velocities from 12 m/s to 24 m/s. At each flow velocity, each AR model is employed to design a GPC control law. The identification and controller parameters are set as follows:L=2000, p=5, hp=5, hc=5, and R and Q are unity matrices. Afterwards, each control law is acted on the aircraft at corresponding test velocity. The results are shown in Table 4.Table 4
Gust response alleviation result for varying frequencies and fixed velocity.
Frequency: Hz
Acceleration at 12 m/s : g
Acceleration at 14 m/s : g
Acceleration at 16 m/s : g
Acceleration at 18 m/s : g
Acceleration at 20 m/s : g
Acceleration at 22 m/s : g
Acceleration at 24 m/s : g
1.0
0.08
0.20
∗
0.06
0.27
∗
0.20
∗
0.17
0.14
1.5
0.06
0.17
0.05
0.22
0.15
0.13
0.12
2.0
0.10
0.32
0.09
0.42
0.26
0.25
0.24
2.5
0.08
0.25
0.07
0.33
0.20
0.20
0.19
3.0
0.14
0.41
0.15
0.61
0.34
0.39
0.37
3.5
0.10
0.21
0.14
0.47
0.26
0.33
0.31
4.0
0.18
0.11
0.31
0.77
0.46
0.74
0.64
4.5
0.13
0.05
0.20
0.41
0.33
0.80
0.61
5.0
0.21
0.21
0.22
0.11
0.49
1.58
∗
0.82
5.5
0.15
0.21
0.12
0.31
0.38
0.89
0.45
6.0
0.24
0.39
0.16
0.77
0.67
1.14
0.70In this case, the deflection of aileron is less than 2 degrees at the whole test condition. The bold star mark denotes that the closed-loop response is larger than that of the open-loop one. From this table, it can be seen that the controller does not alleviate the wing tip response at some low frequencies. It may be caused by the turbulence at low frequencies. From Figure3, in the real wind tunnel test, the PI controller does not alleviate gust response at these frequencies, either. At the low or high frequencies, we have to develop other excellent control methods. Note that the data length for identification may influence the alleviation effect.


## 3.4. Gust Response Alleviation at Varying Velocities

In the above section, we identify a unified AR model at a specified flow velocity. Hence, we need to use 7 different AR models for GPC design on the total 77 test conditions. This will save the switch time in the real test. In this section, we still want to merge the discrete data into one data segment at several test conditions. Other than identifying AR models for varying frequencies, alternatively, we identify the AR model for varying velocities at the same gust frequency. In this case, the amplitude of gust disturbance varies with flow velocities. Pointed out in Section2.2, it is represented by the polynomials of flow velocity, and the gust frequency is fixed for a specific AR model. At a specific gust frequency ω0, the gust input is written as(16)udk=sin⁡ω0tk,Visin⁡ω0tk,Vi2sin⁡ω0tkT,where Vi is the flow velocity; it ranges from 12 m/s to 24 m/s. A significant change on the standard GPC method is to increase the dimension of gust disturbance from one to three. Together with merged experimental response, the specific identified AR model can be adaptive to varying velocities. There are total of 11 different gust frequencies in the test. Hence, we have to identify 11 different AR models. By this approach, at each gust frequency in the wind tunnel test, we use the same AR model for gust alleviation. Therefore, we do not need to change the parameters of T, B, A, Dp, and Df at the same gust frequency. It may save the experimental time of switching controllers at different test conditions. The length of each data segment for identification is chosen as 2000. Simulation results indicate that this data length will result in better alleviation effect. The other parameters are set the same as above. Table 5 shows the gust response alleviation results at each test condition for varying velocities. The bold star infers that the closed loop response is not alleviated by the controller.Table 5
Gust response alleviation result for fixed frequency and varying velocities.
Frequency: Hz
Acceleration at 12 m/s : g
Acceleration at 14 m/s : g
Acceleration at 16 m/s : g
Acceleration at 18 m/s : g
Acceleration at 20 m/s : g
Acceleration at 22 m/s : g
Acceleration at 24 m/s : g
1.0
0.13
∗
0.07
0.03
0.01
0.03
0.07
0.14
1.5
0.13
0.12
0.09
0.06
0.07
0.14
0.24
2.0
0.56
0.16
0.24
0.47
0.62
0.70
0.73
2.5
0.20
0.11
0.16
0.19
0.18
0.13
0.19
3.0
0.41
0.23
0.35
0.45
0.68
1.21
2.02
∗
3.5
0.26
0.05
0.05
0.05
0.14
0.40
0.78
4.0
0.22
0.11
0.27
0.25
0.14
0.40
0.94
4.5
0.28
0.04
0.16
0.21
0.22
0.40
0.77
5.0
0.03
0.19
0.26
0.29
0.40
0.67
1.10
5.5
0.20
0.19
0.32
0.33
0.21
0.20
0.57
6.0
0.30
0.32
0.28
0.19
0.30
0.68
1.22From this table, it can be seen that this controller suppresses most of the response. However, the closed-loop response is larger than the open-loop one at 2 test conditions. The amplified closed-loop response mainly occurs at very low and very high frequencies, similar to the results by the PI controller. At these test conditions, the gust frequency is far away from its frequency of first bending mode, which is about 3 Hz.


## 3.5. Gust Alleviation at Varying Flow Velocities and Gust Frequencies

In this section, we merged the whole input-output data to one data series to identify only one AR model. There are total of seven different flow velocities. And, at each velocity, there are 11 different frequencies. Hence, a total 77 varying open-loop responses are used to identify a unified AR model. The first 0.6 seconds of each data segment is chosen to compose the total 46.2 seconds open-loop response data. A fifth-order AR model with 3 gust inputs is identified by the linear square algorithm. The comparison of experimental open-loop response and identified acceleration is shown in Figure7. Figure 7(a) shows the response comparison for lower frequency at varying flow velocity. Figure 7(b) shows the response comparison for higher frequency. From the comparison, it can be seen that though the order of AR model is only five, the experimental response can also agree well with the identified one under varying flow velocities. The identified model can describe the aircraft model’s aeroelastic behavior well. Noting that there is some disturbance when the gust frequency is only 1 Hz, the response of identified model in the low frequency may have some discrepancies from the experimental one, indicated in Figure 7(a). Noting that alleviation factor may fluctuate with the length of each data segment, 0.6-second data here is enough to get a not bad alleviation result.Figure 7
The comparison between the experimental response and the output of identified AR model.
(a)
Response at the first 5 seconds
(b)
Response at the last 5 secondsAfter a unified AR model is identified using the 77 pieces of test data, the modified GPC controller is applied to gust response alleviation at all the experimental conditions. The amplitude of closed-loop response is shown in Table6. From this table, it can be seen that the modified GPC controller is effective at all the test conditions. Comparing it with the open-loop response shown in Table 1, the alleviation factor can be up to 20% for all the test conditions. The closed-loop response and aileron’s deflection are shown in Figure 8. It is compared with the experimental one at the velocity of 24 m/s and the gust frequency of 6 Hz. Though there is a shock at the time of controller acting on the aircraft model, the deflection does not reach 0.5 degrees all the time. It needs only a little control effort to alleviate the wing tip acceleration.Table 6
Gust response alleviation result for varying frequencies and varying velocities.
Frequency: Hz
Acceleration at 12 m/s : g
Acceleration at 14 m/s : g
Acceleration at 16 m/s : g
Acceleration at 18 m/s : g
Acceleration at 20 m/s : g
Acceleration at 22 m/s : g
Acceleration at 24 m/s : g
1.0
0.03
0.01
0.01
0.02
0.05
0.09
0.14
1.5
0.02
0.01
0.01
0.02
0.03
0.06
0.09
2.0
0.03
0.02
0.02
0.03
0.05
0.09
0.14
2.5
0.02
0.01
0.01
0.02
0.03
0.06
0.09
3.0
0.04
0.03
0.03
0.03
0.05
0.09
0.15
3.5
0.03
0.02
0.02
0.02
0.03
0.07
0.12
4.0
0.05
0.04
0.04
0.03
0.05
0.12
0.23
4.5
0.04
0.03
0.03
0.02
0.04
0.10
0.20
5.0
0.07
0.05
0.06
0.03
0.07
0.20
0.39
5.5
0.05
0.04
0.05
0.02
0.05
0.16
0.31
6.0
0.07
0.06
0.08
0.03
0.08
0.27
0.52Figure 8
The closed-loop response (a) and the deflection of aileron (b).
(a)
(b)


## 4. Conclusions

The gust response alleviation of an aircraft model is developed with the modified GPC design algorithm. The simulated closed-loop response by GPC is compared with experimental results of the sinusoidal gust response test in the wind tunnel. In order to adapt the controller to varying experimental flow velocity, the gust disturbance is represented by second-order polynomials of flow velocities. Consequently, the influence of varying gust input on the AR model is considered. By the wind tunnel test and closed-loop response simulations, some conclusions are made as follows:(1)
The standard GPC can be applied to gust response alleviation in the wind tunnel test. Compared with model-based method, there is no need to construct a theoretical aeroelastic model by the data-based GPC method. It can suppress the wing tip response significantly at high frequencies. However, it does not work well at some low frequencies.(2)
The modified GPC design algorithm can be adapted to varying flow velocities and frequencies. It is advantageous in the wind tunnel test that we do not need to switch the control law in the whole test. It may save the wind tunnel test time. The modified GPC controller is acted on the experimental open-loop acceleration data. Simulation indicates that the unified controller can alleviate the wing tip accelerations up to 20%. Note that the alleviation effect may fluctuate with different parameters chosen in the model identification process.



---

*Source: 348971-2015-05-31.xml*