叶片颤振问题的流固耦合分析.docxVIP

叶片颤振问题的流固耦合分析 叶片机械的叶片振动是影响其可靠性的最严重问题之一。现代发动机的设计应提高压气机和风扇的压缩比，并使用轻新材料。然而，在现代压气机的前级和超级别空气扇中，有强烈的激发波对附板层的干扰，以及通过近速点激发波诱导的附板层的分离。叶片振动预测的经典方法，如单参数试验模型或变形激盘法和半体验法，通常采用线性化气动力模型，预测精度难以提高。考虑到近年来叶片振动引起的发动机故障，提出了一种更详细的叶片振动预测方法。这对提高发动机的安全性和可靠性具有极其重要的意义。 基于目前对流固耦合叶片振动物理机制的认识,气动弹性分析模型大致分为3种:能量法、半耦合方法与全耦合方法.目前广泛使用的气动弹性分析模型是由文献提出的能量法,这种方法认为颤振发作时的频率与结构的固有频率一致.基于文献的行波理论的全环叶片振动方法认为颤振发作时全环叶片以行波节径形式振动(以节径数或叶片间相位角表征),根据结构的循环对称性,取一个循环扇区并给定叶片的振动频率、振型和叶片间相位角进行振荡叶片的非定常流动分析,由一个振动周期内气动力对叶片做功得到气动阻尼;针对每个气动工况,需要遍历所有可能的叶片间相位角(叶片角相位角的数目等于该叶排内的叶片数M),以寻求气动弹性最不稳定叶片间相位角.这种方法的优点在于某特定叶片间相位角的颤振分析效率较高,但颤振分析一般要求对所有的行波模态的稳定性进行分析,所以当考虑所有叶片间相位角和多个模态(N代表模态数)的颤振分析时,该方法需要进行N×M个叶片振动的非定常流动分析,计算量大;而且该方法中需要预先给定叶片振动方式,且不能考虑流体与固体之间的相互作用. 流固耦合法是交替求解流体和结构体运动方程组,流体域和固体域在每一个时间步通过相互传递边界条件建立联系,是3种气动弹性模型中最灵活的一种.针对不同的工程需求,可以耦合不同复杂程度的流体模型和叶片振动模型.文献将二维非定常流N-S方程与简化为特征截面上二维振动形式的动力学方程组耦合,研究表明解除流固耦合约束的叶片表面非定常力作功与流固耦合时相比有显著差异.文献求解三维非定常N-S方程及纯扭转振动的简化方程,分析攻角和背压对三维透平叶片失速颤振的影响.流固耦合方法的主要缺点是三维非定常流场和结构动力学方程求解都使用时间推进法,计算时间较长,对计算资源的需求很大.文献使用模态叠加法解耦结构动力学方程大大减少了流固耦合分析的时间. 本文在自行开发的三维流固耦合程序基础上发展了在叶片振动的气动弹性计算中同时考虑多个叶片振动模态的气动弹性模型,采用该全环气动弹性模型,叶片振动无需预设叶片间相位角和振动模式等参数,通过叶片在非定常流场中振幅随时间的变化历程,判断叶片的气动弹性稳定性.本文采用发展的流固耦合的分析方法,考察了NASA Rotor67跨音速风扇转子在不同工作点时前三阶固有模态下所有行波节径形式的气动阻尼的变化,由此确定气动弹性最不稳定状态对应的叶片振动节径形式. 1 气动弹性模型 1.1 控制体t的状态向量 求解动边界和动态变形网格下的流固耦合问题采用了守恒型积分形式的非定常可压缩雷诺平均Navier-Stokes方程,其流体运动方程为 ??t∫Ω(t)qdΩ+∮?Ω(t)(F?w?q)?n?dA=∫Ω(t)SdΩ(1)??t∫Ω(t)qdΩ+∮?Ω(t)(F-w?q)?n^dA=∫Ω(t)SdΩ(1) 其中,Ω,?Ω和n?n^是某一物理时刻的控制体的体积、表面面积和表面外法向;q是守恒型流体状态向量;F代表了对流和耗散通量项;w是运动控制体表面的网格速度;S为坐标系运动等带来的源项.湍流流动中采用了Sparlart-Allmaras的湍流模型. 1.2 叶片模态的推导 使用广义坐标解耦叶片动力学方程,得到N个线性无关的常微分方程组: Mgζ??+Cgζ˙+Kgζ=Qg(2)Μgζ??+Cgζ˙+Κgζ=Qg(2) 其中,ζ是模态位移向量;ζ˙ζ˙是模态速度向量;ζ??ζ??是模态加速度向量;使用质量归一化质量矩阵、阻尼矩阵和刚度矩阵,得到Mg,Cg和Kg;Qg为广义气动力.详细推导过程参见文献. 广义坐标ζ与物理位移x之间的关系如式(3).Φ为质量归一化的系统模态.首先从有限元程序中获得单叶片的N阶模态,然后,使用式(4)组集第i阶叶片全环模态φi. ζ=Φx(3)φi=Xˉˉˉk=0+∑k=1k[XˉˉˉkAcos(j?1)kα+XˉˉˉkBsin(j?1)kα]0kN/2+(?1)j?1Xk=N2(4)ζ=Φx(3)φi=Xˉk=0+∑k=1k[XˉkAcos(j-1)kα+XˉkBsin(j-1)kα]0kΝ/2+(-1)j-1Xk=Ν2(4) 其中,包括0节径整体对称模态;前行波后行波成对出现的0kN/2节径整体模态;当叶片数为偶数时,N/2节径反对称整体模