表面凹槽对流体动压润滑油膜厚度的影响

利用自行研发的面接触光干涉油膜厚度测量系统,对表面凹槽滑块的流体动压润滑油膜厚度进行了试验测量,试验中以静止的微型凹槽滑块平面和旋转的光学透明圆盘平面构成润滑副,且两润滑平面始终保持平行;在固定的载荷(速度)条件下,对油膜厚度-速度(载荷)曲线进行测量.结果表明：凹槽的宽度,深度,方向和位置等因素对油膜厚度有着重要影响.同时采用经典Reynolds方程对油膜厚度进行了理论计算,结果表明理论值在某些条件下并不能解释试验结果.     

超薄膜磁头滑块气动力特性

采用有限差分法对广义润滑方程进行数值求解,计算出计算机磁头滑块压强场的分布情况。分析、研究了其稳态和动态气动力特性,并将计算结果分别与求解一阶、二阶修正雷诺方程所得到的结果进行了比较,得到如下三个结论:(1)当飞行高度很小,飞行速度较低时,必须采用广义润滑方程进行磁头滑块的气动力计算,(2)与广义润滑方程结果比较,求解一阶修正雷诺方程所得到的计算结果总是偏高,而求解二阶修正雷诺方程所得到的计算结果总是偏低。此外,还解决了大压缩数下数值失稳问题,使得压缩数可以计算到120万,足以适应任何实际工程的需要。     

纳米薄膜润滑物理—数学模型及数值分析

基于润滑剂分子通常具有链状结构的事实，在分析润滑剂分子链长同膜厚关系的基础上，建立了纳米薄膜润滑物理模型，并利用含旋转量的流体力学运动方程得到了相应的Reynolds方程，同时对薄膜润滑Reynolds方程进行了数值计算，以考察特征长度对薄膜润滑状态参数的影响。结果表明，同相应的厚膜解相比，薄膜模型中润滑剂的粘度及承载能力均明显提高，且其提高幅度随着特征长度的增大而增大。根据润滑剂分子链长度确定的薄膜润滑区间以及膜厚－速度关系数值解同相应的试验结果一致。     

Numerical analyses of rectangular micro-textures in hydrodynamic lubrication regime for sliding contacts

Hydrodynamic lubrication of rectangular micro-textures on sliding contact surfaces is investigated using numerical calculation methods. The theoretical models for the slider surface are developed and the film pressure is used to evaluate the hydrodynamic lubrication based on the Reynolds equation. Meanwhile, the geometry and distribution of the rectangular dimples are optimized for maximizing the average film pressure. Results show that the film pressure is dependent on the geometry and distribution of the rectangular micro-dimples. The optimal geometry of the single rectangular dimple is obtained, and the spacing has an important influence on the film pressure. The distribution types of rectangular dimples affect the hydrodynamic lubrication significantly and the interlaced array of the rectangular micro-dimples is beneficial to enhancing the hydrodynamic lubrication. Meanwhile, the rectangular dimples with 72° interlaced angle exhibits the best effectivity.

     

Elasto-hydrodynamic lubrication: A non-linear complementarity problem

The classical problem of elasto-hydrodynamic lubrication of cylinders in line contact is formulated as a non-linear complementarity problem. A direct algorithm is applied to the approximation obtained by finite differences. Implementation considerations are emphasized. The new method provides reliable and automatic location of the previously troublesome lubricant free boundary. Numerical results reveal the qualitative behaviour of the pressure distribution and the lubricant film thickness under variation of key physical parameters.     

Combined effect of viscosity variation and non-Newtonian Rabinowitsch fluid in wide parallel rectangular-porous plate with squeeze-film characteristics

