技术领域
[0001] 本
发明涉及振动故障诊断领域，特别是一种用于滚动
轴承内圈磨损故障诊断及趋势预测的方法。
背景技术
[0002] 随着人们生活
水平的提高，对于物质需求的
质量和数量都在不断提高，工业生产需要有较高的效率，一旦用于工业生产的机械发生故障，将会极大的影响工业企业的生产进度，带来巨大的经济损失甚至危及人身安全。
滚动轴承是各种机械装备中经常使用的零件。滚动轴承的长期工作容易使得滚动轴承的内圈极易磨损。当滚动轴承发生轻微磨损时，其仍可以正常工作，当滚动轴承内圈磨损达到一定程度时，滚动轴承便无法正常工作，此时如果继续工作将会违背安全生产原则，并且会影响产品合格率。因此，及时的对滚动轴承内圈磨损故障做出诊断及趋势预测是十分必要的。这样可以及时了解滚动轴承的状态，使得滚动轴承发挥最大的作用。目前应用于滚动轴承故障诊断的方法有小波包分析(Wavelet packet analysis)方法、
深度学习(Deep learning)方法和
频谱分析(Spectrum analysis)方法等，这些方法存在较大的局限性。小波包分析方法依赖于小波基的选择，深度学习方法对故障的诊断具有很好的效果，但是其具有较大的计算量，实际应用较困难，频谱分析方法可以有效地发现滚动轴承的故障
频率，但是对于故障程度无法做出具体判断，因此目前没有有效的方法可以同时实现滚动轴承内圈磨损故障的诊断和趋势预测。基于计算时域指标，根据时域指标值大小进行故障诊断及故障趋势预测在理论上是可以实现的。通过计算时域指标值找出故障趋势的规律，根据规律判断滚动轴承内圈磨损故障的程度，预测出设备的停机时间，使得在不影响产品质量的前提下，滚动轴承能够发挥最大作用。
发明内容
[0003] 本发明的目的是为了解决上述问题，设计了一种用于滚动轴承内圈磨损故障诊断及趋势预测的方法。
[0004] 实现上述目的本发明的技术方案为，一种用于滚动轴承内圈磨损故障诊断及趋势预测的方法，包括以下步骤：
[0005] 步骤1：基于动态时间规整(DTW)
算法的振动
信号预处理；
[0006] 步骤2：时域指标的计算；
[0007] 步骤3：滚动轴承内圈磨损故障的诊断及趋势预测。
[0008] 进一步地，所述的步骤1中，基于动态时间规整(DTW)算法的振动信号预处理，其具体步骤如下：
[0009] 步骤1.1：将参考信号A与测试信号B同时输入到DTW算法；
[0010] 步骤1.2：通过DTW算法获得最优规整路径WR；
[0011] 步骤1.3：使用重
采样技术将原始残差向量信号OR恢复原长。
[0012] 进一步地，所述的步骤2中，时域指标的计算，其具体步骤如下：
[0013] 步骤2.1：计算传统时域指标值；
[0014] 步骤2.2：计算时域新指标。
[0015] 进一步地，所述的步骤3中，滚动轴承内圈磨损故障的诊断及趋势预测，其具体步骤如下：
[0016] 步骤3.1：做出时域新指标的曲线图；
[0017] 步骤3.2：根据时域新指标曲线图判断滚动轴承内圈磨损故障的程度及其故障趋势。
[0018] 进一步地，所述的步骤1.1中，参考信号A与测试信号B同时输入到DTW算法，其具体为，滚动轴承的振动信号是由
加速度振动检测
传感器测得的，滚动轴承在无故障情况下测得的信号为健康振动信号Xh(t)，滚动轴承在内圈磨损故障情况下测得的信号为故障振动信号Xf(t)，健康振动信号作为DTW算法的参考信号A，故障振动信号作为DTW算法的测试信号B，参考信号A与测试信号B同时输入到DTW算法。
[0019] 进一步地，所述的步骤1.2中，通过DTW算法获得最优规整路径WR，其具体为，通过动态时间规整算法计算参考信号A的数据点与测试信号B的数据点的最小累积距离，获得最优规整路径WR。
[0020] 进一步地，所述的步骤2.1中，计算传统时域指标值，其具体为，健康振动信号Xh(t)和故障振动信号Xf(t)分别需计算的传统时域指标的数量为4个，分别是峰值Peak、均方根值RMS、
波峰因数CF、峭度值Ku。
[0021] 进一步地，所述的步骤2.2中，计算时域新指标，其具体为，在计算信号的4个传统时域指标值之后，需要计算时域新指标TALAP和TRCKT，其具体计算方法如下：
[0022]
[0023]
[0024] 其中，XRMS表示信号的均方根值，RMS0表示健康信号的均方根值，log表示取对数。
[0025] 进一步地，所述的步骤3.1中，做出时域新指标的曲线图，其具体步骤如下：
[0026] 步骤3.1.1：根据滚动轴承健康情况下和内圈不同磨损程度情况下的四个时域指标值，计算健康振动信号Xh(t)和故障振动信号Xf(t)的TALAP值和TRCKT值；
[0027] 步骤3.1.2：拟合出TALAP和TRCKT的曲线图。
[0028] 进一步地，所述的步骤3.2中，根据时域新指标曲线图判断滚动轴承内圈磨损故障的程度及其故障趋势，其具体为，根据TALAP的曲线的变化趋势，TALAP曲线图分为三个区域，早期微弱磨损区域、中度磨损区域、严重故障预警区域与失效区域的
混叠区域，根据TRCKT的曲线变化趋势，TRCKT曲线图分为两个阶段，非完全下降区域与下降区域，下降区域对应TALAP曲线图的失效区域，当滚动轴承的振动信号的TALAP值处在早期微弱故障磨损区域和中度磨损区域时，滚动轴承仍可以继续正常工作，当达到严重故障预警区域与失效区域的混叠区域时，要根据TRCKT曲线图判断TRCKT值是否处于下降区域，若处于下降区域则应立即停止运行。
[0029] 利用本发明的技术方案制作的一种用于滚动轴承内圈磨损故障诊断及趋势预测的方法，通过该方法可以对滚动轴承的振动信号进行预处理，得到包含大量故障信息的振动信号，使之适用于时域指标计算，便于进行故障诊断和趋势预测，TALAP和TRCKT指标综合了传统时域指标的特点和优势，避免了多指标综合判断的复杂性，仅可根据TALAP指标实现故障诊断，根据TRCKT指标实现故障趋势预测，降低了故障诊断和趋势预测的难度，方法简单实用，此方法适用于任何恒定转速条件下的滚动轴承内圈磨损故障的诊断及趋势预测。因此，该发明具有实用价值，能够为维修人员提供准确的滚动轴承运行信息，减少因轴承故障带来的经济损失，满足了故障诊断及趋势预测的需要。
附图说明
[0030] 图1是本发明所述一种用于滚动轴承内圈磨损故障诊断及趋势预测的方法的
流程图。
[0031] 图2是滚动轴承转速1728r/min情况下不同磨损程度的原始残差向量信号时域图。
[0032] 图3是滚动轴承转速1750r/min情况下不同磨损程度的原始残差向量信号时域图。
[0033] 图4是滚动轴承转速1772r/min情况下不同磨损程度的原始残差向量信号时域图。
[0034] 图5是滚动轴承转速1796r/min情况下不同磨损程度的原始残差向量信号时域图。
[0035] 图6是滚动轴承转速1728r/min情况下不同磨损程度的残差向量信号时域图。
[0036] 图7是滚动轴承转速1750r/min情况下不同磨损程度的残差向量信号时域图。
[0037] 图8是滚动轴承转速1772r/min情况下不同磨损程度的残差向量信号时域图。
[0038] 图9是滚动轴承转速1796r/min情况下不同磨损程度的残差向量信号时域图。
[0039] 图10是时域新指标TALAP在不同转速条件下的曲线图。
[0040] 图11是时域新指标TRCKT在不同转速条件下的曲线图。
具体实施方式
[0041] 下面结合附图对本发明进行具体描述，如图1所示，一种用于滚动轴承内圈磨损故障诊断及趋势预测的方法，包括以下步骤：
[0042] 步骤1：基于动态时间规整(DTW)算法的振动信号预处理；
[0043] 步骤2：时域指标的计算；
[0044] 步骤3：滚动轴承内圈磨损故障的诊断及趋势预测。
[0045] 下面通过一个具体
实施例对本发明进行具体阐述。
[0046] 在本实施例中，数据来源于美国凯斯西储大学轴承
数据中心，所使用数据的基本参数如表1所示，实验轴承的基本参数如表2所示，参考信号的基本参数如表3所示。
[0047] 表1 实验所使用数据的基本参数
[0048]
[0049] 表2 轴承基本参数表
[0050]
[0051] 表3 参考信号基本参数
[0052]
[0053] 下面应用动态时间规整算法进行预处理，从原始数据中截取实验所使用的所有信号，长度为2400点。
[0054] 步骤1：基于动态时间规整(DTW)算法的振动信号预处理；
[0055] 步骤1.1：滚动轴承的振动信号是由加速度振动检测传感器测得的，健康振动信号作为DTW算法的参考信号A，相应的故障振动信号作为DTW算法的测试信号B，参考信号A与测试信号B同时输入到DTW算法。
[0056] 步骤1.2：通过DTW算法获得最优规整路径WR；通过DTW算法计算参考信号A的数据点与测试信号B的数据点的最小累积距离，获得最优规整路径WR；DTW算法的具体计算过程如下：
[0057] (1)输入参考信号A和测试信号B；
[0058] (2)DTW算法处理的过程就是将参考信号A的数据点与测试信号B的数据点进行匹配，顺序的以参考信号A中的每一个数据点作为计算欧氏距离的标准数据点，寻找测试信号B中与之欧式距离最小的数据点，该最小数据点便于标准数据点匹配成功，对应的两点形成规整路径的一个坐标，所有坐标按照顺序排列形成完整的规整路径，匹配点之间欧氏距离最小时，累计距离便是最小值。
[0059] 在规整路径中，第i对匹配数据点表示为：wi＝(Am，Bn)，第i对匹配数据点的规整路径坐标表示为：wRi＝(m，n)，匹配结果表示为：W＝w1，w2，w3，……，wi，……，wj，完整的规整路径表示为WR＝wR1，wR2，wR3，……，wRi，……，wRK，
[0060] 累计距离可以表达为如下形式：
[0061]
[0062] 其中，m是来自参考信号A的索引，n是来自测试信号B的索引，K为匹配数据点对的总数。因此动态时间规整必须满足以下3个条件：
[0063] 条件1(边界条件)：必须满足wRi＝(1，1)，wRK＝(K，K)；
[0064] 条件2(连续性条件)：给定wRi＝(m，n)，且wR(i-1)＝(m′，n′)，则m-m′≤1，且n-n′≤1；
[0065] 条件3(单调性条件)：给定wRi＝(m，n)，且wR(i-1)＝(m′，n′)，则m-m′≥0，且n-n′≥0。
[0066] (3)输出规整路径WR；
[0067] 得到规整路径后，可以得到规整后的信号，将参考信号A与测试信号B直接做差得到原始残差向量信号OR，图2所示为转速1728r/min情况下的原始残差向量信号时域图，图3为转速1750r/min情况下的原始残差向量信号时域图，图4为转速1772r/min情况下的原始残差向量信号时域图，图5为转速1796r/min情况下的原始残差向量信号时域图。
[0068] 步骤1.3：使用重采样技术将原始残差向量信号OR恢复原长；重采样技术的具体计算过程如下：
[0069] (1)输入原始残差向量信号OR＝OR1，OR2，……，ORi，……，ORj和规整路径WR，令wRi＝(a，b)，wR(i+1)＝(a′，b′)，将a设置为参考信号A的索引，b设置为测试信号B的索引，m为从1到K的整数，i和m的初始值均为1；
[0070] (2)判断b′是否与b相等。若二者相等，则令Rm取ORi和ORi+1中的较大值；若二者不相等，则令Rm＝ORi；
[0071] (3)令i＝i+1；
[0072] (4)判断i是否等于j，若不相等，令m＝m+1，返回(2)继续循环，若相等则停止循环；
[0073] (5)输出残差向量信号R＝R1，R2，……，Rm，……，RK；
[0074] 在得到残差向量信号后，经过重采样算法的处理，信号恢复原长，图6所示为转速1728r/min情况下的残差向量信号时域图，图7为转速1750r/min情况下的残差向量信号时域图，图8为转速1772r/min情况下的残差向量信号时域图，图9为转速1796r/min情况下的残差向量信号时域图。
[0075] 步骤2：时域指标的计算；
[0076] 步骤2.1：健康振动信号Xh(t)和故障振动信号Xf(t)分别需计算的传统时域指标的数量为4个，分别是峰值Peak、均方根值RMS、波峰因数CF、峭度值Ku；传统时域指标的计算公式如下所示：
[0077]
[0078]
[0079] Peak＝Xmax
[0080]
[0081] XRMS表示信号的均方根值，Xmax表示信号的最大值，xk表示信号中的第k个数据点，表示信号中各个数据点的平均值。不同转速和磨损程度条件下传统时域指标的计算结果如表4所示。
[0082] 表4 不同转速和磨损程度条件下传统时域指标的计算结果
[0083]
[0084] 步骤2.2：在计算信号的4个传统时域指标值之后，需要计算时域新指标TALAP和TRCKT，其具体计算方法如下：
[0085]
[0086]
[0087] 其中，RMS0表示健康信号的均方根值，log表示取对数。不同转速和不同磨损程度条件下TALAP和TRCKT的计算结果如表5和表6所示。
[0088] 表5 不同转速和不同磨损程度条件下TALAP的计算结果
[0089]
[0090] 表6 不同转速和不同磨损程度条件下TRCKT的计算结果
[0091]
[0092] 步骤3：滚动轴承内圈磨损故障的诊断及趋势预测；
[0093] 步骤3.1：做出时域新指标的曲线图；
[0094] 步骤3.1.1：根据滚动轴承健康情况下和内圈不同磨损程度情况下的四个时域指标值，计算健康振动信号Xh(t)和故障振动信号Xf(t)的TALAP值和TRCKT值；
[0095] 步骤3.1.2：拟合出TALAP和TRCKT的曲线图。不同故障程度与不同转速情况下的TALAP曲线图如图10所示，不同故障程度与不同转速情况下TRCKT的曲线图如图11所示。
[0096] 步骤3.2：根据时域新指标曲线图判断滚动轴承内圈磨损故障的程度及其故障趋势。根据TALAP的曲线的变化趋势，TALAP曲线图分为三个区域，早期微弱磨损区域、中度磨损区域、严重故障预警区域与失效区域的混叠区域，根据TRCKT的曲线变化趋势，TRCKT曲线图分为两个阶段，非完全下降区域与下降区域，下降区域对应TALAP曲线图的失效区域，当滚动轴承的振动信号的TALAP值处在早期微弱故障磨损区域和中度磨损区域时，滚动轴承仍可以继续正常工作，当达到严重故障预警区域与失效区域的混叠区域时，要根据TRCKT曲线图判断TRCKT值是否处于下降区域，若处于下降区域则应立即停止运行。
[0097] 综上，一种用于滚动轴承内圈磨损故障诊断及趋势预测的方法能够有效地对滚动轴承内圈磨损故障做出准确的故障诊断和故障趋势预测，并可以确定该种故障所处的恶化阶段。DTW算法与时域指标结合使用，能够充分提取滚动轴承振动信号中的故障信息，为设备维护技术人员提供了一种新的滚动轴承状态监测方法。
[0098] 上述技术方案仅体现了本发明技术方案的优选技术方案，
本技术领域的技术人员对其中某些部分所可能做出的一些变动均体现了本发明的原理，属于本发明的保护范围之内。
