Numerical studies of a non-linear aeroelastic system with plunging and pitching freeplays in supersonic/hypersonic regimes

The flutter and post flutter of a two-dimensional double-wedge lifting surface with combined freeplay and cubic stiffness nonlinearities in both plunging and pitching degrees-of-freedom operating in supersonic/hypersonic flight speed regimes have been analyzed. In addition to the structural nonlinearities, the third-order piston theory aerodynamics is used to evaluate the unsteady non-linear aerodynamic force and moment. Such model accounts for stiffness and damping contributions produced by the aerodynamic loads. Responses involving limit cycle oscillation and chaotic motion are observed over a large number of parameters that characterizes the aeroelastic system. Results of the present study show that the freeplay in the pitching degree-of-freedom and soft/hard cubic stiffness in the pitching and plunging degrees-of-freedom have significant effects on the LCOs exhibited by the aeroelastic system in the supersonic/hypersonic flight speed regimes. The simulations also show that the aeroelastic system behavior is greatly affected by physical structural parameters, such as the radius of gyration and the frequency ratio, especially in post-flutter regimes, when accounting for all system nonlinearities. It has been shown that at high Mach numbers the non-linear aerodynamic stiffness yields detrimental effects from the aeroelastic point of view, while the damping one do not.     

Active Aerothermoelastic Control of Hypersonic Double-wedge Lifting Surface

Designing re-entry space vehicles and high-speed aircraft requires special attention to the nonlinear thermoelastic and aerodynamic instability of their structural components. The thermal effects are important since temperature environment brings dramatic influences on the static and dynamic behaviors of flight structures in supersonic/hypersonic regimes and is likely to cause instability, catastrophic failure and oscillations resulting in structural failure due to fatigue. In order to understand the dynamic behaviors of these ＂hot＂ structures, a double-wedge lifting surface with combining freeplay and cubic structural nonlinearities in both plunging and pitching degrees-of-freedom operating in supersonic/hypersonic flight speed regimes has been analyzed. A third order piston theory aerodynamic is used to estimate the applied nonlinear unsteady aerodynamic loads. Also considered is the loss of torsional stiffness that may be incurred by lifting surfaces subject to axial stresses induced by aerodynamic heating. The aerodynamic heating effects are estimated based on the adiabatic wall temperature due to high speed airstreams. As a recently emerging technology, the active aerothermoelastic control is aimed at providing solutions to a large number of problems involving the aeronautical/aerospace flight vehicle structures. To prevent such damaging phenomena from occurring, an application of linear and nonlinear active control methods on both flutter boundary and post-flutter behavior has been fulfilled. In this paper, modeling issues as well as numerical simulation have been presented and pertinent conclusions outlined. It is evidenced that a serious loss of torsional stiffness may induce the dynamic instability; however active control can be used to expand the flutter boundary and convert unstable limit cycle oscillations （LCO） into the stable LCO and/or to shift the transition between these two states toward higher flight Mach numbers.     

基于非线性试验气动力的飞机静气动弹性响应分析

分别使用线性和非线性静气动弹性分析方法,对某飞机纵向静气动弹性响应特性随飞行动压、攻角、平尾偏度和纵向过载变化的趋势进行了分析,并将分析结果和飞行试验进行了比较。线性方法的气动力计算使用亚音速偶极子格网法。非线性方法由风洞试验提供刚体气动力,并使用线性气动力影响系数矩阵对其进行弹性化处理。通过对两种方法的计算结果进行对比分析可以看出:①使用刚体试验气动力的非线性分析方法能够获得和飞行试验比较一致的结果;②在线性方法所提供的结果中,翼面的弯曲变形及剪力、弯矩和扭矩等部分参数能为初步设计提供大致参考,但气动力系数的弹性增量、翼面的扭转变形、翼面的剪力、弯矩和扭矩的弹性部分可能不能为型号设计提供正确的指导;③不论刚体气动力为线性还是非线性,气动力系数随动压变化的趋势均呈线性。     

亚、超音速机翼非定常广义气动力系数之间的关系

 <正> 研究弹性飞机的动力学特性应计及非定常广义气动力。为了尽可能保留非定常运动历程,应首先计算由阶跃函数形式下洗产生的非定常气动力函数(指示函数或指数函数),再利用卷积形式的迭加原理处理任意函数形式下洗产生的非定常气动力。用振型法建立飞机弹性运动方程,可将飞机上任一质点的弹性位移近似地表示为     

飞机大迎角非定常气动力建模研究进展

