1、 年 月第 卷 第 期润滑与密封 ：文献引用：陈阳，张功学，吴垚多叶动压气体滑动轴承静态特性的有限差分算法润滑与密封，（）：，（）：基金项目：陕西省 年自然科学基础研究计划一般项目（）；陕西科技大学引进人才博士启动基金项目（）；西安科学技术局高校院所人才服务企业项目（）收稿日期：；修回日期：作者简介：陈阳（），男，硕士研究生，研究方向为气体滑动轴承动压润滑理论。：。通信作者：吴垚（），男，博士，讲师，研究方向为流体润滑理论和轴承转子系统动力学。：。多叶动压气体滑动轴承静态特性的有限差分算法陈 阳 张功学 吴 垚（陕西科技大学机电工程学院 陕西西安）摘要：多叶动压气体滑动轴承因其结构简单、摩擦阻
2、力低、旋转精度高和无环境污染等优点，在高速离心分离机、空气压缩机和透平膨胀机等旋转机械中应用广泛。为探究多叶动压气体滑动轴承的静态性能，通过数学变换将三叶动压轴承的气体润滑 方程转化为标准偏微分方程形式，利用有限差分法和超松弛迭代法进行数值求解，研究气膜厚度和气膜压力分布、承载力、摩擦因数和质量流量等静态性能，随偏心率、预负荷系数、轴承数、长径比及瓦块分布位置的变化规律。结果表明：三叶轴承的承载力和轴颈表面摩擦因数随偏心率和长径比的增加而增加，而偏位角和质量流量随偏心率和预负荷系数的增加则呈现出相反的变化趋势；随着轴承数和预负荷系数的增大，承载力和摩擦因数显著提高，偏位角和质量流量则逐渐减小；
3、瓦块分布位置对三叶动压气体滑动轴承的静态性能影响显著，其中瓦上承载方式的承载力、偏位角和质量流量明显高于瓦间承载方式。关键词：多叶气体轴承；预负荷系数；承载能力；静态特性；有限差分法中图分类号：（，）：，：；随着旋转机械向高转速、高精度、高功率和低能耗的趋势发展，滑动轴承对高速转子的润滑性能已成为影响设备高效稳定工作的关键因素。相较于传统滚动轴承和油液润滑轴承，动压气体轴承因其结构简单、易于制造、对环境无污染及服役温度范围宽等优点，广泛应用于空气压缩机、高速离心机和透平膨胀机等旋转机械中。多叶动压气体轴承采用多块相同几何参数的圆弧状轴瓦来支承转子，其轴向槽间隙可有效防止压力扰动沿周向扩散以提高
4、轴承的稳定性，同时轴向槽的存在还起到收集灰尘和杂质粒子、维持瓦块 轴颈表面清洁的作用。目前，国内外学者针对动压气体轴承的性能开展了大量深入的研究，并取得了丰富的研究成果。赵三星等发现轴向槽圆柱气体动压轴承的承载力和偏位角会随着开槽位置的改变呈周期性变化。和通过耦合求解波箔型径向箔片轴承的气膜压力控制方程和能量方程，研究了轴承间隙、载荷大小和轴颈转速对转轴三维温度场分布的影响，结果表明转子的最高温度随轴承间隙的减小而升高，随载荷和转速的增大而急剧升高。田宇忠等发现轴承的承载能力会随着轴向槽宽度与数目的增大而逐渐降低。等 通过有限元法求解微极性流体润滑 方程，发现承载力、偏位角和摩擦力随预负荷系数
5、和微极性流体耦合数的增大而增大，端泄流量随预负荷系数的增大而减小。燕震雷和伍林结合牛顿迭代法和有限差分法求解考虑气体稀薄效应的修正 方程，探讨了不同边界滑移条件对三瓦可倾瓦动压气体轴承承载力的影响，表明气体稀薄效应明显降低了轴承承载能力，且随着克努森数的增大，连续模型、一阶滑移模型和 的新滑移模型计算的承载力依次降低。冯凯等人建立了考虑接触面库仑摩擦的三瓣式气体箔片轴承气弹耦合模型，研究了转速和轴承预载对承载力的影响，发现承载力随轴承预载和转速的增大而显著提高。和 通过数值求解非牛顿流体润滑的量纲一化 方程，对比了牛顿 非牛顿流体润滑剂润滑下二、三、四叶滑动轴承静态特性曲线，发现当非牛顿润滑剂
6、的膨胀系数为 时三叶轴承拥有最大承载力，而牛顿流体润滑下的二叶轴承具有最大的姿态角和流量系数。等建立了考虑轴颈倾斜、轴瓦弹性变形和瓦块数目的轴向槽轴承弹流润滑模型，运用等效支点位置探究了轴颈倾斜角和轴向槽数目对承载能力的影响，发现承载能力随偏斜角和沟槽数的增加而减小。等结合线性摄动法和有限差分法求解静压气体轴承的非定常 方程，得到了不同供气孔径且不均匀分布的气膜厚度和气膜压力分布，指出承载力随孔口直径的减小而增大，随偏心率、转速和供气压力增大而增大，并发现在小偏心率时轴承具有更好的稳定性。化文灿等基于有限差分法对 种波箔型动压气体轴承的定常 方程和气膜厚度方程联合求解，得到了承载力和摩擦力矩随
7、转速与偏心率的变化规律。由于轴颈 轴瓦间隙、气体可压缩性、气膜不连续性及瓦块分布位置联合作用的复杂性，针对多叶动压气体轴承润滑性能的研究文献相对较少。而研究多叶动压气体轴承的静态特性，对提高气体轴承转子系统稳定性具有重要的指导意义。本文作者通过数学变换将多叶动压轴承的气体润滑 方程转化为标准偏微分方程形式，利用有限差分法和超松弛迭代法进行数值求解，得到三叶气体轴承各瓦块的气膜膜厚及压力分布、承载力和摩擦因数等静态性能，并分析了偏心率、预负荷系数、轴承数、长径比和瓦块分布位置对多叶气体轴承静态性能的影响规律。数学模型图 所示为不同承载方式的三叶动压气体滑动轴承的结构示意图。三叶动压气体滑动轴承是
8、由三块相同几何尺寸的轴瓦依次均匀地分布在支点圆上构成的，瓦块分布位置的改变可以影响主要承载区域，主要承载区域集中在瓦 上时为瓦上承载方式，位于瓦 与瓦 之间时为瓦间承载方式。图中 为瓦间角，为瓦包角。为了计算轴承气膜的承载力等静态性能，必须计算其气膜压力分布。根据气体运动方程、连续性方程、气体状态方程和 方程，可得到可压缩气体润滑 方程的一般形式为（）（）（）式中：为气体密度；为气膜厚度；为气体动力黏度；为轴颈表面速度，；为轴颈转速；为时间；为轴承周向坐标；为轴承轴向坐标。对多叶动压气体轴承而言，其气膜厚度与轴颈中心相对于轴承中心的偏心距有关，在忽略轴瓦厚度的前提下，第 块轴瓦与轴颈表面的气膜
9、厚度表达式为 （）（）（）（）式中：为轴承半径间隙；为装配间隙；为第 块轴瓦的支点位置角；为轴颈中心相对轴承润滑与密封第 卷中心的偏心距；为轴承偏位角。图 三叶动压气体滑动轴承结构 ：（）；（）对动压轴承进行分析计算时，常以量纲一化形式进行，可以将问题归纳为最紧凑形式，突出各有关因素的作用。则等温条件下气体润滑 方程的量纲一化形式如下：（）（）式中：为量纲一气膜厚度；为轴承的量纲一周向角坐标，为轴承半径；为量纲一气膜压力，为环境压力，；为轴承的量纲一轴向坐标；为轴承数，（）。量纲一气膜厚度表达式为 （）（）（）式中：为量纲一气膜厚度；表示轴颈中心相对支点圆中心的偏心率，；为轴瓦的预负荷系数，。
10、量纲一压力边界条件为（，）（，），（）式中：表示轴承宽度；为从角起线到轴瓦进气端的角度；为从角起线到轴瓦出气端的角度。数值求解为了快速求解动压气体润滑 方程，需将其化为标准偏微分方程的形式进行求解。通过简单的数学变换，令，则 （），式（）可变为（）（）同理可得：（）则：（）（）通过移项，多叶动压气体轴承的静态润滑 方程（）可整理为（）利用有限差分法求解动压气体润滑轴承压力分布，就是在瓦块的气膜区域按照差分的步长、划分为许多网格。如图 所示，共划分有（）（）个网格点，每个节点的位置用（，）编号表示，节点（，）上的 值用，表示，得到一组近似表达气膜压力分布的离散压力数值，再利用相应的数值积分便可求
11、得承载力等静态性能。图 多叶气体轴承瓦块的网格划分 年第 期陈 阳等：多叶动压气体滑动轴承静态特性的有限差分算法 采用二阶差分法和中心差分法对式（）进行离散：，（），（），（），），（），（）（），（，），（）（）同理可得：，（），（），（），），（）（），（）（），（，），（），（）（）对于项由半步长差商表达式得：（），（）则量纲一化 方程可化为（），（），（），），（），（）（），（），（），），（），（），（，），（），（，），（），（）将上式化简整理可得：（，），（，），（）（，），（）（，），（）（，），（，）（）经整理成：，（）其中各系数为，（）（，），（）（，），（）（，），（，
12、）（）因此，可表示为，（，），（）任意一点的气膜压力，可表示为，（）轴承在垂直和水平方向的量纲一气膜合力可按下式计算 （）（）（）式中：和 分别为轴承在垂直和水平方向上的气膜合力；为静态气膜合力。轴承的量纲一承载力和偏位角为 （）（）由于气膜的剪切流动与压力流动，轴颈会与轴承内表面会发生摩擦，轴颈表面的周向摩擦因数的表达式为 （）润滑与密封第 卷式中：“”表示摩擦力方向与轴承旋转方向相反。由于气体润滑剂的可压缩性，任一区域流入流出的体积流量不一定守恒，所以用质量流量表征其流动性更为合理。单位时间内流进轴瓦周向截面的质量流量可表示为 （）基于上述模型及数值求解方法，采用有限差分法和超松弛迭代法编
13、制程序，对图 所示的三叶动压气体滑动轴承进行静态性能计算，具体计算流程如图 所示。图 静态性能计算流程 结果与讨论为验证上述理论模型的正确性及编制程序的可靠性，文中 计 算 了 当 ，；，和 ，时的轴承承载力，并与文献 的结果进行了对比，如图 所示。两者结果吻合良好，误差在 以内，说明了文中有限差分算法和程序的正确性。气膜厚度与气膜压力分布在 ，和 时，计算得到不同瓦块分布位置的量纲一气膜厚度和气膜压力分布，如图 所示。瓦块分布位置对气膜厚度和压力分布有着明显影响，沿着轴颈旋转方向，气膜厚度先逐渐增大后减小再增大，且在最小气膜厚度位置处出现最大气膜压力。瓦上承载方式的最小气膜厚度（为 ）大 于
14、 瓦 间 承 载 的 最 小 气 膜 厚 度（为），且最大气膜压力高出瓦间承载方式近。这是因为瓦上承载方式的主要承载区域位于主要承载瓦 上，而瓦间承载方式的主要承载区域位于瓦 和瓦 间，在外载荷作用下润滑气体经轴向槽的泄漏量大于瓦上承载。偏心率和预负荷系数对静态性能的影响图 给出了当 ，时，承载力、偏位角、轴颈表面周向摩擦因数和单位时间内流入瓦 的质量流量 随偏心率和预负荷系数（、和 ）的变化曲线。由图 （）（）可知，随着偏心率和预负荷系数的增大，和迅速增大；同时，在不同预负荷系数条件下，瓦上承载方式的承载能力始终优于瓦间承载方式。从图（）中可以观察到，当偏心率小于 时，瓦上承载和瓦间承载方式
15、的摩擦因数相差不大；在预负荷系数为 和 且偏心率大于 时，瓦间承载方式的略微大于瓦上承载方式；当预负荷系数为 时，小偏心率下的瓦上和瓦间承载方式的摩擦因数几乎一样；但当偏心率为 时，出现瓦上承载方式的摩擦因数高于瓦间承载方式的情况，这是由于高预负荷系数和大偏心率的条件使瓦上承载方式的轴颈表面更贴近主要承载瓦。从图 （）（）可知，瓦上承载方式的偏位角 呈现出随预负荷系数的增大而逐渐减小，随偏心率的增大先增大后逐渐减小的趋势；随着预负荷系数和偏心率的增大，逐渐下降，这是因为轴瓦间隙随着预负荷系数的增大而变小，所以质量流量减小。瓦间承载方式的偏位角 和 随偏心率和预负荷系数的增大而迅速减小。图 承载
16、力计算结果对比 年第 期陈 阳等：多叶动压气体滑动轴承静态特性的有限差分算法 图 瓦上、瓦间承载方式的气膜厚度和压力分布（，）（，）：（）；（）；（）；（）图 偏心率和预负荷系数对静态性能的影响（，）（，）：（）；（）；（）；（）润滑与密封第 卷 轴承数和长径比对静态性能的影响图 示出了当 ，分别为、和 时，轴承数对承载力、摩擦因数和质量流量的影响曲线。由图 （）可知，承载力随 和 的增大而呈现增长速率放缓的变化趋势，这是由于轴承数的增加提高了轴颈转速，增强了气体动压效应，从而使轴承承载力变大；同时，类似承载力随偏心率和预负荷系数的变化曲线，瓦上承载方式的承载力始终高于瓦间承载方式。由图 （）
17、可知，随轴承数和长径比的增加而线性增加，瓦上承载方式的摩擦因数与瓦间承载方式相差不大。由图（）可知，瓦上承载方式的 随轴承数的增大先减小而后缓慢增加，且随长径比的增加而快速增大，这是因为承载瓦块 的面积随长径比的增大而变大。瓦间承载方式的 随轴承数的增加而缓慢增大，随长径比的增加而明显增大。由于瓦间承载方式较大的气体端泄量，瓦上承载方式的 始终高于瓦间承载方式。图 轴承数和长径比对静态性能的影响（，）（，）：（）；（）；（）结论（）长径比和轴承数对多叶动压气体轴承的承载能力影响显著，长径比的增大可以明显提高轴承的承载能力。这是由于长径比改变了瓦块承载面积大小。轴承数增加对应的是轴颈转速增大，转
18、速提高增强了气体动压效应，从而提高了多叶轴承的承载能力。（）预负荷系数可为多叶动压气体轴承的轴瓦提供预紧压力，较大的预负荷系数可以提高轴承的承载力，但也会使轴颈表面周向摩擦因数增大，同时减小了轴瓦间隙，使单位时间内流入瓦 的气体质量流量减小。（）瓦块分布方式对多叶动压气体滑动轴承的静态性能影响显著，瓦上承载方式的承载能力和质量流量要高于瓦间承载方式，瓦间承载方式的轴承在高速运行时会出现气体从瓦间沟槽泄漏的现象，从而使承载力和质量流量下降，因此在设计轴承转子系统时应充分考虑瓦块分布位置的影响。参考文献 谢伟松，林鑫，王伟韬，等航空发动机弹性箔片气体动压轴承技术研究及性能评价综述润滑与密封，（）：
19、，（）：侯予，杨山举，陈兴亚，等箔片动压气体轴承在低温透平膨胀机中的应用哈尔滨工程大学学报，（）：，（）：吕昕多叶式波箔气体动压轴承静特性研究润滑与密封，（）：，（）：，：赵三星，孟光，李友荣轴向槽圆柱气体动压轴承的性能研 年第 期陈 阳等：多叶动压气体滑动轴承静态特性的有限差分算法 究现代制造工程，（）：，（）：，（）：田宇忠，刘正林，刘宇，等带轴向槽水润滑轴承润滑性能数值分析与研究船海工程，（）：，（）：，：，（）：燕震雷，伍林稀薄效应对可倾瓦动压气体轴承性能的影响航空动力学报，（）：，（）：冯凯，胡小强，赵雪源，等三瓣式气体箔片径向轴承的静动态特性中国机械工程，（）：，（）：，：，（）：，（）：，（）：化文灿，杨山举，冶文莲，等制冷剂润滑的径向波箔气体轴承承载特性分析低温与超导，（）：，（）：李莎，吕延军，罗宏博，等轴向槽气体轴承转子非线性系统的动力学行为机械工程学报，（）：，（）：张镜洋，赵晓荣，常海萍，等边界滑移对波箔型动压气体轴承静特性的影响推进技术，（）：，（）：，：，（）：吴炎，谢振宇，陈李成，等径向磁气组合轴承系统的动态性能航空动力学报，（）：，（）：杨利花，戚社苗，耿海鹏，等可倾瓦气体动压轴承承载能力的数值分析润滑与密封，（）：，（）：虞烈，戚社苗，耿海鹏可压缩气体润滑与弹性箔片气体轴承技术北京：科学出版社，润滑与密封第 卷
