滚动轴承故障诊断(附MATLAB程序)

第二组实验轴承故障数据：Test2.mat数据打开后应采用X105_DE_time作为分析数据，其他可作为参考，转速1797rpm轴承型号：6205-2RSJEMSKF,深沟球轴承采样频率：12kHz1、确定轴承各项参数并计算各部件的故障特征频率通过以上原始数据可知次轴承的参数为：轴承转速r=1797r/min；滚珠个数n=9；滚动体直径d=7.938mm；轴承节径D=39mm；：滚动体接触角。=0由以上数据计算滚动轴承不同部件故障的特征频率为：外圈故障频率f1=r/60*1/2*n(1-d/D*cosa)=107.34Hz内圈故障频率f2=r/60*1/2*n(1+d/D*cosa)=162.21Hz滚动体故障频率f3=r/60*1/2*D/d*[1-(d/D)八2*cos八2(a)]=70.53Hz保持架外圈故障频率f4=r/60*1/2*(1-d/D*cosa)=11.92Hz.对轴承故障数据进行时域波形分析将轴承数据Test2.mat导入MATLAB中直接做FFT分析得到时域图如下：

并求得时域信号的各项特征：(1)有效值:0.2909； (2)峰值:1.5256；(3)峰值因子:5.2441；(4)峭度:5.2793；(5)脉冲因子:7.2884；(6)裕度因子:9.1083：.包络谱分析对信号做EMD模态分解，分解得到的每一个IMF信号分别和原信号做相关分析，找出相关系数较大的IMF分量并对此IMF分量进行Hilbert变换。EmpiricalModeDecomposition离；岫佃3<|和匕3if,*邢*u寸Am■中1见#34山仲叫储仲/心中冲陛小咖产田的网1—.A1尸仲也阿团为《灼《叫删产Il 椒七3小；,”*咒小心―卡$，>^词［)(^*4/卜—中信队1咻1*一看单1>"^*¥*.*加川《沁咻；gz.Mw—CUE J~~Z'/Q'S"c~~~N\1—CUE J~~Z'/Q'S"c~~~N\1铲7'f~f由图中可以看出经过EMD分解后得到的9个IMF分量和一个残余量。IMF分量分别和原信号做相关分析后得出相关系数如下：EMD分量IMF1IMF2IMF3IMF4IMF5相关系数0.95960.19900.10960.00620.0230EMD分量IMF6IMF7IMF8IMF9IMF10相关系数0.00320.00450.00550.00600.0062由上表得：IMF1的相关系数明显最大，所以选用IMF1做Hilbert包络谱分析。所得Hilbert包络谱图如下:

对包络谱图中幅值较大区域局部放大得到下图由以上包络图的局部放大图中可以看出包络图中前三个峰值最大也最明显，三个峰值频率由小到大排列分别为58.59Hz、105.5Hz、164.1Hz。把这三个频率数值和前文计算所得的理论值进行比较可知:频率值最大为164.1Hz和内圈的故障理论计算特征频率f2=162.21Hz相近，说明此轴承的故障发生在轴承的内圈。clc程序1原始信号时域分析及小波去噪处理轴承诊断幅值原信号时域图小波去噪幅值小波去噪后时域图程序2分解及 包络轴承诊断为一阶矩阵，是分量，为数据点实信号的信号重构及 结果显示函数峭度均值包络包络频率幅值、、一,.m差均方值（方根:计算偏度为峭度