基于多重降噪的轴承故障诊断研究

2023-01-07连雄飞王智冲王宇航

技术与市场 2022年12期

连雄飞，马强，王智冲，王宇航，彭冉

(1.河北工程大学机械与装备工程学院，河北邯郸 056038； 2.天津大学力学系，天津 300354； 3.河北省智能工业装备技术重点实验室(河北工程大学)，河北邯郸 056038； 4.天津非线性动力学与控制重点实验室，天津 300354； 5.邯郸熔淬科技有限公司，河北邯郸 056038)

0 引言

滚动轴承有助于机器平稳高速地旋转，且承载能力好，广泛应用于旋转机械中[1]，但由于运行速度高，易发生故障，导致较大的经济损失，因此及时诊断滚动轴承的早期故障非常必要。

轴承信号中包含的大量状态信息，多表现为非线性和非平稳性信号[2]，这类信号重点在于获得其时域局部性质。处理此类信号时，通常通过时域、频域、时频分析3种方法进行分析[3]，其中时频分析应用更广，本文采用经验模态分解[4]的衍生算法——集合经验模态分解(EEMD)便是一种常用的时频分析方法。

旋转机械装置的复杂结构和环境的多样性，导致轴承故障信号淹没在强背景噪声中，故障信号难以提取。本文针对此问题，提出新的降噪方法，以便提取轴承故障信号。首先利用小波阈值对含噪信号初步去噪，然后利用EEMD将信号分解为多个IMF分量，最后对IMF分量进行SVD分解重构得到最终信号，再通过仿真验证，证明了其较于EEMD-SVD而言能够更好地进行降噪处理。

1 基本理论

1.1 小波阈值去噪

小波变换在信号处理方面具有较好的自适应性[5]，且对于信号处理，小波变换引出多种噪声滤除方法[6]，小波阈值去噪便是其中之一，该方法是由donoho[7]提出的。本文根据轴承振动信号的特征，选择紧支集正交小波(daubechies，db)进行阈值去噪。图1为小波去噪基本原理流程图。

图1 小波去噪流程图

小波阈值降噪(WTD)通过小波变换(WT)对加噪信号s(t)进行分解，并为每个尺度设置小波系数，噪声信号系数相较于待处理信号系数较小，且噪声信号主要分布在高频段[8]，因此选取一个合适的λ值作为阈值对WTD降噪尤为重要。经典阈值函数有2种，其公式分别为：

硬阈值：

(1)

软阈值：

(2)

式中，Wj，k为小波系数，λ为设定阈值，sign(*)为符号函数。

1.2 集合经验模态分解

经验模态分解(EMD)方法处理非线性信号时的一个缺陷就是特征模态函数(IMF)会出现混叠状态，即在分解过程中不同的IMF分量会发生渗透，不同频率会出现在同一IMF分量中，集成经验模态分解算法( EEMD) 则是为了消除、抑制这种混叠状态在EMD的基础上改进而来。已有大量学者将EEMD用于机械振动信号的降噪研究。 EEMD原理是假设在x(t)加入高斯白噪声，对分解得到的IMF分量进行求平均。削减了原信号的间歇性现象，抑制EMD的模态混叠，提高了信号的分解效率和准确性[9]。EEMD 步骤如下。

1)在x(t)中加入高斯白噪声Gi(t)，获得综合信号X1(t) 即:

X1(t)=x(t)+Gi(t)

(3)

2)EMD算法将X1(t)分解得到多个IMF分量。

3)在x(t)中加入不同的噪声，进行(1)(2)步骤N次，求平均，最终得到的IMF 分量:

(4)

式中:xk(t)为对x(t)进行集合经验模态分解得到的第k个IMF分量。

(4)最终，得到原始信号重构信号X(t)。

1.3 奇异值分解

由数学线性理论可知，对于一个实矩阵C，可以构成Hankel矩阵：

(5)

对C进行奇异值分解可得：

C=UΣVT

(6)

U、V为正交矩阵，Σ是奇异值构成的对角矩阵，满足Σ=diag(σ1，σ2，…，σp)、σ1≥σ2≥…≥σp≥0[10]。C为待处理矩阵，左右奇异向量分别由U、V列向量组成。奇异值分解具有良好的尺度和稳定性，也是矩阵固有的特征。U、V满足条件：UUT=I、VVT=I、p=min(m，n)。利用信号能量的可分性，选择合适奇异值，按照式(5)进行SVD逆变换，得到降噪后信号。