准确建立非定常气动力数学模型,是飞机大迎角飞行控制律设计、飞行动力学分析和飞行仿真的基础与前提。鉴于此,对大迎角非定常气动力建模研究进展,包括数学建模方法和人工智能建模方法两类进行了系统综述。其中：数学类建模方法是以对非定常流动现象和机理认识为基础的,主要有气动导数模型、非线性阶跃响应模型、状态空间模型、微分方程模型、非线性阶跃响应与状态空间混合模型以及迎角速率模型等;人工智能方法回避了复杂流动机理,属于黑箱非线性系统建模,主要有神经网络模型、模糊逻辑模型和支持向量机模型等。对于每种气动力模型,阐述了其建模思路和方法,给出了典型应用情况,并对其特点和局限性作了简要评述。最后,指出了当前大迎角非定常气动力建模研究工作存在的问题和未来研究方向。     

低速流场中柔性悬臂板的后颤振响应

建立了一个新的非线性气动弹性模型,对低速流场中柔性悬臂板的后颤振响应特性进行了分析。建模中考虑了结构几何非线性、气动力非线性以及两者之间的强耦合效应。通过实验数据对所建立的气动弹性模型进行了验证。发现在低速流场中柔性悬臂板可能会以周期加倍的方式进入混沌运动。结构几何非线性效应和翼尖涡引起的非定常气动力效应对柔性悬臂板的结构响应有显著影响,而尾涡变形引起的非定常气动力对结构运动的影响较小。还研究了不同耦合算法的差异,给出了小展弦比大柔性结构非线性气动弹性数值仿真时耦合策略的选择依据。     

一种基于亚临界响应的颤振稳定性边界预测新方法

基于地面模态试验事先获取的结构模态参数,通过弹性机翼在气流中的亚临界响应,解算出作用在机翼上的非定常模态气动力系数。以弹性机翼的模态位移作为系统输入,模态气动力系数作为系统输出,通过辨识方法获得弹性机翼振动的气动力模型。在时域内耦合结构运动方程和气动力模型,建立基于试验数据的气动弹性稳定性分析模型。通过分析系统稳定性随动压的变化规律,获得弹性机翼的颤振稳定性特性。与经典颤振边界外推方法的主要区别在于该方法实质上只需要一次亚临界响应试验即可预测颤振临界点,可极大降低颤振试验的风险和成本。该方法既可用于颤振风洞试验,也可用于颤振试飞。     

弹性耦合对复合材料旋翼前飞气弹响应及载荷的影响

 给出一套计算复合材料旋翼前飞气弹响应和桨榖载荷的方法,所用结构模型考虑横向剪切变形、剖面面外翘曲变形和复合材料弹性耦合的影响。气动模型采用准定常升力线理论和Drees线性入流模型,翼型升力、阻力系数来自风洞试验。构造出21自由度梁单元,应用Hamilton 原理推导出桨叶运动的有限元方程。在此基础上,研究弯曲-扭转和拉伸-扭转耦合对复合材料旋翼前飞气弹响应和桨榖振动载荷的影响。结果表明：弹性耦合对扭转方向的影响很大,对挥舞和摆振两个方向的影响很小;各种弹性耦合对桨榖振动载荷有着不同程度的影响,负的摆振弯曲 扭转耦合和正的拉伸 扭转耦合使桨榖垂直方向的振动载荷降低5%左右。     

直升机旋翼瞬态气弹运动有限元模型研究

从能量变分原理出发建立了旋翼的瞬态气弹有限元模型,模型考虑了旋翼转速加速度对桨叶动能变分的影响。采用叶素理论建立了桨叶在复杂流场中的气动载荷模型,模型计入了旋翼主轴纵倾和横滚角的影响,研究了流场风速分量的计算。运用所建模型仿真分析了直升机在护卫舰尾流场内起动时桨叶的挥舞、摆振、扭转(变距)等瞬态运动过程。模型可用于旋翼的振动仿真分析、噪声估计和故障诊断方法的验证。     

舰面直升机旋翼瞬态气弹响应分析

建立了舰面直升机起动和停车过程中的旋翼瞬态气弹响应分析方法,在桨叶结构建模中考虑了桨叶的挥舞、摆振、扭转和拉伸变形并计入了桨叶重量影响和因旋翼速度变化产生的惯性力。桨叶气动力计算采用时域非定常空气动力计算模型,包括非线性分离和动态失速的影响。用实验数据对分析方法进行了验证。      

总距突增时共轴式直升机瞬态操纵响应分析

从共轴双旋翼直升机的工程实际出发,建立了共轴式直升机上下旋翼非定常气动特性的计算模型.引入Leishman-Beddoes指数函数的半经验公式,建立了二维翼型非定常气动模型;分别从固定尾迹和自由尾迹,引入干扰因子到动态入流三种方法出发,建立了反映共轴双旋翼直升机上下旋翼气动干扰的诱导速度模型;从跷跷板式旋翼的挥舞动力学方程出发,利用4阶Runge-Kutta算法求解桨叶刚性挥舞角的数值解.通过计算分析,得到了悬停和前飞状态,总距突增时上下旋翼升力的动态响应特性,以及总距突增时上下旋翼桨叶铰链力矩的响应特性.      

