CN115356108A - 一种调制高阶水平挤压变换机械故障诊断方法与装置 - Google Patents

一种调制高阶水平挤压变换机械故障诊断方法与装置 Download PDF

Info

Publication number
CN115356108A
CN115356108A CN202211231265.3A CN202211231265A CN115356108A CN 115356108 A CN115356108 A CN 115356108A CN 202211231265 A CN202211231265 A CN 202211231265A CN 115356108 A CN115356108 A CN 115356108A
Authority
CN
China
Prior art keywords
time
order
frequency
modulation
signal
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.)
Granted
Application number
CN202211231265.3A
Other languages
English (en)
Other versions
CN115356108B (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.)
Suzhou Apqi Internet Of Things Technology Co ltd
Suzhou Qironggu Technology Co ltd
Chengdu Apuqi Technology Co ltd
Original Assignee
Suzhou Apqi Internet Of Things Technology Co ltd
Suzhou Qironggu Technology Co ltd
Chengdu Apuqi Technology Co ltd
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 Suzhou Apqi Internet Of Things Technology Co ltd, Suzhou Qironggu Technology Co ltd, Chengdu Apuqi Technology Co ltd filed Critical Suzhou Apqi Internet Of Things Technology Co ltd
Priority to CN202211231265.3A priority Critical patent/CN115356108B/zh
Publication of CN115356108A publication Critical patent/CN115356108A/zh
Application granted granted Critical
Publication of CN115356108B publication Critical patent/CN115356108B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明提供一种调制高阶水平挤压变换机械故障诊断方法与装置，步骤如下：S10，获取信号，并对信号进行调制，得到调制信号
Figure 100004_DEST_PATH_IMAGE002
；S20，对调制信号
Figure 100004_DEST_PATH_IMAGE003
进行短时傅里叶变换，根据短时傅里叶变换时频值计算调制高阶群时延；S30，基于短时傅里叶变换和调制高阶群时延生成调制高阶水平挤压变换；S40，基于调制高阶水平挤压变换的时频谱进行机械故障诊断。本申请针对高阶时间重排同步挤压变换的时频聚焦性会受短时傅里叶的结果直接影响的问题，采用解调变换改善基于短时傅里叶变换的高阶时间重排同步挤压变换的时频精度。

Description

