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.     

Time-varying nonlinear dynamics of a deploying piezoelectric laminated composite plate under aerodynamic force

Using Reddy’s high-order shear theory for laminated plates and Hamilton’s principle, a nonlinear partial differential equation for the dynamics of a deploying cantilevered piezoelectric laminated composite plate, under the combined action of aerodynamic load and piezoelectric excitation, is introduced. Two-degree of freedom (DOF) nonlinear dynamic models for the time-varying coefficients describing the transverse vibration of the deploying laminate under the combined actions of a first-order aerodynamic force and piezoelectric excitation were obtained by selecting a suitable time-dependent modal function satisfying the displacement boundary conditions and applying second-order discretization using the Galerkin method. Using a numerical method, the time history curves of the deploying laminate were obtained, and its nonlinear dynamic characteristics, including extension speed and different piezoelectric excitations, were studied. The results suggest that the piezoelectric excitation has a clear effect on the change of the nonlinear dynamic characteristics of such piezoelectric laminated composite plates. The nonlinear vibration of the deploying cantilevered laminate can be effectively suppressed by choosing a suitable voltage and polarity.     

复合材料层合板的颤振特性及振动抑制研究

研究了超声速流中压电复合材料层合板的颤振特性及振动抑制方法。采用一阶活塞理论计算了超声速流场中的气动压力,基于经典层合板理论和Hamilton原理推导了压电复合材料层合板的动力学模型,设计了滑模观测器以减少观测溢出,通过Lyapunov方法证明观测器的稳定性,应用观测状态设计了LQR控制器,讨论了几何参数、铺设角度对压电复合材料层合板颤振特性的影响,利用SIMULINK仿真求解了层合板的脉冲响应,验证了控制器的有效性。结果表明,合理规划层合板的几何参数和铺设角度可提高系统颤振稳定性,滑模观测器能够较为准确地追踪原始系统且具有良好的鲁棒性,LQR控制可以在一定范围内消除层合板的颤振点,并且能够有效地控制压电复合材料层合板在颤振边界处的振动,Q矩阵越大,振动控制效果越好,压电层厚度越大,LQR控制效果越好。     

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.     

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

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

Cyber-physical flexible wing for aeroelastic investigations of stall and classical flutter

The fluid–structure interactions of a finite aspect ratio, cantilevered, flexible wing were investigated using a cyber-physical system to virtually augment the torsional dynamics of the wing. Cyber-physical systems (CPS), which have recently been pursued by a small number of research groups, have proven to be a very useful mechanism to interrogate the fluid–structure interaction parameter space. The premise of a CPS is to use dynamic feedback control to make a system behave according to desired equations of motion. Systems are composed of embedded hardware and software coupled with real-time computing to give the user the flexibility to quickly explore a range of structural parameters. With the advancement of modern control theory, robotics and embedded systems, CPSs integrate both simulation and physical properties into a smart structure which can be used to push the boundaries of research investigations.The CPS in this work allows for the investigation of dynamic aeroelastic instabilities of a three-dimensional, flexible, rectangular planform wing. Two dynamic instability regimes are observed: first, stall flutter, in which the torsional or pitch mode is excited through the dynamic stall process, and second, coupled (or classical) flutter, in which the pitch mode couples with the bending mode. By varying the torsional stiffness and therefore the frequency of torsional versus bending oscillations of the wing, both of these regimes can be attained at the same aerodynamic conditions using the CPS.     

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.     

Reduced order nonlinear aeroelasticity of swept composite wings using compressible indicial unsteady aerodynamics

Nonlinear dynamic aeroelasticity of composite wings in compressible flows is investigated. To provide a reasonable model for the problem, the composite wing is modeled as a thin walled beam (TWB) with circumferentially asymmetric stiffness layup configuration. The structural model considers nonlinear strain displacement relations and a number of non-classical effects, such as transverse shear and warping inhibition. Geometrically nonlinear terms of up to third order are retained in the formulation. Unsteady aerodynamic loads are calculated according to a compressible model, described by indicial function approximations in the time domain. The aeroelastic system of equations is augmented by the differential equations governing the aerodynamics lag states to derive the final explicit form of the coupled fluid-structure equations of motion. The final nonlinear governing aeroelastic system of equations is solved using the eigenvectors of the linear structural equations of motion to approximate the spatial variation of the corresponding degrees of freedom in the Ritz solution method. Direct time integrations of the nonlinear equations of motion representing the full aeroelastic system are conducted using the well-known Runge–Kutta method. A comprehensive insight is provided over the effect of parameters such as the lamination fiber angle and the sweep angle on the stability margins and the limit cycle oscillation behavior of the system. Integration of the interpolation method employed for the evaluation of compressible indicial functions at any Mach number in the subsonic compressible range to the derivation process of the third order nonlinear aeroelastic system of equations based on TWB theory is done for the first time. Results show that flutter speeds obtained by the incompressible unsteady aerodynamics are not conservative and as the backward sweep angle of the wing is increased, post-flutter aeroelastic response of the wing becomes more well-behaved.     

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