Adaptive aeroelastic vibration suppression of a supersonic airfoil with flap

This paper concerns the flutter, post-flutter and adaptive control of a non-linear 2-D wing-flap system operating in supersonic/hypersonic flight speed regimes. An output feedback control law is implemented and its performance toward suppressing flutter and limit cycle oscillations (LCOs) as well as reducing the vibrational level in the subcritical flight speed range is demonstrated. This control law is applicable to minimum phase systems and we provide conditions for stability of the zero dynamics. The control objective is to design a control strategy to stabilize the pitch angle while adaptively compensating for uncertainties in all the aeroelastic model parameters. It is shown that all the states of the closed-loop system are asymptotically stable.     

Transient aeroelastic responses and flutter analysis of a variable-span wing during the morphing process

To investigate the transient aeroelastic responses and flutter characteristics of a variablespan wing during the morphing process,a novel frst-order state-space aeroelastic model is proposed.The time-varying structural model of the morphing wing is established based on the Euler-Bernoulli beam theory with time-dependent boundary conditions.A nondimensionalization method is used to translate the time-dependent boundary conditions to be time-independent.The time-domain aerodynamic forces are calculated by the reduced-order unsteady vortex lattice method.The morphing parameters,i.e.,wing span length and morphing speed,are of particular interest for understanding the fundamental aeroelastic behavior of variable-span wings.A test case is proposed and numerical results indicate that the flutter characteristics are sensitive to both of the two morphing parameters.It could be noticed that the aeroelastic characteristics during the wing extracting process are more serious than those during the extending process at the same morphing speed by transient aeroelastic response analysis.In addition,a faster morphing process can get better aeroelastic performance while the mechanism comlexity will arise.     

压电驱动器的气动弹性应用

 随着压电智能材料与结构的发展,压电驱动器在气动弹性控制领域占据重要地位。使用压电驱动器控制翼面变形,利用而不是抵抗气动弹性效应可以控制升力、力矩以及它们的分布。采用基本相同的智能结构翼面控制系统,根据不同的控制目标需求,使用压电智能材料驱动器可以达到多种目的,包括静态的形状控制与动态的颤振抑制、抖振控制与阵风响应控制。静态控制方面例如改变翼面形状获得附加空气动力以增加升力、提供横滚力矩、改变升力分布以减小诱导阻力或减小翼根弯矩等;动态控制例如利用改变翼面形状产生的附加空气动力作为控制载荷,改变气动弹性系统的耦合程度,根据控制效果要求可作为气动阻尼、气动刚度或气动质量。这种控制方法可以减轻结构重量,提高操纵效率,扩大飞行包线,提高材料利用率,已成为可变形飞行器的重要研究内容。本文主要阐述压电驱动器气动弹性应用的动机与机理、发展与成就以及问题与展望。     

Active aeroelastic control over a four control surface wing model

This paper describes the procedure implemented to design, develop and test an aeroelastic control system installed on a forward swept wing of the aeroelastic demonstrator X-DIA. A control method directly based on Nissim aerodynamic energy concept has been chosen. Two different modeling techniques have been adopted for the calculation of generalized aerodynamic forces, such as doublet lattice method and computational fluid dynamics and the obtained results are finally compared. The latter approach, applied to better estimate the control surfaces effectiveness, requires the capability to correctly model the control surface rotation and the grid deformation, usually addressed as non-trivial problems in CFD based aeroelastic analysis. A genetic algorithm optimization technique has been adopted to state and refine all the control gains.     

基于CFD/CSD耦合方法的旋翼气动弹性载荷计算分析

为提高直升机前飞状态下旋翼非定常气动弹性载荷的预估精度,在旋翼气动弹性综合分析方法中引入旋翼CFD模块,建立了一套基于CFD/CSD松耦合分析的计算方法和程序。为高效解决流固耦合方法中由于桨叶挥舞、扭转等弹性变形带来的旋翼贴体网格变形问题,采用基于代数变换方法的网格变形技术,桨叶运动变形量和旋翼气动力信息通过流固交接面传递。旋翼流场分析方法中,主控方程采用耦合S-A湍流模型的Navier-Stokes方程,围绕旋翼流场的网格采用结构嵌套网格方法生成,无黏通量计算采用Roe格式,时间推进采用双时间法。旋翼结构分析中,考虑旋翼配平,基于Hamilton变分原理和20自由度Timoshenko梁模型求解弹性旋翼非线性运动方程。分别对CSD和CFD方法进行验证,在此基础上,计算了SA349/2旋翼桨叶在前飞状态下的非定常气动力、挥舞弯矩、摆振弯矩和扭转力矩,并与飞行测试数据进行了对比。计算表明:CFD/CSD耦合方法可以显著提高旋翼非定常气动弹性载荷的分析精度,精确捕捉桨叶表面压强峰值、激波位置等,表明本文发展的旋翼CFD/CSD耦合方法可以有效地运用到旋翼气动弹性载荷的预测分析中。     

