CN116625681A - 一种基于短时傅里叶变换的谱幅值调制滚动轴承故障诊断方法 - Google Patents

一种基于短时傅里叶变换的谱幅值调制滚动轴承故障诊断方法 Download PDF

Info

Publication number
CN116625681A
CN116625681A CN202211106295.1A CN202211106295A CN116625681A CN 116625681 A CN116625681 A CN 116625681A CN 202211106295 A CN202211106295 A CN 202211106295A CN 116625681 A CN116625681 A CN 116625681A
Authority
CN
China
Prior art keywords
short
fourier transform
spectrum
time fourier
amplitude modulation
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.)
Pending
Application number
CN202211106295.1A
Other languages
English (en)
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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN202211106295.1A priority Critical patent/CN116625681A/zh
Publication of CN116625681A publication Critical patent/CN116625681A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00—Testing of machine parts
    • G01M13/04—Bearings
    • G01M13/045—Acoustic or vibration analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/72—Wind turbines with rotation axis in wind direction

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明公开了一种基于短时傅里叶变换的谱幅值调制滚动轴承故障诊断新方法，对于采集到的复杂轴承信号，通过短时傅里叶变换得到其时频域中的幅值和相位，保持相位不变，用不同的指数权重改变其幅值并将结果与原相位构成编辑谱，然后进行短时傅里叶逆变换得到修正信号，最后对修正信号进行希尔伯特变换和傅里叶变换得到平方包络谱，从而降低噪声的影响并突出故障特征频率。通过对平方包络谱进行分析，可以提取出故障特征频率，最终实现滚动轴承的故障诊断。

Description

