VMD与FastICA相结合的滚动轴承故障诊断方法技术

本发明专利技术涉及一种VMD与FastICA相结合的滚动轴承故障诊断方法，其步骤：利用现有数据采集设备采集滚动轴承的原始振动信号；将采集到的滚动轴承原始振动信号进行VMD分解；原始振动信号经VMD后分解为k个模态分量，以连续的3个模态分量为一序列组合进行FastICA分析，得到重构故障信号；对重构故障信号进行Hilbert变换，得到重构故障信号的包络谱，在所得包络谱中提取各故障特征下的特征频率、调制特征频率及边频带，进而根据提取的特征频率成分判断滚动轴承是否存在故障以及故障类型。

【技术实现步骤摘要】
VMD与FastICA相结合的滚动轴承故障诊断方法
本专利技术涉及一种滚动轴承故障诊断方法，特别是关于一种VMD与FastICA相结合的滚动轴承故障诊断方法。
技术介绍
滚动轴承作为旋转机械如齿轮箱、涡轮机械的核心部件之一，其故障将直接影响机械设备的稳定运行。受工作环境的影响，滚动轴承原始振动信号一般为非平稳、非线性的，且信号中常常夹杂着噪声，导致振动信号的信噪比降低，给滚动轴承的故障特征提取和故障诊断带来了困难。
技术实现思路
针对上述问题，本专利技术的目的是提供一种VMD与FastICA相结合的滚动轴承故障诊断方法，该方法能够有效的降低滚动轴承原始振动信号所含干扰噪声，凸显故障特征频率，以便更好的进行滚动轴承的特征频率成分的提取和故障诊断。为实现上述目的，本专利技术采取以下技术方案：一种VMD与FastICA相结合的滚动轴承故障诊断方法，其特征在于包括以下步骤：1)利用现有数据采集设备采集滚动轴承的原始振动信号y；2)对振动信号y进行VMD分解；3)原始振动信号经VMD后分解为k个模态分量IMF1,IMF2,‥‥,IMFk，以连续的3个模态分量为一序列组合进行FastICA分析，得到重构故障振动信号；4)对重构故障振动信号进行Hilbert变换，得到重构故障振动信号的包络谱，在所得包络谱中提取各故障特征下的特征频率、调制特征频率及边频带，进而根据提取的特征频率成分判断滚动轴承是否存在故障以及故障类型。进一步，所述步骤2)中，VMD分解过程包括步骤如下：2.1)令m＝0,初始化第k个模态函数uk记为当前模态函数功率谱的中心ωk记为当前m＝0所对应的傅里叶变换λ1，其中m为迭代次数；λ为所对应的傅里叶变换；k∈(1,n)，n为原始振动信号y进行VMD分解后得到的模态个数；2.2)更新第k个模态函数uk、第k个模态的中心频率ωk和λ；2.3)判断更新结果是否满足收敛条件，若满足则停止迭代，否则返回步骤2.2)。进一步，所述步骤2.2)中，根据式更新uk；根据式更新ωk；根据式更新λ；其中uk为第k个模态函数；ui为第i个模态函数，i≠k；为第k个模态的第m-1次迭代的模态函数；相当于第m次迭代时对当前剩余量y-∑i≠kui的维纳滤波；ω为中心频率；ωk为第k个模态的中心频率；为第k个模态的第m-1次迭代的中心频率；为m次迭代模态函数功率谱的中心；λm为第m-1次迭代所对应的傅里叶变换；λm+1为第m次迭代所对应的傅里叶变换；α为惩罚因子。进一步，所述步骤2.3)中，收敛条件为：式中，e预先设定的允许误差。进一步，所述步骤3)中，基于负熵的独立分量分析的具体步骤如下：3.1)将连续的3个模态分量为一序列组合得到的输入矩阵Y，并对输入矩阵Y进行中心化，使其均值为0；3.2)对中心化后的数据进行白化处理，得到白化后的数据z；3.3)选择需要估计的独立成分的个数m，设迭代次数p←1；3.4)随机选取一个具有单位范数的初始权矢量wp；3.5)迭代计算，即更新wp：其中T为矩阵转置，E[·]为均值运算，g[·]为非线性函数,可取g(y)＝tanh(a1y),其中y为随机变量，a1为常数，a1取值范围为1≤a1≤2，故本实施例中取a1＝1；3.6)进行正交：3.7)标准化wp：wp←wp/||wp||；3.8)假如wp尚未收敛，则返回步骤4.5)；3.9)使p←p+1,若p≤m,则返回步骤4.4)，若p＞m,则结束；3.10)令W(:,n)＝wp,则Z＝W′*z，Z为输入矩阵Y经过FastICA处理后的结果，W′为W的转置；3.11)对所有的序列组合均经过3次迭代后终止，得到三个独立分量IC1、IC2和IC3，从中选取一个最优的结果ICr1作为本次FastICA分析的结果；共进行n-2次FastICA分析，从每次结果中选取一个最优分量，共n-2个分量ICr1，ICr2,…，ICr(n-2)，利用所得n-2个分量重构故障振动信号yr。进一步，所述步骤3.5)中，wp更新过程为：其中T为矩阵转置，E[·]为均值运算，g[·]为非线性函数,取g(y)＝tanh(a1y),其中y为随机变量，a1为常数。本专利技术由于采取以上技术方案，其具有以下优点：1、本专利技术利用VMD方法对轴承故障信号进行分解可以有效抑制EMD方法存在的模态混叠，能够将不同频率的信号分量很好的区分开来。2、本专利技术将VMD方法与FastICA相结合可以弥补应用FastICA方法时存在的欠定问题，充分发挥两者在信号处理方面的优势。3、采用本专利技术分析实际滚动轴承内外圈故障信号，可清晰准确提取各故障特征下的特征频率、调制特征频率及边带，且易于实现，适合用于滚动轴承故障诊断。附图说明图1为本专利技术的整体流程示意图。具体实施方式下面结合附图和实施例对本专利技术进行详细的描述。如图1所示，本专利技术提供一种VMD与FastICA相结合的滚动轴承故障诊断方法，其包括以下步骤：1)利用现有数据采集设备采集滚动轴承的原始振动信号y；2)将采集到的滚动轴承原始振动信号y进行VMD分解；其分解过程包括步骤如下：2.1)令m＝0,初始化第k个模态函数uk记为当前模态函数功率谱的中心ωk记为当前m＝0所对应的傅里叶变换λ1，其中m为迭代次数；λ为所对应的傅里叶变换；k∈(1,n)，n为原始振动信号y进行VMD分解后得到的模态个数。2.2)更新第k个模态函数uk、第k个模态的中心频率ωk和傅里叶变换λ；根据式更新uk；根据式更新ωk；根据式更新λ；其中,uk为第k个模态函数；ui为第i个模态函数，i≠k；为第k个模态的第m-1次迭代的模态函数；相当于第m次迭代时对当前剩余量y-∑i≠kui的维纳滤波；ω为中心频率；ωk为第k个模态的中心频率；为第k个模态的第m-1次迭代的中心频率；为m次迭代模态函数功率谱的中心；λm为第m-1次迭代所对应的傅里叶变换；λm+1为第m次迭代所对应的傅里叶变换；α为惩罚因子。2.3)判断更新结果是否满足收敛条件，若满足则停止迭代，否则返回步骤2.2)；收敛条件为：e为预先设定的允许误差，其值取为10-6。3)原始振动信号y经VMD后分解为k个模态分量(IMF1,IMF2,‥‥,IMFk,)，以连续的3个模态分量为一序列组合进行FastICA分析,得到重构故障振动信号；基于负熵的独立分量分析的具体步骤如下：3.1)将连续的3个模态分量为一序列组合得到的输入矩阵Y，并对输入矩阵Y进行中心化，使其均值为0；3.2)对中心化后的数据进行白化处理，得到白化后的数据z；3.3)选择需要估计的独立成分的个数m，设迭代次数p←1；3.4)随机选取一个具有单位范数的初始权矢量wp；3.5)迭代计算，即更新wp：其中T为矩阵转置，E[·]为均值运算，g[·]为非线性函数,可取g(y)＝tanh(a1y),其中y为随机变量，a1为常数，a1取值范围为1≤a1≤2，故本实施例中取a1＝1；3.6)进行正交：3.7)标准化wp：wp←wp/||wp||；3.8)假如wp尚未收敛，则返回步骤4.5)；3.9)使p←p+1,若p≤m,则返回步骤4.4)，若p＞m,则结束；3.10)令W(:,n)＝wp,则Z＝W′*z，Z为输入矩阵Y经过FastICA处理后的结果，W′为W的转置。3.11)对所有的序列本文档来自技高网...

