跨音速风扇叶尖间隙对颤振特性的影响研究.docxVIP

跨音速风扇叶尖间隙对颤振特性的影响研究 杨俊，郑赟 （北京航空航天大学，能源与动力工程学院，北京，100191） 摘 要：为了研究跨声速风扇叶片叶尖间隙大小对叶片颤振特性的影响，以NASA rotor67转子为例建立了单排全环流动模型。求解定常和非定常流场使用使用雷诺平均Navier-Stokes方程（RANS），结构动力学方程的求解使用模态叠加法,利用时间推进法模拟颤振的整个非定常过程直到达到收敛条件,对得到的模态位移曲线进行拟合得到对应状态的气动阻尼值。定常计算和颤振计算使用自行编制的HGFS计算程序计算，间隙范围从0.15%span-1.26%span。数值计算结果表明，叶尖间隙的大小与气动阻尼为非线性关系，存在一个中间叶尖间隙使气动阻尼最大，气弹稳定性最优。 关 键 词：叶尖间隙；泄漏涡；颤振；气动弹性稳定性 中图分类号：V231.3 文献标识码：A 文章编号： 引言 随着现代航空发动机的发展，人们对航空发动机的性能提出了更高的要求。一方面要求发动机具有高推重比、重量轻的特点，使发动机得结构更加紧凑；另一方面要求发动机具有高可靠性、长寿命以及低耗油率等特点，导致流体诱发的叶片振动问题越来越突出和复杂。叶片颤振对应的叶轮机械流动诱发振动现象的力学分支被命名为叶轮机气动弹性力学[[] 周盛。叶轮机气动弹性力学引论[M]。 国防工业出版社，1989。 ][[] 任汝根，周盛。叶轮机械中的非定常流动[J]。燃气涡轮试验与研究，1991，4（4）：1-24。 ]。航空发动机转子叶片颤振是一种气动力引起的自激不稳定的振动。由于流体域作用在固体域上做正功，气动力做功产生的能量转子叶片不能吸收，颤振就会发生。而在转子叶尖附近，以顶隙泄漏流、通道激波、叶片边界层及其相互影响为特征，形成了复杂的非定常三维流动。颤振的发生是一种不良振动，这类振动导致叶片高周疲劳失效，与一系列灾难性的事故都有所关联[[] 胡骏, 吴铁鹰, 曹人靖。航空叶片机原理[M]。国防工业出版社, 2006.3。 ]。近几十年来，人们对颤振现象进行了深入的实验和计算研究。 1数值计算方法 流固耦合法是现阶段对于颤振分析应用较多的一种方法。它考虑了叶片结构变形和绕流运动之间的相互作用，相对于纯流体域和固体域的模拟，更加符合物理实际。在流体域的求解中，本文采用自行编制的程序HGAE[[] 郑赟, 杨慧。跨音速风扇全环叶片颤振特性的流固耦合分析[J]。北京航空航天大学学报, 2013, 5: 011。 ] [[] 郑赟. 基于非结构网格的气动弹性数值方法研究[J]. 航空动力学报, 2009, 24(9). ]，其计算器HGFS 使用雷诺平均 Navier-Stokes 方程（RANS）的方法来求解定常和非定常流场。 图  SEQ 图 \* ARABIC 1 流固耦合方法计算流程 如图1，在每一个物理时间步内，在流体域上求解流体控制方程来得到其作用于固体表面的压力，从而为结构计算提供此时刻的边界条件。在固体域则利用模态叠加法求解结构动力学方程，得到结构体的变形位移。将结构体变形后的位置作为流体域在下一时刻的边界条件，通过网格变形技术对流体域网格进行更新。在此基础上利用时间推进法模拟颤振的整个非定常过程直到达到收敛条件。通过对得到的模态位移曲线进行拟合可以得到对应状态的气动阻尼值。 1.1流体域的求解 在流体域的求解中，本文所用的程序为HGAE，其计算器HGFS 使用雷诺平均 Navier-Stokes 方程（RANS）的方法来求解定常和非定常流场。对于转静交界面，在定常计算时使用了掺混面模型，非定常计算时使用了滑移面模型。 非定常、可压缩的雷诺平均N-S方程在笛卡尔坐标系下的微分形式可写为： ?Q?t+?F-Fv?x+?G-Gv?y+?H-Hv?z=0 （1.1） 雷诺平均N-S方程保留了方程中的粘性项，只对湍流脉动做了时间平均，从而产生了雷诺应力项。为了对雷诺应力进行模拟，从而衍生出了各种湍流模型并且每种都有其使用的工况[[] Srinivasan G R, Ekaterinaris J A, McCroskey W J. Evaluation of Turbulence Models for Unsteady Flows of an Oscillating Airfoil[J]. Computers Fluids, 1995, V24: 833-861 ]。本为所选的湍流模型为计算效率和精度都较好的Spalart-Allmaras（S-A）模型。 如式（1.2），S-A模型是一个关于v的一方程模型。v与湍流粘性系数μT的关系为μT=ρvfυ1。 ?v?t+uj?v?xj=Cb11-ft2Ωv+Cb11-f