滑动轴承模拟器 返回
三体学

滑动轴承模拟器

基于短轴承近似雷诺兹方程实时计算期刊轴承的油膜压力分布、承载能力、最小油膜厚度。直观展示索默费尔德数与轴承性能的关系。

轴承尺寸

运行条件

计算结果
Sommerfeld S
离心率 ε
h_min (μm)
P_max (MPa)
摩擦系数 f
动力损失 (W)
轴承截面图
周向油膜压力分布
索默费尔德数 S vs 最小油膜厚度 h_min/C
理论与主要公式

$$h(x) = h_1 + (h_2 - h_1)\frac{x}{L}$$

楔形油膜厚度(m):$h_1, h_2$ 是两端的膜厚,$L$ 是滑动面长度(m)。

$$\frac{dp}{dx} = 6\mu U \frac{h - \bar{h}}{h^3}$$

雷诺兹方程(一维):$\mu$ 是粘度(Pa·s),$U$ 是相对速度(m/s),$\bar{h}$ 是最大压力位置的膜厚。

$$F_L = \int_0^L p(x)\,dx$$

流体浮力(N/m):油膜压力分布积分后的轴承荷载支持能力。

滑动轴承模拟器简介

🙋
这个模拟器里经常出现"索默费尔德数",这是什么?在设计中怎么用?
🎓
简单地说,这是表示"轴承浮起程度"的无量纲指标。公式是 $S = \frac{\mu N L D}{W}\left(\frac{R}{C}\right)^2$,表示粘性力与荷载的比值。S越大(高速、高粘度、低荷载)油膜越厚越稳定,S越小就容易接触摩擦。实际设计中,通过这个值来预测润滑状态。试试调整上面的"转速"和"轴承荷载"滑块,看看S值怎么变化。
🙋
原来如此!但是如果油膜太薄会怎样?屏幕上显示的"最小油膜厚度"是什么意思?
🎓
这是最关键的!最小油膜厚度 $h_{min}= C(1 - \epsilon)$ 就是轴和轴承最接近的距离。如果这个值不够厚(小于表面粗糙度),金属就会接触,导致烧伤和磨损。实务中的设计准则是让 $h_{min}$ 至少是表面粗糙度的3~5倍。在模拟器里调整"径向间隙"变小,或者"离心率"变大,你就能看到 $h_{min}$ 急剧下降,变得很危险。试试看。
🙋
明白了…。屏幕上五彩的压力分布图在实际设计中怎么看?最红的地方是最危险的吗?
🎓
很好的观察!那个压力分布正是轴承支撑荷载的原理。油被挤进楔形间隙产生高压,这个压力把轴撑起来。红色高压区的位置和大小决定了承载能力和轴的位置。现场最常见的问题就是压力分布偏斜导致"边缘负荷",某处压力过高。在模拟器里把"轴承宽度"调得很小,你会看到压力峰值变得很尖锐,这就是边缘负荷的效果。

物理模型和主要公式

本模拟器的基础是短轴承近似的雷诺兹方程,计算轴承圆周方向的油膜压力分布。

$$ \frac{d}{d\theta}\left( h^3 \frac{dp}{d\theta}\right) = 6 \mu U R \frac{dh}{d\theta}$$

其中,$h = C(1 + \epsilon \cos\theta)$ 是油膜厚度,$p$ 是油膜压力,$\mu$ 是油的粘度,$U$ 是轴表面速度,$R$ 是轴半径,$\theta$ 是圆周角度,$C$ 是径向间隙,$\epsilon$ 是离心率。右边的 $dh/d\theta$ 表示"楔形效应",这是产生压力的来源。

求解雷诺兹方程得到压力分布后,积分可得轴承能支持的荷载(承载能力)和摩擦系数,这些是设计的直接指标。

$$ W = L \int_{0}^{2\pi}p R \cos\theta \, d\theta, \quad f = \frac{F_f}{W} $$

$W$ 是承载能力,$L$ 是轴承宽度,$F_f$ 是粘性阻力产生的摩擦力,$f$ 是摩擦系数。索默费尔德数 $S$ 就是这些结果的无量纲化形式,是重要的设计参数。

常见问题

