在阅读此文前，为了方便您进行讨论和分享，麻烦您点击一下“关注”，可以给您带来不一样的参与感，感谢您的支持。

文 |夙烨

编辑 | 夙烨

前言

滚动轴承作为旋转机械设备的关键部件，其运转状况直接关系着整个设备的安全运行口，对轴承进行实时剩余寿命预测尤为重要。

轴承故障诊断领域中一个重要研究方向，就是通过设备状态监测数据对其性能特征进行分析并进行剩余寿命预测，主要包括两个方面，一是对轴承性能退化特征指标进行预测，另一方面是对轴承剩余工作寿命进行预测。

而对滚动轴承剩余寿命的影响因素有很多：一方面滚动轴承在实际工作中，由于工作环境复杂和运行条件的多变性，其疲劳损伤的发展趋势无固定方向且故障种类多样。

这也造成了滚动轴承的剩余寿命离散性较大且具有较强的随机性。

另一方面轴承工作状态监测数据受工作环境和运行状态等条件影响比较大，也存在诸多随机性，因此通过采用上述方法不能反映滚动轴承的真实工作状态，故很难实现对轴承剩余寿命的精准预测。

基于威布尔分布的轴承比例故障模型对轴承振动数据分布，以及其残差分布无固定要求，且通过威布尔分布对轴承故障数据进行拟合能够实现对轴承故障早期数据进行识别。

此外，该模型能够实现对多个故障特征轴承的剩余寿命计算，因此该模型被广泛应用于轴承寿命预测研究中，针对上述问题，可以提出了一种基于PCA和威布尔分布的轴承剩余寿命预测方法，以实现滚动轴承实时剩余寿命的精准预测。

PCA降维方法

作为一种最常见的降维方法，主成分分析法，将存在一定关系的数据向量通过正交变换降维处理后转化为线性独立的变量。

它的原理是通过投影的方式，将高维数据映射到低维空间来实现数据的降维，其目的是在确保信息丢失降至最低的条件下，用低维数据向量来表示高维数据变量。

假设一个样本数据有n个特征参数，所有的特征参数能够完全表示该样本的信息，但是这n个特征参数存在一定的信息重叠，经过PCA处理后的数据，在保留原数据主要信息的同时删除了特征参数的信息重叠。

因此可以采用PCA对轴承振动数据进行降维处理。

PCA计算方法其实就是在n个特征向量中，定义一个新的空间，将原特征向量进行映射后实现特征集的降维，从n维降到k维。

当计算得到的累计贡献率大于95%时，则表明选取的个特征能够表示原数据95%以上的信息，即特征值选取实现了特征集的降维。

轴承性能退化特征参数的PCA降维

为构建轴承的性能退化指标，我们采用IEEEPHM2012数据挑战赛提供的轴承剩余寿命，预测的数据集该挑战赛将轴承进行分组后，分别对每个轴承进行加载运行试验。

试验停止后记录每个轴承的故障状态，试验结果为17个轴承中10个轴承出现障特征。

为了验证PCA是否能实现对轴承性能退化特征的降维，将试验的17个轴承分组，其中6个轴承为训练对象，11个轴承为测试对象，计算训练轴承的各个特征参数，构建轴承性能退化特征参数集并对构建得到的特征参数集进行PCA降维处理。

随后对6个训练轴承去噪后的信号进行分析，将选取信号的特征参数组成特征参数集，其中选取的特征参数分别是时域特征参数、频域特征参数以及时频特征参数，并且是方均根值、峭度系数、峰值、频谱均值、频谱方均根值和EEMD归一化。

此处选取第一组轴承中的1号轴承为例进行分析，其特征参数变化如下图所示：

从下图中可以看出，轴承在开始工作后的前3h内，6个特征参数均未产生明显的波动，这说明在此时间内轴承在平稳运行，平稳运行期过后，各特征参数均开始产生较大幅度的震荡，说明轴承性能正在进行退化，且越来越剧烈。

