浙江大学学报(工学版), 2024, 58(9): 1970-1980 doi: 10.3785/j.issn.1008-973X.2024.09.021

机械工程

三楔气体轴承的动静压特性分析

章晨昕,, 洪伟荣,, 郑水英

浙江大学 化工机械研究所,浙江 杭州 310027

Hydrostatic and hydrodynamic characteristic analysis of three-wedge gas bearing

ZHANG Chenxin,, HONG Weirong,, ZHENG Shuiying

Institute of Process Equipment, Zhejiang University, Hangzhou 310027, China

通讯作者: 洪伟荣,男,教授. orcid.org/0000-0002-8979-0488. E-mail: hongwr@zju.edu.cn

收稿日期: 2023-07-12  

Received: 2023-07-12  

作者简介 About authors

章晨昕(1993—),女,博士生,从事高速旋转机械研究.orcid.org/0000-0003-0950-7724.E-mail:chenxin_zhang@126.com , E-mail:chenxin_zhang@126.com

摘要

针对氢燃料电池汽车车载空压机的轴承-转子系统,提出兼具静压和动压特性的三楔式气体轴承. 开展轴承动静压耦合作用机制和流场特性、动静态参数研究,并与传统的圆柱气体轴承进行对比. 为了计算轴承的流场特性,在Fluent软件的用户自定义函数(UDF)中编写动网格程序,并提出相应的计算流体力学 (CFD)瞬态模拟方法来求解轴颈任意偏心位置时的轴承非线性气膜力与动态参数. 研究表明,在运行工况变化时,三楔气体轴承动静压效应的强度会不断变化,但始终能通过耦合效应为转子提供有效、稳定的支承. 在低转速范围内,静压效应是影响三楔气体轴承性能的主要因素,表现出与圆柱气体轴承相似的支承效果,保证车载压缩机在启停或升速过程中的稳定运转. 当转速增大时,三楔气体轴承的动压效应显著增强,能在更低的供气压力下为高速旋转的车载压缩机提供支承. 相较于需要外接气源的圆柱气体轴承,三楔气体轴承能在仅利用氢燃料电池汽车供气系统内部气路的情况下为转子提供稳定的支承,适应车载空压机中气体压力有限且须频繁启停、变速的工作环境.

关键词: 气体轴承 ; 动静压耦合 ; 动压效应 ; 轴承特性 ; 氢燃料电池汽车

Abstract

A three-wedge gas bearing, with both hydrostatic and hydrodynamic pressure characteristics, was specially designed for the bearing-rotor system of the air compressor in hydrogen fuel cell vehicle. The mechanism of hydrostatic and hydrodynamic coupling effect, the flow field characteristics, and the dynamic and static parameters of the bearing were investigated. These results were also compared with those of traditional cylindrical gas bearings. To analyze the flow field characteristics of the bearing, a dynamic mesh program was developed using UDF in Fluent software, and a related CFD transient simulation method was proposed to solve the nonlinear gas film forces and the dynamic parameters of the bearing at any eccentric journal position. The analysis showed that the hydrodynamic and hydrostatic pressure effect intensity of the three-wedge bearing varied with the external working condition, but kept on coupling with each other, ensuring that the bearing always provided effective and stable support for the rotor. At low rotational speed, the hydrostatic pressure effect was the main factor affecting bearing performance, and the three-wedge gas bearing exhibited similar working performance to cylindrical bearing. The three-wedge gas bearing can thus provide stable support for on-board compressor rotor during start-stop and acceleration. With the increase of the rotational speed, the hydrodynamic pressure effect of the three-wedge gas bearing was significantly enhanced. The three-wedge gas bearing can thus provide support for high-speed compressors with lower gas supply pressure than the cylindrical bearing. Therefore, compared with the traditional cylindrical gas bearing that required external gas sources, the three-wedge gas bearing can provide a stable support for the rotor only using the internal gas circuit of the hydrogen fuel cell vehicle. It adapted well with the working environment inside the hydrogen fuel cell vehicle where the gas supply pressure was limited, and the rotor underwent frequent start-stop and speed changes.

Keywords: gas bearing ; hydrostatic and hydrodynamic pressure coupling effect ; dynamic pressure effect ; bearing characteristics ; hydrogen fuel cell vehicle

