基于最优秩非负矩阵分解的行星齿轮箱早期故障诊断方法技术

本发明专利技术提供一种基于最优秩非负矩阵分解的行星齿轮箱早期故障诊断方法。首先，采集行星齿轮箱故障振动信号，并计算振动信号的STFT谱；然后，将原始振动信号频谱作为K

【技术实现步骤摘要】
基于最优秩非负矩阵分解的行星齿轮箱早期故障诊断方法


[0001]本专利技术属于旋转机械故障诊断与信号分析
，涉及到一种旋转机械早期故障诊断方法，具体为一种基于最优秩非负矩阵分解(Nonnegative Matrix Factorization)的行星齿轮箱早期故障诊断方法，可用于行星齿轮箱中齿轮与滚动轴承部件的故障诊断。

技术介绍

[0002]行星齿轮传动因其体积小、传动比大和承载能力强等优点在现代机械工业领域被广泛应用。考虑其结构特点，行星齿轮箱通常在重载、冲击、污染等恶劣工况下运行，极易导致其关键部件如齿轮、轴承发生故障，引发设备停机甚至更为严重的后果。因此，对行星齿轮箱展开状态监测，探索高效的故障诊断方法，在损伤发展早期对故障特征进行准确提取与识别，对保证设备安全运行与减少经济损失来说具有重要意义。
[0003]行星齿轮箱内部的运行情况可以通过振动信息表达出来。齿轮箱早期故障特征微弱，易淹没在强噪声中。对包含齿轮箱故障信息的振动信号进行有效分析和处理是实现行星齿轮箱早期故障准确诊断的有力途径。同时，受到传递路径和强环境噪声的影响，行星齿轮箱振动信号可视作由多个信号源混合，具有复杂的频谱结构。为了消除噪声干扰，从原始振动信号中分离出表征齿轮箱故障模式的特征信息，多种信号分析技术如小波变换、经验模态分解、集合经验模态分解、基函数法、多参数联合估计等被广泛应用。但这些方法都具有各自的不足。小波变换需要预先确定小波基与分解层数，缺乏自适应性。经验模态分解虽然具有很好的自适应性，能够将输入信号自动分解成一系列固有模态函数与残差之和。但针对具有微弱特征信息的齿轮箱振动信号来说，分解极易产生模态混叠。集合经验模态分解虽然在一定程度上改善了分解的模态混叠问题，但额外引入了附加白噪声的强度与平均次数两个输入参数，且增加了运算量。不同于上述方法，非负矩阵分解作为一种矩阵分解算法，通过对非负矩阵的聚类分析，能够有效提取信号局部特征，且算法简单，具有可解释性和明确的物理意义，已广泛应用于信号处理，模式识别等领域。但在非负矩阵分解算法的实施过程中，秩k的选取一直是令人困扰的难题。针对机械故障信号的非负矩阵分解，现有的k选取手段大多是基于遍历分解后的实验效果或人为经验选取。前者通常十分耗时，后者的分解质量无法得到保证。此外，对分解产生的分量的精准筛选与基于筛选分量如何提取故障特征也影响着非负矩阵分解在微弱信号特征提取方面的应用。

技术实现思路

