考虑几何非线性的复合材料机翼气动弹性分析

本文研究结构几何非线性与气动力非平面效应对大展弦比复合材料机翼的气动弹性行为的影响．将非线性有限元法与曲面涡格法结合,计算机翼静气动弹性变形;通过曲面偶极子格网法结合静气动弹性平衡位置处的结构切线刚度,建立气动弹性方程并求解得到机翼颤振速度．针对板模型机翼,分析了迎角对机翼几何非线性气动弹性特性的影响．结果表明：本文复合材料板模型机翼的颤振形式不受水平弯曲模态影响,属于经典弯扭颤振;在几何非线性的影响下,机翼扭转频率随结构变形增大而明显减小,颤振速度随迎角增大而减小．     

Nonlinear dynamics of an aeroelastic airfoil with free-play in transonic flow

Nonlinear dynamic behaviors of an aeroelastic airfoil with free-play in transonic air flow are studied. The aeroelastic response is obtained by using time-marching approach with computational fluid dynamics (CFD) and reduced order model (ROM) techniques. Several standardized tests of transonic flutter are presented to validate numerical approaches. It is found that in time-marching approach with CFD technique, the time-step size has a significant effect on the calculated aeroelastic response, especially for cases considering both structural and aerodynamic nonlinearities. The nonlinear dynamic behavior for the present model in transonic air flow is greatly different from that in subsonic regime where only simple harmonic oscillations are observed. Major features of the responses in transonic air flow at different flow speeds can be summarized as follows. The aeroelastic responses with the amplitude near the free-play are dominated by single degree of freedom flutter mechanism, and snap-though phenomenon can be observed when the air speed is low. The bifurcation diagram can be captured by using ROM technique, and it is observed that the route to chaos for the present model is via period-doubling, which is essentially caused by the free-play nonlinearity. When the flow speed approaches the linear flutter speed, the aeroelastic system vibrates with large amplitude, which is dominated by the aerodynamic nonlinearity. Effects of boundary layer and airfoil profile on the nonlinear responses of the aeroelastic system are also discussed.     

Nonlinear aeroservoelastic analysis of a controlled multiple-actuated-wing model with free-play

In this paper, the effects of structural nonlinearity due to free-play in both leading-edge and trailing-edge outboard control surfaces on the linear flutter control system are analyzed for an aeroelastic model of three-dimensional multiple-actuated-wing. The free-play nonlinearities in the control surfaces are modeled theoretically by using the fictitious mass approach. The nonlinear aeroelastic equations of the presented model can be divided into nine sub-linear modal-based aeroelastic equations according to the different combinations of deflections of the leading-edge and trailing-edge outboard control surfaces. The nonlinear aeroelastic responses can be computed based on these sub-linear aeroelastic systems. To demonstrate the effects of nonlinearity on the linear flutter control system, a single-input and single-output controller and a multi-input and multi-output controller are designed based on the unconstrained optimization techniques. The numerical results indicate that the free-play nonlinearity can lead to either limit cycle oscillations or divergent motions when the linear control system is implemented.     

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.     

A theoretical and experimental investigation of the effects of a steady angle of attack on the nonlinear flutter of a delta wing plate model

Limit cycle oscillations (LCO) of wings on certain modern high performance aircraft have been observed in flight and in wind tunnel experiments. Whether the physical mechanism that gives rise to this behavior is a fluid or structural nonlinearity or both is still uncertain. It has been shown that an aeroelastic theoretical model with only a structural nonlinearity can predict accurately the limit cycle behavior at low subsonic flow for a plate-like wing at zero angle of attack. Changes in the limit cycle and flutter behavior as the angle of attack is varied have also been observed in flight. It has been suggested that this sensitivity to angle of attack is due to a fluid nonlinearity. In this investigation, we study the flutter and limit cycle behavior of a wing in low subsonic flow at small steady angles of attack. Experimental results are compared to those predicted using an aeroelastic theoretical model with only a structural nonlinearity. Results from both experiment and theory show a change in flutter speed as the steady angle of attack is varied. Also the LCO magnitude increased at a given velocity as the angle of attack was increased for both the experiment and theory. While not proving that the observed sensitivity to angle of attack of LCO in aircraft is due to a structural nonlinearity, the results do show that a change in the aeroelastic behavior at angles of attack can be caused by a structural nonlinearity as well as a fluid nonlinearity. In this paper, only structural nonlinearities are considered, but an extension to include aerodynamic nonlinearities would be very worthwhile.     

Aeroelastic analysis of an idealized landing gear door in subsonic potential flow