PDF (3556KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

章晨昕, 洪伟荣, 郑水英. 三楔气体轴承的动静压特性分析. 浙江大学学报(工学版)[J], 2024, 58(9): 1970-1980 doi:10.3785/j.issn.1008-973X.2024.09.021

ZHANG Chenxin, HONG Weirong, ZHENG Shuiying. Hydrostatic and hydrodynamic characteristic analysis of three-wedge gas bearing. Journal of Zhejiang University(Engineering Science)[J], 2024, 58(9): 1970-1980 doi:10.3785/j.issn.1008-973X.2024.09.021

氢燃料电池汽车几乎无污染排放、燃料来源广泛、充能快速,被称为绿色新型环保汽车. 空气压缩机作为燃料电池阴极气体供应系统的重要部件,向电堆提供高压空气作为反应原料,有无油、高效、小型化的要求[1].

在燃料电池汽车车载空压机的支承选择方面,对气体润滑、水润滑、磁悬浮等类型轴承都有相关研究[2-4]. 其中,气体轴承具有耐寒耐热、摩擦小、无污染、寿命长等优点,且能在高转速、高精度的要求下工作,成为此类空压机的首选支承部件之一. 根据中国最新研发计划,基于对空压机技术要求和寄生功率的综合考量,燃料电池车车载空压机的压缩比一般选取为2.5,即能输出0.25 MPa的高压气体[5]. 在这一供气压力下,传统的静压圆柱气体轴承无法达到高速空压机转子所需的支承能力和稳定性要求,须外接高压气源.

空气箔片轴承是当前的主流选择,因其在高速下的功率损耗较低,相关的空气回路也相对简单. Daejong[6]利用高速旋转试验台测试2种不同的多瓦箔片轴承的性能以及对相应轴系进行转子动力学分析,选取测试中性能最优的轴承提出单级空压机设计;Zhang等[7] 针对120 kW氢燃料电池系统,设计22 kW离心式空压机及配套的箔片式轴承,搭建相应的试验台对空压机以及轴承性能进行测试验证;张雯等[8]针对野外战车应急发电机组的氢燃料电池空压机提出设计空气动压箔片轴承-转子系统的方法,并对相应转子系统进行实验测试;冯凯等[9]建立气体箔片轴承刚度预测模型,搭建刚度测试试验台验证模型的有效性;王伟等[10]将自主研发设计的鼓泡型箔片全动压气体轴承在微型透平膨胀机上进行试验研究;靳彩妍等[11]对箔片轴承的动态特性进行研究,并对其支承的两级空压机轴系进行转子动力学分析. 然而,箔轴承的启动扭矩高,整体尺寸和重量大[12]. 在建立有效的空气压力支撑之前,箔轴承必须达到最小“起飞”转速,在压缩机的启停与低速运转过程中会产生磨损.

为了解决这一问题,本研究提出适用于氢燃料电池汽车车载高速离心空压机的新型三楔气体轴承. 不管是在启停、变速还是稳定高速运转阶段,该轴承都能在仅利用车载空压机内部气路的前提下,为车载空压机转子提供有效、稳定的支承. 在对该轴承进行性能研究的过程中,为了求解其非线性气膜力,提出由自编的用户自定义函数(user-defined functions, UDF)动网格程序为基础的计算流体力学(computational fluid dynamics, CFD)计算方法,能在轴颈中心移动的过程中不断更新计算模型的网格,保证瞬态计算中的模型网格质量. 应用该方法研究三楔气体轴承在不同的偏心率和转速下的流场特性、动静压耦合效应,以及静、动态参数的变化,深入解析三楔动压气体轴承的工作机理和性能优势.

1. 三楔气体轴承流场建模

1.1. 流场结构

三楔式气体轴承的内表面由沿圆周方向首尾相接的3个均匀的圆弧形楔槽构成,如图1所示(为了示意清楚,图中对气膜间隙进行了放大). 当轴颈中心位于轴承几何中心时,轴颈表面到气楔弧面的最小距离,称为半径间隙c*,在数值上等于设计气膜厚度的最小值. 与此相对的,轴颈表面到气楔弧面的最大距离,称为半径楔隙c,等于设计气膜厚度的最大值[13].

图 1

图 1   三楔气体轴承及其流场结构示意图

Fig.1   Structural diagram of three-wedge gas bearing and its flow field


当轴高速旋转时,气膜中的高压气体就会跟随轴颈以一定的速度通过厚度变化的楔形间隙. 随着间隙大小的变化,气膜流场中会出现收敛和发散的结构. 由于流体黏度的存在,通过收敛楔的流体会在运动法向上产生一个对转子的支持力. 这种现象被称为动压效应[14],而楔槽的存在加强了动压效应的形成.

在高速离心空压机启停过程中,外部气源气体通过微孔产生节流作用后,进入三楔轴承的润滑间隙,形成高压气膜支撑转子,使轴承呈现出静压气体轴承的特性,包括稳定性强、承载力大和刚度大等. 在转速变化的过程中,由于转子始终处于悬浮状态,有效避免了轴与轴承之间的摩擦. 当转子在恒定的高速下运行时,由于气膜收敛楔的存在,动压效应将取代静压效应起主导作用,能在提高轴承工作性能的同时降低系统的功耗[15].

总体而言,该轴承在燃料电池车载高速空压机的使用工况下,既可以实现在压缩机启停、变速时无磨损,又能在压缩机高速运转时,在有限的供气压力下表现出优良的承载力与稳定性,兼具了静压和动压气体轴承的优点.

1.2. 几何建模

本研究的三楔气体轴承几何参数如表1所示.

表 1   三楔气体轴承参数

Tab.1  Three-wedge gas bearing parameters

物理量变量符号数值
供气孔直径dl/mm2.0
节流孔直径d2/mm0.35
平均气膜间隙$\bar c $/mm0.02
轴承半径r/mm12.5
轴承长度l/mm34.5

新窗口打开| 下载CSV


1.3. 网格划分和无关性验证

三楔气体轴承的流场结构与实体结构分别如图12所示,轴承表面均匀分布着4排共36个节流孔. 气体经由这些节流微孔进入到轴与轴承之间的间隙,形成高压气膜并为转子提供支承.

图 2

图 2   三楔气体轴承实体结构

Fig.2   Entity structure of three-wedge gas bearing


利用CFD分析的前处理软件GAMBIT创建包括气膜、供气孔和节流孔在内的气体轴承流场的计算模型,然后将其划分为结构化的六面体单元. 为了选择最合适的网格密度,选择气膜的最大压力pmax和轴承承载力Fc为参数进行网格无关性分析. 如表2所示为不同网格数N的计算模型在相同条件下(偏心率ε = 0.5、转速n = 100 000 r/min)的模拟结果比较. 表中,$\Delta {{D_{\mathrm{P}}}} $$\Delta {{D_{\mathrm{F}}}} $分别为不同网格模型计算得出的最大压力pmax、承载力Fc与参考值的偏差百分比. 结果表明,使用42万和85万网格的计算模型得出的计算结果与参照值的偏差较大,而在网格数量超过114万后,偏差变小. 为了平衡计算精度、计算量和计算时间等,本研究采用具有 114 万个网格的计算模型(见图3).

表 2   网格无关性验证

Tab.2  Mesh independence verification

N/104pmax / PaΔDP / %Fc / NΔDF / %
42511416.64.145124.74.954
85519595.02.617128.32.210
114529237.90.810130.10.838
180532921.80.119130.90.229
250533559.2参考值131.2参考值

新窗口打开| 下载CSV


图 3

图 3   三楔气体轴承流场建模和网格划分

Fig.3   Modeling and meshing of three-wedge gas bearing fluid field


由于气体轴承在空间的3个方向上存在较大的尺寸差异(气膜厚度在几十微米以内),为了得到合理的网格长宽比,在保证计算量合理的前提下能高效收敛至与实际情况最接近的模拟计算结果,对计算模型气膜厚度方向的网格层数进行手动调整. 在压力和速度分布变化更显著的孔口区域周围对网格也进行了细化.

1.3.1. Fluent计算方法

在数值模拟中采用的是Fluent湍流计算模型中的RNG k-$\varepsilon $模型. 这是一种改进的k-$\varepsilon $模型,考虑了涡流的影响,有效改进了精度,也突破了原本的k-$\varepsilon $模型只能适用于高雷诺数模型的限制. 相关研究表明,较其他方法而言,该模型更适用于近壁区的层流及其他低雷诺数的流动计算[16-18],也能更好地处理流线弯曲程度较高或应变率较高的情况[19-20],因此在更广泛的流动中具有可信度和信赖度,能有效地模拟本研究涉及到的低雷诺数流动模型.

1.3.2. 边界条件和瞬态计算设置

在Fluent软件中对该网格模型进行模拟计算时,采用三维的压力基(pressure-based)隐式求解器PISO进行瞬态计算的求解. 流场计算使用进口、出口和壁面3个边界条件. 1)气体轴承的供气孔为压力进口边界条件,在计算中提供进口总压力;2)气体轴承两侧面端为压力出口边界条件,根据实际情况,压力设置为环境大气压力;3)气膜的内表面设置为随轴旋转的移动壁面,剩余边界设置为固定壁面,壁面之间不发生相对滑动.

