轴承接触应力与寿命预测

在机械传动与旋转机械系统中，轴承是实现支撑、减小摩擦及传递载荷的核心部件。轴承的可靠性直接决定了整个装备的运行稳定性。从固体力学角度来看，轴承失效的核心诱因通常是由于长期循环载荷作用下产生的接触疲劳。因此，深入理解轴承的接触应力分布规律，并

在机械传动与旋转机械系统中，轴承是实现支撑、减小摩擦及传递载荷的核心部件。轴承的可靠性直接决定了整个装备的运行稳定性。从固体力学角度来看，轴承失效的核心诱因通常是由于长期循环载荷作用下产生的接触疲劳。因此，深入理解轴承的接触应力分布规律，并建立准确的寿命预测模型，是进行轴承设计与故障诊断的关键。
轴承内部的接触行为（如滚动体与滚道之间的接触）通常被视为非线性弹性接触问题。在工程实践中，最常用的理论工具是赫兹（Hertz）接触理论。

根据赫兹理论，当两个弹性体（如球形滚动体与圆柱形滚道）相互挤压时，接触区域并非一个点，而是一个具有一定几何形状的接触斑（Contact Patch）。对于球与圆柱的接触，接触形状通常近似为椭圆。

接触应力的核心特征包括：

  • 应力分布：接触压力在接触斑中心达到最大值 $p_{max}$，并向边缘逐渐减小至零。
  • 应力集中：最大接触压力位于接触表面，但最大剪应力（$\tau_{max}$）通常位于接触表面下方一定深度处。这解释了为什么许多轴承失效表现为从内部向表面发展的剥落现象。
  • 影响因素：接触应力的大小取决于载荷 $F$、材料的弹性模量 $E$、泊松比 $\nu$ 以及接触几何参数（如滚动体半径 $R$ 和滚道曲率半径 $R_c$）。

对于球形滚动体，最大接触压力 $p_{max}$ 与载荷 $F$ 的关系大致遵循 $p_{max} \propto F^{2/3}$ 的幂函数关系。这意味着载荷的微小增加会导致接触应力的显著上升。

疲劳失效机制

轴承的失效过程是一个从微观损伤到宏观破坏的过程，主要分为以下几个阶段：

  1. 微观损伤阶段：在循环接触应力的作用下，材料内部产生微观塑性变形或位错堆积。
  2. 裂纹萌生阶段：由于最大剪应力位于表面下方，裂纹往往在材料内部的应力集中点萌生。
  3. 裂纹扩展阶段：随着循环载荷的持续，微裂纹逐渐汇合，形成宏观裂纹。
  4. 剥落阶段（Spalling）：裂纹扩展至表面，导致材料块脱落，接触面积减小，应力进一步集中，最终导致轴承失效。

根据失效位置的不同，可以将疲劳分为表面疲劳（Surface Fatigue）和亚表面疲劳（Subsurface Fatigue）。在高速、轻载条件下，表面疲劳较为常见；而在低速、重载条件下，亚表面疲劳是主要的失效模式。

轴承寿命预测模型

在工程设计中，最常用的寿命评价指标是 $L_{10}$ 寿命，即在 90% 的置信度下，轴承能够运行的额定寿命（以转数计）。

1. 基本寿命计算公式

根据 ISO 281 标准，轴承的寿命计算公式如下：

$$L_{10} = \left( \frac{C}{P} \right)^p \times 10^6 \text{ (revolutions)}$$

其中：

  • $C$：轴承的基本额定动载荷（Basic Dynamic Load Rating），由厂家提供。
  • $P$：轴承在实际运行中的等效动载荷（Equivalent Dynamic Load），需综合考虑径向载荷和轴向载荷。
  • $p$：寿命指数。对于球轴承，$p = 3$；对于滚子轴承，$p = 10/3$。

2. 寿命预测的步骤

在实际工程应用中，进行寿命预测通常遵循以下流程：

  • 第一步：确定载荷条件。计算轴承在工作循环中承受的径向载荷 $F_r$ 和轴向载荷 $F_a$。
  • 第二步：计算等效动载荷 $P$。利用公式 $P = X F_r + Y F_a$，其中 $X$ 和 $Y$ 是取决于轴承类型和内/外圈几何参数的系数。
  • 第三步：应用寿命公式。代入 $C$、$P$ 和 $p$ 值计算 $L_{10}$。
  • 第四步：转换为时间。若需得到运行小时数 $L_{10h}$，则需结合转速 $n$ 进行换算：$L_{10h} = \frac{10^6}{60n} \cdot \left( \frac{C}{P} \right)^p$。

工程计算示例

假设某工程应用中需要选用一款深沟球轴承，已知参数如下：

  • 基本额定动载荷 $C = 25,000 \text{ N}$
  • 实际工作中的等效动载荷 $P = 3,000 \text{ N}$
  • 轴承转速 $n = 1,800 \text{ rpm}$

计算过程：

  1. 计算转数寿命 $L_{10}$：
    由于是球轴承，$p = 3$。
    $$L_{10} = \left( \frac{25,000}{3,000} \right)^3 \times 10^6 \approx 8.33^3 \times 10^6 \approx 578.7 \times 10^6 \text{ revolutions}$$

  2. 计算运行小时数 $L_{10h}$：
    $$L_{10h} = \frac{578.7 \times 10^6}{60 \times 1,800} \approx \frac{578,700,000}{108,000} \approx 5,358 \text{ hours}$$

通过该计算，工程师可以评估该轴承是否满足设计要求的服役时长。

前沿研究方向

随着工业 4.0 和智能制造的发展，轴承的接触应力与寿命预测研究已不再局限于传统的解析公式，而是向以下方向演进：

  • 有限元分析 (FEA)：利用非线性有限元方法，模拟复杂几何形状下的三维应力分布，尤其是考虑材料非线性（塑性）和润滑油膜的影响。
  • 概率寿命预测：考虑到载荷的随机性和材料属性的波动，引入 Weibull 分布等概率模型，实现从确定性寿命到概率性寿命的跨越。
  • 数字孪生与状态监测：结合振动传感器数据与物理模型，通过机器学习算法实时评估轴承的剩余使用寿命 (RUL)，实现从“定期维护”向“预测性维护”的转变。