Active flutter suppression in a 2-D airfoil using linear matrix inequalities techniques

Flutter is an in-flight vibration of flexible structures caused by energy in the airstream absorbed by the lifting surface. This aeroelastic phenomenon is a problem of considerable interest in the aeronautic industry, because flutter is a potentially destructive instability resulting from an interaction between aerodynamic, inertial, and elastic forces. To overcome this effect, it is possible to use passive or active methodologies, but passive control adds mass to the structure and it is, therefore, undesirable. Thus, in this paper, the goal is to use linear matrix inequalities (LMIs) techniques to design an active state-feedback control to suppress flutter. Due to unmeasurable aerodynamic-lag states, one needs to use a dynamic observer. So, LMIs also were applied to design a state-estimator. The simulated model, consists of a classical flat plate in a two-dimensional flow. Two regulators were designed, the first one is a non-robust design for parametric variation and the second one is a robust control design, both designed by using LMIs. The parametric uncertainties are modeled through polytopic uncertainties. The paper concludes with numerical simulations for each controller. The open-loop and closed-loop responses are also compared and the results show the flutter suppression. The perfomance for both controllers are compared and discussed. Copyright © 2006 by ABCM.

Investigation of helicopter rotor blade flutter and flapwise bending response in hovering

"Propulsion Laboratory, Contract no. AF33(616)-5210, Project no. 3137, Task no. 33113."

Bibliographical review of panel flutter and effects of aerodynamic noise

"Materials Laboratory, Contract no. AF 33(616)-5426, Project no. 7360."

The use of Gramian matrices for aeroelastic stability analysis

Most of the established procedures for analysis of aeroelastic flutter in the development of aircraft are based on frequency domain methods. Proposing new methodologies in this field is always a challenge, because the new methods need to be validated by many experimental procedures. With the interest for new flight control systems and nonlinear behavior of aeroelastic structures, other strategies may be necessary to complete the analysis of such systems. If the aeroelastic model can be written in time domain, using state-space formulation, for instance, then many of the tools used in stability analysis of dynamic systems may be used to help providing an insight into the aeroelastic phenomenon. In this respect, this paper presents a discussion on the use of Gramian matrices to determine conditions of aeroelastic flutter. The main goal of this work is to introduce how observability gramian matrix can be used to identify the system instability. To explain the approach, the theory is outlined and simulations are carried out on two benchmark problems. Results are compared with classical methods to validate the approach and a reduction of computational time is obtained for the second example. © 2013 Douglas Domingues Bueno et al.

Asymptotic Description of Flutter Amplitude Saturation by Nonlinear Friction Forces

The computation of the non-linear vibration dynamics of an aerodynamically unstable bladed-disk is a formidable numerical task, even for the simplified case of aerodynamic forces assumed to be linear. The nonlinear friction forces effectively couple dif- ferent travelling waves modes and, in order to properly elucidate the dynamics of the system, large time simulations are typically required to reach a final, saturated state. Despite of all the above complications, the output of the system (in the friction microslip regime) is basically a superposition of the linear aeroelastic un- stable travelling waves, which exhibit a slow time modulation that is much longer than the elastic oscillation period. This slow time modulation is due to both, the small aerodynamic effects and the small nonlinear friction forces, and it is crucial to deter- mine the final amplitude of the flutter vibration. In this presenta- tion we apply asymptotic techniques to obtain a new simplified model that captures the slow time dynamics of the amplitudes of the travelling waves. The resulting asymptotic model is very re- duced and extremely cheap to simulate, and it has the advantage that it gives precise information about the characteristics of the nonlinear friction models that actually play a role in the satura- tion of the vibration amplitude.

Asymptotic Description of Flutter Amplitude Saturation by Nonlinear Friction Forces

