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.

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.

Experimental evidence of flutter and divergence instabilities induced by dry friction

Flutter and divergence instabilities have been advocated to be possible in elastic structures with Coulomb friction, but no direct experimental evidence has ever been provided. Moreover, the same types of instability can be induced by tangential follower forces, but these are commonly thought to be of extremely difficult, if not impossible, practical realization. Therefore, a clear experimental basis for flutter and divergence induced by friction or follower-loading is still lacking. This is provided for the first time in the present article, showing how a follower force of tangential type can be realized via Coulomb friction and how this, in full agreement with the theory, can induce a blowing-up vibrational motion of increasing amplitude (flutter) or an exponentially growing motion (divergence). In addition, our results show the limits of a treatment based on the linearized equations, so that nonlinearities yield the initial blowing-up vibration of flutter to reach eventually a steady state. The presented results give full evidence to potential problems in the design of mechanical systems subject to friction, open a new perspective in the realization of follower-loading systems and of innovative structures exhibiting 'unusual' dynamical behaviors. © 2011 Elsevier Ltd.

A finite volume unstructured mesh approach to dynamic fluid-structure interaction: an assessment of the challenge of flutter analysis

Computational modelling of dynamic fluid-structure interaction (DFSI) is problematical since conventionally computational fluid dynamics (CFD) is solved using finite volume (FV) methods and computational structural mechanics (CSM) is based entirely on finite element (FE) methods. Hence, progress in modelling the emerging multi-physics problem of dynamic fluid-structure interaction in a consistent manner is frustrated and significant problems in computation convergence may be encountered in transferring and filtering data from one mesh and solution procedure to another, unless the fluid-structure coupling is either one way, very weak or both. This paper sets out the solution procedure for modelling the multi-physics dynamic fluid-structure interaction problem within a single software framework PHYSICA, using finite volume, unstructured mesh (FV-UM) procedures and will focus upon some of the problems and issues that have to be resolved for time accurate closely coupled dynamic fluid-structure flutter analysis.

A finite volume unstructured mesh approach to dynamic fluid–structure interaction: an assessment of the challenge of predicting the onset of flutter

Computational modelling of dynamic fluid–structure interaction (DFSI) is a considerable challenge. Our approach to this class of problems involves the use of a single software framework for all the phenomena involved, employing finite volume methods on unstructured meshes in three dimensions. This method enables time and space accurate calculations in a consistent manner. One key application of DFSI simulation is the analysis of the onset of flutter in aircraft wings, where the work of Yates et al. [Measured and Calculated Subsonic and Transonic Flutter Characteristics of a 45° degree Sweptback Wing Planform in Air and Freon-12 in the Langley Transonic Dynamic Tunnel. NASA Technical Note D-1616, 1963] on the AGARD 445.6 wing planform still provides the most comprehensive benchmark data available. This paper presents the results of a significant effort to model the onset of flutter for the AGARD 445.6 wing planform geometry. A series of key issues needs to be addressed for this computational approach. • The advantage of using a single mesh, in order to eliminate numerical problems when applying boundary conditions at the fluid-structure interface, is counteracted by the challenge of generating a suitably high quality mesh in both the fluid and structural domains. • The computational effort for this DFSI procedure, in terms of run time and memory requirements, is very significant. Practical simulations require even finer meshes and shorter time steps, requiring parallel implementation for operation on large, high performance parallel systems. • The consistency and completeness of the AGARD data in the public domain is inadequate for use in the validation of DFSI codes when predicting the onset of flutter.

Reducing Uncertainty in Aeroelastic Flutter Boundaries Using Experimental Data

Flutter prediction as currently practiced is usually deterministic, with a single structural model used to represent an aircraft. By using interval analysis to take into account structural variability, recent work has demonstrated that small changes in the structure can lead to very large changes in the altitude at which
utter occurs (Marques, Badcock, et al., J. Aircraft, 2010). In this follow-up work we examine the same phenomenon using probabilistic collocation (PC), an uncertainty quantification technique which can eficiently propagate multivariate stochastic input through a simulation code,
in this case an eigenvalue-based fluid-structure stability code. The resulting analysis predicts the consequences of an uncertain structure on incidence of
utter in probabilistic terms { information that could be useful in planning
flight-tests and assessing the risk of structural failure. The uncertainty in
utter altitude is confirmed to be substantial. Assuming that the structural uncertainty represents a epistemic uncertainty regarding the
structure, it may be reduced with the availability of additional information { for example aeroelastic response data from a flight-test. Such data is used to update the structural uncertainty using Bayes' theorem. The consequent
utter uncertainty is significantly reduced across the entire Mach number range.

An investigation of yarn snarling and balloon flutter instabilities in ring spinning

This article reports theoretical and experimental investigation on yarn snarling and balloon fluttering in ring spinning. Yarn snarling and balloon fluttering affect yarn breakage in ring spinning. The theoretical model has incorporated the tangential component of air drag on a ballooning yarn, which was ignored in previous models. The results show that yarn snarling happens in the balloon when the ratio of yarn length in the balloon to balloon height is greater than a specific value that depends on the yarn type and count. Yarn tension experiences an obvious change before and after yarn snarling. The balloon flutter appears between normal balloons while the balloon loops are changing. Fluttering balloon shapes that oscillate periodically between two and three loop configurations as yarn tension varies periodically have also been observed experimentally.

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.