CN107944199B - 一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法 - Google Patents

一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法 Download PDF

Info

Publication number
CN107944199B
CN107944199B CN201711408267.4A CN201711408267A CN107944199B CN 107944199 B CN107944199 B CN 107944199B CN 201711408267 A CN201711408267 A CN 201711408267A CN 107944199 B CN107944199 B CN 107944199B
Authority
CN
China
Prior art keywords
modal
vibration signal
frequency
frequency spectrum
trend
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
Application number
CN201711408267.4A
Other languages
English (en)
Other versions
CN107944199A (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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201711408267.4A priority Critical patent/CN107944199B/zh
Publication of CN107944199A publication Critical patent/CN107944199A/zh
Application granted granted Critical
Publication of CN107944199B publication Critical patent/CN107944199B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/021Gearings

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Control Of Transmission Device (AREA)

Abstract

本发明公开了一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法,它包括如下步骤:先采集目标齿轮箱的振动信号;再采用基于频谱趋势的方式对变分模态分解方法进行改进,利用改进的变分模态分解方法对齿轮箱的振动信号进行分解,得到振动信号的各个模态分量;接着对获得的模态分量进行解调分析,获得各个模态信号的特征频率;最后将特征频率与齿轮箱的各个齿轮转频进行对比,定位齿轮箱的故障源。本发明将时域范围的趋势概念引入频域范围内,提出了基于频谱趋势的改进变分模态分解的方法,该方法很好的解决了变分模态分解中模态数量需要预设的问题,利用频谱趋势极大值的归一化横坐标来初始化中心频率,使得分解的各个模态更为合理。

Description

一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法
技术领域
本发明属于齿轮箱故障识别技术领域,具体涉及一种基于频谱趋势和变分模态分解(VMD)的齿轮箱故障识别方法,它是一种基于频谱趋势(Spectrum Trend)的多分量个数估计并以频谱趋势为基础来改进变分模态分解(Variational Mode Decomposition,VMD)算法,用于多分量信号模态的提取。
背景技术
齿轮箱振动信号具有复杂多分量和调幅调频(AM-FM)的特点,幅值解调和频率解调方法能够避免传统Fourier频谱中的复杂边带分析,有效识别故障特征频率。但是对于多分量AM-FM信号有效的解调分析的一个前提就是提取有效的模态分量,小波变换(WaveletTransform,WT)和经验模态分解(Empirical Mode Decomposition,EMD)都是运用比较广泛的模态分量提取方法。但是WT存在小波基的选择和小波基确定之后缺乏自适应性的问题,而EMD是一种经验性的方法,缺乏完备的理论基础。
Dragomiretski等人提出VMD多分量信号分解方法,该方法采用频域非递归的迭代求解方式,将信号分解转换成以模态分量的估计带宽之和最小为目标的变分方式。该方法将维纳滤波器内嵌在算法之中,具有较好的鲁棒性。
VMD方法在直接处理齿轮箱振动信号时,模态分量的个数需要预设,分量个数估计的过大或者过小都会造成振动信号分解的不合理,影响故障识别的精度。此外不同的中心频率初始化方式也会造成分解得到的模态分量在频谱上不同的分布,也会对故障识别造成不利影响。这些问题也影响了VMD方法在齿轮箱故障识别领域中的广泛运用。
发明内容
针对现有技术中存在的上述问题,本发明专利的目的在于提供一种简便的、准确度高的基于频谱趋势和变分模态分解的齿轮箱故障识别方法,本发明的基于频谱趋势(Spectrum Trend)的改进VMD方法,命名为ST-VMD,并将其运用于齿轮箱故障诊断领域。
所述的一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法,其特征在于包括如下步骤:
步骤1采集目标齿轮箱的振动信号;
步骤2采用基于频谱趋势的方式对变分模态分解方法进行改进,利用改进的变分模态分解方法对齿轮箱的振动信号进行分解,得到振动信号的各个模态分量;
步骤3对获得的模态分量进行解调分析,获得各个模态信号的特征频率;
步骤4将特征频率与齿轮箱的各个齿轮转频进行对比,定位齿轮箱的故障源。
所述的一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法,其特征在于;步骤2)采用基于频谱趋势的改进变分模态分解方法来对齿轮箱的振动信号进行分解,得到振动信号的各个模态分量,具体步骤如下:
2.1)对振动信号进行Fourier变换,获得振动信号的频谱;
2.2)利用经验模态分解(EMD)算法对振动信号的频谱进行自适应性的分解,获得信号频谱的本征模态函数(IMFs)及其残余;
2.3)对部分IMFs和残余进行重构,获得振动信号的频谱趋势;
2.4)统计频谱趋势中的极大值个数M并将所有极大值的横坐标并归一化为xm
2.5)设置模态个数估计量K=M,初始的中心频率
Figure BDA0001520780430000031
2.6)利用步骤2.5)中设置的参数对振动信号进行变分模态分解,最终得到振动信号的各个模式函数。
所述的一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法,其特征在于步骤2.3)中的部分IMFs一般为最后4个或5个。
所述的一种基于改进经验小波变换的齿轮箱故障识别方法,其特征在于步骤3中的解调分析的方法包括基于Hilbert变换的包络解调方法和基于Teager能量算子的解调方法。
所述的一种基于改进经验小波变换的齿轮箱故障识别方法,其特征在于步骤2.5)中将频谱趋势中极大值的横坐标归一化之后作为初始化的中心频率。
通过采用上述技术,与现有技术相比,本发明的有益效果体现在:
1)本发明通过检测频谱趋势极大值个数,自适应地解决了VMD算法中模态个数预设的问题;
2)本发明通过将频谱趋势中极大值的横坐标归一化之后作为初始化的中心频率,使初始化的中心频率更符合振动信号频谱的分布,从而更具适应性;
3)本发明通过ST-VMD方法分解得到的模态分量在频谱上的分布更加合理,提高了齿轮箱故障识别的精度。
附图说明
图1为本发明的基于改进变分模态分解的齿轮箱故障识别方法流程图;
图2为本发明的基于频谱趋势的改进变分模态分解算法流程图;
图3为本发明的汽车座椅水平驱动器(HDM)齿轮箱振动时域波形图;
图4为本发明的HDM齿轮箱振动信号傅里叶频谱图;
图5为本发明的傅里叶频谱的经验模态分解结果图;
图6为本发明的频谱趋势(幅值被乘以2)和振动信号的频谱图;
图7为本发明的各个模态分量时域图;
图8为本发明的各个模态分量频谱图;
图9为本发明的各个模态分量频谱分布图;
图10为本发明的ST-VMD分解得到模态分量的解调分解结果图。
具体实施方式
以下结合说明书附图对本发明作进一步的描述,但本发明的保护范围并不仅限于此:
本发明基于频谱趋势和变分模态分解的齿轮箱故障识别方法,主要是针对原VMD方法在处理具有复杂多分量AM-FM信号时,模态分量个数难以估计,以及分解得到的各个模态分量不合理的问题,提出了基于频谱趋势的改进变分模态分解(ST-VMD)方法。该方法将时域信号的趋势概念引入到振动信号的频谱范围内,因此齿轮箱振动信号频谱的整体形态变化被考虑进来,并利用EMD算法来提取振动信号的频谱趋势。频谱趋势每个波峰能够很好的反应一个AM-FM信号,波峰的个数(或频谱趋势极大值的个数)即为需要预设的模态分量个数,同时将频谱趋势中极大值的横坐标归一化之后作为初始化的中心频率,提高初始化的中心频率的适应性,也使得分解得到的模态分量在频谱上的分布更加合理。
如图1所示,本发明的基于频谱趋势和变分模态分解的齿轮箱故障识别方法,具体步骤如下:
步骤1,采集齿轮箱的振动信号,对齿轮箱振动信号采用基于频谱趋势的方式对变分模态分解方法(ST-VMD)进行改进,利用改进的变分模态分解方法对齿轮箱的振动信号进行分解,得到各个模态分量,其具体步骤如图2所示:
1.1)对齿轮箱振动信号进行傅里叶变换(Fourier变换),获得振动信号的频谱;
1.2)对步骤1.1)获得的频谱进行EMD分解,获得若干模态分量IMFs和残余;
1.3)将步骤1.2)中部分IMFs和残余进行重构,获得振动信号的频谱趋势;
1.4)统计频谱趋势中的极大值个数M和所有极大值的横坐标并归一化为xm
1.5)设置模态个数估计量K=M,初始的中心频率
Figure BDA0001520780430000051
1.6)利用步骤1.5)中设置的参数对振动信号进行VMD分解,最终得到振动信号的各个模态分量;
步骤2,对步骤1中获得各个模态分量进行解调分析,得到特征频率,解调的方法可以是基于Hilbert变换的包络解调方法,也可以是基于Teager能量算子的解调方法,在本发明中,对于实际的齿轮箱振动信号,当获得的模态分量比较合理时,Hilbert变换解调方法和能量算子解调方法均能获得较好的解调结果;
步骤3,将步骤2中解调分析获得的特征频率与齿轮箱中各个齿轮的转频进行比较,对故障源定位,从而实现齿轮箱的故障诊断。
实施例1
如图1-2所示,基于频谱趋势和变分模态分解的齿轮箱故障识别方法,包括如下步骤:
1)以采集的汽车座椅水平驱动器(Horizontal Driving Machine,HDM)齿轮箱的振动信号为例,其时域波形如图3所示,对其进行Fourier变换,获得其信号频谱如图4所示;
2)利用EMD算法对获得的HDM振动信号的频谱进行分解,获得的IMFs和残余如图5所示;
3)对最后4个IMFs和残余进行重构(即IMF10~IMF13与残余res之和),获得其频谱趋势(为了清楚的显示,频谱趋势的幅值被乘以2)和振动信号的傅里叶频谱如图6所示;
4)检测到频谱趋势的极大值点的个数为7,将7个极大值点的横坐标进行归一化为xm
5)设置模态个数估计量K=7,初始的中心频率
Figure BDA0001520780430000061
6)利用步骤5)中设置的参数对振动信号进行VMD分解,最终得到振动信号的各个模式函数,各个模态分量的时域图如图7所示,各个模态分量的频谱图如图8所示,图9显示的是各个模态分量在频谱中的分布;
7)利用Teager能量算子对各个模态分量进行解调,解调分析的结果如图10所示;
8)本发明图5、7、8、10中的纵坐标均为幅值,从图10可知,解调出来的特征频率约为53.5Hz,因此可以判断转频约为53.5Hz的齿轮存在故障。
本发明与原VMD算法相比,根据频谱趋势的极大值的个数,自适应的解决了模态分量需要预设的问题,此外利用极大值点的横坐标初始化的中心频率更加符合振动信号在频域的分布,分解得到的模态分量也更加合理,从而提高了齿轮箱故障识别的准确度。