1.3.3. 动网格的动态参数求解法

为了计算轴颈径向位移的过程中轴承气膜力的变化,流场的计算模型及其网格须随着轴颈中心的运动进行相应的变化和更新. 气体轴承的流场计算模型在空间的3个方向上尺寸的尺度存在较大差异,且在瞬态计算过程中轴位移的方向和大小在每个时间步都不固定。因此,在使用Fluent自带的3种动网格方法时,可能会出现网格质量下降、网格变形率大甚至负体积网格等问题,难以获得准确的计算结果.

为了避免网格变形,保证瞬态计算的进行,提出基于结构化网格的动网格移动方法. 该方法利用C语言编写了计算每个时间步中网格节点坐标的程序,并通过UDF接口将该动态网格函数加载到Fluent中. 在瞬态计算过程中,当轴颈中心发生移动时,程序将根据轴颈中心的位置计算出所有网格节点的新坐标,并将节点更新至相应位置[21].

当气体轴承流场模型采用全结构化网格划分后,取流场计算模型轴承面上的任意节点P1,其沿网格线投影到轴表面上的节点为P2. 当轴心移动时,在轴承壁面上的节点P1保持不动,而轴颈表面的节点P2移动到$ {P}'_{2} $,其新坐标如下:

$ \left.\begin{array}{l}x_{P_2^{\prime}}=x_{P_2}+x_{\mathrm{d}}, \\y_{P_2^{\prime}}=y_{P_2}+y_{\mathrm{d}}.\end{array}\right\} $

当轴心移动时,在所有计算模型的节点中,轴颈上节点(如P2)的位移距离最大,离轴心越远的节点位移越小. 对于P1P2的连接网格线上的任意节点P,其坐标可以计算如下:

$ \left.\begin{array}{l}x_{P^{\prime}}=x_P+\dfrac{N_i}{N} x_{\mathrm{d}} , \\y_{P^{\prime}}=y_P+\dfrac{N_i}{N} y_{\mathrm{d}} .\end{array}\right\} $

式中:Ni为节点P所在的径向网格层数,N为径向网格的总层数. 如图4所示,径向网格层数的编号从轴承表面(编号为0)开始到轴颈表面(编号为N).

图 4

图 4   动网格方法的计算原理图

Fig.4   Calculation principle of dynamic mesh method


可以看出,在配备该动网格程序后,在网格更新的过程中,网格节点的总数量及节点之间的拓扑关系保持不变. 即使经过多次的网格更新后,网格质量依旧良好. 该方法可以保证轴颈中心运动时,流场计算模型网格单元的光滑性和规律性,有效避免单元体积为负的情况. 如此,在轴颈扰动的瞬态计算过程中,可以通过网格更新后的模型求解得到轴颈在不同位置时的气膜力与轴承参数.

