航空学报 > 2013, Vol. 34 Issue (7): 1501-1511   doi: 10.7527/S1000-6893.2013.0122

基于径向基函数的网格变形及非线性气动弹性时域仿真研究

谢亮1, 徐敏1, 安效民1, 蔡天星2, 陈韦贤1   

  1. 1. 西北工业大学 航天学院, 陕西 西安 710072;
    2. 上海机电工程研究所, 上海 201109
  • 收稿日期:2012-12-19 修回日期:2013-02-02 出版日期:2013-07-25 发布日期:2013-03-04
  • 通讯作者: 徐敏 E-mail:xumin@nwpu.edu.cn
  • 作者简介:谢亮,男,博士研究生。主要研究方向:计算流体力学与计算气动弹性力学。Tel:029-88494614,E-mail:jjccblws02@163.com;徐敏,女,博士,教授,博士生导师。主要研究方向:气动弹性力学、多力学场耦合问题。Tel:029-88494614,E-mail:xumin@nwpu.edu.cn;安效民,男,博士,讲师。主要研究方向:非线性气动弹性力学。Tel:029-88494614,E-mail:frank805@nwpu.edu.cn;蔡天星,男,博士,工程师。主要研究方向:飞行器气动弹性力学Tel:021-24185234,E-mail:star.cai@163.com;陈韦贤,男,硕士研究生。主要研究方向:非线性气动弹性问题。Tel:029-88494614,E-mail:378673941@qq.com
  • 基金资助:

    国家自然科学基金(90816008, 11202165)

Research of Mesh Deforming Method Based on Radial Basis Functions and Nonlinear Aeroelastic Simulation

XIE Liang1, XU Min1, AN Xiaomin1, CAI Tianxing2, CHEN Weixian1   

  1. 1. College of Astronautics, Northwestern Polytechnical University, Xi’an 710072, China;
    2. Shanghai Electro-Mechanical Engineering Institute, Shanghai 201109, China
  • Received:2012-12-19 Revised:2013-02-02 Online:2013-07-25 Published:2013-03-04
  • Supported by:

    National Natural Science Foundation of China (90816008, 11202165)

摘要:

为开展非线性气动弹性研究,基于非线性结构有限元软件NASTRAN和自主研制的多块结构化计算流体力学(CFD)求解器,开发了一套基于计算流体力学/计算结构动力学(CFD/CSD)耦合求解方法的气动弹性时域仿真程序。该程序采用径向基函数(RBF)交换两套求解器之间的数据并进行网格变形。为提高RBF方法的效率,构造了基于多次插值的空间待插值点精简算法。在多次插值过程中,每次插值的对象为上次插值的误差,并同时限制插值区域,以此实现了空间待插值网格数的精简。数个网格变形的算例表明该方法可支持大变形运动,并且具有较高的计算效率。采用此程序开展了AGARD 445.6机翼颤振计算、大展弦比机翼的静气动弹性计算与切尖三角翼极限环振荡(LCO)现象的动气动弹性仿真,结果揭示了当机翼展弦比较大或者响应幅值较大时,结构非线性对于气动弹性有显著影响。

关键词: 径向基函数, 网格变形, 插值, 非线性, 气动弹性, 大展弦比, 三角翼, 极限环振荡

Abstract:

To carry out research of nonlinear aeroelasticity, this paper presents a nonlinear time-domain aeroelastic simulation program which adopts computational fluid dynamics/computational structural dynamics (CFD/CSD) couple method. The program couples the nonlinear structural finite element software NASTRAN with a multi-block structured grid-based CFD solver developed by myself. Radial basis function (RBF) is used to transfer information between the two solvers and deform grid. To improve the efficiency of the RBF method, a space point reduction algorithm based on multi-level interpolation is constructed. For the purpose of reducing the number of space points, the interpolation object of each interpolation step is the error of previous step, and interpolation region is limited simultaneously. Results of some mesh deformation examples show that this way can support the movement of large deformation with high computational efficiency. With the use of this program, flutter computation of AGARD 445.6 wing is completed, static aeroelastic computation of a wing with a large aspect ratio is performed and a dynamic aeroelastic simulation of limit-cycle-oscillation (LCO) of a cropped delta wing is carried out. Results demonstrate that structural nonlinearity has strong influences on aeroelastic phenomenas when the aspect ratio of the wing is large or the amplitude of response is great.

Key words: radial basis function, grid deformation, interpolation, nonlinearity, aeroelasticity, large aspect ratio, delta wing, limit-cycle-oscillation

中图分类号: