汽轮机转子动静碰摩故障诊断技术探究.doc

汽轮机转子动静碰摩故障诊断技术研究 1、相关定义 1.1、离散小波变换定义 与 Fourier 分析中的离散 Fourier 变换类似,离散小波变换也是用离散点的和式来代 替连续点的积分,离散点之间的间隔将决定离散变换对连续变换近似的精度。在实际应 用中一般将尺度参数a以二分法离散化,即只取 2 的整数次幂,而位移参数b 可在 Z 域 (即整数域)上等间隔取值,离散小波正变换(Discrete Wavelet Transform-DWT)的定 义如下: DWTfttdt mnZ mn()mn( ) , ,, ∑ ∈ = ψ (3-10) 写成内积形式有 DWTm,n=f(t) ,ψ m ,n (t) m,n z (3-11) 1.2、频谱分析的基本概念 频谱分析是信号分析处理的一个重要环节。信号频谱是信号在频率域中的表现方 式,反映了信号的幅度、相位、能量、功率等信息在各频率上的分布。谱分析是进行故 障诊断的依据,通过计算有限长度序列的离散傅立叶变换即快速傅立叶变换,得到振动 信号的频谱,在频域内研究振动特性。波形分析是对信号在时间域上,即各种物理量的 幅值在时间坐标上的分布进行分析,如幅值随时间的的变化,起始与结束的时间,波形 的畸变等。 信号的波形分析和频谱分析是相互独立又密切相关的,它们之间可以通过傅立叶变 换相互转换。 波形分析是对信号的时间、幅值和周期的测量,除相关分析以外,一般使用通常的 测量标度,其方式较简单迅速,表现更直观。而频谱分析是对动态信号在频率域内的分 布进行分析,得到的是以频率为横坐标的各种物理量的谱线和曲线,即各种以频率为变 量的频谱函数F(? ) 。频谱分析中可求得幅值谱、相位谱、功率谱和各种谱密度等等。 频谱分析过程较为复杂,它是以傅立叶级数和傅立叶积分为基础的。 频域和时域分析是动态信号分析的两个观察角度。波形分析中以时间为坐标表示动 态信号;频谱分析则是以频率为坐标表示动态信号。时域表示较为直观形象,而频域表 示动态信号则更能反映信号深层次的特性。目前,动态信号分析的趋势是从时域向频域 发展,但同时这两种方法仍是互相联系、相辅相成,被同时应用于信号特性的研究中 。 频谱分析和波形分析既有悠久的历史,同时又随着科技的发展而不断涌现出新方 法。目前,在动态信号的研究日益广泛的情况下,已经从确定性信号发展到随机信号处 理分析的阶段,对于随机数据的处理还涉及到概率论和数理统计方法[16]。 现在频谱分析和波形分析与随机数据处理分析方法,已成了现代科学研究的有力工 具,在各个科学研究领域,例如航空航天、机械、土木、地震与资源勘测等领域的各种 结构的动态实验分析中,都得到了应用。目前,西方国家的理工专业的学生,都必须学 习这门课程,一个受过良好教育的工程师和科学家,需要至少对这方面的概念和方法有 1.3、小波包的定义 设支撑长度为 L 的正交镜像滤波器(Quadrature Mirror Filter-QMF)组h(n ) 、g(n ) 满 足以下条件: (1)∑ 0 (由于h(n ) 的支撑长度 L 是有限的,故其显然是满足 的); (2)∑?? = n h(n2k)h(n2l) δk ,l ; (3)∑ = n h(n) 2 且g(k)=(?1) kh(L?1 ? k ) 。其中,h(n ) 、g(n ) 分别对应低通和高通滤波,定义 l2(z)→l2(2z)的算子 H 和G 如下: { } ∑ =? kz Hs(k)(i)s(k)h(k2i ) (3-34) { } ∑ ∈ =? kz Gs(k)(i)s(k)g(k2 i) (3-35) 其中s(k ) 代表输入序列,H {s(k )}( i) ,G {s(k) } ( i) 表示s(k ) 经过算子作用后的输出。 它们是一对反卷积-下采样的过程。与算子H 和G 相对应的共轭算子H ? 和G ? 定义为 { } ∑ ?= ? jz Hs(j)(k)s(j)h(k2i ) (3-36) { } ∑ ∈ ?= ? jz Gs(j)(k)s(j)g(k2 i) (3-37) 其中算式中的变量含义与上同,这对算子是上采样-卷积过程。有时为了表示方便, 将算子分别写成H 1 和H 2 。算子H 、G 及它们的共轭H ? 、G ? 之间必须满足完备的重构 条件(Perfect Reconstruction Condition-PRC): ?? ??? == == += ?? ?? ?? HHGGI HGGH HHGGI 0 (3-38) 下面利用递归的方式定义函数序列 1.4、关联维数定义 X是重构的相空间,X中的每一行都是m维相空间的一个向量。现将它们任 意两者之间的差值的绝对值记为ri, j =|xi-x j |,xi和x j 为m维