1.4. 基于动网格的刚度阻尼计算方法

刚度和阻尼系数反映了轴颈处于静平衡位置时,在外部载荷和小扰动的作用下,轴承气膜力与轴颈中心位移、速度变化的关系. 两者直接影响转子-轴承系统的稳定性、临界转速、谐响应等动态特性[22].

1.4.1. 理论计算方法

当气体轴承的轴颈中心处于静平衡位置时,在外界的位移或速度扰动下,轴颈受到的气膜力就会发生变化. 该气膜力的变化与外界扰动(位移或速度变化)呈非线性关系,但当扰动为微小量的情况下,为了简化分析可以将其视为线性关系. 因此,当轴在静平衡位置附近发生小扰动时,气膜力$ {\boldsymbol{F}} $可表示为关于轴颈中心坐标xy、轴颈中心的移动速度$ \dot{x} $$ \dot{y} $的函数:

$ \quad{\boldsymbol{F}}={{\boldsymbol{F}}}_{0}+\Delta {\boldsymbol{F}}={{\boldsymbol{F}}}_{0}+{\boldsymbol{k}}\left[\begin{array}{c}\Delta x\\ \Delta y\end{array}\right]+{\boldsymbol{c}}\left[\begin{array}{c}\Delta \dot{x}\\ \Delta \dot{y}\end{array}\right], $

$ {\boldsymbol{k}}=\left[\begin{array}{cc}{k}_{xx}& {k}_{xy}\\ {k}_{yx}& {k}_{yy}\end{array}\right]=\left[\begin{array}{cc}\dfrac{\partial {F}_{x}}{x}& \dfrac{\partial {F}_{y}}{x}\\ \dfrac{\partial {F}_{x}}{y}& \dfrac{\partial {F}_{y}}{y}\end{array}\right], $

$ {\boldsymbol{c}}=\left[\begin{array}{cc}{c}_{xx}& {c}_{xy}\\ {c}_{yx}& {c}_{yy}\end{array}\right]=\left[\begin{array}{cc}\dfrac{\partial {F}_{x}}{\dot{x}}& \dfrac{\partial {F}_{y}}{\dot{x}}\\ \dfrac{\partial {F}_{x}}{\dot{y}}& \dfrac{\partial {F}_{y}}{\dot{y}}\end{array}\right]. $

式中:$ {{\boldsymbol{F}}}_{0} $为静平衡位置上的气膜力;kc为刚度系数和阻尼系数矩阵,矩阵中包含了它们的主参数和交叉参数. 计算模型的坐标方向如图5所示.

图 5

图 5   轴承流场网格模型坐标方向

Fig.5   Coordinate direction of bearing flow field mesh model


1.4.2. 动网格CFD瞬态模拟求解方法

当轴颈匀速运动时,表示气膜力的方程(式(3))可以简化,气膜力的变化量$ \Delta {\boldsymbol{F}} $可在两正交方向xy上分解为

$ \left.\begin{array}{l}\Delta {{{F}}}_{x}={k}_{xx}\Delta x+{k}_{xy}\;\Delta y,\\ \Delta {{{F}}}_{y}={k}_{yx}\Delta x+{k}_{yy}\;\Delta y.\end{array}\right\} $

轴承刚度即可使用CFD模拟分析进行计算. 当轴颈沿x方向匀速移动时,沿y方向的位移$ \Delta y=0 $. 利用位移中产生的气膜力变化,可以计算出刚度系数矩阵中的两参数$ {k}_{xx}={\Delta {F}_{x}}/{\Delta x} $以及$ {k}_{yx}= {\Delta {F}_{y}}/{\Delta x} $. 同理,当轴颈沿y方向运动时,$ \Delta x=0 $,另两参数$ {k}_{yy}={\Delta {F}_{y}}/{\Delta y} $$ {k}_{xy}={\Delta {F}_{x}}/{\Delta y} $可以计算. 运用1.3.3节中提到的动网格方法,对于轴颈的每个微小位移,可在Fluent软件中通过对xy方向的网格压力分布积分得到相应瞬态时间步的气膜力. 由于研究的是轴承径向的动态参数,沿z方向的气膜力变化不在讨论范围之内.

为了计算阻尼系数,须再进行一次轴颈位移相同,但时间步长不同的瞬态模拟分析. 由于时间步长的改变,在此次分析中,轴颈会进行相同位移和方向的运动,但速度与前次模拟中不同. 根据式(5),当$ \Delta y=0 $时,阻尼系数矩阵中的$ {c}_{xx}= {\Delta {F}_{x}}/{\Delta \dot{x}} $$ {c}_{yx}={\Delta {F}_{y}}/{\Delta \dot{x}} $. 同理,当Δx=0时,cyy=$ {\Delta {F}_{y}}/{\Delta \dot{y}} $$ {c}_{xy}= {\Delta {F}_{x}}/{\Delta \dot{y}} $.

该方法已在微孔静压气体轴承的研究中进行过使用. 相应的实验验证表明,该方法能对轴承的非线性气膜力及动态参数进行瞬态模拟计算,且能比传统的稳态计算方法更准确地预测转子-轴承系统的临界转速[23].

2. 三楔气体轴承仿真结果分析

由1.1节中的理论分析可以看出,动静压耦合效应的直接决定因素是流场中气膜的收敛楔结构和气体随转子的转动. 因此本研究主要关注不同轴颈偏心率和转子转速下三楔气体轴承的性能参数变化.

在直接利用车载压缩机气路进行供气的前提下,气体轴承供气压力最高为0.25 MPa,因此,后续研究均在该条件下进行.

2.1. 轴承流场压力分布

图6所示为不同转速n三楔轴承气体流场圆周方向压力p1的分布. 图6(a)示意了流场计算模型中网格节点所在的周向角度位置$ {\theta }_{{\mathrm{node}}} $,由此可知,30°、150°和270°这3处为圆弧形楔槽的中心,是最小气膜间隙所在的区域. 在偏心率为0的情况下,此处气膜厚度数值等于半径间隙c*. 在取下半周气膜进行分析时,将周向气膜高压中心与轴承几何对称中心270°处的夹角称为高压中心偏转角α(见图6(b)).

图 6

图 6   不同转速下的气膜周向压力分布

Fig.6   Circumferential pressure distribution of gas bearing fluid film under different rotational speeds


图6可以看出,当转速为零时,轴承几乎没有受到动压效应的作用,几个小型的压力峰值代表了进气孔带来的静压节流作用,各进气孔之间的压力不同则是气膜厚度不同带来的流场静压力变化. 此时气膜压力关于高压中心对称分布,且最大压力不超过供气压力0.25 MPa. 随着转速升高,动静压耦合效应增强,可以明显观察到流场最大压力的上升. 在轴承转速增加至100 000 r/min的过程中,气膜压力峰值由0.2347 MPa逐渐增大到0.2734 MPa,流场最大压力超过了供气压力,是动压效应才能带来的效果. 其次,在流场中出现了高压中心的偏心现象. 在低转速(小于20 000 r/min)范围内,轴承的静压效应起主导作用,旋转的转子受到气体的黏性阻力而产生偏转. 在转速由0加速到20 000 r/min的过程中,流场的高压中心偏转角α由0°增加至30.33°. 而后,转子转速持续升高产生的动压效应使轴承的动静压耦合效应增强,气体的可压缩性提高且其带来的影响逐渐超过其他因素,因此在转速由20 000 r/min增高到100 000 r/min的过程中,偏转角α从30.33°递减到24.72°. 高压中心偏转角α的变化是轴颈中心静平衡位置的变化带来的,它会造成轴承流场形状和周向气膜厚度的变化,进而导致轴承的楔形效应增加,收敛楔与发散楔间压力差值增大,动压效应增强,最终对轴承的承载力和动态参数产生影响.

图7所示为流场气体密度ρ的分布. 结合图67可以看出,当气体向楔槽中心流动时,由于被压缩而密度增大. 当转速升高时,气体压缩性增强,进一步解释了气膜压力峰值和轴承承载力随转速增加而增大的原因. 气体在被压缩的过程中体现出的压力和密度的变化率也随转速升高而增加.

图 7

图 7   不同转速下气膜周向密度分布

Fig.7   Circumferential density distribution of gas bearing fluid film under different rotational speeds


图8所示为不同转速下流场节流孔附近的周向压力p2分布. 由于轴承微孔带来的节流静压效应,气体在节流孔处达到压力峰值. 在气体流出节流孔并进入楔形间隙的过程中,压力由峰值逐渐减小. 转速越高,节流孔周围由轴承楔形效应带来的动压现象表现得越强.

图 8

图 8   不同转速下节流孔区域周向压力分布

Fig.8   Circumferential pressure distribution around gas bearing orifices under different rotational speeds


从气膜的平铺压力云图(见图9)可以更直观地看到,在低转速下,压力峰值出现在节流孔和3个楔槽的中心处. 随着转速增加,轴颈中心的静平衡位置开始变化,高压区中心从楔槽中心逐渐偏移. 此时,随着楔形间隙内的气体与转子一起高速旋转,气体的可压缩性增大并能提供更高的气膜压力,可观察到云图中的压力峰值从0.2500 MPa增加到0.2745 MPa,从与供气压力相等增加至超越供气压力. 由于动压效应随着转速的增加逐渐增强,节流孔周围由于静压效应产生的压力变化区域逐渐变小,压力分布也由随气膜厚度变化变得更向高压中心集中,说明动压逐渐替代了静压效应成为影响流场压力分布的主导因素.

图 9

图 9   不同转速下的气膜压力云图

Fig.9   Gas film pressure nephogram under different rotational speeds


2.2. 轴承动态参数分析研究

2.1节中研究的是偏心率为0的三楔轴承在转速改变时的动静压耦合强度变化. 当三楔轴承所在的转子系统受到外部载荷作用时,轴颈中心的偏心率会持续增加直至稳定于静平衡位置,轴承流场的周向压力分布情况也会由图6(b)变化至如图10所示的情况. 可以看到,当偏心率增加时,气膜厚度会发生改变并引起周向压力的不对称分布,引起轴承承载力与刚度的上升. 为了进一步探究三楔轴承的动态特性变化,对不同转速与偏心率下的轴承进行了瞬态模拟计算.

图 10

图 10   偏心率为0.3时的气膜周向压力分布

Fig.10   Circumferential pressure distribution of gas bearing fluid film under eccentricity of 0.3


