Quasi-steady prediction of coupled bending–torsion flutter under rotating stall

A method is presented in this paper to predict cascade flutter under subsonic stalled flow condition in a quasi-steady manner. The ability to predict the occurrence of aeroelastic flutter is highly important from the compressor design point of view. In the present work, the well known Moore–Greitzer compression system model is used to evaluate the flow under rotating stall and the linearized aerodynamic theory of Whitehead is used to estimate the blade loading. The cascade stability is then predicted by solving the structural model, which is posed as a complex eigenvalue problem. The possibility of occurrence of flutter in both bending and torsional modes is considered and the latter is found to be the dominant one, under subsonic stalled flow, for a large range of frequency ratios examined. It is also shown that the design of compressor blades at frequency ratios close to unity may result in rapid initiation of torsional flutter in the presence of stalled flow. A frequency ratio of 0.9 is primarily emphasized for most part of the study as many interesting features are revealed and the results are physically interpreted. Roughly a pitchfork pattern of energy distribution appears to occur between bending mode and torsional mode which ensures that only one flutter mode is possible at any instant in time. A bifurcation from bending flutter to torsional flutter is shown to occur during which the frequency of the two vibrating modes appear to coalesce for a very short period of time.     

Response analysis of a pitch–plunge airfoil with structural and aerodynamic nonlinearities subjected to randomly fluctuating flows

This study focuses on numerically investigating the response dynamics of a pitch–plunge airfoil with structural nonlinearity under dynamic stall conditions. The aeroelastic responses are investigated for both deterministic and randomly time varying flow conditions. To that end, a pitch–plunge airfoil under dynamic stall condition is considered and the nonlinear aerodynamic loads are computed using a Leishman–Beddoes formulation. It is shown that the presence of structural nonlinearities can give rise to a variety of dynamical responses in the pre-flutter regime. Next, a response analysis under the presence of a randomly fluctuating wind is carried out. It is demonstrated that the route to flutter occurs via a regime of pre-flutter oscillations called intermittency. Finally, the manifestation of these stochastic responses is characterized by invoking stochastic bifurcation concepts. The route to flutter via intermittency is presented in terms of topological changes occurring in the joint-probability density function of the state variables.     

Numerical investigations on the dynamic behavior of a 2-DOF airfoil in the transitional Re number regime based on fully coupled simulations relying on an eddy-resolving technique

The paper is the numerical counterpart of the experimental investigation on the fluid–structure interaction (FSI) of a wing with two degrees of freedom (DOF), i.e., pitch and heave. Wood et al. (2020) has provided the experimental basis by studying the flutter stability of an elastically mounted straight wing (NACA 0012 airfoil) in a wind tunnel considering the transitional Reynolds number regime. Three different configurations with varying distances between the fixed elastic axis and the variable center of gravity were considered. Additional free-oscillation tests in still air were carried out in order to make the mechanical properties of the setup available for the simulations. The present contribution describes the numerical methodology applied consisting of a partitioned coupled solver combining eddy-resolving large-eddy simulations on the fluid side with a solver for the governing equations of the translation and rotation of the rigid wing. In order to prove the parameters provided by the experiment and to determine the pure material damping coefficients not available from the measurements, simulations of 1-DOF free-oscillation tests in still air are carried out and analyzed. For validation purposes the corresponding 2-DOF free-oscillation tests in still air are assessed and a good agreement with the experimental data is achieved. Finally, the wing exposed to a constant free-stream of varying strength is analyzed leading to the characterization of complex instantaneous FSI phenomena such as limit-cycle oscillations and flutter. Under full utilization of the supplementary measurements the predictions are evaluated in detail. Contrary to the experiments the simulations provide the entire fluid data and unique data for the translatory and rotatory movement allowing to investigate the causes of the observed phenomena. Both limit-cycle oscillations and flutter can be reproduced by the coupled FSI predictions.     

Aeroelastic analysis of swept pre-twisted wings

