-
摘要: 采用球-盘点接触光干涉润滑油膜测量装置,观察了纯滚条件下近接触区油池形态随供油量和卷吸速度的变化过程. 结果表明,油池在低速、中速和高速下分别呈现出闭合态、半开放态和分离态3种状态. 闭合态油池内气穴长度随卷吸速度增加而增加,且存在使尾部油池分离的临界速度. 气穴区表面张力作用使气穴长度预测值和测量值出现差别;卷吸作用和压差作用是油池变形的原因;供油量和卷吸速度决定了润滑剂的供-失平衡态,影响了油池侧宽的变化趋势.Abstract: Using the ball-on-disc optical lubricating film measurement device, the evolution of oil reservoir patterns near pure rolling contact was observed. It was found that three patterns of oil reservoirs, i.e., closed, half-opened, and fully separated were formed under low speed, moderate speed, and elevated speed, respectively. For the closed pattern, the dimensionless cavitation length almost proportionally increased with the increase in the entrainment speed. There existed a critical entrainment speed beyond which the rear reservoir was split by a cavitation bubble. It was the parameter of surface tension at the boundary of the cavitation bubble, that resulted in the deviation of predicted cavitation length from that of measured one. The entrainment effect and the pressure difference cause the deformation of the oil reservoir shape. The lubricant feed-loss balance was determined by the oil charging amount and entrainment speed. Consequently, the trends of the width of the side oil reservoir with entrainment speed were different.
-
油膜润滑将接触副两表面分离,把固体壁面摩擦转化为流体的内摩擦,是减少机械装备零部件(如滚动轴承)摩擦磨损的有效方法. 基于零部件运行可靠性考虑,工程中常采用大量外部供油试图保证接触副间油膜的充分建立,现有的滚动轴承使用寿命估算方法也大都建立在充分供油条件下[1]. 而在滚动轴承内部,滚动体与内外圈接触所产生的压力梯度使润滑剂向滚道两侧迁移(侧泄),在高速[2]、低温[3]和脂润滑[4]等条件下,迁移出的润滑剂难以有效向滚道回填,进而使接触副入口处润滑剂供给不足,产生运动乏油并导致油膜厚度降低. 在限制性供油条件下,例如使用油气润滑方式,定量的润滑剂穿透气帘被输送到轴承内部[5],但由于润滑剂量供给较少,接触区大都处于乏油状态,且润滑状态与接触区周围的油池几何参数、形态及动态演化呈现出较强的相关性,而该方面的研究并不充分,因而有必要对此进行细致的试验观察与分析.
针对乏油润滑机理研究,Wedeven等[6]率先用光干涉技术对乏油润滑现象进行了可视化观察,并以入口距离量化了入口区供油状态与膜厚降低间的关系. 此后四十年,乏油润滑试验与理论研究得到了广泛的开展[7],并用于指导滚动轴承的润滑设计[8]. 在乏油润滑基础试验中,为了便于获取或复现乏油状态,常采用微量定量的供油方式,例如Qian等[9]观察了供油量为20~100 μl条件下的油池形态与润滑状态. 鉴于油气润滑微液滴定量的供油特征,王少朋等[10]初步观察了单个微油滴穿过接触区的动态润滑过程. 刘成龙等[11]探究了油滴大小、卷吸速度和载荷等因素影响下微油滴的润滑行为. 李书义等[12]建立了微油滴供油润滑模型,模拟了微油滴体积、卷吸速度和润滑油黏度等对润滑状态的影响.
多数的乏油润滑研究关注接触区内部润滑状态,而外部润滑介质分布对于入口区供油及油膜厚度存在明显影响. 例如,Chiu等[13]分析了表面张力作用下滚道两侧润滑剂的回填; Li等[14]和李哲等[15]分别用化学和物理的方法制备润湿性梯度表面,有效增强了润滑剂的回填作用;韩兵等[16]研究了毛细力和楔裂压作用对润滑剂自发回填的影响. 此外,江楠等[17]考虑滚动体的滑动,观察了滑滚条件下定量润滑剂的回填效应,在实质上,增强润滑剂回填是为了改变接触区周围的油池形态,对油池形态进行细致研究是实现润滑剂回填有效调控的前提. 已有的研究已经表明油池几何特征受到多种因素影响[18],Emden等[19]发现了“闭合”与“分离”油池形态. 目前,针对油池形态的研究仍然不够充分,油池形态的动态演化过程及其对润滑剂在壁面分布的影响等问题尚需进一步分析. 本文中采用球-盘点接触光干涉润滑油膜测量装置,对供油量为0.1~0.4 μl条件下的油池形态进行观察,通过细致观察油池尾部气穴区变化规律与油池分离临界条件,定量分析油池几何参数的变化,探究油池尾部布油和入口供油对油池形态演化的影响机理.
1. 试验部分
试验在球-盘点接触光干涉润滑油膜测量装置上进行,如图1所示,钢球与蓝宝石盘组成接触副,以蓝宝石盘带动钢球的方式实现纯滚动. 试验所用蓝宝石盘直径为150 mm;钢球为G5精度,直径为25.4 mm,表面粗糙度约为14 nm. 加载装置位于球托下方,可实现钢球对蓝宝石盘的加载. 接触区附近的油池分布经显微镜放大后被CCD捕获并储存,为了采集油池动态变化过程,图像的采集模式设置为连续采集,采集帧率依据试验要求设定.
表1中列出了试验条件,其中ue=(ud + ub)/2,ue为卷吸速度,ud和ub分别为蓝宝石盘和钢球接触点线速度. 为了考察恒定载荷下卷吸速度对油池的影响,试验设定的卷吸速度变化范围为ue=0.2~400 mm/s. 试验采用的润滑油为PAO10,其性质列于表2中. 在采集试验数据之前,先将球盘进行表面清洁,再将定量的润滑剂布置于蓝宝石盘表面,然后将卷吸速度设置为ue=0.2 mm/s使油滴经过接触区,随后将卷吸速度设置为ue=0.5 mm/s,并在该速度下使蓝宝石盘转过4~5圈,待润滑剂均匀地铺设在蓝宝石盘接触轨道上后进行试验数据的采集.
表 1 试验条件Table 1. Experimental conditionsParameters Specifications Entrainment speed, ue/(mm/s) 0.2~400 Load/N 10 Maximum Hertz pressure/GPa 0.234 Lubricant amount/μl 0.1; 0.2; 0.3; 0.4 Ambient temperature/℃ 20±1 表 2 试验用润滑油性质Table 2. Properties of lubricants used in experimentsLubricant Viscosity at 20 ℃
/(Pa·s)Density at 20 ℃
/(kg/m3)PAO10 0.13 835 2. 试验结果
2.1 供油量对气穴长度的影响
在图2(a)中给出了纯滚条件下,供油量分别为0.1、0.2、0.3和0.4 μl时油池外形随速度变化干涉图. 为了量化油池形状,图中定义a为赫兹接触半径,l1为气穴长度,l2为侧油池宽度. 当供油量为0.1 μl时,在卷吸速度ue=1.0 mm/s时,油池呈现为“葫芦”形;随卷吸速度增加,油池的前端变得突出,尾部逐渐收缩并形成尾部液滴,油池外形演变为“蝌蚪”形,此过程中出口气穴长度l1逐渐增加,在卷吸速度ue=2.8 mm/s时达到临界状态;在这一卷吸速度下尾部油池被气穴分离为两部分,定义该速度为油池分离临界速度,用ue*表示. 随着供油量的增加,可以看出油池尺寸整体扩大,且在ue=1.0 mm /s时尾部油池明显大于其前端油池,但随着卷吸速度增加,尾部油池仍然逐渐收缩. 观察气穴长度变化不难发现,尽管油池外形随供油量增加而扩大,但在相同卷吸速度ue下,出口气穴长度基本相等. 所不同的是尾部油池分离临界速度ue*逐渐增大,供油量为0.1、0.2、0.3和0.4 μl时油池分离临界速度分别为2.8、6.2、7.0和8.0 mm/s. 图2(a)表明,在相同的卷吸速度下,供油量几乎未对气穴长度产生影响,但决定了尾部气穴分离的临界速度.
为了进一步观察尾部油池在临界卷吸速度下的分离过程,图2(b)给出了供油量为0.1和0.4 μl时尾部油池分离的动态过程图. 为了便于表示油池的瞬态变化,定义卷吸速度刚到达临界速度的时刻为t0 (参考时间). 油量为0.1 μl条件下,在t=t0+1.333 s时刻,油池的初始尾部液滴断开后脱离主油池,此后初始液滴的残余部分在主油池尾部收缩成次液滴. 但残余部分的快速收缩导致了次液滴不稳定,在t=t0+1.5 s时发生明显的抖动,并在t=t0+1.667 s时刻与主油池脱离,此时尾部油池也被气穴打开,分离开的油池尾部也均形成液滴. 在油量为0.4 μl条件下,当速度达到8 mm/s时,初始尾部液滴迅速脱离主油池,与0.1 μl条件下不同的是残余部分并未形成明显的次液滴(t=t0时刻),而是逐渐融入到主油池中,且主油池尾部边界逐渐逼近气穴边界;在t=t0+2.333 s时刻,残余部分几乎消失;在t=t0+2.367 s时刻,主油池边界被气穴边界击穿,尾部油池分离为两部分. 需要说明的是,尾部液滴起布油作用,闭合状态下的尾部液滴将润滑油布于滚道上,分离状态下尾部的两个液滴将润滑油布于滚道两侧,形成油脊. 两种布油方式对入口区供油产生了显著影响,并决定了接触区润滑状态.
为了定量分析不同供油量下气穴长度随卷吸速度的变化,定义无量纲气穴长度为λ1:
λ1=l1a (1) 式中a为赫兹接触区半径. 图3给出了四种供油量下气穴无量纲长度λ1随卷吸速度ue的变化曲线. 可以看出,气穴无量纲长度λ1随着卷吸速度的增大而近似线性增加,且各曲线也接近重合,即相同卷吸速度、不同供油量下的无量纲气穴长度基本相等,进一步说明了供油量对气穴长度的影响不明显. 图中的4个坐标值为尾部油池分离临界速度ue*及其所对应的气穴长度,可见临界气穴长度随着供油量的增加而增加.
2.2 供油量对侧油池宽度的影响
上述图2中观察了较低速度条件下油池从闭合到尾部油池打开的过程,为了观察速度进一步增大后接触区两侧油池形态的变化,图4给出了PAO10在纯滚动的条件下,供油量分别为0.1、0.2、0.3和0.4 μl时侧油池随速度变化图像,图中侧油池用白色虚线标出. 可以看出,在供油量为0.1 μl的条件下,当卷吸速度为ue=4 mm/s时,尾部油池打开,油池呈现为经典的“蝴蝶”形;在同一卷吸速度下,供油量为0.2、0.3和0.4 μl时所对应的油池仍为“葫芦”形. 当卷吸速度为ue˃16 mm/s时,四种供油量所形成的油池形状均为“蝴蝶”形. 当供油量为0.1 和0.2 μl时,随卷吸速度增加,两侧油池尺寸逐渐减少;而当供油量为0.3和0.4 μl时,随卷吸速度增加,两侧油池的油量基本不变.定义侧油池无量纲宽度λ2用以衡量侧油池的大小.
λ2=l2a (2) 图5给出了四种供油量下侧油池无量纲宽度λ2随卷吸速度ue的变化曲线. 可以看出,当卷吸速度ue≤16 mm/s时,随着卷吸速度的增加,四种供油量曲线均有1个上升的过程;当卷吸速度ue˃16 mm/s,供油量为0.1和0.2 μl时,随着卷吸速度的增加,曲线有下降的趋势;供油量为0.3和0.4 μl时,随着卷吸速度增加,曲线比较平缓.
3. 机理分析
3.1 近接触区油池的形成与变形
近接触区油池的形成是因为润滑剂在毛细力作用下充满钢球与玻璃盘接触区附近的微小楔形间隙,对于球-盘静态接触条件下毛细力[20]可表示为式(3).
Fcap=ΔpA=γr0A=2πRγ(cosθ1+cosθ2) (3) 式中:Δp为油池内外压差,A为油池面积,γ为气液界面张力,r0为球盘间隙内的液桥半径,R为钢球半径,θ1为润滑剂与玻璃盘表面的接触角,θ2为润滑剂与钢球表面的接触角. 尽管该毛细力的大小随油池的形态不同而发生变化,但始终是油池形成的主导机制. 如图6所示,当球、盘两表面运动时,润滑剂在表面卷吸作用下流动,入口区液桥半径减小,而出口区液桥半径增大,即
r0=r22R(cosθ1+cosθ2) (4) 式中:r为接触区中心到油池边界的距离,显然,液桥半径随r的增大而增大. 由(3)式可知,相对于静态接触,油池边界上的压差Δp在入口区域增大,而在出口区域减小,这种差别也导致油池外形沿着卷吸方向伸长.
为了便于分析油池外形随速度的变化,图6中以1 mm/s下的两个相交的(虚线)圆为参考,对比了不同速度下的油池外形. 由于出口气穴区内的压力小于环境压力,引起了图6中M、N和Q处的油池颈缩. 随着速度增加,气穴区逐渐扩展,颈缩处油池宽度明显减小且位置逐渐下移,在颈缩区下方形成尾部液滴. 显然,速度进一步增加,颈缩处油池宽度将减小到难以承受尾部液滴的拉力,导致尾部液滴脱离主油池. 同样受到尾部油池气穴区内外压差作用及其布油的影响,以及入口区的压力梯度作用,导致入口区头部油池凸出.
3.2 近接触区油池形态与气穴长度
油池形态随速度增加可分为图7所示的3种形态:闭合态、半开放态和分离态. 其中,闭合态对应于富油润滑状态,取决于供油量;半开放态(仅尾部油池分离)对应于中等乏油状态或者富油润滑状态;分离态对应于严重乏油状态或干涸润滑状态.
闭合态油池气穴区压力小于大气压力(pcav<pamb,即负压 ),甚至小于饱和压力. 由于压差
Δp′=pamb−pcav 的量级明显小于接触区压力,常在数值计算中将气穴区压力设定为大气压力. 但负压[图7(b)]是气穴产生根本原因,一方面使溶解在润滑剂中的空气释放出来,另一方面也使润滑油内的气体向气穴区聚集,使气穴区域扩大[21]. 闭合态油池内气穴长度受到多种因素影响,Stadler等[22]通过数值计算回归出了预测气穴长度的经验公式,如式(5)所示.la=9.2G0.3(2U)0.2775W−0.27k1k2 (5) 式中:G、U和W分别为无量纲材料、速度和载荷参数[23],
k1=(pgaugepcav)0.32 和k2=(αrefα)0.3 为修正系数. 其中pgauge 和αref 分别为数值计算使用的气穴压力和润滑剂黏压系数,可视为参考量. 但采用式(5)及图4中的参数,计算得到的气穴长度低于测量值,且随着速度增加差距明显增大,如图8所示. Emden等[19,24]的计算结果也明显低于试验测量值,其主要原因是在计算过程中并未考虑气穴区的表面张力作用. 若考虑表面张力作用以及液体的密度,气穴区各参数应符合式(6).F(l,a,u,η,Δp′,ρ,γ)=0 (6) 采用π定理量纲分析方法对(6)式进行量纲分析,可得:
la=F(ρuaη,aΔp′uη,γuη) = F(Re,CaRe2,ReWe) (7) 式中:Re、Ca和We分别为雷诺数、毛细数和韦伯数. 上式可进一步表示为式(8)[19].
la=c1(Re)b(CaRe2)c(ReWe)d + c2 (8) 式中各系数需通过试验确定. 若采用Emden等[19]的数值b=1/4,c=-5/4和d=1/2,可得到如图9所示的气穴长度与无量纲参数组良好线性关系.
闭合态油池尾部位于滚道中间,可将润滑剂布置于滚道上[图7(c) State 1],有效保证了滚动体的入口区供油,使接触副处于富油润滑状态下.
当卷吸速度达到临界卷吸速度时,尾部油池被气穴区分离,此时气穴区压力与大气压力相等[pcav= pamb, 图7(b)中State 2],油池为半开放态. 处于该状态下的油池,分离开的尾部油池将润滑剂布置于滚道两侧形成油脊[图7(c)中State 2],同时,进入接触区的润滑剂在压力梯度作用下流向滚道两侧(侧泄). 两表面油脊在接触区附近再次重叠,并在挤压作用下铺展,铺展的面积与供油量相关. 另一方面,两侧油脊可在表面张力作用下向滚道回填,其回填量
hrep 可表示为式(9)[13].hrep=C0γt2η (9) 式中:C0为回油系数,
C0∝a/ahsidehside ,即回油系数由接触半径和油脊高度(供油量)的比值决定;t为相继运动滚动体的时间间隔,由滚动体运动速度及其数量决定. 两侧油池的铺展与回填量决定了入口区的供油状态,当两侧油脊的回填量不足时,润滑剂不能完全充满入口间隙,使入口油池呈现内凹的外形[图7(a)中State 2],其入口油池边界与接触区边缘距离S决定了接触区的油膜厚度. 由图4中ue=16 mm/s时的干涉图可见,在较大供油量0.4 μl时,因铺展面积和回填量较大,其入口距离大于0.3 μl工况时.当卷吸速度较大时,两侧油脊的回填量进一步减小,入口距离S=0时,入口油池与接触区边缘接触,使入口油池也分离,油池处于分离态[图7(a)中State 3]. 当供油量较小时,随着速度增加,润滑剂的侧泄量大于回填量,使侧油池宽度λ2持续减小;当供油量较大时,润滑剂的侧泄与回填达到平衡,使侧油池宽度λ2趋向定值,进而呈现出图5所示的变化趋势.
4. 结论
试验中观察了纯滚动工况不同供油量时近接触区油池随卷吸速度的演化过程,得到如下结论:
a. 近接触区油池在低速、中速和高速时分别呈现出闭合态、半开放态和分离态三种状态.
b. 油池为闭合态时,出口气穴长度随着速度的增加而增加,尾部油池分离临界速度随着供油量的增加而增加;油池为半开放态时,入口距离随着卷吸速度的增加而减小,随着供油量的增加而增加;油池为分离态时,不同供油量时润滑剂处于不同供-失平衡状态,侧油池宽度呈现出不同变化趋势.
c. 在不考虑气穴区表面张力作用的条件时,气穴长度的计算数值低于试验测量值;而考虑表面张力作用时,气穴长度与无量纲参数组呈现出线性关系.
d. 因球-盘接触压力、接触几何特征及供油参数与真实轴承有所差别,其油池形态和润滑状态也与真实轴承的润滑特征存在差异.
-
表 1 试验条件
Table 1 Experimental conditions
Parameters Specifications Entrainment speed, ue/(mm/s) 0.2~400 Load/N 10 Maximum Hertz pressure/GPa 0.234 Lubricant amount/μl 0.1; 0.2; 0.3; 0.4 Ambient temperature/℃ 20±1 表 2 试验用润滑油性质
Table 2 Properties of lubricants used in experiments
Lubricant Viscosity at 20 ℃
/(Pa·s)Density at 20 ℃
/(kg/m3)PAO10 0.13 835 -
[1] 刘维民, 许俊, 冯大鹏, 等. 合成润滑油的研究现状及发展趋势[J]. 摩擦学学报, 2013, 33(1): 91–104 doi: 10.16078/j.tribology.2013.01.002 Liu Weimin, Xu Jun, Feng Dapeng, et al. The research status and prospect of synthetic lubricating oils[J]. Tribology, 2013, 33(1): 91–104 doi: 10.16078/j.tribology.2013.01.002
[2] Liang He, Guo Dan, Luo Jianbin. Experimental investigation of lubrication film starvation of polyalphaolefin oil at high speeds[J]. Tribology Letters, 2014, 56(3): 491–500. doi: 10.1007/s11249-014-0426-9
[3] Wikström V, Höglund E. Starting and steady-state friction torque of grease-lubricated rolling element bearings at low temperatures-part I: a parameter study[J]. Tribology Transactions, 1996, 39(3): 517–526. doi: 10.1080/10402009608983562
[4] Li X M, Guo F, Poll G, et al. On the Grease Film Evolution in Rolling EHL Contacts[J]. Friction, 2021, 9(1): 179–190. doi: 10.1007/s40544-020-0381-4
[5] Jauhari K. Oil lubrication on high-speed spindle bearing system: a review[J]. Proceeding of Tribology Industry, 2014, 16: 216–231.
[6] Wedeven L D, Evans D, Cameron A. Optical analysis of ball bearing starvation[J]. Journal of Lubrication Technology, 1971, 93(3): 349–361. doi: 10.1115/1.3451591
[7] Poll G, Li X M, Bader N, et al. Starved lubrication in rolling contacts-a review[J]. Bearing World Journal, 2019, 4: 69–81.
[8] Lee S W, Maeda T. Angular contact ball bearings that require a minimum quantity of grease supply for lubrication[J]. NTN Technical Review, 2008, 76: 90–95.
[9] Qian Shanhua, Guo Dan, Liu Shuhai, et al. Experimental investigation of lubricant flow properties under micro oil supply condition[J]. Journal of Tribology, 2012, 134(4): 041501. doi: 10.1115/1.4007107
[10] 王少朋. 微量供油条件下润滑油膜特性实验研究[D]. 青岛: 青岛理工大学, 2014 Wang Shaopeng. Experimental studies on the lubrication films with minimal quantity oil feeding[D]. Qingdao: Qingdao Technological University, 2014 (in Chinese)
[11] 刘成龙, 栗心明, 郭峰, 等. 单个微油滴弹流润滑行为的试验研究[J]. 摩擦学学报, 2017, 37(3): 340–347 doi: 10.16078/j.tribology.2017.03.009 Liu Chenglong, Li Xinming, Guo Feng, et al. Experimental observation of EHL behavior under single oil droplet supply[J]. Tribology, 2017, 37(3): 340–347 doi: 10.16078/j.tribology.2017.03.009
[12] 李书义, 郭峰, 栗心明. 微油滴供油条件下弹流润滑行为的数值分析[J]. 摩擦学学报, 2016, 36(4): 413–420 doi: 10.16078/j.tribology.2016.04.003 Li Shuyi, Guo Feng, Li Xinming. Numerical analysis of elastohydrodynamic lubrication behavior under micro oil droplet supply[J]. Tribology, 2016, 36(4): 413–420 doi: 10.16078/j.tribology.2016.04.003
[13] Chiu Y P. An analysis and prediction of lubricant film starvation in rolling contact systems[J]. ASLE Transactions, 1974, 17(1): 22–35. doi: 10.1080/05698197408981435
[14] Li X M, Guo F, Wong P L, et al. Regulation of lubricant supply by wettability gradient in rolling EHL contacts[J]. Tribology International, 2018, 120: 565–574. doi: 10.1016/j.triboint.2018.01.020
[15] 李哲, 刘成龙, 李超, 等. 条状润湿表面对面接触限量润滑油膜的影响研究[J]. 摩擦学学报, 2020, 40(5): 601–607 doi: 10.16078/j.tribology.2020003 Li Zhe, Liu Chenglong, Li Chao, et al. Influence of strip-shaped wet surface on lubrication films in slider-on-disc conformal contact under limited lubricant supply[J]. Tribology, 2020, 40(5): 601–607 doi: 10.16078/j.tribology.2020003
[16] 韩兵, 王文中, 赵自强. 低速下润滑接触区补充供油机制的研究[J]. 摩擦学学报, 2016, 36(3): 341–347 doi: 10.16078/j.tribology.2016.03.011 Han Bing, Wang Wenzhong, Zhao Ziqiang. Oil replenishment mechanism of lubricated contact at low speed[J]. Tribology, 2016, 36(3): 341–347 doi: 10.16078/j.tribology.2016.03.011
[17] 江楠, 栗心明, 杨萍, 等. 滑滚条件下定量润滑状态与润滑剂回填效应试验观察[J]. 摩擦学学报, 2021, 41(1): 9–16 doi: 10.16078/j.tribology.2020045 Jiang Nan, Li Xinming, Yang Ping, et al. Observation of lubrication state and lubricant replenishment under sliding-rolling motions with single-charging amount[J]. Tribology, 2021, 41(1): 9–16 doi: 10.16078/j.tribology.2020045
[18] 梁鹤, 张宇, 王文中. 轴承内部润滑油分布及回流的试验观察与研究[J]. 摩擦学学报, 2020, 40(4): 450–456 doi: 10.16078/j.tribology.2019180 Liang He, Zhang Yu, Wang Wenzhong. Experimental observation and investigation of oil distribution and replenishment in a rolling bearing model[J]. Tribology, 2020, 40(4): 450–456 doi: 10.16078/j.tribology.2019180
[19] van Emden E, Venner C H, Morales-Espejel G E. Aspects of flow and cavitation around an EHL contact[J]. Tribology International, 2016, 95: 435–448. doi: 10.1016/j.triboint.2015.11.042
[20] Popov V L. Contact mechanics and friction-physical principles and applications[M]. Berlin: Springer, 2010.
[21] Braun M J, Hannon W M. Cavitation formation and modelling for fluid film bearings: a review[J]. Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, 2010, 224(9): 839–863. doi: 10.1243/13506501jet772
[22] Stadler K, Izumi N, Morita T, et al. Estimation of cavity length in EHL rolling point contact[J]. Journal of Tribology, 2008, 130(3): 031502. doi: 10.1115/1.2919780
[23] Chiu Y P. Discussion: “isothermal elastohydrodynamic lubrication of point contacts: part II-ellipticity parameter results[J]. Journal of Lubrication Technology, 1976, 98(3): 381–382. doi: 10.1115/1.3452862
[24] van Emden E, Venner C H, Morales-Espejel G E. A challenge to cavitation modeling in the outlet flow of an EHL contact[J]. Tribology International, 2016, 102: 275–286. doi: 10.1016/j.triboint.2016.05.039