轴承几何参数
运行条件
Sommerfeld数 S
—
偏心率 ε
—
h_min (µm)
—
P_max (MPa)
—
摩擦系数 f
—
功率损失 (W)
—
轴承截面图
周向油膜压力分布
Sommerfeld数 S vs 最小油膜 h_min/C
采用短轴承近似雷诺方程分析滑动轴承的流体润滑。实时计算油膜压力分布、最小油膜厚度和Sommerfeld数。
本模拟器基于“短轴承近似”的雷诺方程，它描述了油膜压力在轴承间隙中的分布。对于长度L远小于直径D的轴承，可以忽略周向压力变化，只考虑轴向变化，使计算大大简化。
$$\frac{\partial}{\partial z}\left( h^3 \frac{\partial p}{\partial z}\right) = 6\mu U \frac{\partial h}{\partial x}$$其中，$p$是油膜压力，$h$是局部油膜厚度，$\mu$是润滑油动力粘度，$U$是轴颈表面速度，$x$是周向坐标，$z$是轴向坐标。这个方程本质是流体质量守恒，左边是压力流动项，右边是剪切（速度）流动项。
轴承性能由两个关键指标衡量：Sommerfeld数和最小油膜厚度。它们直接决定了轴承的运行状态和安全裕度。
$$S = \frac{\mu N L D}{W}\left( \frac{R}{C}\right)^2, \quad h_{min}= C(1 - \varepsilon)$$$S$是Sommerfeld数（无量纲），$N$是轴转速（rps），$W$是载荷，$R$是轴半径，$C$是径向间隙。$h_{min}$是最小油膜厚度，$\varepsilon$是偏心率（轴心偏移量与间隙的比值）。设计目标是保证$h_{min}$大于表面粗糙度的数倍。
大型船舶推进轴系：船舶主机的曲轴轴承和推力轴承承受巨大的交变载荷。CAE模拟用于精确预测不同航速和工况下的油膜压力与厚度，确保在跨洋航行中不会发生因润滑失效导致的“烧瓦”事故。
发电厂汽轮发电机组：汽轮机转子轴承转速极高（3000 rpm以上），要求极其稳定的油膜支撑。模拟分析用于优化轴承的型线（如椭圆轴承、可倾瓦轴承），抑制油膜振荡，保证电网频率的稳定。
内燃机连杆轴承：发动机每个循环中，连杆轴承的载荷大小和方向都在剧烈变化。通过模拟可以分析油膜压力峰值和最小厚度的瞬态变化，为选择机油粘度、设计轴承间隙提供关键依据，减少磨损和摩擦功耗。
高速精密机床主轴：机床主轴轴承要求高精度、低振动。液体动压滑动轴承（如液体静压轴承）能提供极高的刚度和阻尼。模拟用于设计油腔和节流器，确保在高速切削时主轴变形最小，加工精度最高。
开始使用此模拟器时，有几个容易陷入的误区。首先，人们常认为“索末菲数S越大就越安全”，但实际情况并非如此简单。虽然S值增大会使油膜变厚，但发热会变得不可忽视。在高速、高粘度条件下过度增大S值，可能因油液剪切发热导致粘度下降（热失控现象），反而陷入油膜变薄的恶性循环。例如，将转速从1000rpm提升至10000rpm时，发热量几乎增加10倍。冷却设计必不可少。
其次要注意参数输入顺序。若直接设定较大的“偏心率”并点击“计算”，会立即得到显示接触状态的结果。实际工程设计中，应先根据目标载荷与转速确定所需粘度及间隙，再将该状态下的偏心率和最小油膜厚度作为结果进行验证。关键在于理解“偏心率”并非原因，而是载荷与刚度共同作用产生的现象。
最后，切勿忽略此模拟器的重要前提“无限宽度近似”。实际轴承存在侧向油液泄漏的“端泄效应”，因此计算得到的承载能力往往被高估。尤其在宽径比较小（L/D < 0.5左右）的轴承中需特别注意。可在模拟器中对比L/D=0.3与L/D=1.0的情况：宽径比较大时压力分布呈梯形，较小时则呈尖锐三角形——这正是端泄影响的直观体现。
这种滑动轴承流体润滑模拟技术实际上与多领域基础技术直接相关。首当其冲的是摩擦学——研究摩擦、磨损与润滑的科学，轴承正是其典型应用场景。此处涉及的油膜形成原理，可延伸至齿轮齿面、凸轮-从动件等所有“运动机械元件”的接触面设计。
另一关联领域是旋转机械振动工程。滑动轴承油膜不仅支撑轴系，还兼具弹簧与阻尼器功能。油膜的刚度与阻尼特性直接影响轴的临界转速（共振转速）和不稳定振动（油膜振荡）的产生。在模拟器中逐步提升“转速”时，可观察到轴心轨迹的变化，这正是振动特性演变的微观呈现。
此外与CFD（计算流体力学）也密切相关。本工具求解的雷诺方程，是将纳维-斯托克斯方程应用于油膜这种薄层领域后大幅简化的结果。在CFD专家眼中，这堪称“高度专门化的CFD求解器”。因此理解其计算原理，也是学习更复杂流动数值模拟的重要起点。
若对此模拟产生兴趣并希望深入探索，可尝试以下进阶步骤。首先学习“有限长轴承”模型——突破本工具基于的“无限宽”假设，同时考虑轴向压力变化。其控制方程将变为二维偏微分方程：$$ \frac{\partial}{\partial x}\left(\frac{h^3}{\mu}\frac{\partial p}{\partial x}\right) + \frac{\partial}{\partial z}\left(\frac{h^3}{\mu}\frac{\partial p}{\partial z}\right) = 6U \frac{\partial h}{\partial x} $$ 求解需掌握数值分析（如有限差分法）知识，难度显著提升但更贴近实际。
数学基础方面，建议掌握偏微分方程数值解法与无量纲分析思想。为何使用无量纲数（如索末菲数）？其意义在于将物理相似现象归类，整合实验结果与设计数据。例如将模型试验结果缩放至实际设备时，这项技术不可或缺。
具体课题可研究“动压推力轴承”与“静压轴承”。此处涉及的动压效应（楔形效应）可拓展至旋转斜面产生升力的推力轴承；而通过外部供油形成高压油膜支撑的静压轴承，具有发热量小、可实现超精密定位的特点，展现了截然不同的原理与魅力。了解两者将极大拓展轴承设计的视野。