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

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

基于当流活塞理论的气动弹性计算方法研究

发展了一种高效、高精度的超音速、高超音速非定常气动力计算方法--基于定常CFD技术的当地流活塞理论.运用当地流活塞理论计算非定常气动力,耦合结构运动方程,实现超音速、高超音速气动弹性的时域模拟.运用这种方法计算了一系列非定常气动力算例和颤振算例,并和原始活塞理论、非定常Euler方程结果作了比较.由于局部地使用活塞理论假设,这种方法大大地克服了原始活塞理论对飞行马赫数、翼型厚度和飞行迎角的限制.与非定常Euler方程方法相比,当地流活塞理论的效率很高.     

超声速复合材料层合壁板结构的颤振特性分析

对超声速复合材料壁板结构的气动弹性颠振特性进行了分析研究.采用Hamilton原理和假设模态法建立结构的运动方程,采用活塞理论模拟超声速非定常气动力,通过求解本征值问题,得到结构的固有频率和阻尼比等物理量.数值计算了结构无量纲固有频率随气动压力的变化曲线,确定颤振临界气动压力(或颤振速度),并计算了结构的受迫振动时间响应历程曲线,分析比较了不同纤维铺设方式和不同铺设角度对超声速复合材料壁板结构气动弹性稳定性的影响.本文研究结果对超声速飞行器壁板结构的气动弹性稳定性分析和设计具有理论参考价值.     

基于当地流活塞理论的气动弹性计算方法研究

发展了一种高效、高精度的超音速、高超音速非定常气动力计算 方法------基于定常CFD技术的当地流活塞理论. 运用当地流活塞理论计算非定常 气动力，耦合结构运动方程，实现超音速、高超音速气动弹性的时域模拟. 运用这 种方法计算了一系列非定常气动力算例和颤振算例，并和原始活塞理论、非定 常Euler方程结果作了比较. 由于局部地使用活塞理论假设，这种方法大大地克服 了原始活塞理论对飞行马赫数、翼型厚度和飞行迎角的 限制. 与非定常Euler方程方法相比，当地流活塞理论的效率很高.     

风洞实验中翼型大迎角气动力分散性的机理分析

基于$k$-$\omega$的SST两方程湍流模型, 求解雷诺平均N-S方程获得定常和非定常气动力, 耦合 翼型弹性振动方程, 在时间域内模拟了不同的翼型非定常流动, 重点研究了大迎角下的分离 流问题, 研究结果表明: 在百万雷诺数条件下, 由于振动引起分离涡的不规律脱落, 可能导 致气动力平均值的变化; 而厚度大于20\%的翼型在一定大迎角范围内, 会出现分离涡流场平 衡态的转化, 从流体力学稳定性的角度, 解释了风洞实验中大迎角气动力数据的分散性, 为 大迎角气动力风洞实验的重复性和数据分散性给出了一种新解释.     

复合材料机翼/外挂系统颤振特性分析

综合考虑复合材料机翼的材料耦合和几何耦合的特点,利用Hamilton原理建立了非定常气动力作用下复合材料机翼/外挂系统的气动弹性运动方程。采用微分求积法进行离散,运用MATLAB语言编程对系统响应进行数值仿真。计算结果表明:用微分求积法计算复合材料机翼的颤振特性是可行的;随材料耦合刚度的增加,机翼的颤振临界速度出现了双峰值;而对于复合材料机翼/外挂系统,颤振临界速度随材料耦合刚度的变化规律较为复杂,与外挂参数的设置相关。     

弯度翼型大迎角时分离涡流场的多平衡态现象

基于kω的SST两方程湍流模型,在时间域求解雷诺平均Navier-Stokes方程,模拟弯度翼型大迎角时的分离流动。通过给翼型施加一定形式的扰动,重点关注了翼型弯度对大迎角分离涡流场平衡态转移的影响。研究结果表明:与相同厚度20%以上的对称翼型相比,2%弯度的翼型出现分离涡流场平衡态转移的起始迎角变小2°左右,迎角区间变宽约1°;在厚度相对较小的NACA2416翼型上也发现上述分离涡平衡态转移现象。由此说明翼型弯度在一定程度上促使了分离涡平衡态的转移。     

激波主导流动下壁板的热气动弹性稳定性理论分析

