第 39 卷 第 5 期
2019 年 5 月
电 力 自 动 化 设 备
Electric Power Automation Equipment
Vol.39 No.5
May 2019
基于 COT⁃SSD 的变转速滚动轴承微弱故障诊断
王晓龙,唐贵基,何玉灵
(华北电力大学 机械工程系,河北 保定 071000)
摘要:针对背景噪声干扰及转速波动工况下滚动轴承微弱故障识别困难这一问题,提出一种结合计算阶次追
踪( COT)和奇异谱分解( SSD)的新型诊断方法。 利用 COT 算法对采集的原始时域信号进行等角度重采样,
继而利用 SSD 算法对重采样角域信号进行处理,通过自适应构建的轨迹矩阵的奇异值分解重组,将角域信号
从高频至低频分解为若干个奇异谱分量,利用融合峭度指标筛选最佳奇异谱分量,选定最佳分量后对其进行
进一步的包络解调运算,最终通过分析包络阶次谱中幅值突出的成分来准确判定滚动轴承运行状态。 滚动
轴承内外圈故障实测信号分析结果表明,所提方法能够有效提取出变速工况下滚动轴承的微弱故障特征
信息。
关键词:滚动轴承;变转速;微弱故障;计算阶次追踪;奇异谱分解
中图分类号:TM 315;TH 17 文献标识码:A DOI:10.16081 / j.issn.1006
-
6047.2019.05.028
收稿日期:2018
-
04
-
01;修回日期:2019
-
03
-
20
基金项目:国家自然科学基金资助项目(51777074) ;中央高
校基本科研业务费专项资金资助项目(2018MS124,2018YQ03)
Project supported by the National Natural Science Foundation of
China(51777074) and the Fundamental Research Funds for the
Central Universities(2018MS124,2018YQ03)
0 引言
滚动轴承是风力发电机的重要组成零部件,在
其故障萌芽阶段,损伤造成的振动冲击能量较弱,加
之传递路径衰减效应及背景噪声的干扰影响,使得
传感器获取的振动信号的信噪比往往较低。 此外,
风机启停和负载变化时,转速会发生波动,使得采集
的信号源具有明显的非平稳特性,利用定速轴承诊
断方法进行分析无法有效提取变转速过程信号中的
状态特征信息,因此研究适用于滚动轴承变转速工
况下的诊断方法具有重要的现实意义
[1⁃3]
。
滚动轴承故障振动信号为多分量调制信号,在
解调分析前需要将其分解为若干个单分量调制信
号,现有 的 诸 多 非 平 稳 信 号 处 理 手 段, 如 小 波 变
换
[4]
、经验模态分解
[5]
、局部均值分解
[6]
、局部特征
尺度分解
[7]
、变分模态分解
[8]
等方法虽然在滚动轴
承诊断中已得到较为广泛的应用,但是依然存在一
些未能从根本上得以解决的缺陷或不足。 Bonizzi 等
人于 2014 年提出一种新的自适应信号处理方法 ———
奇异谱分解 SSD(Singular Spectrum Decomposition)
[9]
法。 SSD 算法建立在奇异谱分析理论的基础上,其
通过重新定义轨迹矩阵表达形式来增强信号中的振
荡分量,降低重构残余项的能量,在运算过程中能够
根据所分析信号的特点自动选取轨迹矩阵的嵌入维
数,对轨迹矩阵的主成分进行奇异值分解及特征重
组,从而实现信号频带的自适应划分,得到一系列频
率范围由高至低的单分量信号。 SSD 算法的提出为
多分量 旋 转 机 械 故 障 非 平 稳 信 号 处 理 提 供 了 新
思路。
当滚动轴承在变转速工况下运行时,其调制信
号中的调制频率是与转速相关的时变量,直接对解
调后的信号进行频谱分析会出现“频率涂抹” 现象,
将极大地影响故障特征提取效果,而计算阶次跟踪
COT(Computed Order Tracking) 算法则是面向该问
题的一种有效应对措施
[10]
。 COT 算法分析精度高,
无需特定硬件支持,应用较为广泛,通过软件计算即
可追踪参考轴的转速波动并进行等角度增量重采
样,可将时域非平稳信号转换为角域平稳信号,从而
提取信号中与转速有关的信息,同时抑制与转速无
关的干扰成分。
鉴于上述分析,本文尝试将 SSD 算法引入到变
转速轴承故障诊断领域中,并将其与 COT 算法相结
合,提出了基于 COT⁃SSD 的滚动轴承微弱故障诊断
方法,有望从信噪比较低的信号源中提取出敏感特
征成分,实现滚动轴承损伤的精确判定。
1 COT 算法
同时对振动信号 x( t) 和键相信号 s(t) 进行等
时间间隔采样,识别键相信号 s(t)的脉冲时刻 t
zi
,即
参考轴每旋转 1 周对应的时刻,对脉冲时刻 t
zi
进行
差分运算可得到参考轴每旋转 1 周所需时间 T
i
。 假
设参考轴在旋转 1 周的时长内做匀变速运动,并假
定旋转前 2 周的角加速度相等,根据运动学定律可
计算出参考轴前 2 周的初速度 v
10
和角加速度 a
1
[11]
:
v
10
T
1
+
1
2
a
1
T
2
1
=
2π
v
10
(T
1
+
T
2
)
+
1
2
a
1
(T
1
+
T
2
)
2
=
4π
ì
î
í
ï
ï
ï
ï
(1)
设参考轴每旋转 1 周等角度采集 N 个点,即等