一种基于CFD技术的高效气动弹性时域求解方法.pdfVIP

第十届全国空气弹性学术交流会 CARS．200r750002 一种基于CFD技术的高效气动弹性时域求解方法 蒋跃文张伟伟叶正寅 (西北工业大学翼型叶栅空气动力学国防科技重点实验室，陕西西安710072) 摘要针对流场，结构耦合方程，基于Ad锄s隐式方法发展了一种杂交的线性多步法(mJ订S)，广义气动力 项采用插值求解，结构项仍采用隐式格式计算，融合了显式方法耦合简单和隐式方法稳定性好的特点。该时 域耦合迭代格式使用预估一校正法进行时间推进，仅需在预估步调用一次流场求解器，校正步的广义气动力 项采用插值求解。算例运用非定常Eulcr方程计算了一个标准气动弹性模型(Isogaiwing)的颤振结果，并 与几种经典的时域模拟方法进行了比较，证明该方法具有效率高、稳定性好、精度高的优点。 关键词气动弹性颤振时域模拟线性多步杂交的线性多步法学习 ’ 1引言 由于跨声速、大迎角等飞行范围的工程问题需要，以及传统线性气动弹性分析方法的不足，近 年来，非线性流场和结构运动方程的数值耦合求解技术在气动弹性领域得到很大的发展。气动力求 解方法已经由线化理论为基础的方法【l刁】发展到现在的CFD技术陆101，但是求解非定常Euler小．S方 程的CFD技术应用于非线性气动弹性时域分析的计算量非常大，耗费的时间长。因此如何提高气动 弹性计算的效率和准确度一直是该领域研究的重点问题。 文献[4]采用流场和结构进行同步积分的方法，在每个时间步内结构和气动力之间信息不断交换 直到收敛。文献[5]流场和结构耦合是在每个有限控制体内进行。但这两种算法需对流场求解器代码 作较大改动，使用不便，也不便于并行运算。 在进行时域推进时，程序模块性比较好的一种方法是在流体和结构区域分别创建网格和物理模 型，流场和结构参数在每个时间步进行信息交换，同步推进。这种耦合算法最大的好处是在满足精 确耦合的同时能充分利用现有的计算流体力学(cFD)技术和计算结构动力学的方法，只需增加少 量的数据交换模块即可，从而保持了程序的模块化，所以使用方便。 目前，冻结气动力的龙格—库塔法慨使用的也比较多。该方法虽然算法简单，但精度低，时间 推进时需要选取很小的时间步长才能得到较准确的解，计算量较大。文献[8]使用插值来近似求解龙 格一库塔法在每个时间步内变化的模态气动力，虽然算例显示该技术可大大提高计算的精度，但缺乏 数学上的误差分析。 线性多步法在气动弹性时域分析中也得到了广泛利用。对于气动弹性控制方程，Borl锄d和 Rizz眦a【lJ使用中心差分方法求解。而文献【2】使用四阶预估—校正方法求解，文中提到此方法由于每 引入了状态方程表达式，利用显式方法求解，虽然这种方法应用方便但稳定性不高。Robinson等【lo】 使用预估一校正方法求解，即在一个时间步内求解了两次气动力。虽然该方法的精度和稳定性较好， 但是因此增加了计算量，效率不高，而且需要修改流场求解器的代码，程序的模块化较差。 针对以上计算方法在效率、精度、稳定性以及程序的模块性等方面难以统一的问题，发展了一 种杂交的线性多步法(刚S)，利用AId锄s显式方法作为其预估步，使用预估一校正方法来进行 气动弹性时域推进，并通过计算一个跨音速气动弹性模型(Isogaiwing)的颤振结果来与几种经典 ．6． 第十届全国空气弹性学术交流会 CARs．200r750002 的算法进行比较。 2气动弹性控制方程 应用拉格朗日方程，流场中结构的运动方程写成矩阵形式 M蛆+G寸+Kq=Q 其中M为广义质量矩阵，G为广义阻尼矩阵，K为广义刚度系数矩阵，Q为广义气动力矩阵，q为 结构运动的广义位移。 引入状态变量x=嘞，92。…，圣l，圣2，…，，由于Q与x’f有关，把方程(1)改写成状态方程 文=F(x，f)=A·x+B·Q(x，f) (2) 其忆[_二一划， (2)式就是一个一阶常微分方程组。 B=瞄]· 对于二元机翼的颤振问题，设无因次时间。一％一，则有