基于相似原理的环瓣式石墨密封泄漏流动特性求解模型.pdf

年 月第 卷 第 期润滑与密封 ：文献引用：薛庆增，胡振圆，孙丹，等基于相似原理的环瓣式石墨密封泄漏流动特性求解模型润滑与密封，（）：，（）：基金项目：国家自然科学基金项目（）；辽宁省兴辽英才计划青年拔尖人才项目（）收稿日期：；修回日期：作者简介：薛庆增（），男，硕士，工程师。：。通信作者：孙丹（），男，博士，教授，主要研究方向为透平机械先进密封技术。：。基于相似原理的环瓣式石墨密封泄漏流动特性求解模型薛庆增 胡振圆 孙 丹 王 双 赵 欢（海装沈阳局驻沈阳地区某军事代表室 辽宁沈阳；西安交通大学航天航空学院，复杂服役环境重大装备结构强度与寿命全国重点实验室 陕西西安；沈阳航空航天大学航空发动机学院，辽宁省航空推进系统先进测试技术重点实验室 辽宁沈阳）摘要：环瓣式石墨密封因其泄漏通道尺寸微小，导致其建模、网格划分以及计算困难。基于相似原理方法建立环瓣式石墨密封泄漏通道求解模型，采用方程分析法推导环瓣式石墨密封泄漏通道内流体流动相似准则，获得遵循几何相似和力学相似的映射模型，并采用建立的泄漏通道映射模型分析环瓣式石墨密封的泄漏流动特性，并与实际模型计算结果进行比较。研究结果表明：泄漏通道内气体流动相似性可综合采用弗劳德、欧拉、雷诺相似准则表征；映射前后模型相同结构位置处流体压力、速度分布具有较好的一致性。通过映射模型求解的泄漏量与实际模型求解的泄漏量相对误差在误差允许范围内，验证了推导的泄漏通道流体流动相似准则和映射方法的可靠性，为研究环瓣式石墨密封微小泄漏通道泄漏流动特性提供新方法。关键词：环瓣式石墨密封；相似原理；泄漏流动；方程分析法中图分类号：（，；，；，）：，：；环瓣式石墨密封因其优良密封性能，在航空发动机等透平机械中得以广泛应用。其中石墨密封环在工作时因流体动压效应而被径向推开，其与跑道间形成极其微小的泄漏通道，形成非接触式密封。然而因其泄漏通道尺寸微小导致建模及网格划分困难，开展数值计算工作较为困难，因此研究环瓣式石墨密封泄漏流动特性求解方法具有重要意义。目前，国内外学者就环瓣式石墨密封的泄漏流动特性做了深入的研究。理论研究方面，采用热力学原理分析了理想密封模型内流体流动特性并提出密封泄漏量计算公式；和 在 公式的基础上进行修正，进一步研究了石墨密封的 泄 漏 流 动 特 性。数 值 计 算 方 面，、采用数值分析方法研究了转速、几何参数和工况条件对密封泄漏量的影响；通过数值计算方法，研究了石墨环内表面有无槽结构对石墨密封泄漏量和流场压力分布的影响；运睿德等通过建立泄漏量计算模型，结合 软件分析了不同弹簧力下环瓣式石墨密封的泄漏流动特性。此外，其他学者也采用数值计算方法研究了石墨密封泄漏流动特性。在关于密封泄漏流动特性方面，也有少数学者采用其他先进理论来研究密封泄漏流动特性，其中就包括相似原理理论。朱毅征和何枫基于相似原理，通过实验回归分析法归纳出流量系数、迷宫函数曲线与迷宫密封泄漏率之间的关系；崔海龙等基于相似原理建立有限元模型并通过 计算得到微小气膜间隙内的压力场；刘兴旺等采用相似原理中的 定理来研究径向迷宫密封泄漏量的影响因素，并预测了实际模型中泄漏量的变化趋势。综上所述，国内外专家学者对石墨密封泄漏流动特性开展了大量研究，但在研究泄漏流动特性时建立的泄漏通道模型与实际工况存在一定偏差，且石墨密封中关于相似原理的应用较少。本文作者将相似原理应用到环瓣式石墨密封泄漏流动特性研究中，通过相似原理中的方程分析法建立环瓣式石墨密封泄漏通道内气体流动相似准则，基于相似准则建立石墨密封泄漏流动特性求解模型，为环瓣式石墨密封泄漏流动特性研究提供重要基础。环瓣式石墨密封泄漏流动相似原理理论分析 环瓣式石墨密封工作原理环瓣式石墨密封组件由石墨环、密封座、拉伸弹簧、轴向弹簧、防转销、挡板等部件组成。石墨密封沿周向分成若干段相等弧长的石墨环瓣，各石墨环瓣之间通过搭接头和拉伸弹簧连接，各搭接头之间有一定间隙，用来补偿密封环的磨损和制造误差。石墨环瓣内侧开有横槽、环槽和浅槽。石墨环在周向弹簧的作用下紧箍在密封跑道上，形成主密封面，在轴向弹簧的作用下紧压在密封座上，形成次密封面。图 所示为环瓣式石墨密封的整环结构示意图。带有浅槽结构的环瓣式石墨密封在工作时会开启，石墨环内表面与转子跑道之间形成尺寸微小的主泄漏通道，形成非接触式密封。图 所示为环瓣式石墨密封的整体结构示意图。图 环瓣式石墨密封装置示意 ：（）；（）润滑与密封第 卷图 环瓣式石墨密封泄漏通道示意 泄漏流动控制方程物理模型的相似必须以物理本质相同为基本原则，泄漏通道内气体流动的控制方程反映了泄漏通道内气体流动的物理本质，故泄漏通道内气体流动的几何和力学相似应立足于相同的物理流动控制方程。文中将泄漏通道内气体流动简化为不可压缩黏性流体的非定常等温流动过程，暂时忽略温度变化对环瓣式石墨密封泄漏通道内流体流动的影响。泄漏通道内气体流动状态的连续性方程和动量守恒方程可分别用式（）和式（）表示。（）（）式中：、为流体速度矢量在、个方向的速度分量；为流体密度；为流体压力；为流体动力黏度；、为质量力。泄漏流体流动状态相似准则泄漏通道内气体流动相似分析的一个原则是进气压力、封严间隙、转速等工况参数组成的相似准则应相等，想要映射后的泄漏通道内气体流动反映实际泄漏通道内气体流动，则需要二者之间的相似准则数应相等。采用方程分析法从式（）和式（）中推导出相似准则数，使映射后的泄漏通道内气体流动和实际泄漏通道内气体流动相似准则数相等而形成两者相似。以上标“”表示映射后的泄漏通道内流体流动，则映射泄漏通道内流体流动与实际泄漏通道内气体流动几何及物理量之间的对照关系，个相似流动现象的速度、压力、密度、质量力和坐标的对应关系如式（）所示。，（）式中：、和 为流场直角坐标系空间坐标；、和 分别为直角坐标下、和 方向上流体的速度分量；为相似倍数，不同几何及物理量的相似倍数以下标区分。利用式（）对实际泄漏通道内气体流动控制方程进行相似变换，通过方程中的几何量及物理量的对应转化，形成映射泄漏通道内气体流动控制方程，以在满足几何和力学相似的条件下基于映射泄漏通道进行气体流动问题分析。映射后泄漏通道内气体流动控制方程可用式（）和（）表示：（）（）比较相似转化前后的连续方程式（）和式（）及动量方程式（）和式（）。根据量纲理论，相似倍数应相等，即：（）年第 期薛庆增等：基于相似原理的环瓣式石墨密封泄漏流动特性求解模型 常数（）从而可得：（）把式（）代入上面的相似指标关系式中，整理可得到封严间隙气体流动的弗劳德准则、欧拉准则、雷诺准则，如下：记 （）记 （）记 （）式（）（）给出的各项相似准则数中，依据环瓣式石墨密封结构的应用场景，特征速度 分别用进气特征速度、转子表面线速度 表征。环瓣式石墨密封泄漏通道内的气体流动过程类似于同心缝隙管道流动，特征尺寸 可由封严间隙 表示为 （）泄漏通道内进气特征速度 和转子表面线速度定义为（）式中：为封严压差；为封严气体密度；为转子转速；为转子半径。利用式（）得到表征进气状态下的弗劳德数和转子转动状态下的弗劳德数，分别为（）式中：为封严间隙；为重力加速度。利用式（）得到表征进气状态下的欧拉数 和转子转动状态下的欧拉数，分别为（）式中：为进气压力。利用式（）得到表征进气状态下的雷诺数 和转子转动状态下的雷诺数，分别为（）环瓣式石墨密封结构与工况 实际泄漏通道结构与工况文中数值求解模型参考文献的环瓣式石墨密封，以整环中单环瓣为研究对象，建立密封泄漏流动特性计算模型。图 所示为文献中环瓣式石墨环瓣单环瓣的结构示意图，图 所示为依据文献 结构参数建立的泄漏通道模型及泄漏通道部分尺寸示意图，表 及表 分别给出了泄漏通道结构参数及工况参数。图 石墨环瓣结构示意 图 泄漏通道模型及尺寸示意 表 泄漏通道结构参数 密封半径 封严间隙 主密封唇宽度 主泄漏通道长度 环槽宽度 环槽深度 横槽长度 横槽宽度 浅槽长度 浅槽宽度 浅槽厚度 润滑与密封第 卷表 泄漏通道工况参数 工况编号高压侧压力 低压侧压力 转速（）映射泄漏通道结构与工况 泄漏通道结构参数进行映射模型泄漏量分析需要确定映射后泄漏通道的结构和工况条件。在分析过程中应保证实际模型与映射模型的几何相似。可以根据研究内容来确定映射模型与实际模型的比例因子，这样可以保证映射模型是经过实际模型缩小或放大之后得到的，因此可以保证在几何层面映射模型与实际模型的相似性以及结果的准确性。文中根据所研究内容，选取了 个映射模型与实际模型之间的相似倍数，分别为 、。表 所示的是实际模型结构按照上文中所描述的相似方法计算出的相对应的映射模型结构参数。表 映射泄漏通道结构参数 映射模型编号相似倍数密封半径 封严间隙 主密封唇宽度 主泄漏通道长度 环槽宽度 环槽深度 横槽长度 横槽宽度 浅槽长度 浅槽宽度 浅槽厚度 泄漏通道工况参数依据实际模型应与映射模型相似准则数相等的原理，按照式（）（）可以计算出映射模型的进口压力、轴转速、重力加速度 个工况参数。表 所示的是实际模型按照上文中所描述的相似方法计算出的相对应映射模型工况参数。表 映射泄漏通道工况参数 映射模型编号相似倍数工况编号高压侧压力 轴转速（）重力加速度（）环瓣式石墨密封数值建模 数值求解模型在 软件平台利用 模块分别构造实际模型和映射模型（周向 ，实际模型为 环瓣，选取其中之一；轴向位置以主泄漏通道长度 做量纲一化处理并以 表示）。图 所示为实际模型示意图，其中流场边界条件及分析所需圆柱坐标系如图所示，映射后的 种模型也采用同样的处理方式。图 实际模型示意 高压侧和低压侧均施加压力边界条件，转子表面施加运动壁面条件。在 平台中设置迭代步数为 步，收敛标准为，湍流模型为 模型，封严气体物理性质同理想气体一致。网格划分及无关性验证文中采用多区域网格划分方法对实际模型进行网格划分，共划分出 种网格数量模型，在高压侧压力 ，低压侧压力 ，转速 工况条件下对实际模型的泄漏量以及压力分布展开分析。模型的网格数量分别为 、，图 所示为泄漏量随网格数量的变化。可以看出网格数量对计算结果产生的误差不大，综合考虑计算准确度、计算时间及计算资源的因素，文中将模型划分为 个网格。年第 期薛庆增等：基于相似原理的环瓣式石墨密封泄漏流动特性求解模型 图 泄漏量随网格数量的变化 模型准确性验证为了验证文中所建立的实际模型数值计算的准确性，基于文献 给出的泄漏通道结构和工况条件进行了数值仿真计算，并与文献 中的数据进行了比较。表 给出了文中实际模型与文献 中参考模型密封泄漏量对比。文中实际模型泄漏量与文献 中泄漏量结果平均误差为 ，两者结果较吻合，验证了文中所建立求解模型的准确性。表 石墨密封泄漏量对比 进口压力 泄漏量 （）文中文献相对误差 环瓣式石墨密封数值计算分析 流场特性分析 流场压力分布分析在实际模型和 种映射模型周向 、轴向 （）共四处位置提取泄漏通道内气体流动压力值。图 给出的是在工况 条件下映射模型的流场压力分布云图，由于流体动压效应的作用在浅槽位置处压力明显较高。图 表示的是在 种工况条件下，种映射模型与实际气膜模型压力值对比。图 （）和图 （）均表明映射泄漏通道内气体流动与实际泄漏通道内气体流动压力分布变化趋势基本一致。在工况 条件下，实际模型中最大压力为 ，进口压力为 ，压力增高 ；映射模型中最大压力为 ，进口压力为 ，压力增高 ；映射模型中最大压力为 ，进口压力为 ，压力增高 ；映射模型中最大压力为 ，进口压力为 ，压力增高 。种映射模型的压力平均增量为 。以上数据表明映射后的各个模型泄漏通道内压力增量有较好的一致性。图 工况 下映射模型压力分布云图 图 不同工况下各模型泄漏通道内压力对比 ：（）；（）流场速度分布分析在实际模型和 种映射模型周向 、轴向 共四处位置提取泄漏通道内气体流动轴向速度值。图 给出的是在工况 条件下映射模型的流场速度分布云图，图 表示的是在 种工况条件下，种映射模型与实际模型轴向速度值对比。图 （）润滑与密封第 卷和图 （）均表明映射泄漏通道内气体流动与实际泄漏通道内气体流动轴向速度变化趋势基本一致。在工况 条件下，实际模型中最大速度为 ，最小速度为 ，速度增大；映射模型中最大速度为 ，最小速度为 ，速度增大 ；映射模型中最大速度为 ，最小速度为 ，速度增大 ；映射模型 中最大速度为 ，最小速度为 ，速度增大 。种映射模型速度的平均增量为 。以上数据表明映射后的各个模型泄漏通道内速度增量有较好的一致性。图 工况 下映射模型轴向速度分布云图 图 不同工况下各模型泄漏通道内轴向速度对比 ：（）；（）泄漏特性分析图 分别给出了在工况 和工况 条件下实际模型泄漏量与映射模型、数值计算泄漏量对比。从图 （）中可以看出，在工况 条件下，映射模型和的泄漏量与实际模型泄漏量相对误差值较小，分别为 和 ，但映射模型相对误差较大，为 ，种映射模型的平均相对误差为 。从图 （）中可以看出，在工况 条件下，映射模型和的泄漏量与实际模型泄漏量相对误差值较小，分别为 和 ，但映射模型相对误差较大，为 ，种映射模型的平均相对误差为 。产生误差的原因为：映射后求解模型节点及网格疏密程度发生变化，网格划分质量对数值计算结果产生一定影响；文中简化了求解模型的建立过程，将流体流动过程设定为等温过程，暂时忽略了实际情况中温度变化对流体流动的影响，故产生上述误差。图 不同工况下各模型泄漏量对比 ：（）；（）结论（）环瓣式石墨密封泄漏通道内气体流动的相似性可采用弗劳德相似准则、欧拉相似准则、雷诺相似准则 个相似特征数来表征。年第 期薛庆增等：基于相似原理的环瓣式石墨密封泄漏流动特性求解模型 （）基于相似原理建立的映射模型能较好地反映实际模型泄漏通道内流体流动状态，种模型流动状态相吻合。（）映射模型求解的泄漏量与实际模型求解的泄漏量平均相对误差为 ，证明所建立的相似准则可准确计算环瓣式石墨密封泄漏量，为石墨密封泄漏流动特性分析提供重要理论基础。参考文献 张家远，王利恒，顾健石墨密封基础理论及工程应用技术研究进展装备制造技术，（）：，（）：徐凯松石墨密封应用发展现状的进展研究炭素，（）：，（）：，（）：，（）：，（）：，（）：，运睿德，陈志英，刘勇，等周向弹簧力分布对圆周密封装置密封性能的影响推进技术，（）：，（）：，（）：，（）：胡廷勋，周坤，王晓燕，等浮环密封泄漏特性数值计算与试验航空动力学报，（）：，（）：王衍，王达，曹彬，等基于逾渗理论的接触式机械密封泄漏特性分析润滑与密封，（）：，（）：朱毅征，何枫迷宫密封的理论与实验研究航空学报，（）：，（）：崔海龙，岳晓斌，雷大江，等基于 的气体静压轴承径向特性分析润滑与密封，（）：，（）：刘兴旺，张官正，李超，等基于相似理论的涡旋压缩机径向迷宫密封泄漏量研究化工机械，（）：，（）：张堃元，金志光流体动力学北京：科学出版社，：；，：，（）：润滑与密封第 卷