[0004]为了解决上述问题，本专利技术的目的在于提供一种基于最优秩非负矩阵分解的行星齿轮箱早期故障诊断方法。本专利技术能够摆脱对信号处理技术和诊断经验等先验知识的依赖，克服传统非负矩阵分解算法对分解秩的先验条件，改善强噪声背景下行星齿轮箱故障信息微弱诊断难等问题。
[0005]本专利技术的技术方案为：
[0006]一种基于最优秩非负矩阵分解的行星齿轮箱早期故障诊断方法，包括以下步骤：
[0007]步骤S1：沿行星齿轮箱易于放置传感器的壳体正上方布置加速度传感器，采集行星齿轮箱故障原始振动信号s(t)，并进行快速傅里叶变换得到原始信号频谱S(f)；
[0008]步骤S2：将齿轮箱原始信号频谱S(f)作为输入，初始聚类中心数目k＝2，依次加一，遍历进行K
‑
Means聚类，获得以聚类中心数目k为横坐标，频谱K
‑
Means聚类质量(pkmcq)为纵坐标的折线图。横纵坐标均为无量纲参数；
[0009]其中，新指标pkmcq定义为信号频谱K
‑
means聚类的类内平均距离与类间平均距离之比：
[0010][0011]式中，C
i
表示第i个聚类簇，num(C
i
)表示第i个聚类簇的样本数目，c
i
为C
i
的质心，p为C
i
中的样本点。pkmcq代表了聚类效果的好坏。随着聚类数k由小到大逐渐逼近真实聚类数时，pkmcq值的下降速率也会逐渐减小。当聚类数k超过真实聚类数并逐渐远离时，pkmcq值的下降速率呈现出平缓趋势。以聚类数k为横坐标，pkmcq值为纵坐标绘制k
‑
pkmcq曲线。因此，在k
‑
pkmcq曲线中会首次出现一个明显的拐点，对应的k值即为估计的聚类数；
[0012]步骤S3：对原始振动信号s(t)进行短时傅里叶变换，得到短时傅里叶变换谱图矩阵V，将其作为非负矩阵分解的输入矩阵；
[0013]其中，谱图矩阵V按下述公式得到：
[0014][0015]式中，γ(t)表示短时傅里叶变换窗函数。t表示时间，τ表示时延。随时间t的改变，窗函数会在时间轴上位移，并计算s(t)经窗函数截取部分的傅里叶变换；
[0016]步骤S4：将根据k
‑
pkmcq曲线得到的最佳聚类数目k作为非负矩阵分解最优秩，指导谱图矩阵V的非负矩阵分解，得到满足非负性的两个矩阵W＝[w1,w2,
…
,w
k
]∈R
m
×
k
和H＝[h1,h2,
…
,h
k
]T
∈R
k
×
n
，具体分解过程为：
[0017]V
m
×
n
≈W
m
×
k
×
H
k
×
n
(s.t.W≥0,H≥0)
[0018]式中，W称为基矩阵，H为权重矩阵。
[0019]步骤S5：计算基矩阵W中每个基向量的峭度值，选择峭度值最大的基向量w
o
作为最佳分量进行下一步的处理；
[0020]步骤S6：根据最佳基向量w
o
对原始信号频谱S(f)进行滤波处理，得到滤波信号频谱S
w
(f)，接着通过快速傅里叶反变换得到滤波信号s
w
(t)。滤波具体实现过程表示为：
[0021][0022]步骤S7：对滤波信号s
w
(t)进行希尔伯特变换解调处理，从包络谱中提取反映行星齿轮箱故障的特征信息(齿轮局部故障频率、转频及两者的组合频率)，准确判别故障类型。
[0023]本专利技术的有益效果为：
[0024](1)本专利技术发挥了非负矩阵分解算法在自适应分解信号矩阵方面的优势，通过对
分解得到的矩阵内部分量的有效筛选与处理，解决了行星齿轮箱早期微弱故障诊断难的问题。
[0025](2)本专利技术提出了一种新的信号非负矩阵分解最优秩的估计策略，构造聚类中心数目k与频谱K
‑
Means聚类质量的折线图，有效克服了非负矩阵分解算法对分解秩的先验依赖性，使得非负矩阵分解算法能够适用于强噪声背景下。
[0026](3)本专利技术设计了一种新的自适应滤波方法，借助最大峭度准则自适应选取非负矩阵分解算法处理得到的基矩阵中包本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于最优秩非负矩阵分解的行星齿轮箱早期故障诊断方法，其特征在于，包括以下步骤：步骤S1：沿行星齿轮箱易于放置传感器的壳体正上方布置加速度传感器，采集行星齿轮箱故障原始振动信号s(t)，并进行快速傅里叶变换得到原始信号频谱S(f)；步骤S2：将齿轮箱原始信号频谱S(f)作为输入，初始聚类中心数目k＝2，依次加一，遍历进行K
‑
Means聚类，获得以聚类中心数目k为横坐标，频谱K
‑
Means聚类质量pkmcq为纵坐标的折线图；横纵坐标均为无量纲参数；其中，指标pkmcq定义为信号频谱K
‑
means聚类的类内平均距离与类间平均距离之比：式中，C
i
表示第i个聚类簇，num(C
i
)表示第i个聚类簇的样本数目，c
i
为C
i
的质心，p为C
i
中的样本点；pkmcq代表聚类效果的好坏；随着聚类数k由小到大逐渐逼近真实聚类数时，pkmcq值的下降速率也会逐渐减小；当聚类数k超过真实聚类数并逐渐远离时，pkmcq值的下降速率呈现出平缓趋势；以聚类数k为横坐标，pkmcq值为纵坐标绘制k
‑
pkmcq曲线；因此，在k
‑
pkmcq曲线中会首次出现一个明显的拐点，对应的k值即为估计的聚类数；步骤S3：对原始振动信号s(t)进行短时傅里叶变换，得到短时傅里叶变换谱图矩阵V，将其作为非负矩阵分解的输入矩阵；其中，谱图矩阵V按下述公式得到：式中，γ(t)表示短时傅里叶变换窗函数；t表示时...

【专利技术属性】
技术研发人员：李宏坤，杨晨，陈钧，曹顺心，张孔亮，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：