一种调制高阶水平挤压变换机械故障诊断方法与装置
技术领域
本发明涉及机械故障诊断领域，具体而言，涉及一种调制高阶水平挤压变换机械故障诊断方法与装置。
背景技术
随着现代化生产的发展和科学技术的不断进步，机械设备发展日趋大型化、高速化、集成化、自动化、精密化和复杂化，且数量也是迅速增长，随之带来的就是对其要求的高的成本效益、实用性、可靠性。大型机械由于运行工作环境恶劣、持续的高负荷运作及使用寿命的限制，零部件易出现损伤甚至损坏。由于部件间的联系密切，一旦系统某个部位出现故障，将会引起连锁反应，关键设备部件的损坏甚至会导致整个系统瘫痪，造成巨大的经济损失，引起灾难性的后果。
轴承是大型机械中最关键的机械部件，也是最脆弱的机械部件。轴承缺陷会产生周期性机械脉冲，安装在机器上的状态监测传感器可以捕捉到这种脉冲，因此如何从状态监测信号中识别这种故障显得尤为关键。
据数据统计表明，由轴承、转子和齿轮引起的旋转机械故障占到总故障数的70%，而其中仅轴承便占到了30%以上。几十年来，对于旋转机械尤其是滚动轴承的故障诊断取得了很大的进展。其中常用的时频分析方法主要包括短时傅里叶变换（STFT）、连续小波变换（CWT）、S变换（ST）等，这些基于窗函数的方法都受海森伯格（Heisenberg）测不准原理的约束，其时间分辨率和频率分辨率无法同时得到提高，使得信号特征在时频谱上的表征相对较模糊，对轴承的故障诊断精度产生影响。近年来，重排（RM）及同步挤压变换（SST）等新兴后处理技术的广泛应用，极大地提高了传统时频分析的分辨率，使得时频表征趋于理想情况。然而，时频后处理技术很大程度上依赖于处理前的时频表征性能，无法准确表征非平稳信号的时变特征，故如何获得更加精确的时频分布，在提高参数估计精度的同时，获取高聚焦性能的时频表征结果依然是目前机械故障诊断领域的研究重点。
发明内容
有鉴于此，本发明提供一种调制高阶水平挤压变换机械故障诊断方法与装置，通过获取更精确的群时延估计值，提高时频能量聚集性，从而高效地识别机械故障。
第一方面，提出一种调制高阶水平挤压变换机械故障诊断方法，该方法包括以下步骤：
S10，获取信号，并对信号进行调制，得到调制信号
Figure DEST_PATH_IMAGE001
，其中
Figure 623303DEST_PATH_IMAGE002
表示信号频 率，
Figure DEST_PATH_IMAGE003
表示调制频率；
S20，对调制信号
Figure 992973DEST_PATH_IMAGE004
进行短时傅里叶变换，根据短时傅里叶变换时频值计算 调制高阶群时延；
S30，基于短时傅里叶变换和调制高阶群时延生成调制高阶水平挤压变换；
S40，基于调制高阶水平挤压变换的时频谱进行机械故障诊断。
作为优选：对信号进行调制，调制算子定义为
Figure DEST_PATH_IMAGE005
，
Figure 444814DEST_PATH_IMAGE006
为虚 数单位，
Figure 495816DEST_PATH_IMAGE007
表示信号的相位，
Figure 100002_DEST_PATH_IMAGE008
是相位调节算子，
Figure 513450DEST_PATH_IMAGE009
表示对相位调节算子求关 于调制频率
Figure 257284DEST_PATH_IMAGE003
的一阶偏导数。通过相位调制得到调制信号
Figure 100002_DEST_PATH_IMAGE010
，式中
Figure 727580DEST_PATH_IMAGE011
表示信号振幅，
Figure 100002_DEST_PATH_IMAGE012
表示时间。
作为优选：对所述调制信号
Figure 847851DEST_PATH_IMAGE004
进行短时傅里叶变换，得到时频值
Figure 47888DEST_PATH_IMAGE013
，其中
Figure 100002_DEST_PATH_IMAGE014
表示窗函数再利用所述时频值计算调制高阶群时延，计算过程如下：
首先，求得
Figure 244515DEST_PATH_IMAGE013
关于时间
Figure 926075DEST_PATH_IMAGE015
的
Figure 100002_DEST_PATH_IMAGE016
阶偏导数为：
Figure 866349DEST_PATH_IMAGE017
，
再次，根据上式对频率
Figure 100002_DEST_PATH_IMAGE018
求偏导数得到：
Figure 45526DEST_PATH_IMAGE019
，
式中，
Figure 100002_DEST_PATH_IMAGE020
。然后，利用上式 构造高阶方阵
Figure 944212DEST_PATH_IMAGE021
与
Figure 100002_DEST_PATH_IMAGE022
：
Figure 513734DEST_PATH_IMAGE023
，
Figure 100002_DEST_PATH_IMAGE024
，
式中,
Figure 241387DEST_PATH_IMAGE025
；
Figure 100002_DEST_PATH_IMAGE026
为窗函数
Figure 88121DEST_PATH_IMAGE027
中的调节因 子；
最后，利用高阶方阵
Figure 344658DEST_PATH_IMAGE021
与
Figure DEST_PATH_IMAGE028
计算调制高阶群时延
Figure 276842DEST_PATH_IMAGE029
：
Figure DEST_PATH_IMAGE030
,
式中，
Figure 11449DEST_PATH_IMAGE031
表示取复数的实部。
作为优选：基于短时傅里叶变换和调制高阶群时延得到调制高阶水平挤压变换
Figure DEST_PATH_IMAGE032
：
Figure 712689DEST_PATH_IMAGE033
,
式中，
Figure DEST_PATH_IMAGE034
为一时间变量，且为与t取值无关的不同时间变量，函数
Figure 936866DEST_PATH_IMAGE035
为Dirac函 数,
Figure DEST_PATH_IMAGE036
为阈值。
作为优选：基于调制高阶水平挤压变换时频谱进行机械故障诊断，包括：
基于调制高阶水平挤压变换时频谱确定目标脉冲；
基于所述目标脉冲的时频特征确定机械故障。
第二方面，提出一种调制高阶水平挤压变换机械故障诊断方法的装置，该装置包括：
变换模块，用于对机械信号进行调制，得到调制信号，然后进行短时傅里叶变换，生成所述机械信号的短时傅里叶时频谱；
时间模块，用于基于短时傅里叶变换时频谱生成所述机械信号在各个时频位置的群时延；
挤压模块，用于基于短时傅里叶变换和调制高阶群时延生成调制高阶水平挤压变换；
诊断模块，用于基于调制高阶水平挤压变换的时频谱进行机械故障诊断。
作为优选：对信号进行调制，调制算子定义为
Figure 480980DEST_PATH_IMAGE037
，通过 相位调制得到调制信号
Figure DEST_PATH_IMAGE038
。
作为优选：对所述调制信号
Figure 894643DEST_PATH_IMAGE004
进行短时傅里叶变换，得到时频值
Figure 919231DEST_PATH_IMAGE039
，再利用所述时频值计算调制高阶群时延，计算过程如下：
首先，求得
Figure 392938DEST_PATH_IMAGE013
关于时间
Figure 746384DEST_PATH_IMAGE015
的
Figure 698160DEST_PATH_IMAGE016
阶偏导数为：
Figure DEST_PATH_IMAGE040
，
再次，根据上式对频率
Figure 967467DEST_PATH_IMAGE018
求偏导数得到：
Figure 736709DEST_PATH_IMAGE041
，
然后，利用上式构造高阶方阵
Figure 334044DEST_PATH_IMAGE021
与
Figure 89510DEST_PATH_IMAGE028
：
Figure DEST_PATH_IMAGE042
，
Figure 603537DEST_PATH_IMAGE043
，
最后，再利用高阶方阵
Figure 559991DEST_PATH_IMAGE021
与
Figure 972518DEST_PATH_IMAGE028
计算调制高阶群时延
Figure DEST_PATH_IMAGE044
：
Figure 453047DEST_PATH_IMAGE045
,
作为优选：基于短时傅里叶变换和调制高阶群时延得到调制高阶水平挤压变换
Figure 306734DEST_PATH_IMAGE032
：
Figure DEST_PATH_IMAGE046
,
作为优选：基于调制高阶水平挤压变换时频谱进行机械故障诊断，包括：
计算机械信号的调制高阶水平挤压变换时频谱的时频包络；将所述时频包络中最大振幅对应的频率作为目标频率；获取所述目标频率对应的目标脉冲；将所述目标脉冲的时频特征和多个预设时频特征进行匹配；根据匹配结果确定所述机械的故障。
本申请提供一种调制高阶水平挤压变换机械故障诊断方法与装置，采用解调变换改善基于短时傅里叶变换的高阶时间重排同步挤压变换的时频精度，解决了高阶时间重排同步挤压变换的时频聚焦性会受短时傅里叶的结果直接影响的问题。
附图说明
图1为本发明流程图。
图2（a）为模拟信号波形图，（b）为模拟信号加噪后的波形图。
图3为图2经短时傅里叶变换时频谱图。
图4为图2经调制高阶水平挤压变换时频谱图。
图5（a）为轴承外圈故障信号的波形图,（b）为图5（a）中矩形区域内的信号的频谱图。
图6为轴承外圈故障信号进行（a）STFT、（b）SET、（c）HTSST、（d）DHTSST得到的时频谱，右侧（e）图至（h）图为（a）图至（d）图矩形框的局部放大图。
图7（a）为轴承外圈故障信号DHTSST的时频包络谱，（b）为提取的故障特征。
图8为机械故障诊断装置的框图。
具体实施方式
下面将结合附图对本发明作进一步说明。
第一方面，如图1，为本申请一种调制高阶水平挤压变换机械故障诊断方法的流程示意图，包括以下步骤：
S10，获取信号，并对信号进行调制，得到调制信号
Figure 948936DEST_PATH_IMAGE001
，其中
Figure 255284DEST_PATH_IMAGE002
表示信号频 率，
Figure 742766DEST_PATH_IMAGE003
表示调制频率；
S20，对调制信号
Figure 575593DEST_PATH_IMAGE004
进行短时傅里叶变换，根据短时傅里叶变换时频值计算 调制高阶群时延；
S30，基于短时傅里叶变换和调制高阶群时延生成调制高阶水平挤压变换；
S40，基于调制高阶水平挤压变换的时频谱进行机械故障诊断。
方法中对信号进行调制，是为了提高信号的时频表示的能量聚集性，从而能够更准确地研究信号的时频特征。
具体地，本申请的一种实施：
（1）根据定义的调制算子
Figure 608271DEST_PATH_IMAGE047
，式中
Figure 916761DEST_PATH_IMAGE006
为虚数单位，
Figure 100002_DEST_PATH_IMAGE048
表示信号的相位，
Figure 817721DEST_PATH_IMAGE049
是相位调节算子，
Figure 380421DEST_PATH_IMAGE009
表示对相位调节算子求关于调制频率
Figure 708634DEST_PATH_IMAGE003
的一阶偏导数。对信号进行相位调制得到调制信号
Figure 244700DEST_PATH_IMAGE038
，式 中
Figure 824718DEST_PATH_IMAGE011
表示信号振幅，
Figure 632137DEST_PATH_IMAGE012
表示时间。
（2）对所述调制信号
Figure 255885DEST_PATH_IMAGE004
进行短时傅里叶变换，得到时频值
Figure 100002_DEST_PATH_IMAGE050
，其中
Figure 24121DEST_PATH_IMAGE051
表示窗函数
Figure 100002_DEST_PATH_IMAGE052
再利用所述时频值计算调制高阶群时延，计算过程如下：
首先，求得
Figure 188255DEST_PATH_IMAGE013
关于时间
Figure 725546DEST_PATH_IMAGE015
的
Figure 395562DEST_PATH_IMAGE016
阶偏导数为：
Figure 900362DEST_PATH_IMAGE040
，
再次，根据上式对频率
Figure 681236DEST_PATH_IMAGE018
求偏导数得到：
Figure 338613DEST_PATH_IMAGE041
，
式中，
Figure 648372DEST_PATH_IMAGE053
。然后，利用上式 构造高阶方阵
Figure 640467DEST_PATH_IMAGE021
与
Figure 225033DEST_PATH_IMAGE028
：
Figure 736916DEST_PATH_IMAGE042
，
Figure 873369DEST_PATH_IMAGE043
，
式中,
Figure 962547DEST_PATH_IMAGE025
；
Figure 226170DEST_PATH_IMAGE026
为窗函数
Figure 920456DEST_PATH_IMAGE027
中的调节因 子；
最后，利用高阶方阵
Figure 493389DEST_PATH_IMAGE021
与
Figure 69864DEST_PATH_IMAGE022
计算调制高阶群时延
Figure 137177DEST_PATH_IMAGE044
：
Figure 607341DEST_PATH_IMAGE054
,
式中，
Figure 226541DEST_PATH_IMAGE031
表示取复数的实部。
（3）基于短时傅里叶变换和调制高阶群时延得到调制高阶水平挤压变换
Figure 165679DEST_PATH_IMAGE032
：
Figure 895737DEST_PATH_IMAGE046
,
式中，
Figure 237986DEST_PATH_IMAGE034
为一时间变量，且为与t取值无关的不同时间变量，函数
Figure 903454DEST_PATH_IMAGE035
为Dirac函 数，
Figure DEST_PATH_IMAGE055
为阈值。
（4）基于调制高阶水平挤压变换时频谱进行机械故障诊断，包括：
基于调制高阶水平挤压变换时频谱确定目标脉冲；
基于所述目标脉冲的时频特征确定机械故障。
参见图2-图4，其中具体的：
步骤（1）中，本申请原始信号
Figure 110313DEST_PATH_IMAGE056
的时域图如图2，定义调制算子
Figure DEST_PATH_IMAGE057
，然后 通过相位调制得到信号
Figure 519429DEST_PATH_IMAGE004
，其短时傅里叶变换时频谱图如图3。
步骤（3）中，调制高阶水平挤压变换后的时频谱图如图4。
从图中可以看出，经过本发明的方法处理后，图4得到的时频谱图，明显优于图3的时频谱图，也就是说，经过本发明的方法处理后得到的时频谱图明显优于原始时频谱图。
参见图5-图7。在实验中，外圈故障的振动信号由加速度传感器进行采集，传感器被放置在驱动端的位置。此时采样频率为12kHz，电机轴的工作速度为1797rpm。根据轴承参数和转速，通过计算，外圈故障的特征频率为107.4Hz。采集到的外圈故障振动信号及其频谱图绘制到图5中（a）中，选取其中一部分信号来进行分析，选取区域为图5中（a）中的矩形框，图5中（b）展示了所选区域内信号的频谱。可以看出，信号的能量主要分布在频率2.5kHz~4kH之间，在频率3kHz左右频谱的能量发生了明显的变化。STFT的结果显示在图6中（a）中，白色矩形框部分的表征结果的局部放大图展示在图6中（e）中，由于模糊的时频能量，信号的脉冲特征信息分布在较大的区域范围内，甚至难以观察到明显的群时延变化信息。SET的结果显示在图6中（b）中，该方法白色矩形框部分的局部放大图结果显示在图6中（f）中，SET保留了时频特征最明显处的时频能量，因此SET能够大幅提高能量聚集性。但是SET将同一时间区间内的瞬态信号识别为了两个不同的AM—FM信号分量，因此SET不能用于提取瞬态信号的时频特征。HTSST和DHTSST的表征结果分别如图6中（c）、图6中（d）所示，它们的白色矩形框部分的局部放大图如图6中（g）、图6中（h）所示。HTSST与DTSST在提高时频分辨率的同时，保留了信号的脉冲特征。由于噪声干扰，HTSST在群时延脊线上估计出了一条模糊的曲线，而DHTSST在提高聚集性后的时频表示上再次进行挤压，因此很好地估计了群时延脊线。在图7中，由DHTSST提取的相邻脉冲特征的时间间隔为9.33ms，对应的故障频率为107.2HZ，与真实值一致，验证了DHTSST在提取故障特征方面的效果。
第二方面，如图8所示，机械故障诊断装置包括：变换模块，时间模块，挤压模块，诊断模块。图8是机械故障诊断装置的框图。
变换模块用于对机械信号进行调制，然后进行短时傅里叶变换，生成所述机械信号的短时傅里叶时频谱；变换模块还用于基于传感器获取轴承部位的所述机械信号；对所述机械信号进行短时傅里叶变换以生成所述短时傅里叶变换时频谱。
时间模块用于基于短时傅里叶变换时频谱生成所述机械信号在各个时频位置的群时延；时间模块还用于基于所述短时傅里叶变换时频谱中的相位信息生成所述机械信号在各个时频位置的群时延。
挤压模块用于基于短时傅里叶变换和调制高阶群时延生成调制高阶水平挤压变换的时频谱；挤压模块还用于基于短时傅里叶变换和调制高阶群时延生成时频系数；基于所述时频系数生成高阶水平挤压变换时频谱。
诊断模块用于基于高阶水平挤压变换的时频谱进行机械故障诊断。诊断模块还用于基于高阶水平挤压变换时频谱确定目标脉冲；基于所述目标脉冲的时频特征确定机械故障。
根据本发明的机械故障诊断装置，通过对机械信号进行调制，然后进行短时傅里叶变换，生成所述机械信号的短时傅里叶时频谱；基于短时傅里叶变换时频谱生成所述机械信号在各个时频位置的群时延；基于短时傅里叶变换和调制高阶群时延生成高阶水平挤压变换的时频谱；基于高阶水平挤压变换的时频谱进行机械故障诊断的方式，能够更加精细地刻画信号的非平稳特性，提取机械故障信号的脉冲特征，从而高效地识别机械故障。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种调制高阶水平挤压变换机械故障诊断方法，其特征在于，包括：
S10，获取信号，并对信号进行调制，得到调制信号
Figure DEST_PATH_IMAGE002
，其中表示信号频率，
Figure DEST_PATH_IMAGE004
表示调制频率；
S20，对调制信号
Figure 527096DEST_PATH_IMAGE002
进行短时傅里叶变换，根据短时傅里叶变换时频值计算调制高阶群时延；
S30，基于短时傅里叶变换和调制高阶群时延生成调制高阶水平挤压变换；
S40，基于调制高阶水平挤压变换的时频谱进行机械故障诊断。
2.根据权利要求1所述的一种调制高阶水平挤压变换机械故障诊断方法，其特征在于，对信号进行调制，调制算子定义为
Figure DEST_PATH_IMAGE006
，式中
Figure DEST_PATH_IMAGE008
为虚数单位，
Figure DEST_PATH_IMAGE010
表示信号的相位，
Figure DEST_PATH_IMAGE012
是相位调节算子，
Figure DEST_PATH_IMAGE014
表示对相位调节算子求关于调制频率
Figure 745369DEST_PATH_IMAGE004
的一阶偏导数，通过相位调制得到调制信号
Figure DEST_PATH_IMAGE016
，式中
Figure DEST_PATH_IMAGE018
表示信号振幅，
Figure DEST_PATH_IMAGE020
表示时间。
3.根据权利要求2所述的一种调制高阶水平挤压变换机械故障诊断方法，其特征在于，对所述调制信号
Figure 538488DEST_PATH_IMAGE002
进行短时傅里叶变换，得到时频值
Figure DEST_PATH_IMAGE022
，其中
Figure DEST_PATH_IMAGE024
表示窗函数
Figure DEST_PATH_IMAGE026
，再利用所述时频值计算调制高阶群时延，计算过程如下：
首先，求得
Figure DEST_PATH_IMAGE027
关于时间的
Figure DEST_PATH_IMAGE029
阶偏导数为：
Figure DEST_PATH_IMAGE031
，
再次，根据上式对频率
Figure DEST_PATH_IMAGE033
求偏导数得到：
Figure DEST_PATH_IMAGE035
，
式中，
Figure DEST_PATH_IMAGE037
；然后，利用上式构造高阶方阵
Figure DEST_PATH_IMAGE039
与：
Figure DEST_PATH_IMAGE041
，
Figure DEST_PATH_IMAGE043
，
式中,
Figure DEST_PATH_IMAGE045
；
Figure DEST_PATH_IMAGE047
为窗函数
Figure DEST_PATH_IMAGE048
中的调节因子；
最后，利用高阶方阵
Figure 387493DEST_PATH_IMAGE039
与
Figure DEST_PATH_IMAGE050
计算调制高阶群时延：
Figure DEST_PATH_IMAGE052
,
式中，
Figure DEST_PATH_IMAGE054
表示取复数的实部。
4.根据权利要求3所述的一种调制高阶水平挤压变换机械故障诊断方法，其特征在于，基于短时傅里叶变换和调制高阶群时延得到调制高阶水平挤压变换
Figure DEST_PATH_IMAGE056
：
Figure DEST_PATH_IMAGE058
,
式中，
Figure DEST_PATH_IMAGE060
为一时间变量，且为与t取值无关的不同时间变量，函数
Figure DEST_PATH_IMAGE062
为Dirac函数,
Figure DEST_PATH_IMAGE064
为阈值。
5.根据权利要求1-4任意一项所述的一种调制高阶水平挤压变换机械故障诊断方法，其特征在于，基于调制高阶水平挤压变换时频谱进行机械故障诊断，包括：
基于调制高阶水平挤压变换时频谱确定目标脉冲；
基于所述目标脉冲的时频特征确定机械故障。
6.一种调制高阶水平挤压变换机械故障诊断方法的装置，其特征在于，包括：
变换模块，用于对机械信号进行调制，得到调制信号
Figure 762718DEST_PATH_IMAGE002
，然后进行短时傅里叶变换，生成所述机械信号的短时傅里叶时频谱；
时间模块，用于基于短时傅里叶变换时频谱生成所述机械信号在各个时频位置的群时延；
挤压模块，用于基于短时傅里叶变换和调制高阶群时延生成调制高阶水平挤压变换；
诊断模块，用于基于调制高阶水平挤压变换的时频谱进行机械故障诊断。
7.根据权利要求6所述的一种调制高阶水平挤压变换机械故障诊断方法的装置，其特征在于，对机械信号进行调制，调制算子定义为
Figure DEST_PATH_IMAGE066
，通过相位调制得到调制信号
Figure DEST_PATH_IMAGE067
。
8.根据权利要求6所述的一种调制高阶水平挤压变换机械故障诊断方法的装置，其特征在于，对所述调制信号
Figure 381918DEST_PATH_IMAGE002
进行短时傅里叶变换，得到时频值
Figure DEST_PATH_IMAGE068
，再利用所述时频值计算调制高阶群时延，计算过程如下：
首先，求得关于时间
Figure DEST_PATH_IMAGE070
的
Figure 524317DEST_PATH_IMAGE029
阶偏导数为：
Figure DEST_PATH_IMAGE071
，
再次，根据上式对频率
Figure 254376DEST_PATH_IMAGE033
求偏导数得到：
Figure DEST_PATH_IMAGE072
，
然后，利用上式构造高阶方阵
Figure 533042DEST_PATH_IMAGE039
与
Figure 791985DEST_PATH_IMAGE050
：
Figure DEST_PATH_IMAGE073
，
Figure DEST_PATH_IMAGE074
，
最后，再利用高阶方阵
Figure 156101DEST_PATH_IMAGE039
与
Figure 158692DEST_PATH_IMAGE050
计算调制高阶群时延
Figure DEST_PATH_IMAGE076
：
9.根据权利要求8所述的一种调制高阶水平挤压变换机械故障诊断方法的装置，其特征在于，基于短时傅里叶变换和调制高阶群时延得到调制高阶水平挤压变换
Figure DEST_PATH_IMAGE077
：
Figure 825952DEST_PATH_IMAGE058
。
10.根据权利要求6-9任意一项所述的一种调制高阶水平挤压变换机械故障诊断方法的装置，其特征在于，基于调制高阶水平挤压变换时频谱进行机械故障诊断，包括：
计算机械信号的调制高阶水平挤压变换时频谱的时频包络；将所述时频包络中最大振幅对应的频率作为目标频率；获取所述目标频率对应的目标脉冲；将所述目标脉冲的时频特征和多个预设时频特征进行匹配；根据匹配结果确定所述机械的故障。
CN202211231265.3A 2022-10-10 2022-10-10 一种调制高阶水平挤压变换机械故障诊断方法与装置 Active CN115356108B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211231265.3A CN115356108B (zh) 2022-10-10 2022-10-10 一种调制高阶水平挤压变换机械故障诊断方法与装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211231265.3A CN115356108B (zh) 2022-10-10 2022-10-10 一种调制高阶水平挤压变换机械故障诊断方法与装置

