中心频率收敛趋势作用下的故障诊断方法技术

本发明专利技术公开了一种中心频率收敛趋势作用下的故障诊断方法，包括(1)采集旋转机械设备的动态信号x(t)；(2)设置变分模型的初始分解参数；(3)使用设定初始分解参数的变分模型分解动态信号x(t)，在中心频率收敛趋势引导下遍历信号分析频带迭代分解动态信号x(t)，得到优化模态{m1...mn...mN}和相应的中心频率{ω1...ωn...ωN}；(4)搜索故障相关模态mI，以故障相关模态mI的中心频率ωI引导参数优化，提取包含故障信息的最优目标分量

Fault diagnosis method under the effect of convergence trend of central frequency

The invention discloses a fault diagnosis method under the effect of central frequency convergence trend, which includes: (1) collecting dynamic signal x (T) of rotating machinery; (2) setting initial decomposition parameters of variation model; (3) decomposing dynamic signal x (T) by using variation model with initial decomposition parameters, traversing signal analysis frequency band and iteratively decomposing dynamic signal x (T) under the guidance of central frequency convergence trend To the optimal mode {M1... Mn... Mn} and the corresponding center frequency {\u03c9 1... \u03c9 n... \u03c9 n}; (4) search the fault-related mode MI, guide the parameter optimization with the center frequency \u03c9 I of the fault-related mode MI, and extract the optimal target component containing the fault information