由于外加载荷方向的动态参数是轴承性能主要关注的研究对象,因此后文中轴承动态参数的主刚度、阻尼方向为yy,交叉刚度、阻尼方向为xy.

图11所示为不同转速和偏心率下的轴承承载力. 三楔气体轴承的承载力Fc随偏心率ε和转速n的上升而增加. 其中,承载力关于转速线性增大,而关于偏心率呈指数式快速增大,特别是在偏心率大于0.4以及转速大于40 000 r/min之后.

图 11

图 11   不同转速和偏心率下的轴承承载力

Fig.11   Bearing capacity of gas bearing under different rotational speeds and eccentricities


不同偏心率和转速对轴承动态系数的影响如图12所示. 随着偏心率和转速的增加,动压效应变强,主刚度kyy、交叉刚度kxy系数均有增大的趋势,且在偏心率大于0.4后变化更显著. 在大于40 000 r/min的高转速下,交叉刚度系数受转速影响较小,变化趋势相近.

图 12

图 12   不同转速和偏心率下的轴承刚度系数

Fig.12   Stiffness coefficients of gas bearing under different rotational speeds and eccentricities


阻尼系数表征了气膜力随速度的变化率,气膜阻尼的存在对旋涡能量消耗有影响,是抑制旋涡运动的一个因素. 如图13所示为不同转速和偏心率下的轴承阻尼系数. 可以看出,随着转速的增加,转子更容易克服气膜力进行运动,阻尼系数有减小的趋势. 随着偏心率的增大,动压效应增强,导致转子必须施加更多的功来克服气膜力. 因此可看到在偏心率ε< 0.4时,随着偏心率的升高,轴承阻尼系数有增大的趋势,轴承对涡动的抑制加强,转子系统运行更稳定. 但在偏心率大于0.4后,阻尼的主系数cyy与交叉系数cxy均表现出下降的趋势. 说明偏心率的持续增加导致了轴承气膜过薄,无法再提供有效支承.

图 13

图 13   不同转速和偏心率下的轴承阻尼系数

Fig.13   Damping coefficients of gas bearing under different rotational speeds and eccentricities


结合图1213可以发现,当转速为零时,交叉刚度和交叉阻尼系数均接近于零. 同时,刚度与阻尼的交叉系数数值都远小于主系数,这是多楔轴承特有的一种动态特性.

2.3. 与圆柱气体轴承的性能对比

在已知外接气源的圆柱轴承能为转速高达100 000 r/min的氢燃料电池汽车车载空压机转子提供足够支承[23]的情况下,将0.25 MPa供气压力下的三楔气体轴承与传统的0.50 MPa供气压力下的静压圆柱气体轴承进行性能对比,进而确定三楔轴承是否能在仅使用内部气路的情况下满足车载空压机的支承需求. 2种轴承的结构及流场形状对比如图14所示,其中三楔气体轴承的设计平均气膜间隙(c*c的平均值)与圆柱轴承的气膜间隙cr相同.

图 14

图 14   轴承结构及流场形状对比图

Fig.14   Structure and fluid field comparison of gas bearings


图15所示为不同供气压力下,2种轴承在不同转速和偏心率下的轴承承载力对比. 可以看出,三楔轴承与圆柱轴承的承载力Fc随偏心率的增加均呈整体上升的趋势. 在低转速(小于60 000 r/min)与小偏心率(ε<0.4)范围内,三楔轴承与圆柱轴承能达到相似的支承效果. 但在偏心率ε>0.4之后,通过分析曲线斜率可知,三楔轴承承载力的增加速度远超圆柱轴承的,这主要是由于三楔轴承动压效应的增强趋势远大于圆柱轴承的. 这种偏心率达到某个数值之后,承载力会出现非线性陡增的现象,这是动压结构的轴承特有的属性.

图 15

图 15   不同转速和偏心率下三楔与圆柱气体轴承承载力对比

Fig.15   Comparison of load capacity between three-wedge and cylindrical gas bearings under different rotational speeds and eccentricities


从对轴承承载力和动态参数的理论研究中可以看出,轴承的实际工作性能与轴颈中心的静平衡位置密切相关. 为了确定三楔轴承在指定供气压力下能否代替传统的圆柱轴承,须研究2种轴承在实际工况下的平衡位置,并对比相应位置处的动态性能.

在实际使用中,单个轴承须承载的车载离心空压机转子外载荷约为20 N[24]. 在此条件下,计算出不同转速下轴承静平衡位置,如图16所示. 图中,①~⑤分别表示转速为10000080000600004000020000 r/min. 由图可知,当轴承转速在工作范围内变化时,三楔轴承在相同转速下的偏心率始终小于圆柱轴承的.

图 16

图 16   不同转速下三楔与圆柱气体轴承在车载空压机转子上的静平衡位置对比

Fig.16   Comparison of static equilibrium positions between three-wedge and cylindrical gas bearings on FCV compressor rotors under different rotational speeds


基于图16的计算结果,利用1.4节中的动网格瞬态计算方法,可求解出2种轴承在各平衡位置时的动态参数值,如图17所示. 由于每个动态参数计算点对应的轴承转速与偏心率不同,轴承的动态参数不会呈现单调变化. 对比2种轴承的动态参数可知,在低转速下,三楔气体轴承在支承性能方面表现出与传统圆柱轴承相似的静压效应,但当转速升高时,三楔轴承的动静压耦合特征能使其在高速运转时能表现出更好的支承特性和稳定性. 因此可以推论出,在使用燃料电池汽车空压机自身气路作为供气来源的前提下,新设计的三楔轴承能在0.25 MPa的供气压力下满足车载空压机的使用需求.

