CN109883706B - 一种滚动轴承局部损伤微弱故障特征提取方法 - Google Patents
一种滚动轴承局部损伤微弱故障特征提取方法 Download PDFInfo
- Publication number
- CN109883706B CN109883706B CN201910276672.8A CN201910276672A CN109883706B CN 109883706 B CN109883706 B CN 109883706B CN 201910276672 A CN201910276672 A CN 201910276672A CN 109883706 B CN109883706 B CN 109883706B
- Authority
- CN
- China
- Prior art keywords
- rolling bearing
- node
- signal
- fault
- 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.)
- Active
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种滚动轴承局部损伤微弱故障特征提取方法,使用加速度传感器采集故障滚动轴承的振动信号,使用MODWPT分解原信号至不同的节点,接下来计算各节点均方包络自相关的SK,然后选取各层中不小于最大SK一半的节点信号进行MED滤波,再依据层数叠加同层所选节点的归一化频谱信息并平均,最后合并所有层的频谱,提取滚动轴承的局部损伤故障特征;本发明给出了一种强背景噪声干扰下提取滚动轴承局部损伤故障特征的方法,为实现滚动轴承的PHM早期微弱故障特征提取提供了一种有效的诊断工具。
Description
技术领域
本发明涉及故障诊断与信号处理分析技术领域,特别涉及一种滚动轴承局部损伤微弱故障特征提取方法。
背景技术
滚动轴承作为旋转机械的核心零部件之一,其运行状态好坏与机械设备的工作状态正常与否密切相关。存在局部损伤的故障滚动轴承工作时,会生成的冲击信号,根据损伤位置不同故障特征频率也不同可完成滚动轴承的故障诊断。在故障特征提取时,傅里叶变换是一种有效的故障特征提取方法,但受到环境噪声等强干扰的影响,原始信号中的微弱故障特征极易被淹没而难以提取,有必要结合相关的信号处理方法进行微弱故障特征提取。
合并均方包络谱方法是一种有效的故障特征提取方法,其使用最大迭代离散小波包变换(Maximal Overlap Discrete Wavelet Packet Transform,MODWPT)分解信号。MODWPT是一种特殊的小波变换,其各分解层所得信号的时间分辨率与原信号的分辨率相同,没有相位扭曲,且小波系数与尺度系数具有平移不变性。因此,原信号经MODWPT分解后不存在因降采样操作导致的时间长度减半问题,这对于小波分解后应用各小波系数的性能十分重要。在选择有用节点信号时,其通过对瞬态冲击特征明显的峭度进行选择,而峭度指标并没有考虑其频率分布的情况,对于后期通过频谱合并的方法进行特征故障提取是不利的。且需要指出的是,选定的各频带信号内同样存在噪声干扰,不益于微弱故障特征的提取。MED(Minimum Entropy Deconvolution)最早用于地震信号的处理,之后被引入故障诊断领域,对振动信号中的故障特征进行增强。MED通过迭代的方法解卷积振动信号,当达到规定的终止条件后,根据最大峭度准则输出MED最优滤波信号,用于滚动轴承的降噪分析。但滚动轴承的早期故障较微弱,对原信号直接使用MED降噪易受瞬态冲击等干扰信号的影响,难以有效提取故障特征信号。
发明内容
为了克服上述现有技术的缺陷,本发明的目的在于提供了一种滚动轴承局部损伤微弱故障特征提取方法,通过使用MODWPT分解信号,计算各节点均方包络自相关的SK,然后选取各层不小于最大SK的节点信号进行MED滤波,最后将所得信号的均方包络谱归一化并联合分析得到故障特征频谱,解决了滚动轴承早期微弱故障受环境噪声干扰难以提取的问题。
为了达到上述目的,本发明的技术方案是:
一种滚动轴承局部损伤微弱故障特征提取方法,具体步骤如下:
Step1、通过加速度传感器获取加速度振动信号x(t);
Step2、将加速度振动信号x(t)使用MODWPT方法分解,得到不同中心频率fc,不同滤波频带bw的节点信号式中k表示第k层分解层数,k=0,1,…,K,K为确定的最大分解层数,i表示第k层分解的第i个节点信号,i=1,2,…,2k,各节点信号具有相同的长度,保证了一致的时间分辨率;
Step3、对步骤Step2所得分解信号取平方,并根据Hilbert变换获取包络,接下来计算均方包络无偏自相关函数的SK(Spectrum Kurtosis),去除各层中SK值小于该层最大SK值一半的所有节点信号,其它节点信号留作进一步分析;
Step4、使用MED方法对各层保留的节点信号进行滤波,得到故障特征增强的节点信号m为第k层选出的第m个节点信号;
Step5、计算步骤Step4中得到的各节点信号的均方包络,再使用傅里叶变换方法转换至频域,然后对各频谱进行归一化,然后依所属分解层数叠加上述频谱并取平均,得到与分解层数相同数量的频谱;
Step6、将上一步骤所得频谱进行叠加,得出最终的频谱信号,找出对应的故障特征频率,诊断滚动轴承的故障。
本发明所提出的方法通过对滚动轴承局部损伤微弱故障特征的提取,为实现强背景噪声干扰下的滚动轴承故障诊断给出了一种较为有效的方法,为实现滚动轴承PHM(Prognostic and Health Management)中提取早期局部损伤微弱故障特征提供了一种有效的诊断工具。
本发明的有益效果是:通过计算MODWPT分解后各节点均方包络谱自相关的SK,可以反映峭度分布随频率的变化,利用MED对各所选节点信号进行滤波,增强频带内的重复性故障冲击特征,通过叠加各层归一化频谱信息并平均,然后合并所有层的频谱,提取滚动轴承的局部损伤故障特征。采用本发明益于提取强干扰噪声下的微弱故障特征,适用于滚动轴承的早期故障诊断。
附图说明
图1是本发明方法步骤流程图。
图2是本发明实施例中原始振动信号。
图3是本发明图2经本发明方法得到的频谱。
图4是本发明图2经包络分析方法得到的频谱。
具体实施方式
下面结合附图和具体实施方式对本发明做详细叙述。
一种滚动轴承局部损伤微弱故障特征提取方法,具体步骤如下:
具体参数如下:1)选用美国凯斯西储大学的滚动轴承故障数据进行分析,滚动轴承为SKF的6205型,接触角为0度,节径大小为39mm,滚动体直径为7.94mm,滚动体个数为9个;2)故障类型为外滚道单点损伤,直径为0.533mm;3)采集实施例中的采样频率为12KHz,驱动电机的转速为1797r/min。
Step1、通过加速度传感器获取加速度振动信号x(t);由于存在一定的机械松动,该信号受到强烈的调制干扰,原始信号如图2所示;
Step2、将加速度振动信号x(t)使用MODWPT方法分解,得到不同中心频率fc,不同滤波频带bw的节点信号式中k表示第k层分解层数,k=0,1,…,K,设定分解层数K为5,i表示第k层分解的第i个节点信号,i=1,2,…,2k,各节点信号具有相同的长度,保证了一致的时间分辨率;
Step3、对步骤Step2所得分解信号取平方,并根据Hilbert变换获取包络,接下来根据公式(1)计算均方包络无偏自相关函数的SK(Spectrum Kurtosis),去除各层中SK值小于该层最大SK值一半的所有节点信号,留作进一步分析的其它节点信号位置如表1所示;
式中,[·]表示数学期望,X(tn,f)为频率f处信号x(tn)的复包络信号。
表1节点信号位置
Step4、使用MED方法对各层保留的节点信号进行滤波,得到故障特征增强的节点信号m为第k层选出的第m个节点信号;
Step5、计算步骤Step4中得到的各信号的均方包络,再使用傅里叶变换方法转换至频域,然后对各频谱进行归一化,然后依所属分解层数叠加上述频谱并取平均,得到与分解层数相同数量的频谱;
Step6、将上一步骤所得频谱进行叠加,得出最终的频谱,如图3所示。图中,点划线对应上述滚动轴承的故障特征频率162.186Hz及其2阶至5阶倍频成分,可以看到,方框标注的故障特征频率比较明显,完成了滚动轴承的微弱故障诊断。为对比分析,使用当前应用较为广泛的包络分析方法处理同一组信号,所得结果如图4所示。可以看到,方框标注的故障特征频率受干扰严重,不能有效提取,且倍频成分并不明显,故障诊断失效。因此,新方法能够有效提取滚动轴承微弱故障特征。
上面结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。
Claims (1)
1.一种滚动轴承局部损伤微弱故障特征提取方法,其特征在于,具体步骤如下:
Step1、通过加速度传感器获取加速度振动信号x(t);
Step2、将加速度振动信号x(t)使用MODWPT方法分解,得到不同中心频率fc,不同滤波频带bw的节点信号式中k表示第k层分解层数,k=0,1,…,K,K为确定的最大分解层数,i表示第k层分解的第i个节点信号,i=1,2,…,2k,各节点信号具有相同的长度,保证了一致的时间分辨率;
Step3、对步骤Step2所得分解信号取平方,并根据Hilbert变换获取包络,接下来计算均方包络无偏自相关函数的SK(Spectrum Kurtosis),去除各层中SK值小于该层最大SK值一半的所有节点信号,其它节点信号留作进一步分析;
Step4、使用MED方法对各层保留的节点信号进行滤波,得到故障特征增强的节点信号m为第k层选出的第m个节点信号;
Step5、计算步骤Step4中得到的各节点信号的均方包络,再使用傅里叶变换方法转换至频域,然后对各频谱进行归一化,然后依所属分解层数叠加上述频谱并取平均,得到与分解层数相同数量的频谱;
Step6、将上一步骤所得频谱进行叠加,得出最终的频谱信号,找出对应的故障特征频率,诊断滚动轴承的故障。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910276672.8A CN109883706B (zh) | 2019-04-08 | 2019-04-08 | 一种滚动轴承局部损伤微弱故障特征提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910276672.8A CN109883706B (zh) | 2019-04-08 | 2019-04-08 | 一种滚动轴承局部损伤微弱故障特征提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109883706A CN109883706A (zh) | 2019-06-14 |
CN109883706B true CN109883706B (zh) | 2019-12-31 |
Family
ID=66936384
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910276672.8A Active CN109883706B (zh) | 2019-04-08 | 2019-04-08 | 一种滚动轴承局部损伤微弱故障特征提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109883706B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110455530B (zh) * | 2019-09-18 | 2021-08-31 | 福州大学 | 谱峭度结合卷积神经网络的风机齿轮箱复合故障诊断方法 |
CN111024347B (zh) * | 2019-12-26 | 2021-11-02 | 华中科技大学 | 一种电抗器健康状态评估方法 |
CN111238813B (zh) * | 2020-01-19 | 2021-05-07 | 西安交通大学 | 一种强干扰下提取滚动轴承故障特征的方法 |
CN112577722B (zh) * | 2020-11-19 | 2023-06-23 | 江苏科技大学 | 基于平方包络和零频率谐振器的弱故障诊断方法 |
CN112857798B (zh) * | 2021-02-23 | 2023-03-10 | 太原理工大学 | 一种基于频谱分析的多轴系机械传动系统质量评估方法及装置 |
CN113702042B (zh) * | 2021-07-14 | 2022-10-28 | 西南交通大学 | 基于一种增强最小熵解卷积的机械故障诊断方法及系统 |
CN114441174A (zh) * | 2022-02-09 | 2022-05-06 | 上海电气集团股份有限公司 | 滚动轴承复合故障的诊断方法、系统、设备及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104655425A (zh) * | 2015-03-06 | 2015-05-27 | 重庆大学 | 基于稀疏表示和大间隔分布学习的轴承故障分类诊断方法 |
CN106289774A (zh) * | 2016-07-26 | 2017-01-04 | 北京工业大学 | 一种滚动轴承故障识别与趋势预测方法 |
CN107631866A (zh) * | 2017-08-31 | 2018-01-26 | 西安交通大学 | 一种低速运行旋转机械的早期故障识别方法 |
CN107631877A (zh) * | 2017-08-11 | 2018-01-26 | 南京航空航天大学 | 一种针对机匣振动信号的滚动轴承故障协同诊断方法 |
CN108168891A (zh) * | 2018-02-26 | 2018-06-15 | 成都昊铭科技有限公司 | 滚动轴承微弱故障信号特征的提取方法及设备 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10677685B2 (en) * | 2013-04-05 | 2020-06-09 | Aktiebolaget Skf | Method, computer program product and system |
-
2019
- 2019-04-08 CN CN201910276672.8A patent/CN109883706B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104655425A (zh) * | 2015-03-06 | 2015-05-27 | 重庆大学 | 基于稀疏表示和大间隔分布学习的轴承故障分类诊断方法 |
CN106289774A (zh) * | 2016-07-26 | 2017-01-04 | 北京工业大学 | 一种滚动轴承故障识别与趋势预测方法 |
CN107631877A (zh) * | 2017-08-11 | 2018-01-26 | 南京航空航天大学 | 一种针对机匣振动信号的滚动轴承故障协同诊断方法 |
CN107631866A (zh) * | 2017-08-31 | 2018-01-26 | 西安交通大学 | 一种低速运行旋转机械的早期故障识别方法 |
CN108168891A (zh) * | 2018-02-26 | 2018-06-15 | 成都昊铭科技有限公司 | 滚动轴承微弱故障信号特征的提取方法及设备 |
Non-Patent Citations (1)
Title |
---|
基于MODWPT的包络阶次谱在滚动轴承故障诊断中的应用;杨宇等;《哈尔滨工程大学学报》;20101031;第31卷(第10期);第1380-1385页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109883706A (zh) | 2019-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109883706B (zh) | 一种滚动轴承局部损伤微弱故障特征提取方法 | |
CN107505135B (zh) | 一种滚动轴承复合故障提取方法及系统 | |
CN111089726B (zh) | 一种基于最优维数奇异谱分解的滚动轴承故障诊断方法 | |
CN108444709B (zh) | VMD与FastICA相结合的滚动轴承故障诊断方法 | |
CN106404399B (zh) | 基于自适应冗余提升小波包分解树的轴承故障诊断方法 | |
CN108388908B (zh) | 基于k-svd和滑窗降噪的滚动轴承冲击故障诊断方法 | |
CN111238813B (zh) | 一种强干扰下提取滚动轴承故障特征的方法 | |
CN112101245A (zh) | 基于频域窗函数的短时傅里叶变换机械冲击特征提取方法 | |
CN110046476B (zh) | 滚动轴承故障的三元二进分形小波稀疏诊断方法 | |
CN110672327A (zh) | 一种基于多层降噪技术的异步电机轴承故障诊断方法 | |
CN108398260B (zh) | 基于混合概率方法的齿轮箱瞬时角速度的快速评估方法 | |
CN113607415A (zh) | 一种变转速下基于短时随机共振的轴承故障诊断方法 | |
CN111769810A (zh) | 一种基于能量峭度谱的流体机械调制频率提取方法 | |
CN114486263B (zh) | 一种旋转机械滚动轴承振动信号降噪解调方法 | |
CN106908232A (zh) | 一种旋转机械的振动信号故障识别方法 | |
CN112926014A (zh) | 基于rls与rssd的滚动轴承声信号多频带融合故障诊断方法 | |
CN117571316A (zh) | 一种复合故障诊断方法及系统 | |
CN111323233B (zh) | 一种用于低速旋转机械故障诊断的局部均值分解方法 | |
CN107490477B (zh) | 基于频谱核密度函数相关性比较的齿轮箱故障诊断方法 | |
CN116150585A (zh) | 基于乘积包络谱的旋转机械故障诊断方法 | |
CN116361733A (zh) | 一种故障诊断方法、装置、系统以及存储介质 | |
CN114001961B (zh) | 一种滚动轴承故障特征靶向定位及提取方法 | |
CN112345247A (zh) | 一种滚动轴承的故障诊断方法及装置 | |
CN110595751B (zh) | 基尼指数引导的早期故障特征小波重构方法及其应用 | |
CN114923689A (zh) | 一种基于局部特征尺度分解的滚动轴承故障诊断方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |