Rolling element bearing diagnostics using the Case Western Reserve University data-学习笔记

论文对比分析了三种诊断技术在CWRU数据集上的应用,重点在于滚动轴承的故障诊断。关键步骤包括使用平方包络频谱作为诊断工具,以及通过不同的预处理技术提取故障信息。轴承故障诊断依赖于通过幅度解调得到的包络信号,该信号通常包含比原始信号更明显的故障特征。实验表明,诊断效果与故障尺寸、速度/负载无关,而与组件功能有关。倒谱预白化和基准方法(离散/随机分离+频谱峰度)在诊断中表现出较好的性能。
摘要由CSDN通过智能技术生成

Rolling element bearing diagnostics using the Case Western Reserve University data: A benchmark study.

  1. 论文通过将三种已建立的诊断技术应用于整个CWRU数据集来提供这样的基准。所有方法都使用平方包络频谱(即平方包络的频谱)作为最终诊断工具,但是在获取包络信号之前使用了不同的预处理步骤。
  2. 滚动轴承的局部故障会在加速度信号中产生一系列宽带脉冲响应,因为轴承组件会反复出现故障。轴承诊断的关键是,通过幅度解调获得的包络信号通常包含比原始信号更清晰的故障信息。
  3. 大多数轴承诊断技术的主要工具是包络谱,包络谱:对信号进行hilbert变换(信号处理中的一种常用手段,本质上是卷积)之后,然后取极值,然后对取极值之后得到的一维数据取包络,对包络信号进行FFT变换得到的数据。(横坐标为频率,纵坐标为幅值)包络谱对冲击事件的故障比较敏感。
  4. 轴承承受的唯一载荷(理论上)是轴和任何附加组件的重量产生的静态重力载荷,尽管有证据表明可能存在动态载荷叠加在此静态负载上。
  5. 轴承上的唯一径向负载(理论上)是通过6.00点钟位置(而不是3.00点钟位置)作用的静态重力负载。
  6. envelope analysis of the raw signal:仅包括全带宽原始信号的包络分析(包络平方频谱)。
  7. cepstrum prewhitening :1.倒谱预白化,将所有频率分量设置为相同的幅度;2.全带宽信号的包络分析(平方包络频谱)。
  8. benchmark method:1.离散/随机分离(DRS)删除确定性(离散频率)分量;2.频谱峰度以确定最冲动的频带,然后进行带通滤波;3.带通滤波后信号的包络分析(平方包络谱)。
  9. 论文中使用的DRS设置-滤波器长度N和延迟Δ(以样本数为单位)–是通过反复试验在少量数据集上建立的,其中12k数据选择了N = 16384,Δ= 500,N = 8192,Δ= 500用于48k数据。光谱峰度是使用 Antoni’s Fast Kurtogram 。
  10. 功率谱密度(power spectral density):单位频率间隔的光功率或噪声功率。

  1. 新算法可以对P类(P1:数据可能是可诊断的;例如,包络频谱显示了预期故障频率处的离散分量,但它们在频谱中并不占主导地位;P2:潜在可诊断的数据;例如,包络频谱显示拖尾的成分似乎与预期的故障频率一致。)中的数据集进行更全面的诊断或对N类(N1:无法针对指定的轴承故障诊断数据,但存在其他可识别的问题(例如松动);N2:数据不可诊断,并且几乎无法与噪声区分开,包络谱中的轴谐波可能例外。)中的数据集进行成功的诊断。
  2. 诊断结果似乎与故障尺寸或速度/负载无关,而与组件的功能有关,这对于每种故障尺寸大概都是相同的,但是当安装新轴承时,在故障尺寸之间会有所不同。怀疑这种现象可归因于机械松动,其严重性随每个轴承的安装而改变。
  3. 最清晰的数据集之一是数据集3007DE,尽管可以通过方法1进行诊断,但使用方法3(基准)进行处理时却具有非常经典的特性。最强的断球谐波是BSF的2倍和4倍,边带间隔为FTF(最强的二阶边带)。 FTF也至少存在三个谐波。
  4. 许多球故障数据集的另一个有趣特征是,经常有外部和内部种族故障的证据。
  5. 在(48k数据)频谱(直接使用fft,fft_size=1024)的11-14 kHz区域中通常很明显:轴速度的非常高的谐波。
  6. 功率谱是信号自相关函数的傅里叶变换,能量谱是信号本身傅立叶变换幅度的平方。

