留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于经验小波变换与谱峭度的船舶轴系故障特征提取方法

何洋洋 王馨怡 董晶

何洋洋, 王馨怡, 董晶. 基于经验小波变换与谱峭度的船舶轴系故障特征提取方法[J]. 中国舰船研究, 2020, 15(X): 1–9 doi:  10.19693/j.issn.1673-3185.01771
引用本文: 何洋洋, 王馨怡, 董晶. 基于经验小波变换与谱峭度的船舶轴系故障特征提取方法[J]. 中国舰船研究, 2020, 15(X): 1–9 doi:  10.19693/j.issn.1673-3185.01771
He Y Y, Wang X Y, Dong J. Fault feature extraction method for marine shafting based on empirical wavelet transform-spectral kurtosis[J]. Chinese Journal of Ship Research, 2020, 15(0): 1–9 doi:  10.19693/j.issn.1673-3185.01771
Citation: He Y Y, Wang X Y, Dong J. Fault feature extraction method for marine shafting based on empirical wavelet transform-spectral kurtosis[J]. Chinese Journal of Ship Research, 2020, 15(0): 1–9 doi:  10.19693/j.issn.1673-3185.01771

基于经验小波变换与谱峭度的船舶轴系故障特征提取方法

doi: 10.19693/j.issn.1673-3185.01771
详细信息
    作者简介:

    何洋洋,男,1993年生,硕士,助理工程师。研究方向:旋转机械状态监测和故障诊断。E-mail:heyangyang_0507@163.com

    王馨怡,女,1994年生,硕士,助理工程师。研究方向:电磁兼容。E-mail:wxinyi2014@126.com

    董晶,女,1990年生,硕士,工程师。研究方向:电磁兼容。E-mail:dj19900805@zju.edu.cn

    通讯作者:

    何洋洋

  • 中图分类号: U664.21