In this paper, aeroelastic modeling of aircraft wings with variations in sweep angle, taper ratio, and variable pre-twist angle along the span is considered. The wing structure is modeled as a classical beam with torsion and bending flexibility. The governing equations are derived based on Hamilton’s principle. Moreover, Peters’ finite state aerodynamic model which is modified to take into account the effects of the wing finite-span, the wing sweep angle, and the wing pre-twist angle, is used to simulate the aerodynamic loads on the wing. The coupled partially differential equations are discretized to a set of ordinary differential equations using Galerkin’s approach. By solving these equations the aeroelastic instability conditions are derived. The results are compared with some experimental and analytical results of previous published papers and good agreement is attained. Effects of the wing sweep angle, taper ratio, bending to torsional rigidity, and pre-twist angle on the flutter boundary in several cases are studied. Results show that these geometrical and physical parameters have considerable effects on the wing flutter boundary.     

Effect of uncertainty on the bifurcation behavior of pitching airfoil stall flutter

In this paper the effect of system parametric uncertainty on the stall flutter bifurcation behavior of a pitching airfoil is studied. The aerodynamic moment on the two-dimensional rigid airfoil with nonlinear torsional stiffness is computed using the ONERA dynamic stall model. The pitch natural frequency, a cubic structural nonlinearity parameter, and the structural equilibrium angle are assumed to be uncertain. The effect on the amplitude of the response, the bifurcation of the probability distribution, and the flutter boundary is considered. It is demonstrated that the system parametric uncertainty results already in 5% probability of pitching stall flutter at a 12.5% earlier position than the point where a deterministic analysis would predict unstable behavior. Probabilistic collocation is found to be more efficient than the Galerkin polynomial chaos method and Monte Carlo simulation for modeling uncertainty in the post-bifurcation domain.     

Aeroelastic stability of a cantilevered plate in yawed subsonic flow

The aeroelastic stability of cantilevered plates with their clamped edge oriented both parallel and normal to subsonic flow is a classical fluid–structure interaction problem. When the clamped edge is parallel to the flow the system loses stability in a coupled bending and torsion motion known as wing flutter. When the clamped edge is normal to the flow the instability is exclusively bending and is referred to as flapping flag flutter. This paper explores the stability of plates during the transition between these classic aeroelastic configurations. The aeroelastic model couples a classical beam structural model to a three-dimensional vortex lattice aerodynamic model. The aeroelastic stability is evaluated in the frequency domain and the flutter boundary is presented as the plate is rotated from the flapping flag to the wing configuration. The transition between the flag-like and wing-like instability is often abrupt and the yaw angle of the flow for the transition is dependent on the relative spacing of the first torsion and second bending natural frequencies. This paper also includes ground vibration and aeroelastic experiments carried out in the Duke University Wind Tunnel that confirm the theoretical predictions.     

On bimodal flutter behavior of a flexible airfoil

The dynamic aeroelastic behavior of an elastically supported airfoil is studied in order to investigate the possibilities of increasing critical flutter speed by exploiting its chord-wise flexibility. The flexible airfoil concept is implemented using a rigid airfoil-shaped leading edge, and a flexible thin laminated composite plate conformally attached to its trailing edge. The flutter behavior is studied in terms of the number of laminate plies used in the composite plate for a given aeroelastic system configuration. The flutter behavior is predicted by using an eigenfunction expansion approach which is also used to design a laminated plate in order to attain superior flutter characteristics. Such an airfoil is characterized by two types of flutter responses, the classical airfoil flutter and the plate flutter. Analysis shows that a significant increase in the critical flutter speed can be achieved with high plunge and low pitch stiffness in the region where the aeroelastic system exhibits a bimodal flutter behavior, e.g., where the airfoil flutter and the plate flutter occur simultaneously. The predicted flutter behavior of a flexible airfoil is experimentally verified by conducting a series of systematic aeroelastic system configurations wind tunnel flutter campaigns. The experimental investigations provide, for each type of flutter, a measured flutter response, including the one with indicated bimodal behavior.     