针对激波主导流动下弹性壁板的热气动弹性稳定性分析问题,建立了基于当地活塞流理论的分析模型,并用数值仿真方法来验证其正确性. 首先基于Hamilton原理和Von-Karman大变形理论,建立壁板的热气动弹性运动方程,其中假设壁板受热后温度均匀分布,激波前后区域的气动力模型采用当地一阶活塞流理论;利用Galerkin方法将具有连续参数系统的偏微分颤振方程离散为有限个自由度的常微分方程;基于李雅普诺夫间接法将非线性颤振方程组在平衡位置处进行线化,再用Routh-Hurwits判据来判断线性系统的稳定性,从而来推论出非线性颤振系统的气动弹性稳定性. 在时域中采用龙格--库塔法对非线性颤振方程进行数值积分,得到壁板非线性颤振响应的时间历程,与理论分析结果进行对比. 研究结果表明,壁板受到斜激波冲击时,更容易发生颤振失稳,并且激波强度越大,极限环幅值和频率越大;激波主导流场中的壁板失稳边界不同于传统单纯超声速气流中壁板颤振的失稳边界;只有在斜激波前后不同的动压值都满足颤振稳定性边界的条件下,壁板才可能保持其气动弹性稳定性.      

翼型模型振动的实验研究与气动特性分析

通过在NF-3低速风洞专门研制的翼型模型及相应的俯仰和沉浮振动机构,选用NACA0012翼型进行大迎角下不同频率的振动实验,研究了模型振动平均状态下对其气动力特性的影响情况,并在N-S方程基础上对振动流场进行了初步分析。实验与计算研究的结果表明:在临近定常失速迎角的大迎角条件下,翼型的振动可以引起旋涡分离,导致翼型升力减小和失速迎角的提前;就所讨论的两种振动模式而言,俯仰振动的影响大于沉浮振动。所以,为了提高飞机模型,尤其是大展弦比飞机模型的风洞实验精准度,在模型设计和加工时要特别注意加强机翼弦向的扭转刚度。     

一种层流翼型附面层流动的数值模拟和实验研究

论文通过数值模拟方法对NACA6412层流翼型的低速绕流特性进行了研究,采用修正的BaldwinLomax湍流模型对来流速度为25m/s、8°迎角时翼型的流场、背风区的速度型分布、不同弦向位置附面层厚度的变化规律以及转捩位置进行了模拟,数值模拟结果与风洞实验结果具有一定的相关性.     

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.     

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.     

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.     

FLUTTER OF AN AIRFOIL WITH A CUBIC RESTORING FORCE

In this paper, the effect of a cubic structural restoring force on the flutter characteristics of a two-dimensional airfoil placed in an incompressible flow is investigated. The aeroelastic equations of motion are written as a system of eight first-order ordinary differential equations. Given the initial values of plunge and pitch displacements and their velocities, the system of equations is integrated numerically using a fourth order Runge-Kutta scheme. Results for soft and hard springs are presented for a pitch degree-of-freedom nonlinearity. The study shows the dependence of the divergence flutter boundary on initial conditions for a soft spring. For a hard spring, the nonlinear flutter boundary is independent of initial conditions for the spring constants considered. The flutter speed is identical to that for a linear spring. Divergent flutter is not encountered, but instead limit-cycle oscillation occurs for velocities greater than the flutter speed. The behaviour of the airfoil is also analysed using analytical techniques developed for nonlinear dynamical systems. The Hopf bifurcation point is determined analytically and the amplitude of the limit-cycle oscillation in post-Hopf bifurcation for a hard spring is predicted using an asymptotic theory. The frequency of the limit-cycle oscillation is estimated from an approximate method. Comparisons with numerical simulations are carried out and the accuracy of the approximate method is discussed. The analysis can readily be extended to study limit-cycle oscillation of airfoils with nonlinear polynomial spring forces in both plunge and pitch degrees of freedom.     

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.     

Limit-cycle oscillations in unsteady flows dominated by intermittent leading-edge vortex shedding

High-frequency limit-cycle oscillations of an airfoil at low Reynolds number are studied numerically. This regime is characterized by large apparent-mass effects and intermittent shedding of leading-edge vortices. Under these conditions, leading-edge vortex shedding has been shown to result in favorable consequences such as high lift and efficiencies in propulsion/power extraction, thus motivating this study. The aerodynamic model used in the aeroelastic framework is a potential-flow-based discrete-vortex method, augmented with intermittent leading-edge vortex shedding based on a leading-edge suction parameter reaching a critical value. This model has been validated extensively in the regime under consideration and is computationally cheap in comparison with Navier–Stokes solvers. The structural model used has degrees of freedom in pitch and plunge, and allows for large amplitudes and cubic stiffening. The aeroelastic framework developed in this paper is employed to undertake parametric studies which evaluate the impact of different types of nonlinearity. Structural configurations with pitch-to-plunge frequency ratios close to unity are considered, where the flutter speeds are lowest (ideal for power generation) and reduced frequencies are highest. The range of reduced frequencies studied is two to three times higher than most airfoil studies, a virtually unexplored regime. Aerodynamic nonlinearity resulting from intermittent leading-edge vortex shedding always causes a supercritical Hopf bifurcation, where limit-cycle oscillations occur at freestream velocities greater than the linear flutter speed. The variations in amplitude and frequency of limit-cycle oscillations as functions of aerodynamic and structural parameters are presented through the parametric studies. The excellent accuracy/cost balance offered by the methodology presented in this paper suggests that it could be successfully employed to investigate optimum setups for power harvesting in the low-Reynolds-number regime.     

