一种几何大变形下的非线性气动弹性求解方法

非线性气动弹性的时域求解中, 涉及到非线性的流体动力学(CFD)和非线性的结构动力学(CSD)耦合问题. 基于Co-rotational理论, 推导了三维壳单元几何非线性下的切线刚阵和内力公式, 针对推进过程中的能量守恒, 引入预估-校正推进格式, 发展了一种近似能量守恒的非线性动态响应算法; 基于1/2时间步的交错耦合格式, 结合带有几何守恒律的双时间推进求解雷诺平均N-S方程的求解器, 发展了非线性气动弹性求解的高精度耦合格式. 通过结构几何大变形下的静力和动力分析验证了所发展的结构非线性求解器, 并通过AGARD445.6机翼的非线性气动弹性响应分析, 说明了所发展耦合求解方法的实用性.      

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

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

Performance enhancement of wing-based piezoaeroelastic energy harvesting through freeplay nonlinearity

We investigate experimentally how controlled freeplay nonlinearity affects harvesting energy from a wing-based piezoaeroelastic energy harvesting system. This system consisits of a rigid airfoil which is supported by a nonlinear torsional spring (freeplay) in the pitch degree of freedom and a linear flexural spring in the plunge degree of freedom. By attaching a piezoelectric material (PSI-5A4E) to the plunge degree of freedom, we can convert aeroelastic vibrations to electrical energy. The focus of this study is placed on the effects of the freeplay nonlinearity gap on the behavior of the harvester in terms of cut-in speed and level of harvested power. Although the freeplay nonlinearity may result in subcritical Hopf bifurcations (catastrophic for real aircrafts), harvesting energy at low wind speeds is beneficial for designing piezoaeroelastic systems. It is demonstrated that increasing the freeplay nonlinearity gap can decrease the cut-in speed through a subcritical instability and gives the possibility to harvest energy at low wind speeds. The results also demonstrate that an optimum value of the load resistance exists, at which the level of the harvested power is maximized.     

Passive morphing of flying wing aircraft: Z-shaped configuration

High Altitude, Long Endurance (HALE) aircraft can achieve sustained, uninterrupted flight time if they use solar power. Wing morphing of solar powered HALE aircraft can significantly increase solar energy absorbency. An example of the kind of morphing considered in this paper requires the wings to fold so as to orient a solar panel to be hit more directly by the sun's rays at specific times of the day. An example of the kind of morphing considered in this paper requires the wings to fold so as to orient a solar panel that increases the absorption of solar energy by decreasing the angle of incidence of the solar radiation at specific times of the day. In this paper solar powered HALE flying wing aircraft are modeled with three beams with lockable hinge connections. Such aircraft are shown to be capable of morphing passively, following the sun by means of aerodynamic forces and engine thrusts. The analysis underlying NATASHA (Nonlinear Aeroelastic Trim And Stability of HALE Aircraft), a computer program that is based on geometrically exact, fully intrinsic beam equations and a finite-state induced flow model, was extended to include the ability to simulate morphing of the aircraft into a “Z” configuration. Because of the “long endurance” feature of HALE aircraft, such morphing needs to be done without relying on actuators and at as near zero energy cost as possible. The emphasis of this study is to substantially demonstrate the processes required to passively morph a flying wing into a Z-shaped configuration and back again.     

Applications and developments of aeroelasticity of flexible structure in flow controls

柔性结构与空气动力耦合形成的系统呈现出丰富的非定常、非线性流动和结构动力学行为,对其气动弹性效应合理地控制和利用,能够大幅度提高飞机机翼、风力机叶片等结构的气动性能,并使其具有气动自适应能力.本文总结了近年来与气弹效应应用相关的研究进展及存在的问题,具体介绍了薄膜翼型的流动控制特性、柔性壁面减阻技术以及Sinha扰流装置的发展过程、主要成果以及未来发展趋势,着重对相关试验、流固耦合数值分析、Lagrangian拟序结构动力学等理论分析方法进行总结,展示了气弹效应在流动控制方面的巨大潜力和深远的学术意义,以便更多的研究人员开展该领域的研究工作.     

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

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

基于RBF插值技术的CFD/CSD非线性耦合分析方法研究

基于非结构混合网格的N-S方程求解器和结构柔度影响系数法,发展了一种考虑气动、结构非线性的基于RBF插值技术CFD/CSD耦合分析方法,适用于解决现代大展弦比飞机的非线性静气动弹性问题。该方法采用时间相关法（即求解非定常方程组,用长时间的渐近解趋于定常状态）求解静气弹分析时的定常流动。考虑大展弦比飞机结构变形问题为大变形小应力问题,在利用柔度系数法求解结构方程时,假设每次求解结构方程时应力与应变为线性关系,整体静气弹分析过程为非线性关系,因此每次求解结构方程时要更新柔度影响系数矩阵。在非定常N-S方程每求解一个时间步耦合一次结构有限元分析,由于结构有限元分析的时间相对于气动分析时间是很短的,所以这种方法实际上近似使用了一次求解非定常气动力的时间完成了整个静气动弹性分析的过程。对于气动网格与结构有限元网格不一致性,本文采用径向基函数（RBF）插值方法中的TPS方法进行结构弹性变形和气动载荷插值,采用虚功原理完成气动载荷数据交换。为了节省气弹分析时间,采用动网格方法对气动网格进行更新,本文基于RBF插值方法发展一种适用于混合网格（四面体、三棱柱、金字塔和六面体）变形的动网格方法,可以保证附面层网格的质量与分布从而准确模拟其流动。利用该方法对M6机翼、DLR-F6翼身组合体和某大型客机机翼进行了静气动弹性特性分析,结果验证了本文开发的非线性CFD/CSD耦合分析方法的可行性、精确性和高效性。     

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.     

直升机旋翼／机身耦合系统的气弹响应分析（二）方程的求解

在文（一）得到的直升机旋翼系统运动方程的基础上,运用动态子结构的方法对旋翼系统和机身作为相对独立的部分进行分析,通过界布的力平衡和几何协调将子系统耦合系统整体系统,用分割－迭代法求解直至旋翼系统和机身的响应同时收敛到精度要求,并研制相应的计算程序,给出了工程算例。