Publications (2)

Publication Number Publication Date
CN115356108A true CN115356108A (zh) 2022-11-18
CN115356108B CN115356108B (zh) 2023-02-10

Family

ID=84008565

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211231265.3A Active CN115356108B (zh) 2022-10-10 2022-10-10 一种调制高阶水平挤压变换机械故障诊断方法与装置

Country Status (1)

Country Link
CN (1) CN115356108B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117009742A (zh) * 2023-06-20 2023-11-07 南方电网调峰调频发电有限公司储能科研院 一种带宽提取域机械信号的故障特征表征方法

Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5646350A (en) * 1996-01-23 1997-07-08 Computational Systems Inc. Monitoring slow speed machinery using integrator and selective correction of frequency spectrum
WO2002065157A2 (en) * 2001-02-14 2002-08-22 The United States Of America, As Represented By The Aministrator Of The National Aeronautics And Space Administration (Nasa) Empirical mode decomposition for analyzing acoustical signals
US20050102116A1 (en) * 2003-11-07 2005-05-12 National Chiao Tung University High-resolution intelligent rotor machine diagnostic system and method
JP2007282006A (ja) * 2006-04-10 2007-10-25 Megachips Lsi Solutions Inc 画像信号処理装置
CN106053080A (zh) * 2016-08-16 2016-10-26 石家庄铁道大学 基于能量切片小波变换的滚动轴承故障特征提取方法
CN108776031A (zh) * 2018-03-21 2018-11-09 南京航空航天大学 一种基于改进的同步挤压变换的旋转机械故障诊断方法
CN109632312A (zh) * 2019-01-22 2019-04-16 北京化工大学 基于多约束非负矩阵分解算法的轴承复合故障诊断方法
CN111289795A (zh) * 2020-02-12 2020-06-16 成都理工大学 高精度高阶时间重排同步挤压变换时频分析方法
CN111474581A (zh) * 2020-05-14 2020-07-31 中国海洋石油集团有限公司 一种基于非线性时间挤压时频变换的瞬态弱信号检测方法
CN111504640A (zh) * 2020-04-30 2020-08-07 电子科技大学 一种加权滑动窗二阶同步压缩s变换轴承故障诊断方法
CN111856562A (zh) * 2020-07-30 2020-10-30 成都理工大学 一种广义高阶同步挤压地震信号时频分解与重构方法
CN112083495A (zh) * 2020-10-15 2020-12-15 中国石油化工股份有限公司 基于变分模态分解的同步压缩小波变换提高分辨率方法
CN112101245A (zh) * 2020-09-18 2020-12-18 丽水市特种设备检测院 基于频域窗函数的短时傅里叶变换机械冲击特征提取方法
CN112668518A (zh) * 2020-12-31 2021-04-16 中国地质大学（武汉） 一种对振动故障信号的vmsst时频分析方法
CN113607446A (zh) * 2021-05-20 2021-11-05 西安交通大学 一种机械设备早期故障诊断方法、系统、设备及存储介质
CN114240915A (zh) * 2021-12-22 2022-03-25 成都阿普奇科技股份有限公司 一种基于图像三维重构的裂纹检测算法
CN114252915A (zh) * 2021-11-03 2022-03-29 成都理工大学 一种基于二阶水平多重同步挤压变换的油气储层识别方法

