基于ITD和改进形态滤波的轴承故障诊断方法技术

本发明专利技术公开了一种基于ITD和改进形态滤波的轴承故障诊断方法，为基于固有时间尺度分解(Intrinsic time scale decomposition,ITD)和改进形态滤波的轴承故障诊断方法，属于机械故障诊断领域，其内容包括如下步骤：首先对轴承振动信号进行ITD分解得到一系列旋转分量(Proper rotation,PR)，再以峭度为准则筛选出含故障信息丰富的有效PR，对每个有效分量进行改进形态滤波—ACDIF提取冲击成分然后信号重构，最后利用频谱分析提取重构信号中的故障特征。本发明专利技术利用了ITD方法对信号的自适应分解能力和ACDIF的有效滤噪并保留冲击成分的优点，提高了滤波输出的准确度。ACDIF滤波器不仅结构设计简单而且可应用于振动信号的故障诊断中，通过该方法滤波之后的信号冲击特征被完整保留。

【技术实现步骤摘要】
基于ITD和改进形态滤波的轴承故障诊断方法
本专利技术属于机械设备的故障诊断领域，涉及一种基于ITD和改进形态滤波的轴承故障诊断方法。
技术介绍
滚动轴承的状态监测与故障诊断一般以采集的振动信号作为分析对象，通过信号分析得到故障特征实现诊断。由于轴承工况的不稳定和零部件的损伤等引起的非线性振动，采集到的信号大多表现出非线性、非平稳的特征；同时不可避免地受到各种噪声与信号调制干扰的影响。因此，如何有效地抑制噪声等干扰，从成分复杂的振动信号中提取出故障特征成为故障检测的关键。时频分析法是处理非线性、非平稳信号的常用方法，如短时傅里叶变换(ShorttimeFouriertransform,STFT)、Wigner-vill分布(WVD)、小波变换(Wavelettransform,WT)、经验模态分解(Empiricalmodedecomposition,EMD)等能同时从时域和频域揭示信号成分，实现故障检测，在机械故障诊断中取得了广泛应用。短时傅里叶变换方法简单高效，但受Heisenberg不确定原理的限制，时频分辨率较低；WVD能够描述信号的能量分布情况，在分析多分量信号时，存在比较严重的交叉项干扰；小波变换具有可变的时频窗口，但要求选择基函数，有限的时宽也有可能导致能量泄漏；EMD基于信号本身的局部特征时间尺度，虽有很大发展，但还存在端点效应、模态混叠和负频率等问题。固有时间尺度分解(Intrinsictimescaledecomposition,ITD)是由Frei提出的另一种自适应时频分析方法，该方法能将复杂信号分解为若干相互独立且瞬时频率具有物理意义的固有旋转分量(Properrotation,PR)之和。与EMD相比，ITD采用线性插值取代了基于三次样条曲线的基线信号构建方法，得到的PR分量保留了原始信号的极值时刻，不会出现过包络和欠包络现象，同时端点效应被严格限制在两端。原始信号经ITD分解之后，进一步突显了故障信号的局部特征信息，但噪声抑制作用不明显，高频PR分量仍含有大量噪声，直接对高频分量进行故障分析，噪声干扰会严重影响诊断结果的准确性。不同于时频分析法，形态滤波(Morphologicalfilter,MF)是基于数学形态学变换的一种重要的非线性滤波方法，它利用结构元素与信号的形态特征进行匹配，可以有效地提取信号的边缘轮廓和形状特征，其原理简明、运算简单高效，且具有明确的物理意义。
技术实现思路
本专利技术的目的在于提供一种基于ITD和改进形态滤波的轴承故障诊断方法，针对现有技术中单一的ITD分解降噪方法和故障特征提取技术的不足，以及故障信号受谐波和随机噪声干扰严重的特点和形态学基本理论，改进振动信号降噪和特征提取技术，有效滤除谐波干扰,提取强背景噪声下的冲击故障特征，实现设备的故障诊断。为了实现上述目的，本专利技术提供如下技术方案：首先对采集到的振动信号进行ITD分解得到一系列旋转分量(Properrotation,PR)，每个PR分量包含着原始信号部分特征信息，其复杂程度相对简单。初步细化信号后以峭度为准则筛选出含故障信息丰富的有效PR，去除多余分量避免噪声分量的干扰。再对每个有效分量进行ACDIF形态滤波，改进的形态滤波方法可以有效地提取信号的边缘轮廓和信号的形状特征，并同时提取出与故障相关的正负冲击成分；形态滤波过程中，采用TEK作为评价指标自适应筛选最佳结构元素长度,有效提高了滤波处理的效率和精确度。最后将滤波处理过的信号重构避免细节信息的丢失，并利用频谱分析提取重构信号中的故障特征实现故障诊断。具体步骤如下：一种基于ITD和改进形态滤波的轴承故障诊断方法，其特征在于，包括如下步骤：第一，对采样信号x进行ITD分解，得到多个PR分量；第二，利用峭度指标从步骤一中筛选出包含故障信息最多的分量为有效PR分量；第三，以TEK为依据计算出步骤二每个有效PR分量信号的最佳结构元素长度Li；第四，结合步骤三中的最佳结构元素长度Li，分别对步骤二有效分量PR进行ACDIF抑制随机噪声，突出其冲击成分；第五，重构步骤四中形态滤波后的各个分量得到合成信号x′；第六，对步骤五合成信号x′进行快速傅里叶变换(FastFouriertransform,FFT)，提取故障特征频率，实现故障诊断。所述第一步骤中采用ITD对采集到的振动信号进行分解，初步滤除干扰信号，突出信号局部特征。所述第二步骤中利用峭度指标筛选出包含故障信息最多的分量为有效PR分量，峭度的计算公式为：其中:σ和μ分别表示信号的标准差和均值，E(t)表示变量t的数学期望。所述第三步骤中以TEK为依据计算出每个有效分量信号的最佳结构元素长度Li，Teager能量峭度(Teagerenergykurtosis,TEK)正是峭度和Teager能量算子的融合，对于一维离散信号X＝{xt；t＝1,2,…,N}，其TEK定义如下：其中，为离散信号X在采样点t处的Teager能量，是的平均值。TEK结合了峭度和Teager能量算子的优缺点，对信号的冲击成分很敏感，同时考虑了对尖峰的描述及信号瞬时能量的变化可以更加精确有效地反映出故障信号的冲击特征。引入TEK作为筛选结构元素最佳长度的评判指标，TEK会跟随原始信号及形态运算方法的变换而改变，具有一定的自适应性。所述第四步骤中结合步骤三计算出的最佳结构元素，对步骤二筛选出的有效PR分量分别进行ACDIF抑制随机噪声，突出其冲击成分，ACDIF滤波的具体方法为：首先定义膨胀算子和闭算子的结合：膨胀-闭(Dilation-Closing,DC)或闭-膨胀(Closing-Dilation,CD)滤波算子，定义式如下：再定义腐蚀算子和开算子的结合：腐蚀-开(Erosion-Opening,EO)或开-腐蚀(Opening-Erosion,OE)滤波算子，定义式如下：为了同时提取正负冲击，取FDC或FCD与FEO或FOE的差值为新的组合差值形态算子。CDIF的滤波输出形式主要有以下四种：当f(n)满足一定条件时，即FCD-OE(f(n))＝MG(f(n)。)对于同一信号f(n)，每种运算滤波效果各异。通过实验分析比较，发现FCD-EO和FDC-EO的输出信号的冲击幅值较高，相对于其余两种滤波，整体最优。所以，采用两者的平均滤波输出—平均组合差值形态滤波作为最终结果，即整合了FCD-EO和FDC-EO的优缺点，既可以有效抑制随机噪声又能够突出信号冲击特征，兼顾细节信息。所述第五步骤中重构形态滤波后的各个PR分量得到合成信号，避免单一分量分析所导致的细节信息损失，使得故障诊断更加准确。本专利技术的有益效果是：组合四大基本数学形态算子构造出一种新的形态滤波器—平均组合差值形态滤波(ACDIF)，不仅结构设计简单而且可应用于振动信号的故障诊断中，通过该方法滤波之后的信号冲击特征被完整保留，干扰基本滤除，信号幅值谱中故障频率及其倍频谱线清晰显著。利用滤波后信号的TEK值筛选最佳结构元素长度，充分考虑了故障信号的冲击特性和瞬时能量的变化，提高了形态滤波处理的效率，也使得滤波输出更加精确，具有一定的自适应性。将ITD和ACDIF相结合，利用了ITD方法对信号的自适应分解能力和ACDIF的有效滤噪并保留冲击成分的优点，提高了滤波本文档来自技高网...

