基于最大平均窗峭度盲解卷积的机械故障诊断方法及系统技术方案

本发明专利技术公开了基于最大平均窗峭度盲解卷积的机械故障诊断方法及系统，涉及机械故障诊断技术领域，包括以下步骤：S1，输入测量信号，随机初始化滤波器系数；S2，求解滤波信号；S3，对滤波信号进行加窗实现信号均等分割，求得每段滤波信号的峭度，并求其平均值得到该滤波信号的平均窗峭度；S4，更新滤波器系数，获得新的滤波器；S5，重复步骤S2

【技术实现步骤摘要】
基于最大平均窗峭度盲解卷积的机械故障诊断方法及系统


[0001]本专利技术涉及机械故障诊断
，具体为基于最大平均窗峭度盲解卷积的机械故障诊断方法及系统。

技术介绍

[0002]在基于振动信号分析的旋转机械健康监测和故障诊断中，由于存在多个振动源以及受复杂传递路径和强噪声影响，测量振动信号通常成分复杂，可能包含高斯白噪声以及谐波分量、强冲击干扰等非高斯干扰成分，导致周期性故障信号十分微弱。
[0003]为解决最小熵解卷积(Minimum Entropy Deconvolution,MED)方法在用于故障诊断时倾向于恢复少量主导冲击而非周期性故障冲击序列的问题，提出了最大平均窗峭度盲解卷积方法。方法定义了一种新盲解卷积指标——平均窗峭度，该指标通过对信号进行加窗均等分割计算各等分段峭度的均值，以抑制因少量主导冲击造成信号峭度过大的问题。方法能实现在高斯白噪声、强冲击干扰下均能准确恢复周期性故障冲击序列，提取故障特征信息。齿轮植入故障案例分析结果表明方法可以实现齿轮微弱周期性故障冲击序列的增强检测，同时与MED等传统盲解卷积方法的对比研究进一步突显了所提方法在故障诊断中的优势。

技术实现思路