The present article carries out the study of viscosity variation of non-Newtonian fluid with the homogeneous porous wall on wide parallel rectangular-plate based on the Rabinowitsch fluid model. The non-linear modified Reynolds equation is derived for the lubrication of rectangular squeeze film bearing with viscosity variation and porous parameter. Using the Morgan–Cameron approximation, the nonlinear Reynolds-type equation for squeeze-film which governs the film pressure is solved within the fundamentals of small perturbation technique. The characteristic of the wide parallel rectangular-porous plate is numerically computed for different physical quantities such as film pressure, load carrying capacity and response time. Moreover, as limiting cases some of the results from the available literature are recovered also. Further, the findings reveal that the viscosity variation of non-Newtonian fluid and the presence of porous wall lead to reduction in the load capacity and the response time respectively. Here, the porous matrix consists of a system of capillaries of very small radii with the homogeneous porous wall. The impact of porosity is incorporated as a result it acts as self-lubrication on bearing surface. Also, the effect of viscosity variation is one of the most important characteristics of fluid which helps in the design of bearings for lubrication in engineering and industrial applications.

     

计及轴颈倾斜的径向滑动轴承湍流润滑分析

分析了轴颈倾斜状态下,径向滑动轴承的湍流润滑性能. 基于轴颈倾斜的统一Reynolds方程和能量方程,应用有限差分法求解了不同轴颈倾斜方位角、轴颈倾斜度、偏心率和平均雷诺数下的径向滑动轴承湍流润滑性能. 结果表明:轴颈倾斜方位角α=0°时,随着轴颈倾斜度的增大,轴承油膜的压力峰向轴承一端移动,轴承一端的轴向油膜温度梯度增大;α=90°时,随着轴颈倾斜度的增大,轴承油膜压力逐渐出现双峰分布,且向轴承两端移动,轴承两端的轴向温度梯度也不断增大. 在相同轴颈倾斜度增量下,轴承最大油膜压力、最高油膜温度、承载力和稳定工作力矩的增量随轴承中央截面偏心率的增大而增大. 相同轴颈倾斜度增量下,轴承最大油膜压力增量、最高油膜温度增量、承载力增量、摩擦系数减量和稳定工作力矩增量随平均雷诺数的增大而增大. 可见,径向滑动轴承湍流润滑分析中有必要考虑轴颈倾斜因素的影响.      

Convective inertia effects in wall-bounded thin film flows

We apply the boundary layer equations to inertial flow in wall bounded films that might be characterized as ‘thin’, say ɛ ≤ 0.1 where ɛ is the ratio of the characteristic lengths, yet to which the lubrication approximation of Reynolds no longer applies. Two particular flow geometries are investigated, nominally parallel plates and nominally inclined plates, both with and without spatially periodic perturbation of the stationary plate. A Galerkin-B spline formulation of the governing equations is employed, and we rely on parametric continuation to obtain solutions at higher values of the Reynolds number. In particular, we are able to demonstrate that the boundary layer equations yield accurate results for a wide range of Reynolds numbers when the aspect ratio is less than 1/10. We also find that in both nominally parallel and nominally inclined geometries the sign of the inertial force correction is determined by the film contour in the neighborhood of the exit, this result might have implications in the design of MEMS devices.     

斜齿轮弹流润滑下的接触疲劳寿命计算

经典齿轮接触疲劳强度理论是基于光滑表面赫兹干接触理论,而实际齿面具有粗糙度,且啮合轮齿多数处于混合润滑状态.本文基于齿轮润滑接触分析建立了渐开线斜齿轮的接触疲劳寿命计算模型.模型由齿轮润滑接触分析模型和基于次表面应力分布的疲劳寿命模型组成.首先将斜齿圆柱齿轮一对齿的瞬时啮合等效为两反向圆锥的接触问题,建立了齿轮的有限长弹流润滑计算模型,考虑了齿轮啮合周期内瞬时载荷、接触线长、卷吸速度等因素的影响,基于统一雷诺方程方法求得啮合齿对间的润滑压力和油膜厚度分布;在此基础上,计算轮齿接触区次表面的米歇斯应力分布,根据Zaretsky接触疲劳寿命计算模型,对齿轮组的接触疲劳寿命进行模拟预测.针对不同工况参数下接触疲劳寿命计算表明:润滑油黏度、轮齿表面粗糙度等因素对齿面接触疲劳寿命均有显著的影响.     