Experimental investigations on the dynamic behavior of a 2-DOF airfoil in the transitional Re number regime based on digital-image correlation measurements

The present paper investigates the fluid–structure interaction (FSI) of a wing with two degrees of freedom (DOF), i.e., pitch and heave, in the transitional Reynolds number regime. This 2-DOF setup marks a classic configuration in aeroelasticity to demonstrate flutter stability of wings. In the past, mainly analytic approaches have been developed to investigate this challenging problem under simplifying assumptions such as potential flow. Although the classical theory offers satisfying results for certain cases, modern numerical simulations based on fully coupled approaches, which are more generally applicable and powerful, are still rarely found. Thus, the aim of this paper is to provide appropriate experimental reference data for well-defined configurations under clear operating conditions. In a follow-up contribution these will be used to demonstrate the capability of modern simulation techniques to capture instantaneous physical phenomena such as flutter. The measurements in a wind tunnel are carried out based on digital-image correlation (DIC). The investigated setup consists of a straight wing using a symmetric NACA 0012 airfoil. For the experiments the model is mounted into a frame by means of bending and torsional springs imitating the elastic behavior of the wing. Three different configurations of the wing possessing a fixed elastic axis are considered. For this purpose, the center of gravity is shifted along the chord line of the airfoil influencing the flutter stability of the setup. Still air free-oscillation tests are used to determine characteristic properties of the unloaded system (e.g. mass moment of inertia and damping ratios) for one (pitch or heave) and two degrees (pitch and heave) of freedom. The investigations on the coupled 2-DOF system in the wind tunnel are performed in an overall chord Reynolds number range of $9.66\times 10^3\le \text{Re}\le 8.77\times 10^4$. The effect of the fluid-load induced damping is studied for the three configurations. Furthermore, the cases of limit-cycle oscillation (LCO) as well as diverging flutter motion of the wing are characterized in detail. In addition to the DIC measurements, hot-film measurements of the wake flow for the rigid and the oscillating airfoil are presented in order to distinguish effects originating from the flow and the structure.     

大跨度悬索桥钢箱加劲梁中央开槽的颤振控制机理

针对近流线型的闭口钢箱梁这种大跨悬索桥常用的加劲梁断面形式,采用风洞试验与理论分析相结合的方法研究不同宽度中央开槽的颤振控制作用与颤振控制机理.通过对原型断面和5种槽宽开槽断面(开槽宽度分别为原型断面宽度的20%、40%、60%、80%和100%)的试验研究表明,适当宽度的中央开槽可以较大幅度的提高颤振临界风速,但槽宽并非越宽越好,而是存在一个"最优槽宽".采用二维三自由度耦合颤振分析方法从气动阻尼和颤振形态矢量两方面对中央开槽的颤振控制机理进行理论研究,结果表明,适当宽度的中央开槽可以改善气动阻尼的形成和发展规律,并提高颤振发生时竖弯自由度的参与程度和扭弯耦合程度,从而增强结构的颤振稳定性能.     

Direct flutter and limit cycle computations of highly flexible wings for efficient analysis and optimization

The usefulness of flutter as a design metric is diluted for wings with destabilizing (softening) nonlinearities, as a stable high-amplitude limit cycle (subcritical) may exist for flight speeds well below the flutter point. It is thus desired to design aeroelastic structures such that the post-flutter behavior is as benign (i.e., supercritical) as possible, among the other constraints commonly considered in the optimization process. In order to account for these metrics in an accurate and efficient manner, direct tools are utilized to first locate the Hopf-point (flutter speed), and then to obtain a nonlinear perturbation solution via the method of multiple scales. The latter scheme provides a scalar variable whose sign and magnitude dictate the nature of the limit cycle. The accuracy of these methods is demonstrated with a high-aspect-ratio highly flexible wing, modeled with nonlinear beam finite elements and the ONERA dynamic stall tool. Stiffness and inertial design variables are allowed to vary spatially throughout the wing, in order to conduct gradient-based optimization of the limit cycle under flutter and mass constraints. The resulting wing structure demonstrates strongly supercritical behavior, as well as several design conflicts between linear (flutter) and nonlinear (limit cycles) sensitivities, which are not present in the uniform baseline wing.     

