高级检索

  • ISSN 1006-3080
  • CN 31-1691/TQ
引用本文:
Citation:

考虑波箔变形的波箔型气体箔片轴承润滑性能的数值研究

    作者简介: 许浩杰(1994-),男,浙江宁波人,硕士生,研究方向为机械设计及理论。E-mail:13162583657@163.com;
    通讯作者: 安琦, anqi@ecust.edu.cn
  • 中图分类号: TH133.35;TH117.2

Numerical Study on the Lubrication Performance of Bump-Type Gas Foil Bearings with Considering the Deformation of Bump Foil

    Corresponding author: Qi AN, anqi@ecust.edu.cn
  • CLC number: TH133.35;TH117.2

  • 摘要: 以波箔型气体箔片轴承为研究对象,运用弹性力学理论对波箔的力学性能进行了研究,构建了其受力变形模型,并结合可压缩气体Reynold方程,构建了一种流固耦合算法。同时结合具体算例,对比研究了弹性箔片气体轴承和刚性气体轴承的气膜压力分布及气膜厚度分布,发现两者分布形状相似,但弹性箔片气体轴承气膜压力的最大值较小,轴向平分面处的最小气膜厚度较大;偏心率一定时,弹性箔片气体轴承的摩擦力矩小于刚性气体轴承摩擦力矩,且气体端泄流量大于刚性气体轴承,故弹性箔片气体轴承运行时的摩擦力矩较小,发热较小,且散热较快。偏心率一定时,增大波箔厚度tB和减小波箔波距s可减小轴承量纲为一最小气膜厚度,并减小轴承偏位角;减小tB和增大s可减小轴承摩擦力矩,有利于降低摩擦功耗,减少摩擦热,且可增大轴承气体端泄流量,有利于散热。
  • 图 1  波箔型气体箔片轴承结构

    Figure 1.  The structure of bump-type gas foil bearings

    图 2  箔片结构受力分析

    Figure 2.  Force analysis of bump foil structure

    图 3  单波顶点位移计算分析

    Figure 3.  Calculation and analysis of single wave’s peak displacement

    图 4  网格划分及差商示意图

    Figure 4.  The schematic diagram of mesh generation and difference quotient

    图 5  计算流程图

    Figure 5.  Calculation flow chart

    图 6  无量纲气膜压力分布对比图

    Figure 6.  The comparison of dimensionless gas film pressure distribution

    图 7  轴向平分面处无量纲气膜压力与气膜厚度对比图

    Figure 7.  The comparison of dimensionless gas film pressure and thickness at the axial bisector

    图 8  不同偏心率下摩擦力矩与气体端泄流量对比图

    Figure 8.  Comparison of friction torque and side leakage flow of gas with different eccentricities

    图 9  不同偏心率下波箔厚度和波距对气膜承载力的影响

    Figure 9.  Influences of thickness and wave distance of bump foil on gas foil capacity with different eccentricities

    图 10  不同偏心率下波箔厚度和波距对偏位角的影响

    Figure 10.  Influences of thickness and wave distance of bump foil on attitude angle with different eccentricities

    图 11  不同偏心率下波箔厚度和波距对量纲为一般最小气膜厚度的影响

    Figure 11.  Influences of thickness and wave distance of bump foil on the dimensionless minimum gas film thickness with different eccentricities

    图 12  不同偏心率下波箔厚度和波距对摩擦力矩的影响关系图

    Figure 12.  Influences of thickness and wave distance of bump foil on friction torque with different eccentricities

    图 13  不同偏心率下波箔厚度和波距对气体端泄流量的影响关系图

    Figure 13.  Influence of thickness and wave distance of bump foil on side leakage flow of gas with different eccentricities

    表 1  轴承结构及润滑气体参数

    Table 1.  Parameters of bearing structure and lubricating gas

    $R$/m$L$/m$c$/m${t_B}$/m${h_B}$/m$s$/m${R_0}$/m
    ${\rm{5}} \times {10^{ - 2}}$$7.5 \times {10^{ - 2}}$$1.0 \times {10^{ - 4}}$$7.62 \times {10^{ - 5}}$$6.3 \times {10^{ - 4}}$$4.064 \times {10^{ - 3}}$$2.655 \times {10^{ - 3}}$
    ${\rm{2}}\theta $/(°) $2l$/m $E$/Pa $\mu $/(Pa·s) ${p_0}$/Pa $\Re $/(J·kg−1·K−1) T/K
    1.407 ${\rm{3}}{\rm{.434}} \times {10^{ - 3}}$ $2.07 \times {10^11}$ $1.932 \times {10^{ - 5}}$ $1.01325 \times {10^5}$ 287 298
    下载: 导出CSV
  • [1] CHRISTOPHER D. Oil-free shaft support system rotordynamics: past, present and future challenges and opportunities[J]. Mechanical Systems and Signal Processing, 2012, 29: 67-76. doi: 10.1016/j.ymssp.2011.07.024
    [2] BRANAGAN M, GRIFFIN D, GOYNE C, et al. Compliant gas foil bearings and analysis tools[J]. Journal of Engineering for Gas Turbines and Power, 2016, 138(5): 054001.1-054001.8.
    [3] 徐磊, 焦映厚, 陈照波, 等. 波箔型径向气体轴承静态特性分析[J]. 润滑与密封, 2012, 37(8): 27-32.
    [4] WALOWIT J A, ANNO J N. Modern developments in lubrication mechanics[M]. London: Applied Science Publishers Ltd., 1975:214-221.
    [5] HESHMAT H, WALOWIT J A, PINKUS O. Analysis of gas-lubricated foil journal bearings[J]. Journal of Lubrication Technology, 1983, 105(4): 647-655. doi: 10.1115/1.3254697
    [6] PENG Z-C, KHONSARI M M. Hydrodynamic analysis of compliant foil bearings with compressible air flow[J]. Journal of Tribology, 2004, 126(3): 542-546. doi: 10.1115/1.1739242
    [7] LEE D H, KIM Y C, KIM K W. The static performance analysis of foil journal bearings considering three-dimensional shape of the foil structure[J]. Journal of Tribology, 2008, 130(3): 031102-031111. doi: 10.1115/1.2913538
    [8] 龚焕孙, 孙鸿兴, 周国明. 径向波箔轴承的设计计算和试验研究-(Ⅰ)理论分析[J]. 上海理工大学学报, 1991, 13(4): 1-17.
    [9] 许怀锦, 刘占生, 张广辉, 等. 波箔片气体动压径向轴承气膜压力分布特点[J]. 轴承, 2008(6): 23-27. doi: 10.3969/j.issn.1000-3762.2008.06.009
    [10] LE LEZ S, ARGHIR M, FRENE J. A new bump-type foil bearing structure analytical model[J]. Journal of Engineering for Gas Turbines and Power, 2007, 129(4): 1047-1057. doi: 10.1115/1.2747638
    [11] 许怀锦.转子-箔片轴承系统动力学特性理论及试验研究[D].哈尔滨工业大学, 2009.
    [12] 李姗姗, 安琦. 考虑微观形貌时平面型推力滑动轴承润滑性能数值研究[J]. 华东理工大学学报(自然科学版), 2017, 43(5): 733-739.
  • [1] 粟爽格安琦 . 内外圈装配过盈量对圆柱滚子轴承力学性能的影响. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190121002
    [2] 潘斌翁涛杨家鹏安琦 . 汽车发动机皮带系统臂式张紧轮力学性能研究. 华东理工大学学报(自然科学版), doi: 10.14135/j.cmki.1006-3080.20190114006
    [3] 罗雪莉郎美东 . 聚合物前药载药性能的计算机模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190514001
    [4] 罗安王汉奎王建文 . 基于小冲杆试验数据的力学性能的数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180609002
    [5] 徐田刘艺军何颖程起林 . CoNi-MOF/RGO复合物纳米片的制备及其储锂性能. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190329001
    [6] 高源安琦 . 轴承座同心度误差对深沟球轴承-转子系统振动性能的影响. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180415001
    [7] 熊文莉杨家鹏安琦 . 轴承滚子加工精度对高速电主轴动力学性能的影响. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180613002
    [8] 许叶龙刘迎圆惠虎於晔鸿 . 固液混合过程的数值模拟及实验研究. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180709001
    [9] 赖兆林冯翔虞慧群 . 基于逆向学习行为粒子群算法的云计算大规模任务调度. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190218001
    [10] 叶贞成倪泽雨程辉 . 一种改进的分隔壁精馏塔简捷计算方法. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191010001
    [11] 熊祥军邵方明张祖渊 . 直径限定可靠性计算的冗余边的检测算法. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.
    [12] 徐海涛缪宏钢翁涛安琦 . 汽车发动机张紧轮主轴力学分析及其疲劳寿命计算方法. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191107001
    [13] 胡贵华隆建杜文莉 . 考虑汽化影响的乙烯裂解炉对流段数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180913002
    [14] 蔡冬莹王跃林段先建胡丹李玉冰胡彦杰 . 预混合高速射流燃烧器内温度场的数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190221004
    [15] 柏德鸿宗原赵玲 . 扭曲片强化裂解炉管多尺度混合特性分析. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190315003
    [16] 于苗钱锋胡贵华隆建李天越 . 工业级MIP提升管反应器气固两相流动特性的数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180913003
    [17] 贾翔飞钱嘉澍吴幼青陈剑佩 . 鼓泡塔反应器中两相流动CFD-PBM耦合数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191121001
    [18] 陆俊杰吴晨安琦 . 基于微观接触力学的旋塞阀密封面疲劳寿命数值研究. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190103001
    [19] 杨家鹏熊文莉安琦 . 滚子直径误差对双列调心轴承疲劳寿命的影响. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180411001
    [20] 顾烨晨杨家鹏安琦 . 空调贯流风扇转子振动对支承轴承受力及寿命影响. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180604001
  • 加载中