图 17

图 17   不同转速下三楔与圆柱气体轴承在车载空压机转子上的动态参数对比

Fig.17   Comparison of hydrodynamic characteristics between three-wedge and cylindrical gas bearings on compressor rotors in fuel cell vehicle under different rotational speeds


3. 结 论

提出适用于氢燃料电池汽车车载高速离心空压机的新型三楔气体轴承,以克服传统的静压圆柱轴承需要额外供气的缺点. 经由对三楔轴承的动、静压特性分析可知,此类轴承能在车载压缩机启停、变速、稳定运行中始终表现出良好的支承特性. 主要内容如下:

(1)提出基于自编UDF动网格程序的CFD瞬态模拟方法,用以计算轴颈在任意偏心位置时的轴承非线性气膜力. 使用该方法研究三楔气体轴承的流场特性、动静压耦合效应以及静、动态参数的变化,并与圆柱气体轴承的各项参数对比.

(2)气体黏性和压缩性均会影响三楔气体轴承的动静压耦合效应强度. 三楔轴承流场中气体的黏性力和小孔节流的静压效应在转子系统低转速(小于20 000 r/min)时起主导作用;当转速较高(大于20 000 r/min)时,气体的压缩性变化对轴承流场动静压耦合效应强度的影响逐渐增大,动压效应占主导地位.

(3)将传统使用的0.50 MPa供气压力下的圆柱气体轴承与使用车载空压机内部气路的0.25 MPa供气压力下的三楔气体轴承进行性能对比. 结果表明,相较于圆柱轴承,动压效应对三楔轴承特性的增益更显著.

通过对2类轴承实际使用中静平衡位置及相应位置的动态参数计算,可以观察到三楔轴承在低转速与下能表现出与传统圆柱轴承相似的静压效应,而当转速升高时,三楔轴承的动静压耦合特征能使其在更低的供气压力下表现出良好的支承特性和稳定性. 较传统的圆柱轴承,三楔气体轴承因其无须外接气源的优点,能降低燃料电池汽车系统内气体回路的复杂度,同时其也能适应车载空压机需要频繁启停、变速的工作环境.

本研究已经取得了初步的成果,在后续工作中,可设计加工氢燃料电池汽车车载空压机样机,通过实验与理论分析相结合的方式,进一步证明三楔动压气体轴承设计的有效性,以及仅使用空压机内部气路为轴承供气这一研究方向的合理性.

参考文献

付佩, 周紫佳, 兰利波, 等

氢燃料电池汽车发动机关键技术研究现状及趋势展望

[J]. 汽车工程学报, 2022, 12 (4): 11

DOI:10.3969/j.issn.2095-1469.2022.04.05      [本文引用: 1]

FU Pei, ZHOU Zijia, LAN Libo, et al

Research and development trends of key technologies of hydrogen fuel cell engines

[J]. Chinese Journal of Automotive Engineering, 2022, 12 (4): 11

DOI:10.3969/j.issn.2095-1469.2022.04.05      [本文引用: 1]

任天明. 高速水润滑轴承电动离心式空压机关键技术研究[D]. 北京: 北京科技大学, 2017.

[本文引用: 1]

REN Tianming. Study on the high speed motorized centrifugal air compressor with water lubricated journal bearings [D]. Beijing: University of Science and Technology Beijing, 2017.

[本文引用: 1]

徐文章. 双螺杆空压机的水润滑动静压轴承-转子动力学分析[D]. 南京: 东南大学, 2019.

XU Wenzhang. Dynamic analysis of the water-lubricated hybrid bearing-rotor system of twin-screw air compressor [D]. Nanjing: Southeast University, 2019.

王智洋, 张庆源. 一种基于磁悬浮轴承的车载燃料电池气体压缩机: CN110886706A[P]. 2021−03−19.

[本文引用: 1]

WAN Y, GUAN J, XU S

Improved empirical parameters design method for centrifugal compressor in PEM fuel cell vehicle application

[J]. International Journal of Hydrogen Energy, 2016, 42 (8): 5590- 5605

[本文引用: 1]

DAEJONG K

Foil bearing technology for high speed single stage air compressors for fuel cell applications

[J]. Chinese Journal of Turbomachinery, 2021, 63 (4): 54- 59

[本文引用: 1]

ZHANG H, YU W, HUA W

Design and key technology of oil-free centrifugal air compressor for hydrogen fuel cell

[J]. China Electrotechnical Society Transactions on Electrical Machines and Systems, 2022, 6 (1): 11- 19

DOI:10.30941/CESTEMS.2022.00003      [本文引用: 1]

张雯, 董泽达, 徐方程

氢燃料电池空压机动压箔片轴承设计与试验

[J]. 兵器装备工程学报, 2022, 43 (4): 172- 178

DOI:10.11809/bqzbgcxb2022.04.028      [本文引用: 1]

ZHANG Wen, DONG Zeda, XU Fangcheng

Design and test of aerodynamic foil bearings for hydrogen fuel cell air compressors

[J]. Journal of Ordnance Equipment Engineering, 2022, 43 (4): 172- 178

DOI:10.11809/bqzbgcxb2022.04.028      [本文引用: 1]