Flutter and resonances of a flag near a free surface

We investigate the effects of a nearby free surface on the stability of a flexible plate in axial flow. Confinement by rigid boundaries is known to affect flag flutter thresholds and fluttering dynamics significantly, and this work considers the effects of a more general confinement involving a deformable free surface. To this end, a local linear stability is proposed for a flag in axial uniform flow and parallel to a free surface, using one-dimensional beam and potential flow models to revisit this classical fluid–structure interaction problem. The physical behaviour of the confining free surface is characterized by the Froude number, corresponding to the ratio of the incoming flow velocity to that of the gravity waves. After presenting the simplified limit of infinite span (i.e. two-dimensional problem), the results are generalized to include finite-span and lateral confinement effects. In both cases, three unstable regimes are identified for varying Froude number. Rigidly-confined flutter is observed for low Froude number, i.e. when the free surface behaves as a rigid wall, and is equivalent to the classical problem of the confined flag. When the flow and wave velocities are comparable, a new instability is observed before the onset of flutter (i.e. at lower reduced flow speed) and results from the resonance of a structural bending wave and one of the fundamental modes of surface gravity waves. Finally, for large Froude number (low effect of gravity), flutter is observed with significant but passive deformation of the free surface in response of the flag’s displacement.     

A check on the energy method of predicting blade transonic stall flutter

An improved structural dynamic model of an oscillating blade in two degrees of freedom is combined with an unsteady aerodynamic model for the transonic flow about a cascade with separation, which results in a coupled system. The system is solved in an iterative way between the two models. As a check on the current energy methods, the stall flutter boundaries for two real rotors are predicted by using the present method and the results are compared with the experiments and those predicted by using an energy method.     

On the static and dynamic stability of thin beam conveying fluid

In this paper, numerical investigation of the statical and dynamical stability of aligned and misaligned viscoelastic cantilevered beam is performed with a terminal nozzle in the presence of gravity in two cases: (1) effect of fluid velocity on the flutter boundary of beam conveying fluid and (2) effect of gravity on the buckling boundary of beam conveying fluid. The beam is assumed to have a large width-to-thickness ratio, so the out-of-plane bending rigidity is far higher than the in-plane bending and torsional rigidities. Gravity vector is considered in the vertical direction. Thus, deflection of the beam because of the gravity effect couples the in-plane bending and torsional equations. The beam is modeled by Euler–Bernoulli beam theory, with the flow-induced inertia, Coriolis and centrifugal forces along the beam considered as a distributed load along the beam. Furthermore, the end nozzle is regarded as a lumped mass and modeled as a follower axial force. The extended Hamilton’s principle and the Galerkin method are utilized to derive the bending–torsional equations of motion. The coupled equations of motion are solved as eigenvalue problems. Also, several cases are examined to study the impact of gravity, beam inclination angle, mass ratio, nozzle aspect ratio, bending-to-bending rigidity ratio and bending-to-torsional rigidity ratio on flutter and buckling margin of the system.

     

Validation of CFD‐CSD coupling interface methodology using commercial codes

Coupling interface between computational fluid dynamics (CFD) and computational structural dynamics (CSD) is required to provide exchange of information for the simulation of fluid–structure interaction (FSI) phenomena. Accuracy and consistency of information exchanged through coupling interface between the independent CFD and CSD solvers plays a central role in the simulation and prediction of FSI phenomenon, like flutter. In this paper validation of an implemented coupling interface methodology is presented for subsonic, transonic and near supersonic mach regime. The test case chosen for this purpose is the flutter of AGARD445.6 standard I‐wing weakened model configuration for subsonic to near transonic flow regime. Gambit^® and Fluent^® are used for CFD grid generation and solution of fluid dynamic equations, respectively. CSD modeling and simulation are provided by numerical time integration of modal dynamic equations derived through the finite element modeling in ANSYS^® environment. Copyright © 2009 John Wiley & Sons, Ltd.     