Fault feature extraction method for marine shafting based on empirical wavelet transform-spectral kurtosis

  • 摘要:   目的  为了解决船舶轴系振动信号的冗余信息过多、故障特征提取率较低等问题,提出一种基于经验小波变换与谱峭度(EWT-SK)的故障特征提取方法。  方法  首先,利用经验小波变换(EWT)对原始信号进行处理,以分离冗余振动成分,从而解决经验模态分解(EMD)固有的端点效应和模态混叠问题;然后,基于谱峭度和相关系数,筛选模态函数并进行重构,以突出故障信息,从而提高信噪比;最后,利用谱峭度获得最优带通滤波器的参数,并对滤波之后的信号进行包络解调,从而完成故障诊断。  结果  根据实例分析验证结果:在信号分解方面,EWT方法的特征提取稳定性和效率更高,可以保证轴系故障信息的完整性;在去噪效果方面,采用EWT-SK方法之后,故障信号的峭度值为4.761 6,相关系数为0.708 8,信噪比为3.762 4,其表现优于EMD和变分模态分解(VMD)方法。  结论  EWT-SK方法具有良好的特征提取能力与噪声抑制能力,适用于船舶轴系的故障诊断。
  • 图  1  傅里叶轴的分割

    Figure  1.  Partitioning of Fourier axis

    图  2  船舶轴系故障的诊断流程图

    Figure  2.  Flow chart of marine shafting fault diagnosis

    图  3  故障轴承模拟及采集装置图

    Figure  3.  Diagram of fault bearing simulation and acquisition device

    图  4  轴承内圈故障信号的时域、频域波形

    Figure  4.  Waveform of fault signal bearing inner ring in time domain and frequency domain

    图  5  EMD的前4阶模态时域图

    Figure  5.  Time-domain graph of the first fourth order IMFs using EMD

    图  6  EMD的前4阶模态频域图

    Figure  6.  Frequency-domain graph of the first fourth order IMFs using EMD

    图  7  VMD的模态函数时域图

    Figure  7.  Time-domain graph of IMFs using VMD

    图  8  VMD的模态函数频域图

    Figure  8.  Frequency-domain graph of IMFs using VMD

    图  9  故障信号的EWT支撑边界

    Figure  9.  EWT support boundary of fault signal

    图  10  EWT构造的带通滤波器幅频特性曲线

    Figure  10.  Amplitude frequency characteristic curve of bandpass filter constructed by EWT

    图  11  EWT的模态函数时域图

    Figure  11.  Time-domain graph of IMFs using EWT

    图  12  EWT的模态函数频域图

    Figure  12.  Frequency-domain graph of IMFs using EWT

    图  13  EMD-SK滤波信号时域图及包络谱

    Figure  13.  Time domain figure and envelope spectrum of signal filtered by EMD-SK

    图  14  VMD-SK滤波信号时域图及包络谱

    Figure  14.  Time domain figure and envelope spectrum of signal filtered by VMD-SK

    图  15  EWT重构信号的谱峭度图

    Figure  15.  Spectral kurtosis figure of EWT reconstructed signal

    图  16  EWT-SK滤波信号时域图及包络谱

    Figure  16.  Time domain figure and envelope spectrum of signal filtered by EWT-SK

    表  1  EWT分解模态函数的参数

    Table  1.   The IMFs' parameters of EWT

    模态峭度相关系数
    IMF13.588 60.333 1
    IMF24.154 90.543 7
    IMF33.669 80.440 4
    IMF45.297 20.351 7
    下载: 导出CSV

    表  2  3种方法的降噪对比结果

    Table  2.   Comparision of noise reduction results by three methods

    特征提取方法评价指标
    峭度相关系数信噪比
    EMD-SK3.325 50.212 22.349 9
    VMD-SK3.299 40.620 32.457 9
    EWT-SK4.761 60.708 83.762 4
    下载: 导出CSV
  • 董晨辰. 基于数据驱动的船舶旋转机械故障诊断方法研究及应用[D]. 南京: 南京航空航天大学, 2016.

    DONG C C. The research and application of fault diagnosis of ship rotation machinery based on data driven methods[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2016 (in Chinese).
    吕跃刚, 何洋洋. EWT和ICA联合降噪在轴承故障诊断中的应用[J]. 振动与冲击, 2019, 38(16): 42–48, 70.

    LV Y G, HE Y Y. Application of an EWT-ICA combined method in fault diagnosis of rolling bearings[J]. Journal of Vibration and Shock, 2019, 38(16): 42–48, 70 (in Chinese).
    ANTONI J. Fast computation of the kurtogram for the detection of transient faults[J]. Mechanical Systems and Signal Processing, 2007, 21(1): 108–124. doi:  10.1016/j.ymssp.2005.12.002
    苏文胜, 王奉涛, 张志新, 等. EMD降噪和谱峭度法在滚动轴承早期故障诊断中的应用[J]. 振动与冲击, 2010, 29(3): 18–21. doi:  10.3969/j.issn.1000-3835.2010.03.005

    SU W S, WANG F T, ZHANG Z X, et al. Application of EMD denoising and spectral kurtosis in early fault diagnosis of rolling element bearings[J]. Journal of Vibration and Shock, 2010, 29(3): 18–21 (in Chinese). doi:  10.3969/j.issn.1000-3835.2010.03.005
    程军圣, 杨怡, 杨宇. 基于LMD的谱峭度方法在齿轮故障诊断中的应用[J]. 振动与冲击, 2012, 31(18): 20–23, 54.

    CHENG J S, YANG Y, YANG Y. Application of spectral kurtosis approach based on local mean decomposition (LMD) in gear fault diagnosis[J]. Journal of Vibration and Shock, 2012, 31(18): 20–23, 54 (in Chinese).
    彭畅, 柏林, 谢小亮. 基于EEMD、度量因子和快速峭度图的滚动轴承故障诊断方法[J]. 振动与冲击, 2012, 31(20): 143–146.

    PENG C, BO L, XIE X L. Fault diagnosis method of rolling element bearings based on EEMD, measure-factor and fast kurtogram[J]. Journal of Vibration and Shock, 2012, 31(20): 143–146 (in Chinese).
    马增强, 李亚超, 阮婉莹, 等. 基于变分模态分解和谱峭度法的滚动轴承故障特征提取[J]. 电子测量与仪器学报, 2017, 31(11): 1782–1787.

    MA Z Q, LI Y C, RUAN W Y, et al. Rolling bearing fault feature extraction based on VMD and spectral kurtosis[J]. Journal of Electronic Measurement and Instrumentation, 2017, 31(11): 1782–1787 (in Chinese).
    GILLES J. Empirical wavelet transform[J]. IEEE Transactions on Signal Processing, 2013, 61(16): 3999–4010. doi:  10.1109/TSP.2013.2265222
    魏高飞, 张丹瑞. 基于局域均值分解的船舶轴系故障特征分析[J]. 上海船舶运输科学研究所学报, 2019, 42(1): 27–32.

    WEI G F, ZHANG D R. Fault feature analysis of ship shafting based on local mean decomposition[J]. Journal of Shanghai Ship and Shipping Research Institute, 2019, 42(1): 27–32 (in Chinese).
    何洋洋, 吕跃刚, 刘俊承. 基于VMD与粒子滤波的滚动轴承故障诊断[J]. 可再生能源, 2019, 37(1): 126–131. doi:  10.3969/j.issn.1671-5292.2019.01.020

    HE Y Y, LYU Y G, LIU J C. Fault diagnosis of rolling bearing based on VMD and particle filter[J]. Renewable Energy Resources, 2019, 37(1): 126–131 (in Chinese). doi:  10.3969/j.issn.1671-5292.2019.01.020
    唐贵基, 王晓龙. 参数优化变分模态分解方法在滚动轴承早期故障诊断中的应用[J]. 西安交通大学学报, 2015, 49(5): 73–81. doi:  10.7652/xjtuxb201505012

    TANG G J, WANG X L. Parameter optimized variational mode decomposition method with application to incipient fault diagnosis of rolling bearing[J]. Journal of Xi'an Jiaotong University, 2015, 49(5): 73–81 (in Chinese). doi:  10.7652/xjtuxb201505012
    田晶, 王英杰, 刘丽丽, 等. 基于Birge-Massart阈值降噪与EEMD及谱峭度的滚动轴承故障特征提取[J]. 航空动力学报, 2019, 34(6): 1399–1408.

    TIAN J, WANG Y J, LIU L L, et al. Fault feature extraction of rolling bearing using Birge-Massart threshold denoising with EEMD and spectral kurtosis[J]. Journal of Aerospace Power, 2019, 34(6): 1399–1408 (in Chinese).
    胡泽超, 何琳, 徐伟, 等. 船舶推进轴系纵向振动共振转换器的优化设计[J]. 中国舰船研究, 2019, 14(1): 107–113.

    HU Z C, HE L, XU W, et al. Optimzation design of resonance changer for marine propulsion shafting in longitudinal vibration[J]. Chinese Journal of Ship Research, 2019, 14(1): 107–113 (in Chinese).
  • [1] 周海军, 王超伟, 周国敬, 周博2020: 基于随机森林的离心泵滚动轴承故障诊断 , 中国舰船研究, 15, 129-135. doi: 10.19693/j.issn.1673-3185.01646
    [2] 范磊, 李维波, 徐聪, 李巍, 李齐2020: 基于晶闸管的自动调压变压器故障诊断技术 , 中国舰船研究, 15, 74-81,106. doi: 10.19693/j.issn.1673-3185.01532
    [3] 刘国强, 林叶锦, 张志政, 庞水2020: 基于粗糙集和优化DAG-SVM的船舶主机故障诊断方法 , 中国舰船研究, 15, 68-73. doi: 10.19693/j.issn.1673-3185.01650
    [4] 蒋佳炜, 胡以怀, 方云虎, 李方玉2020: 船舶动力装置智能故障诊断技术的应用与展望 , 中国舰船研究, 15, 56-67. doi: 10.19693/j.issn.1673-3185.01679
    [5] 康兴, 李维波, 余万祥, 徐聪, 何凯彦2019: 某舰用变压器在线监测与故障诊断系统研究 , 中国舰船研究, 14, 147-154. doi: 10.19693/j.issn.1673-3185.01434
    [6] 胡红钱, 施伟锋, 兰莹, 卓金宝2018: 基于以太网的船舶电力系统动态电能质量监测与故障诊断系统设计 , 中国舰船研究, 13, 120-126. doi: 10.3969/j.issn.1673-3185.2018.01.018
    [7] 魏一, 张跃文, 李斌2018: 基于DPC-GMM算法的船舶燃油系统故障诊断 , 中国舰船研究, 13, 147-153,165. doi: 10.19693/j.issn.1673-3185.01262
    [8] 任鑫, 孔衍, 周碧松, 张凯, 马丽2013: 改进遗传算法在船用核动力装置概率因果故障诊断中的应用 , 中国舰船研究, 8, 107-111. doi: 10.3969/j.issn.1673-3185.2013.01.017
    [9] 张国华, 郝刘仓, 曹振芳2012: 基于经验模态分解的大型轮船锚机齿轮箱故障诊断 , 中国舰船研究, 7, 75-80. doi: 10.3969/j.issn.1673-3185.2012.01.015
    [10] 田奕洋, 邱金水, 刘伯运2011: 基于图论的舰船水消防系统故障诊断顺序确定方法 , 中国舰船研究, 6, 72-74. doi: 10.3969/j.issn.1673-3185.2011.04.015
    [11] 刘峻华, 孟清正, 杨涛, 张聘亭2011: 船舶动力装置可组态智能故障诊断系统设计 , 中国舰船研究, 6, 77-80. doi: 10.3969/j.issn.1673-3185.2011.02.016
    [12] 崔志超, 阮礽忠, 徐华勇, 陈甜, 郑依妹2011: 基于小波分析的船舶同步发电机定子绕组故障在线诊断仿真研究 , 中国舰船研究, 6, 69-72,92. doi: 10.3969/j.issn.1673-3185.2011.02.014
    [13] 方芳2010: 感应电机转子断条故障诊断方法研究进展 , 中国舰船研究, 5, 77-81. doi: 10.3969/j.issn.1673-3185.2010.05.016
    [14] 张文胜, 龚涛, 胡锴2010: 基于D-S证据理论的直流电机故障诊断研究 , 中国舰船研究, 5, 72-74,80. doi: 10.3969/j.issn.1673-3185.2010.04.016
    [15] 方芳2010: 感应电机定子绕组短路故障诊断方法综述 , 中国舰船研究, 5, 98-101. doi: 10.3969/j.issn.1673-3185.2010.06.020
    [16] 王学杰, 贺国, 陈连树2009: 基于二叉树的某型船柴油机控制系统故障诊断专家系统 , 中国舰船研究, 4, 34-37. doi: 10.3969/j.issn.1673-3185.2009.06.008
    [17] 孙雅囡2008: 基于模糊神经网络和D-S证据理论的故障诊断方法研究 , 中国舰船研究, 3, 55-57. doi: 10.3969/j.issn.1673-3185.2008.01.013
    [18] 张威, 曾凡明, 苏启军, 袁利国, 刘金林2008: 综合层次分类方法在舰船动力装置故障诊断中的应用研究 , 中国舰船研究, 3, 34-37,46. doi: 10.3969/j.issn.1673-3185.2008.03.008
    [19] 黄冰阳, 陈奇, 王焕杰2008: 船舶轴系用液压紧配螺栓连接结构的有限元分析 , 中国舰船研究, 3, 45-49. doi: 10.3969/j.issn.1673-3185.2008.05.010
    [20] 周宗和, 杨自春2008: 基于小波包——时间序列分析的水位调节器故障诊断方法 , 中国舰船研究, 3, 61-64. doi: 10.3969/j.issn.1673-3185.2008.01.015
  • 加载中
图(16) / 表 (2)
计量
  • 文章访问数:  44
  • HTML全文浏览量:  4
  • PDF下载量:  718
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-09-11
  • 修回日期:  2019-12-07
  • 网络出版日期:  2020-12-10

基于经验小波变换与谱峭度的船舶轴系故障特征提取方法

doi: 10.19693/j.issn.1673-3185.01771
    作者简介:

    何洋洋,男,1993年生,硕士,助理工程师。研究方向:旋转机械状态监测和故障诊断。E-mail:heyangyang_0507@163.com

    王馨怡,女,1994年生,硕士,助理工程师。研究方向:电磁兼容。E-mail:wxinyi2014@126.com

    董晶,女,1990年生,硕士,工程师。研究方向:电磁兼容。E-mail:dj19900805@zju.edu.cn

    通讯作者: 何洋洋
  • 中图分类号: U664.21

摘要:   目的  为了解决船舶轴系振动信号的冗余信息过多、故障特征提取率较低等问题,提出一种基于经验小波变换与谱峭度(EWT-SK)的故障特征提取方法。  方法  首先,利用经验小波变换(EWT)对原始信号进行处理,以分离冗余振动成分,从而解决经验模态分解(EMD)固有的端点效应和模态混叠问题;然后,基于谱峭度和相关系数,筛选模态函数并进行重构,以突出故障信息,从而提高信噪比;最后,利用谱峭度获得最优带通滤波器的参数,并对滤波之后的信号进行包络解调,从而完成故障诊断。  结果  根据实例分析验证结果:在信号分解方面,EWT方法的特征提取稳定性和效率更高,可以保证轴系故障信息的完整性;在去噪效果方面,采用EWT-SK方法之后,故障信号的峭度值为4.761 6,相关系数为0.708 8,信噪比为3.762 4,其表现优于EMD和变分模态分解(VMD)方法。  结论  EWT-SK方法具有良好的特征提取能力与噪声抑制能力,适用于船舶轴系的故障诊断。

English Abstract

何洋洋, 王馨怡, 董晶. 基于经验小波变换与谱峭度的船舶轴系故障特征提取方法[J]. 中国舰船研究, 2020, 15(X): 1–9 doi:  10.19693/j.issn.1673-3185.01771
引用本文: 何洋洋, 王馨怡, 董晶. 基于经验小波变换与谱峭度的船舶轴系故障特征提取方法[J]. 中国舰船研究, 2020, 15(X): 1–9 doi:  10.19693/j.issn.1673-3185.01771
He Y Y, Wang X Y, Dong J. Fault feature extraction method for marine shafting based on empirical wavelet transform-spectral kurtosis[J]. Chinese Journal of Ship Research, 2020, 15(0): 1–9 doi:  10.19693/j.issn.1673-3185.01771
Citation: He Y Y, Wang X Y, Dong J. Fault feature extraction method for marine shafting based on empirical wavelet transform-spectral kurtosis[J]. Chinese Journal of Ship Research, 2020, 15(0): 1–9 doi:  10.19693/j.issn.1673-3185.01771
    • 船舶轴系包含动力系统主轴、推力轴、中间轴以及艉轴(螺旋桨轴)等旋转机械系统[1],整体规模庞大且结构复杂,各设备之间的耦合性较强。由于工作环境恶劣,所以轴系易产生各类故障,且诊断和维修的难度一般较高。目前,基于船舶轴系的振动信号时频域分析方法已逐步应用于船舶故障诊断领域,但由于船舶的强背景噪声干扰,复杂的传输路径以及信号采集、传递、衰减等因素影响[2],噪声干扰的幅值和能量将远远超过故障冲击信号。为了实现强背景噪声下的微弱特征信息提取,需对振动信号进行降噪处理。

      传统的共振解调方法一般基于经验和历史数据,不仅存在带通滤波器的参数(例如中心频率、通带宽度)选择问题,还存在效率低、不稳定等问题。Antoni[3]基于短时傅里叶变换(short-time fourier transform,STFT)和有限冲激响应(finite impulse response,FIR)滤波器提出了基于谱峭度(spectral kurtosis,SK)的分析方法,通过采用1/3-二叉树滤波器结构,并以谱峭度划分最优带通滤波器的中心频率和带宽,显著提高了故障检测效率;然而,该方法易受噪声影响,并不适用于信噪比较低的信号或随机脉冲噪声信号。

      在利用谱峭度在强背景噪声下进行特征信息提取方面,相继有学者提出了很多改进方法。苏文胜等[4]提出了基于经验模态分解(empirical mode decomposition,EMD)与谱峭度的分析方法,可以突出高频共振成分,从而改善低频干扰的影响;但由于EMD固有的模态混叠和端点效应问题,其对低信噪比信号的处理效果较差,且易出现故障辨别失效。程军圣等[5]和彭畅等[6]基于谱峭度分别在故障诊断中应用了局部均值分解(local mean decomposition,LMD)和集成经验模态分解(ensemble empirical mode decomposition,EEMD)方法,可以有效降低噪声干扰对谱峭度的影响;但该类递归算法没有从根本上解决模态混叠的问题,故仍未实现有效降噪,且其计算量有所增加。马增强等[7]提出了变分模态分解(variational mode decomposion,VMD)与谱峭度相结合的方法,不仅克服了模态混叠问题,还取得了良好的诊断效果;但VMD存在模态数、惩罚因子以及时间步长等参数的设定问题,所以无法实现自适应滤波降噪。Gilles[8]基于经验小波变换(empirical wavelet transform,EWT)提出了一种新型多分量信号分析方法,通过结合经验模态分解和小波理论的优势,不仅解决了EMD固有的端点效应和模态混叠问题,还具有计算量小、可靠性高、自适应性强等优点。

      鉴于船舶轴系振动信号是包含多个分量的复杂信号,且具有非线性、非平稳性的特点,本文拟提出经验模态分解与谱峭度( EMD-SK)相结合的特征提取方法:首先,利用EWT筛选出故障特征信息丰富的振动成分,完成初步降噪;其次,依据谱峭度获取最佳带通滤波参数;然后,对带通滤波信号进行包络解调,完成特征提取,从而实现故障诊断;最后,对实例轴承信号进行对比分析,用以验证本文方法的有效性。

    • 经验小波变换是基于小波框架和窄带信号分析理论的多分量信号分解方法,可以分离船舶轴系的复杂信号成分,其主要思想是根据频域极值点来自适应划分傅里叶谱,并通过构造正交小波滤波器组进行频谱滤波,进而提取不同尺度的本征模态函数(intrinsic mode function,IMF)。

    • 假设信号角频率的取值范围是[0,π],则N个单分量成分的信号需将傅里叶支撑[0,π]分割为N个连续空间。设各区间边界为ωn,即信号频谱内2个相邻极大值的中点,其中n=1,2,···,N。以ωn为中心,定义如图1中阴影区所示的过渡段,其宽度为Tn=2τn。将区间${\Lambda _n}$=[ωn−1ωn]内的带通滤波器定义为经验小波,则依据迈耶小波(Meyer wavelet)的构造方法,即可得到经验小波的尺度函数${\hat \phi _n}(\omega )$和小波函数${\hat \psi _n}(\omega )$

      图  1  傅里叶轴的分割

      Figure 1.  Partitioning of Fourier axis

    • 设信号$f(t)$的傅里叶变换为$F(\omega )$,其中t为时间,$\omega \in \left[ {0,{\rm{\pi }}} \right]$。采用阈值法确定$\left| {F(\omega )} \right|$的前N个极大值,利用各极值的频域角频率${\varOmega _n}$,即可求得频带边界角频率${\omega _n}$

      $${\omega _n} = \frac{{{\varOmega _{n + 1}} + {\varOmega _n}}}{2}$$ (1)

      因此,在傅里叶支撑[0,π]范围内,EWT通过信号自身的频率特性定义了N个窄带带通滤波器,其上、下截止频率分别为ωn+1ωn,且边界处满足起始点ω0=0,终点ωN=π,所以EWT分解算法具有很强的自适应性。

    • 参考经典小波变换方法,定义经验小波变换的细节系数$W_f^\varepsilon (n,t)$为信号$f(t)$与小波函数${\psi _n}(t)$的内积,且近似系数$W_f^\varepsilon (0,t)$定义为信号$f(t)$与尺度函数${\phi _1}(t)$的内积,则原信号$f(t)$可重构为

      $$\begin{split}& f(t) = W_f^\varepsilon (0,t) * {\phi _1}(t) + \sum\limits_{n = 1}^N {W_f^\varepsilon (n,t)*{\psi _n}(t)} = \\&\;\; {F^{ - 1}}\left[ {\hat W_f^\varepsilon (0,\omega ){{\hat \phi }_1}(\omega ) + \sum\limits_{n = 1}^N {\hat W_f^\varepsilon (n,\omega ){{\hat \psi }_n}(\omega )} } \right] \end{split}$$ (2)

      式中:${\rm{*}}$为卷积运算符;${F^{ - 1}}\left[ {} \right]$为傅里叶逆变换;$\hat W_f^\varepsilon (0,\omega )$${\hat \phi _1}(\omega )$$\hat W_f^\varepsilon (n,\omega )$${\hat \psi _n}(\omega )$分别为$W_f^\varepsilon (0,t)$${\phi _1}(t)$$W_f^\varepsilon (n,t)$${\psi _n}(t)$的傅里叶变换。

      EWT算法所分解的固有模态函数${f_0}(t)$${f_i}\left( t \right)$分别为

      $${f_0}(t) = W_f^\varepsilon (0,t)*{\phi _1}(t)$$ (3)
      $${f_i}(t) = W_f^\varepsilon (n,t)*{\psi _i}(t), i \geqslant 1$$ (4)

      式中:${f_i}\left( t \right)$为所分解的第i个模态函数;${\psi _i}(t)$为相应的经验小波函数。

    • 实船振动信号包含船体振动、结构激振、螺旋桨噪声等多类噪声信号,具有复杂性、随机性等特点。由于包络谱分析方法通常仅适用于窄带信号,所以为了实现船舶轴系故障诊断,本文将引入谱峭度理论,即通过检测最大脉冲分量所在的频带,用以为包络谱分析自适应选取最优的频带。将谱峭度值${K_Y}(f)$定义为某频率下信号概率密度函数的峰值度量,其本质是令谱峭度最大化的带宽与中心频率组合的寻优过程。根据傅里叶变换窗口内各频域的谱线峭度值,基于谱峭度对瞬态冲击的敏感性,即可实现故障冲击信息的检测和定位[5]

      假设故障振动信号$X(t)$基于Wold-Cramer分解的系统激励响应为$Y(t)$,基于四阶谱阶距,其谱峭度${K_Y}(f)$

      $${K_Y}(f) = \frac{{{S_{4Y}}(f)}}{{S_{2Y}^2(f)}} - 2$$ (5)

      式中:${S_{2Y}}(f)$${S_{4Y}}(f)$分别为度量复包络能量的2阶和4阶瞬时矩。

      在谱峭度图中,本文将以频率为横坐标,分解层数k为纵坐标,设定频率分辨率${2^{ - (k + 1)}}$,并将采用不同颜色表示不同分解层、不同频率处的谱峭度。谱峭度最大色块层数所在的频率即为共振解调技术带通滤波的最优中心频率,而横坐标长度即为最佳带宽,因此,基于该理论可以快速、准确地求得带通滤波器的设计参数。

    • 船舶轴系振动信号中包含多个分量的复杂随机信号,其能量弱、信噪比低,所以早期故障的特征信息常常淹没于强噪声干扰之中[9]。由于传统共振解调方法存在滤波器参数选取困难的问题,且其在强背景噪声的干扰下,无法单独提取耦合于轴系振动信号中的准确周期故障信息,因此,本文将提出基于EWT与谱峭度的特征提取方法。首先,对原信号进行EWT分解,基于谱峭度和互相关系数筛选模态分量,用以削弱噪声,提高信噪比,从而避免不相关频段信息的干扰;然后,根据谱峭度准确定位瞬态成分所在的频带,获取客观滤波参数,从而避免人为因素的干扰。结合这2种方法,即可获得更好的噪声抑制效果和故障特征提取信息。如果将其应用于船舶轴系的故障诊断中,即可解决轴系故障特征信息被强背景噪声干扰而难以识别的问题,具体诊断步骤如图2所示。

      图  2  船舶轴系故障的诊断流程图

      Figure 2.  Flow chart of marine shafting fault diagnosis

      1) 对船舶轴系故障的原始采样信号进行EWT分解,从而分离多成分振动信息;

      2) 经EWT分解之后,计算各模态函数的谱峭度值和互相关系数;

      3) 筛选谱峭度值、互相关系数均较大的最优IMF分量,并构成合成信号,实现轴系信号的初步降噪;

      4) 对初步降噪信号进行快速谱峭度计算,获得带通滤波器的中心频率和带宽,并基于该参数对初步降噪的轴系振动信号进行带通滤波;

      5) 对滤波之后的信号进行包络解调,从而提取轴承故障特征频率并开展故障诊断。

    • 为了模拟船舶轴系故障,本文将以美国凯斯西储大学电气工程实验室的故障轴承模拟及采集装置(图3)作为实例分析对象,其加速度传感器分别安装于电机壳体的风扇端及驱动端。进行实船轴系数据采集时,可以在主轴、推力轴、中间轴及艉轴等多处结构布置横向和纵向测点,以便进一步对比分析。本实验轴承采用SKF6205型深沟球轴承,轴承直径D为39.04 mm,滚动体直径d为7.94 mm,滚动体个数Z为9,接触角α${0^ \circ }$。实验过程中,在轴承内圈人工布置了直径0.177 8 mm、深度0.279 4 mm的单点故障。该装置的采样频率为12 000 Hz,电机负荷为0.746 kW,电机转速nc为1 772 r/min。

      图  3  故障轴承模拟及采集装置图

      Figure 3.  Diagram of fault bearing simulation and acquisition device

      该装置轴承的转动频率fr

      $${f_{\rm{r}}} = \frac{{{n_{\rm{c}}}}}{{60}}$$ (6)

      由式(6)可知,fr为29.53 r/s。

      根据轴承故障特征频率的经验公式,内圈故障频率${f_{{\rm{ic}}}}$为159.91 Hz,其计算公式如下:

      $${f_{{\rm{ic}}}} = 0.5Z\left(1 + \frac{d}{D}\cos \alpha \right){f_{\rm{r}}}$$ (7)

      虽然实船轴系传感器的测点位置和能量传递损失与本实验装置存在一定差异,且两者所提取的周期振幅及高谐次频率也略有不同,但其具体差值主要受实船主机缸数、螺旋桨叶片数等参数的影响,因此,本实验装置模拟的船舶轴系故障振动信号仍然具备有效性和参考性。

      图4所示为轴承内圈故障信号的时域和频域波形:在时域图中,清晰地反映了周期性冲击故障,但由于背景噪声的影响,难以判断冲击间隔,且故障信息也模糊不清;在频域图中,干扰谱线较多,所以频率ncficncfr均湮没于背景噪声之中。

      图  4  轴承内圈故障信号的时域、频域波形

      Figure 4.  Waveform of fault signal bearing inner ring in time domain and frequency domain

    • 采用EMD方法对轴承内圈的故障信号进行处理,即可获得12个模态函数(IMF)和1个残余分量,共计13阶模态。为便于描述,本文仅列出了前4阶模态(IMF1,IMF2,IMF3,IMF4)的时域图和频域图,如图5图6所示。由于故障信号包含了大量的噪声成分,在边界效应和其他因素的影响下,进行经验模态分解时易发生模态混叠现象[10]图6中,IMF1几乎涵盖了整个频率范围,IMF3与IMF4的中心频率几乎重合。模态混叠将对各阶模态函数的物理意义判别产生影响,对于复杂的实船振动信号而言,这将导致部分特征信息的混叠、丢失问题。

      图  5  EMD的前4阶模态时域图

      Figure 5.  Time-domain graph of the first fourth order IMFs using EMD

      图  6  EMD的前4阶模态频域图

      Figure 6.  Frequency-domain graph of the first fourth order IMFs using EMD

    • 采用VMD方法对轴承内圈的故障信号进行处理,其脱离了EMD循环筛分的分解方式,可以在变分框架内通过搜寻非约束变分模型最优解[11],从而有效地分离各类频率成分。依据文献[11]的参数设定策略,设定模态数为4,惩罚因子为10 000,时间步长为0。故障信号经VMD方法处理之后,共分解为4阶模态分量,其时域图及频域图分别如图7图8所示。

      图  7  VMD的模态函数时域图

      Figure 7.  Time-domain graph of IMFs using VMD

      图  8  VMD的模态函数频域图

      Figure 8.  Frequency-domain graph of IMFs using VMD

      图7可知,各波形整体无明显的冲击脉冲,即各模态函数的谱峭度指标整体偏小,这将导致轴系故障信号中瞬变成分识别的可靠性有所降低。由图8可知,各模态的频域波形基本消除了模态混叠现象,各频段波形的层次清晰,可以用于下一步的信号特征提取。

    • 利用EWT方法处理故障信号时,其可以利用信号频域的稀疏特性,依据频域极大值自适应划分频带边界(图9),从而保证模态函数为窄带信号。将傅里叶支撑[0,π]分割为如图9所示的4个部分,则经验小波定义在各区间的4个带通滤波器(4种颜色区分的线条)如图10所示。

      图  9  故障信号的EWT支撑边界

      Figure 9.  EWT support boundary of fault signal

      图  10  EWT构造的带通滤波器幅频特性曲线

      Figure 10.  Amplitude frequency characteristic curve of bandpass filter constructed by EWT

      经验小波变换之后,故障信号共分解为4个模态分量,其时域图及频域图分别如图11图12所示。由图11可知,EWT可以自动估算理想的分解层数,从而获得比EMD更合理的模态函数,故其各个模态的时域波形更清晰,且周期性冲击序列的规律也更为明显。由图12可知,该方法避免了经验模态分解方法中欠包络、过包络等引起的模态混叠问题,且有效抑制了端点效应。由此可见,EWT方法可以有效地剔除环境噪声,故适用于船舶轴系故障的特征提取。

      图  11  EWT的模态函数时域图

      Figure 11.  Time-domain graph of IMFs using EWT

      图  12  EWT的模态函数频域图

      Figure 12.  Frequency-domain graph of IMFs using EWT

      对比以上3种分解方法,可以得出如下结论:

      1) 对于多激励源的船舶轴系振动信号,EWT的分解层数明显优于EMD和VMD,故其实际提取稳定性更高,这与EWT根据频域极大值来自适应分割傅里叶谱有关。EWT方法基于小波分析的理论框架,解决了EMD缺乏严格数学基础的问题,可以自适应确定模态函数的数量,从而避免虚假模态干扰。虽然VMD可以通过预设模态分解层数k来进行最优化分解,但还需解决船舶复杂轴系信号的参数设定问题。

      2) EWT可以保证轴系故障信息的完整性。该方法通过构造一组正交小波滤波器组提取调幅调频(AM-FM)成分,克服了EMD固有的端点效应及模态混叠缺陷,避免了相同成分的信号信息干扰模态函数所致的有效信息丢失问题,从而明显提高了船舶轴系有效信号的估计质量和还原精度。

      3) EWT可以快速高效地提取轴系故障特征。VMD需预先确定最佳模态数、二次惩罚因子和双重上升的时间步长,且各模态的中心频率和带宽均需在变分模型中进行迭代求解;而EWT采用自适应局部极大值划分了频带边界[12],通过较少的迭代次数显著降低了计算量,故EWT方法的分解效率优于VMD。

    • 根据EMD分解各模态函数的谱峭度和相关系数指标,选取谱峭度与相关系数最大的IMF1和IMF2来构成筛选合成信号以完成初步降噪,这2个模态含有最多的故障冲击信息与轴系振动信号成分。根据合成信号即可求得谱峭度图,可知其最大峭度频带的中心频率为1 750 Hz,频带宽度为500 Hz,即谱峭度图确定的最优带通滤波器的范围是[1 500 Hz,2 000 Hz]。在该频段范围内,信号的谱峭度最大,信噪比较高,包含了更多的故障冲击特征信息。基于该频段对初步降噪之后的信号进行带通滤波,其信号时域图和包络频谱如图13所示。

      图  13  EMD-SK滤波信号时域图及包络谱

      Figure 13.  Time domain figure and envelope spectrum of signal filtered by EMD-SK

      经EMD-SK滤波之后,虽然在包络频谱中可以辨别轴承二倍转频$2{f_{\rm{r}}}$和内圈故障特征频率${f_{{\rm{ic}}}}$,但由于混杂的噪声干扰,信号频谱的整体分辨率较低,且特征频率${f_{{\rm{ic}}}}$的倍频信息较模糊,所以船舶轴系故障诊断的可信度与准确度也相对较低。

    • 依据峭度及相关系数指标,选取最优分量IMF3和IMF4来构成筛选合成信号。通过对该信号进行谱峭度分析,由谱峭度图即可自适应确定最优带通滤波器的范围是[4 500 Hz,6 000 Hz]。对滤波之后的信号进行包络解调,其时域图和包络频谱如图14所示。

      图  14  VMD-SK滤波信号时域图及包络谱

      Figure 14.  Time domain figure and envelope spectrum of signal filtered by VMD-SK

      对比EMD-SK的分解结果,VMD-SK筛选合成信号的时域波形更清晰,振动规律较明显,且部分噪声也得到了抑制。在包络频谱中可以观察到一倍和二倍转频(${f_{\rm{r}}}$$2{f_{\rm{r}}}$),且轴承内圈的特征频率${f_{{\rm{ic}}}}$及其二倍频幅也更为凸显,但其故障特征的高倍频对比度较低。究其原因,在实际轴系故障诊断中,测点位置及能量传递损失等条件干扰将对强背景噪声的信号诊断结果造成影响。

    • EWT分解模态函数的参数指标如表1所示,根据峭度和相关系数选取了最优分量IMF2和IMF4,以保留更多的故障冲击信息,进而完成振动数据的初步降噪。由于EWT处理轴系故障信号时存在去噪不完全、特征信息识别不稳定等问题,故需基于谱峭度图开展进一步处理。谱峭度图横坐标为频率,纵坐标为分解层数k,如图15所示,当分解层数(窗长)为3时,在[0,750 ]Hz频带(中心频率fc=375 Hz)存在最大谱峭度值Kmax(无量纲参数),据此即可确定共振解调最优带通滤波器的参数。对滤波之后的信号进行包络谱分析,其时域图和包络图如图16所示。

      表 1  EWT分解模态函数的参数

      Table 1.  The IMFs' parameters of EWT

      模态峭度相关系数
      IMF13.588 60.333 1
      IMF24.154 90.543 7
      IMF33.669 80.440 4
      IMF45.297 20.351 7

      图  15  EWT重构信号的谱峭度图

      Figure 15.  Spectral kurtosis figure of EWT reconstructed signal

      图  16  EWT-SK滤波信号时域图及包络谱

      Figure 16.  Time domain figure and envelope spectrum of signal filtered by EWT-SK

      根据时域图和包络谱,可以对3种方法的去噪效果进行定性评估。与EMD-SK和VMD-SK相比,EWT-SK的分解时域波形受噪声干扰的影响较小,故其轴系振动冲击的特征也更为明显。在包络频谱中,较大幅值谱线位于轴承内圈特征频率${f_{{\rm{ic}}}}$及其二倍、三倍频率处,可以有效地诊断轴承内圈故障,且二倍转频($2{f_{\rm{r}}}$)也得到了凸显,所以故障特征的整体幅值略高于EMD与VMD,从而证明了采用EWT-SK方法提取船舶轴系故障信号微弱特征信息的有效性。

      表2所示为3种特征提取方法的降噪对比结果,由表2可知,经EWT-SK方法降噪之后的谱峭度值最大,这反映了滤波信号中含有较为丰富的轴系故障周期性特征信息;EWT-SK方法的相关系数最高,这代表了该信号与原始信号具有较高的相关性[13],故可降低异常非相关噪声对船舶故障诊断的干扰;EWT-SK方法的信噪比最高,这表明了该方法具有较强的噪声抑制能力。因此,本文提出的EWT-SK方法可以剔除船舶的强背景噪声,并且快速、准确、稳定地筛选故障特征信息,具备实船适用性。

      表 2  3种方法的降噪对比结果

      Table 2.  Comparision of noise reduction results by three methods

      特征提取方法评价指标
      峭度相关系数信噪比
      EMD-SK3.325 50.212 22.349 9
      VMD-SK3.299 40.620 32.457 9
      EWT-SK4.761 60.708 83.762 4
    • 针对船舶轴系的多成分复杂信号的故障诊断问题,本文提出了基于经验小波变换与谱峭度的特征提取方法, EWT方法可以分离船舶的冗余振动成分并实现初步降噪,谱峭度可以定位轴系的瞬态冲击分量,从而分离故障特征频段。通过实例验证分析,得到如下结论:

      1) EWT可以在小波框架下自适应划分频谱,根据实船轴系振动信号的特点来分离模态函数,在保留特征信息的同时也具有明确的物理意义,从而解决了EMD固有的端点效应和模态混叠问题。同时,较少的迭代次数也显著降低了计算工作量,不仅提高了分解效率,也避免了难以解释的虚假模态,从而提高了船舶轴系复杂信号的分解可靠性。

      2) 通过结合谱峭度、相关系数指标来筛选本征模态函数,不仅避免了盲目取舍模态函数所导致的故障特征丢失问题,还可以突出船舶故障信息,进而提高了信噪比。

      3) 谱峭度解决了传统共振解调中的参数选取问题,可以快速定位、度量船舶轴系故障信号中的瞬态冲击信息,不仅保留了较多的故障特征,还避免了故障冲击频段之外的噪声干扰,从而提高了诊断精度。

      4) EWT与谱峭度相结合的方法解决了低信噪比条件下谱峭度故障信息提取率低、稳定性差的问题,具有自适应性好、鲁棒性强、模态分解能力高等优点,可为船舶轴系振动信号的故障特征提取提供参考。

参考文献 (13)

目录

    /

    返回文章
    返回