The computation of the non-linear vibration dynamics of an aerodynamically unstable bladed-disk is a formidable numerical task, even for the simplified case of aerodynamic forces assumed to be linear. The nonlinear friction forces effectively couple dif- ferent travelling waves modes and, in order to properly elucidate the dynamics of the system, large time simulations are typically required to reach a final, saturated state. Despite of all the above complications, the output of the system (in the friction microslip regime) is basically a superposition of the linear aeroelastic un- stable travelling waves, which exhibit a slow time modulation that is much longer than the elastic oscillation period. This slow time modulation is due to both, the small aerodynamic effects and the small nonlinear friction forces, and it is crucial to deter- mine the final amplitude of the flutter vibration. In this presenta- tion we apply asymptotic techniques to obtain a new simplified model that captures the slow time dynamics of the amplitudes of the travelling waves. The resulting asymptotic model is very re- duced and extremely cheap to simulate, and it has the advantage that it gives precise information about the characteristics of the nonlinear friction models that actually play a role in the satura- tion of the vibration amplitude.

Supersonic flutter of composite skin panels of repeated-sublaminate layup

In high-speed aerospace vehicles, supersonic flutter is a well-known phenomenon of dynamic instability to which external skin panels are prone. In theory, the instability stage is expressed by the 'flutter critical parameter' Q(crit), which is a function of the stiffness-, and dynamic pressure parameters. For a composite skin panel, Q(crit) can be maximised by lay-up optimisation. Repeated-sublaminate lay-up schemes possess good potential for economical lay-up optimisation because the corresponding effort is limited to a family of sublaminates of few layers only. When Q(crit) is obtained for all sublaminates of a family, and the sublaminates ranked accordingly, the resulting ranking reveals not only the optimum lay-up, but also the near-optimum lay-ups, which are useful design alternatives, and the inferior lay-ups which should be avoided. In this paper, we examine sublaminate-ranking characteristics for a composite panel prone to supersonic flutter. In particular, we consider a simple supported midplane-symmetrical rectangular panel of typical aspect ratio alpha and flow angle psi, and for four-layered sublaminates, obtain the Q(crit)-based rankings for a wide range of the number of repeats, r. From the rankings, we find that an optimum lay-up can exist for which the outermost layer is oriented wide of, rather than along, the flow. Furthermore, for many lay-ups other than the optimum and the inferior, we see that as r increases, Q(crit) undergoes significant change in the course of converging. To reconcile these findings, eigenvalue-coalescence characteristics are discussed in detail for specific cases.

KFVS (Kinetic Flux Vector Splitting) on moving grids for unsteady aerodynamics

Accurate numerical solutions to the problems in fluid-structure (aeroelasticity) interaction are becoming increasingly important in recent years. The methods based on FCD (Fixed Computational Domain) and ALE (Alternate Lagrangian Eulerian) to solve such problems suffer from numerical instability and loss of accuracy. They are not general and can not be extended to the flowsolvers on unstructured meshes. Also, global upwind schemes can not be used in ALE formulation thus leads to the development of flow solvers on moving grids. The KFVS method has been shown to be easily amenable on moving grids required in unsteady aerodynamics. The ability of KFMG (Kinetic Flux vector splitting on Moving Grid) Euler solver in capturing shocks, expansion waves with small and very large pressure ratios and contact discontinuities has been demonstrated.

Aerodynamics of delta wings with flaps at hypersonic speeds

