CN115791159A - 一种基于振动信号的齿轮故障在线实时诊断方法 - Google Patents
一种基于振动信号的齿轮故障在线实时诊断方法 Download PDFInfo
- Publication number
- CN115791159A CN115791159A CN202211669935.XA CN202211669935A CN115791159A CN 115791159 A CN115791159 A CN 115791159A CN 202211669935 A CN202211669935 A CN 202211669935A CN 115791159 A CN115791159 A CN 115791159A
- Authority
- CN
- China
- Prior art keywords
- gear
- fault
- frequency
- signal
- vibration 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.)
- Pending
Links
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明的目的在于提供一种基于振动信号的齿轮故障在线实时诊断方法,包括如下步骤:齿轮系统健康监测,故障类型判断,故障位置判断。针对齿轮断齿、磨损等典型故障,本发明基于齿轮箱体的振动信号,可以实现齿轮传动系统的健康状态在线监测,进而进行故障类型诊断及发生位置确定。本发明提取振动信号的峭度KT作为健康监测特征参数,基于KT的正态分布特点,应用3σ准则对参数的阈值进行了优化。基于EEMD分解提取了IMF分量的能量占比pf和总能量熵HEN作为故障类型判断的特征参数,该特征参数随工况变化很小。
Description
技术领域
本发明涉及的是一种齿轮故障诊断方法,具体地说是齿轮箱中齿轮故障在线诊断方法方法。
背景技术
齿轮箱是机械设备中传递动力和运动的通用零部件,由于其具有传动比固定、结构紧凑、传动精度高等优点,已经在各领域被广泛引用。然而因为轮齿承载的力较大,所处工作环境比较恶劣,所以齿轮容易发生故障,齿轮故障导致设备的力传递过程出现误差,进而导致设备故障。
基于振动信号的方法是齿轮故障诊断常用的方法之一。由于故障发生早期及振动传递路径影响,基于常用的时频域分析方法提取的故障特征不明显。
发明内容
本发明的目的在于提供能通过在时域计算振动信号的峭度,实现齿轮系统在线实时健康状态监测的一种基于振动信号的齿轮故障在线实时诊断方法。
本发明的目的是这样实现的:
本发明一种基于振动信号的齿轮故障在线实时诊断方法,其特征是:
(1)齿轮系统健康监测:
齿轮系统健康状态特征提取:
当齿轮发生故障时,振动信号会包含冲击成分,峭度可以表征信号中的冲击成分,根据振动信号的峭度特征进行齿轮传动系统的运行状态评估;
对于一段数据长度为N振动信号YN=[y(1)…y(N)]T,其峭度KT计算公式为
特征参数的阈值优化与健康状态评估:
分别在不同转速和负载下采集健康、断齿及磨损状态下的振动信号,生成计算样本,提取特征参数KT,经检验在健康、断齿及磨损三种状态下,提取的特征参数KT服从正态分布,表示为
KT~N(μ,σ2)
其中,μ为KT的均值,σ为KT的标准差,σ2为KT的方差;
采用健康状态下的KT作为齿轮健康状态评估的特征参数;
采用3σ准则对KT的阈值进行优化:即特征参数KT的阈值范围为[μ-3σ,μ-3σ];
在线健康状态评估方法:实时测取齿轮箱振动信号,间隔N个数据长度,计算振动信号的峭度特征KT,根据KT阈值范围[μ-3σ,μ-3σ],来判断齿轮系统是否发生故障;
(2)故障类型判断:
齿轮发生故障时会影响到齿轮箱振动的内部激励,从而使振动信号频域特性发生变化,由于齿轮箱振动信号属于非平稳的随机信号,基于单纯的频域分析方法无法得到准确的特征信息,而且提取到的频域特征会受到工况的影响,采用EEMD分解对振动信号进行时频分析,提取两个特征参数作为区分断齿、磨损故障类型的判断依据,在此基础上,应用SVM方法进行故障类型的判断;
(3)故障位置判断:
通过谱峭度及包络解调方法,提取与转频有关的信息来确定齿轮故障发生的位置。
本发明还可以包括:
1、步骤(2)中,针对一对齿轮副系统,包括主动齿和从动齿,其齿数分别为z1、z2,如果主动齿转速n1,则主动齿转频f1为
则,从动齿转频f2为
啮合频率fGMF为
fGMF=f1×z1=f2×z2
对于一段数据长度为N振动信号YN=[y(1)…y(N)]T,进行EEMD分解,得到n个内禀模态分量IMFi(k)(i=1,…,n;k=1,…,N),每个IMFi(k)的能量为Ei
总能量E为:
计算IMFi(k)的能量Ei占总能量E的比值为pi
为了描述振动信号的能量在频带内的分布变化情况,提取所有IMFi(k)分量的总能量熵HEN为
通过EEMD分解选取二个特征参数,一是啮合频率所在频段IMFi(k)的能量占比,定义为pf;二是所有IMFi(k)分量的总能量熵HEN;
针对特征参数pf、HEN,应用“一对一”的二分类方法,设计SVM分类器,对断齿与磨损故障类型进行分类。
2、步骤(3)具体为:
信号预处理:
对信号进行基于线性预测的预处理,去掉振动信号中的线性部分,得到包含故障信息的预测残差信号:
其中,k=1,…,N,a=[a1 a2…am]为AR模型系数;m为模型阶数;
通过BIC方法优化模型的阶数m,并利用最小二乘法确定模型系数a,则模型在k时刻的预测残差信号ε(k)为
确定冲击频段:
通过对ε(k)进行谱峭度分析确定冲击频段:
首先对ε(k)做Wold-Cramér分解,定义白噪声激励信号为X(k),对应的谱增量为dX(f),其中f为频率变量,则ε(k)的Wold-Cramer分解表示为
其中H(k,f)为信号ε(k)在频率f处的复包络;
为了估计ε(k)的能量强度,根据复包络H(k,f)计算ε(k)的2阶矩累积量S2ε(k,f)为
S2ε(k,f)=|H(k,f)|2
再对2阶矩累积量S2ε(k,f)进行时域平均,得到ε(k)的2阶谱累积量S2ε(f)
同理计算得到ε(k)的4阶谱累积量S4ε(f)为
根据以上两个式子,得到ε(k)随频率f变化的谱峭度Ks(f)为
根据预测残差信号的谱峭度,找到最大谱峭度所对应的频率fmax,则fmax附近所对应的频段即为冲击频段;
包络分析:
在冲击频段内,对预测误差信号ε(k)进行带通滤波,得到新的信号g(k)(k=1,…,N),然后对g(k)进行包络谱分析,得到包络谱g(f),确定包络谱中的最大值gmax(fb)所对应频率fb,作为故障位置诊断的特征参数;
故障位置诊断方法:根据特征参数fb,确定fb与主/从动齿转频及倍频之间的相关性,来确定故障发生的位置是主动齿或者从动齿。
本发明的优势在于:
1、针对齿轮断齿、磨损等典型故障,本发明基于齿轮箱体的振动信号,提出了“健康监测-故障类型判断-故障位置判断”的在线实时诊断方法,可以实现齿轮传动系统的健康状态在线监测,进而进行故障类型诊断及发生位置确定。
2、本发明提取振动信号的峭度KT作为健康监测特征参数,基于KT的正态分布特点,应用3σ准则对参数的阈值进行了优化。基于EEMD分解提取了IMF分量的能量占比pf和总能量熵HEN作为故障类型判断的特征参数,该特征参数随工况变化很小。
3、由于故障发生早期及振动传递路径影响,基于常用的时频域分析方法提取的故障特征不明显,本发明提出了基于谱峭度及包络谱分析的齿轮故障定位方法。
附图说明
图1为不同状态下峭度KT的数据分布图;
图2为健康状态下峭度KT的数据分布图;
图3为3000rpm,2N·m工况下齿轮振动信号的EEMD分解图;
图4为基于SVM的故障类型分类结果示意图;
图5为本发明流程图;
图6为断齿情况下的谱峭度;
图7为断齿情况下的包络谱。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
结合图1-7,本发明针对断齿、磨损故障,基于齿轮箱体轴向振动信号,设计了“健康监测-故障类型判断-故障位置判断”的齿轮故障在线诊断方法如下:
1.齿轮系统健康监测
1.1齿轮系统健康状态特征提取
当齿轮发生故障时,振动信号会包含一定的冲击成分。峭度可以有效的表征信号中的冲击成分,而且由于峭度指标是无量纲数,所以其结果不受齿轮运行工况的影响。本专利提出根据振动信号的峭度特征进行齿轮传动系统的运行状态评估。
对于一段数据长度为N振动信号YN=[y(1)…y(N)]T,其峭度KT计算公式为
1.2特征参数的阈值优化与健康状态评估
本专利采用3σ准则对KT的阈值进行优化。
分别在不同转速和负载下采集健康、断齿及磨损状态下的振动信号,生成计算样本,提取特征参数KT。经检验在健康、断齿及磨损三种状态下,本专利针提取的特征参数KT较好的服从正态分布,可以表示为
KT~N(μ,σ2)(2)
其中,μ为KT的均值,σ为KT的标准差,σ2为KT的方差。
由于故障与健康状态下特征参数KT重合的概率很低,因此,本专利采用健康状态下的KT作为齿轮健康状态评估的特征参数。
KT分布在(μ-3σ,μ+3σ)中的概率为0.9974,落在±3σ以外的概率为0.27%。为准确有效确定KT的阈值,本专利提出了基于3σ准则阈值优化方法,即特征参数KT的阈值范围为[μ-3σ,μ-3σ]。
在线健康状态评估方法:实时测取齿轮箱振动信号,间隔N个数据长度,计算振动信号的峭度特征KT,根据KT阈值范围[μ-3σ,μ-3σ],来判断齿轮系统是否发生故障。
2.故障类型判断
齿轮发生故障时会影响到齿轮箱振动的内部激励,从而使振动信号频域特性发生变化。由于齿轮箱振动信号属于非平稳的随机信号,基于单纯的频域分析方法无法得到准确的特征信息,而且提取到的频域特征会受到工况的影响。本发明采用EEMD分解对振动信号进行时频分析,提取了两个特征参数作为区分断齿、磨损故障类型的判断依据。在此基础上,应用SVM方法进行故障类型的判断。
针对一对齿轮副系统,包括主动齿和从动齿,其齿数分别为z1、z2。如果主动齿转速n1,则主动齿转频f1为
则,从动齿转频f2为
啮合频率fGMF为
fGMF=f1×z1=f2×z2 (5)
对于一段数据长度为N振动信号YN=[y(1)…y(N)]T,进行EEMD分解,得到n个内禀模态分量IMFi(k)(i=1,…,n;k=1,…,N),每个IMFi(k)的能量为Ei
总能量E为:
计算IMFi(k)的能量Ei占总能量E的比值为pi
为了描述振动信号的能量在频带内的分布变化情况,提取所有IMFi(k)分量的总能量熵HEN为
考虑到齿轮故障发生时,主要影响啮合频率处的能量,总能量也将发生变化。通过EEMD分解选取二个特征参数,一是啮合频率所在频段IMFi(k)的能量占比,定义为pf;二是所有IMFi(k)分量的总能量熵HEN。
针对特征参数pf、HEN,本专利应用“一对一”的二分类方法,设计SVM分类器,对断齿与磨损故障类型进行分类。
3.故障位置判断
由于齿轮传动系统包括主动齿和从动齿,断齿、磨损等故障可能发生在不同的位置,有必要进行故障位置判断。齿轮故障将导致振动信号出现在边频带调制现象,调制频率与主动齿或从动齿转频及其倍频有关。由于故障发生早期及振动传递路径影响,基于常用的时频域分析方法提取的故障特征不明显,本专利通过谱峭度及包络解调方法,提取与转频有关的信息来确定齿轮故障发生的位置。具体如下:
(1)信号预处理
齿轮箱的振动信号中包含了大量的噪声以及其他部件的振动成分。线性预测是一种基于信号自回归模型的预处理方法,本专利首先对信号进行基于线性预测的预处理,去掉振动信号中的线性部分,得到包含故障信息的预测残差信号。
其中,k=1,…,N,a=[a1 a2…am]为AR模型系数;m为模型阶数。
通过BIC等方法优化模型的阶数m,并利用最小二乘法确定模型系数a。则模型在k时刻的预测残差信号ε(k)为
(2)冲击频段确定
在故障情况下,将对齿轮结构产生冲击,但是由于故障发生早期及振动传递路径影响,这种冲击成分不明显。为了提取残差信号中的冲击成分,本专利通过对ε(k)进行谱峭度分析确定冲击频段。
首先对ε(k)做Wold-Cramér分解,定义白噪声激励信号为X(k),对应的谱增量为dX(f),其中f为频率变量,则ε(k)的Wold-Cramer分解可以表示为
其中H(k,f)为信号ε(k)在频率f处的复包络。
为了估计ε(k)的能量强度,根据复包络H(k,f)计算ε(k)的2阶矩累积量S2ε(k,f)为
式(13)化简为
S2ε(k,f)=|H(k,f)|2 (14)
再对2阶矩累积量S2ε(k,f)进行时域平均,得到ε(k)的2阶谱累积量S2ε(f)
同理计算得到ε(k)的4阶谱累积量S4ε(f)为
根据式(15)、(16),得到ε(k)随频率f变化的谱峭度Ks(f)为
根据预测残差信号的谱峭度,找到最大谱峭度所对应的频率fmax,则fmax附近所对应的频段即为冲击频段,这里取[fmax-200,fmax+200]Hz。
(3)包络分析
在冲击频段[fmax-200,fmax+200]内,对预测误差信号ε(k)进行带通滤波,得到新的信号g(k)(k=1,…,N)。然后对g(k)进行包络谱分析,得到包络谱g(f),确定包络谱中的最大值gmax(fb)所对应频率fb,作为故障位置诊断的特征参数。
故障位置诊断方法:根据特征参数fb,确定fb与主/从动齿转频及倍频之间的相关性,来确定故障发生的位置是主动齿或者从动齿。
实施例:
本发明提出的特征参数阈值优化与分类器设计方法,是在离线状态下进行。利用试验获取的计算样本提取特征参数KT,基于3σ准则对KT的阈值进行优化。然后,利用计算样本提取特征参数pf、HEN,设计“一对一”SVM分类器并进行训练,使分类器准确的区分磨损与断齿故障。具体步骤如下:
步骤1:计算样本生成
计算样本需要足够的长度,才能反应齿轮的运行状态。考虑时频分析对数据长度的需求,计算样本长度应保证主动齿至少要转动50次的振动数据,即计算样本的长度N可按下式计算
式中,n1为主动齿轮转速;fs为振动信号的采样频率。
设定不同转速和负载下,包括转速880、1200、2000、3000rpm,齿轮负载0、0.4、0.8、1.2、1.6、2.0N·m时,模拟断齿与磨损故障,在不同工况下运行。设定采样频率fs=12.8kHz,采集齿轮箱体的振动数据,获得正常、断齿与磨损情况下的振动信号数据,形成计算样本。
本实施例共生成72组计算样本,其中健康、断齿、磨损状态下的计算样本各24组,每组数据长度为N。
步骤2:特征参数KT提取及阈值范围确定
对于一组振动信号YN=[y(1)…y(N)]T,其峭度KT计算公式为
针对72组计算样本,按式(19)计算健康、断齿、磨损三种状态下的峭度,如图1所示。从图中可以看出,三种状态下的峭度都服从正态分布,健康与故障状态下虽然有重合,但是概率很小。健康下的3σ为0.5877,断齿下的3σ为5.4216,磨损下的与3σ为5.7533,故障与健康状态下3σ值相差较远。因此,本专利采用健康状态下的3σ值作为故障与健康状态的判断依据。
健康状态下KT的特征分布图如图2所示,KT均值μ=2.7484,标准差σ=0.1959。应用3σ准则确定KT的阈值范围为[μ-3σ,μ+3σ]=[2.161,3.336]。
步骤3:特征参数pf、HEN提取
对72组计算样本进行EEMD分解,得到n个内禀模态分量IMFi(k)(i=1,…,n;k=1,…,N)。如图3所示为3000rpm,2N·m工况下健康状态下振动信号的IMFi(i=1,…,9)。每个IMFi的能量为Ei
总能量E为
则计算每个IMFi的能量Ei占总能量E的比值,及能量占比pi
为了描述振动信号的能量在频带内的分布变化情况,计算EEMD分解后所有IMFi分量的总能量熵HEN为
考虑到齿轮故障发生时,主要影响啮合频率处的能量,总能量也将发生变化。通过EEMD分解选取二个特征参数,一是啮合频率所在频段的IMF能量占比,定义为pf;二是所有IMFi分量的总能量熵HEN。
以3000rpm时断齿、磨损故障状态下的振动信息为例,提取的不同负载下的特征参数pf、HEN如表1所示。
表1 3000rpm故障状态下的特征参数
从表1可见,在不同负载下,特征参数pf和HEN变化很小;同时,在磨损状态下的pf和HEN都要明显低于断齿状态下的特征参数,在其他工况下也有同样的规律。因此,pf和HEN基本不受工况变化影响。
步骤4:SVM分类器设计与训练
本专利针对特征参数pf、HEN,采用二分类方法,设计“一对一”SVM分类器。然后,利用故障状态下的48组计算样本,提取48组特征参数pf、HEN,并将其分成40组训练样本和8组测试样本。将其中训练样本作为分类器的输入,进行训练,得到断齿与磨损故障分类器。用8组测试样本进行测试,准确率97.2%,其中一组测试样本的分类结果如图4所示。
基于优化后的健康监测阈值和故障类型SVM分类器,以及特征参数KT、pf和HEN的提取方法,本专利提出齿轮传动系统在线诊断方法,诊断过程如图5所示。
齿轮传动系统在线诊断包括三部分:健康状态监测、故障类型判断和故障位置判断。
具体步骤如下:
步骤5:健康状态监测
设置采样频率fs=12.8kHz,实时采集齿轮箱振动信号,根据转速计算振动信号的数据长度N,截取长度为N的振动信号形成数组YN。
根据式(19)计算YN的峭度KT,根据KT的阈值判断齿轮传动系统的健康状态。如此循环往复,每个N个数据长度判断一次。
步骤6:故障类型判断
当KT超过了健康阈值范围,则判定为齿轮发生故障,进入故障类型判断流程。
首先对振动信号YN进行EEMD分解,得到n个内禀模态分量IMFi(k)(i=1,…,n;k=1,…,N)。根据式(20)-(23)提取到啮合频率所在频段的IMF分量的能量占比pf和所有IMFi分量的总能量熵HEN。
将特征参数pf、HEN输入到步骤4设计的SVM分类器中进行分类,从而判断断齿或磨损故障。
步骤7:故障位置判断
确定故障类型后,需要进一步确定故障发生位置,即发生在主动齿或从动齿。
步骤7.1:信号预处理
其中,k=1,…,N,a=[a1 a2…am]为AR模型系数;m为模型阶数。在本实例中模型阶数取10,利用最小二乘法确定模型系数为a=[1.00001.92462.11620.8979-0.6697-1.3201-1.0974-0.5706-0.0738-0.04910.0456]。则模型在k时刻的预测残差信号ε(k)为
步骤7.2:冲击频段确定
首先对ε(k)做Wold-Cramér分解,定义白噪声激励信号为X(k),对应的谱增量为dX(f),其中f为频率变量,则ε(k)的Wold-Cramer分解可以表示为
其中H(k,f)为信号ε(k)在频率f处的复包络。
为了估计ε(k)的能量强度,根据复包络H(k,f)计算ε(k)的2阶矩累积量S2ε(k,f)为
式(27)化简为
S2ε(k,f)=|H(k,f)|2 (28)
再对2阶矩累积量S2ε(k,f)进行时域平均,得到ε(k)的2阶谱累积量S2ε(f)
同理计算得到ε(k)的4阶谱累积量S4ε(f)为
根据式(29)、(30),得到ε(k)在随频率f变化的谱峭度Ks(f)为
如图6所示为3000rpm,2N·m工况下主动齿断齿状态下的谱峭度。根据该谱峭度图,可以确定最大谱峭度值所对应的频率为fmax=1176Hz,则fmax附近所对应的频段即为冲击频段,本实施例中取冲击频段为[976,1376]Hz。
步骤7.3:包络分析
在冲击频段内,对预测误差信号ε(k)进行带通滤波,得到新的信号g(k)(k=1,…,N)。然后对g(k)进行包络谱分析,得到包络谱g(f),确定包络谱中的最大值gmax(fb)所对应频率fb,作为故障位置诊断的特征参数。
如图7所示为断齿情况下的包络谱图,图(a)包络谱图中的最大值所对应频率fb=50Hz,与转速为3000rpm时主动齿转频f1=50Hz相同,所以可以判断为主动齿发生故障。
综上所述,本专利针对断齿、磨损故障,基于齿轮箱体轴向振动信号,设计了“健康监测-故障类型判断-故障位置判断”的齿轮故障在线诊断方法,该方法可以实现对齿轮的在线实时健康监测,并且可以快速、准确的判断故障类型以及故障发生的位置,保证了齿轮传动系统安全、稳定运行。
Claims (3)
1.一种基于振动信号的齿轮故障在线实时诊断方法,其特征是:
(1)齿轮系统健康监测:
齿轮系统健康状态特征提取:
当齿轮发生故障时,振动信号会包含冲击成分,峭度可以表征信号中的冲击成分,根据振动信号的峭度特征进行齿轮传动系统的运行状态评估;
对于一段数据长度为N振动信号YN=[y(1)…y(N)]T,其峭度KT计算公式为
特征参数的阈值优化与健康状态评估:
分别在不同转速和负载下采集健康、断齿及磨损状态下的振动信号,生成计算样本,提取特征参数KT,经检验在健康、断齿及磨损三种状态下,提取的特征参数KT服从正态分布,表示为KT~N(μ,σ2)
其中,μ为KT的均值,σ为KT的标准差,σ2为KT的方差;
采用健康状态下的KT作为齿轮健康状态评估的特征参数;
采用3σ准则对KT的阈值进行优化:即特征参数KT的阈值范围为[μ-3σ,μ-3σ];
在线健康状态评估方法:实时测取齿轮箱振动信号,间隔N个数据长度,计算振动信号的峭度特征KT,根据KT阈值范围[μ-3σ,μ-3σ],来判断齿轮系统是否发生故障;
(2)故障类型判断:
齿轮发生故障时会影响到齿轮箱振动的内部激励,从而使振动信号频域特性发生变化,由于齿轮箱振动信号属于非平稳的随机信号,基于单纯的频域分析方法无法得到准确的特征信息,而且提取到的频域特征会受到工况的影响,采用EEMD分解对振动信号进行时频分析,提取两个特征参数作为区分断齿、磨损故障类型的判断依据,在此基础上,应用SVM方法进行故障类型的判断;
(3)故障位置判断:
通过谱峭度及包络解调方法,提取与转频有关的信息来确定齿轮故障发生的位置。
2.根据权利要求1所述的一种基于振动信号的齿轮故障在线实时诊断方法,其特征是:步骤(2)中,针对一对齿轮副系统,包括主动齿和从动齿,其齿数分别为z1、z2,如果主动齿转速n1,则主动齿转频f1为
则,从动齿转频f2为
啮合频率fGMF为
fGMF=f1×z1=f2×z2
对于一段数据长度为N振动信号YN=[y(1)…y(N)]T,进行EEMD分解,得到n个内禀模态分量IMFi(k)(i=1,…,n;k=1,…,N),每个IMFi(k)的能量为Ei
总能量E为:
计算IMFi(k)的能量Ei占总能量E的比值为pi
为了描述振动信号的能量在频带内的分布变化情况,提取所有IMFi(k)分量的总能量熵HEN为
通过EEMD分解选取二个特征参数,一是啮合频率所在频段IMFi(k)的能量占比,定义为pf;二是所有IMFi(k)分量的总能量熵HEN;
针对特征参数pf、HEN,应用“一对一”的二分类方法,设计SVM分类器,对断齿与磨损故障类型进行分类。
3.根据权利要求1所述的一种基于振动信号的齿轮故障在线实时诊断方法,其特征是:步骤(3)具体为:
信号预处理:
对信号进行基于线性预测的预处理,去掉振动信号中的线性部分,得到包含故障信息的预测残差信号:
其中,k=1,…,N,a=[a1 a2…am]为AR模型系数;m为模型阶数;
通过BIC方法优化模型的阶数m,并利用最小二乘法确定模型系数a,则模型在k时刻的预测残差信号ε(k)为
确定冲击频段:
通过对ε(k)进行谱峭度分析确定冲击频段:
首先对ε(k)做Wold-Cramér分解,定义白噪声激励信号为X(k),对应的谱增量为dX(f),其中f为频率变量,则ε(k)的Wold-Cramer分解表示为
其中H(k,f)为信号ε(k)在频率f处的复包络;
为了估计ε(k)的能量强度,根据复包络H(k,f)计算ε(k)的2阶矩累积量S2ε(k,f)为
S2ε(k,f)=|H(k,f)|2
再对2阶矩累积量S2ε(k,f)进行时域平均,得到ε(k)的2阶谱累积量S2ε(f)
同理计算得到ε(k)的4阶谱累积量S4ε(f)为
根据以上两个式子,得到ε(k)随频率f变化的谱峭度Ks(f)为
根据预测残差信号的谱峭度,找到最大谱峭度所对应的频率fmax,则fmax附近所对应的频段即为冲击频段;
包络分析:
在冲击频段内,对预测误差信号ε(k)进行带通滤波,得到新的信号g(k)(k=1,…,N),然后对g(k)进行包络谱分析,得到包络谱g(f),确定包络谱中的最大值gmax(fb)所对应频率fb,作为故障位置诊断的特征参数;
故障位置诊断方法:根据特征参数fb,确定fb与主/从动齿转频及倍频之间的相关性,来确定故障发生的位置是主动齿或者从动齿。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211669935.XA CN115791159A (zh) | 2022-12-25 | 2022-12-25 | 一种基于振动信号的齿轮故障在线实时诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211669935.XA CN115791159A (zh) | 2022-12-25 | 2022-12-25 | 一种基于振动信号的齿轮故障在线实时诊断方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115791159A true CN115791159A (zh) | 2023-03-14 |
Family
ID=85426654
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211669935.XA Pending CN115791159A (zh) | 2022-12-25 | 2022-12-25 | 一种基于振动信号的齿轮故障在线实时诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115791159A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117150349A (zh) * | 2023-10-31 | 2023-12-01 | 济南嘉宏科技有限责任公司 | 一种智能化设备基础故障自主定位与量化评估方法及系统 |
-
2022
- 2022-12-25 CN CN202211669935.XA patent/CN115791159A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117150349A (zh) * | 2023-10-31 | 2023-12-01 | 济南嘉宏科技有限责任公司 | 一种智能化设备基础故障自主定位与量化评估方法及系统 |
CN117150349B (zh) * | 2023-10-31 | 2024-02-02 | 济南嘉宏科技有限责任公司 | 一种智能化设备基础故障自主定位与量化评估方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Saufi et al. | Challenges and opportunities of deep learning models for machinery fault detection and diagnosis: A review | |
Wang et al. | Support vector data description for fusion of multiple health indicators for enhancing gearbox fault diagnosis and prognosis | |
CN105716857B (zh) | 一种行星齿轮箱健康状态评估方法 | |
Yang et al. | Vibration feature extraction techniques for fault diagnosis of rotating machinery: a literature survey | |
CN107003665B (zh) | 用于评估连接到马达的旋转机械的状态的方法 | |
CN109374293B (zh) | 一种齿轮故障诊断方法 | |
CN113532848A (zh) | 一种行星齿轮箱故障诊断系统 | |
Ramteke et al. | Automated gear fault detection of micron level wear in bevel gears using variational mode decomposition | |
CN113553930A (zh) | 变压器有载分接开关机械故障诊断方法 | |
CN115791159A (zh) | 一种基于振动信号的齿轮故障在线实时诊断方法 | |
US20240053225A1 (en) | Method and Apparatus for Identifying an Abnormality in Mechanical Apparatus or Mechanical Component | |
Chen et al. | A visualized classification method via t-distributed stochastic neighbor embedding and various diagnostic parameters for planetary gearbox fault identification from raw mechanical data | |
CN109635306A (zh) | 基于小波分解和谱峭度的旋转机械故障诊断方法 | |
CN106198009B (zh) | 一种滚动轴承的emd、谱峭度和平滑迭代包络分析方法 | |
CN112686279B (zh) | 一种基于k均值聚类和证据融合的齿轮箱故障诊断方法 | |
Giurgiutiu et al. | Review of vibration-based helicopters health and usage monitoring methods | |
CN106153333B (zh) | 一种基于小波分解滤波的包络分析方法 | |
Bartelmus et al. | Gearbox vibration signal pre-processing and input values choice for neural network training | |
CN112597607A (zh) | 风力发电机高速轴承的预测性维修方法及系统 | |
CN116499742A (zh) | 一种基于集成学习的旋转机械故障监测方法 | |
CN109580218B (zh) | 一种基于似然学习机的风机齿轮箱状态识别方法 | |
Zhang et al. | A novel hybrid compound fault pattern identification method for gearbox based on NIC, MFDFA and WOASVM | |
CN107490477A (zh) | 基于频谱核密度函数相关性比较的齿轮箱故障诊断方法 | |
Mukane et al. | Labview based implementation of fuzzy logic for vibration analysis to identify machinery faults | |
CN113686577B (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 |