【技术保护点】
1.一种VMD与FastICA相结合的滚动轴承故障诊断方法，其特征在于包括以下步骤：1)利用现有数据采集设备采集滚动轴承的原始振动信号y；2)对振动信号y进行VMD分解；3)原始振动信号经VMD后分解为k个模态分量IMF1,IMF2,‥‥,IMFk，以连续的3个模态分量为一序列组合进行FastICA分析，得到重构故障振动信号；4)对重构故障振动信号进行Hilbert变换，得到重构故障振动信号的包络谱，在所得包络谱中提取各故障特征下的特征频率、调制特征频率及边频带，进而根据提取的特征频率成分判断滚动轴承是否存在故障以及故障类型。

【技术特征摘要】
1.一种VMD与FastICA相结合的滚动轴承故障诊断方法，其特征在于包括以下步骤：1)利用现有数据采集设备采集滚动轴承的原始振动信号y；2)对振动信号y进行VMD分解；3)原始振动信号经VMD后分解为k个模态分量IMF1,IMF2,‥‥,IMFk，以连续的3个模态分量为一序列组合进行FastICA分析，得到重构故障振动信号；4)对重构故障振动信号进行Hilbert变换，得到重构故障振动信号的包络谱，在所得包络谱中提取各故障特征下的特征频率、调制特征频率及边频带，进而根据提取的特征频率成分判断滚动轴承是否存在故障以及故障类型。2.如权利要求1所述方法，其特征在于：所述步骤2)中，VMD分解过程包括步骤如下：2.1)令m＝0,初始化第k个模态函数uk记为当前模态函数功率谱的中心ωk记为当前m＝0所对应的傅里叶变换λ1，其中m为迭代次数；λ为所对应的傅里叶变换；k∈(1,n)，n为原始振动信号y进行VMD分解后得到的模态个数；2.2)更新第k个模态函数uk、第k个模态的中心频率ωk和λ；2.3)判断更新结果是否满足收敛条件，若满足则停止迭代，否则返回步骤2.2)。3.如权利要求2所述方法，其特征在于：所述步骤2.2)中，根据式更新uk；根据式更新ωk；根据式更新λ；其中uk为第k个模态函数；ui为第i个模态函数，i≠k；为第k个模态的第m-1次迭代的模态函数；相当于第m次迭代时对当前剩余量y-∑i≠kui的维纳滤波；ω为中心频率；ωk为第k个模态的中心频率；为第k个模态的第m-1次迭代的中心频率；为m次迭代模态函数功率谱的中心；λm为第m-1次迭代所对应的傅里叶变换；λm+1为第m次迭代所对应的傅里叶变换；α为惩罚因子。4.如权利...

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