Patent Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5646350A (en) * 1996-01-23 1997-07-08 Computational Systems Inc. Monitoring slow speed machinery using integrator and selective correction of frequency spectrum
WO2002065157A2 (en) * 2001-02-14 2002-08-22 The United States Of America, As Represented By The Aministrator Of The National Aeronautics And Space Administration (Nasa) Empirical mode decomposition for analyzing acoustical signals
US20050102116A1 (en) * 2003-11-07 2005-05-12 National Chiao Tung University High-resolution intelligent rotor machine diagnostic system and method
JP2007282006A (ja) * 2006-04-10 2007-10-25 Megachips Lsi Solutions Inc 画像信号処理装置
CN106053080A (zh) * 2016-08-16 2016-10-26 石家庄铁道大学 基于能量切片小波变换的滚动轴承故障特征提取方法
CN108776031A (zh) * 2018-03-21 2018-11-09 南京航空航天大学 一种基于改进的同步挤压变换的旋转机械故障诊断方法
CN109632312A (zh) * 2019-01-22 2019-04-16 北京化工大学 基于多约束非负矩阵分解算法的轴承复合故障诊断方法
CN111289795A (zh) * 2020-02-12 2020-06-16 成都理工大学 高精度高阶时间重排同步挤压变换时频分析方法
CN111504640A (zh) * 2020-04-30 2020-08-07 电子科技大学 一种加权滑动窗二阶同步压缩s变换轴承故障诊断方法
CN111474581A (zh) * 2020-05-14 2020-07-31 中国海洋石油集团有限公司 一种基于非线性时间挤压时频变换的瞬态弱信号检测方法
CN111856562A (zh) * 2020-07-30 2020-10-30 成都理工大学 一种广义高阶同步挤压地震信号时频分解与重构方法
CN112101245A (zh) * 2020-09-18 2020-12-18 丽水市特种设备检测院 基于频域窗函数的短时傅里叶变换机械冲击特征提取方法
CN112083495A (zh) * 2020-10-15 2020-12-15 中国石油化工股份有限公司 基于变分模态分解的同步压缩小波变换提高分辨率方法
CN112668518A (zh) * 2020-12-31 2021-04-16 中国地质大学（武汉） 一种对振动故障信号的vmsst时频分析方法
CN113607446A (zh) * 2021-05-20 2021-11-05 西安交通大学 一种机械设备早期故障诊断方法、系统、设备及存储介质
CN114252915A (zh) * 2021-11-03 2022-03-29 成都理工大学 一种基于二阶水平多重同步挤压变换的油气储层识别方法
CN114240915A (zh) * 2021-12-22 2022-03-25 成都阿普奇科技股份有限公司 一种基于图像三维重构的裂纹检测算法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
包文杰等: "参数化的短时傅里叶变换及齿轮箱故障诊断", 《振动.测试与诊断》 *
张云鹏等: "S变换在滚动轴承故障诊断上的应用", 《应用科技》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117009742A (zh) * 2023-06-20 2023-11-07 南方电网调峰调频发电有限公司储能科研院 一种带宽提取域机械信号的故障特征表征方法