Claims (3)

1.一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法,其特征在于包括如下步骤:
步骤1 采集目标齿轮箱的振动信号;
步骤2 采用基于频谱趋势的方式对变分模态分解方法进行改进,利用改进的变分模态分解方法对齿轮箱的振动信号进行分解,得到振动信号的各个模态分量,具体步骤如下:
2.1)对振动信号进行Fourier变换,获得振动信号的频谱;
2.2)利用经验模态分解EMD算法对振动信号的频谱进行自适应性的分解,获得信号频谱的本征模态函数IMFs及其残余;
2.3)对部分IMFs和残余进行重构,获得振动信号的频谱趋势;
2.4)统计频谱趋势中的极大值个数M并将所有极大值的横坐标并归一化为
Figure DEST_PATH_IMAGE002
2.5)设置模态个数估计量K=M,初始的中心频率
Figure DEST_PATH_IMAGE004
2.6)利用步骤2.5)中设置的参数对振动信号进行变分模态分解,最终得到振动信号的各个模式函数;
步骤3 对获得的模态分量进行解调分析,获得各个模态信号的特征频率,解调分析的方法包括基于Hilbert变换的包络解调方法和基于Teager能量算子的解调方法;
步骤4 将特征频率与齿轮箱的各个齿轮转频进行对比,定位齿轮箱的故障源。
2.根据权利要求1所述的一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法,其特征在于步骤2.3)中的部分IMFs为最后4个或5个。
3.根据权利要求2所述的一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法,其特征在于步骤2.5)中将频谱趋势中极大值的横坐标归一化之后作为初始化的中心频率。
CN201711408267.4A 2017-12-22 2017-12-22 一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法 Active CN107944199B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711408267.4A CN107944199B (zh) 2017-12-22 2017-12-22 一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711408267.4A CN107944199B (zh) 2017-12-22 2017-12-22 一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法

