迷宫密封不稳定力三维CFD数值分析方法研究.pdfVIP

迷宫密封不稳定力三维CFD 数值分析方法研究 刘晓锋，陆颂元 东南大学振动控制与信息系统研究所，南京（210096 ） E-mail ：yourado@163.com 摘 要：迷宫密封中的泄漏流体与转子相互作用产生的气动力能诱发轴系失稳，密封转子动 特性系数描述了这一不稳定力的重要特性。论文建立了涡动转子－迷宫密封三维CFD 数学 模型，并以压缩机入口迷宫密封为计算实例，利用通用CFD 软件Fluent 求解三维CFD 模型， 获得转子表面气动力，利用涡动转子－迷宫密封动力学模型预测了迷宫密封动特性系数。计 算结果与通用 CFD 软件 TASCflow 和基于整体流动理论的迷宫密封动特性计算程序 DYNLAB 的计算结果进行比较，表明三维CFD 方法对于漏气量和动特性的预测精度高于整 体流动理论。此外，论文进一步计算分析了密封入口预旋速度和间隙变化以及轴向锥形间隙 对迷宫密封稳定性的影响，为消除汽流激振提供理论依据。 关键词：迷宫密封；计算流体动力学；整体流动理论；转子动特性系数 中图分类号：TK26 TH113.1 文章标识码：A 1 引言 旋转机械转子系统中的非接触密封在防止流体泄漏的同时，还会产生重要的流体激振 力，影响转子的振动稳定性。与动压滑动轴承的油膜类似，密封中流体膜 （汽膜、水膜或油 膜）动力特性中的交叉刚度k 是促使转子做非同步低频涡动的激振力来源；同时，动特性中 的直接阻尼C 可以生成对这种低频涡动的抑制力。多年来，研究人员对密封的动特性、流 体激振力和轴系失稳，从理论分析、数值计算、试验测定等各方面做了深入研究；研究对象 以迷宫密封为主，近年延伸到新型密封，如蜂窝密封中；传统的理论研究分析方法采用以整 体流动理论为基础的控制容积法，近年开始了采用计算流体动力学（Computational Fluid Dynamics，CFD ）数值计算和分析的研究。 1987 年，Scharrer和Childs [1]将回流速度比例系数作常数处理，计入了腔内回流速度对 壁面剪应力的影响，建立了较完善的整体流两控制容积模型。至今，两容积模型一直是西方 国家对迷宫密封动特性分析研究的主要方法，虽然在工程实际中得到广泛应用，但它最大的 缺陷在于不能正确模拟密封中的回流，因此，整体流动模型预测泄漏量有较高的精度，但预 测动特性的精度不理想。 Dietzen和Nordman[2] 1987 年首先提出了二维CFD摄动模型，用有限差分法求解不可压缩 流体N-S方程，获得了环状密封的动特性系数；其后，Tam[3] 、Rhode[4]和Athevale[5]进一步 采用三维CFD方法分析了环状密封和迷宫密封动特性； 2003 年，Jeffrey Moore[6]利用商业 CFD程序SCISEAL成功计算出了静子齿的 8 齿迷宫密封动特性；同年，Kirk和Toshio[7]等应 用另一通用CFD程序TASCflow计算分析了压缩机和汽轮机的迷宫密封动特性。 国内 1991 年，鲁周勋、谢有柏[8]应用二维有限差分法和摄动法计算了环状紊流密封的 动特性，用流线分离法划分两控制容积的界限计算了迷宫密封的动特性；2003 年，李松涛、 黄典贵等[9]利用自编程序，对迷宫密封流场进行三维非定常数值求解，得到密封汽流激振力 和动特性系数。 本文的研究是建立涡动转子－迷宫密封三维CFD 数学模型；应用当前国内流通的CFD 软件 Fluent 计算压缩机叶轮入口迷宫密封流场以及转子表面气动力；根据涡动转子－迷宫 密封动力学模型计算密封动特性系数；最后，本文计算分析了密封入口预旋速度和间隙变化 以及轴向锥形间隙对迷宫密封稳定性的影响。 -1- 2 涡动转子－迷宫密封动力学模型 准确的密封动特性系数是转子－轴承系统稳定性分析的基本计算条件。如果转子在静子 中的某一偏心位置做小轨迹涡动，可用密封气动力与转子位移、速度的线性模型定