专利名称：基于自相关的大型轴流风机不平衡量的识别方法
技术领域：
本发明涉及转动件的动平衡技术，尤其涉及一种基于自相关的大型轴流风机不平
衡量的识别方法。
背景技术：
大型轴流风机在使用过程中容易受到不平衡引起的振动困扰，严重影响机组的使用寿命和周围工作环境。大型轴流风机在生产过程中，虽然在动平衡机上对叶轮进行了动平衡，但往往因为装配误差，和叶片安装运输过程中产生碰撞或变形导致新的不平衡振动。大型轴流风机电机被封闭在机壳内，平衡所需的键相信号只能通过在叶轮上贴反光片用光电头来获得，但是轴流风机工作时，轴向风力很大，实际上很难正确获得键相信号，且十分危险。从而很难采用常规的测相测幅法进行现场动平衡测试，如何解决该类旋转机械的不平衡振动问题具有十分重要的现实意义。
发明内容
本发明的目的是克服现有技术的不足，提供一种基于自相关的大型轴流风机不平衡量的识别方法。 基于自相关的大型轴流风机不平衡量的识别方法包括如下步骤 1)假设轴流风机振动故障由不平衡量引起，即轴流风机振动信号中工频成分占轴
流风机振动信号总能量的一半以上，定义轴流风机振动信号中工频信号的幅值与其他频率
信号的幅值之比大于等于2，采样长度大于转速信号周期，其中轴流风机振动信号包括工频
成分、半倍频成分、二倍频成分、其他频率成分及噪声信号，不平衡量引起的振动对应轴流
风机振动信号的工频成分； 2)对轴流风机振动信号进行低通滤波，抑制高频信号； 3)对经过低通滤波后的轴流风机振动信号进行二重自相关，从轴流风机振动信号众多频率成分中分离出工频成分； 4)对二重自相关后的轴流风机振动信号进行有效值的计算，还原出轴流风机振动信号中工频成分的幅值； 5)用测幅法进行动平衡，确定实际轴流风机不平衡质量的大小和位置。 2.根据权利要求1所述的一种基于自相关的大型轴流风机不平衡量的识别方法，
其特征在于所述的假设轴流风机振动故障由不平衡量引起，即轴流风机振动信号中工频成
分占轴流风机振动信号总能量的一半以上，定义轴流风机振动信号工频信号的幅值与其他
频率信号的幅值之比大于等于2，采样长度大于转速信号周期步骤为 大型轴流风机引起振动信号用数学式表示式为 <formula>formula see original document page 4</formula>
式中，A-工频振动信号的振幅；f_工频；13 _工频振动信号的相位； Bi-其他频率振动信号的振幅；Vi-其他振动信号的频率；n i_其他频率振动信号
的相位；n(t)-标准的高斯白噪声， 假设轴流风机振动故障由不平衡量引起，即轴流风机振动信号工频成分占轴流风 机振动信号总能量的一半以上，定义轴流风机振动信号工频信号的幅值与其他频率信号的 幅值之比A/Bi > 2 ;其次，采样长度远大于转速信号的周期。 3.根据权利要求1所述的一种基于自相关的大型轴流风机不平衡量的识别方法，
其特征在于所述的对轴流风机振动信号进行低通滤波，抑制高频信号步骤为 轴流风机的转速小于6000转/分，即不平衡引起的振动信号的工频分量频率在
100Hz以下，引入二阶RC低通滤波器，将其截止频率设置在lOOHz，复杂振动信号经过上述
低通滤波器后可抑制振动信号的中的高频分量，而让不平衡引起的工频分量通过。 4.根据权利要求1所述的一种基于自相关的大型轴流风机不平衡量的识别方法，
其特征在于所述的对经过低通滤波后的轴流风机振动信号进行二重自相关，从轴流风机振
动信号众多频率成分中分离出工频成分步骤为 轴流风机振动信号经过低通滤波后，衰减了高频信号分量，依然存在跟转速非同 频的其他振动分量和干扰分量，为提取工频分量，对低通滤波后的振动信号进行自相关运 算，这样工频信号的幅值与其他频率信号幅值之比增大到原来的2次方，一次自相关的后 轴流风机振动信号表示式为
Rxl00 = y"r戶丁sin(2;r/V) + S十sin(2;n^) + (2)
2 i=i 2 式中，A-工频振动信号的振幅；f-工频；Bi-其他频率振动信号的振幅；Vi-其他频 率频率t )-标准的高斯白噪声自相关表达式， 5.根据权利要求1所述的一种基于自相关的大型轴流风机不平衡量的识别方法， 其特征在于所述的对二重自相关后的轴流风机振动信号进行有效值的计算，还原出轴流风 机振动信号中工频成分的幅值步骤为 进一步把yi ( t )当做输入信号，进行二重自相关，得到轴流风机振动信号的二重
自相关函数表达式
j4 m B 4
Rx2 0)= 丁 cos(2;r/r) + X!十cos(2呻)+ "2 (r) (3) 8 i=i 8 式中n2( O-标准的高斯白噪声二重自相关函数， 忽略噪声信号，轴流风机振动信号二重自相关函数表示式为
爿4
^m二重自相关0) = i cos(^;r/V) (斗) 对轴流风机振动信号的数字采样值进行二重自相关并进行有效值的计算，得到轴 流风机原始振动信号中工频信号的幅值A :
5
式中，U
_轴流风机振动信号的数字采样值的二重自相关表达式。 6.根据权利要求1所述的一种基于自相关的大型轴流风机不平衡量的识别方
法，其特征在于所述的用测幅法进行动平衡，确定转子的原始不平衡量为M。，不平衡相位为
e。，步骤为 将一个已知试重M安装在转子校正面的角度9i处，i = 1，2，3，开动轴流风机，测
得相应的幅值为Ai，i = 1，2，3，使用余弦定理和线性系统假设推导出原始不平衡量M。及相
位9 。的表达式为
u M
0 4 、o 其中，S == —1z
2
3< ° 本发明首次提出将自相关算法应用到对轴流风机振动信号的提取上，通过二重自 相关法有效的提取出了信号中的工频部分，通过算术平均值的方法计算出二重自相关后信 号的幅值，最后还原出原始信号工频部分的幅值，采用侧幅测动平衡的方法，确定所需加平 衡质量及角度。有效地解决了大型轴流风机无法提取转速信号的难题，避免了需要通过测 得振动信号的相位信息而推导出不平衡质量及角度，并通过实验数据验证了该方法的可行 性，有较大的工程应用价值。
图1是轴流风机原始振动信号时频域图；
图2是一次自相关振动信号时频域图；
图3是二重自相关振动信号时频域图；
图4是测幅测振识别模型图。
具体实施例方式
基于自相关的大型轴流风机不平衡量的识别方法包括如下步骤 1)假设轴流风机振动故障由不平衡量引起，即轴流风机振动信号中工频成分占轴
流风机振动信号总能量的一半以上，定义轴流风机振动信号中工频信号的幅值与其他频率
信号的幅值之比大于等于2，采样长度大于转速信号周期，其中轴流风机振动信号包括工频
成分、半倍频成分、二倍频成分、其他频率成分及噪声信号，不平衡量引起的振动对应轴流
风机振动信号的工频成分； 2)对轴流风机振动信号进行低通滤波，抑制高频信号；
3)对经过低通滤波后的轴流风机振动信号进行二重自相关，从轴流风机振动信号 众多频率成分中分离出工频成分； 4)对二重自相关后的轴流风机振动信号进行有效值的计算，还原出轴流风机振动 信号中工频成分的幅值； 5)用测幅法进行动平衡，确定实际轴流风机不平衡质量的大小和位置。 所述的假设轴流风机振动故障由不平衡量引起，即轴流风机振动信号中工频成分
占轴流风机振动信号总能量的一半以上，定义轴流风机振动信号工频信号的幅值与其他频
率信号的幅值之比大于等于2，采样长度大于转速信号周期步骤为 大型轴流风机引起振动信号用数学式表示式为
附
:c[t] = j sin(2;r力+ 〃) + sin(2;rv,t + /7,.) + "(t) (1)
!■=1 式中，A-工频振动信号的振幅；f-工频；13 _工频振动信号的相位；B「其他频率 振动信号的振幅；Vi-其他振动信号的频率；ni-其他频率振动信号的相位；n(t)-标准的 高斯白噪声， 假设轴流风机振动故障由不平衡量引起，即轴流风机振动信号工频成分占轴流风
机振动信号总能量的一半以上，定义轴流风机振动信号工频信号的幅值与其他频率信号的
幅值之比A/Bi > 2 ;其次，采样长度远大于转速信号的周期。 所述的对轴流风机振动信号进行低通滤波，抑制高频信号步骤为 轴流风机的转速小于6000转/分，即不平衡引起的振动信号的工频分量频率在
100Hz以下，引入二阶RC低通滤波器，将其截止频率设置在100Hz，复杂振动信号经过上述
低通滤波器后可抑制振动信号的中的高频分量，而让不平衡引起的工频分量通过。 所述的对经过低通滤波后的轴流风机振动信号进行二重自相关，从轴流风机振动
信号众多频率成分中分离出工频成分步骤为 轴流风机振动信号经过低通滤波后，衰减了高频信号分量，依然存在跟转速非同
频的其他振动分量和干扰分量，为提取工频分量，对低通滤波后的振动信号进行自相关运
算，振动信号的自相关函数为
1^0")=五00)"(> + "] =- fx昨("r)(i/ (2)
71 — r J J& o 假设x (t) = s (t) +n (t)，其中s(t)=爿sin(2;r々+ y ) + Z A sin(2;rv,. +《.)，
! =1
信号的自相关函数变为 Rx1(t) =E[s(t) *s(t+t)] =Rss(T)+E[s(t) *n (t+t ) ]+E [s (t+t ) *n(t)]+Rm
(t ) (3) 上式第一项为信号的自相关函数Rss( t )，第二、三项为信号和噪声的互相关函数
E[s(t) .n(t+"]和E[s(t+" "(t)]，第四项为噪声的自相关函数IC(0。 利用信号和噪声、噪声和噪声的不相关特性，当采样长度远大于周期信号的周期
7时，E[s(t) .n(t+0]，E[s(t+0 .n(t)]和Rm("的值均可忽略不计。 根据不同频不相关的原理可知，信号的自相关函数Rss( O相当于各自频率的自
相关函数的相加，工频信号的自相关函数可以表示为
《,工频(0 = ;^" f'x(t)x(/ +力^ = A f' Asin(2;r// + 〃)Asin[2;7"/(/ + r) +
1 —T "° i 一T (4)
」21 J2
=—cos(2;r/V)-- —[ cos[2;r/(2"r) + 2州f 由于采样长度T远大于转动信号的周期，^,4" f—、os[2;r/(2, + r) + 2々]^的
T一r 2力
值可忽略不计，故可以得出&sIS(r) = f cos(2;r/r)，同理其它频率的自相关函数表达式
为
、m' b 2
"威他频率0") = X!十cos—if) (5)
i=l Z —次自相关后，轴流风机振动信号表示式变为
m b 2
r xi 0) = y t ( r)= 7 cos(2;r/V) + Z十cos(2，) + (r) (6)
2 i=i 2 式中，A-工频振动信号的振幅；f-工频；Bi-其他频率振动信号的振幅；Vi-其他频 率频率；^(t)-标准的高斯白噪声自相关函数， 这样，可以进一步把yji )当做输入信号，再做自相关运算，得到轴流风机振动信
号的二重自相关函数表达式
rx2 0)=T cos(2;r/r) + Z t cos(2，) + "2 (r) (7) 8 i=i 8 式中n2( O-标准的高斯白噪声二重自相关函数， 由于工频信号的幅值与其他频率信号的幅值之比A/Bi > 2，当经过二次自 相关运算后，工频信号的幅值与其他频率信号的幅值之比为AVB^ > 16，因此可将
S"^cos(2;n^)忽略，同时标准的高斯白噪声的二重自相关函数112")也可以忽略，这 i=i 8
样轴流风机振动信号的二重自相关函数表达式可表示为
^二重自被")^ + cos(2;rA) (8) 令t = 0， (8)式变为
凡二廳关(0) = + (9) 所述的对二重自相关后的轴流风机振动信号进行有效值的计算，还原出轴流风机 振动信号中工频成分的幅值步骤为 由于实际测试中，只能得到信号x[t]的N个采样值x[n]，所以必须通过这N个采 样值来估计信号的一次自相关函数r[m]和二重自相关函数U[m]。计算r[m]的直接估计算
法为
<formula>formula see original document page 9</formula>
(10)
<formula>formula see original document page 9</formula> 令m二 o，式(io)变为r
二^J]x(n)2 ，则轴流风机振动信号的数字采样值 的二重自相关表达式为
<formula>formula see original document page 9</formula> (9) (11)两式联立求得轴流风机原始振动信号中工频信号的幅值A :
<formula>formula see original document page 9</formula>(12) 采用虚拟仪器Labview进行二重自相关仿真实验，假设原始振动信号为 Asin (2 Ji f。t) +Bsin (2 ji f\t) +C sin (2 ji f2t) +n (t)，主要由半倍频、工频、二倍频及噪声信 号构成。其中工频f。 = 10Hz，幅值A = IV ;半倍频分量^ = 5Hz，B = 0. 5V ;二倍频分量f2 =20Hz， C = 0. 3V ;n(t)为高斯白噪声。信号的采样频率为1000Hz，采样点数为10000个点。 对仿真结果进行分析，从图1可以看出原始信号由于夹杂噪声及其他倍频信号的 干扰，由不平衡引起的工频信号很难从中提取出来。从频域图也可以看出，工频成分的幅值 与半倍频信号的幅值之比为2，与二倍频的幅值之比为10/3。经过自相关后，虽然相位信 息丢失了，但是频率并未发生变化，此时的幅值也和原来的幅值之间存在解析关系，通过频 谱分析和相应的计算也可以得到原信号的频率信息和幅值信息。从图2可以看出，一次自 相关后工频成分的幅值与半倍频信号的幅值之比为4，与二倍频的幅值之比约等于ll，有 用工频信号的分量相对于其他干扰分量的信噪比有了明显提高。从图3可以看出，经过二 重自相关后的工频成分的幅值与半倍频信号的幅值之比为16，与二倍频的幅值之比约等于 123，不管从时域还是频域都可以明显的看出信号主要由工频成分构成，然后根据解析关系 就可以非常容易的拾取出原始信号中不平衡引起的工频分量，仿真实验也很好的验证此方 法的现实可行性。 所述的一种基于自相关的大型轴流风机不平衡量的识别方法，其特征在于所述的用测幅法进行动平衡，确定实际轴流风机不平衡质量的大小和位置步骤为
如
中图4所示，首先设轴流风机转子的原始不平衡量为M。，不平衡相位为 9 。，由不平衡重引起的振动信号的幅值为A。，为了确定M。和e 。的值，需使用一个已知试重 M，将它安装在转子校正面的角度9 i处，i = 1，2，3，开动轴流风机，测得相应的幅值为Ai，i =1，2，3，一般来说ei可取任意值，但若取2Ji圆周上的等分点，可大大简化计算，这里取 ei= (i-l)2Ji/3。针对图4，使用余弦定理得 A/ = A02+A2+2AA0cos ( 9 「 9 0) (i = 1， 2， 3) (13)
并根据线性系统假设可得
<formula>formula see original document page 10</formula>
根据(14)，推导出原始不平衡j M
3
M。及相位e 0的表达式<formula>formula see original document page 10</formula>
通过以上计算，理论上准确得出了风机叶轮的不平衡质量M。及其所处的角度位置
权利要求
一种基于自相关的大型轴流风机不平衡量的识别方法，其特征在于包括如下步骤1)假设轴流风机振动故障由不平衡量引起，即轴流风机振动信号中工频成分占轴流风机振动信号总能量的一半以上，定义轴流风机振动信号中工频信号的幅值与其他频率信号的幅值之比大于等于2，采样长度大于转速信号周期，其中轴流风机振动信号包括工频成分、半倍频成分、二倍频成分、其他频率成分及噪声信号，不平衡量引起的振动对应轴流风机振动信号的工频成分；2)对轴流风机振动信号进行低通滤波，抑制高频信号；3)对经过低通滤波后的轴流风机振动信号进行二重自相关，从轴流风机振动信号众多频率成分中分离出工频成分；4)对二重自相关后的轴流风机振动信号进行有效值的计算，还原出轴流风机振动信号中工频成分的幅值；5)用测幅法进行动平衡，确定实际轴流风机不平衡质量的大小和位置。
2. 根据权利要求1所述的一种基于自相关的大型轴流风机不平衡量的识别方法，其特征在于所述的假设轴流风机振动故障由不平衡量引起，即轴流风机振动信号中工频成分占轴流风机振动信号总能量的一半以上，定义轴流风机振动信号工频信号的幅值与其他频率信号的幅值之比大于等于2，采样长度大于转速信号周期步骤为大型轴流风机引起振动信号用数学式表示式为加 <formula>formula see original document page 2</formula>式中，A-工频振动信号的振幅；f_工频；13 _工频振动信号的相位；Bi-其他频率振动信号的振幅；Vi-其他振动信号的频率；ni-其他频率振动信号的相位；n(t)-标准的高斯白噪声，假设轴流风机振动故障由不平衡量引起，即轴流风机振动信号工频成分占轴流风机振动信号总能量的一半以上，定义轴流风机振动信号工频信号的幅值与其他频率信号的幅值之比A/Bi > 2 ;其次，采样长度远大于转速信号的周期。
3. 根据权利要求1所述的一种基于自相关的大型轴流风机不平衡量的识别方法，其特征在于所述的对轴流风机振动信号进行低通滤波，抑制高频信号步骤为轴流风机的转速小于6000转/分，即不平衡引起的振动信号的工频分量频率在100Hz以下，引入二阶RC低通滤波器，将其截止频率设置在100Hz，复杂振动信号经过上述低通滤波器后可抑制振动信号的中的高频分量，而让不平衡引起的工频分量通过。
4. 根据权利要求1所述的一种基于自相关的大型轴流风机不平衡量的识别方法，其特征在于所述的对经过低通滤波后的轴流风机振动信号进行二重自相关，从轴流风机振动信号众多频率成分中分离出工频成分步骤为轴流风机振动信号经过低通滤波后，衰减了高频信号分量，依然存在跟转速非同频的其他振动分量和干扰分量，为提取工频分量，对低通滤波后的振动信号进行自相关运算，这样工频信号的幅值与其他频率信号幅值之比增大到原来的2次方，一次自相关的后轴流风机振动信号表示式为<formula>formula see original document page 3</formula>式中，A-工频振动信号的振幅；f-工频；Bi-其他频率振动信号的振幅；Vi-其他频率频率( T )-标准的高斯白噪声自相关表达式。
5.根据权利要求1所述的一种基于自相关的大型轴流风机不平衡量的识别方法，其特征在于所述的对二重自相关后的轴流风机振动信号进行有效值的计算，还原出轴流风机振动信号中工频成分的幅值步骤为进一步把yi ( T )当做输入信号，进行二重自相关，得到轴流风机振动信号的二重自相关函数表达式 <formula>formula see original document page 3</formula>式中n2( O -标准的高斯白噪声二重自相关函数，忽略噪声信号和其它频率信号，轴流风机振动信号二重自相关函数表示式为 <formula>formula see original document page 3</formula>4)对轴流风机振动信号的数字采样值进行二重自相关并进行有效值的计算，得到轴流风机原始振动信号中工频信号的幅值A : <formula>formula see original document page 3</formula>)式中，u [o]-轴流风机振动信号的数字采样值的二重自相关表达式。
6.根据权利要求1所述的一种基于自相关的大型轴流风机不平衡量的识别方法，其特征在于所述的用测幅法进行动平衡，确定转子的原始不平衡量为l，不平衡相位为e 。，步骤为将一个已知试重M安装在转子校正面的角度9i处，i = 1，2，3，开动轴流风机，测得相应的幅值为Ai，i二l，2，3，使用余弦定理和线性系统假设推导出原始不平衡量M。及相位e。的表达式为 <formula>formula see original document page 3</formula>其中，《=
全文摘要
本发明公开了一种基于自相关的大型轴流风机不平衡量的识别方法。方法1)假设轴流风机振动由不平衡量引起，定义轴流风机振动信号工频信号的幅值与其他频率信号的幅值之比大于等于二，采样长度大于转速信号周期；2)对轴流风机振动信号进行低通滤波；3)进行二重自相关，从轴流风机振动信号众多频率成分中分离出工频成分；4)进行有效值计算，还原出轴流风机振动信号中工频成分幅值；5)用测幅法进行动平衡，确定轴流风机不平衡质量的大小和位置。本发明可以在不提取键相信号的前提下，从复杂振动信号中有效分离出不平衡振动幅值分量，最终通过测幅平衡法找出轴流风机的不平衡大小和位置，有较大的工程应用价值。
文档编号F04D19/00GK101749256SQ20101003960
公开日2010年6月23日 申请日期2010年1月8日 优先权日2010年1月8日
发明者张倩, 张志新, 钟越波 申请人:浙江大学