但是在轴承性能退化过程中，各参数并不能清晰地描绘出轴承具体退化特征，即不能区分轴承退化阶段，并且各个特征参数之间存在着信息混叠。

因此，可将轴承的6个特征参数构建成性能退化特征参数集，然后对其进行PCA处理后得到轴承性能退化指标。

需要特别注意的是，由于各个轴承安装工艺以及试验环境的不同，会导致不同轴承的相同性能特征参数存在差异。

因不同轴承在退化同一时间段内其方均根值不尽相同，故在对性能退化参数集进行PCA处理之前，须先对计算得到的特征参数进行标准化。

首先计算所有采样时间点的不同特征参数变量的平均值，然后用各个采样点的特征参数与计算得到的参数变量平均值相除，得到的结果就是该特征参数的标准化数列，计算公式如下：

随后对训练样本轴承的各特征参数变量实现标准化后，选取每个轴承的试验采样点，试验点的选取根据轴承退化阶段，分别选取20个退化试验初始时的采样点、50个退化中期的采样点和130个退化后期采样点，共200个采样点。

6个训练轴承每个轴承6个特征参数共构成12006的特征参数集，采用PCA方法对构建的轴承特征参数集进行降维处理，从PCA降维处理后的结果得到特征参数集的前两个主元，累计贡献率达到了原始数据的98.23%。

计算结果表明：轴承原始特征参数集98.23%的信息可以通过前两个主元表示出，即此两个主元可以反映出轴承的整个性能退化过程。

轴承性能退化特征集的前两个主元输出如下图所示，从图中可以看出此两个主元可以清晰地将轴承退化阶段区分开，即可以反映出轴承的性能退化过程。

轴承威布尔分布故障模型

轴承在使用寿命期间，其发生故障的概率与工作时间存在一定的概率分布关系，不同轴承的寿命以及全寿命期间出现故障的次数不同，但无论处在什么阶段轴承的剩余寿命与其故障发生概率的关系都是一致的。

轴承故障率与工作时间之间的关系曲线称为“浴盆曲线”轴承在t时刻的故障率计算公式如下：

式中：F(t)为累计故障分布函数；f(t)为故障密度分布函数；R(t)为可靠度函数。

威布尔分布利用概率值可以很容易计算出它的分布参数，通过参数的调整，就可以实现产品全寿命周期的表示。

即浴盆曲线的早期失效、随机失效和老化失效的3个阶段，因此它被广泛应用于各种寿命估算工程中，轴承在威布尔分布中t时刻的故障密度分布函数为：

则其故障分布函数为：

可靠度函数为：

代入到式中计算故障率函数为：

平均寿命函数为：

式中：η为尺度参数；m为形状参数，表示函数走势，y为位置参数，Г函数可以通过查表求得，因轴承在安装后就伴随着故障的发生，其位置参数y为0，故轴承威布尔分布比例故障率模型为：

该模型的可靠度函数为：

剩余寿命函数为：

从比例故障模型的可靠度函数和剩余寿命函数公式来看，在计算之前需先确定尺度参数η、形状参数m和位置参数γ，极大似然估计法常用在模型参数估算，采用此方法对故障模型参数进行估计具体步骤为：

构造比例故障模型的似然函数为：

对似然函数取对数分别对η、m和γ求偏导，并令：

此处采用牛顿迭代法对方程进行求解，求得参数似然估计值计算结果如表1所示：

轴承性能退化指标的灰色预测

在对轴承进行寿命预测之前，如何对轴承性能退化指标进行预测是研究的关键，此处选用灰色模型GM(1，1)来对性能退化指标进行预测。

以第3个训练轴承为例，对其去噪处理后的振动信号进行特征集提取后做PCA降维，得到该轴承性能退化指标变化如下图所示：

从图中可以看出，轴承开始运行的前4h性能退化指标变化较平稳，第5h后开始剧烈波动，即轴承开始退化加剧。

而选取轴承开始工作第5h后的连续9个采样节点作为灰色模型的原始序列，并通过建立好的灰色模型，对第10个时间节点的性能特征集前两个主元进行预测，预测结果如下图所示：