The post-Hopf-bifurcation response of an airfoil in incompressible two-dimensional flow

A bifurcation analysis of a two-dimensional airfoil with a structural nonlinearity in the pitch direction and subject to incompressible flow is presented. The nonlinearity is an analytical third-order rational curve fitted to a structural freeplay. The aeroelastic equations-of-motion are reformulated into a system of eight first-order ordinary differential equations. An eigenvalue analysis of the linearized equations is used to give the linear flutter speed. The nonlinear equations of motion are either integrated numerically using a fourth-order Runge-Kutta method or analyzed using the AUTO software package. Fixed points of the system are found analytically and regions of limit cycle oscillations are detected for velocities well below the divergent flutter boundary. Bifurcation diagrams showing both stable and unstable periodic solutions are calculated, and the types of bifurcations are assessed by evaluating the Floquet multipliers. In cases where the structural preload is small, regions of chaotic motion are obtained, as demonstrated by bifurcation diagrams, power spectral densities, phase-plane plots and Poincaré sections of the airfoil motion; the existence of chaos is also confirmed via calculation of the Lyapunov exponents. The general behaviour of the system is explained by the effectiveness of the freeplay part of the nonlinearity in a complete cycle of oscillation. Results obtained using this reformulated set of equations and the analytical nonlinearity are in good agreement with previously obtained finite difference results for a freeplay nonlinearity.     

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.     

Analytical formulation of 2-D aeroelastic model in weak ground effect

This paper deals with the aeroelastic modeling and analysis of a 2-D oscillating airfoil in ground effect, elastically constrained by linear and torsional springs and immersed in an incompressible potential flow (typical section) at a finite distance from the ground. This work aims to extend Theodorsen theory, valid in an unbounded flow domain, to the case of weak ground effect, i.e., for clearances above half the airfoil chord. The key point is the determination of the aerodynamic loads, first in the frequency domain and then in the time domain, accounting for their dependence on the ground distance. The method of images is exploited in order to comply with the impermeability condition on the ground. The new integral equation in the unknown vortex distribution along the chord and the wake is solved using asymptotic expansions in the perturbation parameter defined as the inverse of the non-dimensional ground clearance of the airfoil. The mathematical model describing the aeroelastic system is transformed from the frequency domain into the time domain and then in a pure differential form using a finite-state aerodynamic approximation (augmented states). The typical section, which the developed theory is applied to, is obtained as a reduced model of a wing box finite element representation, thus allowing comparison with the corresponding aeroelastic analysis carried out by a commercial solver based on a 3-D lifting surface aerodynamic model. Stability (flutter margins) and response of the airfoil both in frequency and time domains are then investigated. In particular, within the developed theory, the solution of the Wagner problem can be directly achieved confirming an asymptotic trend of the aerodynamic coefficients toward the steady-state conditions different from that relative to the unbounded domain case. The dependence of flutter speed and the frequency response functions on ground clearance is highlighted, showing the usefulness of this approach in efficiently and robustly accounting for the presence of the ground when unsteady analysis of elastic lifting surfaces in weak ground effect is required.     

翼型颤振压电俘能器的输出特性研究

压电俘能器能够为自然界中低功率的微机电系统持续供能. 为了模拟机翼的沉浮?俯仰二自由度运动和有效俘获气动弹性振动能量, 本文提出一种新颖的翼型颤振压电俘能器. 基于非定常气动力模型, 推导翼型颤振压电俘能器流?固?电耦合场的数学模型. 建立有限元模型, 模拟机翼的沉浮?俯仰二自由度运动, 获得机翼附近的涡旋脱落和流场特性. 搭建风洞实验系统, 制作压电俘能器样机. 利用实验验证理论和仿真模型的正确性, 仿真分析压电俘能器结构参数对其气动弹性振动响应和俘获性能的影响. 结果表明: 理论分析、仿真模拟和实验研究获得的输出电压具有较好的一致性, 验证建立数学和仿真模型的正确性. 仿真分析获得机翼附近的压力场变化云图, 表明交替的压力差驱动机翼发生二自由度沉浮?俯仰运动. 当风速超过颤振起始速度时, 压电俘能器发生颤振, 并表现为极限环振荡. 当偏心距为0.3和风速为16 m/s时, 可获得最大输出电压为17.88 V和输出功率为1.278 mW. 功率密度为7.99 mW/cm3, 相比较于其他压电俘能器, 能实现优越的俘获性能. 研究结果对设计更高效的翼型颤振压电俘能器提供重要的指导意义.