索默费尔德数是描述轴承性能的无量纲数,定义为 S = (μN/P) × (R/C)²,其中 μ 是粘度,N 是转速,P 是面压,R 是轴半径,C 是间隙。这个值越大,油膜越容易形成,是评估轴承负荷容量和稳定性的基本指标。
离心率ε(0~1)表示轴的偏心程度。增大ε会使最小油膜厚度减少,油膜压力峰值升高。承载能力增加,但油膜变薄,烧伤风险升高。ε接近0时油膜均匀稳定。
短轴承近似适用于轴承宽度相对直径较小(L/D < 0.5)的情况。此时可忽略周向压力梯度,只考虑轴向流动。对于宽度大的轴承(L/D > 1),需用长轴承近似或二维分析。本工具采用短轴承近似。
增大粘度μ或转速N会增大雷诺兹方程右侧的楔形效应,油膜压力升高,承载能力改善。但过度增加会导致发热和气穴风险,需在适当范围内调整。

实际应用

大型发电用汽轮机和发电机:巨大的高速转子需要滑动轴承支撑。通过模拟计算索默费尔德数和最小油膜厚度,确保在启动、停止和全部运行工况下都保持流体润滑状态,是长期可靠性的关键。

船舶推进轴轴承:在海水腐蚀环境下承受推力的滑动轴承。根据推力和螺旋桨转速预测油膜压力分布,优化轴承宽度和间隙,是设计的核心内容。

工业电动机和泵:工厂常用中等规模电动机的轴承设计需平衡成本和性能。通过对比不同润滑油粘度下的摩擦系数和承载能力,既降低能耗又保证耐久性。

汽车发动机曲轴轴承:燃烧压力冲击加高速旋转的极端工况。调整轴承宽度和间隙,评估油膜压力峰值和最小油膜厚度,防止疲劳和烧伤。

常见误区和注意事项

开始使用这个模拟器时,容易踩几个坑。第一个是"索默费尔德数S越大越安全"的误解。虽然S大时油膜厚,但发热变得严重。高速高粘度下,剪切摩擦产生的热导致油粘度下降(热失控),反而油膜变薄。比如转速从1000转升到10000转,发热量增加约10倍。必须配套冷却设计。

第二个是参数输入顺序的问题。不能随意调大"离心率"就算结果,实际工程中离心率是荷载和刚度的结果,而非原因。正确的设计流程是:先根据目标荷载和转速确定粘度和间隙,然后把结果中的离心率和最小油膜厚度当作验证指标来看。

第三是这个模拟器的大前提"无限宽度近似"容易被忽视。现实轴承有侧面漏油的"端部泄漏",所以计算的承载能力通常偏高。特别是宽度很窄(L/D < 0.5)的轴承要特别留意。在模拟器里分别试试L/D=0.3和1.0,你会看到压力分布从尖三角形变成平台形,这就是端部泄漏的表现。

使用指南

  1. 输入轴承半径R(mm)、轴承宽度L(mm)、径向间隙C(μm)。钢制期刊轴承的典型值为R=25mm、L=50mm、C=50μm
  2. 设置转速N(rpm)。电动机轴承通常3000转,机床主轴6000转以上。自动计算索默费尔德数So=(μ·N·R²)/(P·C²)
  3. 用有限差分法离散短轴承雷诺兹方程,实时同步计算油膜压力分布、最小油膜厚度δmin、承载能力W。显示压力分布二维等高线和轴承性能曲线

具体计算例

设R=20mm、L=40mm、C=45μm、N=2800rpm、ISO VG 32矿物油(μ=32mPa·s)的期刊轴承,荷载P=5kN时,索默费尔德数So≈0.95。通过雷诺兹方程分析得最小油膜厚δmin=8.2μm、最大压力pmax=2.8MPa、承载能力W=6.2kN。转速增至5600转时So=1.90,δmin=15.1μm、pmax=1.5MPa。展示了同一荷载下转速升高如何增大油膜厚度的机理

实务注意事项

  1. 间隙C需考虑制造公差和热膨胀。铝轴承座在温度升高时间隙增大,So下降,油膜破坏风险急增,必须在额定荷载基础上留20%安全系数
  2. 粘度输入值应为轴承运行温度下的动粘度。常温32cSt油在80℃时降到约8cSt,必须预测稳态运行温度后再设定模拟条件
  3. 最小油膜厚度小于间隙的0.2倍(C=50μm时低于10μm)会进入边界润滑领域,磨损加速。超过警告值必须更换轴承或调整运行条件