齿轮箱润滑：仿真与实验对比

齿轮箱中的飞溅润滑看似简单，但若要精确仿真，却绝非易事。齿轮高速旋转，润滑油四处飞溅，每一毫秒需要润滑的表面都在不断变化。自由液面持续碎裂、几何体运动以及流体物性的强敏感性三者叠加，正是传统基于网格的CFD方法力不从心之处，也是仿真方法的选择远比大多数工程师预期更为关键之处。

引言

本文将阐述其中的原因，分析飞溅润滑的物理机制对仿真方法的实际要求，并通过两项shonDy齿轮箱润滑研究的验证结果，展示在方法极限条件下的表现——包括润滑剂在25°C与-40°C之间的行为对比。

飞溅润滑难以仿真的物理根源

在飞溅润滑中，旋转齿轮浸入油池，通过将润滑油甩向轴承、油道和壳体表面来实现润滑。系统中没有专用泵送回路或受控射流，润滑油的分布完全取决于转速、油位、壳体几何形状以及工作温度下润滑剂的黏度。

仿真难题直接源于上述物理过程。齿轮与油池相互作用时，油面持续破碎。液滴形成后撞击壳体壁面，汇聚成油膜，再流回油底壳。因此，润滑油的相界面是一个高度瞬态、拓扑结构复杂的自由液面，在整个仿真过程中不断经历碎裂与重新汇合。与此同时，齿轮几何体处于旋转状态，计算域始终处于动态变化之中。

仿真需要具体预测的是：润滑剂能否到达轴承。这要求精确捕捉自由液面动力学过程——正是这一过程将润滑油输送至宽度仅数毫米的轴承腔入口通道。几何体的持续变化与流体相分布的持续变化相叠加，正是shonDy所采用的粒子法大显身手之处。

基于网格的CFD方法在此场景中的局限

有限体积法（FVM）CFD在固定域或缓慢变形域中的流动仿真已有成熟应用。当几何体发生旋转、自由液面经历大变形（包括碎裂）时，网格必须持续变形或定期重构以保持有效性。两种策略均会带来额外的计算开销，并在界面处引入近似误差。

齿轮箱飞溅润滑的具体难点在于：润滑油界面并非简单的平移或平滑变形，而是在齿轮齿面离开油池时撕裂，在润滑油落回油底壳时重新汇合。基于网格的流体体积（VOF）方法在固定网格上追踪该界面，对于适度的界面变形效果尚可，但随着碎裂程度加剧，近似误差也随之增大。通往轴承腔的狭窄通道（其中的流动由上方振荡油柱驱动）需要精细的局部分辨率，并需仔细处理通道入口附近的界面重构。

旋转几何体进一步加剧了上述难题。标准做法是在旋转区域与静止区域之间使用动网格或滑移界面，这不仅增加了设置复杂度，还在区域边界引入了插值误差。对于具有多个旋转轴及不同转速齿轮副的齿轮箱而言，这将成为一项不容忽视的建模决策。

飞溅润滑同时呈现上述三重挑战：持续碎裂的界面、旋转多体几何结构以及强黏度依赖性。关于粒子法与有限体积法CFD的通用对比，请参阅专题文章；本文其余部分将聚焦于上述组合对齿轮箱润滑的具体影响。

粒子法CFD的不同处理方式

在shonDy这样采用MPS（Moving Particle Semi-implicit）方法的粒子法求解器中，流体由一组携带质量、速度和热力学属性的离散粒子表示。无需对网格进行变形或重构。自由液面由粒子分布本身定义，碎裂、汇聚和成膜过程均可自然处理，无需额外的界面追踪逻辑。

旋转几何体的处理同样直接。在shonDy中，齿轮几何体被定义为运动固壁边界，粒子在每个时间步根据预设的旋转运动与之相互作用。齿轮、轴和壳体均属于同一多体系统，其运动直接施加于仿真，无需滑移网格区域或边界插值。这并非针对运动几何体问题的变通方案，而是无网格方法的天然特性——在粒子-壁面相互作用的计算层面，固定边界与运动边界并无本质区别。