Control surface effectiveness is an important parameter for any aeroplane. For a hypersonic aircraft, though the power required to operate the flaps is determined by low speed flying conditions, it is imperative to know the effect of flaps at hypersonic speeds. Hence, studies have been done on this topic by aerodynamicists for over 40 years. In spite of this, only a limited data is available in the literature on this subject. This paper discusses the experimental study of the effect of sweep on the aerodynamic characteristics of thin slab delta wings with flaps at hypersonic speeds. For the purpose of this investigation, a novel special thin six-component balance, which has a thickness of 4mm and can be housed inside wings with 8mm thickness, has been designed. The wings had a sweep of 76degrees, 70degrees and 65degrees, t/c of 0.053 and flaps with 12% of wing area and 12% of wing chord. Testing were done at Mach 8.2, Re number of 2.13 x 10(6) (based on chord), from alpha = -12degrees to 12degrees and flap angle of 20degrees, 30degrees and 40degrees. Separation lengths, measured from Schlieren pictures, clearly show that there is 'no appreciable' effect of sweep on them. Also, using a simple local flow field calculation, the separation has been identified to be transitional in nature. These features of separation reflect in the force data. Because of the small separation length, the flaps (inspite of their small size) were very effective in generating additional C-N, C-M and C-l, which increased with increase in flap angle. In general, the C-N, C-M and X-CP were unaffected by sweep for symmetric flap deflection at positive incidences and asymmetric flap case, For symmetric flap case at negative incidences, only C-N was not influenced by the sweep but C-M decreased and X-CP moved upstream as the sweep is decreased, The wing with lower sweep produces higher CA and lower (L/D)(max) for both symmetric and asymmetric flaps. The rolling moment and adverse yaw increased with decrease in sweep for asymmetric flap deflection. Newtonian theory is shown to be incapable of predicting the effect of sweep on C-l, C-n and on the incremental values of C-N, C-M and C-A. In conclusion, it can be said that a small flap is generally adequate for hypersonic aeroplanes provided they operate at altitudes where transitional and turbulent separation can be expected to occur. This would make the flaps effective and thus enable ample control authority.

Stall flutter of NACA 0012 airfoil at low Reynolds numbers

In the present work, we experimentally study and demarcate the stall flutter boundaries of a NACA 0012 airfoil at low Reynolds numbers (Re similar to 10(4)) by measuring the forces and flow fields around the airfoil when it is forced to oscillate. The airfoil is placed at large mean angle of attack (alpha(m)), and is forced to undergo small amplitude pitch oscillations, the amplitude (Delta alpha) and frequency (f) of which are systematically varied. The unsteady loads on the oscillating airfoil are directly measured, and are used to calculate the energy transfer to the airfoil from the flow. These measurements indicate that for large mean angles of attack of the airfoil (alpha(m)), there is positive energy transfer to the airfoil over a range of reduced frequencies (k=pi fc/U), indicating that there is a possibility of airfoil excitation or stall flutter even at these low Re (c=chord length). Outside this range of reduced frequencies, the energy transfer is negative and under these conditions the oscillations would be damped. Particle Image Velocimetry (PIV) measurements of the flow around the oscillating airfoil show that the shear layer separates from the leading edge and forms a leading edge vortex, although it is not very clear and distinct due to the low oscillation amplitudes. On the other hand, the shear layer formed after separation is found to clearly move periodically away from the airfoil suction surface and towards it with a phase lag to the airfoil oscillations. The phase of the shear layer motion with respect to the airfoil motions shows a clear difference between the exciting and the damping case.

Shock-boundary layer interaction and transonic flutter

The transonic flutter dip of an aeroelastic system is primarily caused by compressibility of the flowing fluid. Viscous effects are not dominant in the pre-transonic dip region. In fact, an Euler solver can predict this flutter boundary with considerable accuracy. However with an increase in Mach number the shock moves towards the trailing edge causing shock induced separation. This shock-boundary layer interaction changes the flutter boundary in the transonic and post-transonic dip region significantly. We discuss the effect of viscosity in changing the flutter boundary in the post-transonic dip region using a RANS solver coupled to a two-degree of freedom model of the structural dynamics of a wing.

Transonic flutter, shocks, and parameter effects

There is a drop in the flutter boundary of an aeroelastic system placed in a transonic flow due to compressibility effects and is known as the transonic dip. Viscous effects can shift the lo-cation of the shock and depending on the shock strength the boundary layer may separate leading to changes in the flutter speed. An unsteady Euler flow solver coupled with the structural dynamic equations is used to understand the effect of shock on the transonic dip. The effect of various system parameters such as mass ratio, location of the center of mass, position of the elastic axis, ratio of uncoupled natural frequencies in heave and pitch are also studied. Steady turbulent flow results are presented to demonstrate the effect of viscosity on the location and strength of the shock.