【技术保护点】
一种基于ITD和改进形态滤波的轴承故障诊断方法，其特征在于，包括如下步骤：第一步，对采样信号x进行ITD分解，得到多个PR分量；第二步，利用峭度指标筛选出包含故障信息最多的分量为有效PR分量；第三步，以TEK为依据计算出每个分量信号的最佳结构元素长度Li；第四步，结合最佳结构元素，对有效PR分量分别进行ACDIF抑制随机噪声，突出其冲击成分；第五步，重构形态滤波后的各个分量得到合成信号x′；第六步，对合成信号x′进行快速傅里叶变换，提取故障特征频率，实现故障诊断。

【技术特征摘要】
1.一种基于ITD和改进形态滤波的轴承故障诊断方法，其特征在于，包括如下步骤：第一步，对采样信号x进行ITD分解，得到多个PR分量；第二步，利用峭度指标筛选出包含故障信息最多的分量为有效PR分量；第三步，以TEK为依据计算出每个分量信号的最佳结构元素长度Li；第四步，结合最佳结构元素，对有效PR分量分别进行ACDIF抑制随机噪声，突出其冲击成分；第五步，重构形态滤波后的各个分量得到合成信号x′；第六步，对合成信号x′进行快速傅里叶变换，提取故障特征频率，实现故障诊断。2.根据权利要求1所述的基于ITD和改进形态滤波的轴承故障诊断方法，其特征在于，所述第一步中，采用ITD对采集到的振动信号进行分解，初步滤除干扰信号，突出信号局部特征。3.根据权利要求1所述的基于ITD和改进形态滤波的轴承故障诊断方法，其特征在于，所述第三步中，采用TEK指标筛选形态滤波的最佳结构元素长度，一维离散信号X＝{xt；t＝1,2,…,N}的TEK定义为：其中，为离散信号X在采样点t处的Teager能量，是的平均值。4.根据权利要求1所述的基于ITD和改进形态滤波的轴承故障诊断方法，其特征在于，所述第四步中，对有效分量分别进行ACDIF抑制随机噪声，突出其冲击成分，ACDIF滤波的具体方法为：首先定义膨胀算子和闭算子的结合：膨胀-闭(Dilation-Closing,DC)或...

【专利技术属性】
技术研发人员：余建波，吕靖香，程辉，孙晓东，晁晓娜，
申请(专利权)人：同济大学，上海航天设备制造总厂，
类型：发明
国别省市：上海,31