基于计算流体动力学的迷宫密封流场计算.docxVIP

基于计算流体动力学的迷宫密封流场计算 密封动力学特性 由于结构简单、无摩擦、油耗低、使用寿命长、维护方便等优点，广泛应用于汽轮机、离心分离器、鼓风机、涡轮机等旋转机械的轴端和级间密封。随着科学技术的不断进步,为了满足封严特性的要求,密封间隙的尺寸不断减小;另一方面,为了符合某些设计要求,转子的柔性不断增大,这些均会导致转子—密封系统自激振动易于发生。研究表明,压气机中的气流激振力主要来自于其中的迷宫密封;而对于汽轮发电机组而言,蒸汽诱发的转子不稳定振动有时会成为限制机组出力的重要因素。 国内外研究人员对密封的动力学特性从理论分析、数值计算和试验测定等各方面做了深入研究。在迷宫密封动力学特性的预测方面,传统的理论分析方法主要是采用以整体流(bulk flow)理论为基础的控制容积法。1980年,日本的Iwatsubo T首先应用整体流理论建立密封的单控制体方程,分析转子泄漏量和转子动力学特性。Scharrer J K将回流速度比例系数作常数处理,计入空腔回流速度对壁面剪应力的影响,建立比较完善的整体流两控制容积模型。国内鲁周勋和谢友柏等在1991年用流线分离法划分两控制容积的界限,计算迷宫密封动力学特性。两控制容积方法虽然在预测密封泄漏量上具有较高的精度,但在预测转子动力学特性方面不太理想。随着计算机技术和计算流体动力学的发展,流体计算动力学方法在分析研究迷宫密封动力学特性方面得到相应的应用。Dietzed F J和Nordman R首先提出二维CFD(computational fluid dynamics)摄动模型,用有限差分法求解不可压缩流体的N—S方程,得到环状密封的动特性系数;其后,Tam L T等采用三维CFD方法分析环状密封和迷宫密封的动特性;Jeffery Moore J利用商用软件Sciseal计算8个静子齿的迷宫密封动力学特性;Toshio H等820-826利用商用软件CFX-TASCflow计算和分析空气压缩机和汽轮机的迷宫密封动力学特性。近年来,国内很多专家学者也都应用流体计算软件对迷宫密封流场进行深入的研究,如孙婷梅、叶建槐、李军、王炜哲。 本文根据文献820-822中的模型,应用Gambit软件建立转子—密封三维CFD模型,并应用CFD计算软件Fluent进行密封流场计算。Fluent计算结果与文献822-826中CFX(computational fluid dynamiX)的计算结果进行比较。文献中也对相关模型进行了分析。本文为了更好地研究迷宫密封动力学特性,进一步讨论出口与进口压比对密封的泄漏量和转子动力学系数的影响,研究表明气流的轴向速度在跨音速前后有明显的不同。 1 密封流场的定常求解方法 当密封中的转子在偏心位置运动时,根据文献,在直角坐标系(x,y)下,气流激振力与密封转子位移(x,y)、速度(˙xx˙,˙yy˙)和加速度(¨xx¨,¨yy¨)关系可表示为 -[FxFy]=[ΚxxΚxyΚyxΚyy][xy]+[DxxDxyDyxDyy][˙x˙y]+[ΜxxΜxyΜyxΜyy][¨x¨y](1)?[FxFy]=[KxxKyxKxyKyy][xy]+[DxxDyxDxyDyy][x˙y˙]+[MxxMyxMxyMyy][x¨y¨](1) 式中,Fx、Fy为气流激振力在直角坐标中x、y方向上的力,Kxx、Kyy分别表示气体密封力的主刚度,Dxx、Dyy分别表示气体密封力的主阻尼,Mxx、Myy分别表示气体密封力的主惯性效应,Kxy、Kyx分别表示气流激振力的交叉耦合刚度;Dxy、Dyx分别表示气流激振力的交叉耦合阻尼,Mxy、Myx分别表示气流激振力的交叉耦合惯性效应。当密封转子在静子中作小幅振动时,气流激振力可以表示成比较简单的斜对称矩阵形式167-168 -[FxFy]=[Κˉk-ˉkΚ][xy]+[Dd-dD][˙x˙y]+[ΜmmΜ][¨x¨y](2)?[FxFy]=[K?kˉkˉK][xy]+[D?ddD][x˙y˙]+[MmmM][x¨y¨](2) 主刚度K反映气流对转子主刚度项的影响,可以改变系统的临界转速;交叉刚度项k是产生切向力的主要因素,即促使转子作非同步低频涡动的主要因素,当k增大时,转子的稳定性下降;主阻尼D反映流场对低频涡动的抑制力大小,一般增加主阻尼会让转子稳定;交叉阻尼项d对系统影响较小;质量M、m分别表示气流的主惯性效应和交叉惯性效应,有的文献中认为其数值较小可以忽略不计822-824。 当转子在密封中的转速为ω、涡动速度为Ω时,在固定坐标系(X、Y)下,转子轴心的运动轨迹可近似认为是以e为半径、以密封腔室几何中心为圆心的圆,如图1a所示。如果在固定坐标系下求解,流场是随时间变化的,则需要采用动网格技术进行非定常求解分析密封流场,这样增加了模型设置的难度