Landing gear doors on aircraft have experienced flutter during preliminary flight testing. While designs vary widely, landing gear doors are typically plate-like structures with a relatively rigid actuator attached to their inside surface. To better understand the aeroelasticity of landing gear doors, this study investigates the aeroelastic stability of an idealized model. The model consists of a hinged plate with an interior constraint approximating the actuator attachment. The plate is subject to uniform flow, and an unsteady vortex lattice model is coupled to the structural model to predict critical flow velocities. The location and footprint area of the internal constraint, along with plate aspect and mass ratios, are varied to investigate a large parameter space. Results reveal that the critical flow speed and instability mechanism are sensitive to the postulated actuator placement. In general, flutter is the dominant mode of instability when the actuator is postulated in the leading quarter of the plate. In other postulated locations, divergence dominates. However, the exact shape and location of the boundary between flutter and divergence is configuration dependent and found to be especially sensitive to changes in aspect ratio.     

Effects of multiple structural nonlinearities on limit cycle oscillation of missile control fin

Aeroelastic analyses are performed for a 2-D typical section model with multiple nonlinearities. The differences between a system with multiple nonlinearities in its pitch and plunge spring and a system with a single nonlinearity in its pitch are thoroughly investigated. The unsteady supersonic aerodynamic forces are calculated by the doublet point method (DPM). The iterative V-g method is used for a multiple-nonlinear aeroelastic analysis in the frequency domain and the freeplay nonlinearity is linearized using a describing function method. In the time domain, the DPM unsteady aerodynamic forces, which are based on a function of the reduced frequency, are approximated by the minimum state approximation method. Consequently, multiple structural nonlinearities in the 2-D typical wing section model are influenced by the pitch to plunge frequency ratio. This result is important in that it demonstrates that the flutter speed is closely connected with the frequency ratio, considering that both pitch and plunge nonlinearities result in a higher flutter speed boundary than a conventional aeroelastic system with only one pitch nonlinearity. Furthermore, the gap size of the freeplay affects the amplitude of the limit cycle oscillation (LCO) to gap size ratio.     

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.     

CHAOTIC MOTIONS AND LIMIT CYCLE FLUTTER OF TWO-DIMENSIONAL WING IN SUPERSONIC FLOW

Based on the piston theory of supersonic flow and the energy method, the flutter motion equations of a two-dimensional wing with cubic stiffness in the pitching direction are established. The aeroelastic system contains both structural and aerodynamic nonlinearities. Hopf bifurcation theory is used to analyze the flutter speed of the system. The effects of system parameters on the flutter speed are studied. The 4th order Runge-Kutta method is used to calculate the stable limit cycle responses and chaotic motions of the aeroelastic system. Results show that the number and the stability of equilibrium points of the system vary with the increase of flow speed. Besides the simple limit cycle response of period 1, there are also period-doubling responses and chaotic motions in the flutter system. The route leading to chaos in the aeroelastic model used here is the period-doubling bifurcation. The chaotic motions in the system occur only when the flow speed is higher than the linear divergent speed and the initial condition is very small. Moreover, the flow speed regions in which the system behaves chaos axe very narrow.     

Post-flutter response of a flexible cantilever plate in low subsonic flows

Numerical simulations on the post-flutter response of a flexible cantilever plate are carried out by establishing a nonlinear aeroelastic model. The present study shows that chaotic movements may exist in the three-dimensional panel flutter problems in case of low subsonic flows. In the analysis, time traces, phase-plane plots, Poincare maps as well as power spectral densities are employed to identify the dynamic behavior of the system. It is observed that the plate undergoes period-1, period-3 and non-periodic motions with the increase of inflow velocity. The post-flutter behavior is dominated by both geometric and aerodynamic nonlinearities. Numerical results show that wingtip vortexes are in fact an important source of aerodynamic nonlinearities, which have not been fully studied before. The study also provides a criterion on how to choose a coupling strategy in the nonlinear aeroelastic simulation of a low-aspect-ratio flexible structure in low subsonic flows when the dominant nonlinear effect is different in the post-flutter response.     

Unsteady aerodynamic reduced-order modeling of an aeroelastic wing using arbitrary mode shapes