变体飞行器的气动弹性力学建模是当前先进飞行器设计的研究热点和难点. 然而传统的气动弹性动力学建模方法对于具有结构参变特性的变体飞行器气动弹性力学研究存在建模效率低、计算复杂等问题. 本研究提出了一种基于流形切空间插值的可折叠式变体机翼参数化气动弹性建模方法. 首先, 该方法建立若干个典型折叠角下的折叠翼结构有限元模型, 通过流形切空间插值方法建立折叠翼参数化结构动力学模型. 其次, 采用偶极子网格法得到参数化非定常气动力模型, 进而建立气动和结构相互耦合的折叠翼参数化气动弹性模型. 为了验证该参数化建模方法在折叠翼气动弹性分析中的准确性, 本文以一小展弦比折叠翼为研究对象, 从折叠翼自由振动时的参变模态特性、颤振边界预测两方面进行了算例验证, 并与直接计算方法进行了对比, 进一步验证了参数化气动弹性模型的有效性. 研究结果表明, 该参数化气动弹性模型对上述两类问题的计算精度与直接计算方法一致, 并且有着计算效率更高的优势.      

Effects of actuator nonlinearity on aeroelastic characteristics of a control fin

The influences of actuator nonlinearities on actuator dynamics and the aeroelastic characteristics of a control fin were investigated by using iterative V-g methods in subsonic flows; in addition, the doublet-hybrid method (DHM) was used to calculate unsteady aerodynamic forces. The changes of actuator dynamics induced by nonlinearities, such as backlash or freeplay, and the variations of flutter boundaries due to the changes of actuator dynamics were observed. Results show that the aeroelastic characteristics can be significantly dependent on actuator dynamics. Thus, the actuator nonlinearities may play an important role in the nonlinear aeroelastic characteristics of an aeroelastic system. The present results also indicate that it is necessary to seriously consider the influence of actuator dynamics on the flutter characteristics at the design stage of actuators to prevent aeroelastic instabilities of aircraft or missiles.     

Camber effects in the dynamic aeroelasticity of compliant airfoils

This paper numerically investigates the effect of chordwise flexibility on the dynamic stability of compliant airfoils. A classical two-dimensional aeroelastic model is expanded with an additional degree of freedom to capture time-varying camber deformations, defined by a parabolic bending profile of the mean aerodynamic chord. Aerodynamic forces are obtained from unsteady thin airfoil theory and the corresponding compliant-airfoil inertia and stiffness from finite-element analysis. V–g and state-space stability methods have been implemented in order to compute flutter speeds. The study looks at physical realizations with an increasing number of degrees of freedom, starting with a camber-alone system. It is shown that single camber leads to flutter, which occurs at a constant reduced frequency and is due to the lock in between the shed wake and the camber motion. The different combinations of camber deformations with pitch and plunge motions are also studied, including parametric analyses of their aeroelastic stability characteristics. A number of situations are identified in which the flutter boundary of the compliant airfoil exhibits a significant dip with respect to the rigid airfoil models. These results can be used as a first estimation of the aeroelastic stability boundaries of membrane-wing micro air vehicles.     

CFD／CSD紧耦合及新型动网格方法在气动弹性模拟中的应用

研发出一套基于紧耦合的CFD／CSD耦合方法和程序。非定常流场求解采用混合网格有限体积方法,时间离散采用基于LU-SGS隐式格式的双时间步长法。通过求解雷诺平均Navier—Stokes方程模拟了三维机翼的跨音速气动弹性现象。得到了其颤振边界,与风洞实验结果吻合较好,验证了方法和程序的有效性和实用意义。同时将Delau...     

基于气动力降阶模型的跨音速气动弹性稳定性分析

基于离散型输入输出差分模型,运用非定常CFD方法训练信号,然后运用最小二乘方法进行参数辨识,得到降阶的非定常气动力模型,再将该离散差分模型转换为连续时间域内的状态方程。耦合气动状态方程和结构状态方程,得到耦合系统的气动弹性状态方程。求解不同动压下状态矩阵的特征值,根据根轨迹图分析系统的稳定性特性。分析结果与直接耦合CFD/CSD方法结果相吻合,可以计算跨音速非线性气动弹性问题。其计算效率比直接耦合CFD/CSD方法提高1~2个数量级。针对Isogai wing在跨音速出现的S型颤振边界进行了较为细致的分析,阐述了该现象是由于系统诱发颤振的分支随着速度(来流动压)的提高而发生转移所导致的。     

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

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

