专利名称：一种基于三维空间轴心轨迹的压缩机轴系故障诊断方法
技术领域：
本发明属于机械设备故障诊断等领域，涉及机械振动信号的分析及处理，具体涉 及一种基于三维空间轴心轨迹的压缩机轴系故障诊断方法。
背景技术：
转子是压缩机的主要部件，其运行速度高、故障危害大，一直是压缩机监测诊断技 术研究的重点对象。转子在轴承中高速旋转时不仅围绕自身中心旋转，还环绕某一中心作 涡动运动。在涡动的过程中，由于受到周期性的离心力的激励作用，转子将产生横向强迫振 动，从而影响转子的平稳运转，甚至会引发严重的设备事故。转子的回转中心涡动运动的轨 迹则称为轴心轨迹，因此常规的轴心轨迹是一个二维的图像。监测轴心轨迹并提取其特征 是压缩机故障诊断的重要方法。轴心轨迹的获取一般采用两个互成90°安置的非接触式电 涡流传感器，在各自的方向上测量转轴组件相对机座的振动，在复平面对其进行组合后即 可得到轴心轨迹。然而在实际中，转子运转中除了有X和Y两个方向的横向振动外，常常还包括轴向 的振动。轴向的振动往往与横向振动相关，相互之间存在一定的耦合。事实上，现场机组除 了在X和Y两个方向上安装了相应的传感器外，还沿轴向(Z方向)安装了传感器来检测转 子运行过程中轴向位置的变化。本发明中公开的基于三维空间轴心轨迹的压缩机轴系故障 诊断方法将三个传感器信号(X方向、Y方向和Z方向)合成得到三维空间轴心轨迹，三维 空间轴心轨迹考虑了转子轴向振动这个重要诊断信息，将三个方向传感器获取的振动信号 集成起来，能够更为全面的反映转子振动状态，可为检测和诊断提供更加丰富的信息。同时 充分利用多源信息的功能，可以大大提高故障诊断的准确度和效率。为了能让获得的振动信号更加清晰的表现压缩机轴系故障，我们将原始轴心轨迹 进行了处理。即利用信号处理方法、对构成轴心轨迹的振动信号分别进行滤波和提纯的基 础上，提出了三维滤波轴心轨迹和三维合成轴心轨迹。利用三维滤波轴心轨迹和三维合成 轴心轨迹能够消弱信号中噪声的干扰，更好的突现压缩机轴系的实际振动特征，为故障诊 断提供了更为全面的信息。
发明内容
在全息谱技术的基础上，提出了三维振动分析方法，该方法包括合成三维振动分 析与滤波三维振动分析，能够更加清楚、生动地反映设备的运行状态和故障特征。转子的振 动包含水平、垂直以及轴向的三个方向振动。我们将三个方向的振动信号加以合成，构造出 三维轴心轨迹，包含三维滤波轴心轨迹以及三维合成轴心轨迹。这种方法能够真实的反映 转子轴心的运动情况，为故障诊断提供了更为全面的信息。本发明的基本操作步骤如下(1)在压缩机转子轴颈处某一截面或同时在某几个截面上安装两支相互垂直的非 接触式电涡流传感器，然后在压缩机转子轴轴向安装相同的电涡流传感器，得到一组振动信号 x(t)、y(t)、z(t)；(2)按照轴心轨迹的合成方法，将三个传感器所获得信号的瞬时幅值一一对应地 作为X-Y-Z空间上的坐标点(X，y，ζ)。随着转子的转动，该点即形成X-Y-Z空间中的三维 轴心轨迹；(3)将三个方向的振动信号x(t)、y(t)和z(t)分别进行滤波，然后再将滤波后的 信号χ' (n)、y' (η)和ζ' (η)在时域内重新构造，得到三维滤波心轨迹；(4)利用频谱分析原理，将振动信号先分解为一系列的简谐振动分量，再将故障 诊断时所需的频率分量复合起来，然后按照轴心轨迹的合成方法进行合成得到合成轴心轨 迹；合成轴心轨迹中的频率分量不是固定的，是随着条件变化的，需要由诊断人员在诊断时 根据条件自己确定。具体操作步骤如下(1)根据需要在压缩机转子轴颈处某一截面或同时在某几个截面上安装非接触式 电涡流传感器，得到一组振动信号X (t)、y (t)；(2)在压缩机转子轴轴向安装相同的电涡流传感器，得到振动信号z(t)；(3)将步骤(1)和(2)中得到的振动信号x(t)、y(t)、z(t)进行低通滤波剔除高 频干扰信号；(4)设w(t)为在X-Y-Z三维空间进行组合后的三维轴心轨迹，则其表达式为 w(t) = ix(t)+jy(t)+kz (t)；(5)按照轴心轨迹的合成方法，将三个传感器所获得信号的瞬时幅值一一对应地 作为X-Y-Z空间上的坐标点(X，y，ζ)。随着转子的转动，该点即形成X-Y-Z空间中的三维 轴心轨迹；(6)将三个方向的振动信号x(t)、y(t)和z(t)分别进行滤波，然后再将滤波后的 信号χ' (n)、y' (η)和ζ' (η)在时域内按照三维滤波轴心轨迹的构造方法重新构造，得 到三维滤波心轨迹；(7)利用频谱分析原理，将振动信号x(t)、y(t)和z(t)分别先分解为一系列的简 谐振动分量，再将故障诊断时所需的频率分量复合起来，然后按照三维合成轴心轨迹的构 造方法进行构造得到三维合成轴心轨迹。所述的三维滤波轴心轨迹的构造方法为(1)设信号χ (t)、y⑴和ζ⑴经低通滤波后的采样量化值为χ (n)、y (η)和ζ (η)， 滤波器的冲激响应为h (η)；(2)对步骤⑴中采样量化后的信号χ (n)、y (η)和ζ (η)进行滤波，滤波后的信号 为χ' (n)、y' (η)和ζ' (η)，其表达式为: (3)由此可得滤波后的轴心轨迹表达式为W' (n) = ix' (n)+jy' (η) +kz' (η)(2)(4)滤波后轴心轨迹的向量表达式为
(5)按照轴心轨迹的合成方法，将滤波后的三组信号χ' (n)、y' (η)和ζ' (η)
随着转子的转动，该点即形
的瞬时幅值一一对应地作为X-Y-Z空间上的坐标点(X，y，ζ) 成X-Y-Z空间中的三维轴心轨迹；所述的三维合成轴心轨迹的构造方法为(1)设轴心轨迹Χ、Υ和Z方向的振动信号分别为x(t)、y(t)和ζ (t)；(2)将上述三个方向上的信号分解为各频率分量下的一个个分量x(t) = A1Sin (ω 0t+Φ》+A2Sin (2 ω 0t+Φ 2)+...+Ansin (η ω 0t+Φ n)= X1 (t) +X2 (t) +... +χη(t)y (t) = B1Sin (ω 0t+ Ψ i) +B2Sin (2 ω 0t+ Ψ2) +··· +Bnsin (η ω 0t+ Ψη) =Yiζ (t) = C1Sin (ω 0t+ ξ ^ +C2sin (2 ω 0t+ ξ 2) +... +Cnsin (η ω 0t+ ξ η) =z^O+z^O+'-'+z^t)其中，An，Bn，Cn，Φη，Ψη, ξη——水平、垂直和轴向信号的幅值与相位；ω0——信号的旋转角频率； η——频率分量的最高阶数；(3)分别把分量Xi (t)、y, (t)和Zi (t)合成为方程组的形式 (4)根据信号的各个振动分量信息，可以将各分量合成，从而得到合成的轴心轨迹 的表达式 (5)使上述合成轴心轨迹的表达式在Χ、Υ、Ζ三维空间展现即可得到三维合成轴心 轨迹。本发明所提出的基于三维空间轴心轨迹的压缩机轴系故障诊断方法，包括三维滤 波轴心轨迹和三维合成轴心轨迹，考虑了转子轴向振动这个重要诊断信息，将三个方向传 感器获取的振动信号集成起来，能够更为全面、生动的反映转子振动状态，可为检测和诊断 提供更加丰富的信息。同时充分利用多源信息的功能，可以大大提高故障诊断的准确度和 效率。
图1为轴心轨迹测试装置示意图。
图2为轴向位移测量时传感器安装方式示意图。图3为单截面三维轴心轨迹图。图4为二截面三维轴心轨迹图。图5为三维滤波轴心轨迹的实现示意图。图6为三维原始轴心轨迹与滤波轴心轨迹。图7为三维原始轴心轨迹与三维合成轴心轨迹。图8为三维空间轴心轨迹分析方法示意图。
具体实施例方式下面结合附图对本发明的内容作进一步详细说明参见图1所示，利用安装在同一截面内相互垂直的两个电涡流传感器(X方向和Y 方向)和安装在轴向的一个电涡流传感器(Z方向)对轴颈振动进行测量后，得到压缩机转 轴的一组振动信号x(t)、y(t)、ζ (t)，仿照轴心轨迹的构成方法，将三个传感器信号(X方 向、Y方向和Z方向)进行合成即可得到三维空间轴心轨迹。参见图2所示，测量轴向振动信号时传感器的安装方式。在实际的压缩机状态监 测和故障诊断中，不仅要测量压缩机转轴径向的振动位移，其轴向位移也是非常重要的监 测参量之一。因为转子的横向振动往往伴随着轴向的窜动，当机器的状态发生变化时，轴向 位移也会发生变化。所以，在诊断分析中引入轴向振动位移量，对故障的确诊有重要的作用。参见图3所示，设传感器测量得到振动信号分别为x(t)、y(t)和z(t)。设w(t)为 在X-Y-Z三维空间进行组合后的三维轴心轨迹，其表达式为w(t) = ix(t)+jy(t)+kz(t), 表示为向量的形式r =乏。按照轴心轨迹的合成方法，将三个传感器所获得信号的 瞬时幅值一一对应地作为X-Y-Z空间上的坐标点(X，y, z)，随着转子的转动，该点即形成 X-Y-Z空间中的三维轴心轨迹。由三维轴心轨迹的向量表达式，三维轴心轨迹包含着转子三 个方向上的振动信息，形成了三维的空间图形。由于不同的故障在不同的测量方向上有着 不同的反映，所以我们将三维轴心轨迹投影到二维平面上，就能对转子的三维轴心轨迹有 更清晰的认识。图4所示，为了考察转子整体的振动情况，可以显示出转子上多个支承面上的三 维轴心轨迹图，其构造的方法和单截面三维轴心轨迹的方法类似，比如，二截面三维轴心轨
迹可以用向量的形式表示为
式中^工^与^^分别对应于两个
截面上水平方向、垂直方向和轴向方向的信号向量。由此，可以通过观察多个截面的三维轴 心轨迹的合成，来分析转子系统的故障机理。参见图5所示，由于信号采集和传输的过程中周围环境的干扰，模拟信号在传送 的过程中往往夹杂着各种干扰，使采样后的数据产生一定的失真，这影响了振动的检测以 及分析的正确性，因此必须采取一定的滤波手段去除噪声干扰。三维滤波轴心轨迹在构造 上，首先将三个方向的振动信号x(t)、y(t)和z(t)分别进行滤波，然后再将滤波后的信号 x' (n)、y' (η)和ζ' (η)在时域内重新构造出三维轴心轨迹。设信号χ(t)、y (t)和ζ (t) 经低通滤波后的采样量化值为χ (n)、y (η)和ζ (η)，滤波器的冲激响应为h (η)，滤波后的信号为χ' (n)、y' (η)和ζ' (η)，其表达式为
由此可得滤波后的轴心 ζ (n) = z(n) * h(n)
轨迹表达式为
根据轴心轨迹表达式即可构造出三 维滤波轴心轨迹。参见图6所示，左边为某转子的三维原始轴心轨迹，图中的线条繁杂混乱，图形特 征不明显，给轴心轨迹的识别和故障诊断带来困难。而右边则为在低频区带通滤波后的三 维滤波轴心轨迹，从图中可以明显的看到含有噪声的轨迹基本恢复了原貌，滤波后的轴心 轨迹是由很多椭圆叠加而成，线条清晰。参见图7、8所示，一般情况下，转子故障中存在着主次之分，其振动信号中除包含 由不平衡引起的同步振动分量外，还存在频率低于转子转速的亚同步振动分量和频率为转 子转速整数倍的高次谐波振动分量。这就使得转子轴心轨迹形状比较复杂、混乱，给轴心轨 迹的识别带来了很大困难。为了克服原始轴心轨迹难以分析的缺点，我们可以利用频谱分 析原理，将振动信号先分解为一系列的简谐振动分量，再将故障诊断时所需的频率分量复 合起来，就可以通过合成的轴心轨迹为故障识别提供更为突出的特征信息。设轴心轨迹X、 Y和ζ方向的振动信号为x(t)、y(t)和ζ (t)，将这三个方向上的信号分解为各频率分量下 的一个个分量χ (t) = A1Sin ( ω 0t + Φ i) +A2Sin (2 ω 0t + Φ 2) + ··· +Ansin (η ω 0t+ Φ n)=
Xl(t)+X2(t)+-+Xn(t)y (t) = B1Sin ( ω 0t + Ψ i) +B2Sin (2 ω 0t + Ψ 2) + ··· +Bnsin (η ω 0t+ Ψη)=ζ (t) = C1Sin ( ω 0t + ξ i) +C2Sin (2 ω 0t + ξ 2) + ··· +Cnsin (η ω 0t + ξ η)=
Zl(t)+Z2(t)+-+Zn(t)其中，An，Bn，Cn，Φη，Ψη, ξη——水平、垂直和轴向信号的幅值与相位，ω0—信号的旋转角频率，η——频率分量的最高阶数；分别把分量Xi(t)、yi(t)和Zi (t)合成为方程组的形式
X1 (f) = A1 sin(0。i + φ)
(O = A sinK^+^l)
^) = Οι η(ω0 + ξι)
x2(t) = A2 η{2ω0 + φ2) γ2( ) = Β2 η(2ω0 + ψ2)· z2(t) = C2 ^ιη(2ω0 + ξ2)
Χη( ) = Αη^η(ηω0 +Φη) yn(t)=Bn^<nω0t+Ψn)
= sin( cV +么)可以看到先前复杂的轴心轨迹是由不同频率、不同相位以及不同幅值的分量叠加 而成。有了信号的各个振动分量信息，就可以很容易的将各分量合成，从而得到三维合成的 轴心轨迹的表达式 式中，1代表需要合成的频率分量。如图7显示的是某转子三维原始轴心轨迹以 及按1X、2X和4X倍频分量合成后的三维合成轴心轨迹。从图中可以看到，合成后的图形线 条非常清晰，可以利用轴心轨迹的形状识别方法来判断其故障特征。以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施方式
仅限于此，对于本发明所属技术领域的普通技术人员来说，在不脱 离本发明构思的前提下，还可以做出若干简单的推演或替换，都应当视为属于本发明由所 提交的权利要求书确定专利保护范围。
权利要求
一种基于三维空间轴心轨迹的压缩机轴系故障诊断方法，其特征在于(1)根据需要在压缩机转子轴颈处某一截面或同时在某几个截面上安装非接触式电涡流传感器，得到一组振动信号x(t)、y(t)；(2)在压缩机转子轴轴向安装相同的非接触式电涡流传感器，得到振动信号z(t)；(3)将步骤(1)和步骤(2)中得到的振动信号x(t)、y(t)、z(t)进行低通滤波剔除高频干扰信号；(4)设w(t)为x(t)、y(t)、z(t)在X Y Z三维空间进行组合后的三维轴心轨迹，则其表达式为w(t)＝ix(t)+jy(t)+kz(t)；(5)按照轴心轨迹的合成方法，将三个传感器所获得信号的瞬时幅值一一对应地作为X Y Z空间上的坐标点(x，y，z)；随着转子的转动，该点即形成X Y Z空间中的三维轴心轨迹；(6)将三个方向的振动信号x(t)、y(t)和z(t)分别进行滤波，然后再将滤波后的信号x′(n)、y′(n)和z′(n)在时域内按照三维滤波轴心轨迹的构造方法重新构造，得到三维滤波心轨迹；(7)利用频谱分析原理，将振动信号x(t)、y(t)和z(t)分别先分解为一系列的简谐振动分量，再将故障诊断时所需的频率分量复合起来，然后按照三维合成轴心轨迹的构造方法进行构造得到三维合成轴心轨迹。
2.如权利要求1所述基于三维空间轴心轨迹的压缩机轴系故障诊断方法，其特征在 于，所述的三维滤波轴心轨迹的构造方法是(1)设信号χ⑴、y⑴和ζ⑴经低通滤波后的采样量化值为χ(n)、y (η)和ζ (η)，滤波 器的冲激响应为h (η)；(2)对步骤(1)中采样量化后的信号x(n)、y(n)和ζ(η)进行滤波，滤波后的信号为 x' (n)、y' (η)和ζ' (η)，其表达式为χ (n) = x(n) * h{n)\y {η) - y{n) * h{n)(1)ζ (η) = z(n)^h(n)(3)由此可得滤波后的轴心轨迹表达式为W' (n) = ix' (n)+jy' (η) +kz' (η)(2)(4)滤波后轴心轨迹的向量表达式为wF = r+F+zT(3)(5)按照轴心轨迹的合成方法，将滤波后的三组信号χ'(n).y' (η)和ζ' (η)的瞬时 幅值一一对应地作为X-Y-Z空间上的坐标点(X，y, ζ)，随着转子的转动，该点即形成X-Y-Z 空间中的三维轴心轨迹。
3.如权利要求1所述基于三维空间轴心轨迹的压缩机轴系故障诊断方法，其特征在 于，所述的三维合成轴心轨迹的构造方法是(1)设压缩机转轴Χ、γ和Z方向的振动信号分别为x(t)、y(t)和z(t)；(2)将上述三个方向上的信号分解为各频率分量下的一个个分量χ (t) = A1Sin (ω0 +Φ ^ +A2Sin (2 ω 0t+ φ 2) +... +Ansin (η ω 0t+ φ η) 其中，An，Bn，Cn，Φη，Ψη, ξη~水平、垂直和轴向信号的幅值与相位； ω0——信号的旋转角频率； η——频率分量的最高阶数；(3)分别把分量Xi (t)、y, (t)和Zi (t)合成为方程组的形式 (4)根据信号的各个振动分量信息，可以将各分量合成，从而得到合成的轴心轨迹的表达式 (5)使上述合成轴心轨迹的表达式在X、Y、Z三维空间展现即可得到三维合成轴心轨迹。
全文摘要
本发明公开了一种基于三维空间轴心轨迹的压缩机故障诊断方法，(1)检测压缩机转子的径向振动信号x(t)、y(t)；(2)检测压缩机转子的轴向振动信号z(t)；(3)将振动信号进行低通滤波剔除高频干扰信号；(4)获得三维轴心轨迹w(t)；(5)形成三维轴心轨迹；(6)将三个方向的振动信号x(t)、y(t)和z(t)分别进行滤波后在时域内重新构造，得到三维滤波心轨迹；(7)利用频谱分析原理，进行合成得到合成轴心轨迹。本发明的基于三维空间轴心轨迹的压缩机故障诊断方法在全息谱技术的基础上，提出了三维振动分析方法，该方法包括合成三维振动分析与滤波三维振动分析，能够更加清楚、生动地反映设备的运行状态和故障特征。转子的振动包含水平、垂直以及轴向的三个方向振动。将三个方向的振动信号加以合成，构造出三维轴心轨迹，包含三维滤波轴心轨迹以及三维合成轴心轨迹。这种方法能够真实的反映转子轴心的运动情况，为故障诊断提供了更为全面的信息。
文档编号F04B51/00GK101907089SQ201010257659
公开日2010年12月8日 申请日期2010年8月20日 优先权日2010年8月20日
发明者吴婷婷, 张西宁, 温广瑞, 田利凯, 臧廷朋 申请人:西安交通大学