Computational fluid dynamics (CFD) based unsteady aerodynamic reduced-order model (ROM) can offer significant improvements to the efficiency of transonic aeroelastic analysis. To construct a ROM based on mode shapes, one run of CFD solver is needed to compute aerodynamic responses corresponding to mode excitations. When mode shapes change with structure, another run of the CFD solver is required to construct the new ROM. The typically large computational cost associated with repeated runs of the CFD solver impedes the application of existing unsteady aerodynamic reduced-order modeling methods to transonic aeroelastic design optimization and aeroelastic uncertainty analysis. This paper demonstrates a method that can replace the CFD solver used in the process of existing unsteady aerodynamic reduced-order modeling. It can produce aerodynamic responses corresponding to mode excitations for arbitrary mode shapes within a few seconds. Computational cost can be reduced by two orders of magnitude using the mode excitations and the corresponding aerodynamic responses computed by the method to construct the ROMs used for flutter analyses in aeroelastic design optimization or aeroelastic uncertainty analysis in transonic regime compared with the existing unsteady aerodynamic reduced-order modeling methods. Results show that the method can accurately produce the aerodynamic responses corresponding to the mode excitations and predict the flutter characteristics of AGARD 445.6 wings root-attached in three different ways.     

Aeroelastic analysis and active flutter control of nonlinear lattice sandwich beams

Nonlinear aeroelastic characteristics of sandwich beams with pyramidal lattice core are investigated, and the active flutter control of the nonlinear structural system is also studied using the piezoelectric actuator/sensor pair. In the structural modeling, Reddy’s third-order shear deformation theory is applied. Aerodynamic pressure is evaluated by the supersonic piston theory. Hamilton’s principle and the assumed mode method are used to derive the equation of motion. The proportional feedback and the optimal H _∞ control methods are performed to design the controller. In the robust control, the uncertainty caused by omitting the nonlinear terms of the control equation is taken into account, and the mixed sensitivity method is used to solve the problem. The nonlinear aeroelastic property of the sandwich beam is analyzed and is compared with that of the equivalent isotropic beam with the same weight to show the superior aeroelastic characteristics of the lattice sandwich beam. Controlled vibration responses under the two different controllers are calculated and compared. Simulation results show that the robust controller is much more effective than the proportional feedback controller in the flutter suppression of the nonlinear sandwich beam.     

Transonic limit cycle oscillation behavior of an aeroelastic airfoil with free-play

The limit cycle oscillation (LCO) behaviors of an aeroelastic airfoil with free-play for different Mach numbers are studied. Euler equations are adopted to obtain the unsteady aerodynamic forces. Aerodynamic and structural describing functions are employed to deal with aerodynamic and structural nonlinearities, respectively. Then the flutter speed and flutter frequency are obtained by V-g method. The LCO solutions for the aeroelastic airfoil obtained by using dynamically linear aerodynamics agree well with those obtained directly by using nonlinear aerodynamics. Subsequently, the dynamically linear aerodynamics is assumed, and results show that the LCOs behave variously in different Mach number ranges. A subcritical bifurcation, consisting of both stable and unstable branches, is firstly observed in subsonic and high subsonic regime. Then in a narrow Mach number range, the unstable LCOs with small amplitudes turn to be stable ones dominated by the single degree of freedom flutter. Meanwhile, these LCOs can persist down to very low flutter speeds. When the Mach number is increased further, the stable branch turns back to be unstable. To address the reason of the stability variation for different Mach numbers at small amplitude LCOs, we find that the Mach number freeze phenomenon provides a physics-based explanation and the phase reversal of the aerodynamic forces will trigger the single degree of freedom flutter in the narrow Mach number range between the low and high Mach numbers of the chimney region. The high Mach number can be predicted by the freeze Mach number, and the low one can be estimated by the Mach number at which the aerodynamic center of the airfoil lies near its elastic axis. Influence of angle of attack and viscous effects on the LCO behavior is also discussed.     

Passive aeroelastic control of a suspension bridge during erection

This study presents a system based on passively controlled leading- and trailing-edge flaps that is designed to suppress wind-induced instabilities such as flutter and torsional divergence. The utility of the approach is demonstrated on a three-dimensional bridge model. Particular emphasis is placed on the early stages of the deck erection process when the bridge is particularly vulnerable to flutter. The flaps are activated by the deck's movements though passive phase-compensating mechanisms comprising of springs, dampers and inerters. It is demonstrated that optimised compensator parameters, and optimum hinge locations, result in a substantially improved deck aerodynamic performance. Particular importance is given to ensuring that the controlled system has good closed-loop ‘robustness’ properties, or in other words, that the controlled system has a high tolerance to parameter variations and uncertainties in the system dynamics. The practical use of a nonlinear optimisation algorithm with a FE bridge aeroelastic model, which includes the flap dynamics, necessitates the use of reduced-order models. A novel model reduction procedure that is based on the retention of dominant poles is introduced into the aeroelastic modelling framework. Multimodal interactions are observed at the various erection stages and conclusions are drawn with regard to the contributions of various modes of vibration to aeroelastic instabilities. The main advantage of this approach lies on the passive system's simplicity and its ability to simultaneously increase the flutter and torsional divergence boundaries. The Humber Bridge in the U.K. is chosen as a study example for numerical simulations.     

