CN114781466B - 基于旋转机械振动信号谐波基频的故障诊断方法及系统 - Google Patents
基于旋转机械振动信号谐波基频的故障诊断方法及系统 Download PDFInfo
- Publication number
- CN114781466B CN114781466B CN202210700846.0A CN202210700846A CN114781466B CN 114781466 B CN114781466 B CN 114781466B CN 202210700846 A CN202210700846 A CN 202210700846A CN 114781466 B CN114781466 B CN 114781466B
- Authority
- CN
- China
- Prior art keywords
- frequency
- spectrum
- harmonic
- original
- fundamental frequency
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
- G01H1/003—Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines
-
- 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
-
- 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/02—Gearings; Transmission mechanisms
- G01M13/021—Gearings
-
- 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/02—Gearings; Transmission mechanisms
- G01M13/028—Acoustic or vibration analysis
-
- 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
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Acoustics & Sound (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (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
技术领域
本发明涉及故障诊断分析技术领域,特别是涉及一种基于旋转机械振动信号谐波基频的故障诊断方法及系统。
背景技术
在进行旋转机械设备状态监测时,采集振动信号后(通常为:加速度信号、速度信号和位移信号)需要分析设备故障特征频率,而大部分故障信号在频谱中都表现为大量丰富的谐波现象,这些谐波的基础频率可能为设备旋转工作频率,设备轴承部件的内圈、外圈、滚动体、保持架故障特征频率,齿轮箱的啮合频率,或者是和电气故障相关的工频特征频率。通过检测识别振动信号的谐波基础频率,辨别谐波基础频率的类别进而确定故障类型是故障诊断中的常规诊断方法。由此可见谐波基础频率识别是故障诊断分析中的基础。
从旋转机械设备的振动信号提取谐波基础频率目前主要有以下几种做法:1)直接识别频谱中最大幅值对应的频率作为基础频率,这种方法计算简单高效,但是往往频谱幅值最大的频率往往并不一定是谐波基础频率,可能为任一高次谐波对应的频率值甚至是其他无关频率,因此这种方法易导致识别错误;2)识别频谱中前几大幅值对应的频率值作为谐波基础频率候选集,通过判断这些频率的倍数关系,确定谐波基础频率;3)通过谐波求和计算各次谐波幅值加权求和,以谐波加权求和的最大值所对应的基础频率作为谐波基频结果,这种方式需设定谐波次数,同时谐波求和的物理含义不明确导致对应的加权求和最大值的基础频率并不一定为实际的基础频率;4)发明专利CN109410980A针对音频谐波信号,使用常Q变换,得到常Q变换谱,并使用求和方法获取基础频率的初始候选集,最后根据最大后验概率准则提取谐波信号基频;5)通过其他信号处理转换,包括自相关法、倒谱法、小波变换法等提取信号的中最大幅值的频率作为基础频率。上述这些方法中本质核心其实都是对谐波频率幅值的判断,进而识别基础频率。但是由于计算的物理意义不明确,以及谐波基础频率确定准则不一致,导致各个方法的识别效果千差万别,并且各个方法并不能直接应用于旋转机械设备的振动信号谐波基础频率提取中。
可见,由于旋转机械振动信号谐波基础频率现有识别准则不一致,识别准确率依赖于各个方法或对信号的先验参数设定,导致了实际在旋转机械振动信号中谐波基础频率识别普遍存在准确率低的问题。
发明内容
为解决或至少缓解上述问题,本发明提出一种基于旋转机械振动信号谐波基频的故障诊断方法及系统,能够准确识别旋转机械振动信号谐波基频,从而实现旋转机械设备的精准故障诊断。
为实现上述目的,本发明提供了如下方案:
一种基于旋转机械振动信号谐波基频的故障诊断方法,包括:
获取旋转机械振动信号并计算所述旋转机械振动信号的原始频谱;
计算所述原始频谱在各个频率下的噪声水平作为噪声频谱;
根据所述原始频谱和所述噪声频谱计算归一化谐波乘积谱;
提取所述归一化谐波乘积谱中最大幅值对应的频率作为谐波基频基准值;
基于所述谐波基频基准值获取实际的谐波基础频率;
根据所述谐波基础频率进行旋转机械设备的故障诊断,得到故障诊断结果。
可选地,所述获取旋转机械振动信号并计算所述旋转机械振动信号的原始频谱,具体包括:
获取旋转机械振动信号;所述旋转机械振动信号包括旋转机械设备的加速度信号、速度信号和位移信号;
采用频谱分析方法、包络解调分析方法或小波分析方法中的一种或几种方法计算所述旋转机械振动信号的原始频谱;所述原始频谱包括频率以及与频率对应的幅值;其中为所述原始频谱中的第个频率数据,为所述原始频谱中的第个幅值数据。
可选地,所述计算所述原始频谱在各个频率下的噪声水平作为噪声频谱,具体包括:
可选地,所述根据所述原始频谱和所述噪声频谱计算归一化谐波乘积谱,具体包括:
计算所述原始频谱幅值谐波乘积谱与所述噪声频谱幅值乘积谱对应值的比值,作为所述归一化谐波乘积谱。
可选地,所述基于所述谐波基频基准值获取实际的谐波基础频率,具体包括:
获取所述原始频谱中最大幅值对应的频率作为实际疑似谐波基频;
若所述谐波基频基准值大于等于所述实际疑似谐波基频,判断所述谐波基频基准值除以所述实际疑似谐波基频的结果的小数部分是否小于等于0.1,若满足小于等于0.1,则确定实际的谐波基础频率为所述实际疑似谐波基频,若满足大于0.1,则确定实际的谐波基础频率为所述谐波基频基准值;
若所述谐波基频基准值小于所述实际疑似谐波基频,则确定实际的谐波基础频率为所述谐波基频基准值。
一种基于旋转机械振动信号谐波基频的故障诊断系统,包括:
原始频谱计算模块,用于获取旋转机械振动信号并计算所述旋转机械振动信号的原始频谱;
噪声频谱计算模块,用于计算所述原始频谱在各个频率下的噪声水平作为噪声频谱;
归一化谐波乘积谱计算模块,用于根据所述原始频谱和所述噪声频谱计算归一化谐波乘积谱;
谐波基频基准值提取模块,用于提取所述归一化谐波乘积谱中最大幅值对应的频率作为谐波基频基准值;
谐波基础频率获取模块,用于基于所述谐波基频基准值获取实际的谐波基础频率;
设备故障诊断模块,用于根据所述谐波基础频率进行旋转机械设备的故障诊断,得到故障诊断结果。
可选地,所述原始频谱计算模块具体包括:
旋转机械振动信号获取单元,用于获取旋转机械振动信号;所述旋转机械振动信号包括旋转机械设备的加速度信号、速度信号和位移信号;
原始频谱计算单元,用于采用频谱分析方法、包络解调分析方法或小波分析方法中的一种或几种方法计算所述旋转机械振动信号的原始频谱;所述原始频谱包括频率以及与频率对应的幅值;其中为所述原始频谱中的第个频率数据,为所述原始频谱中的第个幅值数据。
可选地,所述噪声频谱计算模块具体包括:
可选地,所述归一化谐波乘积谱计算模块具体包括:
归一化谐波乘积谱计算单元,用于计算所述原始频谱幅值谐波乘积谱与所述噪声频谱幅值乘积谱对应值的比值,作为所述归一化谐波乘积谱。
可选地,所述谐波基础频率获取模块具体包括:
实际疑似谐波基频获取单元,用于获取所述原始频谱中最大幅值对应的频率作为实际疑似谐波基频;
第一谐波基础频率获取单元,用于若所述谐波基频基准值大于等于所述实际疑似谐波基频,判断所述谐波基频基准值除以所述实际疑似谐波基频的结果的小数部分是否小于等于0.1,若满足小于等于0.1,则确定实际的谐波基础频率为所述实际疑似谐波基频,若满足大于0.1,则确定实际的谐波基础频率为所述谐波基频基准值;
第二谐波基础频率获取单元,用于若所述谐波基频基准值小于所述实际疑似谐波基频,则确定实际的谐波基础频率为所述谐波基频基准值。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种基于旋转机械振动信号谐波基频的故障诊断方法及系统,所述方法包括:获取旋转机械振动信号并计算所述旋转机械振动信号的原始频谱;计算所述原始频谱在各个频率下的噪声水平作为噪声频谱;根据所述原始频谱和所述噪声频谱计算归一化谐波乘积谱;提取所述归一化谐波乘积谱中最大幅值对应的频率作为谐波基频基准值;基于所述谐波基频基准值获取实际的谐波基础频率;根据所述谐波基础频率进行旋转机械设备的故障诊断,得到故障诊断结果。采用本发明方法能够准确识别旋转机械振动信号谐波基频,从而实现旋转机械设备的精准故障诊断。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一种基于旋转机械振动信号谐波基频的故障诊断方法的流程图;
图2为本发明方法基于旋转机械振动信号进行谐波基频识别的过程示意图;
图3为本发明实施例提供的某旋转机械速度振动信号的原始频谱和噪声频谱示意图;
图4为本发明实施例提供的某旋转机械速度振动信号的归一化谐波乘积谱示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种基于旋转机械振动信号谐波基频的故障诊断方法及系统,能够准确识别旋转机械振动信号谐波基频,从而实现旋转机械设备的精准故障诊断。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明一种基于旋转机械振动信号谐波基频的故障诊断方法的流程图,图2为本发明方法基于旋转机械振动信号进行谐波基频识别的过程示意图。参见图2,本发明方法基于旋转机械振动信号的频谱数据,频谱数据来源于直接的FFT分析或其他信号提取分析方法所得,通过获取频谱在各个频率下的噪声水平作为噪声频谱;之后计算归一化的谐波乘积谱,归一化主要体现在采用原始频谱幅值与噪声频谱的比值,同时在计算频谱幅值乘积之后经过开方计算明确乘积谱的物理含义,为基于幅值水平提取谐波基础频率提供了基础;提取谐波乘积谱中的最大幅值的频率作为基础频率的基准值(称为谐波基频基准值),最后基于实际频谱(也称为原始频谱)中最大幅值对应的频率值(称为实际疑似谐波基频)进行校准,进而获取最终的谐波基础频率(简称谐波基频)。通过本发明的应用,最终可以准确识别旋转机械振动信号谐波基频,为旋转机械实现精准的故障诊断提供了技术基础。
参见图1,本发明一种基于旋转机械振动信号谐波基频的故障诊断方法,具体包括:
步骤1:获取旋转机械振动信号并计算所述旋转机械振动信号的原始频谱。
该步骤1主要计算旋转机械振动信号的频谱。本发明对旋转机械振动信号(包括加速度信号、速度信号、位移信号),采用频谱分析方法(包括直接FFT(快速傅立叶变换)分析),获取频谱;或采用包络解调分析方法获取加速度信号的包络解调分析频谱;或采用小波分析或其他分析方法得到频谱,频谱包括频率数据与对应频率的幅值数据。
因此,所述步骤1获取旋转机械振动信号并计算所述旋转机械振动信号的原始频谱,具体包括:
获取旋转机械振动信号;所述旋转机械振动信号包括旋转机械设备的加速度信号、速度信号和位移信号;
采用频谱分析方法、包络解调分析方法或小波分析方法中的一种或几种方法计算所述旋转机械振动信号的原始频谱;所述原始频谱包括频率以及与频率对应的幅值;其中为所述原始频谱中的第个频率数据,为所述原始频谱中的第个幅值数据;则为的集合,为的集合。
步骤2:计算所述原始频谱在各个频率下的噪声水平作为噪声频谱。
该步骤2用于计算频谱各个频率下的噪声水平。获取频谱各个频率下的噪声,记为noise_y具体过程如下:
计算第个噪声值,表示第个噪声值,具体计算方法为:在中取第个幅值的前W个幅值和后W个幅值组成的幅值数组,计算该幅值数组的中位数作为。若前面不足W个数据,则直接取下标之前所有数据;若后面不足W个数据,则直接取下标之后所有数据。其中为幅值数据的下标值,从1开始;W取值范围为:1<W<10。
因此,所述步骤2计算所述原始频谱在各个频率下的噪声水平作为噪声频谱,具体包括:
步骤3:根据所述原始频谱和所述噪声频谱计算归一化谐波乘积谱。
图3为本发明实施例提供的某旋转机械速度振动信号的原始频谱和噪声频谱示意图。从图3中可发现该速度振动信号中明显包含丰富的谐波频率,因此各个谐波的幅值大小不一致,无法确保谐波基础频率一定是幅值最大所对应的频率,因此需要更精准和通用的谐波基频提取方法解决该问题。本发明计算原始频谱的谐波乘积谱与噪声频谱的谐波乘积谱的比值作为归一化谐波乘积谱,同时在计算频谱幅值乘积之后经过开方计算将乘积谱的物理含义明确,为提取谐波基础频率提供了基础。
所述步骤3根据所述原始频谱和所述噪声频谱计算归一化谐波乘积谱,具体包括:
步骤3.2:确定谐波乘积谱的频率上限:计算所有频率中最大值除以谐波乘积次数N的结果,并对结果向下取整数作为频率上限候选值,取步骤3.1中计算的所述初始频率上限和本步骤3.2计算的所述频率上限候选值中的最小值作为谐波乘积谱的频率上限;
步骤3.5:计算所述原始频谱幅值谐波乘积谱与所述噪声频谱幅值乘积谱对应值的比值,作为所述归一化谐波乘积谱。
步骤4:提取所述归一化谐波乘积谱中最大幅值对应的频率作为谐波基频基准值。
本发明中将谐波基频基准值记为base_freq。
步骤5:基于所述谐波基频基准值获取实际的谐波基础频率。
该步骤5基于谐波基频基准值base_freq获取实际频谱的基础频率值。所述步骤5基于所述谐波基频基准值获取实际的谐波基础频率,具体包括:
步骤5.1:获取所述原始频谱中最大幅值对应的频率作为实际疑似谐波基频,记为suspect_freq;
步骤5.2:若所述谐波基频基准值base_freq大于等于所述实际疑似谐波基频suspect_freq,判断所述谐波基频基准值base_freq除以所述实际疑似谐波基频suspect_freq的结果的小数部分是否小于等于0.1,若满足小于等于0.1,则确定实际的谐波基础频率为所述实际疑似谐波基频suspect_freq,若满足大于0.1,则确定实际的谐波基础频率为所述谐波基频基准值base_freq;
步骤5.3:若所述谐波基频基准值base_freq小于所述实际疑似谐波基频suspect_freq,则确定实际的谐波基础频率为所述谐波基频基准值base_freq。
步骤6:根据所述谐波基础频率进行旋转机械设备的故障诊断,得到故障诊断结果。
本发明方法实现了旋转机械振动信号谐波基频的准确识别。在实际应用中,通过辨别谐波基础频率的类别可以确定旋转机械设备的故障类型,作为旋转机械设备的故障诊断结果。
本发明方法基于旋转机械振动信号的频谱数据(频谱数据来源于直接的FFT分析或其他信号提取分析方法所得),通过获取频谱在各个频率下的噪声水平作为噪声频谱;之后计算原始频谱的谐波乘积谱(即原始频谱幅值乘积谱)与噪声频谱的谐波乘积谱(即噪声频谱幅值乘积谱)的比值作为归一化谐波乘积谱;提取谐波乘积谱中的最大频率作为基础频率的基准值,最后基于实际频谱最大幅值对应的频率值校准进而获取最终的谐波基础频率。通过本发明的应用,谐波提取过程被简化,同时应用乘积谱,并通过归一化和赋予乘积谱幅值实际物理含义,最终直接提取最大频率作为基础频率的基准值,整个过程无需依赖设备部件转速或者设备部件故障频率系数等信息,实现过程简单方便,具有非常高的通用性,针对所有频率数据皆可通过本发明方法计算频谱中隐藏的谐波基础频率,为旋转机械设备的故障诊断提供了可靠有效的谐波基频识别方式,从而提高了旋转机械设备的故障诊断准确性。
图4为本发明实施例提供的某旋转机械速度振动信号的归一化谐波乘积谱示意图。在该实施例中,通过转速测量装置获取该组信号对应的设备转速为935r/min,对应频率为15.58Hz,采用本发明方法,通过归一化谐波乘积谱最终识别的谐波基频为15.5Hz,与实际设备的转速频率误差低于1%。可见本发明方法在无需依赖外部参数输入的情况下实现了旋转机械设备谐波基础频率的准确识别,并且识别准确度高,为旋转机械设备的故障诊断分析提供了坚实基础,同时本发明方法识别谐波基频具备较高的通用性,具有广泛的应用前景。
基于本发明提供的方法,本发明还提供一种基于旋转机械振动信号谐波基频的故障诊断系统,所述系统包括:
原始频谱计算模块,用于获取旋转机械振动信号并计算所述旋转机械振动信号的原始频谱;
噪声频谱计算模块,用于计算所述原始频谱在各个频率下的噪声水平作为噪声频谱;
归一化谐波乘积谱计算模块,用于根据所述原始频谱和所述噪声频谱计算归一化谐波乘积谱;
谐波基频基准值提取模块,用于提取所述归一化谐波乘积谱中最大幅值对应的频率作为谐波基频基准值;
谐波基础频率获取模块,用于基于所述谐波基频基准值获取实际的谐波基础频率;
设备故障诊断模块,用于根据所述谐波基础频率进行旋转机械设备的故障诊断,得到故障诊断结果。
其中,所述原始频谱计算模块具体包括:
旋转机械振动信号获取单元,用于获取旋转机械振动信号;所述旋转机械振动信号包括旋转机械设备的加速度信号、速度信号和位移信号;
原始频谱计算单元,用于采用频谱分析方法、包络解调分析方法或小波分析方法中的一种或几种方法计算所述旋转机械振动信号的原始频谱;所述原始频谱包括频率以及与频率对应的幅值;其中为所述原始频谱中的第个频率数据,为所述原始频谱中的第个幅值数据。
所述噪声频谱计算模块具体包括:
所述归一化谐波乘积谱计算模块具体包括:
归一化谐波乘积谱计算单元,用于计算所述原始频谱幅值谐波乘积谱与所述噪声频谱幅值乘积谱对应值的比值,作为所述归一化谐波乘积谱。
所述谐波基础频率获取模块具体包括:
实际疑似谐波基频获取单元,用于获取所述原始频谱中最大幅值对应的频率作为实际疑似谐波基频;
第一谐波基础频率获取单元,用于若所述谐波基频基准值大于等于所述实际疑似谐波基频,判断所述谐波基频基准值除以所述实际疑似谐波基频的结果的小数部分是否小于等于0.1,若满足小于等于0.1,则确定实际的谐波基础频率为所述实际疑似谐波基频,若满足大于0.1,则确定实际的谐波基础频率为所述谐波基频基准值;
第二谐波基础频率获取单元,用于若所述谐波基频基准值小于所述实际疑似谐波基频,则确定实际的谐波基础频率为所述谐波基频基准值。
对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (4)
1.一种基于旋转机械振动信号谐波基频的故障诊断方法,其特征在于,包括:
获取旋转机械振动信号并计算所述旋转机械振动信号的原始频谱;
所述获取旋转机械振动信号并计算所述旋转机械振动信号的原始频谱,具体包括:
获取旋转机械振动信号;所述旋转机械振动信号包括旋转机械设备的加速度信号、速度信号和位移信号;
采用频谱分析方法、包络解调分析方法或小波分析方法中的一种或几种方法计算所述旋转机械振动信号的原始频谱;所述原始频谱包括频率以及与频率对应的幅值;其中为所述原始频谱中的第个频率数据,为所述原始频谱中的第个幅值数据;
计算所述原始频谱在各个频率下的噪声水平作为噪声频谱;
所述计算所述原始频谱在各个频率下的噪声水平作为噪声频谱,具体包括:
根据所述原始频谱和所述噪声频谱计算归一化谐波乘积谱;
所述根据所述原始频谱和所述噪声频谱计算归一化谐波乘积谱,具体包括:
计算所述原始频谱幅值谐波乘积谱与所述噪声频谱幅值乘积谱对应值的比值,作为所述归一化谐波乘积谱;
提取所述归一化谐波乘积谱中最大幅值对应的频率作为谐波基频基准值;
基于所述谐波基频基准值获取实际的谐波基础频率;
根据所述谐波基础频率进行旋转机械设备的故障诊断,得到故障诊断结果。
2.根据权利要求1所述的基于旋转机械振动信号谐波基频的故障诊断方法,其特征在于,所述基于所述谐波基频基准值获取实际的谐波基础频率,具体包括:
获取所述原始频谱中最大幅值对应的频率作为实际疑似谐波基频;
若所述谐波基频基准值大于等于所述实际疑似谐波基频,判断所述谐波基频基准值除以所述实际疑似谐波基频的结果的小数部分是否小于等于0.1,若满足小于等于0.1,则确定实际的谐波基础频率为所述实际疑似谐波基频,若满足大于0.1,则确定实际的谐波基础频率为所述谐波基频基准值;
若所述谐波基频基准值小于所述实际疑似谐波基频,则确定实际的谐波基础频率为所述谐波基频基准值。
3.一种基于旋转机械振动信号谐波基频的故障诊断系统,其特征在于,包括:
原始频谱计算模块,用于获取旋转机械振动信号并计算所述旋转机械振动信号的原始频谱;
所述原始频谱计算模块具体包括:
旋转机械振动信号获取单元,用于获取旋转机械振动信号;所述旋转机械振动信号包括旋转机械设备的加速度信号、速度信号和位移信号;
原始频谱计算单元,用于采用频谱分析方法、包络解调分析方法或小波分析方法中的一种或几种方法计算所述旋转机械振动信号的原始频谱;所述原始频谱包括频率以及与频率对应的幅值;其中为所述原始频谱中的第个频率数据,为所述原始频谱中的第个幅值数据;
噪声频谱计算模块,用于计算所述原始频谱在各个频率下的噪声水平作为噪声频谱;
所述噪声频谱计算模块具体包括:
归一化谐波乘积谱计算模块,用于根据所述原始频谱和所述噪声频谱计算归一化谐波乘积谱;
所述归一化谐波乘积谱计算模块具体包括:
归一化谐波乘积谱计算单元,用于计算所述原始频谱幅值谐波乘积谱与所述噪声频谱幅值乘积谱对应值的比值,作为所述归一化谐波乘积谱;
谐波基频基准值提取模块,用于提取所述归一化谐波乘积谱中最大幅值对应的频率作为谐波基频基准值;
谐波基础频率获取模块,用于基于所述谐波基频基准值获取实际的谐波基础频率;
设备故障诊断模块,用于根据所述谐波基础频率进行旋转机械设备的故障诊断,得到故障诊断结果。
4.根据权利要求3所述的基于旋转机械振动信号谐波基频的故障诊断系统,其特征在于,所述谐波基础频率获取模块具体包括:
实际疑似谐波基频获取单元,用于获取所述原始频谱中最大幅值对应的频率作为实际疑似谐波基频;
第一谐波基础频率获取单元,用于若所述谐波基频基准值大于等于所述实际疑似谐波基频,判断所述谐波基频基准值除以所述实际疑似谐波基频的结果的小数部分是否小于等于0.1,若满足小于等于0.1,则确定实际的谐波基础频率为所述实际疑似谐波基频,若满足大于0.1,则确定实际的谐波基础频率为所述谐波基频基准值;
第二谐波基础频率获取单元,用于若所述谐波基频基准值小于所述实际疑似谐波基频,则确定实际的谐波基础频率为所述谐波基频基准值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210700846.0A CN114781466B (zh) | 2022-06-21 | 2022-06-21 | 基于旋转机械振动信号谐波基频的故障诊断方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210700846.0A CN114781466B (zh) | 2022-06-21 | 2022-06-21 | 基于旋转机械振动信号谐波基频的故障诊断方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114781466A CN114781466A (zh) | 2022-07-22 |
CN114781466B true CN114781466B (zh) | 2022-09-27 |
Family
ID=82421578
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210700846.0A Active CN114781466B (zh) | 2022-06-21 | 2022-06-21 | 基于旋转机械振动信号谐波基频的故障诊断方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114781466B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116358864B (zh) * | 2023-06-01 | 2023-08-29 | 西安因联信息科技有限公司 | 一种旋转机械设备故障类型诊断方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109034277A (zh) * | 2018-09-17 | 2018-12-18 | 吉林大学 | 基于多特征融合的电能质量扰动分类方法及系统 |
CN109030628A (zh) * | 2018-06-25 | 2018-12-18 | 北京交通大学 | 基于振动噪声频谱特征匹配的轨道结构力学性能评价方法 |
CN109404285A (zh) * | 2018-09-13 | 2019-03-01 | 温州大学 | 一种螺杆式压缩机故障诊断的改进混洗蛙跳算法增强自适应带通滤波方法 |
CN111368129A (zh) * | 2018-12-25 | 2020-07-03 | 天津大学青岛海洋技术研究院 | 基于深度神经网络的哼唱检索法 |
CN113686576A (zh) * | 2021-08-17 | 2021-11-23 | 山东科技大学 | 基于非线性卷积稀疏滤波的旋转机械早期故障诊断方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20200302187A1 (en) * | 2015-07-17 | 2020-09-24 | Origin Wireless, Inc. | Method, apparatus, and system for people counting and recognition based on rhythmic motion monitoring |
CN106124043A (zh) * | 2016-06-17 | 2016-11-16 | 北京中元瑞讯科技有限公司 | 基于mcu和gprs的可组态输出智能传感系统 |
-
2022
- 2022-06-21 CN CN202210700846.0A patent/CN114781466B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109030628A (zh) * | 2018-06-25 | 2018-12-18 | 北京交通大学 | 基于振动噪声频谱特征匹配的轨道结构力学性能评价方法 |
CN109404285A (zh) * | 2018-09-13 | 2019-03-01 | 温州大学 | 一种螺杆式压缩机故障诊断的改进混洗蛙跳算法增强自适应带通滤波方法 |
CN109034277A (zh) * | 2018-09-17 | 2018-12-18 | 吉林大学 | 基于多特征融合的电能质量扰动分类方法及系统 |
CN111368129A (zh) * | 2018-12-25 | 2020-07-03 | 天津大学青岛海洋技术研究院 | 基于深度神经网络的哼唱检索法 |
CN113686576A (zh) * | 2021-08-17 | 2021-11-23 | 山东科技大学 | 基于非线性卷积稀疏滤波的旋转机械早期故障诊断方法 |
Non-Patent Citations (3)
Title |
---|
《Detection and Recovery of Fault Impulses via Improved Harmonic Product Spectrum and Its Application in Defect Size Estimation of Train Bearings》;Zhao M et al;《ScienceDirect》;20160705;全文 * |
《基于改进增强包络谱的齿轮箱轴承故障诊断方法研究》;唐若瑜;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20210715;正文第4、5章 * |
《船舶电力系统逆变模块故障检测与诊断研究》;田维;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20200715;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114781466A (zh) | 2022-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109668733B (zh) | 变分非线性模式分解变转速轴承故障诊断方法 | |
CN111238843B (zh) | 一种基于快速谱峭度分析的风机健康评价方法 | |
Sheng et al. | Applications in bearing fault diagnosis of an improved Kurtogram algorithm based on flexible frequency slice wavelet transform filter bank | |
CN112101174A (zh) | 一种基于LOF-Kurtogram的机械故障诊断方法 | |
CN102840907B (zh) | 早期故障状态下滚动轴承振动信号特征提取和分析方法 | |
CN111307438B (zh) | 一种基于信息熵的旋转机械振动故障诊断方法及其系统 | |
CN105628381A (zh) | 一种基于改进的局部均值分解的往复压缩机轴承故障诊断方法 | |
CN110160765A (zh) | 一种基于声音或振动信号的冲击特征识别方法及系统 | |
CN113086799B (zh) | 一种基于图像识别的电梯故障检测及预警方法 | |
Zhao et al. | Vibration health monitoring of rolling bearings under variable speed conditions by novel demodulation technique | |
CN114781466B (zh) | 基于旋转机械振动信号谐波基频的故障诊断方法及系统 | |
CN112364706A (zh) | 一种基于类不平衡的小样本轴承故障诊断方法 | |
CN108398260B (zh) | 基于混合概率方法的齿轮箱瞬时角速度的快速评估方法 | |
CN114486263B (zh) | 一种旋转机械滚动轴承振动信号降噪解调方法 | |
CN107679013A (zh) | 基于eemd‑hht与时频重排结合的转速曲线估计方法 | |
CN112052712A (zh) | 一种电力设备状态监测与故障识别方法及系统 | |
Chen et al. | Rolling bearing fault feature extraction method using adaptive maximum cyclostationarity blind deconvolution | |
CN112033656A (zh) | 一种基于宽带谱处理的机械系统故障检测方法 | |
CN113345443A (zh) | 基于梅尔频率倒谱系数的海洋哺乳动物发声检测识别方法 | |
Xu et al. | Online detection method for variable load conditions and anomalous sound of hydro turbines using correlation analysis and PCA-adaptive-K-means | |
CN115064182A (zh) | 自适应梅尔滤波器在强噪声环境下风机故障特征识别方法 | |
CN113219333B (zh) | 一种电机故障诊断时的频谱参数处理方法 | |
CN112345247A (zh) | 一种滚动轴承的故障诊断方法及装置 | |
CN113889147A (zh) | 一种基于动态时间弯曲算法的机械系统声纹信号匹配方法 | |
CN107271184A (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 |