[0004]本专利技术的目的在于克服现有技术的不足，提供基于最大平均窗峭度盲解卷积的机械故障诊断方法及系统。
[0005]本专利技术的目的是通过以下技术方案来实现的：
[0006]基于最大平均窗峭度盲解卷积的机械故障诊断方法及系统，包括以下步骤：
[0007]S1，输入测量信号，随机初始化滤波器系数，执行步骤S2；
[0008]S2，求解滤波信号，执行步骤S3；
[0009]S3，对将滤波信号进行加窗实现信号均等分割，求得每段滤波信号的峭度，并求其平均值得到该滤波信号的平均窗峭度，执行步骤S4；
[0010]S4，更新滤波器系数，获得新的滤波器，执行步骤S5；
[0011]S5，重复步骤S2
‑
S4，使得滤波信号的平均窗峭度达到最大，执行步骤S6；
[0012]S6，选择滤波信号的平均窗峭度达到最大时对应的滤波器作为最优滤波器，执行步骤S7；
[0013]S7，对滤波信号进行时域分析，据此准确诊断齿轮故障。
[0014]进一步的，所述步骤S3中，平均窗峭度的获得的表达式为：
[0015][0016]式中，AK表示滤波信号的平均窗峭度，M为分割段数，y
m
表示第m段信号，K(y
m
)为其
对应的峭度，m取1到M。
[0017]进一步的，所述K(y
m
)的计算公式为：
[0018][0019]式中：<
·
>表示时域平均算子。
[0020]进一步的，所述步骤S5中的重复步骤S2
‑
S4，使得滤波信号的平均窗峭度达到最大，选择最大平均窗峭度对应的滤波器参数，描述此过程公式如下：
[0021][0022]式中，为当前滤波器系数值，f为滤波器。
[0023]进一步的，所述步骤S4中，新的滤波器的系数获得的公式为：
[0024][0025]式中，X0为一矩阵，N为信号长度，y为滤波信号。
[0026]进一步的，所述步骤4中新的滤波器的系数获得的公式中，所述矩阵式中，L为滤波器长度。
[0027]进一步的，所述测量信号的卷积形式为：
[0028]x＝e*h
e
+n*h
n
[0029]式中，*代表卷积运算，x为测量信号，e为周期性故障信号，n为干扰成分，h
e
、h
n
分别为e和n对应的传递函数。
[0030]进一步的，所述干扰成分包括:杂乱冲击干扰、非高斯谐波分量、高斯白噪声。
[0031]本专利技术的有益效果是：
[0032]为解决MED倾向于恢复主导冲击的问题，本专利技术提出了一种新的盲解卷积指标——平均窗峭度，该指标通过对信号进行加窗均等分割计算各等分段峭度的均值，以抑制因少量主导冲击造成信号峭度过大的问题。基于此提出了一种新的盲解卷积方法——最大平均窗峭度盲解卷积方法，该方法解决了传统盲解卷积方法倾向于恢复少量主导冲击而非周期性故障冲击的问题，方法能实现在高斯白噪声、强冲击干扰下均能准确恢复周期性故障冲击序列，提取故障特征信息。
附图说明
[0033]图1为本专利技术的流程示意图；
[0034]图2为本专利技术的齿轮植入故障测量信号；
[0035]图3分别为MED方法、MEDA方法、OMEDA方法、本专利技术方法得到的故障齿轮测量信号
的滤波信号图(a
‑
MED方法、b
‑
MEDA方法、c
‑
OMEDA方法、d
‑
本专利技术方法)。
具体实施方式
[0036]下面结合本专利技术的附图1～3，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术的一部分实施例，而不是全部的实施例，基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施。
[0037]在本专利技术的描述中，需要理解的是，术语“逆时针”、“顺时针”“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本专利技术，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本专利技术的限制。
[0038]基于最大平均窗峭度盲解卷积的机械故障诊断方法及系统，包括以下步骤：
[0039]S1，输入测量信号，随机初始化滤波器系数，执行步骤S2；
[0040]S2，求解滤波信号，执行步骤S3；
[0041]S3，对将滤波信号进行加窗实现信号均等分割，求得每段滤波信号的峭度，并求其平均值得到该滤波信号的平均窗峭度，执行步骤S4；
[0042]S4，更新滤波器系数，获得新的滤波器，执行步骤S5；
[0043]S5，重复步骤S2
‑
S4，使得滤波信号的平均窗峭度达到最大，执行步骤S6；
[0044]S6，选择滤波信号的平均窗峭度达到最大时对应的滤波器作为最优滤波器，执行步骤S7；
[0045]S7，对滤波信号进行时域分析，据此准确诊断齿轮故障。
[0046]实施例：
[0047]本实施例进行了齿轮植入故障案例分析以验证方法对实际信号分析的有效性。测量信号采自某机械故障试验台，该故障试验台主要由同步皮带轮、三相交流变频电机、传动轴支座、连轴器、磁粉扭力器、齿轮箱组成。在齿轮箱(减速比55/75)故障模拟实验中，在小齿轮某齿根处通过电火花加工植入了局部点蚀故障，并将加速度计安装在齿轮箱输入轴顶部的壳体上。信号采样频本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于最大平均窗峭度盲解卷积的机械故障诊断方法及系统，其特征在于，包括以下步骤：S1，输入测量信号，随机初始化滤波器系数，执行步骤S2；S2，求解滤波信号，执行步骤S3；S3，对将滤波信号进行加窗实现信号均等分割，求得每段滤波信号的峭度，并求其平均值得到该滤波信号的平均窗峭度，执行步骤S4；S4，更新滤波器系数，获得新的滤波器，执行步骤S5；S5，重复步骤S2
‑
S4，使得滤波信号的平均窗峭度达到最大，执行步骤S6；S6，选择滤波信号的平均窗峭度达到最大时对应的滤波器作为最优滤波器，执行步骤S7；S7，对滤波信号进行时域分析，据此准确诊断齿轮故障。2.根据权利要求1所述的基于最大平均窗峭度盲解卷积的机械故障诊断方法及系统，其特征在于，所述步骤S3中，平均窗峭度的获得的表达式为：式中，AK表示滤波信号的平均窗峭度，M为分割段数，y
m
表示第m段信号，K(y
m
)为其对应的峭度，m取1到M。3.根据权利要求2所述的基于最大平均窗峭度盲解卷积的机械故障诊断方法及系统，其特征在于，所述K(y
m
)的计算公式为：式中：<
·
>表示时域平均算子。4.根据权力1所...

【专利技术属性】
技术研发人员：张新，赵艺珂，王家序，吴磊，张忠强，何劲峰，
申请(专利权)人：西南交通大学，
类型：发明
国别省市：