两种湍流模型时域颤振计算方法研究

采用时域计算分析方法进行了机翼跨音速颤振特性研究。在结构运动网格的基础上,采用格点格式有限体积方法进行空间离散和双时间全隐式方法进行时间推进求解雷诺平均N-S方程。针对流动粘性分别应用了SST湍流模型和SSG雷诺应力模型,通过对跨音速标模算例AGARD445.6机翼的计算结果与实验值的对比分析,其中应用SST湍流模型得到的颤振速度与实验值最为接近,特别是在跨音速段平均相对误差在3%以内;并且计算结果整体上反映了跨音速颤振"凹坑"物理特性,验证了方法的有效性。     

Implications of cubic physical/aerodynamic non-linearities on the character of the flutter instability boundary

The present paper deals with a study of the benign and catastrophic characters of the flutter instability boundary of 2-D lifting surfaces in a supersonic flow field. The objectives of this work are: (i) to contribute to a better understanding of the implications of aerodynamic and physical non-linearities on the character of the flutter boundary and (ii), to outline the effects exerted in the same respect by some important parameters of the aeroelastic system. With the aim of addressing this problem, the method based on the First Liapunov Quantity is used to study the bifurcational behavior of the aeroelastic system in the vicinity of the flutter boundary. The expected outcomes of this study are: (a) to greatly enhance the scope and reliability of the aeroelastic analysis and design criteria of advanced aircraft and, (b) to provide a theoretical basis for the analysis of more complex non-linear aeroelastic systems.     

Parametric aeroelastic modeling based on component modal synthesis and stability analysis for horizontally folding wing with hinge joints

To investigate the aeroelastic stability of a folding wing effectively, a parametric aeroelastic analysis approach is proposed. First, the fixed interface component modal synthesis is used to derive the structural dynamic equation for a folding wing, in which the elastic connection is considered. The unsteady aerodynamic model is established by the doublet lattice method (DLM), and the aeroelastic model is achieved from integration of the DLM with the component modal analysis. The generalized aerodynamic influence coefficient matrix is established by modes kept and constraint modes of each component. The aeroelastic stability of a folding wing is investigated based on the Gram matrix in control theory. The effectiveness of the proposed method is verified via comparison with traditional flutter eigenvalue analysis for both extended and folded configurations. The proposed method identifies coupled modes and improves computational efficiency when compared to classical aeroelastic stability analysis methods, such as the p–k method.     

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.     

Identification of Nonlinear Aeroelastic Systems Based on the Volterra Theory: Progress and Opportunities

The identification of nonlinear aeroelastic systems based on the Volterra theory of nonlinear systems is presented. Recent applications of the theory to problems in computational and experimental aeroelasticity are reviewed. Computational results include the development of computationally efficient reduced-order models (ROMs) using an Euler/Navier–Stokes flow solver and the analytical derivation of Volterra kernels for a nonlinear aeroelastic system. Experimental results include the identification of aerodynamic impulse responses, the application of higher-order spectra (HOS) to wind-tunnel flutter data, and the identification of nonlinear aeroelastic phenomena from flight flutter test data of the active aeroelastic wing (AAW) aircraft.     

热过屈曲功能梯度壁板的气动弹性颤振

功能梯度材料的宏观物理性能随空间位置连续变化,能充分减少不同组份材料结合部位界面性能的不匹配因素.功能梯度壁板用作高速飞行器的热防护结构,能有效消除气动加热带来的壁板内部热应力集中.本文考虑热过屈曲变形引入的结构几何非线性,分析功能梯度壁板的气动弹性颤振边界.基于幂函数材料分布假设,采用混合定律计算功能梯度材料的等效力学性能.根据一阶剪切变形板理论、冯·卡门应变-位移关系和一阶活塞理论,基于虚功原理建立超声速气流中受热功能梯度壁板的非线性气动弹性有限元方程.采用牛顿-拉弗森迭代法数值求解壁板的热屈曲变形,分析超声速气流对热屈曲变形的影响机理.在壁板热过屈曲的静力平衡位置分析动态稳定性,确定了壁板的颤振边界.研究表明,当陶瓷-金属功能梯度壁板的组份材料沿厚度方向梯度分布时,会破坏结构的对称性导致壁板在面内热应力作用下发生指向金属侧的热屈曲变形.超声速气流中壁板热屈曲变形最大的位置随气流速压增大向下游推移,并伴随屈曲变形量的减小.热过屈曲壁板的几何非线性效应会提高壁板的颤振边界,这种影响在高温、低无量纲速压且壁板发生大挠度热屈曲变形时表现显著.较高无量纲气流速压下由于壁板的热屈曲变形被气动力限定在小挠度范围,几何非线性效应不明显.