一种基于短时傅里叶变换的谱幅值调制滚动轴承故障诊断 方法
技术领域
本发明属于旋转机械故障诊断领域，具体涉及一种利用短时傅里叶变换改进谱幅值调制算法的滚动轴承故障诊断方法。
背景技术
滚动轴承是旋转机械中使用最广泛的零部件，同时也是最容易产生故障和损坏的元件之一，因此，为了避免重大安全事故的发生，减少财产损失和人身伤亡，需要对滚动轴承进行状态检测和故障诊断。由于其恶劣的工作环境和复杂的背景噪声，传统的故障诊断方法往往难以识别设备的早期故障，因此，如何识别和提取复杂条件下设备的微弱故障是目前研究的主要方向。
在目前提出的众多故障诊断方法中，包络分析，即共振解调技术，是应用最为广泛且有效的方法之一。这一方法先对由激励产生的共振频带进行带通滤波处理，然后对滤波所得的信号进行包络解调，计算轴承的故障特征频率并与包络谱中的故障特征频率进行对比，从而判断是否存在相应类型的故障。
然而，共振频带的选取是目前面临的一个主要问题，基于谱峭度的峭度图方法可以找到一个最佳的解调频带，但如果仅从一个频带恢复原信号，就会忽略很多微弱的故障信息，所以有必要对原信号采用多频段滤波的方法来完全提取出故障信号。
为了解决这一问题，人们提出了一种新的非线性滤波信号分离方法谱幅值调制(SAM)，它通过对傅里叶变换得到的幅值谱赋予不同的指数权重后进行傅里叶逆变换得到修正信号，再通过修正信号的平方包络谱来识别故障特征。这一方法凭借其效率高、自适应强、计算简便，无需依赖于评价指标的优势，取得了显著的效果。但由于该方法是对整个时域信号进行傅里叶变换，在频域中存在明显的平均效应，导致其在复杂噪声影响下的修正信号平方包络谱模糊，难以识别故障特征频率及其倍频成分。
发明内容
为了弥补现有算法中的不足，本发明提供了一种基于短时傅里叶变换的谱幅值调制滚动轴承故障诊断方法。具体来说，该方法通过对振动信号进行短时傅里叶变换获得时频域中的幅值和相位，然后对幅值赋予不同的指数权重，并结合信号的原始相位进行短时傅里叶逆变换，从而得到一系列的修正信号，最后通过计算和分析修正信号的平方包络谱得到故障特征信息。
本发明包括如下步骤：
第一步：采集轴承振动信号x(τ)，其中采样频率为fs，采样点数为N，信号的短时傅里叶变换可以表达为
式中j为虚数单位，STFT表示短时傅里叶变换，g(τ-t)是移动的窗。将信号短时傅里叶变换的结果表达为幅值A(t,f)和相位的形式如下
第二步：保持相位不变，在-0.5≤MO≤1.5范围内每隔0.1取一个MO值作为指数权重并计算A(t,f)MO，得到21组不同幅值；
第三步：将得到的不同指数权重下的A(t,f)MO与原信号的相位组成新的编辑谱；
第四步：对得到的编辑谱利用短时傅里叶逆变换算法，得到一系列不同的修正信号，其公式为
这一公式也可以表达为
式中xm代表修正信号，iSTFT表示短时傅里叶逆变换。
第五步：对修正信号进行希尔伯特变换和傅里叶变换得到平方包络谱，进一步降低噪声的影响并突出故障特征频率；
第六步：对平方包络谱进行归一化处理，得到基于短时傅里叶变换的谱幅值调制；
第七步：根据转速计算相应的故障特征频率，并判断轴承故障类型。
相比于使用傅里叶变换计算幅值和相位，基于短时傅里叶变换的谱幅值调制能够同时获得信号时间和频率两方面的信息，在幅值调制过程中能够得到更加准确和全面的结果。运算结果显示，在复杂噪声环境下，这一方法相比于传统包络谱能够获得更加清晰的故障信息，极大程度地减小了噪声对故障特征频率识别的影响；相比于现有的使用傅里叶变换的谱幅值调制，基于短时傅里叶变换的谱幅值调制能够获得更加清晰的最大平方包络谱，从而更容易识别故障特征频率及其倍频成分。
附图说明
图1为本发明的滚动轴承故障诊断方法流程图；
图2为本发明中振动信号的波形和频谱；
图3为本发明中振动信号的包络谱；
图4为本发明中振动信号基于短时傅里叶变换的谱幅值调制结果的三维视图；
图5为本发明中振动信号基于短时傅里叶变换的谱幅值调制结果的二维视图；
图6为本发明中振动信号基于短时傅里叶变换的谱幅值调制结果的最大平方包络谱；
图7为本发明中振动信号基于短时傅里叶变换的谱幅值调制结果MO＝0.6时的平方包络谱。
具体实施方式
本发明方法以6307型号轴承的内圈故障信号为例，电机转速为1496r/min，采样频率为15360Hz。经计算后，该轴承的内圈故障特征频率fi＝122.74Hz。
首先选择带有内圈故障的滚动轴承箱，将电机启动起来，利用数据采集器对轴承振动信号进行数据采集，并将采集到的数据传输到电脑中，利用本方法进行后续的数据处理分析。
绘制原始波形，并利用快速傅里叶变换算法得到采集信号的频谱。图2为采集到的故障信号波形及其频谱。从信号的原始波形中无法观察到明显的周期性故障冲击，且频谱中的边频带成分也较难分辨。
然后利用包络解调的方法对振动信号进行处理，得到信号的包络谱，如图3所示。从图中可以看出，受到复杂噪声的影响，包络谱中的有效信息被淹没，很难进行准确的故障诊断。
对原始信号进行短时傅里叶变换，对结果分别求绝对值和相角获取信号时频域中的幅值A(t,f)和相位其形式用下式表示
式中j为虚数单位，STFT表示短时傅里叶变换。
保持求得的相位不变，设置指数权重MO范围为-0.5≤MO≤1.5且每隔0.1取一个MO值计算A(t,f)MO，得到21组不同幅值，再将得到的不同指数权重下的A(t,f)MO与原信号的相位/>组成21组新的编辑谱，对不同MO得到的编辑谱进行短时傅里叶逆变换，得到一系列不同的修正信号，其形式用下式表示
式中xm代表修正信号，iSTFT表示短时傅里叶逆变换。
分别对每组修正信号进行希尔伯特变换和傅里叶变换得到21组平方包络谱，然后进行归一化处理得到归一化后的平方包络谱，设置x轴为频率，z轴为归一化后的幅值，将不同MO值下的结果沿y轴排列，得到振动信号基于短时傅里叶变换的谱幅值调制结果的三维视图如图4所示，同时得到图5为基于短时傅里叶变换的谱幅值调制结果的二维视图，即图4的俯视结果，图中可以明显观察到6307型号轴承内圈故障特征频率fi及其倍频所在峰值。
将沿MO轴方向观察所得的结果画在如图6所示的最大平方包络谱中，图中每个频率对应的幅值均为该频率下的最大幅值，可以视为最优MO值对应的结果，可以观察到图中故障特征频率fi及其倍频较为清晰。
单独提取MO＝0.6时的平方包络谱如图7所示，图中低频部分出现了明显的峰值(图7中标出的频率成分)，其大小与6307型号轴承内圈故障特征频率fi及其倍频相似，故可以认为该轴承内圈出现故障。
因此，本发明针对SAM易受复杂噪声影响的不足，提出了一种基于短时傅里叶变换的谱幅值调制方法，该方法操作简单，计算效率高，自适应性强，可以有效地对滚动轴承故障信号进行故障特征提取，且结果中故障特征频率突出，诊断效果较好。

Claims (1)