图(13)表(1)
计量
  • 文章访问数:  1244
  • HTML全文浏览量:  133
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-01-08
  • 网络出版日期:  2020-01-21

考虑波箔变形的波箔型气体箔片轴承润滑性能的数值研究

    作者简介:许浩杰(1994-),男,浙江宁波人,硕士生,研究方向为机械设计及理论。E-mail:13162583657@163.com
    通讯作者: 安琦, anqi@ecust.edu.cn
  • 1. 华东理工大学机械与动力工程学院,上海 200237
  • 2. 长城汽车股份有限公司技术中心,保定 071000

摘要: 以波箔型气体箔片轴承为研究对象,运用弹性力学理论对波箔的力学性能进行了研究,构建了其受力变形模型,并结合可压缩气体Reynold方程,构建了一种流固耦合算法。同时结合具体算例,对比研究了弹性箔片气体轴承和刚性气体轴承的气膜压力分布及气膜厚度分布,发现两者分布形状相似,但弹性箔片气体轴承气膜压力的最大值较小,轴向平分面处的最小气膜厚度较大;偏心率一定时,弹性箔片气体轴承的摩擦力矩小于刚性气体轴承摩擦力矩,且气体端泄流量大于刚性气体轴承,故弹性箔片气体轴承运行时的摩擦力矩较小,发热较小,且散热较快。偏心率一定时,增大波箔厚度tB和减小波箔波距s可减小轴承量纲为一最小气膜厚度,并减小轴承偏位角;减小tB和增大s可减小轴承摩擦力矩,有利于降低摩擦功耗,减少摩擦热,且可增大轴承气体端泄流量,有利于散热。

