基于奇异值差分谱与改进包络分析的轴承故障特征提取

本资料包含pdf文件1个，下载需要1积分

Bearing Fault Feature Extraction Based on Diference Spectrum of SingularValue and Improved Envelope AnalysisYang Wang-can，Zhang Pei-lin，Wang Huai-guang，Chen Yan-long(Seventh Department，Ordnance Engineering Colege，Shijiazhuang 050003，China)Abstract：To overcome the dificulty of extracting fault feature from roling bearing vibration signal，a method is putforward based on diference spectrum of singular value and improved envelope analysis.Firstly，the original bearing vi-bration signal is decomposed into a series of component signals which is able to be added linearly by singular value de-composition(SVD).Due to the diference of singular values between the component signals of fault feature and compo-nent signals of noise，the efective singular values are selected according to the propey of difference spectrum of singu-lar value and then the component reconstruction signals containing fault feature are selected.Aiming at residual noiseafter denoised by singular value decomposition，the improved envelope analysis is adopted to fuher remove residualnoise in reconstructed signal in frequency domain.Finally，actual inspected bearing signal is analyzed，the bearingfault signal with obvious fault feature and prominent fault frequency is extracted accurately and then fault diagnosis iscompleted。

Key words：roling bearing；diference spectrum of singular value；improved envelope analysis；feature extraction滚动轴承是旋转机械中应用广泛但也易损坏的零件，据不完全统计，旋转机械的故障约有30%是因滚动轴承发生故障引起的 J。因此，对轴承的状态进行监测和故障诊断显得极为重要，如何有效地从混有强烈噪声的原始振动信号中提取到轴承故障特征是轴承故障诊断的关键。

收稿日期：2012-08-27；修回日期：2012-11-25作者简介：杨望灿(1988-)，男，硕士研究生，主要研究方向为信号处理与故障诊断。

奇异值分解(Singular Value Decomposition，SVD)能够将信号分解为-系列线性分量的叠加，具有零相移、波形失真小等优点 J，在信号处理和图像去噪等方面得到了广泛的应用 J。奇异值分解在进行故障信号分析时关键在于如何合理地选择有效奇异值来重构信号，并提取故障特征r 。

奇异值分解的特点决定了原始信号经奇异值分解重构后，重构信号中仍不可避免地存在残余噪声信号，造成对重构信号包络解调分析时仍存在干扰频率，影响轴承的故障特征提韧故障诊断。

《轴承)2013.No.5针对此问题 ，提出了-种基于奇异值差分谱和改进包络分析的轴承故障特征提取方法。首先采用奇异值分解方法对信号进行分解，然后根据奇异值差分谱图中最大峰值的位置确定含有轴承故障特征的分量所对应的奇异值来重构信号，最后利用改进包络分析在频域里将重构信号中的残余噪声消除，凸显轴承故障频率，完成轴承的故障特征提取，诊断轴承故障。

