Research Article  Open Access
Xingzhi Xu, Yakui Gao, Weiguo Zhang, "Iterative Learning Control of a Nonlinear Aeroelastic System despite Gust Load", International Journal of Aerospace Engineering, vol. 2015, Article ID 237804, 13 pages, 2015. https://doi.org/10.1155/2015/237804
Iterative Learning Control of a Nonlinear Aeroelastic System despite Gust Load
Abstract
The development of a control strategy appropriate for the suppression of aeroelastic vibration of a twodimensional nonlinear wing section based on iterative learning control (ILC) theory is described. Structural stiffness in pitch degree of freedom is represented by nonlinear polynomials. The uncontrolled aeroelastic model exhibits limit cycle oscillations beyond a critical value of the freestream velocity. Using a single trailingedge control surface as the control input, a ILC law under alignment condition is developed to ensure convergence of state tracking error. A novel Barrier Lyapunov Function (BLF) is incorporated in the proposed Barrier Composite Energy Function (BCEF) approach. Numerical simulation results clearly demonstrate the effectiveness of the control strategy toward suppressing aeroelastic vibration in the presence of parameter uncertainties and triangular, sinusoidal, and graded gust loads.
1. Introduction
Aeroelasticity is the field of study that deals with the interaction of structural, inertial, and aerodynamic loads. When combined, these loads may encounter adverse instabilities, such as flutter. Within classic aeroelastic theories, linear assumptions are made for the aerodynamics and the structures, and the aeroelastic problem is usually reduced to the straightforward solution of linear coupled equations [1]. However, the assumption of a linear and known structure, as well as the often oversimplified aerodynamic model, usually leads to inaccurate results.
In many cases, transonic dip, flow separation, and dynamic stall can introduce aerodynamic phenomena that classical aeroelasticity is unable to handle. One example is the transonic dip that linear aerodynamics cannot detect. In addition, flow separation and shock oscillations are beyond the capability of classic aeroelasticity [2]. Nonlinear aerodynamic effects are more difficult to analyze in a theoretical fashion, because the dynamics of flow is governed by the full potential, Euler, and NavierStokes formulations where analytical solutions are practically nonexistent [2, 3]. By contrast, structural nonlinearities arising from worn hinges of control surfaces, loose control linkages, and material behavior, as well as various other sources, can significantly complicate the wing dynamics. With structural nonlinearities, the aeroelastic system may exhibit various phenomena, including instability, limit cycle oscillations (LCOs), and chaotic vibration [4, 5]. This topic has been extensively analyzed and reviewed [2]. Wind tunnel experiments, primarily for a typical airfoil section and beamlike wings, can validate the results from numerical or theoretical schemes. In particular, Dowell and colleagues at Duke University [6, 7] and Ko and colleagues at Texas A&M University [8] have made several significant contributions. These researchers designed and installed an aeroelastic modeling with control surface freeplay and even created a periodic or a linear frequency sweep gust excitation in some experiments. Through these experimental facilities, a series of theoretical and experimental studies involving flutter and LCOs, as well as gust responses and alleviation, have been completed.
Several studies have focused on developing strategies to suppress flutter by active control. Control strategies to suppress flutter rely on the use of control surfaces, particularly for the twodimensional wing section with structural stiffness nonlinearity. For example, a classical linear fullstate feedback controller was developed, which guaranteed the wing section system with nonlinear stiffness to stabilize in some circumstances [9]. Bhoir and Singh designed a backsteppingbased output feedback nonlinear control strategy for flutter suppression in [10]. Lee and Singh [11] proposed nonlinear controller for the aeroelastic system using adaptive feedback technique. The statedependent Riccati equation method was developed for nonlinear control problems, which was used to design suboptimal control laws of nonlinear aeroelastic systems considering a quasisteady and unsteady aeroelastic model [12, 13].
A considerable number of studies have dealt with the influence of uncertainty on aeroelastic response prediction. Pettit [14] briefly described general sources of uncertainty that complicate airframe design and testing. In particular, parametric uncertainty in nonlinear pitch stiffness has been modeled in the third and fifthorder coefficients. Generally, uncertainties are specified in the cubic coefficient of the torsional spring and in the initial pitch angle of the airfoil. Beran et al. [15] applied standard probability concepts and Monte Carlo simulation to the study of the LCO of an airfoil with a nonlinear pitch spring. When the uncertainty was considered, adaptive controllers based on partial or full feedback linearization were derived, which were effective for flutter suppression [16]. Experimental results obtained using the adaptive control system were presented in [17], which verified the validity of the proposed method in [16]. Also, for a class of uncertain nonlinear multivariable systems, a higherorder slidingmode control law for the finitetime control has been developed [18]. Motivated by the limited effectiveness of using a single control surface, improvements in the performance of the adaptive controller were investigated through multiple control surfaces [19]. An output feedback control law has been implemented to suppress flutter and adaptively compensate for uncertainties in all the aeroelastic model parameters [20]. In most of the available literature, only parametric uncertainty in pitch stiffness has been considered. The damping uncertainty in the airframe structure and control system is inevitable [21]. Li et al. [22] designed an adaptive control law for flutter suppression of a nonlinear aeroelastic system with damping uncertainty.
Recently, a new iterative learning control (ILC) theory has been developed for the control of uncertain nonlinear systems [23]. This design uses a Barrier Composite Energy Function (BCEF) method with a novel Barrier Lyapunov Function (BLF). The prerequisite for the ILC presented in the design is based on the alignment condition. This is different from Repetitive Control (RC) and merely requires the time to be reset. Nonetheless, in conventional ILC, a constant initial condition must be met in conventional ILC, that is, the time and state must be reset at the beginning of each iteration. This design approach has been applied for robotic manipulators and the other industrial control [24–28]. Several control systems for the prototypical wing section of Block and Strganac [9] have been designed in the past, but the application of ILC theory for this model has not been attempted. As such, it is of interest to develop an ILC flutter controller for the prototypical plungepitch 2D aeroelastic system in the presence of gust loads.
In this paper, the design of a control system for the stabilization of an aeroelastic system in the presence of parameter uncertainties and gust loads is considered. The design is based on ILC control theory. The model has pitch polynomialtype structural nonlinearities and describes the plunge and pitch motion of a wing section with a single control surface. The aeroelastic model has quasisteady linear aerodynamics. A new ILC scheme is derived for a nonlinear system with both parametric and nonparametric uncertainties under alignment condition. A BLF is incorporated in the proposed BCEF approach. Under the proposed control scheme, uniform state tracking error convergence is guaranteed. Simulation results show that the controller suppresses the oscillatory motion of the system, despite uncertainties and triangular, exponential, and sinusoidal gust loads. The main contributions of this work are the following: (1) the authors present rigorous mathematical proof and simulation results utilizing the recently developed iterative learning control (ILC) theory as introduced in [23], (2) BLF is incorporated with BCEF to handle ILC problems under alignment condition, (3) both parametric and nonparametric nonlinear uncertainties can be handled by the proposed control scheme, and (4) compared to the conventional sliding mode control design, the proposed method has good performance in tracking the pitch angle and plunge displacement trajectories.
2. Aeroelastic Model and Control Problem
Figure 1 shows the aeroelastic model. A prototype of this model has been developed in [17]. The secondorder differential equations governing the evolution of the pitch angle (take the clockwise direction as positive) and the plunge displacement (take the downward direction as positive) including the gust load are given bywhere is the mass of the wing section, is the total mass, and is the semichord of the wing. The parameter is the moment of inertia, is the nondimensionalized distance of the center of mass from the elastic axis, and and are the pitch and plunge damping coefficients, respectively. and are the aerodynamic moment and lift, respectively. Assuming a quasisteady aerodynamic model, the aerodynamic lift and moment coefficients are expressed bywhere is the nondimensionalized distance from the midchord to the elastic axis; is the span; and are the lift and moment coefficients per angle of attack, respectively; and are the lift and moment coefficients per control surface deflection (take the clockwise direction as positive), respectively; and and are the effective dynamic and control moment derivatives, respectively.
The aerodynamic force and moment due to the wind gust is modelled as [29]where denotes the disturbance velocity and is a dimensionless time variable defined as .
The nonlinear function is considered a polynomial in given bywhere is constant, .
Equations (1) and (2) are combined to form the following compact statespace form:whereMoreover, .
In this study, with allowance for the uncertainty of stiffness parameters, formula (5) can be expressed as follows:where denotes the part containing known parameters; indicates the matrix related to nonlinear stiffness; and represents the state vector of uncertain parameters, which is expressed asand the system variables are as follows:
Based on the knowledge of the structure of the aeroelastic model, the control objective is to design a control strategy that drives the pitch angle and plunge displacement to zero while adaptively compensating for uncertainties and external disturbances in the system model. In the next section, we show through the theoretical analysis that, under the perturbation of system uncertainties, convergence of state tracking error can be guaranteed. An aeroelastic state space representation, as expressed in (7) above, requires angle against the trailing edge, as control input for suppressing wing section flutter. Let be the controlled state vector. The tracking error for the aeroelastic system can be defined as . In this study, is the desired state vector. Given that the control objective is to suppress the aeroelastic vibrations, will be zero for all time periods.
3. Control Design and Convergence Analysis
The equation of motion can be expressed aswhere denotes the number of iterations; is a vector of system states; is a continuous vector function such that , where is the state tracking error; is known control input distribution matrix; is nonparametric uncertainty; is the parametric uncertainty function, where is a known statedependent matrix function and is a vector of unknown constants; corresponds to the control input signal; and , where is the operation time in each iteration.
Remark 1. For the twodimensional wing section with structural stiffness nonlinearity, parametric uncertainty was modeled in the third and fifthorder stiffness coefficients of the pitch spring. Hence, in this work is parametric uncertainty. In contrast, in [24], is parametric uncertainty.
Assumption 2. The nonparametric uncertainty is locally Lipschitz continuous in ; that is, where and are known bounding functions and is a Euclidean norm for vector .
Assumption 3. There exists a Barrier Lyapunov Function (BLF) and a nonnegative class function , such that, for a vector , as , andThe state tracking error at the th iteration is defined as . Notice that, under alignment condition, , we haveIn order to achieve state tracking error convergence, the control law with ILC scheme in each iteration is proposed aswhere , are positive ILC gains, , , and , represent projection operation such thatwhere is the sign function defined as
Assumption 4. The projection bound satisfies the following condition: , and , ; .
We introduce a nonnegative Barrier Composite Energy Function (BCEF) to facilitate the analysisand is any BLF satisfying Assumption 2.
Remark 5. In this work, we adapt a BLF in [24] aswhich is also positive and will approach to infinity as , where is the predefined bound. From analysis that will be presented later, by BCEF, and keeping it bounded in close loop, we will guarantee convergence of state tracking error.
Theorem 6. For system (10) under alignment condition and control law proposed as (14) and two ILC laws as (15) and (16), the uniform convergence of the state tracking error is guaranteed over , as iteration index .
Proof. Consider the difference of BCEF (19) between two consecutive iterations at time :we haveNotice that, in the case of , we haveWhen , owing to the relationships and , the following can be obtained:Therefore, for the second term of (24), substituting the control law (14) into the integrandwe obtainFor and , applying the property for vectors , , with the same dimension, and , as well as , we havewhere cancels in (30) and (27).
Hence, combining all terms from (27)–(30) yieldsSince , therefore,which indicates that the BCEF defined at is monotonically decreasing over iteration domain.
To prove the finiteness of the time derivative of BCEF for any iteration, for any iteration index , we haveSimilar to (28), for we havewhere we substitute updating law (15) and (16) with .
For and , we havehenceand since and are finite with respect to finite and , (36) indicates that is finite.
For any iterations , we haveNotice that, similar to (28), we can deriveFor and , we haveNotice is a finite term and denote its trace by , and is also a finite term. Hence we getSince in (38)therefore, we can obtainThe initial value of BCEF at th iteration . Since , and can be proved to be bounded by showing that BCEF is bounded in the th iteration. The boundedness of and implies the boundedness of , which guarantees that the boundedness of will be guaranteed; therefore, will be ensured for any time in all iterations.
From (42), we haveSince is positive and is bounded, converges, as ; namely,Therefore, we can conclude that state tracking error converges to zero uniformly; that is,
4. Simulation Results
This section describes the numerical simulations performed to validate the above control law design. The model parameters taken from [17] are collected in Table 1.