基于流形切空间插值的折叠翼参数化气动弹性建模

变体飞行器的气动弹性力学建模是当前先进飞行器设计的研究热点和难点.然而传统的气动弹性动力学建模方法对于具有结构参变特性的变体飞行器气动弹性力学研究存在建模效率低、计算复杂等问题.本研究提出了一种基于流形切空间插值的可折叠式变体机翼参数化气动弹性建模方法.首先,该方法建立若干个典型折叠角下的折叠翼结构有限元模型,通过流形...     

Nonlinear Analysis of Store-Induced Limit Cycle Oscillations

Flight tests of modern high-performance fighter aircraft reveal the presence of limit cycle oscillation (LCO) responses for aircraft with certain external store configurations. Conventional linear aeroelastic analysis predicts flutter for conditions well beyond the operational envelope, yet these store-induced LCO responses occur at flight conditions within the flight envelope. Several nonlinear sources may be present, including aerodynamic effects such as flow separation and shock-boundary layer interaction and structural effects such as stiffening, damping, and system kinematics. No complete theory has been forwarded to accurately explain the mechanisms responsible. This research examines a two degree-of-freedom aeroelastic system which possesses kinematic nonlinearities and a strong nonlinearity in pitch stiffness. Nonlinear analysis techniques are used to gain insight into the characteristics of the behavior of the system. Numerical simulation is used to verify and validate the analysis. It is found that when system damping is low, the system clearly exhibits nonlinear interaction between aeroelastic modes. It is also shown that although certain applied forcing conditions may appear negligible, these same forces produce large amplitude LCOs under specific realizable circumstances.     

Aeroelastic analysis of a rotating wind turbine blade using a geometrically exact formulation

In this paper, an aeroelastic analysis of a rotating wind turbine blade is performed by considering the effects of geometrical nonlinearities associated with large deflection of the blade produced during wind turbine operation. This source of nonlinearity has become more important in the dynamic analysis of flexible blades used in more recent multi-megawatt wind turbines. The structural modeling, involving the coupled edgewise, flapwise and torsional DOFs, has been performed by using a nonlinear geometrically exact beam formulation. The aerodynamic model is presented based on the strip theory, by applying the principles of quasi-steady and unsteady airfoil aerodynamics. Compared to the conventional steady aerodynamic model, the presented model offers a more realistic consideration of fluid–structure interactions. The resulting governing equation, expanded up to the third-order terms, is analyzed by using the reduced-order model (ROM). The ROM is developed by employing the coupled mode shapes of a cantilever blade under free loading condition. The specifications of the 5MW-NREL wind turbine are used in the simulation study. After verifying the ROM results by comparing them with those of the full FEM model, the model is used in additional static, modal and transient dynamics analyses. The results indicate the important effect of geometrical nonlinearity, especially for larger structural deformations. Moreover, nonlinear analyses reveal the important effects of torsion induced by lateral deformations. It is also found that the governing equation is more efficient, and sufficiently accurate, when it is developed by using the second-order kinetic terms, third-order potential terms and the second-order aerodynamic terms together with third-order damping. Finally, the effects of nonlinearities on the flutter characteristics of wind turbine blades are evaluated through frequency and dynamic analyses.     

LCO flutter of cantilevered woven glass/epoxy laminate in subsonic flow

The paper presents a cantilevered composite wing, aeroelastic characteristics of idealized as a composite flat plate laminate. The composite laminate was made from woven glass fibers with epoxy matrix. The elastic and dynamic properties of the laminate were determined experimentally for aeroelastic calculations. Aeroelastic wind tunnel testing of the laminate was performed and the result showed that flutter, a dynamic instability occurred. The cantilevered laminate also displayed limit cycle amplitude, post-flutter oscillation. The experimental flutter velocity and frequency were verified by our computational analysis.     

翼型厚度对气动弹性特性的影响

结合基于$k$-$\omega$的SST两方程湍流模型,求解雷诺平 均Navier-Stokes方程获得定常和非定常气动力,耦合翼型弹性运动方程,在时间 域内模拟了不同厚度对称翼型在不同迎角下的气动弹性动态过程, 并重点研究了较大迎角下的不同厚度翼型流场特征和气动弹性的性质,研究结果表明：在论 文所涉及的参数情况下,对于迎角从零到大迎角范围,翼型颤振临界速度随迎角的变化不是 单调的. 翼型颤振临界速度迅速下降的起始迎角比最大升力系数对应的迎角小很多.