01
引言
刷式密封是一种广泛应用于航空发动机等透平机械的优良接触式动密封[]。近年来,随着透平机械逐渐向高参数方向发展,由刷式密封引起的泄漏损失越来越大,直接影响透平机械的工作效率。因此开展刷式密封的泄漏流动特性研究具有重要理论和实际意义。
国内外学者主要从理论分析与实验研究两方面对刷式密封的泄漏流动特性进行研究。在理论研究方面,国外Bayley和Long[4]首次采用多孔介质模型,在仅考虑粘性阻力的条件下数值计算了刷式密封的泄漏量;Chew等[]同时考虑了粘性阻力系数和惯性阻力系数的non-Darcian多孔介质模型,理论研究了刷式密封泄漏流动特性并与实验结果相对比。国内王之栎、丁水汀、黄学民等[]建立了多孔介质模型,理论分析了刷式密封泄漏流动特性的影响因素;邱波等[10]采用Ergun方程推导得出了刷丝束区域的粘性阻力系数表达式,并数值验证了刷式密封的泄漏量;张元桥等[]建立多孔介质传热模型,在考虑摩擦热效应的条件下,数值分析了不同工况参数对刷式密封泄漏量的影响;孙丹等[]建立了三维瞬态的流固耦合刷式密封求解模型,在考虑刷丝弯曲变形的条件下,分析了小压比条件下刷式密封泄漏特性的影响规律。
研究结果表明,考虑刷丝变形的瞬态流固耦合求解模型计算得出的泄漏量比多孔介质模型更精确。刷式密封刷丝在流体的作用下会产生变形,刷丝变形进一步影响流体分布,刷丝与流体之间的相互作用是典型的瞬态双向流固耦合问题[]。而传统的刷式密封流固耦合方法在计算时由于网格畸变,导致负体积问题,因此难以计算大压比条件下刷丝接触变形对刷式密封泄漏特性的影响。
在泄漏特性实验研究方面。国外学者Ferguson[18]首先建立了密封实验台,通过实验验证了刷式密封的密封性能远优于传统迷宫密封;Bayley等[19]通过建立静态刷式密封实验台,实验结果表明,随着压比的增大刷式密封泄漏量不断增大,同时测出了刷丝束区域的压力分布;Turner等[20]对Bayley和Long的刷式密封实验台进行改进,实验测量了两种间隙条件下的刷式密封的泄漏量和压力分布;Chupp等[21]进行刷式密封泄漏特性实验,验证了刷式密封泄漏量远低于迷宫密封,进出口压比的加载过程也会影响刷式密封的泄漏量。国内孙晓萍[22]进行刷式密封磨损与泄漏特性实验,发现随着磨损的增大,刷丝与转子之间的泄漏量会增大;曹广州等[23]实验研究了低转速条件下,刷式密封的泄漏量不断减少,后趋于稳定;胡广阳等[24]进行刷式密封磨损特性与泄漏特性实验。
结果表明,转速对泄漏量影响较少,刷式密封的密封性能远优于迷宫密封;杜春华等[25]研究了刷式密封在间隙、过渡和过盈三种条件下磨损特性对泄漏量的影响,在实验的初始阶段,过盈配合时磨损对泄漏影响最大,过盈配合时次之,间隙配合时最小。综上,现有文献对刷式密封泄漏特性实验研究方法比较单一,鲜有文献考虑刷丝接触变形对泄漏特性影响的论文公开发表。
本文提出基于ALE流固耦合方法的刷式密封三维瞬态双向流固耦合求解模型。考虑了刷丝接触变形对刷式密封泄漏特性的影响,在验证数值模型准确性的基础上,分析了刷式密封的流场分布特性,研究了压比、刷丝间隙、刷丝直径、刷丝排数等参数对刷式密封泄漏的影响规律,在验证ALE流固耦合方法准确性基础上,研究了不同结构参数对刷式密封泄漏特性的影响规律。
02
基于ALE流固耦合的刷式密封理论研究
2.1 刷式密封流固耦合特性分析
刷式密封三维模型如图1所示,刷丝束由柔软纤细的刷丝交错层叠构成,并沿着转子旋转方向呈一定角度排列,使刷丝对转子的瞬间径向变形或偏心运动具有更好的自适应性,保持良好的密封性能。研究表明,刷丝在气流力的作用下产生变形,刷丝变形进一步影响流场分布,刷丝与流体之间的相互作用是典型的流固耦合问题。当气流力较大时,刷丝的弯曲变形量也随之增大,且刷丝与刷丝之间产生接触碰撞,此时刷丝的运动状态对刷式密封的泄漏量将产生较大的影响。
Fig. 1 Schematic diagram of brush seal structure
2.2 基于ALE方法的刷式密封流体动力学模型
刷式密封流体采用ALE算法,ALE算法与Euler算法、Lagrange算法相比的优点在于其坐标系不固定,在空间也不依附于物体节点,网格可以做任意的运动,该方法可以克服大变形问题难题,主要用于求解流固耦合问题[26]。ALE算法可以计算刷丝接触和大变形,物质质点A在t时刻的空间坐标可表示为
(1) |
式中A表示物质坐标,也称为拉格朗日坐标,t表示时刻。
网格运动可描述为
(2) |
式中为参考坐标系,也称为ALE坐标,它是指网格点的标记。
由式(1)和(2)可知,质点A在t时刻所对应ALE坐标为
(3) |
质点A在参考坐标系中的运动速度(即质点参考速度)为
(4) |
网格点的位移、速度和加速度分别为
(5) |
(6) |
(7) |
利用链锁法则,可将物质速度表示为
(8) |
式中右端第一项为,右端第二项中的是映射的雅可比矩阵,它将质点的参考速度