前段做了个滑动轴承的支撑性能核算轴承不大但转速和载荷条件都比较极限需要把油膜压力分布、空化区域、承载力这些关键量全部算清楚。整套流程用Fluent跑下来最磨人的不是建模反而是空化模型怎么设、油膜力怎么准确提取这两个环节。这篇文章就把我从几何处理到后处理提取力的完整链路理一遍包括一些常规文档里不会写的判断经验和排错记录给做旋转机械润滑分析的朋友做个参考。1. 滑动轴承仿真先想清楚要什么性能指标、计算路线与两个核心难点1.1 滑动轴承的“关键性能”到底指什么滑动轴承的仿真目标说到底就是回答几个工程问题这个转速和载荷下油膜能不能撑住转子、最小油膜厚度够不够安全、摩擦功耗会不会导致温升超标、转子有没有失稳风险。对应到仿真输出核心指标一般有四类承载力也就是油膜对轴颈的总支撑力由压力场在整个轴颈表面的积分决定单位是牛顿最小油膜厚度由给定偏心率直接算出评估边界润滑和磨损风险摩擦功耗与温升壁面剪切应力积分得到摩擦力矩再乘角速度得到功耗刚度与阻尼进阶对动态仿真而言给转子动力学分析提供油膜力对位移和速度的偏导数。很多刚接触轴承CFD的人容易犯一个错误一上来就把网格和模型堆得非常复杂结果算了几个星期提取出来的力却连量级都不对。我个人的习惯是动手之前先把性能指标写清楚再倒推需要哪些物理模型和后处理变量。否则Fluent里可看的结果太多你反而抓不住重点。1.2 为什么空化和力提取是整个流程的两大命门这两个环节之所以关键是因为它们在物理上强耦合、在操作上又容易出错。空化直接影响压力场——油膜发散区一旦出现空化局部压力被限制在空化压力附近不能再往负压区走这会把压力分布的形态从“完整的半正弦波”截掉一块。压力场是力提取的原始输入压力分布错了承载力、偏位角、功耗全部跟着错。而力提取本身也不是Fluent里点一下Report就结束那么简单。它涉及参考压力的设置、面积分的实现、坐标系的转换以及与理论解的交叉验证。有意思的是这两个环节的坑往往是连环的空化压力设得不准力提取环节看起来一切正常但结果和窄轴承近似公式差了一大截你还以为是网格问题。这篇文的顺序就按实际项目推进来排先做前处理和物性再做求解设置最后重点讲力的提取与验证。2. 前处理阶段最容易埋雷的三个地方油膜网格、狭间隙几何与润滑油物性2.1 油膜厚度很小网格该怎么切滑动轴承的特点就是间隙特别小。以我这次核算的某中型轴承为例轴颈直径D50mm轴承宽度L40mm半径间隙c30微米相对间隙c/R≈1.2‰。在偏心率ε0.5的情况下最小油膜厚度只有c(1-ε)15微米最大油膜厚度45微米。15微米对40mm宽的轴承来说尺度比快到3000比1了这决定了网格不能随便用自动划分拍脑袋生成。我的做法是把油膜区域理解成一个“沿着圆周展开的非常薄的楔形通道”用结构化六面体网格去切。径向方向也就是油膜厚度方向至少分12层推荐16到20层因为压力梯度和速度梯度在这个方向上最剧烈。周向方向按角度均匀分段我一般取360到720个控制体。轴向方向根据轴承宽径比调整L/D0.8左右取60到80个节点就够了要关注的是两个端面的出流边界。这里有个容易忽视的细节当偏心率不为零时油膜厚度随周向变化很剧烈最小膜厚处如果网格径向层数不够速度梯度和压力梯度都捕捉不准。我在偏心率超过0.6的工况里会把最小膜厚一侧再加密一到两层保证该区域的首层网格高度不超过最小膜厚的十分之一。网格长宽比虽然难看得吓人但在间隙流动里这是允许的因为物理上压力沿膜厚方向基本不变主要梯度都在径向。2.2 轴颈壁面速度与坐标系对齐前处理的另一个关键问题是边界条件的几何方向。Fluent里处理旋转最简单的方式是把轴颈表面设成Moving Wall并给切向速度ωR不需要真的把整个轴颈固体域建出来。旋转中心要对准轴承的几何中心否则看起来转速对实际速度分解全是歪的。我这次采用的坐标系是这样的Z轴为轴承轴向X-Y平面为径向平面主轴颈绕Z轴旋转。轴颈表面速度方向就是圆周的切线方向大小等于角速度乘以半径。3000rpm对应的角速度是314.16rad/s轴颈半径25mm所以表面速度是7.85m/s。这个数值后续验雷诺数的时候还要用到。还有一个容易被忽略的设定压力参考位置。Fluent默认参考压力是一个标准大气压所有压力都以它为参考点。多相流空化模型里用的饱和蒸汽压通常是绝对压力这两个参考体系一旦搞混空化区域大小会完全失真。我一般把操作压力设成0让整个计算都用绝对压力这样后处理提取力的时候不用反复做加减法。2.3 润滑油物性粘度、密度与“等效空化压力”润滑油参数看起来简单实际上坑也不少。粘度对轴承性能的影响非常敏感我用的某牌号透平油在40摄氏度下动力粘度约0.02Pa·s密度约870kg/m³。如果做常温等温分析这个数据可以接受但高速重载工况下油膜温升可能到十几度甚至几十度粘度随温度下降非常明显这时候等温假设会导致承载力明显高估。物性里面最麻烦的是空化压力。严格意义上的饱和蒸汽压对润滑油来说很低可能只有几十帕。但工程滑动轴承里出现的“空化”很大程度是溶解在油里的空气析出不是油本身的沸腾。空气分离压力比饱和蒸汽压高得多可能到几千帕甚至更高。所以这里要设置的是一个工程等效空化压力而不是查表查润滑油的饱和蒸汽压。实操中我习惯先用一个相对保守的等效空化压力做首轮计算比如绝对压力1000Pa再看空化区域形状是否和理论判断一致空化应该出现在发散楔形区沿着圆周在最大膜厚附近偏下游的位置呈月牙形分布。如果空化区域过大或根本没有再反向调整这个压力值。这个调参过程听起来不“优雅”却是工程仿真里最实在的做法。3. 求解器设置空化模型、层流判断和时间步长的选择逻辑3.1 间隙流是层流还是湍流先算雷诺数再决定很多人一听到高速旋转就上湍流模型其实滑动轴承的间隙流动大部分都在层流范围。雷诺数怎么取对轴承间隙流动常用的定义是Reρ·U·c/μU是轴颈表面线速度c是半径间隙。代入我的工况870×7.85×3e-5/0.02≈102这个数值离层流到湍流的转捩雷诺数一般在1000以上还很远。所以在Fluent里我直接用了层流模型没有开任何湍流模型。这不是偷懒是物理判断。只有当转速特别高、轴承直径大、或者间隙异常大的时候才需要考虑低雷诺数湍流模型比如SST k-omega配合低雷诺数修正。判断依据永远是雷诺数而不是直觉。这一点值得反复强调因为很多人宁可算很贵的湍流模型也不肯花两分钟算一下雷诺数。3.2 多相流与空化模型的参数怎么设空化模拟需要用到多相流。Fluent里滑动轴承空化我推荐用Mixture模型而不是VOF。原因很简单气穴在油膜里以小气泡团簇的形式存在相界面不清晰VOF对这种弥散型两相流既不稳定又算得慢Mixture模型把两相当作混合流体处理配合空化源项在工程尺度上精度和稳定性都够。空化模型我用默认的Zwart-Gerber-Belamri这是工程计算里最常用的一个。主要设置参数包括饱和蒸汽压、气泡半径、成核位点密度和蒸发/冷凝系数。气泡半径默认1e-6米成核位点密度默认5e13蒸发系数默认50冷凝系数0.01。这些默认值对水空化验证过很多但润滑油场景不见得最优。我的经验是优先调整饱和蒸汽压其他参数只有在空化区域明显失真时才动。有个操作小技巧空化源项在压力场剧烈变化时容易让残差猛地抬升。解决办法是把欠松弛因子调低尤其是压力项从默认的0.3降到0.1到0.15先让空化区域稳定下来再慢慢回升。如果你一上来就希望用默认松弛因子拿到收敛解多半会被残差曲线教做人。3.3 边界条件、求解器与时间步的工程取法边界条件其实不复杂。进油口设质量流量入口或压力入口端面出流口设压力出口。关键在压力出口的数值如果参考压力是0端面出口压力应设为油箱或环境压力对应的绝对压力如果计算域包含外部环境就直接用大气压。求解器我习惯用压力基求解器压力速度耦合选Coupled比SIMPLEC更稳代价是内存占用高一点。时间步长的确定原则是每个时间步让轴颈旋转0.5到2度。3000rpm下轴颈转一圈需要0.02秒如果每步旋转1度时间步长就是0.02/360≈5.6e-5秒。这样设置的好处是旋转诱发的压力波前在每个时间步内的移动距离可控空化区演化也能分辨清楚。如果只是求稳态工况下的油膜力可以先用定常计算得到一个初场再切换瞬态继续算这样能节省不少时间。但最终结果一定要基于瞬态稳定后的时均数据因为滑动轴承在高速运转下油膜力本身就有一定波动。4. 油膜力提取的完整操作从Fluent报告到坐标变换与承载力计算4.1 参考压力是力提取的第一大坑先说最坑的地方Fluent的Report→Forces直接输出的力单位是牛顿但它算的是壁面表压的积分而表压是相对参考压力的。如果你操作压力设的是101325Pa那么一个绝对压力等于大气压的面表压是0它对力的积分没有任何贡献。这本身没错但如果你不记得这回事直接用绝对压力场去做面积分误差就是一个大气压乘以整个轴颈表面积算出来会是一个离谱的大数。我的做法是Mesh→Boundary Conditions里把操作压力Operating Pressure改成0让所有压力以绝对压力形式参与计算和后处理。这样Report→Forces读出来的表压就等于绝对压力物理含义清楚也省得后面换算。实际操作中为了让受力分量更可控我给轴颈壁面单独命名了一个边界组比如叫journal-wall后处理时只看这一组。千万不要在整个域上做积分那会把端面出流边界上的压力也算进来。4.2 力的面积分实现与坐标变换虽然Report→Forces很方便但如果你需要做更精细的处理比如分离压力贡献和粘性力贡献或者按特定方向分解力我会把压力场和网格数据导出后用脚本做离线积分。核心公式就是F ∫(p - p_ref) * n * dA其中n是壁面外法线向量dA是面积微元。示意代码大致长这样F_x 0.0 F_y 0.0 for element in wall_elements: A element.area nx, ny element.normal_x, element.normal_y p element.pressure_absolute F_x (p - p_ref) * A * nx F_y (p - p_ref) * A * ny F_total (F_x**2 F_y**2) ** 0.5这里的p_ref就是操作压力如果设成了0那括号里就是p本身。把Fluent导出的壁面压力数据按这个逻辑处理一遍和软件自带报告交叉验证两者应该在数值上非常接近。粘性剪切力对承载力的贡献通常可以忽略但对摩擦功耗很重要积分公式里换成粘性剪应力τ再乘力臂就能得到扭矩。4.3 数值实例从分量合成到偏位角判断假设我这次计算稳定后得到轴颈壁面的力分量是F_x-523NF_y6480N。那么总承载力就是F sqrt(523² 6480²) ≈ 6500N如果外载荷是竖直向下的6500N那么油膜支撑力应该竖直向上Y分量为正就是对的。横向的523N说明轴心位置和外载荷方向之间存在一个小偏位角这个角度可以估算为α atan(523 / 6480) ≈ 4.6度在滑动轴承理论里这个角对应载荷线与偏心线之间的夹角称为偏位角。它的合理范围一般随着偏心率变化偏心率越大偏位角越小。如果你算出的偏位角和经验值差得离谱比如十几二十度那多半不是物理现象而是网格周向划分不对称、坐标轴没对齐或者空化压力设错了导致压力分布畸变。4.4 一个方便交叉验证的量纲估算在Fluent结果出来之前最好先用窄轴承近似公式粗算一遍承载力量级。我的工况数据D50mmL0.04mc3e-5mε0.5ω314.16rad/sμ0.02Pa·s。按无限短轴承近似估算特征力级在μ * ω * R * L³ / c² 0.02 * 314.16 * 0.025 * (0.04³) / (3e-5²) ≈ 1.1e4 N再乘一个随偏心率变化的系数ε0.5时大致在0.6量级粗估在6到7kN的区间。这和我的Fluent结果6500N能对上说明模型设置和边界条件方向正确。如果CFD结果比这个量级高出一倍或低一半先回头查空化、参考压力、端面出流边界基本都能找到问题。5. 收敛性排查与结果可信度几张表讲清楚我踩过的坑5.1 发散和不收敛的排查清单滑动轴承CFD最常见的发散原因不是物理模型不对而是网格和前处理细节。我把排查问题按出现频率整理成一张表方便对照。现象常见原因处理方式残差震荡不下降压力欠松弛系数过大压力松弛降到0.1~0.15空化区域过大等效空化压力设得太高降低到合理范围没有空化区域操作压力参考体系错误检查绝对压力与饱和蒸汽压一致性Fx/Fy方向混乱坐标系或转速方向设置错误逐一核对旋转轴方向计算中途发散网格出现负体积检查偏心率大时最小膜厚处网格质量力报告值跳变瞬态未稳定气穴周期性脱落取稳定段时均值在这些坑里我花时间最多的是“没有空化区域”。当时把润滑油的饱和蒸汽压填成几十帕理论上没问题但实际油中溶解空气析出压力远高于这个值导致空化一直不出现。后来把等效空化压力设成1000Pa绝对压力空化区域才正常出现压力场的负压段也被正确地削平了。5.2 残差、质量守恒与网格无关性检查残差收敛只是最低要求更重要的是物理量守恒。每轮计算结束后我都会单独检查端面泄漏流量和进油流量的平衡泄漏量由进油口流量和空化区引起的质量源共同决定但总体质量流率偏差应该控制在0.1%以内。如果质量守恒差很多回到网格质量排查。网格无关性验证我做了三套基准网格径向16层周向480加密网格径向20层周向720粗网格径向12层周向360。结果对比如下网格方案径向层数周向控制体承载力(N)与基准偏差粗网格123606210-4.5%基准网格1648065000加密网格2072065300.5%粗网格和周向加密网格偏差约4%可以接受但偏心率如果升到0.7以上粗网格偏差会迅速扩大。所以建议偏心率高的工况至少用基准网格起步再视算力决定是否加密。5.3 瞬态力波动与时均值取法高速轴承的油膜力并不是一条平稳直线。空化气泡的产生和溃灭、供油压力的扰动都会让Fx和Fy出现周期性波动。从Report→Forces里看到力在几百牛顿范围内来回跳是很正常的不要急着怀疑模型。我的取值方法是瞬态计算稳定后记录至少五个旋转周期内的力时间历程然后对时间做平均。如果波动的幅值在平均值的10%以内用平均值做承载力判断完全够用如果波动大到15%到20%就要怀疑是否存在气穴周期脱落导致的压力脉动这时除了取平均还要把这个脉动幅值报给下游转子动力学分析因为它可能是一个潜在的激励源。对摩擦功耗也一样处理。积分得到摩擦力矩后功耗PT·ω。我的算例里摩擦力矩在0.35N·m左右乘以角速度314.16后功耗约110W。这个数看起来不大但如果轴承润滑冷却不足全部变成油膜温升效果就不可忽视了。6. 从稳态分析到更复杂的工况留几个扩展方向6.1 热效应与弹性变形的叠加我做的是常温等温分析但工程上很多轴承工况不允许忽略热效应。油膜温升会改变粘度粘度分布又反过来改变压力场和剪切应力场这是双向耦合。Fluent里可以用能量方程加粘性耗散项来做流热耦合润滑油粘度按温度的经验关系式设置。代价是计算量明显增加收敛难度也上一个台阶。弹性变形是另一个方向。高比压轴承的轴瓦表面在油膜压力下会发生弹性变形反过来改变油膜间隙形状。这种情况建议走双向流固耦合或者至少做一次压力场→变形→几何更新的迭代。如果只是粗略评估可以先忽略弹性变形但要把这个假设写进报告结论里提醒决策方注意边界。6.2 从给定偏心率到给定载荷的平衡位置计算最后想说一个容易被误解的点本文这套流程的前提是给定轴心位置也就是给定偏心率和偏位角求油膜力。但在真实工程里往往是给定外载荷问轴心会浮在哪里。这个逆向问题需要动网格或6自由度模型让轴心在外载荷和油膜力作用下逐渐移动到平衡位置。Fluent的动网格加6DOF求解器可以做到但计算成本和时间成本都不是稳态分析能比的。我的建议是先跑通本文的稳态正向流程拿到力和压力分布再往动态平衡走。否则一上来就做全耦合的瞬态动网格分析出问题了你都不知道是空化模型的事还是刚体运动求解的锅。按这个顺序走一遍滑动轴承的CFD分析基本上就不会出现“算了很久结果却没法用”的尴尬情况了。哪怕只是想把Fluent当辅助校核工具从空化到力的提取这条主线理顺了后面的扩展才有底子。
