《滚动轴承故障诊断(附MATLAB程序)》由会员分享,可在线阅读,更多相关《滚动轴承故障诊断(附MATLAB程序)(6页珍藏版)》请在金锄头文库上搜索。
1、第二组实验轴承故障数据:Test2.mat 数据打开后应采用X105_DE_time 作为分析数据,其他可作为参考,转速 1797rpm 轴承型号: 6205-2RS JEM SKF, 深沟球轴承采样频率: 12k Hz 1、确定轴承各项参数并计算各部件的故障特征频率通过以上原始数据可知次轴承的参数为:轴承转速 r=1797r/min;滚珠个数 n=9;滚动体直径 d=7.938mm;轴承节径 D=39mm;:滚动体接触角 =0 由以上数据计算滚动轴承不同部件故障的特征频率为:外圈故障频率 f1=r/60 * 1/2 * n(1-d/D *cos )=107.34Hz 内圈故障频率 f2=r/
2、60 * 1/2 * n(1+d/D *cos )=162.21Hz 滚动体故障频率 f3=r/60*1/2*D/d*1-(d/D)2* cos2( )=70.53Hz 保持架外圈故障频率f4=r/60 * 1/2 * (1-d/D *cos )=11.92Hz 2.对轴承故障数据进行时域波形分析将轴承数据Test2.mat导入 MATLAB 中直接做 FFT 分析得到时域图如下:并求得时域信号的各项特征:(1)有效值: 0.2909;(2)峰值: 1.5256;(3)峰值因子: 5.2441;(4)峭度: 5.2793;(5)脉冲因子: 7.2884;(6)裕度因子: 9.1083:3.包络
3、谱分析对信号做 EMD 模态分解,分解得到的每一个IMF 信号分别和原信号做相关分析,找出相关系数较大的IMF 分量并对此 IMF 分量进行 Hilbert 变换。signalEmpirical Mode Decompositionimf1imf2imf3imf4imf5imf6imf7imf8res.由图中可以看出经过EMD 分解后得到的个IMF 分量和一个残余量。 IMF 分量分别和原信号做相关分析后得出相关系数如下:EMD 分量IMF1 IMF2 IMF3 IMF4 IMF5 相关系数0.9596 0.1990 0.1096 0.0062 0.0230 EMD 分量IMF6 IMF7 I
4、MF8 IMF9 IMF10 相关系数0.0032 0.0045 0.0055 0.0060 0.0062 由上表得:IMF1 的相关系数明显最大, 所以选用 IMF1 做 Hilbert包络谱分析。所得Hilbert 包络谱图如下:对包络谱图中幅值较大区域局部放大得到下图由以上包络图的局部放大图中可以看出包络图中前三个峰值最大也最明显,三个峰值频率由小到大排列分别为58.59Hz、105.5Hz、164.1Hz。把这三个频率数值和前文计算所得的理论值进行比较可知:频率值最大为164.1Hz 和内圈的故障理论计算特征频率f2=162.21Hz相近,说明此轴承的故障发生在轴承的内圈。clc程序
5、1:原始信号时域分析及小波去噪处理clear all 轴承诊断 test2.mat); x1=z.X105_DE_time(1:4096); clear z; N=4096; fs=12000; n=0:N-1; t=n/fs; f=n*fs/N; figure(1); plot(t,x1); xlabel(t); ylabel(幅值); title(原信号时域图 ) % 小波去噪 thr,sorh,keepapp=ddencmp(den,wv,x1); xd=wdencmp(gbl,x1,db3,2,thr,sorh,keepapp); figure(2); plot(t,xd); xlab
6、el(t); ylabel(幅值); title(小波去噪后时域图 ) 程序 2:EMD 分解及 Hilbert包络clc clear all 轴承诊断 test2.mat); x=z.X105_DE_time(1:1024); N=1024; fs=12000; n=0:N-1; f=n*fs/N; lag=N; n=0:N-1; t=n/fs; imf=emd(x); m,n=size(imf); %imf为一 m*n阶矩阵, m是 imf 分量, n 为数据点 emd_visu(x,1:length(x),imf,m); %实信号的信号重构及emd结果显示函数 for i=1:m a(i
7、)=kurtosis(imf(i,:);%峭度b(i)=mean(imf(i,:); %均值; c(i)=var(imf(i,:); %方差; d(i)=std(imf(i,:); %均方值e(i)=std(imf(i,:).0.5; %均方根值f(i)=skewness(imf(i,:); %计算偏度 end k,c=max(a); %k为峭度最大值, c 为最大元素在数组中的位置 r,lags=xcorr(x,lag,unbiased); %计算序列的自相关函数 for i=1:m R,lags=xcorr(imf(i,:),lag,unbiased); %计算序列的自相关函数 a=corrcoef(R(1:N/2),r(1:N/2); %相关系数矩阵【对称】,主对角元素为1 xg(i)=abs(a(1,2); %相关系数 end R,C=max(xg); %R为最大值, C为最大元素在数组中的位置figure(4); y = hilbert(imf(C,:); a = abs(y);%包络 b=fft(a); mag1=abs(b); mag=mag1*2/N; f1=(0:N-1)*fs/N; plot(f1(1:N/2),mag(1:N/2); %set(gca,xlim,0,.400); title(包络); xlabel(频率); ylabel(幅值);