基于LOD-ICA的旋转机械故障诊断方法技术

本发明专利技术涉及旋转机械故障信号处理和故障诊断领域，公开了一种基于LOD

【技术实现步骤摘要】
基于LOD
‑
ICA的旋转机械故障诊断方法


[0001]本专利技术涉及旋转机械故障信号处理和故障诊断领域，特别涉及基于盲源分离的旋转机械故障诊断方法。

技术介绍

[0002]旋转机械的故障诊断技术的研究对于提高工作效率甚至社会生产力、生命安全以及国家经济建设都具有重要意义。
[0003]大型旋转机械出现故障的时候，很多情况下无法预先知道故障信号特征，需要采用盲源分离的方法进行诊断。独立分量分析(independent component analysis，ICA)方法在处理源信号的盲源分离问题和在图形信号处理与分析领域应用广泛，而快速ICA方法以其更快的收敛速度要优于ICA方法，但快速ICA方法与ICA方法一样都存在观测数据通道不足的缺陷。在观测通道数不足的情况下，独立分量分析，需要利用其他分解方法构建虚拟通道数据，以实现各独立信号源的分离。

技术实现思路

[0004]本专利技术要解决的技术问题是：通过LOD(Local Oscillatory
‑
Characteristic Decomposition)方法对故障信号进行预处理，为ICA方法提供足够的虚拟通道数据；在拥有足够多观测通道情况下，利用ICA算法解决LOD方法存在的模态混肴问题，以实现信号的盲源分离。
[0005]本专利技术采用的技术方案是：基于LOD
‑
ICA的旋转机械故障诊断方法，通过基于噪声辅助分析的LOD方法对故障振动信号进行预处理，为快速ICA方法提供足够的虚拟通道数据；其次，在拥有足够多观测通道情况下，以快速ICA算法可以解决基于白噪声辅助分析的LOD方法存在的模态混肴问题，以实现信号的盲源分离；最后，成功提取出混合信号的故障特征。
[0006]LOD方法基于信号的局部特征时标，并且通过微分、坐标域变换和分段线性变换等操作，将信号分解为一系列单振荡分量MOCs，其瞬时频率具有物理意义和残基。该方法计算效率很高，更适合实时在线处理和数据分析。
[0007]进一步的，本专利技术基于快速ICA方法和基于白噪声辅助分析的LOD方法的结合，充分发挥两种方法彼此处理信号的优势，增强了混合故障振动信号特征提取的精确性。
[0008]进一步的，本专利技术实现故障振动信号特征提取的具体步骤包括：
[0009]步骤1：保存采集的旋转机械故障混合信号；
[0010]步骤2：建立任意两相邻极值点区间的分段线性函数，通过坐标域变换将信号从原始数据域坐变换为锯齿域坐标；
[0011]步骤3：对一阶微分信号添加白噪声，求得锯齿域的初始均值函数，进而求得高频波动函数，再通过坐标域变换运算，将高频波动函数转换为原始域高频函数；
[0012]步骤4：鉴别分离出来的信号单分量c1(t)是否含有不同的波动特征，即其瞬时频
(t)，因此判定原始信号x(t)的第1个MOC分量MOC1(t)为信号c1(t)；实际情况下，假使分离出来的信号分量c1(t)是含有不同的波动特征，即其瞬时频率不具有物理意义的一个单分量信号，则循环运算以上步骤m次，最终分离出来的第一个瞬时频率具有物理意义的单分量信号c
m
将被判定为原始信号中的第1个MOC分量MOC1(t)。
[0031]其中，所述步骤4的具体过程为：通过剥离出原始信号x(t)里的单分量信号MOC1(t)来创造出一个全新的信号函数r1(t)，将函数r1(t)视为新的原始信号重新带入LOD算法中，通过以上步骤的循环运算，n次后会得到有且仅有一个单调函数r
n
(t)，最终通过LOD方法对信号的分解，会以n个MOC分量及残差r
n
(t)相加总和构成原始信号x(t)的形式存在。
[0032]进一步的，所述进行实现故障振动信号特征提取关键步骤在于，利用快速ICA算法，依据负熵来建立起目标函数，通过经典的牛顿迭代进行优化，实现混合信号盲源分离。
[0033]本专利技术的有益效果是：对于LMD算法，使用移动平均值来导出局部均值函数和包络估计函数会导致大量计算。相比于LMD算法，LOD算法具有更高的数据处理速度，更快的计算效率，更加精确的分解结果。基于EMD方法而言，其自身计算量，甚至存在的末端效应，过冲和下冲，以及由希尔伯特变换产生的无法解释的负频率等问题并没有在ICA算法的加入后而取得良好的改善效果。所以根据上述内容，LOD算法在信号处理方面较之EMD、LMD都有较大优势。添加了白噪声的LOD方法与快速ICA方法相互结合产生的新方法
‑
LOD
‑
ICA方法可以将原始的混合信号波形很好的分离还原出来，成功实现混合信号的盲源分离。且与添加白噪声的LOD方法相比较，在混合信号故障特征提取方面，LOD
‑
ICA新方法具有更好的效果。因此提出了一种新的时频分析方法一LOD
‑
ICA算法，结合添加了白噪声的LOD和快速ICA两种算法的优点，克服各自的不足，实现旋转机械设备的故障监测与诊断。
附图说明
[0034]图1为算法的流程图；
[0035]图2为实施例中故障信号时域图；
[0036]图3为实施例中故障信号频谱图；
[0037]图4为实施例中故障信号基于噪声辅助分析的LOD分解时频图；
[0038]图5为实施例中故障信号基于噪声辅助的LOD分解频谱图；
[0039]图6为实施例中混合信号分量快速ICA分解时频图；
[0040]图7为实施例中混合信号分量快速ICA频谱图。
具体实施方式
[0041]本专利技术将添加了白噪声的LOD方法与快速ICA方法结合，提出了可自适应分解混合信号的时频分析新方法
‑
LOD
‑
ICA方法，介绍了该方法以添加了白噪声的LOD方法对故障混合信号进行预处理，为快速ICA方法提供了充足的观测通道数据，进而使用快速ICA方法实现混合信号盲源分离，完成故障信号特征提取。
[0042]算法例提供了一种故障信号特征提取方法，算法流程如图1所示，具体步骤如下：
[0043]步骤1：保存以风力机齿轮箱为代表大型旋转机械设备在实际运行时发生故障用接收机采集的混合信号x(t)；
[0044]步骤2：找出信号x(t)中所有极值X
k
(k＝1，2，...，m)及其对应的时刻值τ
k
，建立任
意两个相邻极值点[X
k
，X
k+1
)取值区间的分段线性函数s1(t)，计算方法如式(1)：
[0045][0046]其中，X
k
(k＝1，2，...，m)为x(t)中所有的极值；t表示时间；X
k+1
为X
k
的相邻极值点；τ
k
、τ
k+1
分别X
k
、X
k+1
对应的时刻值；对原始信号x(t)本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于LOD
‑
ICA的旋转机械故障诊断方法，其特征在于，包括以下步骤：步骤1：保存采集的旋转机械故障混合信号；步骤2：建立任意两相邻极值点区间的分段线性函数，通过坐标域变换将信号从原始数据域坐变换为锯齿域坐标；步骤3：对一阶微分信号添加白噪声，求得锯齿域的初始均值函数，进而求得高频波动函数，再通过坐标域变换运算，将高频波动函数转换为原始域高频函数；步骤4：鉴别分离出来的信号单分量c1(t)是否含有不同的波动特征，即其瞬时频率物理意义是否具有实质性，若是，则c1(t)＝MOC1(t)；若否，最终分离出来首次瞬时频率物理意义具有实质性的单分量信号c
m
(t)，c
m
(t)＝MOC1(t)；步骤5：将原始信号拆分为n个MOC(t)分量和余量r(t)；步骤6：对比分析分离出来的信号频谱图，将仍然存在信号模态混肴问题的MOC分量进行标准化(均值为0，方差为1)处理，得到标准化向量Y，通过白化处理实现数据向量Y
→
Z的转变；步骤7：对随机生成的初始化单位范数向量w进行归一化处理；步骤8：根据推导出的迭代公式进行更新迭代计算；步骤9：根据迭代公式w＝w
*
/||w
*
||对向量w进行标准化处理；步骤10：判断向量矩阵w是否收敛，若否，则重复步骤8
→
步骤10；若是则进行下一步；步骤11：重新确立向量w进行对称正交化处理，保证分离出来的每个单位向量w互相独立正交，且线性无关；步骤12：分离得到信号。2.如权利要求1所述的基于LOD
‑
ICA的旋转机械故障诊断方法，其特征在于，步骤2中，找出通过振动信号传感器获取混合信号x(t)中所有极值及其对应的时刻值，建立任意两个相邻极值点取值区间的分段线性函数：其中，X
k
(k＝1,2,
…
,m)为x(t)中所有的极值；t表示时间；X
k+1
为X
k
的相邻极值点；τ
k
、τ
k+1
分别X
k
、X
k+1
对应的时刻值；对原始信号x(t)在任意两个相邻极值点[X
k
，X
k+1
)的信号区间进行分段线性变换运算，以获得分段线性函数s1(t)；通过原始信号x(t)的坐标域变换运算，对原始数据坐标域坐标(t,x)进行坐标转换坐标变换公式表示如下：s(u)＝x(t)
ꢀꢀꢀꢀ
(3)转换后的锯齿域坐标(u,s)极大的减小了分解误差。3.如权利要求1所述的基于LOD
‑
ICA的旋转机械故障诊断方法，...

【专利技术属性】
技术研发人员：廖力达，叶捍，
申请(专利权)人：长沙理工大学，
类型：发明
国别省市：