一种基于VMD的齿轮振动信号降噪及故障诊断方法技术

本发明专利技术涉及一种基于VMD的齿轮振动信号降噪及故障诊断方法，其步骤：采集齿轮在正常运行状态、点蚀故障、磨损故障以及断齿故障下的振动信号xw(n)＝{xw(1),...,xw(N)}，其中，N代表每组数据个数，w代表数据组别，w＝1、2、3、4，分别代表正常运行状态、点蚀故障状态、磨损故障状态、断齿故障状态；对xw(n)＝{xw(1),...,xw(N)}中正常运行状态下的原始振动信号x1(n)进行VMD‑模平方阈值处理；将所有齿轮振动信号xw(n)按照步骤2)进行VMD‑模平方阈值处理，得到处理后的信号

A Noise Reduction and Fault Diagnosis Method of Gear Vibration Signal Based on VMD

The present invention relates to a method for noise reduction and fault diagnosis of gear vibration signal based on VMD. The steps are as follows: acquisition of vibration signal xw(n)={xw(1),...,xw(N)} of gear under normal operation, pitting failure, wear failure and broken tooth failure, where N represents the number of data sets, w represents data groups, w = 1, 2, 3, 4, respectively, representing normal operation status and pitting failure status. The original vibration signal x1(n) in the normal operation state of xw(n)={xw(1),...,xw(N)} is processed by VMD_mode square threshold, and all gear vibration signals xw(n) are processed by VMD_mode square threshold according to step 2 to get the processed signal.