冯凯, 黄明, 李成勤, 等

径向气体箔片轴承高速重载测试及实验台搭建

[J]. 湖南大学学报: 自然科学版, 2020, 47 (2): 1- 5

[本文引用: 1]

FENG Kai, HUANG Ming, LI Chengqin, et al

Test of high speed and heavy load gas foil bearing and set up of test rig

[J]. Journal of Hunan University: Natural Sciences, 2020, 47 (2): 1- 5

[本文引用: 1]

王伟, 李晓疆, 曾强, 等

高速透平机械全金属鼓泡箔片动压气体轴承稳定性研究

[J]. 西安交通大学学报, 2017, 51 (8): 84- 89

[本文引用: 1]

WANG Wei, LI Xiaojiang, ZENG Qiang, et al

Stability analysis for fully hydrodynamic gas-lubricated protuberant foil bearings in high speed turbomachinery

[J]. Journal of Xi'an Jiaotong University, 2017, 51 (8): 84- 89

[本文引用: 1]

靳彩妍, 徐方程, 吴孟龙, 等

空气轴承支撑的燃料电池离心空压机转子动力学分析

[J]. 汽车零部件, 2020, (11): 6- 9

[本文引用: 1]

JIN Caiyan, XU Fangcheng, WU Menglong, et al

Calculation and analysis of air bearing centrifugal compressor rotor dynamics for fuel cells

[J]. Automobile Parts, 2020, (11): 6- 9

[本文引用: 1]

彭万欢, 郑越青, 徐刚

空气动压箔片轴承启停性能的实验研究

[J]. 润滑与密封, 2011, 36 (7): 52- 55

DOI:10.3969/j.issn.0254-0150.2011.07.013      [本文引用: 1]

PENG Wanhuan, ZHENG Yueqing, XU Gang

Experimental study on the start-stop characteristic of foil air bearing

[J]. Lubrication Engineering, 2011, 36 (7): 52- 55

DOI:10.3969/j.issn.0254-0150.2011.07.013      [本文引用: 1]

王云飞. 气体润滑理论与气体轴承设计[M]. 北京: 机械工业出版社, 1999: 29–31.

[本文引用: 1]

十合晋一. 气体轴承设计指南[M]. 东京: 共立出版株式会社, 2002: 22.

[本文引用: 1]

刘暾, 刘育华, 陈世杰. 静压气体润滑[M]. 哈尔滨: 哈尔滨工业大学出版社, 1990.

[本文引用: 1]

李鑫. 高速艇喷水推进器流场特性数值模拟[D]. 镇江: 江苏科技大学, 2014.

[本文引用: 1]

LI Xin. Numerical simulation of flow field characteristics of the high-speed craft waterjet propulsion [D]. Zhenjiang: Jiangsu University of Science and Technology, 2014.

[本文引用: 1]

陈伟威

基于RNG k-$\varepsilon $紊流模型的溢洪道水力特性数值模拟

[J]. 水利科学与寒区工程, 2023, 6 (1): 10- 13

DOI:10.3969/j.issn.2096-5419.2023.01.004     

CHEN Weiwei

Numerical simulation of hydraulic characteristics of spillway based on RNG k-$\varepsilon $ turbulence model

[J]. Hydro Science and Cold Zone Engineering, 2023, 6 (1): 10- 13

DOI:10.3969/j.issn.2096-5419.2023.01.004     

郭晓娟. 脊状表面减阻特性及其结构优化设计研究[D]. 西安: 西北工业大学, 2007.

[本文引用: 1]

GUO Xiaojuan. Research on drag-reduction characteristics of groove surface and its structural optimization design [D]. Xi’an: Northwestern Polytechnical University, 2007.

[本文引用: 1]

杜金威, 王均星, 李辉成

阳升观台阶溢洪道水流数值模拟研究

[J]. 武汉大学学报: 工学版, 2020, 53 (2): 110- 116

[本文引用: 1]

DU Jinwei, WANG Junxing, LI Huicheng

Numerical simulation of the water flow in the stepped spillway of Yangshengguan reservoir

[J]. Engineering Journal of Wuhan University, 2020, 53 (2): 110- 116

[本文引用: 1]

王福军. 计算流体动力学分析: CFD软件原理与应用[M]. 北京: 清华大学出版社, 2004: 124−125.

[本文引用: 1]

LI Q, LIU S, PAN X, , et al

A new method for studying the 3D transient flow of misaligned journal bearings in flexible rotor-bearing systems

[J]. Journal of Zhejiang University: Science A, 2012, 13 (4): 293- 310

[本文引用: 1]

HEINDEL S, MÜLLER P C, RINDERKNECHT S

Unbalance and resonance elimination with active bearings on general rotors

[J]. Journal of Sound and Vibration, 2018, (431): 422- 440

[本文引用: 1]

ZHANG C, HONG W, ZHENG S

The dynamic pressure effect of gas microbearings under high rotational speed

[J]. Engineering Applications of Computational Fluid Mechanics, 2022, 16: 1,1264- 1278

[本文引用: 2]

李宏钦. 氢燃料电池车用高速离心式压缩机的设计及研究[D]. 杭州: 浙江大学, 2019.

[本文引用: 1]

LI Hongqin. Design and research of highspeed centrifugal compressor for hydrogen fuel cell vehicle [D]. Hangzhou: Zhejiang University, 2019.

[本文引用: 1]

/