2. 中国航发四川燃气涡轮院, 成都 610500
2. China Gas Turbine Establishment, Chengdu 610500, China
作为航空发动机高低压转子间的机械耦合点,中介轴承以其在缩短支点跨距和减轻发动机重量等方面的作用获得了广泛应用。随着航空发动机转子转速日益提高,中介轴承热负荷愈发增大,对于轴承润滑方式及润滑效率提出了更高要求。相较于传统的喷油润滑方式[1-2],通过环下供油孔输油的润滑方式更能高效实现润滑功能,提高轴承润滑效率,因而成为发动机主轴轴承的主要润滑方式。中介轴承润滑系统由复杂的环下供油流道组成,也涉及到喷嘴喷油和流道入口收油关联关系,以及滑油输出时变性影响等问题,由于分析设计较为困难而会影响到中介轴承润滑性能和效率提升。开展环下供油流道中滑油流动特性分析,籍以提升发动机主轴轴承润滑性能,确是一项重要的研究工作。
关于中介轴承环下润滑研究较为鲜见,且多数研究工作更常见于国内以及其他工程领域,一些相关的典型研究工作综述如下。陈军[3]针对中介轴承环下供油方式进行了润滑机理层面的分析,研究发现与传统主轴承润滑方式相比,采用环下润滑方式的轴承具有较强的抗断油能力,且滑油可以直接进入有效润滑区,从而减少了滑油损失。雷燕妮等[4]采用CFD仿真方法研究发现,在中介轴承润滑系统进油口处安装外环挡板有利于增加轴承供油流量。张兆会等[5]探究了中介轴承润滑系统中导流板孔数目对导流板孔、回油槽和回油孔处滑油流量的影响。张永红等[6-7]提出了轴承润滑性能更佳的叶扇/环下油孔组合新结构。喷射-收纳滑油研究多以试验分析为手段,所得结果也更为可靠。Krug等[8]构建了轴承喷油润滑试验平台,借此探究了喷油流量和喷油角度等工况与几何因素对收油效率的影响,通过试验发现随着喷油角度增加,喷射滑油与流道壁面碰撞后的飞溅程度增大,导致收油效率降低。Fuller等[9]通过试验分析了喷油角度对喷射油柱在空气扰动作用下破裂过程的影响。Zhang等[10]开展了喷油速度对轴承内部流场中滑油体积分数的影响研究,所得结果表明滑油体积分数随着喷油速度增加而增大。Wu等[11-12]基于轴承喷油润滑试验研究了滑油体积分数对轴承内部温升的影响,研究结果表明随着喷嘴数目和喷油流量增加,轴承内部流场中滑油体积分数增大,轴承温度随之降低。
通过上述分析可以看出存在以下问题,一是有关滑油流动仿真与喷射-收纳滑油的分析工作是相割裂的,并未进行集成性的研究;二来滑油流动仿真是稳态分析,未考虑环下供油孔与滚动体之间相对位置变化而产生的滑油输出时变性影响。这些问题的存在使得相关仿真分析工作与工程实际有较大差异,掣肘着分析方法和成果在工程中的应用。
本文提出了考虑滑油输出时变性影响的喷油-收油与滑油流动集成分析思路,同步形成了环下供油孔流通状态时序分析方法、喷油-收油效率与轴承润滑效率计算公式,利用嵌入了相关方法和公式的数值仿真模型,开展了若干结构和工况参数条件下的喷油-收油与滑油流动数值仿真。论文提出了较为系统的中介轴承环下供油仿真方法,对于实现中介轴承润滑系统的精确设计有明显的工程意义。
1 中介轴承环下供油流道结构和喷油-收油及滑油流动过程描述图 1示意性地给出了包含喷嘴在内的中介轴承环下供油流道结构。分析中涉及的主要结构参数有:喷嘴喷射角度θ;喷嘴与收油孔距离l;收油环上的收油孔直径d1;环下供油孔直径d2;引流孔直径d3;篦齿密封节流间隙δ1;隙缝密封节流间隙δ2。采用1个喷嘴对轴承进行供油,收油孔个数m1、环下供油孔数目m2和引流孔数目m3在分析中视需要而有变化。密封座上开设引流孔的作用是减少通过篦齿密封的滑油泄漏[13-15]。
|
| 图 1 喷嘴及中介轴承环下供油流道结构 Fig. 1 Structures of jet and intershaft bearing innerring oil-flow tube |
在低压转子转动过程中,当收油孔处于喷嘴喷射范围时,喷射的滑油被收油孔收纳进入环下供油流道;当收油孔不在喷嘴喷射范围时,喷射出来的滑油与高速旋转的收油环壁面碰撞,部分滑油因飞溅而流失,还有一部分沉积形成铺展油膜[16-18],其中沿着圆周方向铺展的滑油因收油环上周向布置有收油孔而最终流进环下供油流道,与直接喷射进入流道中的滑油流量共同构成了收油流量。
经收油孔进入环下供油流道中的滑油,在高低压转子旋转离心作用下向着供油流道各出口流动,流道出口有篦齿密封节流间隙出口,隙缝密封节流间隙出口,环下供油孔出口及引流孔出口,其中环下供油孔出口的滑油是对轴承元件进行润滑的有效供给滑油。由于环下供油孔与轴承滚动体瞬时相对位置变化导致环下供油孔流通状态在开放和封闭之间不断转换,使得环下供油孔出口滑油流量是瞬态变化的,因此也影响到流道其他出口滑油流量的变化。
2 喷油-收油数值分析及环下供油孔流通状态时序分析 2.1 喷油-收油数值分析如前所述,喷嘴喷射出的润滑油一部分直接进入收油孔,还有一部分则通过形成壁面油膜而周向流进收油孔。对于前者,需要通过有限元方法计算出收油孔入口处滑油体积分数及滑油流速,对于后者同样需要在计算出壁面沉积油膜厚度及油膜流动速度,最后,进行相应处理获得进入环下供油流道的总滑油流量(收油流量)。
2.1.1 喷油-收油数值模型喷油-收油有限元分析是借助于ANSYS WORKBENCH 17.0软件平台开展的,分析中涉及的油气两相流体流动区域几何模型如图 2所示,其中图 2(b)是剖面图。图 2(b)所示的流体流动边界条件为:1是喷嘴所在截面(自由截面),边界条件为大气压力,除了喷嘴孔外油气两相流体在此截面自由流通;2为喷嘴孔,其边界条件为喷油流量;收油环和低压转子与流体接触界面(如3所示各界面)的边界条件设为无滑移壁面条件,即流体接触界面速度为收油环(低压转子)速度;4为低压转子油孔出口,边界条件为轴承腔腔压。同时对油膜流动区域的网格进行加密处理,以提高油膜厚度及其周向流速的计算精度。
|
| 图 2 几何模型和油气两相流体流动边界 Fig. 2 Geometry model and flowing boundary ofoil/air two-phase flow |
当低压转子转动至收油孔处于喷嘴喷射范围时,通过有限元计算收油孔自进入到移出喷射范围这个时间段中,各时间点的收油孔进油区域上网格滑油含量及滑油流速,平均处理后得到各时间点收油孔入口滑油相体积分数和滑油流速。
当低压转子转动至收油孔位于喷嘴喷射范围外时,喷射滑油在收油环壁面上形成铺展油膜并沿圆周方向流动,在周向上选取若干点,通过有限元计算出各点的铺展油膜厚度及其周向流动速度,作平均处理后获得铺展油膜平均厚度和平均周向流速。
2.1.2 喷射进入收油环的滑油流量在收油环转动一周的时间T0内,将收油孔i处于喷射范围的时间段划分为m个时间间隔δt,形成m+1个时间点(t1, …, tj, …, tm+1)。通过有限元计算tj时刻收油孔入口滑油相体积分数αitj和滑油流速Vitj,tj时刻喷射进入收油孔i的滑油流量[19]为
| $ q_{i{\rm{d}}}^{{t_j}} = \left( {A\alpha _i^{{t_j}}} \right)V_i^{{t_j}} $ | (1) |
式中:Aαitj为收油孔i在tj时刻的进油区域面积。
在收油环转动一周的时间里,喷射进入收油孔i的滑油量为
| $ {Q_{i{\rm{d}}}} = \sum\limits_{j = 1}^m {\left( {q_{i{\rm{d}}}^{{t_j}}{\rm{ \mathit{ δ} }}t} \right)} $ | (2) |
由此得到喷射进入收油孔i的滑油流量为
| $ {q_{i{\rm{d}}}} = {Q_{i{\rm{d}}}}/{T_0} $ | (3) |
若收油环上有n个收油孔,则喷射进入收油环的滑油流量为
| $ {q_{\rm{d}}} = \sum\limits_{i = 1}^n {{q_{i{\rm{d}}}}} $ | (4) |
通过有限元计算出沉积于收油孔之间壁面的铺展油膜平均厚度h及油膜平均周向流速v,沿周向流入收油孔i的滑油流量为
| $ {q_{i{\rm{c}}}} = \frac{{h{d_1}v\Delta t}}{{{T_0}}} $ | (5) |
式中:Δt为低压转子转过两收油孔夹角的时间;铺展油膜流动的轴向宽度范围取为收油孔直径d1。
沿周向流入收油环的滑油流量为
| $ {q_{\rm{c}}} = \sum\limits_{i = 1}^n {{q_{i{\rm{c}}}}} $ | (6) |
收油流量q为qd与qc之和,收油效率ηcatch为q与喷油流量qinject之比,即
| $ q = {q_{\rm{d}}} + {q_{\rm{c}}} $ | (7) |
| $ {\eta _{{\rm{catch}}}} = \frac{q}{{{q_{{\rm{inject}}}}}} $ | (8) |
由于高低压转子的反向转动,环下供油孔与轴承滚动体的瞬时相对位置发生变化而导致供油孔流通状态在开放和封闭之间不断转换,使得供油孔出口的边界条件瞬态变化。时序分析的目的就在于通过揭示环下供油孔流通状态时变规律,为供油孔出口提供瞬态边界条件。结合高低压转子转速,可以确定流通状态变化周期,即为同一环下供油孔出现两次相同流通状态所经历的时间。
2.2.1 环下供油孔流通状态变化周期如图 3所示,将环下供油孔1被滚动体1完全遮挡作为轴承运行初始时刻(t=0 s),根据流通状态变化周期的定义,当轴承运行到周期T时,环下供油孔1被相邻滚动体(2或N)再次遮挡。此时环下供油孔1和滚动体1从各自初始位置转过角度的差值为相邻滚动体夹角φ,即
| $ \varphi = \left| {\theta _1^T - \psi _1^T} \right| = \frac{{{{360}^\circ }}}{N} $ | (9) |
|
| 图 3 t=0 s时环下供油孔与滚动体的相对位置 Fig. 3 Relative position of inner ring holes and rollers(t=0 s) |
式中:N为滚动体数目;θ1T和ψ1T分别是环下供油孔1和滚动体1从各自初始位置转过的角度,分别表示为
| $ \theta _1^T = {n_1}T \times {6^\circ } $ | (10) |
| $ \psi _1^T = \left[ {{n_1}\left( {1 - \gamma } \right)/2 + {n_2}\left( {1 + \gamma } \right)/2} \right]T \times {6^\circ } $ | (11) |
式中:n1和n2分别为高压转子和低压转子的转速,n1(1-γ)/2+n2(1+γ)/2为滚动体公转转速[20-22],在图 3中以顺时针转向为正向,单位为r/min;T为流通状态变化周期,单位为s;γ为滚动体直径和轴承节圆直径之比。
将式(10)和式(11)代入式(9)中,求得流通状态变化周期为
| $ T = \frac{1}{{\left| {{n_1} - {n_2}} \right|}} \cdot \frac{1}{{1 + \gamma }} \cdot \frac{{120}}{N} $ | (12) |
环下供油孔的流通状态有:供油孔被滚动体完全遮挡的封闭态、未被遮挡的开放态和介于其间的半封闭态。根据轴承运行t时刻的环下供油孔与其相邻滚动体的周向坐标,确定二者相对位置关系,由此得到环下供油孔的瞬时流通状态。
t时刻环下供油孔k的周向坐标为
| $ \theta _k^t = \left( {k - 1} \right)\frac{{{{360}^\circ }}}{{{m_2}}} + {n_1}t \times {6^\circ } $ | (13) |
式中:m2为环下供油孔数目。
供油孔k两侧滚动体的编号为k1、k2,它们在t时刻的周向坐标分别为
| $ \varphi _{{k_1}}^t = \left\lfloor {\theta _k^t/\varphi } \right\rfloor \varphi $ | (14) |
| $ \varphi _{{k_2}}^t = \left\lceil {\theta _k^t/\varphi } \right\rceil \varphi $ | (15) |
式中:
分别计算φtk1和φtk2与θkt之间的差值:
| $ {A_1} = \theta _k^t - \varphi _{{k_1}}^{{t_1}} $ | (16) |
| $ {A_2} = \theta _k^t - \varphi _{{k_2}}^t $ | (17) |
若A1 < A2,说明供油孔k更靠近滚动体k1,此时通过A1的大小判断供油孔k的流通状态,反之则以A2判断流通状态,具体判断准则如下:A1或A2 < 1.25°,供油孔为封闭态;1.25° < A1或A2 < 2.5°,供油孔为半封闭态;A1或A2>2.5°,供油孔为开放态。
将流通状态变化周期T分为若干时间节点,对每一个时间节点t均做上述判断,以决定此时间节点的环下供油孔流通状态,进而在滑油流动瞬态数值分析中施加时变边界条件。
3 环下供油流道中滑油流动数值分析及轴承润滑效率计算 3.1 滑油流动数值分析环下供油流道中滑油流动数值分析涉及的计算区域如图 4(a)所示,滑油流动边界如图 4(b)所示。
|
| 图 4 几何模型和滑油流动边界 Fig. 4 Geometry model and oil flowing boundary |
将流通状态变化周期T均分为k个时间间隔δT,形成k+1个时间节点(t1, …, tj, …, tk+1),各时间节点上的边界条件有时不变和时变两种类型,前者结合图 4(b)如表 1所示;后者则根据时序分析方法判断此时间节点上供油孔的流通状态,并依据封闭态、开放态、半封闭态,分别赋予出口流量为0、出口压力为轴承腔压力、开放区域部分与封闭区域部分各取对应边界条件。
| 边界编号 | 边界条件 | 参数符号/单位 |
| 1 | 收油流量 | q/(L·min-1) |
| 2、3 | 轴承腔腔压 | p1/MPa |
| 4 | 大气压 | p2/MPa |
| 5 | 低压转子转速 | n2/(r·min-1) |
| 6 | 高压转子转速 | n1/(r·min-1) |
环下供油流道出口有篦齿密封出口,隙缝密封出口,环下供油孔出口及引流孔出口,其中只有环下供油孔出口的滑油流量有助于轴承的润滑。在一个流通状态变化周期中进行滑油流动数值分析,可获得上述4个出口的滑油流量。
3.2 轴承润滑效率计算计算出的环下供油孔出口瞬时滑油流量
| $ \overline {{q_{\rm{B}}}} = \frac{1}{T}\sum\limits_{j = 1}^k {\left( {q_{\rm{B}}^{{t_j}}\delta T} \right)} $ | (18) |
则轴承润滑效率为
| $ {\eta _{\rm{B}}} = \frac{{\overline {{q_{\rm{B}}}} }}{{{q_{{\rm{inject}}}}}} $ | (19) |
结构参数影响分析中采用的工况条件为:喷油流量qinject=6 L/min;高低压转子反向转动,高压转子转速n1=11 000 r/min、低压转子转速n2=-8 000 r/min;轴承腔压力p1=0.27 MPa。结构参数及其取值如表 2所示。
| 结构参数 | 参数符号/单位 | 数值 |
| 喷嘴喷射角度 | θ/(°) | 30, 45, 60, 75 |
| 环下供油孔数目 | m2 | 4, 6, 8, 10 |
| 滚动体数目 | N | 26, 28, 30, 32 |
进行喷嘴喷射角度影响分析时,环下供油孔数和滚动体数分别为6和30。分析中环下供油孔流通状态变化周期为T=0.189 ms,环下供油孔在11个时间节点上的流通状态如图 5所示。
|
| 图 5 环下供油孔在11个时间节点上的流通状态 Fig. 5 Open-shut modes of inner ring holes on 11moments |
图 6给出了直接喷射进入收油环滑油流量qd、沿周向流入收油环滑油流量qc、收油流量q和收油效率ηcatch随着喷嘴喷射角度θ的变化规律。从图中可以看出,随着θ增加,qd呈现出先增大后降低的趋势,在θ=45°时qd达到最大。原因在于当低压转子转动至收油孔处于喷嘴喷射范围时,θ=45°下滑油喷射路径未与收油环壁面发生干涉,喷射的滑油全部进入收油孔,但当θ=30°、θ=60°和θ=75°时,滑油喷射路径与收油环壁面发生干涉,即喷射的滑油有一部分被收油环壁面反弹,只有部分滑油喷射进入收油孔。
|
| 图 6 喷嘴喷射角度对收油流量和收油效率的影响 Fig. 6 Effect of impingement angles on collectedoil flow rate and catch efficiency |
当收油孔处于喷嘴喷射范围外时,随着θ增大,喷射滑油与收油环壁面碰撞后的飞溅程度减小,使得润滑油更易于沉积于收油环壁面上而导致铺展油膜厚度增大,故沿周向流入收油孔的滑油流量qc也随之增大。作为qd与qc的叠加结果,q随θ的增加也呈现出先增大后降低的趋势,同样在θ=45°时达到最大。在喷油流量一定的条件下,ηcatch与q的变化趋势显然是一致的。
图 7给出了环下供油流道各出口瞬时滑油流量(环下供油孔出口瞬时流量qB、引流孔出口瞬时流量qS、篦齿密封和隙缝密封出口瞬时流量qL和qG)随喷嘴喷射角度θ的变化规律。
|
| 图 7 喷嘴喷射角度对供油流道各出口瞬时滑油流量的影响 Fig. 7 Effect of impingement angles on transientoil flow rate of each outlet |
从图 7中可以看出,由于收油流量(环下供油流道入口流量)随θ增加而呈现先增大后降低的趋势,因此环下供油流道各出口的瞬时滑油流量亦随θ如此变化。由于引流孔的分流效应,篦齿密封出口瞬时滑油流量qL低于隙缝密封出口瞬时滑油流量qG,不过由于二者的直通性,其流量都远大于引流孔出口瞬时流量qS,后者流量较小除与结构有关外,也与高低压转子转速差有关系。在给定的流通状态变化周期中,环下供油孔流通状态经历了封闭-开放-封闭的过程,受其直接影响的环下供油孔出口瞬时滑油流量呈现出先增大后降低的趋势,根据质量守恒原理,其他3个出口的滑油流量时变规律则与之相反。
图 8给出了喷嘴喷射角度θ对供油流道各出口平均滑油流量(轴承供油流量
|
| 图 8 喷嘴喷射角度对供油流道各出口平均滑油流量和轴承润滑效率的影响 Fig. 8 Effect of impingement angles on average oil flow rate of each outlet and bearing lubrication efficiency |
由于
进行环下供油孔数目影响分析时,喷嘴喷射角度为45°,滚动体数目为30,环下供油孔流通状态变化周期为T=0.189 ms。分析中当m2为4时,同一时刻下环下供油孔1和3、2和4的流通状态相同;当m2为8时,同一时刻下环下供油孔1和5、2和6、3和7、4和8的流通状态相同,故图 9给出了m2=4时环下供油孔1和2,m2=8时环下供油孔1、2、3和4的时变流通状态。m2为6和10时,同一时刻下各环下供油孔的流通状态相同,具体如图 5所示。
|
| 图 9 m2=4和m2=8时各环下供油孔的时变流通状态 Fig. 9 Variation of open-shut mode of each inner ring hole over time (m2=4, 8) |
分析中发现环下供油孔数目变化对收油流量和收油效率影响不大,其原因可以归结为:虽然环下供油孔数目变化影响了其出口流量,但由于其他3个出口排放流量的调剂,“稀释”了供油孔数目变化对收油流量和收油效率的影响。实际上在后面分析中,滚动体数目变化对收油流量和收油效率影响也不大,原因是相同的。
图 10给出了环下供油流道各出口瞬时滑油流量随环下供油孔数目m2的变化规律。从图中可以看出当m2=6和m2=10时,由于环下供油孔流通状态呈现出封闭-开放-封闭的规律,故而其出口瞬时滑油流量qB也表现出先增大再减小的趋势,其他3个出口瞬时滑油流量(qL、qS、qG)的时变规律与其相反;而当m2=4和m2=8时,由于各环下供油孔封闭-开放的转换规律都不同,故而环下供油孔出口瞬时滑油流量也发生对应的波动,同样影响到其他3个出口瞬时滑油流量也出现波动变化。
|
| 图 10 环下供油孔数目对各出口瞬时滑油流量的影响 Fig. 10 Effect of inner ring holes' number on transient oil flow rate of each outlet |
图 11给出了环下供油孔数目对供油流道各出口平均滑油流量和轴承润滑效率的影响。容易解释的是,随着m2的增加,环下供油孔出口平均滑油流量
|
| 图 11 环下供油孔数目对各出口平均滑油流量及轴承润滑效率的影响 Fig. 11 Effect of inner ring holes' number on average oil flow rate of each outlet and bearing lubrication efficiency |
在进行滚动体数目的影响分析时,喷嘴喷射角度为45°,环下供油孔数目为6。分析中当滚动体数目N分别为26、28、30和32时,环下供油孔流通状态变化周期为0.219 ms、0.203 ms、0.189 ms和0.178 ms,即周期随滚动体数目增加而变短;当N为26、28和32时,滚动体数与环下供油孔数不成倍数关系,因此同一时刻下各环下供油孔的流通状态不一致,而N=30时,滚动体数为环下供油孔数的整数倍,故同一时刻下各环下供油孔的流通状态相同。
图 12给出了环下供油流道各出口瞬时滑油流量随滚动体数目N的变化规律。图中可以看出当N=26、28和32时,同一时刻下各环下供油孔的流通状态不同,导致环下供油孔出口及其他3个出口的瞬时滑油流量发生波动变化;而N=30时因各环下供油孔开放——封闭状态的转换规律相同,故供油流道各出口的瞬时滑油流量变化趋势较为规律。
|
| 图 12 滚动体数目对各出口瞬时滑油流量的影响 Fig. 12 Effect of rollers' number on transient oil flow rate of each outlet |
图 13给出了滚动体数目对环下供油流道各出口平均滑油流量和轴承润滑效率的影响。从图中可以发现,由于滚动体数目增加使得环下供油孔流通状态变化周期变短,亦即环下供油孔处于开放状态的时长变短,其出口的平均滑油流量和轴承润滑效率也随之减小,自然其他3个出口的平均滑油流量随之增大。
|
| 图 13 滚动体数目对各出口平均滑油流量和轴承润滑效率的影响 Fig. 13 Effect of rollers' number on average oil flow of each outlet |
工况参数影响分析时,喷嘴喷射角度为45°,环下供油孔数目为6,滚动体数目为30;工况参数及其取值如表 3所示。
| 工况参数 | 参数符号/单位 | 数值 |
| 喷油流量 | qinject/(L·min-1) | 6 |
| 低压转子转速 | n2/(r·min-1) | 5 000, 6 500, 8 000, 9 500 |
| 高压转子转速 | n1/(r·min-1) | 11 000 |
| 轴承腔压力 | p1/MPa | 0.27 |
当高低压转子转向相反时,高压转子转速n1为11 000 r/min,低压转子转速n2分别为-5 000、-6 500、-8 000、-9 500 r/min,环下供油孔流通状态变化周期对应为0.225、0.206、0.189、0.176 ms。需要说明的是:低压转子转速中的负号表示低压转子转向为图 3中的逆时针方向,与高压转子转向相反。
当n2=-8 000 r/min时各环下供油孔在11个时间节点上的流通状态见图 5,n2取其他数值时环下供油孔在各时间点的流通状态与之相同。
图 14给出了直接喷射进入收油环滑油流量qd、沿周向流入收油环滑油流量qc、收油流量q和收油效率ηcatch随低压转子转速n2的变化规律。
|
| 图 14 低压转子转速对收油流量和收油效率的影响 Fig. 14 Effect of low-pressure rotor speeds on collectedoil flow rate and catch efficiency |
从图 14中可以看出,随着n2增加,qd减小而qc增大,原因在于当低压转子转速较高时,收油孔处于喷嘴喷射范围的时间较短,这导致直接喷射进入收油环的滑油流量qd较少;此外随着n2增加,收油环壁面上铺展油膜的周向流速增大,故而周向流入收油环的滑油流量qc也随之增大。作为qd和qc的叠加结果,q随n2增大而减小,在喷油流量一定的前提下,ηcatch与q的变化趋势显然是一致的。
图 15给出了环下供油流道各出口瞬时滑油流量随低压转子转速n2的变化规律。从图中可以看出,由于收油流量q随n2增加而减小,故而环下供油流道各出口的瞬时滑油流量亦随n2如此变化。
|
| 图 15 高低压转子转向相反时低压转子转速对各出口瞬时滑油流量的影响 Fig. 15 Effect of low-pressure rotor speeds on transientoil flow rate of each outlet |
从图 15中还可以发现,随着n2增加,环下供油孔出口瞬时滑油流量的峰值减小。其原因是:由于高低压转子转向相反,因此二者转速之差随低压转子转速增加而增大,环下供油孔流通状态变化周期随之变短,亦即环下供油孔处于开放状态的时长减少,其出口瞬时滑油流量的峰值也随之减小。
图 16给出了低压转子转速对环下供油流道各出口平均滑油流量和轴承润滑效率的影响。从图 16中可以发现,由于环下供油流道各出口的瞬时滑油流量随n2增加而减小,因此各出口的平均滑油流量和轴承润滑效率也随之减小。
|
| 图 16 低压转子转速对轴承润滑效率的影响 Fig. 16 Effect of low-pressure rotor speeds on bearinglubrication efficiency |
当高低压转子转向相同时,低压转子转速n2分别为5 000、6 500、8 000、9 500 r/min,环下供油孔流通状态变化周期对应为0.601 1、0.801 5、1.202、2.404 ms。各环下供油孔在11个时间节点的流通状态与n2=-8 000 r/min条件下的流通状态相同。
由于收油环与低压转子连接而共同转动,故而高低压转子转向的异同不影响收油流量和收油效率,当高低压转子转向相同时低压转子转速对收油流量和收油效率的影响亦如图 14所示。
图 17给出了低压转子转速对环下供油流道各出口瞬时滑油流量的影响,从图中可以发现,随着低压转子转速增加,环下供油孔出口瞬时滑油流量的峰值增大。其原因是:由于高低压转子转向相同,故而高低压转子转速之差随低压转子转速增加而减小,环下供油孔流通状态变化周期随之变长,即环下供油孔处于开放状态的时长增加,其出口滑油流量的最大值也随之增大。
|
| 图 17 高低压转子转向相同时低压转子转速对各出口瞬时滑油流量的影响 Fig. 17 Effect of low-pressure rotor speeds on transient oil flow rate of each outlet |
图 18给出了低压转子转速对各出口平均滑油流量和轴承润滑效率的影响,从图中可以看出,由于环下供油孔的开放时长随低压转子转速增大而增加,故而其出口的平均滑油流量和轴承润滑效率亦随之增大。如前所述收油流量(环下供油流道各出口平均滑油流量之和)随低压转子转速增加而减小,因此其他出口(除环下供油孔)的平均滑油流量也随之减小。
|
| 图 18 低压转子转速对各出口平均滑油流量和轴承润滑效率的影响 Fig. 18 Effect of low-pressure rotor speeds on average oil flow rate of each outlet and bearing lubrication efficiency |
从图 18中还可以看出,相较于高低压转子转向相反情况下的轴承润滑效率ηB,二者转向相同时的ηB有较大程度的提升。原因在于高低压转子转向相同时,环下供油孔处于开放态的时间更长,其出口的平均滑油流量也更大,故而轴承润滑效率更高。
4.3 本文分析方法与目前常用分析方法的对比为了体现本文提出的滑油流动瞬态分析方法的先进性,本节在3组结构和工况条件下(见表 4),进行了本文的瞬态分析方法和目前常用的稳态分析方法的对比分析。由于滑油流动稳态分析未考虑因滚动体与环下供油孔相对位置变化而产生的滑油输出时变影响,故在稳态分析中,将轴承滚动体与环下供油孔的相对位置固定(这与实际情况相悖),即环下供油孔始终处于不变的流通状态,因此这里的稳态分析是恒定开放态和半封闭态下的滑油流动稳态分析。
| 参数 | 单位 | 数值 | ||
| 组1 | 组2 | 组3 | ||
| 喷嘴喷射角度 | (°) | 45 | 60 | 45 |
| 环下供油孔数目 | 4 | 6 | 6 | |
| 滚动体数目 | 30 | 26 | 30 | |
| 喷油流量 | L·min-1 | 6 | 6 | 6 |
| 低压转子转速 | r·min-1 | -8 000 | 5 000 | -5 000 |
| 高压转子转速 | r·min-1 | 11 000 | 12 000 | 11 000 |
| 轴承腔压力 | MPa | 0.27 | 0.27 | 0.27 |
图 19(a)、19(b)、19(c)分别给出了参数组1、参数组2和参数组3条件下通过瞬态和稳态分析获得的环下供油孔出口瞬时滑油流量qB。在参数组1和参数组2条件下由于瞬态分析考虑了环下供油孔流通状态时变性影响,且滚动体数与环下供油孔数不成倍数关系,故同一时刻下各环下供油孔的流通状态不一致,导致qB随时间波动变化;对于参数组3,由于滚动体数为供油孔数的整数倍,因此同一时刻下各供油孔的流通状态相同,qB随时间的变化呈现出先增大后减小的趋势;此外,参数组3条件下各供油孔出现了同时开放的情况,而参数组1和2条件下未出现各供油孔同时开放的状态,故参数组3条件下qB的峰值更高些。
|
| 图 19 3组参数下瞬态和稳态分析对应的环下供油孔出口瞬时滑油流量 Fig. 19 Transient oil flow rate of inner ring holes outlets obtained by transient and steady analysis under three parameter groups |
对于参数组1、2、3,通过稳态分析方法获得的恒定流通状态下供油孔出口滑油流量不随时间变化,未能反映出供油孔出口瞬时滑油流量的时变特点,这显然是不符合工程实际的。
从图 19给出的对比结果还可以看出,本文提出的环下供油孔瞬态分析方法能够体现因滚动体与环下供油孔相对位置变化而形成的滑油输出时变性影响,因其更符合工程实际而具有先进性。
5 结论1) 由喷嘴喷射进入中介轴承环下供油流道的滑油由两部分组成,即直接喷射进入流道的滑油和收油环壁面上沿周向流入流道的滑油;环下供油孔流通状态变化周期与高低压转子转速、轴承滚动体数目相关。
2) 随着喷嘴喷射角度增大,进入中介轴承环下供油流道的滑油流量和收油效率先增加后减小,在喷射角为45°时达到最大;环下供油流道各出口的滑油流量和轴承润滑效率亦随喷射角如此变化。
3) 环下供油孔数目对收油效率的影响不大;随着环下供油孔数目增加,环下供油孔出口的滑油流量和轴承润滑效率均增大而其他出口的滑油流量减小。
4) 滚动体数目对收油效率的影响不大;随着滚动体数目增加,环下供油孔流通状态变化周期变短,环下供油孔出口的滑油流量和轴承润滑效率均减小而其他出口的滑油流量增大。
5) 收油效率随低压转子转速增大而减小,且高低压转子转向的异同不影响收油效率。当高低压转子转向相反时,环下供油流道各出口的滑油流量和轴承润滑效率均随低压转子转速增大而降低;当转向相同时,环下供油孔出口的滑油流量和轴承润滑效率随低压转子转速增加而增大,而其他出口的流量随之减少,且相较于转向相反的情况,转向相同条件下的轴承润滑效率有较大程度的提升。
| [1] | PINEL S I, SIGNER H R, ZARETSKY E V. Comparison between oil-mist and oil-jet lubrication of high-speed, small-bore, angular-contact ball bearings[J]. Tribology Transactions, 2001, 44(3): 327-338. |
| Click to display the text | |
| [2] | JIANG S, MAO H. Investigation of the high speed rolling bearing temperature rise with oil-air lubrication[J]. Journal of Tribology, 2011, 133(2): 021101. |
| Click to display the text | |
| [3] |
陈军.涡扇发动机主轴承环下供油方法分析研究[D].西安: 西北工业大学, 2007. CHEN J. Research on the inner ring oil supplying method of mainshaft bearing in turbofan engine[D]. Xi'an: Northwestern Polytechnical University, 2007(in Chinese). |
| Cited By in Cnki (4) | Click to display the text | |
| [4] |
雷燕妮, 苏华, 陈国定, 等. 进油口结构型式对轴承润滑性能影响分析[J]. 润滑与密封, 2005(1): 82-83. LEI Y N, SU H, CHEN G D, et al. Analysis of the influences of lubricating oil entrance's configuration on bearing's lubrication performance[J]. Lubrication and Seal, 2005(1): 82-83. (in Chinese) |
| Cited By in Cnki (3) | Click to display the text | |
| [5] |
张兆会, 张永红, 陈国定, 等. 航空发动机中介轴承润滑系统的流场分析[J]. 轴承, 2004(12): 22-24. ZHANG Z H, ZHANG Y H, CHEN G D, et al. Analysis of oilflow field in lubrication system serving for an aero engine intershaft bearing[J]. Bearing, 2004(12): 22-24. (in Chinese) |
| Cited By in Cnki (1) | Click to display the text | |
| [6] |
张永红, 宋乐民, 陈国定, 等. 基于中介轴承润滑性能改善的新结构研究[J]. 润滑与密封, 2005(3): 124-125. ZHANG Y H, SONG L M, CHEN G D, et al. The research of new structure aimed to improve the lubrication performance[J]. Lubrication and Seal, 2005(3): 124-125. (in Chinese) |
| Cited By in Cnki (1) | Click to display the text | |
| [7] |
张永红, 陈国定, 宋乐民, 等. 叶扇/环下孔组合结构参数对轴承润滑性能的影响分析[J]. 机械传动, 2005, 29(2): 55-57. ZHANG Y H, CHEN G D, SONG L M, et al. Analysis on the bearing lubrication performance influenced by fan-like/inner-ring hole structure parameter[J]. Journal of Mechanical Transmission, 2005, 29(2): 55-57. (in Chinese) |
| Cited By in Cnki | Click to display the text | |
| [8] | KRUG M B, PEDUTO D, KURZ W, et al. Experimental investigation into the efficiency of an aero engine oil jet supply system[J]. Journal of Engineering for Gas Turbines and Power, 2015, 137: 011505. |
| Click to display the text | |
| [9] | FULLER R, WU P K, KIRKENDALL K, et al. Effects of injection angle on the breakup processes of liquid jets in subsonic crossflows[C]//33rd Joint Propulsion Conference and Exhibit, 1997: 2966. |
| [10] | ZHANG R, WEI C, WU W, et al. CFD investigation on the influence of jet velocity of oil-jet lubricated ball bearing on the characteristics of lubrication flow field[C]//2015 International Conference on Fluid Power and Mechatronics (FPM), 2015: 1324-1328. |
| [11] | WU W, HU J, YUAN S, et al. Numerical and experimental investigation of the stratified air-oil flow inside ball bearings[J]. International Journal of Heat & Mass Transfer, 2016, 103: 619-626. |
| Click to display the text | |
| [12] | WU W, HU C, HU J, et al. Jet cooling for rolling bearings:Flow visualization and temperature distribution[J]. Applied Thermal Engineering, 2016, 105: 217-224. |
| Click to display the text | |
| [13] |
杨淋麟, 张丽, 朱惠人. 篦齿结构对封严性能的影响[J]. 科学技术与工程, 2013, 13(33): 10069-10073. YANG L L, ZHANG L, ZHU H R. Effect of labyrinth structure on seal performance[J]. Science Technology and Engineering, 2013, 13(33): 10069-10073. (in Chinese) |
| Cited By in Cnki | Click to display the text | |
| [14] |
张留祥, 陈俊东, 耿旭辉, 等. 直通式篦齿密封性能的数值模拟与试验研究[J]. 润滑与密封, 2013(10): 51-54. ZHANG L X, CHEN J D, GENG X H, et al. Numerical simulation and experimental investigation on straight-through labyrinth seal performance[J]. Lubrication and Seal, 2013(10): 51-54. (in Chinese) |
| Cited By in Cnki (1) | Click to display the text | |
| [15] |
杨军.直通型封严篦齿流动特性的试验研究及数值计算[D].成都: 电子科技大学, 2009. YANG J. Numerical calculation and experimental investigation on straight-through labyrinth seal flow characteristic[D]. Chengdu: University of Electronic Science and Technology of China, 2009(in Chinese). |
| Cited By in Cnki (7) | Click to display the text | |
| [16] |
王莉娜, 陈国定, 孙恒超. 轴承腔油滴沉积特性及油膜流动特征分析[J]. 航空学报, 2016, 37(10): 3159-3169. WANG L N, CHEN G D, SUN H C. Characteristics analysis of oil droplet deposition and oil film flow in a bearing chamber[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(10): 3159-3169. (in Chinese) |
| Cited By in Cnki (1) | Click to display the text | |
| [17] |
陈薄, 陈国定, 王涛. 考虑油滴变形和二次油滴效应的轴承腔壁面油膜流动分析[J]. 航空学报, 2013, 34(8): 1980-1989. CHEN B, CHEN G D, WANG T. Flow characteristics analysis of wall oil film with consideration of oil droplet deformation and secondary oil droplet[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(8): 1980-1989. (in Chinese) |
| Cited By in Cnki (7) | Click to display the text | |
| [18] | CHEN B, CHEN G D, SUN H C, et al. Effect of oil droplet deformation on its deposited characteristics in an aeroengine bearing chamber[J]. Proceedings of the Institution of Mechanical Engineers, Part G:Journal of Aerospace Engineering, 2014, 228(2): 206-218. |
| Click to display the text | |
| [19] | MUNSON B R, OKIISHI T H, ROTHMAYER A P, et al. Fundamentals of fluid mechanics[M]. New York: John Wiley & Sons, 2014. |
| [20] |
徐建东, 陈国定. 高速圆柱滚子轴承运动的简化分析[J]. 机械科学与技术, 1993(4): 57-60. XU J D, CHEN G D. Simplified motion analysis of high speed roller bearings[J]. Mechanical Science and Technology for Aerospace Engineering, 1993(4): 57-60. (in Chinese) |
| Cited By in Cnki | Click to display the text | |
| [21] |
张成铁, 陈国定, 李建华. 高速滚动轴承的动力学分析[J]. 机械科学与技术, 1997, 16(1): 136-139. ZHANG C T, CHEN G D, LI J H. Dynamic analysis of high speed roller bearings[J]. Mechanical Science and Technology for Aerospace Engineering, 1997, 16(1): 136-139. (in Chinese) |
| Cited By in Cnki (42) | Click to display the text | |
| [22] |
刘秀海, 邓四二, 滕弘飞. 高速圆柱滚子轴承保持架运动分析[J]. 航空发动机, 2013, 39(2): 31-38. LIU X H, DENG S E, TENG H F. Kinematics analysis of cages in high-Speed cylindrical roller bearings[J]. Aeroengine, 2013, 39(2): 31-38. (in Chinese) |
| Cited By in Cnki (3) | Click to display the text |