【技术实现步骤摘要】
中心频率收敛趋势作用下的故障诊断方法
本专利技术属于机械微弱故障诊断
，涉及一种中心频率收敛趋势作用下的故障诊断方法。
技术介绍
旋转机械设备在工业生产中得到了广泛的应用，机械零部件的状态直接影响到机械设备的运行状态及其安全状况。当机械零部件出现故障时，会产生周期性的瞬态冲击响应，如何对其进行有效提取和准确评估是轴承故障诊断的关键。但是，由于实际运行环境的复杂性，从设备现场采集的动态信号是包含大量噪声的，信号的微弱故障特征往往被噪声所淹没，从而严重影响了故障特征信号的识别。因此，开展机械微弱故障瞬态特征提取与判别具有实际意义。目前已发展了许多机械故障诊断方法，例如时频域分析方法、经验模态分解、局部均值分解等传统微弱故障诊断方法。可是，这些传统方法因为其自身的局限性，如停止准则选取困难或抗噪性能较差等问题，导致其应用范围受限。变分模态分解方法是一种基于变分模型的信号自适应分解方法，具有较强的抗噪性，采用非筛选分解方式分解信号减少传递误差。近年来，逐渐有学者将变分模态分解方法引入到机械信号处理领域，发展出了基于变分模态分解和分类模型的结合的轴承故障诊断方法；利用变分模态分解方法在多级离心泵滚动轴承故障诊断中的应用；扩展了变分模态分解在结构系统识别中的适用性。但，目前利用变分模态分解方法处理机械信号时，很难预知设备原始动态信号中的实际中心频率、模态分量的个数以及很难完整提取出相应目标分量的最佳平衡参数。
技术实现思路
本专利技术实施例提供一种中心频率收敛趋势作用下的故障诊断方法，在变分模态分解方法的基础上，采用中心频率收敛趋势引导的分解方式实现诊断目标设备原始动态信号的智能分解，克服传统变分模态分解方法初始参数设置的难题，能够对采集的设备动态信号自适应地分析，降低了技术人员使用变分模态分解方法进行机械故障诊断的难度。为了解决上述技术问题，本专利技术提供了一种中心频率收敛趋势作用下的故障诊断方法，包括以下步骤，(1)以fs为采样频率采集诊断目标的动态信号x(t)；(2)设置变分模型的初始分解参数：初始中心频率ω0为0、初始中心频率的增长步长Δω为100Hz、初始计步数z为1、平衡参数α为[1000，4000]和模态分量个数K为1；(3)使用设定初始分解参数的变分模型对所述动态信号x(t)进行一次分解，判断中心频率的收敛趋势，在中心频率收敛趋势引导下遍历信号分析频带迭代分解动态信号x(t)，得到优化模态{m1...mn...mN}和相应的中心频率{ω1...ωn...ωN}；(4)从获取的优化模态{m1...mn...mN}中搜索故障相关模态mI，以故障相关模态mI的中心频率ωI引导参数优化，提取包含故障信息的最优目标分量(5)对提取的最优目标分量进行包络分析，根据其包络谱诊断所述旋转机械设备。本专利技术一个较佳实施例中，进一步包括步骤(3)中，使用交替方向乘子法求解变分模型中的约束模型：式中，x(t)为动态信号，*表示卷积算子，表示对时间t求偏导数，δ(t)为狄利克雷分布函数，指数调节项用于平移各分量的频谱；信号x(t)被分解为K个模态分量mk(k＝1，2，3......K)，每个模态分量mk都围绕其中心频率ωk。本专利技术一个较佳实施例中，进一步包括所述动态信号x(t)在中心频率收敛趋势引导下迭代分解过程为：(S31)使用设定初始分解参数的变分模型对所述动态信号x(t)进行一次分解，得到更新后的中心频率ω1；(S32)判断中心频率的收敛趋势e＝ω1-ω0：若收敛趋势e＝ω1-ω0为上升趋势，则输出相应的模态分量作为优化模态mn，相应的中心频率ωn为提取的最优中心频率；若收敛趋势e＝ω1-ω0为下降趋势，则令ω0＝ω0+zΔω，并同时判断是否遍历整个频带，若ω0＝(ω0+zΔω)＜fs/2，返回步骤(S31)，否则停止迭代分解；(S33)以提取的最优中心频率ωn更新初始中心频率ω0，若新的中心频率ω0＜fs/2，返回步骤(S31)，否则停止迭代分解。本专利技术一个较佳实施例中，进一步包括步骤(4)中，从获取的优化模态{m1...mn...mN}中搜索故障相关模态mI时，通过计算优化模态{m1...mn...mN}的基尼指数值来确定。本专利技术一个较佳实施例中，进一步包括步骤(4)中，以故障相关模态mI的中心频率ωI引导参数优化，提取包含故障信息的最优目标分量的过程为：(S51)设置两组初始分解参数：平衡参数α＝α0+Δα、模态分量个数K＝1、初始中心频率为ωI；平衡参数α＝α0-Δα、模态分量个数K＝1、初始中心频率为ωI；其中，Δα为平衡参数α变化的步长；(S52)分别使用步骤(S51)中设置的两组初始分解参数分解原始动态信号x(t)，得到两组模态分量Ur1、Ul1；(S53)分别计算模态分量Ur1、Ul1的基尼指数值Gnir1、Gnil1；(S54)判断Gnir1、Gnil1的大小：若Gnir1＞Gnil1，执行递增平衡参数优化方案；否则，执行递减平衡参数优化方案。本专利技术一个较佳实施例中，进一步包括所述递增平衡参数优化方案包括，(S61)设置分解参数：平衡参数α＝α0+iΔα(i＝2)，模态分量个数K＝1，初始中心频率为ωI；(S62)使用步骤(S61)中设置的分解参数分解原始动态信号x(t)，得到模态分量Uri，并计算模态分量Uri的基尼指数值Gniri；(S63)判断Gniri、Gniri-1的大小，若Gniri＞Gniri-1，则令i＝i+1，并返回步骤(S61)；否则，令本专利技术一个较佳实施例中，进一步包括所述递减平衡参数优化方案包括，(S71)设置分解参数：平衡参数α＝α0-iΔα(i＝2)，模态分量个数K＝1，初始中心频率为ωI；(S72)使用步骤(S71)中设置的分解参数分解原始动态信号x(t)，得到模态分量Uli，并计算模态分量Uli的基尼指数值Gnili；(S73)判断Gnili、Gnili-1的大小，若Gnili＞Gnili-1，则令i＝i+1，并返回步骤(S71)；否则，令本专利技术的有益效果：其一、本专利技术实施例的中心频率收敛趋势作用下的故障诊断方法，在变分模态分解方法的基础上，采用中心频率收敛趋势引导的分解方式实现诊断目标设备原始动态信号的智能分解，克服传统变分模态分解方法初始参数设置的难题，能够对采集的设备动态信号自适应地分析，降低了技术人员使用变分模态分解方法进行机械故障诊断的难度。其二、本专利技术实施例的中心频率收敛趋势作用下的故障诊断方法，在变分模态分解方法的基础上，采用中心频率收敛趋势引导的分解方式，能够加速分解算法的收敛过程，同时避免了现有分解方法中因预设的不恰当分解模态分量个数而导致的模态混叠和虚假分量问题。其三、本专利技术实施例的中心频率收敛趋势作用下的故障诊断方法，使用中心频率引导平衡参数自适应优化，可以使得最终得到分量的带宽最大限度的匹配实际故障分量的带宽，同时减少计算量。附图说明图1是本专利技术实施例中故障诊断方法的流程图；图2是本专利技术实施例中中心频率收敛趋势引导动态信号分解过程的流程图；图3是本专利技术实施例中中心频率引导参数优化提取包含故障信息的最优目标分量的流程图；图4是采集的一组齿轮箱损伤动态信号波形图；图5是使用本专利技术实施例中故障诊断方法智能分解图4动态信号的四个分量波形图；图6是利本文档来自技高网...