【技术实现步骤摘要】
一种基于VMD的齿轮振动信号降噪及故障诊断方法
本专利技术涉及机械振动信号降噪及故障诊断
，特别是关于一种基于VMD的齿轮振动信号降噪及故障诊断方法。
技术介绍
作为机械传动形式的基本组成部分，齿轮运行情况的好坏直接影响到机械设备的运行状态。由于齿轮故障信号常呈现出强耦合性、非线性、非平稳性的特点，传统降噪与故障诊断方法难以对故障进行有效的诊断，所以有效的信号降噪和故障诊断方法可以显著提高齿轮故障诊断的准确率。
技术实现思路
针对上述问题，本专利技术的目的是提供一种基于VMD的齿轮振动信号降噪及故障诊断方法，该方法能够有效的提取非平稳齿轮运行信号的故障特征，并进行故障诊断，能够对齿轮故障进行有效的识别。为实现上述目的，本专利技术采取以下技术方案：一种基于VMD的齿轮振动信号降噪及故障诊断方法，其包括以下步骤：1)采集齿轮在正常运行状态、点蚀故障、磨损故障以及断齿故障下的振动信号xw(n)＝{xw(1),...,xw(N)}，其中，N代表每组数据个数，w代表数据组别，w＝1、2、3、4，分别代表正常运行状态、点蚀故障状态、磨损故障状态、断齿故障状态；；n为第w组数据中第n个数据点，n∈(1,N)；2)对xw(n)＝{xw(1),...,xw(N)}中正常运行状态下的原始振动信号x1(n)进行VMD-模平方阈值处理；3)将所有齿轮振动信号xw(n)按照步骤2)进行VMD-模平方阈值处理，得到处理后的信号4)提取处理后的信号的每组信号的峭度和均方根值组成特征向量K和R；5)利用PNN进行故障诊断。进一步，所述步骤2)中，VMD-模平方阈值处理步骤如下：2.1)首先对信号x1(n)进行VMD分解；2.2)原始振动信号x1(n)经VMD后分解为k个模态分量IMF1,IMF2,‥‥,IMFk，对每个IMF分量进行模平方阈值处理；2.3)将模平方阈值处理后的分量重构，得到重构后的信号进一步，所述步骤2.1)中，VMD分解的步骤为：2.1.1)令m＝0,初始化第k个模态函数uk记为当前模态函数功率谱的中心ωk记为和m＝0为所对应的傅里叶变换λ1，其中m为迭代次数；λ为所对应的傅里叶变换；k为原始振动信号x1(n)进行VMD分解后得到的模态个数；2.1.2)根据式更新uk；根据式更新ωk；根据式更新λ；其中ui为第i个模态函数，i≠k；为第k个模态的第m-1次迭代的模态函数；相当于第m次迭代时对当前剩余量x1-∑i≠kui的维纳滤波；ω为中心频率；ωk为第k个模态的中心频率；为第k个模态的第m-1次迭代的中心频率；为m次迭代模态函数功率谱的中心；λm为第m-1次迭代所对应的傅里叶变换；λm-1为第m次迭代所对应的傅里叶变换；α为惩罚因子；2.1.3)根据收敛条件判断结果是否满足；e为预先设定的允许误差，若满足则停止迭代，否则返回步骤2.1.2)。进一步，所述步骤2.2)中，对每个IMF分量进行模平方阈值处理的步骤为：2.2.1)对IMFt分量进行小波阈值降噪，首先选择小波分解的分解层数为3层与小波基函数db5，对滚动轴承原始振动信号采用小波分解函数wavedec进行分解；其中，t∈(1,k)；2.2.2)提取小波分解后振动信号的每层的细节系数cd1，cd2，cd3与最后一层的近似系数ca3，运用模平方阈值方法对每层的细节系数进行处理，得到新的细节系数2.2.3)将步骤2.2.2)得到的新的细节系数与最后一层的近似系数ca3组成新的小波分解结构，利用小波重构函数waverec进行小波重构，重构后得到的信号即为模平方阈值降噪后的信号。进一步，所述步骤2.2.2)中，运用模平方阈值方法对每层的细节系数进行处理步骤如下：(1)利用公式调整每一层的阈值，其中λj为每一层的阈值，j为分解尺度，j＝1,2,3，N为信号的长度，cdj,y为第j层第y个细节系数，median()为返回第一层细节系数的中值；经公式调整后得到不同分解尺度的阈值λ1,λ2,λ3；(2)利用得到的阈值对每一层的细节系数进行处理；当细节系数cdj的绝对值小于该层阈值λj时将其置0，当细节系数cdj的绝对值大于或等于该层阈值λj时则将其代入公式得到新的细节系数sign()为符号函数。进一步，所述步骤5)中，PNN故障诊断步骤如下：5.1)利用net＝newpnn(P,T,SPREAD)函数建立PNN神经网络模型，P为输入向量组成的矩阵，T为目标向量组成的矩阵，SPREAD为径向基函数的分布密度；5.2)将提取的特征向量分为训练样本KP,RP和测试样本KT,RT，将训练样本输入已建立的神经网络进行训练；5.3)将测试样本KT,RT输入训练好的神经网络，输出故障分类结果，即可判断齿轮的运行状态，对齿轮进行故障诊断。进一步，所述步骤5.1)中，newpnn算法为：PNN分4层，分别是输入层，模式层，求和层和决策层；输入层单元不参与计算，模式层神经元的个数与各类别训练的样本数之和相等，该模式层各模式单元的输出为：式中，Wi为连接两层之间的权值；δ为平滑系数；求和层神经元通过与同类的各模式层神经元的输出求平均值，计算该样本属于该类别的最大概率，以确定各故障模式的估计概率密度函数，计算为：式中，Xai为输入层神经元向量；d属于某类别的样本数目；决策层利用求和层各模式的估计概率密度，采用Bayes分类规则，选取最大后验概率密度的神经元作为系统输出。本专利技术由于采取以上技术方案，其具有以下优点：1、本专利技术采用了VMD与模平方阈值相结合的信号降噪方法对齿轮的故障振动信号进行处理，有效地减少了齿轮振动信号中的噪声对故障信息的干扰，能有效的提取故障特征，提高故障诊断的效率。2、本专利技术提取降噪后信号的峭度和均方根值并将其组成特征向量，可以更好的反映故障状态，从而获得准确的故障诊断结果。3、本专利技术将PNN应用于故障诊断，并与VMD-模平方阈值结合，使故障诊断的准确率明显提高。附图说明图1是本专利技术的整体流程图。具体实施方式下面结合附图和实施例对本专利技术进行详细的描述。如图1所示，本专利技术提供一种基于VMD(变分模态分解)的齿轮振动信号降噪及故障诊断方法，该方法包括以下步骤：1)采集齿轮在正常运行状态、点蚀故障、磨损故障以及断齿故障下的振动信号xw(n)＝{xw(1),...,xw(N)}，其中，N代表每组数据个数，w代表数据组别，w＝1、2、3、4，分别代表正常运行状态、点蚀故障状态、磨损故障状态、断齿故障状态；n为第w组数据中第n个数据点，n∈(1,N)。2)以xw(n)＝{xw(1),...,xw(N)}中正常运行状态下的一组原始振动信号x1(n)为例，对该原始振动信号x1(n)进行VMD-模平方阈值处理，具体步骤如下：2.1)首先对信号x1(n)进行VMD分解；VMD分解的步骤为：2.1.1)令m＝0,初始化第k个模态函数uk记为当前模态函数功率谱的中心ωk记为和m＝0为所对应的傅里叶变换λ1，其中m为迭代次数；λ为所对应的傅里叶变换；k为原始振动信号x1(n)进行VMD分解后得到的模态个数；2.1.2)根据式更新uk；根据式更新ωk；根据式更新λ；其中ui为第i个模态函数，i≠k；为第k个模态的第m-1次迭代的模态函数；相当于第m次迭代时对当前剩余量x1-∑i≠kui的维纳滤波；ω为中本文档来自技高网...

【技术保护点】
1.一种基于VMD的齿轮振动信号降噪及故障诊断方法，其特征在于，包括以下步骤：1)采集齿轮在正常运行状态、点蚀故障、磨损故障以及断齿故障下的振动信号xw(n)＝{xw(1),...,xw(N)}，其中，N代表每组数据个数，w代表数据组别，w＝1、2、3、4，分别代表正常运行状态、点蚀故障状态、磨损故障状态、断齿故障状态；；n为第w组数据中第n个数据点，n∈(1,N)；2)对xw(n)＝{xw(1),...,xw(N)}中正常运行状态下的原始振动信号x1(n)进行VMD‑模平方阈值处理；3)将所有齿轮振动信号xw(n)按照步骤2)进行VMD‑模平方阈值处理，得到处理后的信号

【技术特征摘要】
1.一种基于VMD的齿轮振动信号降噪及故障诊断方法，其特征在于，包括以下步骤：1)采集齿轮在正常运行状态、点蚀故障、磨损故障以及断齿故障下的振动信号xw(n)＝{xw(1),...,xw(N)}，其中，N代表每组数据个数，w代表数据组别，w＝1、2、3、4，分别代表正常运行状态、点蚀故障状态、磨损故障状态、断齿故障状态；；n为第w组数据中第n个数据点，n∈(1,N)；2)对xw(n)＝{xw(1),...,xw(N)}中正常运行状态下的原始振动信号x1(n)进行VMD-模平方阈值处理；3)将所有齿轮振动信号xw(n)按照步骤2)进行VMD-模平方阈值处理，得到处理后的信号4)提取处理后的信号的每组信号的峭度和均方根值组成特征向量K和R；5)利用PNN进行故障诊断。2.如权利要求1所述的降噪及故障诊断方法，其特征在于：所述步骤2)中，VMD-模平方阈值处理步骤如下：2.1)首先对信号x1(n)进行VMD分解；2.2)原始振动信号x1(n)经VMD后分解为k个模态分量IMF1,IMF2,‥‥,IMFk，对每个IMF分量进行模平方阈值处理；2.3)将模平方阈值处理后的分量重构，得到重构后的信号3.如权利要求2所述的降噪及故障诊断方法，其特征在于：所述步骤2.1)中，VMD分解的步骤为：2.1.1)令m＝0,初始化第k个模态函数uk记为当前模态函数功率谱的中心ωk记为和m＝0为所对应的傅里叶变换λ1，其中m为迭代次数；λ为所对应的傅里叶变换；k为原始振动信号x1(n)进行VMD分解后得到的模态个数；2.1.2)根据式更新uk；根据式更新ωk；根据式更新λ；其中ui为第i个模态函数，i≠k；为第k个模态的第m-1次迭代的模态函数；相当于第m次迭代时对当前剩余量x1-∑i≠kui的维纳滤波；ω为中心频率；ωk为第k个模态的中心频率；为第k个模态的第m-1次迭代的中心频率；为m次迭代模态函数功率谱的中心；λm为第m-1次迭代所对应的傅里叶变换；λm-1为第m次迭代所对应的傅里叶变换；a为惩罚因子；2.1.3)根据收敛条件判断结果是否满足；e为预先设定的允许误差，若满足则停止迭代，否则返回步骤2.1.2)。4.如权利要求2所述的降噪及故障诊断方法，其特征在于：所述步骤2.2)中，对每个IMF分量进行模平方阈值处理的步骤为：2.2.1)对IMFt分量进行小波阈值降噪，首先选择小波分解的分解层数为3层与小波基函数db5，对滚动轴...

【专利技术属性】
技术研发人员：徐小力，刘秀丽，吴国新，蒋章雷，张雪英，
申请(专利权)人：北京信息科技大学，
类型：发明
国别省市：北京,11