粗糙度纹理对有限长线接触混合润滑影响

采用统一Reynolds方程建立有限长线接触混合润滑模型,研究横向、纵向和二维规则表面粗糙度的波长、幅值及工况变化对润滑影响.结果表明:波长、幅值与工况对三种表面粗糙度接触副的润滑影响类似;随着载荷增大,平均膜厚降低,摩擦系数、接触载荷比与接触面积比均增大;随着转速升高,平均膜厚增大,摩擦系数、接触载荷比与面积比均降低,其中摩擦系数随转速进一步增大而小幅升高.在润滑状态转换区域润滑特征参数变化显著,而其他润滑区域变化平缓.沿卷吸速度方向的压力与膜厚波动分布存在相位差,垂直方向则同相位;相同的工况和粗糙度参数时,纵向粗糙度分布更有利于接触润滑.     

球轴承启停过程的瞬态热混合润滑分析

建立了角接触球轴承的几何和数学模型,通过求解考虑了热效应和时变效应的Reynolds方程,对启动和制动过程中的球轴承瞬态热混合润滑问题进行了分析,考虑了不同加速度启动工况下的瞬态热混合润滑情况.结果表明:启动过程中,随转动速度的增大,最小膜厚增大,轴承逐渐由边界润滑进入弹流润滑状态;不同滑滚比下进入弹流润滑状态的时间有所不同,随着滑滚比的增大,进入弹流润滑的时刻有所推迟,轴承处于同一转速条件下的油膜厚度变小;随着转速的增大,油膜温度升高,最高油膜温度增长幅度减小;加速度的增大使边界润滑消失的时间提前,随着转速的增加,油膜温度增大,且在同一时刻加速度越大油膜温度越高;油膜减小过程中的挤压膜作用导致轴承制动过程中的油膜厚度大于启动过程中的油膜厚度;由于在相同转速下轴承在启动时处于边界润滑状态,而在制动时处于弹流润滑状态,润滑状态的不同导致制动过程中的最高油膜温度较启动过程较小.     

粗糙表面滑动轴承非牛顿介质润滑的计算

分析了轴承润滑中非牛顿介质的流变特性以及轴承表面形貌对润滑效果的影响,推导出了修正的雷诺方程,在方程中用差分粘度和第一正应力差函数表征非牛顿介质的流变特性,采用流量因子表征表面形貌的作用.结果表明:差分粘度的变化是影响润滑结果的主要因素,它取决于不同介质的动态参数和剪切频率范围,非牛顿介质润滑的承载力并非总高于或低于牛顿介质;在动载荷条件下,第一正应力差效应使油膜压力明显增大;影响流量因子的2个因素为表面形貌的粗糙度和纹理,但其作用远小于非牛顿介质流变特性的影响.     

供油量对点接触表面微织构润滑性能的影响

采用球盘接触模型,以入口油层厚度为输入变量,基于统一Reynolds方程法建立润滑分析模型,研究了供油量对点接触表面微织构润滑性能的影响.数值计算结果显示:供油量显著影响接触区油膜厚度;微织构能改善摩擦副的润滑摩擦性能;严重乏油时接触区内微坑会产生气穴现象;瞬态下,微织构在乏油时较充分供油时改善润滑的效果明显.     

考虑滚道表面油层分布的滚动轴承润滑分析

