已阅读5页,还剩1页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
第二组实验轴承故障数据: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 内圈故障频率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.包络谱分析对信号做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包络谱图如下:对包络谱图中幅值较大区域局部放大得到下图由以上包络图的局部放大图中可以看出包络图中前三个峰值最大也最明显,三个峰值频率由小到大排列分别为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=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.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,:); %计算偏度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(
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 业务员雇佣合同
- 长期股权投资应用案例-关于持有期间采用权益法核算的会计处理
- 历史丨浙江省2023学年第二学期精诚联盟适应性联考历史试卷及答案
- 书中岁月:借阅制度
- 传染病预防:信息上报篇
- 创意城物业品质内容及管理办法
- 景县广川浩瀚钢结构加工厂年产1000吨新上铁塔生产线项目环评报告
- 房地产投资多元化发展策略
- 油松病虫害防治的最佳时期
- 办公室防疫措施及流程图
- 招标文件(风口风阀、阀门、消声器静压箱)
- 广东省建筑与装饰工程综合定额(说明及工程量计
- 《黑猫》英文版 BY爱伦坡
- 对硝基苯甲醛的技术配方和生产工艺
- (完整版)上海牛津版二年级下册英语全册教材内容
- 六年级语文总复习课《修改病句》修改课件(精)
- 饭店等级四种评定模式
- 霍曼门限位设定
- 初中音乐-《中国人》课件PPT课件
- 肺叶切除术后护理全肺切除术后护理.ppt
- 香港的营养资料标签制度食物环境卫生署
评论
0/150
提交评论