一种旋转叶片振动测量方法与系统的制作方法
【专利摘要】本发明公开了一种旋转叶片振动测量方法与系统，主要包括叶片振动测量方案和基于谱分析的振动参数识别方法，通过单个激光传感器和2?4个光纤传感器分别同步并行采集转速和振动位移，完成振动测量方案；利用最小方差谱分析法对非均匀叶片振动数据进行分析，再通过不同转速位置处的振动峰谱图对比确定叶片振动阶次和频率；最后通过具体叶片转子试验台验证振动测量系统，方案简单易行，方法正确可靠。
【专利说明】
_种旋转叶片振动测量方法与系统
技术领域
[0001] 本发明属于机械测试领域，具体涉及一种旋转叶片振动测量方法与系统。
【背景技术】
[0002] 目前，旋转机械广泛应用于各行各业，如航空发动机、烟气机、叶轮机，它是主要由 涡轮机及其驱动工作机组成的动力机械，其中叶片作为关键部件之一，具有数量多、易损坏 等特点。保证叶片的安全运行能够很大程度降低旋转机械发生故障的概率，因此叶片的安 全受到越来越多的关注。目前测量叶片振动的方法采用非接触测量的叶尖定时方法，叶尖 定时测振技术是将叶尖定时传感器沿径向安装在旋转机械相对静止的壳体上，根据叶片到 达传感器的脉冲时间t与理论到达时间的不同，获得叶片振动的欠采样序列，因此准确得到 叶片振动数据是叶尖定时技术的关键部分，研究更加精确、更加方便的测振方案有重要意 义。而目前对同步和异步叶尖定时参数识别方法中，谱分析的方法发展较慢，然而谱分析在 快速确定振动频率和阶次方面具有优势，因此有必要对谱分析方法进行研究。
【发明内容】
[0003] 本发明所要解决的技术问题是针对上述【背景技术】的不足，提供一种简单易行、准 确可靠的旋转叶片振动测量方法与系统。
[0004] 本发明提供的旋转叶片振动测量方法，包括以下步骤：
[0005] 1)将激光转速传感器安装于旋转机械叶片转子转轴的一侧，用于采集叶盘的转 速;将若干光纤传感器固定在机匣的径向位置，用于采集叶片到来的时间；
[0006] 2)采用同步并行采集卡接收转速与叶片到来时间的信号；
[0007] 3)将采集到的信号传输到控制器中，控制器实现在线数据采集、叶片角度标定、振 动数据位移计算、数据存储和基本分析功能；
[0008] 4)基于最小方差滤波及自相关矩阵的迭代完成对振动峰位移数据的谱分析，对比 不同转速处各振动峰的谱分析结果，确定振动阶次和频率。
[0009] 进一步的，所述光纤传感器设有两至四个。
[0010] 其中，步骤3)具体为：
[0011] 3.1通过占空比得到某一圈叶片的振动位移；
[0012] 3.2在转子低速状态下采集叶片的相对于转速同步信号的初始角度值亦即叶片不 产生振动时到达光纤传感器的位置。
[0013]更进一步的，步骤3.1具体指：
[0014] a.采用划分区间的方法，在两个转速脉冲之间划分若干区间，并保证每个区间中 只有一个叶尖定时脉冲信号;然后利用波形最大值最小值函数对该区间进行最大值点的求 取，同时将该最大值点所对应的X轴的坐标时间求取出来，再求得相对于转速信号的时间 tnbs ;
[0015] b.光纤传感器在第η圈输出的叶片的信号相对转速同步信号转过的角度为anbs = 2 3itnbs/Tn，Tn表示第η圈周期;该叶片的振动角度Aab = anbs-anbsQ，anbsQ为叶片没有振动时该光 纤传感器输出的该叶片的信号相对转速同步信号转过的角度;假设转子的旋转半径为R，该 叶片的振动位移yb= A ab · R，转子旋转η圈后，得到该叶片的振动位移序列{yb}，且该位移 序列对应时间序{tb}。
[0016] 更进一步的，步骤3.2具体指：
[0017] a.测定五个等差转速下叶片相对于转速同步传感器的角度值；
[0018] b.使用线性拟合算法，得出各传感器各叶片角度值随转速变化的线性函数；
[0019] c.通过计算每圈的转速，根据不同叶片不同传感器载入线性方程各自的斜率值和 截距值，运算得出各圈叶片理论角度，即理论条件下叶片无振动时相对于转速同步传感器 的角度值。
[0020]其中，步骤4)具体指：
[0021] 4.1对叶片在不同转速下某两处振动峰数据序列做最小方差谱分析，并加入自相 关矩阵迭代，得到两处振动峰数据序列幅值谱；
[0022] 输入向量:x(ti) = [x(to)，x(ti)，···，x(tk-1)]1"，滤波器权系数向量:h= [ho，hi,… hk-i]T
[0023] 最小方差法权系数的解：
式中…，RG cmxm 为输入向量的k维自相关矩阵，ω为频率点，t为各采样点采样时间，下面加入自相关矩阵的 迭代。
[0024] 自相关矩阵
％幅值谱
[0025]式中i为迭代次数，W( 1} = I，W(i+1 > = diag ( | S(1 |2); I为单位矩阵;N为频点个数。
[0026] 4.2从两幅谱图中找出同一频率位置处的谱峰，该谱峰为当前叶片固有频率，即共 振频率；
[0027] 4.3通过振动峰处的转频及共振频率求得振动阶次。
[0028] 步骤4.2中如果存在多个相同谱峰，通过坎贝尔图加以分别。
[0029] -种旋转叶片振动测量系统，该系统包括传感检测装置、旋转叶片转子、控制器以 及连接线路;其特征在于，所述传感检测装置包括激光转速传感器以及若干光纤传感器，所 述激光转速传感器设置在叶片转子的转轴的一侧，所述光纤传感器固定在机匣的径向位 置;各传感检测装置与控制器相连。
[0030] 本发明提供的旋转叶片振动测量方法与系统针对振动数据非均匀、无规则的特 点，使用了最小方差谱分析方法，相较于传统谱分析，充分考虑非均匀采样的影响。特别适 用于转速变化快，多个叶尖定时传感器数据同时分析的情况。
[0031] 本发明可以使得叶尖定时传感器的安装便捷，不受结构的限制，便于现场测试应 用。
【附图说明】
[0032] 图1是本发明旋转叶片振动测量系统的结构示意图；
[0033] 图2是叶片振动位移获取方式的流程图；
[0034]图3是叶片角度校准方式的流程图；
[0035]图4是叶片振动位移折线图；
[0036]图5是实施例中某一振动峰处振动数据选取及分析结果图；
[0037] 图6是实施例中某一振动峰处振动数据选取及分析结果图；
[0038] 其中，1-激光转速传感器、2-光纤传感器。
【具体实施方式】
[0039] 本发明提供一种旋转叶片振动测量方法与系统，为使本发明的目的，技术方案及 效果更加清楚，明确，以及参照附图并举实例对本发明进一步详细说明。应当理解，此处所 描述的具体实施仅用以解释本发明，并不用于限定本发明。
[0040] 本文重点解决的问题是制定准确快捷叶片测振方案和谱分析参数识别方法，为 此，我们采用一个安装于转轴的激光转速传感器1、二到四个光纤传感器2、PXI并行采集卡 和高性能采集分析软件实现叶片振动采集。方案实施的具体内容如下：
[0041] (1)通过安装于转轴一侧的单个激光转速传感器确定叶盘转速，2至4个光纤传感 器固定于机匣径向位置采集叶片到来时间。
[0042] (2)PXI同步并行采集卡同时接受转速脉冲和叶尖定时脉冲。该采集卡和传感器的 相互配合，具有传感器安装数量少且定时精度高，转速时间和叶片到来时间误差小。
[0043] (3)采集软件完成了在线数据采集、叶片角度标定、振动数据位移计算、数据存储 和基本分析功能。
[0044] (4)离线谱分析软件基于最小方差滤波及自相关矩阵的迭代完成对位移数据的谱 分析，通过不同转速的振动峰谱图对比确定。
[0045] 对上述方案的各部分的关键内容，在下面进行详细描述。
[0046] 1、数据采集
[0047]在线数据采集模块硬件设备基于NI PXIe-4499高速数据采集卡，由于叶尖定时测 振法需要准确测得叶片到来时间，因而只有高采样率才能保证叶片到达时间准确性，该数 据采集卡最高采样率可以达到204.8ks/s。同时单端采集卡提供的16通道模拟信号并行输 入能够使得各个通道采集的信号同时到达PXI系统，使得各传感器信号之间不存在延迟。由 于激光转速传感器和光纤传感器所测得的信号存在严格的对应关系，非并行输入会造成信 号之间延迟误差，这样就需要其他传感器辅助，增加了系统复杂程度且降低采集精度。采用 单个激光转速传感器搭配PXI并行采集卡方法满足了信号获取既方便又准确的要求。采集 卡采用的ADC模数转换分辨率为24位。
[0048]基于labview的采集软件主要在以下模块有以下特点：
[0049] (1)振动位移计算模块
[0050] a.采用划分区间的方法，即使得整圈数内该区间内时域信号只存在单个叶片在单 个传感器上响应的模拟信号。然后利用LabVIEW波形函数工具包中的波形最大值最小值函 数对该区间进行最大值点的求取，同时将该最大值点所对应的X轴的坐标时间求取出来，再 求得相对于转速信号的时间t nbs。
[0051 ] b.传感器s在第η圈输出的叶片b的信号相对转速同步信号转过的角度为anbs = 23T tnbs/Tn，Tn表示第η圈周期。因此叶片b的振动角度Aab = anbs-anbsQ，anbsQ为叶片没有振动时传 感器S输出的叶片b的信号相对转速同步信号转过的角度。假设转子的旋转半径为R，叶片b 的振动位移yb= Aab · R，转子旋转η圈后，得到叶片b的振动位移序列{yb}，且该位移序列对 应时间序{tb}。叶片振动位移获取模块流程图如图2所示。
[0052] (2)叶片角度标定模块
[0053] 该模块实现的功能是叶尖定时计算算法中叶片角度校准功能，即在转子低速状态 下系统采集叶片的相对于转速同步信号的初始角度值亦即叶片不产生振动时到达光纤传 感器的位置。
[0054] a.测定200印111、400印111、600印111、800印111、1000印1]1五个转速下叶片相对于转速同步 传感器的角度值。
[0055] b.使用内置线性拟合模块，得出各传感器各叶片角度值随转速变化的线性函数。
[0056] c.系统通过计算每圈的转速，根据不同叶片不同传感器载入线性方程各自的斜率 值和截距值，通过运算得出各圈叶片理论角度，即理论条件下叶片无振动的相对于转速同 步传感器的角度值。该方法流程图如图3所示。
[0057] 2、参数识别
[0058] 对叶尖定时数据参数辨识方法的研究是本文另一个重点，针对数据不规则、非均 匀的特点，我们希望采用谱分析方法进行参数识别。目前针对叶片振动数据参数辨识谱分 析方法较少，而传统谱分析忽视非均匀采样的问题，我们将最小方差非均匀谱分析应用于 叶片振动数据，并编制相应分析软件。该方法主要包括以下步骤：
[0059] a.对不同转速下某两处振动峰数据序列做最小方差谱分析，并加入自相关矩阵迭 代，得到各自幅值谱。
[0060]
[0061]
[0062] b.从两幅谱图中找出同一频率位置处的谱峰，该谱峰为当前叶片固有频率。如果 存在多个相同谱峰，需要借助坎贝尔图加以分别。
[0063] c.通过振动峰处的转频及共振频率求得振动阶次。
[0064] 3、具体实施
[0065] 3.1操作流程
[0066] (1)检查叶片转子实验台，安装激光转速传感器和4支光纤传感器，确保叶片转动 过程中不会和传感器探头发生干扰，叶片连接处螺栓紧固。
[0067] (2)打开在线采集分析软件，采样率5000Hz，采样点数50000，传感器灵敏度lmv/ mm，并设置存储路径。并保证采集卡正常工作。
[0068] (3)接通传感器电源，启动试验台，完成200印111、400印111、600印111、800印111、1000印111 叶片角度标定工作。
[0069] (4)完成叶片转子Orpm至1500rpm升降速扫频过程，得到振动位移数据。如图4为试 验现场采集数据。
[0070] (5)试验结束后，切断电源。
[0071] (6)使用离线谱分析软件对振动位移数据进行分析，得出结果。
[0072] 3.2谱分析方法
[0073] 测振方案如上节所示，下面对振动数据进行谱分析的过程。4个光纤传感器相对于 叶片1的安装位置分别为7°、22°、37°、52°，整个测试过程持续100s，转速由0Hz升至20Hz再 降为零。图5为724圈处振动数据谱分析结果，图6为759圈处振动数据谱分析结果，两处振动 峰的转频分别为9.7Hz和7.1Hz，分别对应两谱图中第一个谱峰，也说明谱图分析结果存在 转频成分。
[0074] 下一步寻找共振频率，图5中各谱峰为9.6Hz、20.0Hz、30.12Hz，图6中各谱峰为 7.2Hz、14.5Hz、21.7Hz、28.5Hz，其中除去9.6Hz、7.2Hz两处转频成分的谱峰，其余频率成分 中相近的为20.0Hz和21.7Hz、30.12Hz和28.5Hz ;此时我们需要借助一些先验知识，通过坎 贝尔图我们可以预先知道〇. 5mm叶片在转频为9.6Hz至7.2Hz范围内，动频处于30Hz左右，因 此我们认为30.12Hz和28.5Hz分别为转频9.6Hz、7.2Hz处共振频率。需要指出的是，在试验 中由于叶片动频随转速产生轻微改变，不同转速下叶片的共振频率并不完全相等。
[0075]在已知转频和共振频率的情况下，两处振动峰的振动阶次求得为3.1和4.96，表1 列出了最小二乘谱分析法和其他方法对振动峰的分析结果，并与坎贝尔图比较。从中可以 看出，谱分析方法阶次辨识结果与坎贝尔图相比误差在± 5 %以内；频率分析结果经其他方 法检验，结果正确。
[0076] 表1 724圈和757圈两出振动峰各方法分析结果
[0077]
[0078]上述工作原理的说明仅为本发明的一个特例。因此凡是按照本发明提出的技术思 想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。
【主权项】
1. 一种旋转叶片振动测量方法，其特征在于，该方法包括W下步骤： 1) 将激光转速传感器安装于旋转机械叶片转子转轴的一侧，用于采集叶盘的转速;将 若干光纤传感器固定在机厘的径向位置，用于采集叶片到来的时间； 2) 采用同步并行采集卡接收转速与叶片到来时间的信号； 3) 将采集到的信号传输到控制器中，控制器实现在线数据采集、叶片角度标定、振动数 据位移计算、数据存储和基本分析功能； 4) 基于最小方差滤波及自相关矩阵的迭代完成对振动峰位移数据的谱分析，对比不同 转速处各振动峰的谱分析结果，确定振动阶次和频率。2. 根据权利要求1所述的一种旋转叶片振动测量方法，其特征在于，所述光纤传感器设 有两至四个。3. 根据权利要求1所述的一种旋转叶片振动测量方法，其特征在于，步骤3)具体为： 3.1通过占空比方法得到每圈叶片的振动位移； 3.2在转子低速状态下采集叶片的相对于转速同步信号的初始角度值亦即叶片不产生 振动时相对于转速同步信号的初始角度值。4. 根据权利要求3所述的一种旋转叶片振动测量方法，其特征在于，步骤3.1具体指： a. 采用划分区间的方法，在两个转速脉冲之间划分若干区间，并保证每个区间中只有 一个叶尖定时脉冲信号；然后利用波形最大值最小值函数对该区间进行最大值点的求取， 同时将该最大值点所对应的X轴的坐标时间取出，再求得相对于转速信号的时间tnbs; b. 光纤传感器在第η圈输出的叶片的信号相对转速同步信号转过的角度为anbs = 23i tnbs/Tn，Tn表示第η圈周期；该叶片的振动角度Aab = anbs-anbsO，anbs日为叶片没有振动时该光 纤传感器输出的该叶片的信号相对转速同步信号转过的角度;假设转子的旋转半径为R，该 叶片的振动位移yb= Aab . R，转子旋转η圈后，得到该叶片的振动位移序列{yb}，且该位移 序列对应时间序{tb}。5. 根据权利要求3所述的一种旋转叶片振动测量方法，其特征在于，步骤3.2具体指： a. 测定五个等差转速下叶片相对于转速同步传感器的角度值； b. 使用线性拟合算法，得出各传感器各叶片角度值随转速变化的线性函数； C.通过计算每圈的转速，根据不同叶片不同传感器载入线性方程各自的斜率值和截距 值，运算得出各圈叶片理论角度，即理论条件下叶片无振动的相对于转速同步传感器的角 度值。6. 根据权利要求1所述的一种旋转叶片振动测量方法，其特征在于，步骤4)具体指： 4.1对叶片不同转速位置处的两处振动峰数据序列做最小方差谱分析，并加入自相关 矩阵迭代，得到两处振动峰数据序列幅值谱； 输入向量:x(ti) = [x(to)，x(ti)，…，x(tk-i)]T;滤波器权系数向量:h=比o，hi，一hk-i]T 最小方差法权系数的解：氏中史。=[1，6-'心^-'心，RecMXM为输入 向量的k维自相关矩阵，ω为频率点，t为各采样点采样时间，下面加入自相关矩阵的迭代。 自相关矩阵：式中i为迭代次数，W(i> = I，W心i) = diag( I SW 12); I为单位矩阵;N为频点个数。 4.2从两幅谱图中找出同一频率位置处的谱峰，该谱峰为当前叶片固有频率，即共振频 率. 4.3通过振动峰处的转频及共振频率求得振动阶次。7. 根据权利要求6所述的一种旋转叶片振动测量方法，其特征在于，步骤4.2中如果存 在多个相同谱峰，通过坎贝尔图加 W分别。8. -种旋转叶片振动测量系统，该系统包括传感检测装置、旋转叶片转子、控制器W及 连接线路;其特征在于，所述传感检测装置包括激光转速传感器W及若干光纤传感器，所述 激光转速传感器设置在叶片转子的转轴的一侧，所述光纤传感器固定在机厘的径向位置； 各传感检测装置与控制器相连。
【文档编号】G01H9/00GK105973448SQ201610073172
【公开日】2016年9月28日
【申请日】2016年2月2日
【发明人】岳林, 杨彬, 梁安阳, 时辰, 刘豪
【申请人】南京航空航天大学