H(t)为Hilbert变换后的时域信号,f(t)为原始时域信号。那么其包络为:Envelop = sqrt(H^2(t)+f^2(t))。

???

???


注:

Method 1 .envelope analysis of the raw signal:仅包括全带宽原始信号的包络分析(包络平方频谱)。

Method 2.cepstrum prewhitening (!!!):

平方包络信号的仿真 帮助大家学者仿真平方包络信号 % squaredd1 % 调制函数呈衰减趋势 % 滤波法 % 每阶啮合谐波加入两阶调制谐波 clear all % 冲击形状按照文章中的表达式 N=5120; t=1:N; fs=1000; t=t/fs; u=0.925; deta=0.05; d0=exp(-(log(t)-u).^2./(deta.^2))./(t.*deta.*sqrt(2*pi)); a11=0.12.*exp(-0.15.*t).*cos(2*pi.*t); a12=0.08.*exp(-0.15.*t).*cos(2*pi.*t); a1=a11+a12; b11=0.12.*exp(-0.15.*t).*cos(2*pi.*t); b12=0.08.*exp(-0.15.*t).*cos(2*pi.*t); b1=b11+b12; d=1.45.*d0./max(abs(d0)); a1=a1+d/121; b1=b1-d/121; a21=0.12.*exp(-0.15.*t).*cos(4*pi.*t); a22=0.08.*exp(-0.15.*t).*cos(4*pi.*t); a2=a21+a22; b21=0.12.*exp(-0.15.*t).*cos(4*pi.*t); b22=0.08.*exp(-0.15.*t).*cos(4*pi.*t); b2=b21+b22; a2=a2+d/121; b2=b2-d/121; a31=0.12.*exp(-0.15.*t).*cos(6*pi.*t); a32=0.08.*exp(-0.15.*t).*cos(6*pi.*t); a3=a31+a32; b31=0.12.*exp(-0.15.*t).*cos(6*pi.*t); b32=0.08.*exp(-0.15.*t).*cos(6*pi.*t); b3=b31+b32; a41=0.12.*exp(-0.15.*t).*cos(8*pi.*t); a42=0.08.*exp(-0.15.*t).*cos(8*pi.*t); a4=a41+a42; b41=0.12.*exp(-0.15.*t).*cos(8*pi.*t); b42=0.08.*exp(-0.15.*t).*cos(8*pi.*t); b4=b41+b42; a51=0.12.*exp(-0.15.*t).*cos(10*pi.*t); a52=0.08.*exp(-0.15.*t).*cos(10*pi.*t); a5=a51+a52; b51=0.12.*exp(-0.15.*t).*cos(10*pi.*t); b52=0.08.*exp(-0.15.*t).*cos(10*pi.*t); b5=b51+b52; r1=10*(a1.*cos(54*pi*t)-b1.*sin(54*pi*t)); r2=8*(a2.*cos(108*pi*t)-b2.*sin(108*pi*t)); r3=6*(a3.*cos(162*pi*t)-b3.*sin(162*pi*t)); r4=5*(a4.*cos(216*pi*t)-b4.*sin(216*pi*t)); r5=2*(a5.*cos(270*pi*t)-b5.*sin(270*pi*t)); yr=r1+r2+r3+r4+r5+d.*cos(242*pi*t); Y=fft(yr); fil=[zeros(1,610) Y(611:624) zeros(1,4496)];%一共5120个点 fy1=ifft(fil); n1=10;wn1=[110 135]/(fs/2); % 带通滤波 [b,a]=butter(n1,wn1); fy1=filter(b,a,yr); y=fy1.*fy1; Y=fft(y); fil=[Y(1:50) zeros(1,5070)];%一共5120个点 y=abs(ifft(fil)); t=t*360*1000/5120; plot(t,y); xlabel('轴转角/度'); ylabel('Am/mm'); title('平方包络信号的时域波形');
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值
>