Also Published As

Publication number Publication date
CN115356108B (zh) 2023-02-10

Similar Documents

Publication Publication Date Title
Hu et al. An adaptive and tacholess order analysis method based on enhanced empirical wavelet transform for fault detection of bearings with varying speeds
Delgado-Arredondo et al. Methodology for fault detection in induction motors via sound and vibration signals
CN108168891B (zh) 滚动轴承微弱故障信号特征的提取方法及设备
CN110617964A (zh) 用于滚动轴承故障诊断的同步压缩变换阶比分析法
CN109668733B (zh) 变分非线性模式分解变转速轴承故障诊断方法
CN109883706B (zh) 一种滚动轴承局部损伤微弱故障特征提取方法
CN110044610A (zh) 齿轮故障诊断方法
CN111089726A (zh) 一种基于最优维数奇异谱分解的滚动轴承故障诊断方法
CN110763462A (zh) 一种基于同步压缩算子的时变振动信号故障诊断方法
CN109029999B (zh) 基于增强调制双谱分析的滚动轴承故障诊断方法
CN114004256A (zh) 基于数字孪生体的制造装备主轴承的故障诊断方法
CN110987438A (zh) 水轮发电机变转速过程周期性振动冲击信号检测的方法
Mishra et al. Multi-domain bearing fault diagnosis using support vector machine
CN115356108B (zh) 一种调制高阶水平挤压变换机械故障诊断方法与装置
CN108398260B (zh) 基于混合概率方法的齿轮箱瞬时角速度的快速评估方法
CN114486263B (zh) 一种旋转机械滚动轴承振动信号降噪解调方法
CN111307426A (zh) 一种基于FrFT-EWT原理的旋转机械故障特征提取方法
CN113537044B (zh) 基于STFT与改进DenseNet的航空发动机故障诊断方法
CN116358873B (zh) 滚动轴承故障诊断方法、装置、设备及介质
CN102680080B (zh) 一种基于改进的自适应形态滤波的非稳态信号检测方法
CN116361733A (zh) 一种故障诊断方法、装置、系统以及存储介质
CN116256158A (zh) 一种基于深度信号分离的旋转机械瞬时相位自适应提取方法
CN116625681A (zh) 一种基于短时傅里叶变换的谱幅值调制滚动轴承故障诊断方法
CN113138337B (zh) 基于状态划分与频段同步校正的三相电机故障诊断方法
CN112781723B (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
CP02 Change in the address of a patent holder
CP02 Change in the address of a patent holder

Address after: No. 2, Floor 31, Building 1, No. 50, Aerospace Road, Chenghua District, Chengdu, Sichuan 610052

Patentee after: CHENGDU APUQI TECHNOLOGY CO.,LTD.

Patentee after: Suzhou apqi Internet of things Technology Co.,Ltd.

Patentee after: Suzhou Qironggu Technology Co.,Ltd.

Address before: 610052 building B, No.6 Chengjia Road, Longtan Industrial Park, Chenghua District, Chengdu City, Sichuan Province

Patentee before: CHENGDU APUQI TECHNOLOGY CO.,LTD.

Patentee before: Suzhou apqi Internet of things Technology Co.,Ltd.

Patentee before: Suzhou Qironggu Technology Co.,Ltd.