基于EMD方法对往复式压缩机激励源定位

本资料包含pdf文件1个，下载需要1积分

往复式压缩机是生产企业中广泛使用的动力设备，应用于冶金、矿山、化工、船舶、机械制造等行业。尤其是在化工生产上的应用更加普遍.往复式压缩机的故障诊断中使用较多的是振动法。由于其机械结构复杂、运动部件多、工作时振动激励源较多，发生的故障也多种多样，发生故障时，振动信号表现出非平稳、非线性特征，利用传统的时域或频域方法很难从此类振动信号中提取故障信息，诊断难度较大[1-3]。

往复式压缩机机体的主要受到机械力和气体力的作用，其中机械力主要表现为曲轴旋转时产生的旋转惯性力和活塞组件往复运动时产生的往复惯性力，这两种惯性力最终将通过曲轴与机体连接收稿日期：2013-01-14作者简介： (1988- )，男，硕士研究生，主要从事现代信号处理理论、方法及其在故障诊断中的应用等方面的研究。E-mail：zznsurmy###grtmil.com。

通信作者：迟毅林(1953- )，男，教授，博士生导师，主要从事故障监测与诊断、机械设备的集成与创新等方面的研究.E-mail：yilinchi###sina.com。

76 江西理工大学学报 2013年6月的轴承座传递到机体上.气体力主要表现为气缸缸内的压力，这部分的压力将通过汽缸内壁最终传递到机体上.所以在往复式压缩机机体的振动。最终是由机械力和气体力这两种激励联合起来产生的振动[41.其中机械力这-部分产生的故障.称之为动力性能故障；气体力引起的故障，称之为热力性能故障.对这两种激励联合产生的激励信号，利用传统的信号和处理算法，难以准确地对其进行特征分析，因而如果产生了故障，难以准确地判断故障属于哪-类型[5-63.由于压缩机的机械力的振动是通过-个曲柄连杆机构传递的到压缩机机体上的。

先对压缩机机体的振动信号进行 EMD分解，得到信号在不同频率段内的 IMF分量.然后在对分解后的各 IMF分量做频域分析，并对压缩机十字头上采集到的振动信号进行频域分析，从而分离出压缩机机械力振动引起的振动分量。

1 EMD分解经验模态分解法(EMD分解)是黄锷教授等人提出的，用于处理非平稳、非线性信号数据的-种方法.相较于以前的-些方法，该方法具有直观、直接、自适应性等特点，分解方法是直接基于数据本身得出的[71.分解是根据简单的假设，即任何数据由不同的简单振荡固有模式叠加而成.每-个固有模式都代表了-种有着相同数 目极点和零点的简单的振荡.在任何特定时间内，数据可能有许多不同的振荡模式并存，因为-个振荡在另-个振荡上叠加，最后的结果还是复杂的数据.这些振荡模式。每个代表-个固有模态函数(IMF)有如下定义：①在整个数据集中，极值点的数目与零点值的数 目必须相等或最多相差-个；②在任意瞬时，由局部最大值和局部最小值确定的包络的均值为零[8-O1。

各个 IMF分量突出信号的各个局部特征，对其进行分析可以更准确把握原信号的特征信息.EMD分解方法的本质是通过特征时间尺度获得固有振动模式，然后由固有振动模式来分解时间序列数据.运用 EMD方法分解信号的步骤如下：(1)确定信号x(t)的所有局部极大值点和局部极小值点，采用三次样条线将所有的局部极大值和局部极小值连接起来形成上下包络线.将上下包络线的均值记为 m ，原始信号减去 m 可得到-个新的数据序列 h .即：hi (t)-ml (1)(2)如果 h。满足 IMF分量的两个定义 ，则 h为第-个 IMF分量 ，记为 c ；如果不满足，则重复步骤(1)，得到上下包络线的均值 m 并继续判断 -m 。是否满足IMF分量的两个定义.若不满足，则继续重复循环步骤(1)，直至第 k次循环的h。 满足 IMF的两个条件时，将 h 记为c 的第-个 IMF分量。

(3)将第-个 IMF分量 C。从原始信号想 (t)中分离出来，得到：rl- (t)-c1 (2)并将 r 作为新的信号，用步骤(1)、(2)进行分解，得到c：.这样重复 n次循环，直到最后得到 rn为单调函数时停止.所以有：∑ cirn (3)其中：C 表示第 个 IMF分量， 为分解的剩余分量。

分量 C。，C：，，c 分别包含了信号从高到低不同频率段的成分，而且不是等带宽的㈣。

2 压缩机振动信号的EMD分析对某压缩机厂所使用 DW-47/1-2.2型氢气循环压缩机的机体进行振动信号提取.将加速度传感器安装在压缩机气缸的端盖上，采样频率为25.6 kHz，采样点数为 25 600.图 1是往复式压缩机气缸的振动信号.由图 1可知，在往复式压缩机运转的时候会产生周期性的较为规律的冲击.这是由于往复式压缩机的机构特性产生的.是不可避免的冲击.对该原始信号进行频域分析(如图2)。

2墨0- 2山- lI I I-lIl 1.山4 IlI-J I IIJI叩 Pl l1 -l1lf1 唧0 0.2 0.4 0.6 0.8t/s图 1 气缸振动信号225 -岛 ～ 。

频 率 /Hz图2 气缸振动信号的幅值谱第34卷第 3期 张赵宁.等：基于ENID方法对往复式压缩机激励源定位 77从信号的幅值谱中可以看到，该段信号有三个比较明显的峰值频率 225 Hz、.厂2541 Hz、jS783 Hz.但由于往复式压缩机机体上的振动信号的振源只有两个：-种是曲轴旋转时产生的旋转惯性力与活塞组件往复运动时产生的往复惯性力的合力，另-种是气缸缸内的压力产生的气体力[1 2。

那么上述的三个峰值频率中，只有两个峰值频率是机械力和气体力振动的频率，剩下的-个峰值频率，推测其为其他噪声的干扰频率。

对压缩机气缸的振动的原始信号进行 EMD分解，能得到-共 14阶的IMF分量.图 3从上到下依次为截取的气缸振动信号的前八阶固有模态函数 c。～c 由图 3可知，1-6阶的 IMF分量都趋于-个简单振荡信号.7、8阶的 IMF分量趋于正弦衰减.对于趋于正弦衰减固有模态函数 。由于其无明显的振荡峰值 ，故不做研究.因此.对第 8阶以后IMF分量做进-步的研究意义不大.故选择对前 8阶的 IMF分量做进-步的分析，找出其中有用的IMF分量。

5o- 520- 220 o- 2EMD分解前 8阶IMF分量0 0.1 0.2 O.3 0.4 0.5 0.6 0.7 0.8 0.9 1IS 0- 10.50- 0.50 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1f/s图 3 气缸振动信号的EMD分解结果对经过EMD分解后的前 8阶IMF分量进行幅值谱的分析，从图4可以看出，c 这个 IMF分量有- 个 783 Hz的峰值频率.在 c 、c 这两个 IMF分量上，都能找到225 Hz这-峰值频率.而A541 Hz这-峰值频率，并未出现在 IMF分量里面，即说明前面的推测正确 ，≥54l Hz为噪声的干扰频率。

前 8阶IMF分量的幅值谱0.O4广----T-----广----T-----T----]-----r---] L蛐 I-J 10000 0O 2 000 4 000 6 0o0 8 000 10 ooO 12 ooO 14 0o00 2 000 4 000 6 000 8 000 10 oo0 12 00014 00000O0000O 2 000 4 000 6 000 8 000 10 0oo l2 0o0 14 0oo0 2 000 4 000 6 000 8 000 10 0oO 12 00o 14 000频率，Hz图 4 EMD 分解 结果的幅值谱对压缩机十字头的振动信号进行频域分析，如图5、图6所示.我们可以看到有-个明显峰值频率f783 Hz.这-峰值频率与图4中c 这个 IMF分量的峰值频率相对应.由此我可以推测。具有783 Hz这-峰值频率的 IMF分量是机械力所产生的振动分量，其反应的故障就是动力学故障；那么另外两组，具有 225 Hz这个峰值频率的IMF分量就可能是气体力所产生的振动分量.其反应的故障就是热力学故障.由于实验条件有限。对于气体力的振动测试的实验，现阶段还没法完成 ，需要后续的实验进-步论证。

78 江西理工大学学报 2013年6月201O0馨- 10- 20理鏖 l - j I. J . J- Ir 1 1 1 T TrT rf 1r - lOt/s图 5 十字头振动信号I J -。

I Il JIliI 且 -频率 ，H图 6 十字头振动信号的幅值谱本实验进-步论证 EMD分解方法具有将-个复杂的非线性非平稳的信号根据其自适应的特点分解为不同品频段上的 IMF分量的特点.根据往复式压缩机两种激励源的频率分布的差异，将气缸振动的混合信号通过 EMD方法分解成为不同频段上的IMF分量，进而通过监测不同 IMF分量，来确定往复式压缩机的不同故障来源。

3 结束语与旋转机械相比。往复压缩机结构复杂。激励源多，振动信号包含着太多的非平稳性和非线性的成分，因而针对往复机械的精确诊断的实现要困难得多.文中从往复式压缩机的结构出发 .根据其运动特性.分析出气缸振动信号的激励源来自两部分，即机械力和气体压力两部分.由于 EMD分解方法是根据被分析信号本身的特点，自适应地将信号分解为不同频段上的 IMF分量之和，从而降低了信号的非平稳性、非线性的特性[D-I41。

所以根据气缸振动信号的两个振动来源频率的高低不同，先通过 EMD分解的方式，将气缸的振动信号分解成若干个 IMF分量的形式.再对每个IMF分量进行频域分析。找出其中比较敏感的几个IMF分量，与测得的机械力的振动信号的频谱图做对比，找出属于机械力激励的那-部分的IMF分量。

然后通过最初气缸振动信号的频谱图的分析，进而推断出属于气体力激励的那-部分 IMF分量.从而为以后压缩机故障诊断提供-个理论上的依据：对测得气缸振动信号做 EMD分解.通过监测不同的IMF分量，获取往复式压缩机不同的故障信息。