一种新型高鲁棒性动网格技术及其应用

首先对四元数进行李代数空间指数映射,解决了多个四元数插值问题,并结合距离倒数插值方法实现网格边界扰动向空间网格的传播,建立了新型高鲁棒性的四元数变形网格技术. 针对该型动网格技术中由于大型矩阵运算量引起的运算效率低问题,同时利用四元数方法在动网格变形中具备与物面边界高阶一致性的特点,提出了分层次变形策略,避免了面向全流场网格节点的大型矩阵运算;进一步基于无限插值技术较强的逻辑保持能力,建立了面向结构网格分层混合变形方法. 充分利用多区域重叠、对接网格变形技术中隐含的并行性,基于对等式编程思想及MPI 库函数对动网格程序进行并行化编程,建立了高效高鲁棒性的变形网格技术. 以某型客机翼身组合体气动弹性分析为范例,研究了不同方法之间的计算效率以及鲁棒性,进一步将分层混合变形网格技术应用于某型支线客机全机型架外形设计与修正,验证了所建立的新型动网格技术的高效性与鲁棒性.      

直升机气动弹性力学发展现状(续)

Ⅲ.单片桨叶气动弹性问题的求解1.桨叶离散化方法求解旋翼桨叶气动弹性力学问题的第一步是将连续桨叶离散化,即把一个具有无限多个自由度的连续参数系统离散化为具有有限个自由度的离散系统。常用方法有三种: 1)整体模态方法在直升机旋翼气动弹性力学中,以往用得多的离散方法是整体模态法,或称为整体伽辽金方法。此方法的实质是利用桨叶自由振动振型是线性独立      

STATIC AEROELASTICITY METHOD BASED ON EXTERNAL AERODYNAMIC FORCE AND STRIP METHOD

本文通过对载荷设计中的静气动弹性分析方法进行研究,发展了一种基于外部刚性气动力数据和改良片条理论修正的弹性载荷修正方法。其中结构变形通过工程梁方法提取部件刚度阵进行计算,气动力通过网格化模型基于外部试验或CFD气动力数据库插值得到。结构和气动之间位移和力的数据传递分别利用曲面样条和形函数面积坐标加权法插值,弹性载荷修正通过改良后的片条理论计算,由此迭代循环直至结构变形收敛。同时通过对相关工程实例进行分析计算并与成熟方法对比,验证了该方法的可行性和高效性。     

流固耦合力学概述

本文简要介绍了流固耦合力学及其特点、研究分支、一些进展及进一步发展的趋势     

Numerical Continuation Applied to Panel Flutter

The panel flutter formulation studied here includes a geometricnonlinearity. In a physical sense, the nonlinearity represents couplingbetween out-of-plane bending and in-plane stretching. To investigate theinitial condition response, the pseudo-arclength continuation method isapplied to this problem. The continuation method is an alternatenumerical method that complements direct numerical integration. Thelatter procedure is not necessarily computationally efficient when theparameter space under consideration is large. One strength ofcontinuation methods then is the ability to efficiently and accuratelydisplay the solution space, including certain bifurcations, for dynamicsystems. Results from the procedure are in substantial agreement withearlier studies, where applicable, where the panel undergoes cylindricalbending.     

Aeroelastic inverse: Estimation of aerodynamic loads during large amplitude limit cycle oscillations

This paper presents an algorithm to compute the aerodynamic forces and moments of an aeroelastic wing undergoing large amplitude heave and pitch limit cycle oscillations. The technique is based on inverting the equations of motion to solve for the lift and moment experienced by the wing. Bayesian inferencing is used to estimate the structural parameters of the system and generate credible intervals on the lift and moment calculations. The inversion technique is applied to study the affect of mass coupling on limit cycle oscillation amplitude. Examining the force, power, and energy of the system, the reasons for amplitude growth with wind speed can be determined. The results demonstrate that the influence of mass coupling on the pitch–heave difference is the driving factor in amplitude variation. The pitch–heave phase difference not only controls how much aerodynamic energy is transferred into the system but also how the aerodynamic energy is distributed between the degrees of freedom.     

Numerical investigation of structural geometric nonlinearity effect in high-aspect-ratio wing using CFD/CSD coupled approach

An efficient and robust fluid–structure coupled methodology has been developed to investigate the linear and non-linear static aeroelastic behavior of flexible high-aspect-ratio wing. A three-dimensional open source finite element solver has been loosely coupled with an in-house Reynolds-averaged Navier–Stokes solver, designed for hybrid-unstructured meshes, to perform aero-structural coupled simulations. For volume mesh deformation and two-way data interpolation over non-matching grids interface, a radial basis function methodology combined with a data reduction algorithm has been used. This technique is efficient in handling large deflections and provides high-quality deformed meshes. Structural geometric nonlinearity has been considered to predict the deformations in the vertical and torsional directions caused by gravitational and aerodynamic loading. A multi-material finite element model has been generated to match the experimental configuration. Computational aeroelastic simulations were performed on an experimental high-aspect-ratio aeroelastic wing model with a slender body at the tip to get non-linear static deflections, twist and structure natural frequencies. The effect of the geometric nonlinearity is significant for large deformation analysis and has been highlighted in the predicted maximum tip deflection and twist. Good qualitative and quantitative agreement has been achieved between the predicted results and the available experimental data.     

基于边界元方法的气动弹性降阶模型及仿真

传统气动弹性的时域计算耗费了大量时间,为了提高计算效率,本文发展了基于边界元方法的降阶模型技术。首先基于边界元方法建立非定常流场的求解模型,结合特征值分析技术建立了非定常气动力的低阶模型;然后,利用边界元方法建立了气动网格和结构网格之间的信息转换矩阵;最后将非定常气动力降阶模型和结构动力学方程联合,建立了气动弹性系统的低阶状态空间模型。将所发展的降阶模型方法应用于NACA0012翼型的非定常气动力求解中,结果表明降阶模型可以在保证原系统计算精度的同时提高了计算效率;将降阶模型技术应用到三维机翼的气动弹性响应计算中,在系统阶数仅为12阶的情况下可以得到与原系统一致的极限环响应,说明降阶模型技术在求解气动弹性问题中的巨大优势。     

A new frequency domain based approach for a nonlinear aeroelasticity analysis

In this paper, a new frequency-domain based approach for the investigation of aeroelasticity problems is introduced that is capable of handling both linear and nonlinear problems. This approach is based on coupling a conceptual method used in a structural dynamic analysis and an optimum equivalent linear frequency response function (OELF). Additionally, a new criterion for determining the flutter speed and the instability of nonlinear systems is introduced that is based on the condition number of the aeroelastic matrices. Due to the global nature of the condition number, the new criterion proves to be efficient and simple to use. To examine the efficiency of the new technique, a two-dimensional nonlinear airfoil with an unsteady aerodynamic model is considered.