DEVELOPMENT OF A THREE-DIMENSIONAL VISCOUS AEROELASTIC SOLVER FOR NONLINEAR PANEL FLUTTER

A new three-dimensional (3-D) viscous aeroelastic solver for nonlinear panel flutter is developed in this paper. A well-validated full Navier–Stokes code is coupled with a finite-difference procedure for the von Karman plate equations. A subiteration strategy is employed to eliminate lagging errors between the fluid and structural solvers. This approach eliminates the need for the development of a specialized, tightly coupled algorithm for the fluid/structure interaction problem. The new computational scheme is applied to the solution of inviscid two-dimensional panel flutter problems for subsonic and supersonic Mach numbers. Supersonic results are shown to be consistent with the work of previous researchers. Multiple solutions at subsonic Mach numbers are discussed. Viscous effects are shown to raise the flutter dynamic pressure for the supersonic case. For the subsonic viscous case, a different type of flutter behavior occurs for the downward deflected solution with oscillations occurring about a mean deflected position of the panel. This flutter phenomenon results from a true fluid/structure interaction between the flexible panel and the viscous flow above the surface. Initial computations have also been performed for inviscid, 3-D panel flutter for both supersonic and subsonic Mach numbers.     

Flutter suppression of a plate-like wing via parametric excitation

The present work is motivated by the well known stabilizing effect of parametric excitation of some dynamical systems such as the inverted pendulum. The possibility of suppressing wing flutter via parametric excitation along the plane of highest rigidity in the neighborhood of combination resonance is explored. The nonlinear equations of motion in the presence of incompressible fluid flow are derived using Hamilton's principle and Theodorsen's theory for modeling aerodynamic forces. In the presence of air flow, the bending and torsion modes possess nearly the same frequency. Under parametric excitation and in the absence of air flow, each mode oscillates at its own natural frequency. In the neighborhood of combination resonance, the nonlinear response is determined using the multiple scales method at the critical flutter speed and at slightly higher airflow speed. The domains of attraction and bifurcation diagrams are obtained to reveal the conditions under which the parametric excitation can provide stabilizing effect. The basins of attraction for different values of excitation amplitude reveal the stabilizing effect that takes place above a critical excitation level. Below that level, the response experiences limit cycle oscillations, cascade of period doubling, and chaos. For flow speed slightly higher than the critical flutter speed, the response experiences a train of spikes, known as ‘firing,’ a term that is borrowed from neuroscience, followed by ‘refractory’ or recovery effect, up to an excitation level above which the wing is stabilized. The results of the multiple scales method are verified using numerical simulation of the original nonlinear differential equations.     

Nonlinear aeroelasticity of high-aspect-ratio wings excited by time-dependent thrust

Effects of engine placement on flutter characteristics of a very flexible high-aspect-ratio wing are investigated using the code NATASHA (Nonlinear Aeroelastic Trim And Stability of HALE Aircraft). Gravity for this class of wings plays an important role in flutter characteristics. In the absence of aerodynamic and gravitational forces and without an engine, the kinetic energy of the first two modes are calculated. Maximum and minimum flutter speed locations coincide with the area of minimum and maximum kinetic energy of the second bending and torsion modes. Time-dependent dynamic behavior of a turboshaft engine (JetCat SP5) is simulated with a transient engine model and the nonlinear aeroelastic response of the wing to the engine’s time-dependent thrust and dynamic excitation is presented. Below the flutter speed, at the wing tip and behind the elastic axis, the impulse engine excitation leads to a stable limit cycle oscillation; and for the ramp kind of excitation, beyond the flutter speed, at 75 % span, behind the elastic axis, it produces chaotic oscillation in the wing. Both the excitations above the flutter speed are stabilized, inboard of the wing.     

