滚动轴承故障诊断(附MATLAB程序)
更新时间:2024-01-31 05:42:01 阅读量: 教育文库 文档下载
- 滚动轴承的故障诊断推荐度:
- 相关推荐
第二组实验 轴承故障数据:
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/60 * 1/2 * n(1+d/D *cosα)=162.21Hz 滚动体故障频率f3=r/60*1/2*D/d*[1-(d/D)^2* cos^2(α)]=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.包络谱分析
对信号做EMD模态分解,分解得到的每一个IMF信号分别和原信号做相关分析,找出相关系数较大的IMF分量并对此IMF分量进行Hilbert变换。
Empirical Mode Decompositionres.imf8imf7imf6imf5imf4imf3imf2imf1signal
由图中可以看出经过EMD分解后得到的9个IMF分量和一个残余量。IMF分量分别和原信号做相关分析后得出相关系数如下:
EMD分量 相关系数 EMD分量 相关系数 IMF1 0.9596 IMF6 0.0032 IMF2 0.1990 IMF7 0.0045 IMF3 0.1096 IMF8 0.0055 IMF4 0.0062 IMF9 0.0060 IMF5 0.0230 IMF10 0.0062 由上表得:IMF1的相关系数明显最大,所以选用IMF1做Hilbert包络谱分析。所得Hilbert包络谱图如下:
对包络谱图中幅值较大区域局部放大得到下图
由以上包络图的局部放大图中可以看出包络图中前三个峰值最大也最明显,三个峰值频率由小到大排列分别为58.59Hz、105.5Hz、164.1Hz。把这三个频率数值和前文计算所得的理论值进行比较可知:频率值最大为164.1Hz和内圈的故障理论计算特征频率f2=162.21Hz相近,说明此轴承的故障发生在轴承的内圈。
clc
程序1:原始信号时域分析及小波去噪处理
clear all
z=importdata('C:\\Users\\wangkun\\Desktop\\轴承诊断\\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); xlabel('t');
ylabel('幅值');
title('小波去噪后时域图')
程序2:EMD分解及Hilbert包络
clc
clear all
z=load('C:\\Users\\wangkun\\Desktop\\轴承诊断\\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)=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('幅值');
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('幅值');
正在阅读:
滚动轴承故障诊断(附MATLAB程序)01-31
机电系统微机接口技术大作业12-17
三角恒等变换讲义01-27
德育校本课程研究方案140506-01
2015最新苏教版一年级数学上册全册导学案(45页)05-31
利器:HLSL起步教程-完整篇05-19
打油诗_告别小磨蹭02-23
压迫与反抗的结晶 朱仁裕 - 图文11-01
高频电子线路实验报告10-17
第八章单元测验03-18
- exercise2
- 铅锌矿详查地质设计 - 图文
- 厨余垃圾、餐厨垃圾堆肥系统设计方案
- 陈明珠开题报告
- 化工原理精选例题
- 政府形象宣传册营销案例
- 小学一至三年级语文阅读专项练习题
- 2014.民诉 期末考试 复习题
- 巅峰智业 - 做好顶层设计对建设城市的重要意义
- (三起)冀教版三年级英语上册Unit4 Lesson24练习题及答案
- 2017年实心轮胎现状及发展趋势分析(目录)
- 基于GIS的农用地定级技术研究定稿
- 2017-2022年中国医疗保健市场调查与市场前景预测报告(目录) - 图文
- 作业
- OFDM技术仿真(MATLAB代码) - 图文
- Android工程师笔试题及答案
- 生命密码联合密码
- 空间地上权若干法律问题探究
- 江苏学业水平测试《机械基础》模拟试题
- 选课走班实施方案
- 滚动轴承
- 诊断
- 故障
- 程序
- MATLAB
- 自然科学基金 F.信息科学部
- 桂花
- 第七篇 矿井防治水技术管理规定(最终)
- 《普通遗传学》各章试题及答案
- 合规测试试卷(答案)
- 合肥市建筑工程安全风险评估综合表
- 计算机网络复习资料 - 图文
- 舒立平 终端店铺培训20 - 图文
- 广元女儿节
- 探讨波的衍射在生活中的应用(波的衍射专题实验科技论文)
- 安徽2000定额解释
- 台山市水利局行政执法职权核准界定表
- 安庆师范学院2013—2014学年获得奖学金学生名单
- 干部驻村帮扶要做到“三到位”
- 2012高考语文一轮复习(人教版)专题11 考点过关1
- 人事干部个人工作总结(三篇)与人事干部个人年终总结汇编
- 环保认证证书和标志管理办法
- 地质雷达检测报告 - 图文
- 2018-2019年合肥市小升初数学模拟试题(共10套)附详细答案
- 一辩稿-伦理是市场的禁区