基于emd旋转机械裂纹松动耦合故障诊断方法探究

《基于emd旋转机械裂纹松动耦合故障诊断方法探究》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、基于EMD旋转机械裂纹松动耦合故障诊断方法探究摘要：针对旋转机械耦合故障的诊断问题，建立了含有裂纹-松动耦合故障的转子-轴承系统动力学模型，并用龙格库塔法求出故障模型振动信号。利用EMD(EmpiricalModeDecomposition)方法对振动信号进行分解，得到含有故障特征的本征模式函数(IntrinsicModeFunction,简称IMF)。对IMF做希尔伯特变换得到振动信号边界谱，通过分析边界谱的倍频情况并与单一故障信号作比较，诊断出故障信号同时含有裂纹和松动故障特征，说明该故障系统存在裂纹松动耦合故障，并证明EMD方法在旋转机械耦合故障诊断方面的

2、有效性。关键词：故障诊断经验模态分解裂纹松动中图分类号：TH133文献标识码：A文章编号：1672-3791(2012)10(b)-0049-03转轴裂纹和支撑部件松动是旋转机械的常见故障，也是导致机械系统失效甚至造成严重事故的主要原因。目前国内外学者对裂纹、松动故障单独存在的系统研究较多，提出了许多诊断方法。文献［1］研究斜裂纹的动力特性，指出随着裂纹深度的增加，横向响应的组合频率增多。文献［2］结合物理模型与灰色理论，提出行星轮系齿根疲劳裂纹故障预测的新思路，对试验中的疲劳裂纹进行定量检测和故障预测。文献［3］针对工程中出现的支座松动故障，建立了多盘悬臂转子

3、的松动有限元模型，对单支座和双支座松动故障进行动力学特性研究。但在实际转子中，常常出现两种故障同时存在的情况。这种耦合故障转子的动力学行为较单一故障转子更加复杂，而且相互影响，不容易诊断。文献［4］利用求解非线性非自治系统周期解的延拓打靶方法，研究了松动裂纹耦合故障转子轴承系统周期运动的稳定性及其失稳规律。文献［5］建立了带有裂纹-支承松动耦合故障的双跨弹性转子系统动力学模型，利用数值仿真对故障非线性响应进行研究。EMDE6］是近年来发展起来的处理非平稳、非线性信号的时频分析方法。该方法克服了传统时频分析方法中的不足，具有很强的自适应性，并在机械故障诊断领域得到

4、了广泛应用［7〜10］。针对耦合故障信号复杂，具有强非线性的特点,本文提出一种基于EMD的耦合故障诊断方法。该方法先利用EMD将故障信号分解，然后求得有效IMF的边界谱，通过对边界谱分析判断系统状态，达到故障诊断的目的。1系统力学模型和运动微分方程如图1所示，建立含有裂纹-松动耦合故障的刚性支承转子-轴承系统模型，转子圆盘与轴承之间为无质量的弹性轴。模型左端发生松动，轴承座与基础之间的松动最大间隙为。转子圆盘左侧有一弓形横向裂纹，其深度为3。图1中01为轴瓦几何中心；02为转子几何中心；03为转子质心，k为弹性轴刚度；ml为两端轴承处的转子集中质量；m2为转子圆

5、盘的等效集中质量；m3为轴承支座的等效集中质量。模型还考虑了左端滑动轴承作用在转轴上的非线性油膜力，为别为Fx、Fy。设转子右端的径向位移为xl,yl；转子圆盘的径向位移为x2,y2；松动端轴心位移为x3,y3；轴承支座在竖直方向位移为y4,则具有裂纹松动耦合故障的转子-轴承系统运动微分方程为：式中U为转子的偏心量;cl为转子在轴承处的阻尼系数;c2为转子圆盘的阻尼系数；cs为支座松动阻尼系数；ks为支承刚度。为转子转速；g为重力加速的；、为仅与裂纹深度a有关的相对刚度参数。为裂纹开闭函数，本文采用余弦波模型来表示裂纹开闭过程，粗略地考虑裂纹半开半闭的过渡过程，

6、忽略了裂纹的全闭和全开是一个持续过程。余弦波模型的数学表达式为：图2所示，式中为初相位；为裂纹方向与偏心之间的夹角；x,y为转子初始位置松动故障等效成刚度和阻尼的变化;支承间隙系统在位移条件下ks、cs为分段性，其表达式为:式（1）中油膜力沿x和y两个方向的分量为：式（4）中为润滑油粘度；为转子转速；c为轴承径向间隙；R为轴承半径；L为轴承长度。2经验模式分解经验模式分解EMD是一种自适应分解方法，可以把复杂的信号分解为有限个IMF分量。IMF信号一般满足两个条件:(1)从全局特性上看，极值点数必须和过零点数一致或者至多相差一个。(2)在某个局部点，极大值包络和

7、极小值包络在该点的算术平均值是零，即两条包络线关于时间轴对称。我们可以把任何信号按下面步骤分解。(1)用三次样条线将所有的局部极大值点连接起来形成上包络线。(2)用三次样条线将所有的局部极小值点连接起来形成下包络线。(3)上下包络线的平均值记为，求出：理想地，如果是一个IMF,那么就是的第一分量。(4)如果不满足IMF的条件，把作为原始据，重复(1)、(2)、(3),得到上下包络线的平均值再判断是否满足IMF的条件，如不满足，重复循环k次，得到，使得满足IMF条件。记，则为信号的第一个满足IMF条件的分量。(5)将从中分离出来，得到：将作为原始数据重复以上过程，

8、得到的第二个满足IMF条