关于多体动力学： 本文所述齿轮箱仿真包含两个齿轮副，齿数比分别为0.4468和0.3607，总传动比为0.1611，各轴以不同转速旋转。两个齿轮副均包含在同一仿真中，无需对计算域进行分区。

作为基础齿轮箱案例的一部分所进行的粒子尺寸对比研究证实，该案例可在单台工作站上以计算可行的分辨率运行。将约500,000个粒子（粒子半径0.42 mm）的基础案例与约1,000,000个粒子（半径0.33 mm）的高分辨率案例进行比较，两种分辨率下流入轴承腔的体积流量表现出相近的规律。这表明500,000个粒子的设置足以捕捉该几何结构中关键的润滑动力学特征。

决定飞溅润滑效果的关键因素

决定飞溅润滑齿轮箱中润滑油能否到达轴承的因素有三：转速、润滑剂黏度和油底壳油量。本文两项研究分别改变转速和黏度，同时保持油量不变，在两种情况下，仿真结果均再现了物理实验录像中观察到的行为。油量是第三个可调参数，但受到约束：较多的油量有助于轴承获得润滑，但同时也会增加搅油损失，而降低这种损失是齿轮箱设计的持续目标。在实物样机制造前，通过仿真找到仍能可靠润滑的最小注油量，正是此类权衡取舍的典型应用场景。

转速决定润滑油能否到达轴承。 低于某一临界转速时，齿轮无法将润滑油甩得足够高以充满壳体顶部的油池（该油池为轴承油道供油），轴承将因缺油而失效，与其他因素无关。在第一项研究中（矿物油，16 cSt，300 ml），480 RPM时油池无法充满，轴承处于干燥状态；980 RPM时油池充满，左侧和中间轴承获得振荡供油，右侧轴承仍无油供给；1453 RPM时分布规律相近，但达到准稳态的速度大幅加快，约为1.8秒，而前者超过5秒。仿真在每个转速下均与实验吻合，包括480 RPM时润滑不足的情况。

黏度决定该转速是否仍然足够。 对热油足够的转速，对冷油可能不够。第二项研究在固定转速1543 RPM下，对同一齿轮箱在两个润滑剂温度下进行了测试。

在25°C（33.6 cSt）时，三个轴承腔均获得润滑供给。在-40°C时，运动黏度约为前者的97倍（3247 cSt），仅左侧轴承获得润滑油，中间和右侧轴承处于干运转状态。相同转速，相同几何结构：仅黏度升高就导致三个轴承中两个失去润滑，仿真精确再现了实验所示的这一变化。

表1：温度研究中的润滑剂物性参数。

结论

飞溅润滑仿真在技术上要求较高，并非因为工程问题晦涩难懂，而是因为碎裂自由液面、旋转几何体与流体物性强敏感性的组合，恰好对基于网格的CFD方法所依赖的核心假设构成压力。本文所述两项齿轮箱案例研究展示了粒子法CFD在此场景中的优势：通过多体动力学框架自然处理运动几何体，无需界面重构即可处理自由液面，并能直接测试在实际温度范围内黏度变化的敏感性。

上述研究证实，该方法能够在多种转速和温度条件下捕捉润滑行为的定性规律，并在大多数情况下实现定量吻合，以直接实验对比作为参考基准。

对于在被动飞溅润滑作为冷却与润滑策略的一部分的齿轮箱或电驱动单元设计工程师而言，仿真方法的选择将直接影响预测结果的可靠性。

如需了解同一方法在完整电驱动单元上的应用——仿真揭示了一个传动轴承完全处于干运转状态——请参阅电驱动单元飞溅润滑仿真。

如果您正在评估用于齿轮箱或电驱动单元润滑的仿真方法，并希望将该方法应用于您自己的设计，欢迎申请试用shonDy。