磁场对不同温度场中输流悬臂碳纳米管动态特性的影响

本文在采用经典欧拉-伯努利梁模型的基础上,引入考虑小尺度效应的非局部弹性理论,着重研究不同温度场中输流悬臂单层碳纳米管系统（SWCNT）在外加纵向磁场作用下的颤振失稳问题。基于哈密顿原理获得了该流固耦合系统的振动控制方程及相应的边界条件,应用微分变换法（DTM法）求解此高阶偏微分方程,通过数值计算研究了不同温度场中施加纵向磁场对系统动力学特性的影响。结果表明：施加纵向磁场在不同温度场中都将增强输流悬臂碳纳米管的动态稳定性。然而,这种增强程度却与温度场的变化量有关,在不同温度变化量下,磁场对系统稳定性的增强程度有一个峰值,这意味着,实际应用中,为了提高这类流固耦合系统的动态稳定性,一味提高纵向磁场强度并不可取。     

Experiments on transonic limit-cycle-flutter of a flexible swept wing

The aeroelastic behavior of wing models is nonlinear particularly in the transonic speed range. The interaction between aerodynamic and structural forces can lead to the occurrence of Limit-Cycle Oscillations (LCOs). If in addition the wing model is flexible and backward swept, the kinematic coupling between bending and torsion makes the situation even more complex.In the research project “Aerostabil” such a wing was investigated, which was equipped with pressure transducers in three sections and accelerometers. The experiments were performed in the adaptive test section of the transonic wind tunnel TWG in Göttingen. Already Dietz et al. (2003) have reported about experimental details and preliminary results. Based on these data Bendiksen (2008) studied numerically LCO-flutter behavior using a very similar, theoretical model (G-wing) and Stickan et al. (2014) used the original data as a LCO flutter test case. The influence of flexibility on the steady aerodynamics of the wing was described in Schewe & Mai (2018). In this paper now the flutter experiments with the same flexible model were analyzed systematically in the transonic range 0.84 <Ma <0.89 and for six angles of attack from 1.46°to 2.7°. Maps of stability, LCO amplitudes and instantaneous pressure distributions are presented. It was found that unstable regions are islands, whose extent depends on the angle of attack. A LCO test case, already treated in the literature is examined in more detail. The analysis of the time functions showed that during LCO-flutter the motion induced aerodynamic sectional lift forces particularly in the outer wing are asymmetric and thus acting as amplitude limiter. The reason for the asymmetry lies in the shock/boundary layer interaction. The test case, containing the stages of built-up and the transition to the limit cycle provides an excellent opportunity for improving our knowledge about LCOs and for code validation purposes.     

Panel flutter prediction in two dimensional flow with enhanced piston theory

Piston theory may be used in the high Mach number supersonic flow region and/or in very high frequency subsonic or supersonic flow. In this flow model, the pressure at a point on the fluid-solid interface only depends on the downwash at the same point. However the classical piston theory may not be sufficient for some phenomena in aeroelasticity and aeroacoustics (far field prediction). Dowell and Bliss have created an extension of piston theory that allows for higher order effects that take into account the effect the distribution of downwash on pressure at any point. For simple harmonic motion, expansions in reduced frequency, inverse reduced frequency and/or inverse (square of) Mach number have all been created; The effects of higher order terms in these several expansion in creating an enhanced piston theory was illustrated for plunge and pitch motion of an airfoil (discrete system) by Ganji and Dowell. In the present paper, flutter prediction for a flexible panel in two –dimensional flow is investigated using enhanced piston theory. The goal of the present paper is to demonstrate that an enhance version of piston theory can analyze single degree of freedom flutter of a panel as compared to the classical piston theory and quasi-steady aerodynamic models which can only treat coupled mode flutter.