English Abstract

  • 气体箔片轴承是一种较为新型的轴承,特别适用于高速透平机这类机械[1],它的结构有多种不同的型式,其中波箔型气体箔片轴承应用最为广泛[2]。波箔型气体箔片轴承主要由壳体、波箔片(底箔)、平箔片(顶箔)3部分组成[3],工作过程中,气膜压力的产生将使波箔产生弹性变形,目前在箔片轴承的性能计算中鲜有对这个弹性变形进行深入的研究。

    Walowit和Anno[4]于1975年首次从理论上分析了波箔型气体箔片轴承的结构,建立了箔片弹性模型,将箔片的弹性系数代入一维雷诺方程,研究了轴承的静态性能。Heshmat等[5]将波箔片简化为线性弹簧并引入波箔片变形系数α,构建了箔片结构等效模型,利用数值方法研究了波箔型气体箔片轴承的静态特性。Peng等[6]借鉴Heshmat提出的波箔弹簧模型,建立了用来预测波箔型气体径向轴承流体动力学特性的模型,耦合求解了可压缩气体雷诺方程和气膜厚度方程,所得的轴承在一系列转速下的承载能力预测值与已有实验数据吻合较好。Lee等[7]建立了考虑波箔片三维形状的有限元结构模型,比较了相互连接的波箔的变形与离散波箔的变形,发现顶箔和底箔厚度过小可导致轴承承载能力显著降低。龚焕孙等[8]在假设波箔局部刚度相同的前提下,利用圆柱薄壳线性理论推导出了单波的刚度公式,给出了波箔整体刚度计算公式。许怀锦等[9]对比了不考虑箔片变形、考虑其弹性变形时的波箔型气体动压轴承的不同特性,分析了不同转速和偏心率对气膜厚度分布的影响规律。Sébastien等[10]将箔片轴承结构看成波箔各拱起之间互相作用的多自由度系统,每个拱起均包含3个与基本弹簧相连的自由度,由此提出了一种新模型。许怀锦[11]比较了4种基于箔片结构刚度和阻尼的分析模型及相对应的气膜厚度方程,给出了每种模型的适用情况及优缺点。

    综上所述,目前的研究中一般将波箔简化为线性弹簧以构建箔片弹性模型,且多数将波箔的周向变形忽略,这是一种近似研究方法,精度受到限制。实际中,波箔的变形十分复杂,既有径向弹性变形,也有沿周向的弹性变形,采用简化的方法对其力学性能进行描述难以达到较高的精度。为此,本文运用弹性力学理论对波箔的力学性能进行研究,构建其力学模型,结合可压缩流体Reynold方程,构建一种流固耦合算法,同时采用数值方法研究了波箔变形对该轴承润滑性能的影响。

    • 波箔型气体箔片轴承结构如图1所示,它主要由壳体、波箔(底箔)、平箔(顶箔)3部分组成,两箔片的一端共同固定于壳体上,另一端自由搭接。转速较低的启停阶段,轴颈与顶箔直接接触,顶箔上的减摩材料能减少两者之间的摩擦;转速较高的正常运转阶段,顶箔和轴颈之间形成的楔形润滑气膜,为轴颈提供支撑力。

      图  1  波箔型气体箔片轴承结构

      Figure 1.  The structure of bump-type gas foil bearings

      箔片结构受力分析如图2所示,气膜压力通过平箔传递给波箔,波箔在此压力作用下发生弹性变形,虚线部分为变形前的状态,实线部分为变形后的状态,波纹原半长从l伸长为$l + \Delta x$,平箔随波箔顶点而发生的竖直变形量为$\Delta y$。由于波箔沿轴承圆周方向的结构相同,且一端固定,一端自由,因此对单个波纹进行的受力分析可推广到整个波箔。

      图  2  箔片结构受力分析

      Figure 2.  Force analysis of bump foil structure

      在进行力学分析时假设如下:

      (1) 波箔的顶部和平箔之间始终保持接触,单波所受载荷为作用于波峰的集中力;

      (2) 波箔受力沿轴向(宽度方向)均匀,各处径向变形沿轴向均匀;

      (3) 由于波箔和平箔、波箔和壳体之间的摩擦很小,故将其忽略;

      (4) 波箔在受到平箔的压力变形后,只是半径增大,仍然保持为圆弧形。

      单箔顶点位移计算分析如图3所示,单波横截面尺寸为${t_B} \times L$L为轴承宽度),其顶部受到竖直向下的载荷F,由于其横截面高度${t_B}$远小于其轴线半径${R_{\rm{0}}}$,因此可利用卡氏定理计算其变形。首先计算单波顶点的竖直位移$\Delta y$,若规定使单波曲率增大的弯矩为正,则单波任意截面m-m上的弯矩为

      图  3  单波顶点位移计算分析

      Figure 3.  Calculation and analysis of single wave’s peak displacement

      其关于F的偏导数为

      对单波左半段进行积分,并将式(1)、(2)代入可得单波顶点的竖直位移为

      采用附加力法计算单波顶点的水平位移$\Delta x$,在顶点上增加一水平向右的附加力$F{}_a$,于是任意截面m-m上的弯矩为

      其关于$F{}_a$的偏导数为

      式中,l为波箔的波半场;F为单波顶部受到的竖直向下的载荷;α为截面m−m所对应的角度;2θ为波箔的尺寸参数。表示单波圆弧的角度。单波顶点的水平位移为

      单波横截面惯性矩可表示为

      将式(7)代入式(3)、(6),可得

    • 将可压缩流体的运动方程、连续性方程和能量守恒方程联立,得到气体动压润滑Reynolds方程:

      式中,x为轴承的周向坐标,z为轴承的轴向坐标,$\mu $为气体的动力黏度,h为气膜厚度,p为气膜压力,$\rho $为气体密度,t为运行时间,U为轴颈的周向速度。

      理想气体状态方程为:

      上式中,$\Re $为气体常数,其值取决于气体种类,T为热力学温度。

      假设在轴承工作过程中气体温度不变,且忽略压强对气体黏度的影响,将式(11)代入上述Reynolds方程,并将气体视为定常,消去时间项,得到用来描述平箔与轴颈之间气体流动规律的可压缩气体等温定常Reynolds方程:

      将式(12)转换为柱坐标系下的形式,并进行量纲为一化,得

      式中,$x = R\varphi $$z = \lambda \dfrac{L}{2}$$h = Hc$$P = \dfrac p{p_0}$$\varLambda = \dfrac{{6U\mu R}}{{{c^2}{p_0}}}$,其中P为量纲为一的气膜压力,H为量纲为一的气膜厚度,$R$为轴颈半径,$\varphi $为圆周角度,$L$为轴承宽度,$c$为轴承半径间隙,${p_0}$为环境压力,$\varLambda $为轴承系数,$ \lambda $为轴承宽度方向的量纲一的长度。

    • 在工作过程中,波箔在气膜压力下的径向变形会影响气膜厚度分布,从而影响气膜压力分布,气膜压力分布又会反过来影响波箔的径向变形量,这是一流固耦合过程。由于平箔的变形相对于波箔变形来说很小,为便于分析,假设两单波之间的平箔自身不发生变形,而是随波箔的变形而整体发生位移。考虑波箔弹性变形的气膜厚度方程可表示为:

      式中,e为轴承偏心距,将式(8)代入,可得:

      量纲为一化得:

      式中,$\varepsilon $为轴承偏心率。

    • 利用有限差分法求解式(13)所示的量纲为一Reynolds方程,可得气膜压力P[12]。网格划分及差商示意图如图4所示,轴承气膜承载区以波箔固定端和自由端为周向边界(本文假设该区域为$\varphi \in [0,1.8 {\text{π}} ]$),以轴承两端部为轴向边界。将轴承气膜区域离散为$m \times n$的网格,用$i$表示$\varphi $方向的列数,用$j$表示$\lambda $方向的行数,用$(i,j)$表示每个节点的编号,压力值用${P_{i,j}}$表示。沿$\varphi $方向共$m$格,$i$的编号为1~(m+1),每格初始步长为:

      图  4  网格划分及差商示意图

      Figure 4.  The schematic diagram of mesh generation and difference quotient

      $\lambda $方向共$n$格,$j$的编号为1~(n+1),每格步长相等:

      本文采用中差商表示导数:

      将式(19)~(21)代入式(13),化简得:

      式中:${A_{i,j}} = \dfrac{{{P_{i + 1/2,j}}H_{_{i + 1/2,j}}^3}}{{\Delta {\varphi _{i + 1/2}}}}$${B_{i,j}} = \dfrac{{{P_{i - 1/2,j}}H_{_{i - 1/2,j}}^3}}{{\Delta {\varphi _{i - 1/2}}}}, {C_{i,j}} = $${\left( {\dfrac{{2R}}{{L\Delta \lambda }}} \right)^2}\Delta {\varphi _i}{P_{i,j + 1/2}}H_{_{i,j + 1/2}}^3$${D_{i,j}} = {\left( {\dfrac{{2R}}{{L\Delta \lambda }}} \right)^2}\Delta {\varphi _i}{P_{i,j - 1/2}}H_{_{i,j - 1/2}}^3$${E_{i,j}} = {A_{i,j}} + {B_{i,j}} + {C_{i,j}} + {D_{i,j}}\; ,$${F_{i,j}} = \Lambda \left( {{P_{i + 1/2,j}}{H_{_{i + 1/2,j}}} - {P_{i - 1/2,j}}{H_{_{i - 1/2,j}}}} \right)$ 

      求解过程中引入边界条件,可压缩气体润滑中不存在气膜破裂现象,轴承端部及箔片固定端和自由端均与大气环境连通,其压力可取为环境压力,其边界条件量纲为一的形式为:

      引入如下相对收敛准则:

      本文的Eerror${\rm{1}}{{\rm{0}}^{{\rm{ - 5}}}}$,表示迭代收敛允许的误差。式中k为迭代次数。

    • 量纲为一的气膜承载力W

      轴承偏位角$\beta $

      式中,${W_x}$为垂直于偏位线的分力,${W_y}$为平行于偏位线的分力,可表示为:

      作用于轴颈表面的摩擦力矩为:

      气体端泄流量为:

      本文流固耦合的数值计算过程如图5所示:首先对气膜压力P和波箔径向弹性变形量V赋初值,对气膜压力承载区进行网格划分,采用有限差分法迭代求解Reynolds方程,得到一组新的压力分布,计算该压力分布下的波箔径向和周向弹性变形量,波箔沿圆周方向伸长$R\Delta \varphi $,计算区域从原来的$\varphi \in [0,1.8\text{π} ]$变为$\varphi \in [0,1.8\text{π} + \Delta \varphi ]$,据此更新气膜厚度并重新进行网格划分,在此基础上求解又可得一组气膜压力分布,如此不断迭代直至收敛最终得到第k次迭代所得各点气膜压力绝对值之和T和各点气膜压力相对于上一次迭代结果的变化量之和S

      图  5  计算流程图

      Figure 5.  Calculation flow chart

    • 根据上述计算方法,以某一具体波箔型气体箔片轴承为算例进行研究,轴承结构及润滑气体参数如表1所示。

      $R$/m$L$/m$c$/m${t_B}$/m${h_B}$/m$s$/m${R_0}$/m
      ${\rm{5}} \times {10^{ - 2}}$$7.5 \times {10^{ - 2}}$$1.0 \times {10^{ - 4}}$$7.62 \times {10^{ - 5}}$$6.3 \times {10^{ - 4}}$$4.064 \times {10^{ - 3}}$$2.655 \times {10^{ - 3}}$
      ${\rm{2}}\theta $/(°) $2l$/m $E$/Pa $\mu $/(Pa·s) ${p_0}$/Pa $\Re $/(J·kg−1·K−1) T/K
      1.407 ${\rm{3}}{\rm{.434}} \times {10^{ - 3}}$ $2.07 \times {10^11}$ $1.932 \times {10^{ - 5}}$ $1.01325 \times {10^5}$ 287 298

      表 1  轴承结构及润滑气体参数

      Table 1.  Parameters of bearing structure and lubricating gas

      图6所示为偏心率$\varepsilon = 0.{\rm{6}}$、转速为100000 r/min的工况下,不考虑波箔变形(刚性轴承)及考虑波箔变形(弹性箔片轴承)时的量纲为一的气膜压力分布对比图;图7为轴承轴向平分面处量纲为一的气膜压力与气膜厚度对比图;图8为不同偏心率下轴承摩擦力矩与气体端泄流量对比图。

      图  6  无量纲气膜压力分布对比图

      Figure 6.  The comparison of dimensionless gas film pressure distribution

      图  7  轴向平分面处无量纲气膜压力与气膜厚度对比图

      Figure 7.  The comparison of dimensionless gas film pressure and thickness at the axial bisector

      图  8  不同偏心率下摩擦力矩与气体端泄流量对比图

      Figure 8.  Comparison of friction torque and side leakage flow of gas with different eccentricities

      对比图6(a)图6(b)并结合图7(a)可知,不考虑箔片弹性变形与考虑箔片弹性变形时的量纲为一的气膜压力分布形状相似,但前者气膜压力最大值大于后者,说明弹性箔片气体轴承的承载力小于刚性气体轴承。由图7(b)可知,两者轴向平分面处的无量纲气膜厚度变化趋势相似,但前者最小气膜厚度值小于后者。由图8可知,偏心率一定时,弹性箔片轴承的摩擦力矩小于刚性气体轴承,且气体端泄流量大于刚性气体轴承,说明弹性箔片轴承运行时的摩擦力矩较小,摩擦发热小,且散热较快。

      图9图13为其他参数均不改变时,不同偏心率$\varepsilon $下波箔厚度${t_B}$和波箔波距$s$对轴承各润滑性能参数的影响关系图。由图可知:当偏心率一定时,量纲为一的气膜承载力随着波箔厚度的增大而增大,随着波箔波距的增大而减小,且当偏心率较大时,影响较明显;偏位角随着波箔厚度增大而缓慢减小,随着波箔波距增大而缓慢增大,且当偏心率较大时,影响较明显;量纲为一的最小气膜厚度随着波箔厚度增大而减小,随着波箔波距增大而增大;摩擦力矩随着波箔厚度的增大而增大,随着波箔波距的增大而减小,且当偏心率较大时,影响较明显;气体端泄流量随着波箔厚度增大而缓慢减小,随着波箔波距增大而缓慢增大。

      图  9  不同偏心率下波箔厚度和波距对气膜承载力的影响

      Figure 9.  Influences of thickness and wave distance of bump foil on gas foil capacity with different eccentricities

      图  10  不同偏心率下波箔厚度和波距对偏位角的影响

      Figure 10.  Influences of thickness and wave distance of bump foil on attitude angle with different eccentricities

      图  11  不同偏心率下波箔厚度和波距对量纲为一般最小气膜厚度的影响

      Figure 11.  Influences of thickness and wave distance of bump foil on the dimensionless minimum gas film thickness with different eccentricities

      图  12  不同偏心率下波箔厚度和波距对摩擦力矩的影响关系图

      Figure 12.  Influences of thickness and wave distance of bump foil on friction torque with different eccentricities

      图  13  不同偏心率下波箔厚度和波距对气体端泄流量的影响关系图

      Figure 13.  Influence of thickness and wave distance of bump foil on side leakage flow of gas with different eccentricities

    • (1) 以波箔型气体箔片轴承为研究对象,运用弹性力学理论对波箔的力学性能进行了研究,构建了其受力变形模型,结合可压缩流体Reynold方程,构建了一种流固耦合算法,借助Matlab编程实现了在考虑波箔变形的情况下对轴承润滑性能的精确计算。

      (2) 结合一具体算例,对比研究了刚性气体轴承和弹性箔片气体轴承的气膜压力分布及气膜厚度分布。发现两者量纲为一的气膜压力分布形状相似,但前者气膜压力最大值大于后者,说明弹性箔片气体轴承的承载力小于刚性气体轴承。两者轴向平分面处的量纲为一的气膜厚度变化趋势相似,但前者最小气膜厚度值小于后者。偏心率一定时,弹性箔片轴承的摩擦力矩小于刚性气体轴承,且气体端泄流量大于刚性气体轴承,说明弹性箔片轴承运行时的摩擦力矩较小,效率较高,且散热较快,有利于保持工作温度的稳定。

      (3) 结合算例进一步研究了波箔厚度${t_B}$、波箔波距s对波箔型气体箔片轴承各润滑性能参数的影响规律。结果表明:偏心率一定时,增大${t_B}$、减小s会减小轴承量纲为一的最小气膜厚度,减小轴承偏位角,有利于提高轴承承载力,提高承载稳定性;减小${t_B}$、增大s可减小轴承摩擦力矩,有利于降低摩擦功耗,减少摩擦热,且可增大轴承气体端泄流量,有利于散热。

(13)  表(1) 参考文献 (12) 相关文章 (20)

目录

    /

    返回文章