版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、滚动轴承故障诊断(附MATLAB程序)第二组实验轴承故障数据:Test2.mat 数据打开后应采用X105_DE_time作为分析数据,其他可作为参考,转速1797rpm轴承型号:6205-2RS JEM SKF, 深沟球轴承采样频率:12k Hz1、确定轴承各项参数并计算各部件的故障特征频率通过以上原始数据可知次轴承的参数为:轴承转速r=1797r/min;滚珠个数n=9;滚动体直径d=7.938mm; 轴承节径D=39mm;:滚动体接触角=0由以上数据计算滚动轴承不同部件故障的特征频率为: 外圈故障频率f1=r/60 * 1/2 * n(1-d/D *cos)=107.34Hz 内圈故障频
2、率f2=r/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.92Hz2. 对轴承故障数据进行时域波形分析 将轴承数据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变换。由图中可以看出经过EMD分解后得到的个IMF分量和一个残余量。IMF分量分别和原信号做相关分析后得出相关系数如下:EMD分量IMF1IMF2IMF3IMF4IMF5相关系数0.95960.19900.10960.00620.0230EMD分量IMF6IMF7IMF8IMF9IMF10相关系数0.00320.00450.00550.00600.0062由上表得:IMF1的相关系数明显最大,所以选用IMF1做Hilbert包络谱分析。所得Hilbert包络谱图如下:
4、对包络谱图中幅值较大区域局部放大得到下图由以上包络图的局部放大图中可以看出包络图中前三个峰值最大也最明显,三个峰值频率由小到大排列分别为58.59Hz、105.5Hz、164.1Hz。把这三个频率数值和前文计算所得的理论值进行比较可知:频率值最大为164.1Hz和内圈的故障理论计算特征频率f2=162.21Hz相近,说明此轴承的故障发生在轴承的内圈。clc程序1:原始信号时域分析及小波去噪处理clear allz=importdata(C:UserswangkunDesktop轴承诊断test2.mat);x1=z.X105_DE_time(1:4096);clear z;N=4096;fs=
5、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包络clcclear allz=load(C:UserswangkunDesktop轴承诊断test2.m
6、at);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)
7、=std(imf(i,:).0.5; %均方根值 f(i)=skewness(imf(i,:); %计算偏度endk,c=max(a); %k为峭度最大值,c为最大元素在数组中的位置r,lags=xcorr(x,lag,unbiased); %计算序列的自相关函数for i=1:mR,lags=xcorr(imf(i,:),lag,unbiased); %计算序列的自相关函数a=corrcoef(R(1:N/2),r(1:N/2); %相关系数矩阵【对称】,主对角元素为1xg(i)=abs(a(1,2); %相关系数endR,C=max(xg); %R为最大值,C为最大元素在数组中的位置 figure(4);y = hilbert(imf(C,:);a = abs(y);%包络b=fft(a);mag1=abs(b);m
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 做文明有礼的学生演讲稿
- 考前易错小题狂做02:重难单词短语句型强化练100道-2023-2024学年八年级英语下学期期末考点大串讲(牛津译林版) (解析版)
- 中学生读林徽因传有感
- 幼儿英语故事活动教学教案11篇
- 数学实践报告
- 淋巴水肿课程考题及答案
- 智能制造装备单元系统集成 课件 6.4 机器人系统规划与调整-
- 4.2 透镜(考点解读)八年级物理上册同步考点解读与专题训练(苏科版)(解析版)
- 五年级数学一二单元备课教案
- XX市关于进一步加快再生资源回收体系建设的实施方案
- 第六单元综合性学习《以和为贵》课件(共16张PPT) 部编版语文八年级下册
- 2023年济南历下控股集团有限公司招聘笔试题库及答案解析
- 《变形记》卡夫卡课件
- 2023年大连国际机场集团有限公司招聘笔试模拟试题及答案解析
- 一元一次不等式单元测试试题(含答案)
- NY∕T 3349-2021 畜禽屠宰加工人员岗位技能要求
- 张夫也外国工艺设计史
- 管制刀具认定标准
- 11、交付验收标准及方案
- 小学语文教师评职称专业技术工作总结4篇
- 一年级下册语文阅读试题-《小老虎历险记》课外阅读题-人教部编版(含答案)
评论
0/150
提交评论