我们分别采用参差大小检验法和后验差检验法，对所建立的灰色模型精度进行检验，通过预测结果对该模型的平均相对误差Δ均方差比C和小概率误差p进行计算。

计算结果为Δ=0.94，C=0.33，p=0.91，各指标均符合模型精度要求，故该灰色模型可以实现对轴承性能退化指标的精确预测。

剩余寿命预测

提取第3个测试轴承开始工作5h后连续2h的轴承退化特征参数并构建灰色模型，通过建立好的灰色模型，分别对轴承性能退化特征集的第一主元和第二主元进行灰色预测。

并将预测结果代入威布尔比例故障函数中，对测试轴承进行剩余寿命计算，计算结果如下图所示：

计算得到轴承剩余寿命各时间节点预测误差如表2所示：

从上图和表2中可以看出，随着轴承工作时间的推进，轴承剩余寿命预测结果越来越精准，分别对其他测试轴承进行剩余寿命计算，并与轴承真实寿命做对比，结果如表3所示：

轴承工作初期剩余寿命预测误差值较大，是因为轴承安装后运行初期其状态不稳定导致的故障偶发，以及灰色模型的线性变化特点，随着工作时间的推进，该模型的预测精度越来越高。

综上所述，采用基于PCA和威布尔分布的剩余寿命预测方法，可以有效地对轴承剩余寿命进行预测。

结论

针对滚动轴承的剩余寿命随工作时间的推进越来越难预测的问题，通过对滚动轴承的故障特征参数集的分析，提出了一种基于PCA和威布尔分布的滚动轴承剩余寿命预测方法。

首先对轴承的特征参数进行计算，构建故障特征参数集，经过PCA降维处理后得到故障特征集的前两个主元即可清晰反映出轴承的性能退化过程。

然后构建轴承威布尔分布比例故障模型，对轴承性能退化指标进行灰色模型预测后将预测结果代入比例故障模型中，并通过对故障模型进行计算得到轴承的剩余寿命预测值。

最后对比测试轴承的实际寿命和预测寿命，对比结果证明，本研究方法可以实现对轴承剩余寿命的精准预测。

该方法虽然在一定程度上可以实现对轴承剩余寿命的准确预测，具有一定的应用价值，但是研究方法仅仅局限于轴承振动信号特征的分析，未能结合轴承油液分析和温度监测等轴承故障监测方法。

因此，下一步研究方向可以将振动信息参数与油液分析以及温度监测等参数进行融合，并通过改进PCA降维方法对复杂的多参数轴承故障特征进行降维处理，以实现对轴承剩余寿命的更加精准的预测。

基于雨流模型构建的循环数统计原理，提出了一种道路损伤的计算方法，实现了不同道路损伤值的对比分析，为车辆耐久及可靠性试验不同路谱下损伤分析提供了一种新的方法。

所提供的模型分析方法同样可应用到其他耐久及疲劳强度分析中，具有广泛的适用性及实用价值。

依据上述结论，制定了针对该车型的试验场测试规范，其总里程为5000km，实际测试时间为目标规范的60%。

试验结果表明：新规范下的汽车测试能够有效缩短项目开发周期，提高可靠性试验准确性，增强产品竞争力。

参考文献

1.韩清凯，云向河，李宁，等． 大型滚动轴承故障诊断及寿命评估技术进展［J］． 轴承，2021( 9) : 1-13．

2.李苗苗，杨阳，王宇，等． 考虑温度影响的角接触球轴承疲劳寿命 仿 真 研 究［J］． 机 械 制 造 与 自 动 化，2020，49 ( 2 ) :113-116．

3.赵嘉芳，杜杰，等． 基于振动信号的轴承寿命预测［J］． 国外铁道车辆，2021，58( 1) : 35-40．

4.郑小霞，蒋海生，刘静，等． 基于变分模态分解与灰狼算法优化极限学习机的滚动轴承故障诊断［J］． 轴 承，2021 ( 9) :48-53．