研究表明供油量对弹流润滑性能产生显著影响.滚动轴承中由于离心力和滚动体的反复滚压,滚道表面上的润滑剂呈现出非均匀分布的特点.大多数润滑剂被推挤到滚道的两侧,致使接触区的入口间隙不能被完全充满,导致乏油润滑,滚动体与滚道间接触压力接近于赫兹压力分布,膜厚较全膜润滑有明显的减小.本文基于润滑剂的流量连续建立滚道表面油层厚度分布模型,考虑润滑接触压力的影响,计算滚道上的侧流量以预测轴承滚道上补给油层厚度及形状随时间的变化规律;进而以此作为滚动体和滚道接触区的入口油层厚度,采用统一Reynolds方程法数值模拟计算每个时刻轴承滚道与滚动体之间的润滑油膜厚度,压力分布等参数,分析轴承在点接触乏油条件下运行的润滑性能.     

渐开线斜齿轮非稳态弹流润滑数值模拟研究

建立了渐开线斜齿轮啮合的弹流润滑计算模型,将斜齿圆柱齿轮啮合的齿面接触等效为有限长线接触的弹流润滑问题.考虑斜齿轮啮合的实际因素,将斜齿轮啮合过程中的等效曲率半径和齿面载荷的变化反映到弹流润滑计算模型中,应用统一Reynolds方程方法求得轮齿在1个完整啮合周期内的瞬时弹流润滑数值解.结果表明：斜齿轮啮合线上各点处的膜厚、压力均有较大不同,各接触点处的油膜厚度受综合曲率半径的影响较大;斜齿轮传动非稳态效应相对较弱;小齿轮齿根附近和节点位置处润滑状态较差;适当增大压力角可以改善齿轮的润滑.     

聚合物注射成型流动残余应力的数值分析

建立了可压缩黏弹性聚合物熔体在薄壁型腔中充模/保压过程中非等温、非稳态流动 的数学模型,用数值方法实现了注射成型过程中流动应力和取向建立及松弛过程的模拟,研 究了熔体温度、模具温度和注射速率等工艺条件对分子冻结取向的影响,取得了与实验相符 的结果.     

Operability windows in viscoelastic slot coating flows using a simplified viscoelastic-capillary model

The viscoelastic-capillary model to predict approximately coating windows for the stable operations of viscoelastic coating liquids is derived using a lubrication approximation in slot coating processes. Pressure distributions and velocity profiles for viscoelastic liquids based on the Oldroyd-B and Phan-Thien and Tanner (PTT) models are solved in the coating bead region considering the Couette-Poiseuille flow feature and the pressure jumps at upstream and downstream menisci. Practical operating limits for the uniform coating of rheologically different liquids that are free from leaking and bead break-up defects are constructed under various conditions, incorporating the position of the upstream meniscus as an important indicator while determining limits. The shift of the uniform operating range shows different patterns for the Oldroyd-B liquid with a constant shear viscosity and the PTT liquid with a shear-thinning nature in comparison with the Newtonian case. The windows predicted by the simplified model are corroborated with experimental observations for one Newtonian and two viscoelastic liquids.     

变阶梯结构自适应径向滑动轴承的研究

分析了变阶梯结构自适应径向滑动轴承膜压力的形成机理,并建立了二段油膜形耦合变形阶梯结构的流量控制方程,采用有限差分法交叉循环迭代求解了油压力的雷诺方程和变阶梯结构的流量控制方程,比较了设计参数对这咎径向滑动轴承的最小油膜厚度、压力分布、承载能力、摩擦阻力等性能和温升的影响,研究结果表明;合理地选择设计参数可以使得这种滑动轴承具有较好的润滑性能和承载特性。     

A solution to the hydrodynamic lubrication of a circular point contact sliding over a flat surface with cavitation

This letter presents an analytical solution to the hydrodynamic lubrication of a circular point contact sliding over a flat surface with cavitation. The solution is found by solving the Reynolds equation with Reynolds boundary condition for cavitation. The cavitation boundary is shown to be straight lines directed 108.4° against the sliding direction. The result is experimentally verified in the limit of large values of viscosity, sliding velocity and radius of a spherical ball. The solution raises questions about the coupling between cavitation and film rupture and can be used as an independent check on the validity of numerical solutions.