CN102680080B - 一种基于改进的自适应形态滤波的非稳态信号检测方法 - Google Patents

一种基于改进的自适应形态滤波的非稳态信号检测方法 Download PDF

Info

Publication number
CN102680080B
CN102680080B CN201210135037.6A CN201210135037A CN102680080B CN 102680080 B CN102680080 B CN 102680080B CN 201210135037 A CN201210135037 A CN 201210135037A CN 102680080 B CN102680080 B CN 102680080B
Authority
CN
China
Prior art keywords
signal
length
operator
area
signals
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201210135037.6A
Other languages
English (en)
Other versions
CN102680080A (zh
Inventor
沈长青
孔凡让
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN201210135037.6A priority Critical patent/CN102680080B/zh
Publication of CN102680080A publication Critical patent/CN102680080A/zh
Application granted granted Critical
Publication of CN102680080B publication Critical patent/CN102680080B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

一种基于改进的自适应形态滤波的非稳态信号检测方法,用于各类旋转机械设备故障信号特征提取,步骤为:(1)绘制原始信号时域图,从原始信号中寻找并标定出所有信号局部极大值,计算由局部极大值与时域图中时间轴围成的轮廓面积,定为参考面积;(2)利用不同长度的结构元素和闭-开形态算子对原始信号进行形态滤波,并计算信号经不同结构元素形态滤波后得到信号的局部极值与时间轴围成的轮廓面积,并与步骤(1)中得到的参考面积求差,确定最小偏差对应的结构元素长度;(3)利用步骤(2)中得到的结构元素对信号进行形态滤波,分析处理后信号的频谱,判断故障。本发明提高了旋转设备故障判断的效率和准确性。

Description

一种基于改进的自适应形态滤波的非稳态信号检测方法
技术领域
本发明涉及一种非稳态信号的分析检测方法,具体涉及一种自适应形态滤波的非稳态信号检测方法,用于对旋转机械振动信号中反映故障特征的非稳态信号进行检测。
背景技术
对于反映机械设备故障的信号中非稳态成分的检测,在机械设备的故障诊断、生物医学信号的检测等领域具有广泛的应用。常见的非稳态成分检测方法是在时域信号中观察是否存在非稳态成分,但是由于检测过程中不可避免的会混入噪声,其中表征故障的非稳态成分亦会被噪声污染,直接观察的方式准确性和效率都很低。另外一种常用的方法是通过频域分析信号中的周期特征,但是对于信号中持续时间较短的非稳态信号,在频谱中表现为较小的幅值,甚至同样被噪声淹没,因此通过频域分析检测往往不能得到显著的特征。此外,对信号进行滤波分析也是去除信号中噪声干扰,从而提取非稳态成分的一种有效方法,但是其效果受设计的滤波器的关键参数,如截止频率、带宽、中心频率等因素影响很大。随着研究的深入,在实践中,常用的检测方法是对信号进行小波包层层分解,然后对小波节点信号进行解调,分析频域特征,进而得到信号中非平稳成分的周期,但这些技术存在需要先验知识或者耗时长的缺点。因此,总体来说现有的检测方法存在对旋转设备故障判断的效率低和准确性不高的缺点。
发明内容
本发明的技术解决问题:克服现有技术的不足,提供一种基于改进的自适应形态滤波的非稳态信号检测方法,提高了旋转设备故障判断的效率和准确性。
本发明技术解决方案:一种基于改进的自适应形态滤波的非稳态信号检测方法。在待检测设备的壳体上安装加速度传感器,检测设备的振动加速度信号,作为检测信号。步骤如下:
(1)绘制原始信号的时域图,找出信号各局部极大值的位置,并计算其在时域图中与时间轴围成的轮廓区域面积为参考面积:
S ref = Σ i = 1 k - 1 S local _ i
式中,Sref为原始信号极值与时间轴围成的轮廓区域面积,Slocal_i=0.5×lengthi×(maxi+maxi+1)为相邻两个局部极值与时间轴围成的轮廓区域面积,k为局部极值个数,i=1,2,...k-1,lengthi为局部极值间的距离,maxi为第i个局部极值。该步骤以信号局部极值与时间轴围成的轮廓区域面积为形态滤波结构元素设计参考标准,力求通过保留极值信号来实现移除尽可能多的噪声信号,保留尽可能多故障信号;
(2)设定结构元素基本步长step,step为信号采样周期,即采样频率倒数,依次增加扫描各倍率长度,从1×step,2×step,…,m×step(m为采样周期倍数,取值范围为1-30)分别对原始信号进行形态滤波结果与时间轴围成的轮廓区域面积,并与步骤(1)中所得参考面积对比,找出两者面积最接近情况下采用的结构元素的步长,定为最佳结构元素;
(3)根据步骤(2)中得到的最佳结构元素长度,根据下述计算公式利用闭-开算子对信号进行形态滤波并分析频谱:
上式中,Θ,ο和·分别表示腐蚀算子、膨胀算子、开算子和闭算子,其中,基本腐蚀算子、膨胀算子、开算子、闭算子分别表达如下:
腐蚀算子:
(fΘg)(n)=min[f(n+m)-g(m)]
膨胀算子:
( f ⊕ g ) ( n ) = max [ f ( n - m ) + g ( m ) ]
开算子:
闭算子:
( f · g ) ( n ) = ( f ⊕ gΘg ) ( n )
f(n)为待处理信号,长度为N,g(m)为结构元素,长度为M,m=0,1,2...M-1,n=0,1,2...N-1。
本发明与现有技术相比的优点在于:
(1)本发明所提供的非稳态信号检测方法在待检测设备的壳体上安装加速度传感器,检测设备的振动加速度信号,作为检测信号x(t),通过权利要求1所述步骤,即改进的自适应形态滤波器可以自行确定合适的形态滤波器结构元素长度,从而进行滤波得到滤波信号,然后再计算滤波所得的信号的频谱,即可得出故障相关的非稳态信号频率进而判断该非稳态信号频率所对应的故障,本发明能根据信号自身的特点设计形态滤波器,从而能更好地获取非稳态信号,该检测标准具有更高的准确性。
(2)根据权利要求1步骤(1)-(2),本发明在结构元素设计上,完全由信号自身的特性决定,不需要先验知识,因而具有自适应的优点。
(3)根据权利要求1步骤(3),本发明通过对机械设备的振动信号的非稳态成分对应频谱分析,可以实现对机械设备的故障的自动识别。
附图说明
图1为本发明方法实现流程图;
图2为本发明的仿真信号及其频谱;
图3为本发明的改进形态滤波方法处理后信号及其频谱;
图4为本发明的外圈故障信号及其频谱;
图5为本发明中形态滤波算法用于外圈故障信号处理结果及其频谱。
具体实施方式
实施例1:
仿真信号的时域波形如图2(a)所示,采样频率为12KHz,即采样间隔为1/12000s。为验证其降噪性能,对仿真信号叠加幅值为0.5的白噪声,时域波形如图2(b)所示。图2(c)为其频谱,可以看出:不能从中看出有效的频率成分。
利用提出的改进的形态滤波方法对信号进行处理,具体步骤为:
1、根据发明内容的步骤(1),寻找图2(b)所示时域信号极值,并计算其与时间轴围成的面积。
2、根据发明内容的步骤(2),以采样间隔为基准长度,成倍增加结构元素长度,其滤波后信号与横轴围成的面积与参考面积之差随结构元素长度变化如图3(a)所示。可见,随结构元素长度的增加,滤波后信号与横轴围成的面积与参考面积之差呈现先变小,后变大的趋势。当结构元素长度选择为7个采样间隔对应时间长度,即7/12000s时,面积差最小,故选择为最佳结构元素长度。
3、根据发明内容的步骤(3),最佳结构元素处理后的滤波信号如图3(b)所示,滤波结果能很好的保留脉冲信号,噪声信号也得到有效抑制,图3(c)绘制了滤波结果的频谱,可以看出,信号特征频率100Hz及其二次、三次谐波频率成分幅度均较为清晰。从对含有噪声的仿真信号分析结果来看,分析结果与信号实际特征频率吻合,自适应的形态滤波方法能有效的提取脉冲成分,抑制噪声成分,有助于提取信号中周期瞬态成分的时间间隔。
实施例2:
采用实际的轴承数据进行处理。轴承型号为6205-2RS JEM SKF。参数如表1所示。
表16205-2RS JEM SKF轴承参数(单位:英寸)
信号采样频率为12KHz,轴承在1797rpm转速时外圈故障特征频率为107.3Hz。图4(a)为该轴承在1797rpm转速下外圈故障信号,可以从中看出故障脉冲但是无法定性判断故障类别,图4(b)为该信号对应频谱,从信号频谱图中亦无法识别故障相关的频率成分。利用本发明提出的形态滤波方法对该信号处理,具体步骤为:
1、根据发明内容的步骤(1),寻找图4(a)所示时域信号极值,并计算其与时间轴围成的面积。
2、根据发明内容的步骤(2),以采样间隔为基准长度,成倍增加结构元素长度,其滤波后信号与横轴围成的面积与参考面积之差随结构元素长度变化如图5(a)所示。可见,随结构元素长度的增加,滤波后信号与横轴围成的面积与参考面积之差呈现先变小,后变大的趋势。当结构元素长度选择为12个采样间隔对应时间长度,即1/1000s时,面积差最小,故选择为最佳结构元素长度。
3、根据发明内容的步骤(3),最佳结构元素处理后的滤波信号如图5(b)所示,滤波结果能很好的保留脉冲信号,噪声信号也得到有效抑制,图5(c)绘制了滤波结果的频谱,可以看出,外圈故障频率及其二次、三次谐波频率成分幅度均较为清晰。分析结果与实际故障频率吻合,从对外圈故障信号分析结果来看,自适应的形态滤波方法能有效的提取脉冲成分,抑制噪声成分,有助于提取信号中周期瞬态成分的时间间隔。
本发明未详细阐述部分属于本领域技术人员的公知技术。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。

Claims (1)

1.一种基于改进的自适应形态滤波的非稳态信号检测方法,在待检测设备的壳体上安装加速度传感器,检测设备的振动加速度信号,作为检测信号,其特征在于实现步骤如下:
(1)绘制原始信号的时域图,找出信号各局部极大值的位置,并计算信号局部极大值在时域图中与时间轴围成的轮廓区域面积为参考面积:
S ref = Σ i = 1 k - 1 S local _ i
式中,Sref为原始信号极大值与时间轴围成的轮廓区域面积,Slocal_i=0.5×lengthi×(maxi+maxi+1)为相邻两个局部极大值与时间轴围成的轮廓区域面积,k为局部极大值个数,i=1,2,...k-1,lengthi为局部极大值间的距离,maxi为第i个局部极大值;
(2)设定结构元素基本步长step,step为信号采样周期,即采样频率倒数,依次增加扫描各倍率长度,从1×step,2×step,…,z×step分别对原始信号进行形态滤波结果与时间轴围成的轮廓区域面积,并与步骤(1)中所得参考面积对比,找出两者面积最接近情况下采用的结构元素的步长,定为最佳结构元素;z为采样周期倍数;
(3)根据步骤(2)中得到的最佳结构元素长度,根据下述计算公式利用闭-开算子对信号进行滤波,并分析处理后信号的频谱,判断故障类别:
上式中,Θ,ο和·分别表示腐蚀算子、膨胀算子、开算子和闭算子,其中,基本腐蚀算子、膨胀算子、开算子、闭算子分别表达如下:
腐蚀算子:
( fΘg ) ( n ) = min [ f ( n + m ) - g ( m ) ]
膨胀算子:
( f ⊕ g ) ( n ) = max [ f ( n - m ) + g ( m ) ]
开算子:
闭算子:
( f · g ) ( n ) = ( f ⊕ gΘg ) ( n )
f(n)为待处理信号,长度为N,g(m)为结构元素,长度为M,m=0,1,2...M-1,n=0,1,2...N-1。
CN201210135037.6A 2012-05-03 2012-05-03 一种基于改进的自适应形态滤波的非稳态信号检测方法 Expired - Fee Related CN102680080B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210135037.6A CN102680080B (zh) 2012-05-03 2012-05-03 一种基于改进的自适应形态滤波的非稳态信号检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210135037.6A CN102680080B (zh) 2012-05-03 2012-05-03 一种基于改进的自适应形态滤波的非稳态信号检测方法

Publications (2)

Publication Number Publication Date
CN102680080A CN102680080A (zh) 2012-09-19
CN102680080B true CN102680080B (zh) 2014-09-10

Family

ID=46812352

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210135037.6A Expired - Fee Related CN102680080B (zh) 2012-05-03 2012-05-03 一种基于改进的自适应形态滤波的非稳态信号检测方法

Country Status (1)

Country Link
CN (1) CN102680080B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102988041B (zh) * 2012-11-16 2018-04-06 中国科学院上海微系统与信息技术研究所 心磁信号噪声抑制中的信号选择性平均方法
CN105207645B (zh) * 2015-08-25 2018-09-07 苏州汇川技术有限公司 基于数学形态学的滤波方法以及滤波系统
CN108760037B (zh) * 2018-06-15 2020-03-17 西安交通大学 一种基于频谱分析的风力发电机叶片结构损伤检测方法
CN110470476A (zh) * 2019-09-11 2019-11-19 天津浪腾科技有限公司 一种用于滚动轴承故障诊断的脉冲提取方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0969283A1 (en) * 1998-06-25 2000-01-05 Hewlett-Packard Company A method for processing measuring values
CN101435799A (zh) * 2008-12-19 2009-05-20 清华大学 基于声发射技术的水轮机故障诊断方法及装置
CN101644623A (zh) * 2009-06-19 2010-02-10 湖南大学 基于多尺度形态学分析的齿轮故障诊断方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6863791B1 (en) * 2000-09-11 2005-03-08 Spectrumedix Llc Method for in-situ calibration of electrophoretic analysis systems

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0969283A1 (en) * 1998-06-25 2000-01-05 Hewlett-Packard Company A method for processing measuring values
CN101435799A (zh) * 2008-12-19 2009-05-20 清华大学 基于声发射技术的水轮机故障诊断方法及装置
CN101644623A (zh) * 2009-06-19 2010-02-10 湖南大学 基于多尺度形态学分析的齿轮故障诊断方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于ASEGMF的旋转机械振动信号降噪方法研究;张文斌等;《振动与冲击》;20111231;第30卷(第9期);26-29 *
张文斌等.基于ASEGMF的旋转机械振动信号降噪方法研究.《振动与冲击》.2011,第30卷(第9期),26-29.

Also Published As

Publication number Publication date
CN102680080A (zh) 2012-09-19

Similar Documents

Publication Publication Date Title
Yan et al. Fault diagnosis of rolling element bearing using a new optimal scale morphology analysis method
CN108168891B (zh) 滚动轴承微弱故障信号特征的提取方法及设备
Zheng et al. Incipient fault detection of rolling bearing using maximum autocorrelation impulse harmonic to noise deconvolution and parameter optimized fast EEMD
CN109668733B (zh) 变分非线性模式分解变转速轴承故障诊断方法
Shi et al. Bearing fault diagnosis under variable rotational speed via the joint application of windowed fractal dimension transform and generalized demodulation: A method free from prefiltering and resampling
Lei et al. Fault diagnosis of rotating machinery using an improved HHT based on EEMD and sensitive IMFs
Li et al. Research on test bench bearing fault diagnosis of improved EEMD based on improved adaptive resonance technology
CN109855874B (zh) 一种声音辅助振动微弱信号增强检测的随机共振滤波器
CN112101174A (zh) 一种基于LOF-Kurtogram的机械故障诊断方法
CN109883703B (zh) 一种基于振动信号相干倒谱分析的风机轴承健康监测诊断方法
Yu et al. A new morphological filter for fault feature extraction of vibration signals
Yan et al. A feature selection framework-based multiscale morphological analysis algorithm for fault diagnosis of rolling element bearing
He et al. A joint adaptive wavelet filter and morphological signal processing method for weak mechanical impulse extraction
CN102680080B (zh) 一种基于改进的自适应形态滤波的非稳态信号检测方法
Sheng et al. Applications in bearing fault diagnosis of an improved Kurtogram algorithm based on flexible frequency slice wavelet transform filter bank
CN110160765A (zh) 一种基于声音或振动信号的冲击特征识别方法及系统
Liu et al. An online bearing fault diagnosis technique via improved demodulation spectrum analysis under variable speed conditions
Williams et al. Helicopter transmission fault detection via time-frequency, scale and spectral methods
Zhang et al. Improved local cepstrum and its applications for gearbox and rolling bearing fault detection
Berntsen et al. Enhanced demodulation band selection based on Operational Modal Analysis (OMA) for bearing diagnostics
Bastami et al. Estimating the size of naturally generated defects in the outer ring and roller of a tapered roller bearing based on autoregressive model combined with envelope analysis and discrete wavelet transform
CN110231117A (zh) 一种基于s变换的拉索基频特征辨识方法
Du et al. Optimum IMFs selection based envelope analysis of bearing fault diagnosis in plunger pump
Wang et al. Information interval spectrum: a novel methodology for rolling-element bearing diagnosis
CN116361733A (zh) 一种故障诊断方法、装置、系统以及存储介质

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140910

Termination date: 20200503