1 奇异值差分谱理论若A是 m X n的矩阵(假设 m>n)，A的秩为r(ro-，>0，U和y的列向量 lf 和 l 分别为 A的左、右奇异向量。

所以(1)式又可以表示为A∑ U l， 。 (2)由(2)式可得，矩阵A也可以说是以非零奇异值为权重，相对应的左、右奇异向量做外积后的加权和，可见奇异值的选取是矩阵A重构的关键。

假设采集到的离散信号X[ (1)， (2)，，(N)]， (i)s( ) (i)，S(i)为有用信号，u(i)为噪声信号，利用信号 可以构造 Hankel矩阵A (1)(2)：(2)X(3) ：(N-n1) (N-r/，2) pc(npc(n1 (Ⅳ)(3)其中，1

观察 Hankel矩阵，下-行的数据比上-行滞后-个数据点。对于无噪声信号，由于相邻两行具有相关性，这种信号构造的 Hankel矩阵是-种病态矩阵，前面 个奇异值较大，后面的接近于零，奇异值有明显的突变，而对于噪声信号，相邻两行数据没有相关性，其构造的 Hankel矩阵的奇异值没有明显的突变，所以用奇异值差分谱 描述含噪信号的奇异值突变，从而选取有效的奇异值重构信号。

设c为奇异值由大到小形成的序列，C[ 1， 2，， ]。定义b f， - (i1，2，，-1)， (7)则所有的b 构成序列 B[b ，b ，，b ]，序列日即为奇异值的差分谱。差分谱图反应了相邻奇异值之间的变化情况，最大峰值处代表了奇异值序列的最大突变，由于有用信号和噪声的相关性不同而导致2种信号在奇异值上表现出最大的差异，这个最大突变点掖映了有用信号和噪声信号的分界，所以，利用奇异值差分谱能够自动地确定有用分量的个数。

同时注意到在构造上述 Hankel矩阵时，有用信号和噪声信号同时包含其中。由于噪声信号的Hankel矩阵为良态满秩矩阵，经过上述奇异值分解后，噪声信号分布在每个奇异值对应的分量中，提取到的有用分量中仍然不可避免地存在着噪声信号，这也是由奇异值分解的特点决定的，因此需要对重构后的信号进-步做去噪处理，提取信号特征。

(n J 、( 1)l 2 改进包络分析； l 在工程实际中，滚动轴承的故障振动信号多，Xi(Ⅳ) 为受到背景噪声和其他振动频率干扰的幅度调制杨望灿，等：基于奇异值差分谱与改进包络分析的轴承故障特征提取得到的重构信号中仍含有残余噪声，而包络解调分析方法本身不具有去噪的功能，而且易受噪声干扰，所以采用改进包络分析，在频域中进-步去除残余噪声，减小干扰频率的影响，这样对处理后的信号再进行包络解调，可以使得轴承的故障频率突出。

对重构信号进行频域变换后残余噪声在幅值谱上会表现为幅值较小的干扰频率，通过对 Shan-non熵的定义的修改 ，利用改进 Shannon熵函数消除这些干扰频率，再变换到时域后，残余噪声即可被消除。

定义改进的Shannon熵函数为( )P( )l og：P( )l，其中，i1，2，，m，m为幅值谱的点数。

功率谱 P( )定义为- - P(f)P( ) ，(8)归-化(9)u[∑(1 P( )- 1∑P( )) /I 1 E 1(m-1)] 。 (10)上述处理完成后，重构信号的幅值谱中幅值高的主要频率成分突出，幅值低的干扰频率被去除。这样，频域去噪后的信号变换成时域信号再进行包络解调，包络谱中就能够突出轴承故障频率，避免残余噪声信号的影响，清晰有效地提取轴承的故障特征。

3 轴承故障信号分析在-个齿轮箱上采集得到实际的滚动轴承信号，轴承型号为SKF 6205，在测试轴承内圈上利用电蚀加工单点故障，加速度传感器安装在驱动端的轴承座上。采样数据是在电动机转速为1 797r/min下采集到的，采样频率为 12 kHz，采样点数为 1 024。根据内圈故障计算公式可得内圈单点故障频率为162.2 Hz。实测内圈故障信号的时域波形和频谱图如图1和图2所示。

对于轴承而言，当其滚道有损伤时，其时域信号会出现周期性的冲击信号，由于受到轴承共振调制的影响，其特征频率主要集中在 中、高频率段。观察图 1和图2，无法辨别出轴承是否发生故障，图2中虽然中频段幅值较大，但整个频带都存在频率成分，轴承故障信号受噪声干扰严重，无法获取轴承振动情况的准确信息。利用采集到的信号构造Hankel矩阵，进行奇异值分解。采集到的0 0.0l 0.02 0.o3 0.04 0.05 0.05 0.07 0.08时间t/s图l 内圈故障的时域波形0 l 00o 2 000 3 O0o 4 000 5 000 6 00o频f/Hz图2 内圈故障的频谱 图离散数据点为 1 024个，为了充分利用数据，构造行数为 513，列数为 512的最大阶数的 Hankel矩阵，得到512个奇异值。奇异值序列和奇异值差分谱如图3和图4所示。

12108曼 6蚤42图3 奇异值序列0 100 200 300 400 500 600坐标序列i//。

图4 奇异值差分谱0 - s.II)/《 罂 - 量墨 坚 - s.II- 肇《轴承2013.№.5从图3可以看出奇异值 的变化趋势，从图 4可以看出奇异值的突变主要发生在奇异值序列的前部，为了清楚显示，将奇异值差分谱前60个数据局部放大，如图5所示。

图5 奇异值差分谱 的前 60个数据 点从图5可以清晰地看到第4个坐标点的峰值最大，这样根据奇异值差分谱理论可以自动地确定重构信号的有用分量的个数为4，完成对原始信号的重构，去除噪声，提取故障轴承的调幅特征。

O.60.4O.20、基-O.2《 西-0.4坚- 0.6- OI80 0.Ol 0. 0.O3 0.o4 0.)5 0.06 0.07 0.08时间t/s图6 重构信号的时域波形 ，. 。

1 0o0 2 000 3 000 4 000 5 000 6 00o频f/Hz图7 重构信号的频谱图从图6可以看到故障轴承信号的调幅特征，由于内圈存在凹坑，使得轴承运转过程中滚动体受到凹坑的撞击产生周期性冲击，产生了调幅特征。在图7所示的频谱图中重构信号的频率主要集中在3 kHz附近，对照图2可以看到噪声信号产生的其他频率段的大部分干扰频率被剔除，主要频率成分比较突出。但是由于奇异值分解的特点，重构信号中仍保留了部分残余噪声。从图 7中可以看到变换到频域后重构信号的主要频率成分附近还存在-些低幅值的干扰频率。残余噪声在对信号做 Hilbert包络解调时会产生较大影响，因此需要进-步去除。对奇异值分解后的重构信号直接做 Hilbert包络解调，结果如图 8所示。采用改进包络分析方法进行频域去噪，突出幅值高的主要成分，削减低幅值的干扰成分，去除残余噪声 ，然后变换为时域信号做 Hilbert包络解调，结果如图 9所示 。

图8 重构信号的包络谱图9 重构信号的改进包络谱图8中除了故障频率外，仍然受到杂频的干扰。从图9中可以看到 162 Hz的频率非常突出，而且几乎没有干扰频率，与内圈的理论故障频率162.2 Hz非秤近，可以认为提取到的特征频率是内圈的故障频率。

4 结束语研究了基于奇异值差分谱与改进包络分析的轴承故障特征提取方法，奇异值差分谱理论克服了以经验为主选取奇异值重构信号的弊端，利用奇异值差分谱的最大峰值自动选取奇异值重构原始信号，既去除了实测信号中的大量噪声，又保留- s.曼 迥馨加 ∞ 如 ∞ 如 加 m 0- s.量 覃馨 二 !CN41-1148/TH轴承 201 3年5期 Bearing 2013.No.5 53-54转盘轴承安装孔位置的判定张亚平 ，高聪颖 ，李云峰 ，李媛媛 ，王艳(1.洛阳轴研科技股份有限公司，河南 洛阳 471039；2.河南科技大学，河南 洛阳 471003)摘要：测量转盘轴承上安装孑L的孔径及2个相邻孔之间的壁厚可以间接算出弦长，通过弦长误差及弦长误差累积可判断安装孔是否满足使用要求，并结合实例说明了安装孔位置的判定方法。

关键词：转盘轴承；安装孔 ；位置中图分类号：TH133.33 文献标志码：B 文章编号：1000-3762(2013)05-0053-02转盘轴承的安装-般通过安装孑L实现，而产品设计对安装孔都有位置度要求，大尺寸(如直径5 m以上)套圈不能采用数控钻床加工，只能采用传统划线钻孔，这种加工方式带来的误差相对较大，在孔加工完之后，通过测量孔径及 2个相邻孔之间的壁厚可以间接算出弦长，通过弦长可判断安装孔是否满足使用要求。

假设2个相邻孔的孔径分别为A和B，两者之间的壁厚为 C，则弦长为 ( B)/2c。某转盘轴承均布 6O个安装孔，安装孔 中心径为 5 600收稿日期：2012-11-05mm，安装孔直径为45 mm，螺栓为 M42。通过卡尺测量各尺寸，结果见表 1。理论 弦长为 5 600×sin(180。/60)293.08 mm。

常规处理方法是将弦长误差(实际弦长与理论弦长的差值)分别描点绘出，只要最大点和最小点落在-定范围内就认为满足使用要求。该范围- 般炔装孔孔径与螺栓直径差值的-半 ，为保险起见，需再乘以-个小于 1的系数，假定系数为0.8，则本例中的阈值为 1.2。弦长误差如图 1所示 。

从图中可以看出，各点均位于阈值范围内，可了与轴承故障有关的有用信息。针对奇异值分解方法中仍会存留部分残余噪声的问题，利用改进包络分析的方法在频域范围内进-步去除残余噪声，提取到的轴承故障特征频率清晰突出，消除了干扰频率的影响，效果非常明显。