高雷诺数下翼型的气动弹性特性研究.pdfVIP

第六届两北地区计算物理学术会议论文集 2008 高雷诺数下翼型的气动弹性特性研究去 叶正寅张陈安 (西北T业大学，翼型、rtl4栅空气动力学国防科技重点实验室，西安，邮编710072) 摘要通过求解雷诺甲均Navier-Stokes方程并结合修正的B—L代数湍流模型获得定常和非定常气动力，耦合翼 型弹性振动方程，在时间域内数值模拟丫NACA0012翼型和不同迎角下的气动弹性静、动态演化过程，重点研 究r大迎角情况下翼型流场及其气动弹性性质。研究结果表明，在线性气动力范围内，迎角对颤振临界速度的影 响越小；但在人迎角情况下，动态响应对来流速度变得更加敏感，可能在来流速度或初始扰动稍微变化的条件下， 动态响应从收敛突然出现发散，从流场特征看，引起此种突变的根本原冈在于翼型人分离涡的{lj现。 关键词：气动弹性，分离流，失速，颤振 引言 翼硝是构成机翼的基本元素，对于大型飞机而言，机翼一般具有大展弦比、较小后掠角的特 点，所以翼型流动特征在此类飞机的机翼流场中具有代表性，在空气动力学分析中，片条理论等 方法就是以翼璎为基础分析此类机翼的气动性质。另一方面，大型飞机的机翼雷诺数很高，一般 可以达到107以上。冈此，研究高雷诺数下翼型的气动弹性性质，对于认识大型l毛机的气动弹性 问题具有较大的参考价值。 从巡航飞行的角度讲，大型飞机的飞行迎角一般较小，但是，由于大型lS机、尤其是民用飞 机对飞行安全性的要求很高，如比用飞机就要求进行大迎角深失速的气动性质分析，而颤振性质 对飞机的安全性影响更加突出，因此，了解高雷诺数下翼犁的气动弹性性质是一项重要的研究工 作。对于不同的问题，大迎角下的翼璎气动弹性问题一直受到人们的关注ll圳，由于传统的气动 弹性问题研究中，一般使用势流理论下的非定常气动力，适应的迎角范嗣和马赫数范围都有限， 随着计算流体力学(CFD)技术和计算机的发展，近年来采用CFD方法研究气动弹性的T作成 为主要研究途径降‘0l，但是，由于三维气动弹性分析、尤其是涉及到大迎角分离流等非线性气动 弹性问题，时间历程需要很长的数值模拟才能真正了解动态性质，计算量巨大，所以根据大型电 机的机翼特点，选择翼型问题研究具有良好的代表性。本文目的就是系统分析高雷诺数翼型在不 同迎角下的气动弹性性质，尤其是人迎角下翼型的特点。希望深入认识高霄诺数翼型气动弹性特 性的基本规律。 l、气动力计算方法 如上所述，由于计算量的原因，本文以翼型的气动弹性问题代表人型飞机机翼的气动弹性特 征。对于翼型的简单外形而言，为了保证粘性流动特征得到很好模拟，同时考虑到非定常重新生 成网格的效率问题，采用改进的无限插值理论‘111生成绕翼型的C型网格，该方法可以很灵活地 控制翼型表面附近的止交性和网格垂直方向网格尺度变化。 流动基本方程为雷诺平均的可压缩、非定常N—S方程，其积分形式为： 昙归·咖+P面ds=0 (1) ·本文得到国家自然科学基金10572120)资助 ．231． 第六届两北地区计算物理学术会议论文集 2008 其中： Q=[∥，∥，∥，P】。 P，”，1，，e分别为空气密度、墨y方向的速度分量和单位体积的总内能，元为线积分的法向单位向 量，V为积分域，aV为包含体积分域的边界，F为通量项，它包括无粘项R和粘性项F，两部 分： 户：戽+F。 (2) 为了在非定常情况下，采用全隐式双时间方法求解非定常N—S方程，伪时间采用LU--SGS 隐式方法以提高计算效率，离散后的N—S方程在时间和空间上均具有二阶精度。此外，为了更 好地计算低速流动，对空气动力学基本方程运用了预处理方法[121。 由于目前湍流转捩位置没有较好的确定方法，本文以流场为全湍流情况进行计算，湍流模型 采用考虑分离模拟能力的修止B—L代数湍流模型。为了检验数值方法和程序的可靠性，以 NACA0015翼型的过失速1F定常试验【13】情况为验证对象，翼型绕前缘四分