基于ceemd能量熵和阶次跟踪的滚动轴承损伤状态识别.docxVIP

基于ceemd能量熵和阶次跟踪的滚动轴承损伤状态识别 0 振动信号的故障监测 旋转机中的滚动轴承是一个必不可少的部件。近1.3个故障是由滚动轴的承载引起的。因此,以滚动轴承为研究对象,开展相应的故障损伤分析具有重要的工程价值。文献[1] 的研究表明,滚动轴承的振动信号与轴承损伤的类型、程度、部位及其原因等密切相关,目前利用振动信号的故障损伤监测方法应用广泛。 近年来,针对滚动轴承的损伤状态评估,学者们开展了大量的研究。Batista等 目前大多数的旋转机械状态识别研究成果局限于稳定旋转速度的工况条件。但通常情况下, 机械旋转系统由于不同的工况条件,会在不同的旋转速度下工作。在变速度的工况条件下,旋转机械故障的特征参数如特征频率会随着旋转速度的变化而改变。根据Loughlin和Tacer 近年来,通过学者的研究,小波熵和EMD熵能够较好地描述设备的退化损伤状态 1 ceemd方法和经验的标准化过程 1.1 ceemd方法,克服模态混叠问题 经验模态分解( empirical mode decomposition, EMD) 是一种适于非线性、非平稳信号的自适应信号处理方法。其本质是对非平稳信号进行平稳化处理,将信号中不同尺度的波动或趋势逐级分解出来。 同时,EMD方法也存在一些不足如模态混叠现象, 为了克服EMD的模态混叠问题,文献[12] 提出了一种集合经验模态分解( ensemble EMD,EEMD) 方法。EEMD利用白噪声具有频率均匀分布的统计特性,当信号加入高斯白噪声后,使信号在不同尺度上具有连续性,改变信号极值点的特性,促进抗混叠分解,有效地避免了模态混叠现象。但同时在信号中加入白噪声,使分解的结果带来了重构误差。 为了克服噪声对分解结果的影响,文献[13]提出了CEEMD方法,能够在保证分解效果与EEMD相当的情况下,克服所添加的噪声对分解结果的影响。 其分解步骤如下: 1) 向待分解信号中加入n组正、负成对的辅助白噪声,从而生成2套集合IMF。 式中: S为原信号; N为辅助噪声; M 2) 对集合中的每个信号进行EMD分解,每个信号得到1组IMF分量,其中第i个信号的第j个IMF分量用c 3) 对多组IMF分量进行组合,得到分解结果如下: 其中c 1.2 eemd分解计算 为了得到CEEMD能量熵,首先将振动信号进行CEEMD分解,计算各个IMF分量的总能量,计算公式如下: 其中N为每个IMF分量的数据点数。 能量熵的计算方法为: 其中p 归一化后得到的特征向量为: 1.3 经验归一化过程 将原始信号经过CEEMD分解后,各IMF分量c 信号c( t) 的解析信号z( t) 可定义为: 式中: 频率便可定义为相位 θ( t) 对时间t的求导: 运用3次样条曲线将IMF分量c 利用y 当y 定义y 式中F( t) 为统一幅值下的纯调频函数,则调幅部分A( t) 为: 根据式( 4) ,可以得到: 此时即完成了经验频率调制和幅度调制过程, 可以发现,经验幅度调制较为平滑而且有效避免了高频波动和幅度调制系数过调。 计算频率,对相位函数微分: 振动信号的时间- 频率- 幅度谱可表示为: 则包含p个IMF分量信号的边际谱为: 2 滚动轴承损伤状态检测和识别 采用美国凯斯西储大学( Case Western Reserve University) 电气工程实验室的滚动轴承试验数据对所提的方法进行验证 采用转速和载荷分别为1 730 r/min和3 HP, 1 750 r / min和2 HP,1 772 r / min和1 HP,1 797 r / min和0 HP四种工况下,故障直径分别为0. 18 mm和0. 54 mm两种状态下的轴承振动信号对所提的方法进行验证。如表1所示,本次试验包括4种不同工况条件和7种轴承损伤状态,针对各个工作状态条件下的每个损伤状态采集20组数据,一共560组数据,( 每种损伤状态80个,其中40个用来训练,40个用作测试) 。 试验中所采用的滚动轴承主要参数为,轴承直径39. 04 mm,滚珠直径7. 94 mm,滚珠数量为9个,则滚动轴承各故障状态下的特征阶次为: 式中: O 将振动信号进行CEEMD分解后,选取包含主要损伤特征的IMF分量,分别进行经验归一化求得频率调制和幅度调制,并计算得到这些IMF分量的时间- 频率- 幅度谱,进一步计算这些IMF分量组成的边际谱。而特征阶次可以较好地消除转速不同带来的影响。根据边际谱和特征阶次,得到各特征阶次对应的边际谱幅值。并且能量熵能够较好地表征各个损伤状态,由式( 6) 计算每组数据的CEEMD能量熵。最后结合特征阶次边际谱值和CEEMD能量熵构建滚动轴承损伤状态识别的特征向量。将训练和测试各280组特征向量进行归一化后,分别对遗传算法优