专利名称：旋转机械启停车过程波德图/奈奎斯特图的快速计算方法
技术领域：
本发明属于旋转机械状态监测和故障诊断领域中的非平稳振动信号处理领域；特别涉及一种用于旋转机械启停车非平稳过程波德图/奈奎斯特图的快速计算方法。
背景技术：
旋转机械启停车是机组运行过程中的事故多发阶段。在这个过程中机组经历了一个非稳态的运行过程，这个过程包含了机组丰富的状态信息，而这些信息是常规稳态状态下所无法获得的。研究旋转机械启停车过程振动信息的提取和利用，对了解旋转机械的运行状态和旋转机械隐含的故障具有重要的指导作用。对旋转机械启停车过程振动信息的提取和利用寻找一种好的分析方法同样具有重要的意义，本发明意义在于快速、准确的得到 旋转机械启停车过程振动信号的一些特征信息，为旋转机械的动平衡提供指导。旋转机械启停车故障诊断的核心是对其进行故障特征的提取，它直接关系到故障诊断的准确性和故障预报的可靠性。据相关资料分析认为，目前常常采用的旋转机械故障诊断方法大都是在恒转速状稳态下进行的，而旋转机械启停车过程是机组运行过程中的事故多发阶段。因此，旋转机械稳态工作状态下，运用快速傅立叶(FFT)变换可直接将振动信号变换到频域，得到各个频率下的幅值分量和相位分量，进而得到频谱图、功率谱图或轴心轨迹图等。以上方法不能用于非稳态过程的振动信号分析，因此，寻找针对旋转机械启动停车过程的非稳态振动信号的快速计算方法具有重要意义。常用的旋转机械非平稳过程振动信号监测诊断方法有波德图、奈奎斯特图和瀑布图等，而波德图、奈奎斯特图对旋转机械启停车过程的分析尤为重要。波德图和奈奎斯特图的提出最早是在自动控制领域，用来描述控制系统不同频率下系统输入与系统输出之间的关系。后来这一方法被移植和引入到旋转机械非平稳起停车过程的监测和诊断中。如果将旋转机械失衡质量在转动过程中所引起的离心力作为系统的输入，将转子的失衡振动响应作为输出，那么转子的振动可认为是一典型的二阶质量-弹簧一阻尼系统。在旋转机械启停车这一非平稳过程中，转子的转速经历了有由小到大或相反的变化过程，由失衡质量产生的离心力的频率也同样经历了由小到大或相反的变化过程。与旋转机械启停车这一非平稳过程相对应的振动信号，正是失衡质量产生离心力所激发的振动。因此振动信号的波德图和奈奎斯特图可以很好地描述由失衡质量产生离心力所激发的振动。旋转机械启停车过程振动信号波德图和奈奎斯特图是目前旋转机械非平稳振动信号处理中重要和常用的方法。目前常用的旋转机械启停车过程振动信号波德图和奈奎斯特图计算方法有两种
O分段采集启停车过程振动信号并进行后处理得到波德图和奈奎斯特图；2)连续采集启停车过程振动信号，通过跟踪滤波的方法计算波德图和奈奎斯特图。第一种方法需要对采集到的每一段振动信号做快速傅立叶(FFT)变换，获取每一段的IX振动峰值和相位，然后绘制振动信号的波德图和奈奎斯特图，但是绘制的波形失真度比较大，并且曲线不光滑，段与段之间的转速间隔不同，波德图和奈奎斯特图的准确度也不同，间隔越小准确度越高，临界转速越接近真实值，但是间隔小计算量就会变得非常大；第二种方法，对采集到的启停车非平稳振动信号用跟踪滤波的方法除去低频和高频干扰得到转频振动信号，然后根据键相信号的特点寻找转频信号的振动周期，求得每个周期内的振动峰值以及初相位，最后绘制振动信号的波德图和奈奎斯特图。该方法中计算每个周期内的振动峰值和初相位有两种方法一种是通过对该周期为中心、左右一段时间的振动信号进行傅立叶变换，通过傅立叶频谱的到幅值和相位信息。该方法需要重复多次进行傅立叶变换，计算量大。另外，进行傅立叶变换的信号只能近似认为是幅值和相位稳定的信号，但实际上信号的幅值和相位是不断变化的，由此导致幅值和相位计算出现误差；另一种是直接对每个周期内振动信号进行幅值搜索，得到振动幅值。根据每个周期内振动信号采样总点数及极大值点的位置来估算相位。这种相位计算方法对低转速时的相位估计比较准确，在高转速时相位计算的误差较大。
发明内容
本发明的目的在于提供一种用于旋转机械启停车过程非平稳振动信号波德图(Bode plot)、奈奎斯特图(Nyquist plot)的快速计算方法。该方法基于由电润流传感器测量到的旋转机械启停车过程转子振动信号和键相信号，快速计算得到启停车过程的波德图和奈奎斯特图。整个计算过程基于对以上启停车过程中振动信号的Hilbert (希尔伯特)变换处理，并结合振动信号的保相滤波处理、键相信号中周期信息的提取，实现了对旋转机械启停车振动信号波德图和奈奎斯特图的快速计算。旋转机械启停车过程波德图/奈奎斯特图的快速计算方法，包括以下步骤I)采集旋转机械启停车过程的键相信号和振动信号；2)对步骤I)采集的振动信号进行带通滤波，获得旋转机械启停车过程转频振动信号;3)用键相信号相邻的峰值去截取启停车过程转频振动信号，得到振动信号的每一个周期点值；4)对启停车过程转频振动信号作希尔伯特变换，绘制变换后的信号幅值包络波形曲线；5)根据启停车过程转频振动信号的周期点值，得到信号每个周期内的最大值和最小值，然后对最大值和最小值分别作插值、平均处理，即得到振动峰值曲线的波德图；6)对启停车过程转频振动信号作希尔伯特变换，绘制变换后的希尔伯特相位波形曲线，根据相位波形曲线寻找周期点处的相位值，作差值处理后，得到真实的傅立叶初相位变化曲线的波德图。本发明进一步的改进在于根据步骤5)获得的振动峰值曲线和步骤6)获得的傅 立叶相位变化曲线，绘制出奈奎斯特图。本发明进一步的改进在于步骤I)中通过电涡流位移传感器采集旋转机械启停车过程的键相信号和振动信号。本发明进一步的改进在于所述电涡流位移传感器采样频率为Fs=24000Hz。本发明进一步的改进在于测量键相信号的电涡流位移传感器与水平面平行；测量振动信号的电涡流位移传感器与水平面面呈45度夹角。本发明进一步的改进在于步骤2)中对步骤I)采集的启停车振动信号进行保相滤波，在每段信号转频左边3Hz处，右边3Hz处各选取一点，对信号作带通滤波，然后合并滤波后的每一段信号，得到启停车过程转频振动信号。旋转机械启 停车过程波德图/奈奎斯特图的快速计算方法，包括以下步骤I)采集旋转机械启停车过程的键相信号和振动信号；2)用键相信号相邻的峰值去截取启停车过程原始振动信号，得到振动信号的每一个周期点值；3)对启停车过程原始振动信号作希尔伯特变换，绘制变换后的信号幅值包络波形曲线；4)根据启停车过程原始振动信号的周期点值，得到信号每个周期内的最大值和最小值，然后对最大值和最小值分别作插值、平均处理，即得到振动峰值曲线的波德图；5)对启停车过程原始振动信号作希尔伯特变换，绘制变换后的希尔伯特相位波形曲线，根据相位波形曲线寻找周期点处的相位值，作差值处理后，得到真实的傅立叶初相位变化曲线的波德图。本发明进一步的改进在于根据步骤5)获得的振动峰值曲线和步骤6)获得的傅立叶相位变化曲线，绘制出奈奎斯特图。本发明提出的快速计算方法的主要特征有(I)如果要计算通频带振动信号的波德图和奈奎斯特图可跳过第(2)步，直接进行第⑶步处理。(2)对旋转机械启停车振动信号进行保相滤波处理，得到单一的启停车过程转频信号。(3)对⑴或⑵中的启停车过程振动信号进行Hilbert变换，得到旋转机械启停车过程中振动信号各个时刻点处振动幅值的变化曲线和振动信号各个时刻点处的相位变化曲线。(4)对(3)中得到的启停车过程振动幅值变化曲线特征点进行提取、插值等处理，得到旋转机械启停车过程振动的幅值变化曲线。(5)提取启停车过程中键相信号的特征点，结合旋转机械启停车过程中振动信号各个时刻点处的相位变化曲线，获取键相信号的特征点处振动信号的相位信息并进行插值处理，得到启停车过程中振动信号各个时刻点处的真实相位变化曲线。(6)根据(4)中计算的幅值变化曲线以及(5)中计算的真实相位变化曲线，绘制启停车振动信号波德图和奈奎斯特图。本发明中的快速计算方法，具体来讲有以下几个步骤I)根据实验平台自身对振动信号的分析要求，本发明中所使用的数据是由电涡流传感器采集得到，并且要确定合适的键相传感器的安装角度、振动传感器的安装角度(见
图1-1和图1-2)，采集转子在启停车过程中振动数据和转子的键相数据，该算法只针对其中任一个振动数据(见图2);键相数据见图3。2)根据键相信号的特点，获取键相信号峰值对应点的数据，见图4。3)对旋转机械启停车振动信号进行保相滤波处理，得到单一的启停车过程转频信号x(t)，见图5。4)对保相滤波后的转频振动信号作Hilbert变换后，得到信号y (t)，构造解析函数Z(t)，满足:
权利要求
1.旋转机械启停车过程波德图/奈奎斯特图的快速计算方法，其特征在于，包括以下步骤 1)采集旋转机械启停车过程的键相信号和振动信号； 2)对步骤I)采集的振动信号进行带通滤波，获得旋转机械启停车过程转频振动信号； 3)用键相信号相邻的峰值去截取启停车过程转频振动信号，得到振动信号的每一个周期点值； 4)对启停车过程转频振动信号作希尔伯特变换，绘制变换后的信号幅值包络波形曲线. 5)根据启停车过程转频振动信号的周期点值，得到信号每个周期内的最大值和最小值，然后对最大值和最小值分别作插值、平均处理，即得到振动峰值曲线的波德图； 6)对启停车过程转频振动信号作希尔伯特变换，绘制变换后的希尔伯特相位波形曲线，根据相位波形曲线寻找周期点处的相位值，作差值处理后，得到真实的傅立叶初相位变化曲线的波德图。
2.根据权利要求I所述的旋转机械启停车过程波德图/奈奎斯特图的快速计算方法，其特征在于，根据步骤5)获得的振动峰值曲线和步骤6)获得的傅立叶相位变化曲线，绘制出奈奎斯特图。
3.根据权利要求I所述的旋转机械启停车过程波德图/奈奎斯特图的快速计算方法，其特征在于，步骤I)中通过电涡流位移传感器采集旋转机械启停车过程的键相信号和振动信号。
4.根据权利要求3所述的旋转机械启停车过程波德图/奈奎斯特图的快速计算方法，其特征在于，所述电涡流位移传感器采样频率为Fs=24000Hz。
5.根据权利要求3所述的旋转机械启停车过程波德图/奈奎斯特图的快速计算方法，其特征在于，测量键相信号的电涡流位移传感器与水平面平行；测量振动信号的电涡流位移传感器与水平面面呈45度夹角。
6.根据权利要求I所述的旋转机械启停车过程波德图/奈奎斯特图的快速计算方法，其特征在于，步骤2)中对步骤I)采集的启停车振动信号进行保相滤波，在每段信号转频左边3Hz处，右边3Hz处各选取一点，对信号作带通滤波，然后合并滤波后的每一段信号，得到启停车过程转频振动信号。
7.旋转机械启停车过程波德图/奈奎斯特图的快速计算方法，其特征在于，包括以下步骤 1)采集旋转机械启停车过程的键相信号和振动信号； 2)用键相信号相邻的峰值去截取启停车过程原始振动信号，得到振动信号的每一个周期点值； 3)对启停车过程原始振动信号作希尔伯特变换，绘制变换后的信号幅值包络波形曲线. 4)根据启停车过程原始振动信号的周期点值，得到信号每个周期内的最大值和最小值，然后对最大值和最小值分别作插值、平均处理，即得到振动峰值曲线的波德图；5)对启停车过程原始振动信号作希尔伯特变换，绘制变换后的希尔伯特相位波形曲线，根据相位波形曲线寻找周期点处的相位值，作差值处理后，得到真实的傅立叶初相位变化曲线的波德图。
8.根据权利要求7所述的旋转机械启停车过程波德图/奈奎斯特图的快速计算方法，其特征在于，根据步骤5)获得的振动峰值曲线和步骤6)获得的傅立叶相位变化曲线，绘制出奈奎斯特图。
全文摘要
本发明公开了一种旋转机械启停车过程波德图/奈奎斯特图的快速计算方法，该方法基于由电涡流传感器测量到的旋转机械启停车过程转子振动信号和键相信号，快速计算得到启停车过程的波德图和奈奎斯特图；整个计算过程基于对以上启停车过程中振动信号的Hilbert变换处理，并结合振动信号的保相滤波处理、键相信号中周期信息的提取，实现了对旋转机械启停车振动信号波德图和奈奎斯特图的快速计算。实验数据分析结果表明该方法计算的波德图和奈奎斯特图曲线比较光滑，而且得到的振动峰值和振动相位更加准确，相位计算精度不受转速高低的影响，是一种比较好的旋转机械启停车非稳态信号波德图和奈奎斯特图计算方法。
文档编号G01H11/06GK102721462SQ20121019656
公开日2012年10月10日 申请日期2012年6月14日 优先权日2012年6月14日
发明者吴婷婷, 张西宁, 徐进杰, 李兵, 郭金良 申请人:西安交通大学