刚弹耦合飞行器的机动载荷减缓

针对飞行器刚体/弹性模态耦合作用下的机动载荷减缓问题,运用亚音速非定常空气动力学理论以及分析力学原理,建立了柔性飞机机动载荷分析及控制的气动弹性模型,分析了飞机作纵向机动飞行时结构弹性对机动飞行参数的影响。数值分析结果表明,结构弹性模态对飞机机动过程的影响不容忽视。为了有效减缓机动载荷,运用最优控制理论,采用多组控制面联合偏转的作动方案,设计机动载荷减缓最优控制律,实现了对飞机纵向机动载荷的控制。研究结果表明,提出的机动载荷减缓控制方案是有效的。     

小型共轴式直升机旋翼桨叶铰链力矩研究

从工程实际出发,建立了共轴双旋翼直升机桨叶铰链力矩计算模型.引入Leishman-Beddoes指数函数的半经验公式,建立了二维翼型非定常气动模型;引入干扰因子到Pitt-Peters动态入流模型,建立了反映共轴双旋翼直升机上下旋翼气动干扰的诱导速度模型;从桨叶的挥舞动力学模型出发,利用四阶Runge-Kutta算法求解桨叶刚性挥舞角的数值解,并利用模态叠加法计算桨叶的挥舞向变形.最后,通过对在研共轴直升机的计算得知,上下旋翼桨叶铰链力矩虽然在分量上有所差异,但是总的铰链力矩基本相同;中高速飞行时,上下旋翼桨叶铰链力矩在下倾斜盘累积成为使下倾斜盘前倾的力矩.      

Unsteady aerodynamics and flow control for flapping wing flyers

The creation of micro air vehicles (MAVs) of the same general sizes and weight as natural fliers has spawned renewed interest in flapping wing flight. With a wingspan of approximately 15 cm and a flight speed of a few meters per second, MAVs experience the same low Reynolds number (10⁴–10⁵) flight conditions as their biological counterparts. In this flow regime, rigid fixed wings drop dramatically in aerodynamic performance while flexible flapping wings gain efficacy and are the preferred propulsion method for small natural fliers. Researchers have long realized that steady-state aerodynamics does not properly capture the physical phenomena or forces present in flapping flight at this scale. Hence, unsteady flow mechanisms must dominate this regime. Furthermore, due to the low flight speeds, any disturbance such as gusts or wind will dramatically change the aerodynamic conditions around the MAV. In response, a suitable feedback control system and actuation technology must be developed so that the wing can maintain its aerodynamic efficiency in this extremely dynamic situation; one where the unsteady separated flow field and wing structure are tightly coupled and interact nonlinearly. For instance, birds and bats control their flexible wings with muscle tissue to successfully deal with rapid changes in the flow environment. Drawing from their example, perhaps MAVs can use lightweight actuators in conjunction with adaptive feedback control to shape the wing and achieve active flow control. This article first reviews the scaling laws and unsteady flow regime constraining both biological and man-made fliers. Then a summary of vortex dominated unsteady aerodynamics follows. Next, aeroelastic coupling and its effect on lift and thrust are discussed. Afterwards, flow control strategies found in nature and devised by man to deal with separated flows are examined. Recent work is also presented in using microelectromechanical systems (MEMS) actuators and angular speed variation to achieve active flow control for MAVs. Finally, an explanation for aerodynamic gains seen in flexible versus rigid membrane wings, derived from an unsteady three-dimensional computational fluid dynamics model with an integrated distributed control algorithm, is presented.     

Aeroelastic system identification using transonic CFD data for a wing/store configuration

This paper is part of a study investigating the prediction of the aeroelastic behaviour of aircraft subjected to non-linear aerodynamic forces. The main objective of the work is the characterization of the dynamic response of aeroelastic models resulting from coupled Computational Fluid Dynamic and Finite Element calculations. Of interest here is the identification of the flight condition at which the response bifurcates to limited or divergent amplitude self-sustained oscillations without carrying out a comprehensive set of full, computationally expensive, time-marching calculations. The model treated in this work is a three-dimensional wing in a transonic flowfield. Short datasets of pre-bifurcation behaviour are analysed to determine the system's stability and degree of non-linearity. It is found that the calculated responses on the run-up to a transonic Limit Cycle Oscillation show little or no evidence of non-linearity. The non-linearity appears abruptly at the bifurcation flight condition. The variation of the local Mach number over the wing's surface in the steady-state case is used to demonstrate that the non-linearity is due to a shock wave that can move along the surface. At Mach numbers where this is not possible the system behaves in a linear manner and its stability can be analysed using linear methods.