-
气体箔片推力轴承具有摩擦功耗低、稳定性好、耐高温等优点,特别适用于高速透平机械[1-2]。目前常用的的气体箔片推力轴承中,波箔型推力轴承的工程应用最为广泛,人们对其展开了一系列理论和实验研究。在其工作过程中波箔和顶箔都将发生弹性变形,并且该变形与气膜压力互相影响,为动态过程,导致轴承润滑性能的计算十分困难,成为实际研究和应用的一个难点。
1975年,Walowit等[3]首次对波箔型气体箔片轴承的结构进行了理论分析,通过建立箔片弹性模型计算其弹性系数,结合Reynolds方程,研究了气膜压力分布及轴承的一系列静态特性。Heshmat等[4]将弹性箔片气体推力轴承的波箔简化为一系列线性弹簧,采用有限差分法求解Reynolds方程,通过分析轴承几何结构参数对其静态特性的影响,对轴承结构参数进行了优化。Roger等[5]在考虑波箔与轴承座或顶箔之间摩擦力、波箔局部相互作用力及其几何形状的基础上,建立了一个新的箔片变形模型。结果表明,波箔的刚度从固定端到自由端递减,增大波箔和顶箔之间的摩擦系数可有效增大箔片结构刚度,且载荷分布可极大地影响波箔刚度。Iordanoff[6-7]借助弹性空气动力学3D建模,提出了一种针对波箔型气体箔片推力轴承的快速设计方法,根据所需润滑性能通过反向求解得到对应的轴承性能参数。Heshmat等[8]结合有限元法和有限差分法研究波箔型气体箔片推力轴承润滑特性,结果表明该计算方法收敛速度较快,且阐明了顶箔变形对轴承性能的显著影响。戚社苗等[9]通过联立求解压力控制Reynolds方程和箔片结构变形方程,得到了弹性箔片动压气体推力轴承的气弹耦合解,并计算分析了轴承的承载性能。闫佳佳[10]将波箔型动压气体推力轴承的箔片结构简化为线性弹簧,理论研究了其静动特性,并建立了推力轴承-转子系统有限元模型,分析了其动力学特性。Park等[11]针对气体箔片推力轴承推力盘的倾斜情况,理论研究了不同倾斜角和倾斜方向对轴承静态特性和动态特性的影响规律。Lee等[12]利用有限差分法求解柱坐标下的二维薄板方程,以预测混合型气体箔片推力轴承顶箔的变形,并研究了不同结构参数对轴承静动态特性的影响。
文献研究表明,准确计算箔片结构的弹性变形是研究气体箔片推力轴承润滑性能的关键。但过去的研究中一般以线性弹簧模型描述波箔的支承作用,并忽略顶箔在气膜压力作用下的变形,这种近似研究方法的精度有限。实际工况中,波箔的变形十分复杂,且顶箔在波箔支承力和气膜压力共同作用下,将产生不可忽略的变形,采用简化的方法对箔片结构的变形进行描述无法达到较高的精度。为此,本文利用弹性力学理论分别对波箔和顶箔的力学性能进行研究,构建其受力变形模型,结合Reynold方程,构建一种流固耦合算法,采用数值方法研究了轴承结构参数对波箔型气体箔片推力轴承润滑性能的影响。
-
图1示出了波箔型气体箔片推力轴承结构示意图,轴承主要由轴承座、波箔、顶箔及推力盘组成,顶箔位于波箔上方,波箔为顶箔提供弹性支承,两箔片的一端共同固定于轴承座上,由此构成若干扇形瓦块。扇形瓦块结构参数如图2(a)所示,其中
${R_1}$ 为轴承内径,${R_2}$ 为轴承外径,$\beta $ 为瓦块张角,b为节距比($b = \gamma /\beta $ )。图2(b)为扇形瓦块剖面图,推力盘转速较高时,推力盘将气体由进气口带入扇形瓦块区域,然后通过出气口排出,初始气膜间隙由${h_1}$ 减小为${h_2}$ ,顶箔和推力盘之间形成楔形润滑气膜,为推力盘提供支承力,实现气体动压润滑。 -
波箔在工作中的变形如图3所示,气膜压力通过顶箔的传递,间接作用于波箔,使波箔发生弹性变形。由于波箔结构在轴承圆周方向上均匀分布,因此对单波进行的受力变形分析适用于整块波箔。
为便于计算分析,作如下假设:(1) 波箔的顶部和顶箔之间始终保持接触,单波所受载荷为作用于波峰的集中力;(2) 由于波箔和顶箔、波箔和轴承座之间的摩擦很小,故可忽略;(3) 波箔受力后的弹性变形量远小于原始尺寸,属于小变形;(4) 如图3(b)所示,由于单波周向变形量很小,故可认为波箔变形后仍保持为一个扇形,且内外径不变。
波箔的单波在工作过程中的变形如图4所示,单波横截面尺寸为
${t_{\rm{B}}} \times L$ (L为波箔内外半径之差),波拱半角为$\varphi $ ,波拱半长为l,其顶部受到竖直向下的载荷F,由于其横截面高度${t_{\rm{B}}}$ 远小于其轴线半径${R_{\rm{0}}}$ ,因此可利用卡氏定理计算其变形。相较于单波顶点的竖直位移$\Delta y$ ,其水平位移$\Delta x$ 对气膜厚度的影响非常小,故本文不予考虑。下面计算$\Delta y$ ,若规定使单波曲率增大的弯矩为正,则单波任意截面m-m上的弯矩M为:M关于F的偏导数为:
对单波左半段进行积分,并将式(1)、(2)代入可得单波顶点的竖直位移为:
式中:E为箔片弹性模量。
单波横截面惯性矩可表示为
${I_{\rm{B}}} = \dfrac{{L{t_B}^3}}{{12}}$ ,代入式(3)可得: -
轴承在工作时,顶箔在气膜压力和波箔支承力的作用下将发生一定的弹性变形。为计算该弹性变形量,将顶箔视为二维梁进行分析。如图5所示,选取位于两单波顶点之间的某一段顶箔作为分析对象,由于其长度s(即波距)较小,故可将该区间内的气膜压力视为均布载荷q。
变形方程可表示为:
式中:
${w_1}$ 为顶箔竖直变形量,x为距左顶点的距离,${I_{\rm{T}}}$ 为顶箔横截面惯性矩。对式(5)积分可得:
暂时仅考虑顶箔变形,假设波箔不变形,可得以下边界条件:
将式(7)代入式(6),可求得常数项,得到:
-
箔片结构总体变形由顶箔和波箔共同决定,其值可表示为:
式中:
${w_1}$ 为顶箔自身产生的挠度,可由式(8)求得,${w_2}$ 为顶箔在波箔变形影响下所产生的竖直位移。如图6所示,波箔两个顶点的竖直位移分别为${w_{{\rm{B1}}}}$ 和${w_{{\rm{B2}}}}$ ,水平位移分别为$\Delta {x_1}$ 和$\Delta {x_2}$ ,则由几何关系可得,距左顶点x处的位移${w_2}$ 为: -
联立可压缩气体的运动方程、连续性方程和能量守恒方程,并结合理想气体状态方程,可得用来描述顶箔与推力盘之间气体流动规律的等温定常Reynolds方程:
式中:r为径向坐标,
$\theta $ 为周向坐标,p为气膜压力,h为气膜厚度,$\mu $ 为气体的动力黏度,$\omega $ 为转轴角速度。将式(11)量纲为-化得:
式中:
$\bar r = \dfrac{r}{{{R_2}}}$ ,$P = \dfrac{p}{{{p_0}}}$ ,$H = \dfrac{h}{{{\delta _h}}}$ ,$\Lambda = \dfrac{{6\mu \omega {R_2}^2}}{{{p_0}{\delta _h}^2}}$ ,其中${p_0}$ 为环境压力,${\delta _{\rm{h}}}$ 为进出口气膜厚度之差,$\Lambda $ 为轴承系数。根据图2的几何关系,可得推力盘无倾斜时的气膜厚度表达式:
式中:
${w_t}(r,\theta )$ 为箔片结构的总体变形量,$g(\theta )$ 可表示为:故量纲为-气膜厚度表达式为:
-
式(12)所示的Reynolds方程为非线性偏微分方程,本文采用有限差分法及Newton-Raphson迭代法进行求解[13]。当轴承稳定工作时,其每个扇形瓦块内的气膜压力分布及气膜厚度分布相同,故取任一个瓦块作为计算域即可。如图7所示,轴承气膜承载区以波箔固定端和自由端为周向边界,以轴承内外径处为径向边界。将其离散为
$m \times n$ 的网格,压力值${P_{i,j}}$ 表示各节点上气膜压力值,每格步长为:$\Delta \theta = \dfrac{\beta }{m}$ ,$\Delta \bar r = \dfrac{{{R_2} - {R_1}}}{{n{R_2}}}$ 。为提高计算精度,采用中差商表示偏微分项:
将式(16)~(18)代入式(12),化简得:
式中:
${A_{i,j}} = \dfrac{{{{\left( {P{H^3}\bar r} \right)}_{i,j + 1/2}}}}{{{{\bar r}_{i,j}}\Delta {{\bar r}^2}}}$ ,${B_{i,j}} = \dfrac{{{{\left( {P{H^3}\bar r} \right)}_{i,j - 1/2}}}}{{{{\bar r}_{i,j}}\Delta {{\bar r}^2}}}$ ,${C_{i,j}} = \dfrac{{{{\left( {P{H^3}} \right)}_{i + 1/2,j}}}}{{{{\bar r}_{i,j}}^2\Delta {\theta ^2}}}$ ,${D_{i,j}} = \dfrac{{{{\left( {P{H^3}} \right)}_{i - 1/2,j}}}}{{{{\bar r}_{i,j}}^2\Delta {\theta ^2}}}$ ,${E_{i,j}} = {A_{i,j}} + {B_{i,j}} + {C_{i,j}} + {D_{i,j}}$ ,${F_{i,j}} = \Lambda \dfrac{{{{\left( {PH} \right)}_{i + 1/2,j}} - {{\left( {PH} \right)}_{i - 1/2,j}}}}{{\Delta \theta }}$ 求解过程中引入边界条件,可压缩气体润滑中不存在气膜破裂现象,轴承承载区域各边界均与大气环境连通,故其压力与环境压力相同,无量纲形式的边界条件为:
$P\left( {0,\bar r} \right) = 1{\rm{ }},P\left( {\beta ,\bar r} \right) = 1{\rm{ }},P\left( {\theta ,\dfrac{{{R_1}}}{{{R_2}}}} \right) = 1{\rm{ }}, P\left( {\theta ,1} \right) = 1$ 。引入下列相对收敛准则:
图8示出了本文流固耦合的数值计算过程。
-
为验证上述数值解法的正确性,选用与文献[4]相同的刚性气体推力轴承参数,按上述算法借助MATLAB编程计算其单个瓦块的承载力,计算结果对比如表1所示。由对比结果可知,对于不同气膜间隙比和轴承数下的单个瓦块承载力,本文计算所得结果与参考文献结果的误差均小于3%,具有较高的吻合度。
$\frac{{{h_1}}}{{{h_2}}}$ $\Lambda $ Dimensionless bearing capacity $\bar W \times {10^2}$ Deviation/% Results in this paper Results from literature 2 0.1 0.0137 0.014 −2.14 1 0.1402 0.144 −2.64 10 1.6000 1.612 −0.74 20 3.3276 3.33 −0.07 40 6.2657 6.246 0.32 5 0.1 0.1866 0.192 −0.03 1 2.3580 2.343 0.64 10 22.7503 22.2 2.48 20 35.5801 34.66 2.65 40 47.7257 46.4 2.86 表 1 计算结果对比
Table 1. Comparison of calculation results
-
在刚性气体推力轴承的基础上,以某一个具体波箔型气体箔片推力轴承为算例进行研究,计算时取
$m = 60$ ,$n = 50$ 。轴承瓦块数目为6,各瓦块波纹数目均为12,其他轴承结构及运行参数如表2所示[14]。图9至图11分别示出了该波箔型气体箔片推力轴承的量纲为-气膜压力分布图、量纲为-气膜厚度分布图及箔片结构变形量分布图。图12示出了箔片推力轴承和刚性推力轴承在扇形瓦块径向中点处(即
$\overline r = 0.5$ )的量纲为-气膜厚度对比图及其局部放大图。由图12(a)可知,两者气膜厚度的大小相近,且沿周向的变化趋势相似;由图12(b)可知,箔片推力轴承的气膜厚度略大于刚性推力轴承,且有小幅波动。图13示出了不同节距比下箔片推力轴承和刚性推力轴承的承载力及摩擦力矩对比图。由图13(a)可知,箔片推力轴承的承载力略小于刚性推力轴承,且当节距比b在0.6附近时,两者承载力均达到最大值。图13(b)可知,箔片推力轴承的摩擦力矩略小于刚性推力轴承,两者均随节距比b的增大而减小,且近似于线性变化。${R_1}$ /m${R_2}$ /m$\beta $ /(°)b ${t_{\rm{B}}}$ /m${t_{\rm{T} }}$ /m${R_0}$ /m$\varphi $ /rad$2.54 \times {10^{ - 2}}$ $5.08 \times {10^{ - 2}}$ 60 0.5 $1.016 \times {10^{ - 4}}$ $1.524 \times {10^{ - 4}}$ ${\rm{1}}{\rm{.2125}} \times {10^{ - 3}}$ 0.8364 l/m E/Pa μ/(Pa·s) p0/Pa ${h_1}$ /m${h_2}$ /mn/ $(r \cdot {\min ^{ - 1}})$ ${\rm{9}} \times {10^{ - {\rm{4}}}}$ $2.14 \times {10^{11}}$ $1.932 \times {10^{ - 5}}$ $1.01325 \times {10^5}$ $5.2 \times {10^{ - 5}}$ $2.6 \times {10^{ - 5}}$ 30000 表 2 轴承结构及运行参数
Table 2. Parameters of bearing structure and operation
图 12 量纲为-气膜厚度对比及其局部放大图
Figure 12. Comparison of dimensionless gas film thickness and its partial enlargement drawing
图 13 箔片轴承与刚性轴承润滑性能对比
Figure 13. Comparison of lubrication performance of foil bearing and that of rigid bearing
图14至图17为其他参数均不变时,各节距比b下波箔厚度
${t_{\rm{B}} }$ 、波拱半长l、波拱半径${R_0}$ 和顶箔厚度${t_{\rm{T}}}$ 对轴承润滑性能的影响关系图。由图可知:轴承的承载力和摩擦力矩先随波箔厚度增大而快速增大,当波箔厚度达到一定值后,两者增大速度均趋于平缓;承载力和摩擦力矩均随波拱半长增大而缓慢减小;波拱半径较小时,承载力和摩擦力矩受其变化影响较小,波拱半径较大时,两者均随其增大而减小;随着顶箔厚度的增大,承载力和摩擦力矩先快速增大,之后几乎不再变化。因此,可通过增大波箔厚度${t_{\rm{B}}}$ 、顶箔厚度${t_{\rm{T}}}$ 或减小波拱半长l、波拱半径${R_0}$ ,来提高波箔型气体箔片推力轴承的承载力,反之可降低轴承摩擦力矩。图 14 不同节距比下波箔厚度对润滑性能的影响
Figure 14. Influences of bump foil thickness on lubrication performance with different pitch ratios
图 15 不同节距比下波拱半长对润滑性能的影响
Figure 15. Influences of half bump length on lubrication performance with different pitch ratios
-
(1) 以波箔型气体箔片推力轴承为研究对象,运用弹性力学理论构建了波箔和顶箔的受力变形模型,得到了考虑箔片弹性变形的气膜厚度方程,结合可压缩气体Reynold方程,构建了一种流固耦合算法,并利用有限差分法和Newton-Raphson迭代法进行求解,借助MATLAB编程实现了对轴承润滑性能的精确数值计算。结合具体文献对所建立的计算模型进行了对比验证,所得结果与文献结果有良好的吻合度,证明了本文算法的可靠性。
(2) 通过对具体算例的研究,得到了轴承的量纲为-气膜压力分布、量纲为-气膜厚度分布与箔片结构变形量分布。在此基础上对比研究了箔片推力轴承和刚性推力轴承的润滑性能,结果表明,箔片轴承的气膜厚度较大,承载力和摩擦力矩较小。结合算例进一步研究了波箔厚度
${t_{\rm{B}}}$ 、波拱半长l、波拱半径${R_0}$ 和顶箔厚度${t_{\rm{T}}}$ 对轴承润滑性能的影响规律,得出具体的影响规律曲线图。结果表明:轴承的承载力和摩擦力矩随${t_{\rm{B}}}$ 或${t_{\rm{T}}}$ 的增大而增大,随l或${R_0}$ 的增大而减小。故较大的${t_{\rm{B}}}$ 、${t_{\rm{T}}}$ 或较小的l、${R_0}$ 可提高轴承的承载力,反之可减小轴承的摩擦力矩,有利于降低摩擦功耗,减少发热量。
波箔型气体箔片推力轴承润滑性能的数值研究
Numerical Study on the Lubrication Performance of Bump-type Gas Foil Thrust Bearing
-
摘要: 以波箔型气体箔片推力轴承为研究对象,运用弹性力学理论构建了波箔和顶箔的受力变形模型,计算箔片结构变形量,得到气膜厚度方程。在此基础上,结合可压缩气体Reynold方程,构建了一种能够对波箔型气体箔片推力轴承润滑性能进行计算的流固耦合模型。结合具体的文献对所建立的流固耦合计算模型进行了验证,证明了算法的可靠性。通过具体的算例,用数值方法研究了轴承的量纲为-气膜压力分布、量纲为-气膜厚度分布与箔片结构变形量分布。对比研究了箔片推力轴承和刚性推力轴承的润滑性能,发现箔片轴承的气膜厚度较大,承载力和摩擦力矩较小。深入研究了波箔厚度tB、波拱半长l、波拱半径R0和顶箔厚度tT对轴承润滑性能的影响规律,得到了相关曲线图,结果表明:轴承的承载力和摩擦力矩随tB或tT的增大而增大,随l或R0的增大而减小。Abstract: Taking bump-type gas foil thrust bearing as the research object, a force and deformation model of bump foil and top foil were established by elastic mechanics analysis. The foil structure deformation was calculated and the gas film thickness equation was obtained. On this basis, a fluid-solid coupling model used for calculating the lubrication performance of bump-type gas foil thrust bearings was established with use of Reynolds equation for compressible gas. The established fluid-solid coupling calculation model was verified with specific literature, proving the reliability of the algorithm. With a specific example, the dimensionless gas film pressure distribution, dimensionless gas film thickness distribution and foil structure deformation distribution of the bearing were calculated numerically. The lubrication performance of foil thrust bearing and that of rigid thrust bearing were compared. It was found that the gas film thickness of foil bearing was bigger, while the bearing capacity and friction torque of it were smaller. The influences of bump foil thickness tB, half bump length l, bump radius R0 and top foil thickness tT on the lubrication performance of the bearing were studied, and relevant curves were obtained, which showed that bearing capacity and friction torque will increase as tB or tT increase, and decrease as l or R0 increase.
-
Key words:
- bump foil /
- thrust bearing /
- lubrication performance /
- fluid-solid coupling /
- numerical study
-
表 1 计算结果对比
Table 1. Comparison of calculation results
$\frac{{{h_1}}}{{{h_2}}}$ $\Lambda $ Dimensionless bearing capacity $\bar W \times {10^2}$ Deviation/% Results in this paper Results from literature 2 0.1 0.0137 0.014 −2.14 1 0.1402 0.144 −2.64 10 1.6000 1.612 −0.74 20 3.3276 3.33 −0.07 40 6.2657 6.246 0.32 5 0.1 0.1866 0.192 −0.03 1 2.3580 2.343 0.64 10 22.7503 22.2 2.48 20 35.5801 34.66 2.65 40 47.7257 46.4 2.86 表 2 轴承结构及运行参数
Table 2. Parameters of bearing structure and operation
${R_1}$ /m${R_2}$ /m$\beta $ /(°)b ${t_{\rm{B}}}$ /m${t_{\rm{T} }}$ /m${R_0}$ /m$\varphi $ /rad$2.54 \times {10^{ - 2}}$ $5.08 \times {10^{ - 2}}$ 60 0.5 $1.016 \times {10^{ - 4}}$ $1.524 \times {10^{ - 4}}$ ${\rm{1}}{\rm{.2125}} \times {10^{ - 3}}$ 0.8364 l/m E/Pa μ/(Pa·s) p0/Pa ${h_1}$ /m${h_2}$ /mn/ $(r \cdot {\min ^{ - 1}})$ ${\rm{9}} \times {10^{ - {\rm{4}}}}$ $2.14 \times {10^{11}}$ $1.932 \times {10^{ - 5}}$ $1.01325 \times {10^5}$ $5.2 \times {10^{ - 5}}$ $2.6 \times {10^{ - 5}}$ 30000 -
[1] GUO Z, FENG K, LIU T, et al. Nonlinear dynamic analysis of rigid rotor supported by gas foil bearings: Effects of gas film and foil structure on subsynchronous vibrations[J]. Mechanical Systems and Signal Processing, 2018, 107: 549-566. doi: 10.1016/j.ymssp.2018.02.005 [2] DELLACORTE C. 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 [3] WALOWIT J A, ANNO J N. Modern developments in lubrication mechanics[M]. London: Applied Science Publishers Ltd., 1975: 214-221. [4] HESHMAT H, WALOWIT J A, PINKUS O. Analysis of gas lubricated compliant thrust bearings[J]. Journal of Lubrication Technology, 1983, 105(4): 638-646. doi: 10.1115/1.3254696 [5] ROGER KU C P, HESHMAT H. Compliant foil bearing structural stiffness analysis: part ITheoretical model, including strip and variable bump foil geometry[J]. Journal of Tribology, 1992, 114(2): 394-400. doi: 10.1115/1.2920898 [6] IORDANOFF I. Maximum load capacity profiles for gas thrust bearings working under high compressibility number conditions[J]. Journal of Tribology, 1998, 120(3): 571-576. doi: 10.1115/1.2834589 [7] IORDANOFF I. Analysis of an aerodynamic compliant foil thrust bearing: method for a rapid design[J]. Journal of Tribology, 1999, 121(4): 816-822. doi: 10.1115/1.2834140 [8] HESHMAT C A, XU D S, HESHMAT H. Analysis of gas lubricated foil thrust bearings using coupled finite element and finite difference methods[J]. Journal of Tribology, 2000, 122(1): 199-204. doi: 10.1115/1.555343 [9] 戚社苗, 耿海鹏, 郭海刚, 等. 弹性箔片动压气体推力轴承承载性能研究[J]. 润滑与密封, 2006(5): 1-4. doi: 10.3969/j.issn.0254-0150.2006.05.001
[10] 闫佳佳. 弹性箔片动压气体推力轴承-转子系统动力学特性研究[D]. 哈尔滨: 哈尔滨工业大学, 2012.
[11] PARK D J, KIM C H, JANG G H, et al. Theoretical considerations of static and dynamic characteristics of air foil thrust bearing with tilt and slip flow[J]. Tribology International, 2007, 41(4): 282-295. [12] LEE D, KIM D. Design and performance prediction of hybrid air foil thrust bearings[J]. Journal of Engineering for Gas Turbines and Power, 2010, 133(4): 042501. [13] 李姗姗, 安琦. 考虑微观形貌时平面型推力滑动轴承润滑性能数值研究[J]. 华东理工大学学报(自然科学版), 2017, 43(5): 733-739.
[14] FENG K, LIU L J, GUO Z Y, ZHAO X Y. Parametric study on static and dynamic characteristics of bump-type gas foil thrust bearing for oil-free turbomachinery[J]. Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, 2016, 230(8): 944-961. doi: 10.1177/1350650115621015 -