基于集合固有时间尺度分解算法的滚动轴承故障诊断方法
【技术领域】
[0001] 本发明涉及一种滚动轴承故障诊断方法。特别是涉及一种基于集合固有时间尺度 分解算法的滚动轴承故障诊断方法。
【背景技术】
[0002] 振动分析是对滚动轴承进行故障诊断最简单最直接的手段，典型的振动分析方法 包括：小波变换、魏格纳分布、经验模态分解等，但它们都有各自的缺点。小波变换不具有 自适应性且基函数的选择过于依赖使用者的经验。魏格纳分布具有很高的时频分辨率，但 是交叉项的出现限制了它的应用。经验模态分解是一种自适应的时频分解方法，在动力机 械故障诊断中得到了较广泛的应用，但该方法存在着过包络、模态混淆、端点效应以及由 Hilbert变换带来的无法解释的负频率等问题。固有时间尺度分解算法是一种新的非平稳 信号分析方法，它弥补了经验模态分解的很多不足。但是它仍然存在着比较严重的模态混 叠问题。
[0003] 为了解决经验模态分解的模态混叠问题，wu等人提出了集合经验模态分解算法， 该方法通过向信号多次添加相同幅值的白噪声来引导信号向对应的尺度分解，最后通过集 总平均消除加入的噪声。集合经验模态分解算法很好的解决了经验模态分解的模态混叠问 题，但它同时带来了其它两个问题。第一个问题是添加的白噪声不能完全抵消，虽然可以 通过选择较大的添加噪声次数来降低最终的白噪声残余，但是相应的计算时间也会大大增 加，这有悖于故障诊断对算法速度的要求。另外一个问题是不同的白噪声的加入可能会导 致每次分解产生的旋转分量个数不相等，从而给最后的集总平均带来困难。简单的套用集 总经验模态分解的框架来解决固有时间尺度分解算法的模态混叠同样会带来上述两个问 题，因此为了解决固有时间尺度分解算法的模态混叠问题，有必要提出一种新的噪声辅助 分析框架。此外，信号分解之后会产生许多旋转分量，如何从这些旋转分量选取对噪声敏感 分量并进行故障特征的提取也是滚动轴承故障诊断有待解决的问题。
【发明内容】
[0004] 本发明所要解决的技术问题是，提供一种能够准确分析非平稳振动信号，提取具 有代表性故障特征的基于集合固有时间尺度分解算法的滚动轴承故障诊断方法。
[0005] 本发明所采用的技术方案是：一种基于集合固有时间尺度分解算法的滚动轴承故 障诊断方法，包括以下步骤：
[0006] 1)利用位移传感器采集滚动轴承的振动信号x(t);
[0007] 2)利用集合固有时间尺度分解算法对采集到的振动信号x(t)进行分解，生成若 干个旋转分量PRCW)和残差信号e (t);
[0009] 其中，K为旋转分量:PRCjicCO的总数，k为旋转分量标号；
[0010] 3)从所有旋转分量中选取能够反映故障信息的敏感旋转分量；
[0011] 4)对敏感旋转分量进行包络谱分析；
[0012] 5)通过分析故障特征频率对应的包络谱幅值判断故障类型。
[0013] 步骤2)中所述的集合固有时间尺度分解算法，包括首先定义算子Hk( ·)为利用 固有时间尺度分解算法获取信号的第k阶旋转分量，进行如下步骤：
[0014] (1)初始化M个幅值为A的白噪声信号ni(t)，i = 1，···，Μ，其中M为偶数，ni(t) =_nM/2+i ⑴；
[0015] (2)将白噪声信号ni(t)分别添加到振动信号x(t)中生成M个混合信号 x(t)+ni(t)；
[0016] (3)利用固有时间尺度分解算法计算各混合信号X(t)+ni(t)的第一阶旋转分量 PRCii(〇,并对所有第一阶旋转分量PRC>：)求平均值，将所述平均值定义为的集合固有时间 尺度分解算法的第一阶旋转分量
[0020] (5)利用固有时间尺度分解算法计算各ei(t) + a i · Hjr^a))信号的第一阶旋转 分量H1 (ei (t) + a i · H1 (Ii1⑴）），并对所有第一阶旋转分量H1 (ei (t) + a i · H1 (Ii1⑴））求平均 值，将所述平均值定义为的集合固有时间尺度分解算法的第二阶旋转分量
[0022] (6)对于k = 2, 3,…，K，计算第k阶残差信号
[0024] (7)利用固有时间尺度分解算法计算各+ 信号的第一阶旋转 分量H1 (ek (t) + a i · Hk (Ii1⑴）），并对所有第一阶旋转分量H1 (ek (t) + a i · Hk (Ii1⑴））求平均 值，将所述平均值定义为的集合固有时间尺度分解算法的第k+Ι阶旋转分量
[0026] (8)令k = k+l，重复步骤（5)~（7)直到第k阶残差信号e(t)满足固有时间尺 度分解算法的终止条件，最终振动信号x(t)表示为
[0028] 其中的系数a k是采用如下公式得到：
[0029] a k= ε · std (e k (t))/std (Hk (n； (t)))
[0030] 其中，ε为白噪声幅值A与信号标准偏差的比值，std(ek(t))为信号e k(t)的标 准偏差。
[0031] 步骤3)中所述的敏感旋转分量，是所有旋转分量中峭度值最高的分量。
[0032] 本发明的基于集合固有时间尺度分解算法的滚动轴承故障诊断方法，解决了固有 时间尺度分解算法的模态混叠问题，为特征提取奠定了良好的基础，通过计算峭度，选取了 对故障较为敏感的旋转分量，最后通过分析故障特征频率对应的敏感分量包络谱幅值判断 故障类型。本发明可以准确的识别滚动轴承故障，适用于滚动轴承故障诊断。
【附图说明】
[0033] 图1是本发明基于集合固有时间尺度分解算法的滚动轴承故障诊断方法流程图；
[0034] 图2是仿真信号y(t)及其两个组成部分Xl(t)和知(〇时域波形图；
[0035] 图3是利用固有时间尺度分解算法对仿真信号y(t)进行分解的结果图；
[0036] 图4是利用集合固有时间尺度分解算法对仿真信号y(t)进行分解的结果图；
[0037] 图5是故障实验系统示意图；
[0038] 其中
[0039] 1 :电动机；2 :联轴器；3 :轴承座；4 :圆盘；5 :转速传感器；6 :位移传感器；7 :计算 机；8 :数据采集卡；9 :试验台基座
[0040] 图6是滚动轴承滚动体故障振动信号波形图；
[0041] 图7是滚动轴承滚动体故障振动信号经过集合固有时间尺度分解算法分解后得 到的旋转分量和残差信号波形图；
[0042] 图8是敏感旋转分量的包络谱图。
【具体实施方式】
[0043] 下面结合实施例和附图对本发明的一种基于集合固有时间尺度分解算法的滚动 轴承故障诊断方法做出详细说明。
[0044] 本发明的一种基于集合固有时间尺度分解算法的滚动轴承故障诊断方法，如图1 所示，包括以下步骤：
[0045] 1)利用位移传感器采集滚动轴承的振动信号x(t);
[0046] 2)滚动轴承振动信号属于非平稳信号，因此选择最先进的非平稳信号分析方 法-固有时间尺度分解算法对轴承振动信号进行分析。和经验模态分解一样，固有时间尺 度分解算法存在着模态混叠问题。因此，本发明提出了一种集合固有时间尺度分解算法，并 利用该算法对采集到的振动信号x(t)进行分解，生成若干个旋转分量和残差信号 e(t)；
[0048] 其中，K为旋转分量PRC少)的总数，k旋转分量标号；
[0049] 所述的集合固有时间尺度分解算法，包括首先定义算子Hk( ·)为利用固有时间尺 度分解算法获取信号的第k阶旋转分量，进行如下步骤：
[0050] (1)初始化M个幅值为A的白噪声信号ni(t)，i = 1，···，Μ，其中M为偶数，叫⑴ =_nM/2+i ⑴：
[0051] (2)将白噪声信号叫⑴分别添加到振动信号x(t)中生成M个混合信号 x(t)+ni(t)；
[0052] (3)利用固有时间尺度分解算法计算各混合信号x(t)+Ii1 (t)的第一阶旋转分量 PRC')，并对所有第一阶旋转分量PRC^)求平均值，将所述平均值定义为的集合固有时间 尺度分解算法的第一阶旋转分量
[0054] (4)计算一阶残差信号
[0056] (5)利用固有时间尺度分解算法计算各ei(t) + a i ·氏㈨⑴）信号的第一阶旋转