Publications (2)

Publication Number Publication Date
CN107944199A CN107944199A (zh) 2018-04-20
CN107944199B true CN107944199B (zh) 2020-12-01

Family

ID=61939715

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711408267.4A Active CN107944199B (zh) 2017-12-22 2017-12-22 一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法

Country Status (1)

Country Link
CN (1) CN107944199B (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110610115A (zh) * 2018-06-14 2019-12-24 西门子(中国)有限公司 机械故障检测方法、装置及计算机可读存储介质
CN108444698B (zh) * 2018-06-15 2019-07-09 福州大学 基于teo解调和随机共振的行星齿轮箱早期故障诊断方法
CN109269629A (zh) * 2018-08-03 2019-01-25 河海大学 基于经验小波变换的特高压并联电抗器振动信号分析方法
CN109238712A (zh) * 2018-09-12 2019-01-18 温州大学 一种自回归模型增强变分模态分解的风力发电机滚动轴承故障诊断方法
CN109443752B (zh) * 2018-10-17 2020-11-27 北京信息科技大学 一种基于vmd的齿轮振动信号降噪及故障诊断方法
CN109443719B (zh) * 2018-11-01 2020-05-19 河南理工大学 一种钻头振动信号在线虚拟测试方法及其系统
CN109238717A (zh) * 2018-11-09 2019-01-18 福州大学 一种基于vmd-omeda的齿轮箱复合故障诊断方法
CN110333076B (zh) * 2019-06-19 2021-01-26 电子科技大学 基于CNN-Stacking的轴承故障诊断方法
CN110542828A (zh) * 2019-07-17 2019-12-06 佛山电力设计院有限公司 一种配网故障定位方法及系统
CN110427916B (zh) * 2019-08-14 2020-05-01 苏州大学 中心频率收敛趋势作用下的故障诊断方法
CN110514294A (zh) * 2019-08-30 2019-11-29 鞍钢矿业爆破有限公司 一种基于emd与vmd的爆破振动信号降噪方法
CN111178318B (zh) * 2020-01-06 2023-07-11 东南大学 一种基于渐进式vmd的滚动轴承早期复合故障特征提取方法
CN111307277B (zh) * 2020-03-20 2021-10-01 北京工业大学 基于变分模态分解和预测性能的单模态子信号选择方法
CN113049252B (zh) * 2021-03-25 2023-04-14 成都天佑路航轨道交通科技有限公司 一种列车轴承箱的故障检测方法
CN117725539B (zh) * 2024-02-07 2024-04-26 东北大学 一种复杂工况下风机旋转部件的故障特征提取及分析方法
CN117962968B (zh) * 2024-03-29 2024-06-07 浙江省轨道交通运营管理集团有限公司 一种轨道交通故障检测方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015048078A1 (en) * 2013-09-24 2015-04-02 Sikorsky Aircraft Corporation Gear fault detection
CN106096198A (zh) * 2016-06-29 2016-11-09 潍坊学院 一种基于变分模式分解和谱峭度的包络分析方法
CN106153339A (zh) * 2016-06-29 2016-11-23 潍坊学院 一种基于变分模式分解滤波的包络分析方法
CN107192554A (zh) * 2017-05-19 2017-09-22 西安理工大学 一种风电机组滚动轴承的振动故障诊断方法
CN107316653A (zh) * 2016-04-27 2017-11-03 南京理工大学 一种基于改进的经验小波变换的基频检测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015048078A1 (en) * 2013-09-24 2015-04-02 Sikorsky Aircraft Corporation Gear fault detection
CN107316653A (zh) * 2016-04-27 2017-11-03 南京理工大学 一种基于改进的经验小波变换的基频检测方法
CN106096198A (zh) * 2016-06-29 2016-11-09 潍坊学院 一种基于变分模式分解和谱峭度的包络分析方法
CN106153339A (zh) * 2016-06-29 2016-11-23 潍坊学院 一种基于变分模式分解滤波的包络分析方法
CN107192554A (zh) * 2017-05-19 2017-09-22 西安理工大学 一种风电机组滚动轴承的振动故障诊断方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
变分模态分解在齿轮箱故障诊断中的应用;闫洪波 等;《内蒙古科技与经济》;20161115(第21(2016)期);第111页摘要,第111页第1栏第3段至112页第2栏第2段 *
基于改进变分模态分解的旋转机械故障时频分析;刘尚坤 等;《振动工程学报》;20161231;第29卷(第6期);第1119-1126页 *

Also Published As

Publication number Publication date
CN107944199A (zh) 2018-04-20

Similar Documents

Publication Publication Date Title
CN107944199B (zh) 一种基于频谱趋势和变分模态分解的齿轮箱故障识别方法
CN108458871B (zh) 一种基于改进经验小波变换的齿轮箱故障识别方法
Feng et al. Time-varying demodulation analysis for rolling bearing fault diagnosis under variable speed conditions
CN107505135B (zh) 一种滚动轴承复合故障提取方法及系统
CN111238808B (zh) 一种基于经验模态分解与改进变分模态分解的齿轮箱复合故障诊断方法
CN110186682B (zh) 基于分数阶变分模态分解的滚动轴承故障诊断方法
Wang et al. Fault diagnosis of diesel engine based on adaptive wavelet packets and EEMD-fractal dimension
CN107832525B (zh) 一种信息熵优化vmd提取轴承故障特征频率的方法及其应用
CN111238813B (zh) 一种强干扰下提取滚动轴承故障特征的方法
CN104236905A (zh) 一种轴承故障诊断方法
CN103499445A (zh) 一种基于时频切片分析的滚动轴承故障诊断方法
CN110806315B (zh) 一种基于倒位编辑的齿轮箱复合故障诊断方法
CN111257872B (zh) 基于Radon变换与拉普拉斯算子的微多普勒抑制方法
CN111769810B (zh) 一种基于能量峭度谱的流体机械调制频率提取方法
CN109883706A (zh) 一种滚动轴承局部损伤微弱故障特征提取方法
Zheng et al. Faults diagnosis of rolling bearings based on shift invariant K-singular value decomposition with sensitive atom nonlocal means enhancement
CN112287796A (zh) 基于VMD-Teager能量算子的辐射源识别方法
CN112945546A (zh) 一种齿轮箱复杂故障精确诊断方法
Sun et al. Application of wavelet soft threshold denoising algorithm based on EMD decomposition in vibration signals
CN114486263A (zh) 一种旋转机械滚动轴承振动信号降噪解调方法
CN106815434B (zh) 基于小波和emd的噪声量时频分析去噪效果评价方法
Xu et al. Rolling bearing fault feature extraction via improved SSD and a singular-value energy autocorrelation coefficient spectrum
CN111323233B (zh) 一种用于低速旋转机械故障诊断的局部均值分解方法
CN113449630A (zh) 一种改进调制双谱的轴承故障诊断方法、系统及介质
CN114880627B (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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20180420

Assignee: Yunhe County Lixin Arts & Crafts Co.,Ltd.

Assignor: JIANG University OF TECHNOLOGY

Contract record no.: X2023980037396

Denomination of invention: A Gearbox Fault Identification Method Based on Spectrum Trend and Variational Modal Decomposition

Granted publication date: 20201201

License type: Common License

Record date: 20230704

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20180420

Assignee: Wuzhou Wanshun Forging Machine Tool Factory Co.,Ltd.

Assignor: JIANG University OF TECHNOLOGY

Contract record no.: X2023980054147

Denomination of invention: A Gearbox Fault Identification Method Based on Spectral Trend and Variational Mode Decomposition

Granted publication date: 20201201

License type: Common License

Record date: 20231226

EE01 Entry into force of recordation of patent licensing contract