【技术保护点】
1.一种中心频率收敛趋势作用下的故障诊断方法，其特征在于：包括以下步骤，(1)以fs为采样频率采集诊断目标的动态信号x(t)；(2)设置变分模型的初始分解参数：初始中心频率ω0为0、初始中心频率的增长步长Δω为100Hz、初始计步数z为1、平衡参数α为[1000，4000]和模态分量个数K为1；(3)使用设定初始分解参数的变分模型对所述动态信号x(t)进行一次分解，判断中心频率的收敛趋势，在中心频率收敛趋势引导下遍历信号分析频带迭代分解动态信号x(t)，得到优化模态{m1...mn...mN}和相应的中心频率{ω1...ωn...ωN}；(4)从获取的优化模态{m1...mn...mN}中搜索故障相关模态mI，以故障相关模态mI的中心频率ωI引导参数优化，提取包含故障信息的最优目标分量

【技术特征摘要】
1.一种中心频率收敛趋势作用下的故障诊断方法，其特征在于：包括以下步骤，(1)以fs为采样频率采集诊断目标的动态信号x(t)；(2)设置变分模型的初始分解参数：初始中心频率ω0为0、初始中心频率的增长步长Δω为100Hz、初始计步数z为1、平衡参数α为[1000，4000]和模态分量个数K为1；(3)使用设定初始分解参数的变分模型对所述动态信号x(t)进行一次分解，判断中心频率的收敛趋势，在中心频率收敛趋势引导下遍历信号分析频带迭代分解动态信号x(t)，得到优化模态{m1...mn...mN}和相应的中心频率{ω1...ωn...ωN}；(4)从获取的优化模态{m1...mn...mN}中搜索故障相关模态mI，以故障相关模态mI的中心频率ωI引导参数优化，提取包含故障信息的最优目标分量(5)对提取的最优目标分量进行包络分析，根据其包络谱诊断所述旋转机械设备。2.如权利要求1所述的中心频率收敛趋势作用下的故障诊断方法，其特征在于：步骤(3)中，使用交替方向乘子法求解变分模型中的约束模型：式中，x(t)为动态信号，*表示卷积算子，表示对时间t求偏导数，δ(t)为狄利克雷分布函数，指数调节项用于平移各分量的频谱；信号x(t)被分解为K个模态分量mk(k＝1，2，3......K)，每个模态分量mk都围绕其中心频率ωk。3.如权利要求1所述的中心频率收敛趋势作用下的故障诊断方法，其特征在于：所述动态信号x(t)在中心频率收敛趋势引导下迭代分解过程为：(S31)使用设定初始分解参数的变分模型对所述动态信号x(t)进行一次分解，得到更新后的中心频率ω1；(S32)判断中心频率的收敛趋势e＝ω1-ω0：若收敛趋势e＝ω1-ω0为上升趋势，则输出相应的模态分量作为优化模态mn，相应的中心频率ωn为提取的最优中心频率；若收敛趋势e＝ω1-ω0为下降趋势，则令ω0＝ω0+zΔω，并同时判断是否遍历整个频带，若ω0＝(ω0+zΔω)＜fs/2，返回步骤(S31)，否则停止迭代分解；(S33)以提取的最优中心频率ωn更新初始中心频率ω0，若新的中心频率ω0＜fs/2，返回步骤(S31)，否则停止迭代分解。4.如...

【专利技术属性】
技术研发人员：江星星，沈长青，周建芹，宋冬淼，郭文军，杜贵府，王俊，石娟娟，黄伟国，朱忠奎，
申请(专利权)人：苏州大学，
类型：发明
国别省市：江苏,32