The velocity distributions of for simulation are assumed to be (1) a triangular gust, (2) a sinusoidal gust, and (3) a graded gust. For the triangular disturbance input, one haswhere denotes the unit step function, , , and . The sinusoidal and graded velocity distributions arewhere , , and .
In this example, the desired trajectory variable is simply selected as zero. The initial conditions for pitch angle and plunge displacement are set to be and , respectively. All other state variables are also selected as zero. Here the nominal values of the unknown are given as . The upper bounds for the uncertain parameter are arbitrarily selected as . The design parameters are given as follows: , , , , and . The poles of the linearized uncontrolled system without external disturbance for and are (). For and , the poles of the linearized system are (). Therefore, the openloop system is stable for and unstable for . The selected velocity distributions () for sinusoidal gust, triangular gust, and graded gust are shown in Figure 2. The sinusoidal and triangular gusts have finite duration. The openloop responses for the velocity distributions of Figure 2 are shown in Figure 3. For triangular , Figures 3(a) and 3(b) indicate that the pitch angle and plunge displacement converge to zero for . However, for , the system exhibits LCOs after initial transience in all cases (Figures 3(c)–3(f)). Apparently, these undesirable oscillations in the plunge and pitch motion must be suppressed.
(a)
(b)
(c)
(a)
(b)
(c)
(d)
(e)
(f)
First, the responses of the closedloop system, in the presence of a triangular gust load with intensity and are obtained. For a realistic simulation, control surface deflections are limited to . Figure 4 demonstrates the response of the closedloop system to triangular gust disturbance when and . With the addition of control, the pitch angle and plunge displacement converge to zero. The response time for the closedloop is of an order less than 3 s. Compared with the openloop response in Figures 3(a) and 3(b), the closedloop system has a shorter convergent time. Thus, the closedloop system responds quickly after control is implemented. Figure 5 shows the response time histories of the closedloop system with triangular gust of larger intensity and . Evidently, in spite of the stronger gust, the pitch angle and plunge displacement converge to zero in 2 s. This result indicates that the gain of the control input increases with the velocity of the freestream velocity; that is, the efficiency of the control surface increases at this time.
(a)
(b)
(c)
(a)
(b)
(c)
The closedloop system perturbed by a sinusoidal disturbance input (shown in Figure 6) is simulated. The oscillations in the pitch and plunge directions are suppressed within 2 and 3 s, respectively. Figure 7 shows the response time histories of the closedloop system under graded gust of the intensity and . Again the oscillations in the system are suppressed. The convergence time for the pitch angle and plunge trajectories is within 3 s. When the value of increases to 1 while the other parameters of graded gust disturbance remain unchanged, system response remains convergent to a stable state instead of to zero, as shown in Figure 8. The equilibrium values of the stable response convergences in the pitch and plunge directions are approximately 0.065 rad and −0.008 m, respectively.
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
Extensive simulation has been performed for various kinds of wind gust. These results reveal that the control strategy proposed effectively suppresses the oscillation of the plunging displacement and the pitch angle, in spite of the uncertainties in the parameters and wind gusts of different shapes.
Next, for a reference tracking signal, the plunge displacement and pitch angle trajectories are and , respectively. Figures 9 and 10 depicted the corresponding tracking responses. From the results shown in Figure 9, it can be observed that the conventional SMC is much faster than the proposed ILC in tracking the desired pitch trajectory, whereas it is incapable of tracking the plunge displacement trajectory. By contrast, it can be observed from Figures 9 and 10 that the desired tracking objective can be achieved when using the proposed ILC. To further evaluate the tracking performance of the proposed method, the desired tracking signals are set to be and , respectively. In Figures 10 and 11, it is observed that the response curves using conventional SMC have been unable to track the reference signals. However, for more complex reference signals shown in Figure 12, the proposed approach is still effective. From the analysis results, it was also found that the proposed method possesses a remarkably better tracking performance than that using conventional SMC.
5. Conclusions
In this paper, a ILC controller was proposed to suppress aeroelastic vibrations for a nonlinear wing section in the presence of parametric uncertainties and gust load. The control strategy was implemented via a single trailingedge control surface. The proposed scheme was developed based on alignment condition. Both parametric and nonparametric uncertainties can be effectively handled. A Barrier Composite Energy Function (BCEF) method with a novel Barrier Lyapunov Function has been used to guarantee the uniform convergence result for state tracking error. Simulation results for various gust loads were obtained. The controller accomplished suppression of the plunge displacement and pitch angle trajectory oscillations, despite uncertainties in the model parameters and triangular, sinusoidal, and graded gust loads. Unlike the traditional SMC, the proposed method has good tracking performance.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The financial support of this research by National Natural Science Foundation of China, under Grant no. 61374032, is greatly appreciated.
References
 D. Cox and H. C. Curtiss Jr., A Modern Course in Aeroelasticity, Sijthoff and Noordhoff, Alphen aan den Rijn, The Netherlands, 4th edition, 2004.
 B. H. K. Lee, S. J. Price, and Y. S. Wong, “Nonlinear aeroelastic analysis of airfoils: bifurcation and chaos,” Progress in Aerospace Sciences, vol. 35, no. 3, pp. 205–334, 1999. View at: Publisher Site  Google Scholar
 E. Dowell, J. Edwards, and T. Strganac, “Nonlinear aeroelasticity,” Journal of Aircraft, vol. 40, no. 5, pp. 857–874, 2003. View at: Publisher Site  Google Scholar
 Y. Fung, An Introduction to the Theory of Aeroelasticity, Courier Corporation, 1955. View at: MathSciNet
 L. C. Zhao and Z. C. Yang, “Chaotic motions of an airfoil with nonlinear stiffness in incompressible flow,” Journal of Sound and Vibration, vol. 138, no. 2, pp. 245–254, 1990. View at: Publisher Site  Google Scholar  MathSciNet
 D. Tang and E. H. Dowell, “Aeroelastic airfoil with free play at angle of attack with gust excitation,” AIAA Journal, vol. 48, no. 2, pp. 427–442, 2010. View at: Publisher Site  Google Scholar
 M. D. Conner, D. M. Tang, E. H. Dowell, and L. N. Virgin, “Nonlinear behavior of a typical airfoil section with control surface freeplay: a numerical and experimental study,” Journal of Fluids and Structures, vol. 11, no. 1, pp. 89–109, 1997. View at: Publisher Site  Google Scholar
 J. Ko, T. W. Strganac, and A. J. Kurdila, “Stability and control of a structurally nonlinear aeroelastic system,” Journal of Guidance, Control, and Dynamics, vol. 21, no. 5, pp. 718–725, 1998. View at: Publisher Site  Google Scholar
 J. J. Block and T. W. Strganac, “Applied active control for a nonlinear aeroelastic structure,” Journal of Guidance, Control, and Dynamics, vol. 21, no. 6, pp. 838–845, 1998. View at: Publisher Site  Google Scholar
 N. Bhoir and S. N. Singh, “Output feedback nonlinear control of an aeroelastic system with unsteady aerodynamics,” Aerospace Science and Technology, vol. 8, no. 3, pp. 195–205, 2004. View at: Publisher Site  Google Scholar
 K. W. Lee and S. N. Singh, “Control of a wing section using leadingand trailingedge flaps by L1 adaptive feedback despite disturbances,” in Proceedings of the 51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, AIAA Paper 20130333, Dallas, Tex, USA, January 2013. View at: Google Scholar
 S. N. Singh and W. Yim, “State feedback control of an aeroelastic system with structural nonlinearity,” Aerospace Science and Technology, vol. 7, no. 1, pp. 23–31, 2003. View at: Publisher Site  Google Scholar
 D. Li, S. Guo, and J. Xiang, “Aeroelastic dynamic response and control of an airfoil section with control surface nonlinearities,” Journal of Sound and Vibration, vol. 329, no. 22, pp. 4756–4771, 2010. View at: Publisher Site  Google Scholar
 C. L. Pettit, “Uncertainty quantification in aeroelasticity: recent results and research challenges,” Journal of Aircraft, vol. 41, no. 5, pp. 1217–1229, 2004. View at: Publisher Site  Google Scholar
 P. S. Beran, C. L. Pettit, and D. R. Millman, “Uncertainty quantification of limitcycle oscillations,” Journal of Computational Physics, vol. 217, no. 1, pp. 217–247, 2006. View at: Publisher Site  Google Scholar
 J. Ko, T. W. Strganac, and A. J. Kurdila, “Adaptive feedback linearization for the control of a typical wing section with structural nonlinearity,” Nonlinear Dynamics, vol. 18, no. 3, pp. 289–301, 1999. View at: Publisher Site  Google Scholar
 T. W. Strganac, J. Ko, and D. E. Thompson, “Identification and control of limit cycle oscillations in aeroelastic systems,” Journal of Guidance, Control, and Dynamics, vol. 23, no. 6, pp. 1127–1133, 2000. View at: Publisher Site  Google Scholar
 M. Defoort, T. Floquet, A. Kokosy, and W. Perruquetti, “A novel higher order sliding mode control scheme,” Systems & Control Letters, vol. 58, no. 2, pp. 102–108, 2009. View at: Publisher Site  Google Scholar
 S. Gujjula, S. N. Singh, and W. Yim, “Adaptive and neural control of a wing section using leading and trailingedge surfaces,” Aerospace Science and Technology, vol. 9, no. 2, pp. 161–171, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Behal, P. Marzocca, V. M. Rao, and A. Gnann, “Nonlinear adaptive control of an aeroelastic twodimensional lifting surface,” Journal of Guidance, Control, and Dynamics, vol. 29, no. 2, pp. 382–390, 2006. View at: Publisher Site  Google Scholar
 C. C. Marsden and S. J. Price, “Transient and limit cycle simulation of a nonlinear aeroelastic system,” Journal of Aircraft, vol. 44, no. 1, pp. 60–70, 2007. View at: Publisher Site  Google Scholar
 D. Li, J. Xiang, and S. Guo, “Adaptive control of a nonlinear aeroelastic system,” Aerospace Science and Technology, vol. 15, no. 5, pp. 343–352, 2011. View at: Publisher Site  Google Scholar
 J.X. Xu and X. Jin, “Stateconstrained iterative learning control for a class of MIMO systems,” IEEE Transactions on Automatic Control, vol. 58, no. 5, pp. 1322–1327, 2013. View at: Publisher Site  Google Scholar
 P. Sampson, C. Freeman, S. Coote et al., “Using functional electrical stimulation mediated by iterative learning control and robotics to improve arm movement for people with Multiple Sclerosis,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, no. 99, 1 pages, 2015. View at: Publisher Site  Google Scholar
 Y. Q. Chen, K. L. Moore, J. Yu, and T. Zhang, “Iterative learning control and repetitive control in hard disk drive industry—a tutorial,” in Proceedings of the 45th IEEE Conference on Decision and Control (CDC '06), pp. 2338–2351, IEEE, San Diego, Calif, USA, December 2006. View at: Publisher Site  Google Scholar
 Y. M. Zhao, Y. Lin, F. Xi, and S. Guo, “Calibrationbased iterative learning control for path tracking of industrial robots,” IEEE Transactions on Industrial Electronics, vol. 62, no. 5, pp. 2921–2929, 2015. View at: Publisher Site  Google Scholar
 S. Basler, P. Dunow, M. Marquardt, and A. Daasch, “Application of iterative learning control methods for a service robot with multibody kinematics,” in Proceedings of the 20th IEEE International Conference on Methods and Models in Automation and Robotics (MMAR '15), pp. 465–470, Miedzyzdroje, Poland, August 2015. View at: Publisher Site  Google Scholar
 N. R. Kapania and J. C. Gerdes, “Path tracking of highly dynamic autonomous vehicle trajectories via iterative learning control,” in Proceedings of the American Control Conference (ACC '15), pp. 2753–2758, IEEE, Chicago, Ill, USA, July 2015. View at: Publisher Site  Google Scholar
 P. Marzocca, L. Librescu, and G. Chiocchia, “Aeroelastic response of 2D lifting surfaces to gust and arbitrary explosive loading signatures,” International Journal of Impact Engineering, vol. 25, no. 1, pp. 41–65, 2001. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Xingzhi Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.