1.一种基于短时傅里叶变换的谱幅值调制滚动轴承故障诊断方法，其特征在于，该方法由以下步骤组成：
第一步：采集轴承振动信号x(t)，其中采样频率为fs，采样点数为N，利用短时傅里叶变换算法得到其时频域中的幅值A(t,f)和相位其形式为
式中j为虚数单位，STFT表示短时傅里叶变换；
第二步：保持相位不变，在-0.5≤MO≤1.5范围内每隔0.1取一个MO值作为指数权重并计算A(t,f)MO，得到21组不同幅值；
第三步：将得到的不同指数权重下的A(t,f)MO与原信号的相位组成新的编辑谱；
第四步：对得到的编辑谱利用短时傅里叶逆变换算法，得到一系列不同的修正信号，其形式为
式中xm代表修正信号，iSTFT表示短时傅里叶逆变换；
第五步：对修正信号进行希尔伯特变换和傅里叶变换得到平方包络谱，进一步降低噪声的影响并突出故障特征频率；
第六步：对平方包络谱进行归一化处理，得到基于短时傅里叶变换的谱幅值调制；
第七步：根据转速计算相应的故障特征频率，并判断轴承故障类型。
CN202211106295.1A 2022-09-11 2022-09-11 一种基于短时傅里叶变换的谱幅值调制滚动轴承故障诊断方法 Pending CN116625681A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211106295.1A CN116625681A (zh) 2022-09-11 2022-09-11 一种基于短时傅里叶变换的谱幅值调制滚动轴承故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211106295.1A CN116625681A (zh) 2022-09-11 2022-09-11 一种基于短时傅里叶变换的谱幅值调制滚动轴承故障诊断方法

Publications (1)

Publication Number Publication Date
CN116625681A true CN116625681A (zh) 2023-08-22

Family

ID=87615814

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211106295.1A Pending CN116625681A (zh) 2022-09-11 2022-09-11 一种基于短时傅里叶变换的谱幅值调制滚动轴承故障诊断方法

Country Status (1)

Country Link
CN (1) CN116625681A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117786607A (zh) * 2024-02-28 2024-03-29 昆明理工大学 基于时频熵谱的变工况振动信号故障诊断方法、系统
CN117949487A (zh) * 2024-02-01 2024-04-30 中国科学院精密测量科学与技术创新研究院 基于希尔伯特变换和傅里叶基追踪谱的nmr检测方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117949487A (zh) * 2024-02-01 2024-04-30 中国科学院精密测量科学与技术创新研究院 基于希尔伯特变换和傅里叶基追踪谱的nmr检测方法
CN117786607A (zh) * 2024-02-28 2024-03-29 昆明理工大学 基于时频熵谱的变工况振动信号故障诊断方法、系统
CN117786607B (zh) * 2024-02-28 2024-05-17 昆明理工大学 基于时频熵谱的变工况振动信号故障诊断方法、系统

Similar Documents

Publication Publication Date Title
CN116625681A (zh) 一种基于短时傅里叶变换的谱幅值调制滚动轴承故障诊断方法
CN109855874B (zh) 一种声音辅助振动微弱信号增强检测的随机共振滤波器
CN108388908B (zh) 基于k-svd和滑窗降噪的滚动轴承冲击故障诊断方法
CN107506710A (zh) 一种滚动轴承复合故障提取方法
CN111238813B (zh) 一种强干扰下提取滚动轴承故障特征的方法
CN108398260B (zh) 基于混合概率方法的齿轮箱瞬时角速度的快速评估方法
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
CN111504640B (zh) 一种加权滑动窗二阶同步压缩s变换轴承故障诊断方法
CN116361733A (zh) 一种故障诊断方法、装置、系统以及存储介质
CN114486263B (zh) 一种旋转机械滚动轴承振动信号降噪解调方法
CN106908232A (zh) 一种旋转机械的振动信号故障识别方法
CN117928951B (zh) 基于改进经验小波变换和包络谱能量比的故障诊断方法
CN112597958B (zh) 一种滚动轴承故障自动识别方法及系统
CN112903296B (zh) 一种滚动轴承故障检测方法及系统
CN111623968B (zh) 一种基于自适应形态滤波的故障特征提取方法
CN112597969A (zh) 一种滚动轴承故障诊断方法、系统及介质
CN109948487B (zh) 基于时频谱相关性分析的旋转机械故障特征提取方法
CN116304648A (zh) 基于优化脉冲增强与包络同步平均的齿轮故障识别方法
CN112747925B (zh) 一种基于复合形态学滤波的滚动轴承故障诊断方法
CN115356108A (zh) 一种调制高阶水平挤压变换机械故障诊断方法与装置
CN112781723B (zh) 一种基于频谱方差的谐波成分检测方法
Zhang et al. Fault diagnosis of rolling element bearing using ACYCBD based cross correlation spectrum
Shi et al. An auto instantaneous frequency order extraction method for bearing fault diagnosis under time-varying speed operation
CN112149534A (zh) 一种基于pchip-lcd滚动轴承故障特征提取方法
CN110263711B (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