[发明专利]变速激励下叶片异步振动参数识别方法

权利要求书

1.一种变速激励下叶片异步振动参数识别方法，其特征在于：包括以下步骤：

1)对叶片的相邻两个振动信号进行插值处理；

2)通过全相位傅里叶变换获得叶片振动信号的混叠频率Δω和实测相位差ΔΦ^*，求得振动阶次的非整数部分Δm；

3)利用测量振动信号相位与传感器安装角的关系，通过阶次跟踪获得振动阶次N_e，N_e＝(m+Δm)，m是振动阶次的整数部分，Δm是振动阶次的非整数部分；

4)得到叶片的振动频率ω；

所述步骤1)振动信号的插值处理如下：

以第一传感器测量的叶片的振动数据TIP(i)＝[x_i1,x_i2,x_i3,…x_in]^T为基准将第二传感器测量的该叶片数据顺序插入该叶片的振动数据序列中：第二传感器采集的数据表示为:TIP(i+1)＝[x_(i+1)1,x_(i+1)2,x_(i+1)3,…x_(i+1)n]^T,将第二传感器采集的数据插入到第一传感器采集的数据中，插值完成后该传感器采集的数据变为TIP(i)＝[x_i1,x_(i+1)1,x_i2,x_(i+1)2,x_i3,x_(i+1)3,…x_in,x_(i+1)n]^T。

2.根据权利要求1所述的变速激励下叶片异步振动参数识别方法，其特征在于：所述步骤3)通过阶次跟踪获得振动阶次包括以下步骤：

31)获得跟踪相位差ΔΦ_k；

32)求得实测相位差ΔΦ^*和跟踪相位差ΔΦ_k的残值e_k，e_k＝ΔΦ_k-ΔΦ^*；

33)求得残值e_k的均方根值S_k，

34)分别计算Δm为“+”、“-”时的S_k，比较大小，并取S_k为最小时对应的m以及对应的Δm前面的“+”、“-”。

3.根据权利要求1所述的变速激励下叶片异步振动参数识别方法，其特征在于：所述4)中得到叶片的振动频率ω的过程为：将步骤3)获得的振动阶次(m+Δm)以及相应的振动倍频Ω_n代入ω＝(m±Δm)Ω_n求得振动频率ω。