CN111256636B - 一种叶片扭转量的测量方法 - Google Patents
一种叶片扭转量的测量方法 Download PDFInfo
- Publication number
- CN111256636B CN111256636B CN201811450604.0A CN201811450604A CN111256636B CN 111256636 B CN111256636 B CN 111256636B CN 201811450604 A CN201811450604 A CN 201811450604A CN 111256636 B CN111256636 B CN 111256636B
- Authority
- CN
- China
- Prior art keywords
- blade
- value
- point
- measured
- displacement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明涉及一种本叶片扭转量的测量方法,包括以下步骤:S1、获得叶片位移测量值:利用传感器同时对被测级叶片的侧面和转子进行测量,获得被测级叶片的叶片位移测量值Y;S2、确定极值:获取步骤S1得到的叶片位移测量值Y中的所有极大值、得到一组极大值数列,以及叶片位移测量值Y中的所有极小值、得到一组极小值数列;S3、确定有效极值:排除步骤S2得到的极大值数列中的异常极大值、得到一组有效极大值数列{YAm},以及排除极小值数列{YBm}中的异常极小值、得到一组有效极小值数列{YBm};S4、确定叶片扭转量:根据有效极大值数列{YAm}中的YA(j+1)、以及有效极小值数列{YBm}中的YBj,确定所测量的第j+1个被测级叶片的叶片扭转量S(j+1)。
Description
技术领域
本发明涉及汽轮机领域,具体涉及一种叶片扭转量的测量方法。
背景技术
传统的测量汽轮机叶片扭转量的方法是使用专用的夹具夹持百分表,手工沿着转子上的基准面移动百分表及其夹具,逐个叶片测量、记录和计算,汽轮机转子上的叶片多达数千片,这样手工测量和计算的效率非常低。此外,在两级叶片轴向距离较小的情况下,普通的百分表体积过大,无法放入进行测量。
发明内容
鉴于以上所述现有技术的缺点,本发明要解决的技术问题在于提供一种叶片扭转量的测量方法,能够提高叶片扭转量的测量效率和精度。
为实现上述目的,本发明提供一种叶片扭转量的测量方法,包括以下步骤:S1、获得叶片位移测量值:使转子转动,所述转子上安装有至少一圈被测级叶片,利用传感器同时对被测级叶片的侧面和转子进行测量,获得被测级叶片的叶片实际位移量U和转子的转子轴向窜动量V;根据同时刻测得的叶片实际位移量U和转子轴向窜动量V获得被测级叶片的叶片位移测量值Y,所述叶片实际位移量U、转子轴向窜动量V和叶片位移测量值Y都是一系列值;
S2、确定极值:获取步骤S1得到的叶片位移测量值Y中的所有极大值、得到一组极大值数列,获取步骤S1得到的叶片位移测量值Y中的所有极小值、得到一组极小值数列;
S3、确定有效极值:排除步骤S2得到的极大值数列中的异常极大值、得到一组有效极大值数列{YAm},排除步骤S2得到的极小值数列中的异常极小值、得到一组有效极小值数列{YBm},其中有效极大值数列{YAm}中的各项和有效极小值数列{YBm}中的各项都按照测量的先后顺序排列,m为正整数;S4、确定叶片扭转量:取所述有效极大值数列{YAm}中的第j+1项YA(j+1)、以及所述有效极小值数列{YBm}中的第j项YBj,根据YBj和YA(j+1)确定所测量的第j+1个被测级叶片的叶片扭转量S(j+1),j为自然数。
进一步地,所述步骤S1中,所述传感器为位移传感器,所述传感器包括垂直于被测级叶片侧面的第一传感器、以及垂直于转子端面的第二传感器,所述第一传感器用于测量获取叶片实际位移量U,所述第二传感器用于测量获取转子轴向窜动量V;所测得的叶片实际位移量U和转子轴向窜动量V为离散的数值点;根据叶片实际位移量U和转子轴向窜动量V得到离散的叶片位移测量值、构成叶片位移测量值数列{Yn},叶片位移测量值数列{Yn}中的各项按照其测量的序号排列。
优选地,所述步骤S1中叶片位移测量值Y通过下述方式确定:(a)所述第一传感器和第二传感器安装在被测级叶片的同侧时,叶片位移测量值Y等于同时刻的叶片实际位移量U和转子轴向窜动量V的差值;(b)所述第一传感器和第二传感器分别安装在被测级叶片的两侧时,叶片位移测量值Y等于同时刻的叶片实际位移量U和转子轴向窜动量V的和值。
优选地,所述步骤S2中,所述叶片位移测量值数列{Yn}中极大值和极小值的判断标准为:所述极大值大于其前一个值、且大于或等于其后一个值,所述极小值小于其后一个值、且小于或等于其前一个值。
进一步地,所述步骤S2中得到极大值和极小值数量相同,得到所述极大值数列中的各项和极小值数列中的各项都按照测量的先后顺序排列。
进一步地,所述步骤S3中,所述异常极大值和异常极小值的判断方法为:算出极大值和极小值的极值平均值,将极大值数列和极小值数列中相同序号的极大值和极小值进行求差、得到极值差值,当该极值差值小于极值平均值的30%时,则该序号的极大值为异常极大值、该序号的极小值为异常极小值。
进一步地,所述步骤S1中,所述第一传感器测量的被测级叶片侧面为平面;根据所述步骤S1中所获得的叶片位移测量值数列{Yn}得到叶片位移测量值的波形图,所述波形图的横坐标为转子转动距离X,所述波形图的纵坐标为叶片位移测量值Y,所述有效极大值数列{YAm}中各项在波形图中对应的点为Am,所述有效极小值数列{YBm}中各项在波形图中对应的点为Bm;所述异常极大值在波形图中对应的点为Ci,所述异常极小值在波形图中对应的点为Di;所述步骤S4中,第j+1个被测级叶片的叶片扭转量S(j+1)的计算方式为:
其中,XAj和XA(j+1)分别为点Aj和A(j+1)对应的横坐标值;XBj为点Bj对应的横坐标值;e为异常点补偿量。
进一步地,所述异常点补偿量e的确定方式为:(a)当点A(j+1)和点Bj之间无点Ci和点Di时,e=0;(b)当点A(j+1)和点Bj之间有点Ci和点Di、且点Ci和点Di的纵坐标值Y不同时,根据点Ci和点Di在叶片位移测量值数列{Yn}中对应的序号确定补偿点Ei,e为补偿点Ei的横坐标值和点Ci横坐标值的差值;(c)当点A(j+1)和点Bj之间有点Ci和点Di、且点Ci和点Di的纵坐标值Y相同时,e为点Ci的横坐标值和点Di的横坐标值的差值。
进一步地,所述步骤S1中,所述一圈被测级叶片中末叶片的围带顶部安装有感应块,所述传感器还包括第三传感器,当感应块随末叶片转动到第三传感器处时,所述第三传感器能够感应到感应块,并开始记录第一传感器和第二传感器的测量数据。
进一步地,所述被测级叶片具有依次相连的叶根、叶身和围带;将所述第一传感器设置在被测级叶片的叶根侧面,通过所述步骤S1至S4获得被测级叶片的叶根扭转量;将所述第一传感器设置在被测级叶片的围带侧面,通过所述步骤S1至S4获得被测级叶片的围带扭转量;将被测级叶片的叶根扭转量与围带扭转量相减,得到被测级叶片的相对扭转量。
如上所述,本发明涉及的叶片扭转量的测量方法,具有以下有益效果:
利用传感器自动连续地对一圈被测级叶片的侧面进行测量,通过测量得到的叶片实际位移量U和转子的转子轴向窜动量V,经过计算即得到一圈被测级叶片中多个被测级叶片的叶片扭转量,同时可以有效地排除转子短时间停顿或反向转动对叶片扭转量测量的影响。采用本发明的叶片扭转量的测量方法,操作方便,极大地提高了叶片扭转量的测量效率,并有效地提高了测量的准确性。
附图说明
图1为本发明的安装和工作示意图。
图2为图1中K方向视图,该图省略了转子和第三传感器。
图3为本发明中的叶片位移测量值的波形图。
图4为图3中a圈的放大图。
图5为图3中b圈的放大图。
元件标号说明
1 转子
2 被测级叶片
21 叶根
22 叶身
23 围带
3 第一传感器
4 第二传感器
5 第三传感器
具体实施方式
以下由特定的具体实施例说明本发明的实施方式,熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效。
须知,本说明书附图所绘的结构、比例、大小等,均仅用以配合说明书所揭示的内容,以供熟悉此技术的人士了解与阅读,并非用以限定本发明可实施的限定条件,故不具技术上的实质意义,任何结构的修饰、比例关系的改变或大小的调整,在不影响本发明所能产生的功效及所能达成的目的下,均应仍落在本发明所揭示的技术内容能涵盖的范围内。同时,本说明书中所引用的如“上”、“下”、“左”、“右”、“中间”等的用语,亦仅为便于叙述明了,而非用以限定本发明可实施的范围,其相对关系的改变或调整,在无实质变更技术内容下,当亦视为本发明可实施的范畴。
如图1至图5所示,本发明提供一种叶片扭转量的测量方法,包括以下步骤:
S1、获得叶片位移测量值:安装好转子1、叶片以及传感器,使转子1转动,转子1上安装有至少一圈被测级叶片2,利用传感器同时对被测级叶片2的侧面和转子1进行测量,获得被测级叶片2的叶片实际位移量U和转子1的转子轴向窜动量V;根据同时刻测得的叶片实际位移量U和转子轴向窜动量V获得被测级叶片2的叶片位移测量值Y,叶片实际位移量U、转子轴向窜动量V和叶片位移测量值Y都是一系列值。利用传感器测得的叶片实际位移量U和转子1的转子轴向窜动量V可以为连续值,也可以是多个离散的数值。
作为优选设计,如图1和图2所示,在本实施例中,传感器为位移传感器,传感器测量数据的频率约为25Hz,即每个传感器每秒测得约25个数据,精度要高于1μm。传感器包括垂直于被测级叶片2侧面的第一传感器3、以及垂直于转子1端面的第二传感器4,第一传感器3用于测量获取叶片实际位移量U,第二传感器4用于测量获取转子轴向窜动量V。在一圈被测级叶片2中的末叶片的围带顶部安装有感应块,传感器还包括第三传感器5,当感应块随末叶片转动到第三传感器5处时,第三传感器5能够感应到感应块,并开始记录第一传感器3和第二传感器4的测量数据,此时第一传感器3对应的被测级叶片2为第1个被测级叶片2。所测得的叶片实际位移量U和转子轴向窜动量V为离散的数值点,根据叶片实际位移量U和转子轴向窜动量V得到离散的叶片位移测量值Y、构成叶片位移测量值数列{Yn},其中,叶片位移测量值数列{Yn}中的各项按照其测量的序号排列,n为自然数,n具体大小根据实际的测量情况确定,测量过程中保证第一传感器3至少完整地扫过一个被测级叶片2。优选地,可将叶片位移测量值数列{Yn}用列表形式表示,得到下述表1:
表1叶片位移测量值
测量序号 | 0 | 1 | 2 | 3 | …… | n |
叶片位移测量值 | Y<sub>0</sub> | Y<sub>1</sub> | Y<sub>2</sub> | Y<sub>3</sub> | …… | Y<sub>n</sub> |
叶片位移测量值Y的计算根据第一传感器3、以及第二传感器4的位置确定,具体地计算方式如下:(a)当第一传感器3和第二传感器4安装在被测级叶片2的同侧时,如图1所示,叶片位移测量值Y等于同时刻的叶片实际位移量U和转子轴向窜动量V的差值;(b)当第一传感器3和第二传感器4分别安装在被测级叶片2的两侧时,叶片位移测量值Y等于同时刻的叶片实际位移量U和转子轴向窜动量V的和值。
S2、确定极值:获取步骤S1得到的叶片位移测量值Y中的所有极大值、得到一组极大值数列,获取步骤S1得到的叶片位移测量值Y中的所有极小值、得到一组极小值数列。在本实施例中,叶片位移测量值Y为离散型的,构成叶片位移测量值数列{Yn},参考上述表1,当叶片位移测量值数列{Yn}中的某一项数值大于其前一个值、且大于或等于其后一个值,则该数值为极大值,当叶片位移测量值数列{Yn}中的某一项数值小于其后一个值、且小于或等于其前一个值,则该数值为极小值。本实施例中,叶片位移测量值数列{Yn}中的极大值和极小值数量相同,都为k+1个,k为极大值数列中各项和极小值数列中各项的序号,k为自然数。极大值数列中的各项和极小值数列中的各项都按照其原测量序号的先后进行排列。优选地,可将极大值数列和极小值数列分别用列表形式表示,如下表2和表3所示:
表2极大值数列
序号 | 0 | 1 | 2 | 3 | 4 | …… | k |
原测量序号 | 10 | 55 | 100 | 150 | 194 | …… | |
极大值 | Y<sub>10</sub> | Y<sub>55</sub> | Y<sub>100</sub> | Y<sub>150</sub> | Y<sub>194</sub> | …… |
其中,表2的第一行为极大值数列中的各项极大值的序号,第二行为各项极大值在叶片位移测量值数列{Yn}中的原序号,第三行为各项极大值的数值。
表3极小值数列
序号 | 0 | 1 | 2 | 3 | 4 | …… | k |
原测量序号 | 14 | 60 | 104 | 155 | 199 | …… | |
极小值 | Y<sub>14</sub> | Y<sub>60</sub> | Y<sub>104</sub> | Y<sub>155</sub> | Y<sub>199</sub> | …… |
其中,表3的第一行为极小值数列中的各项极小值的序号,第二行为各项极小值在叶片位移测量值数列{Yn}中的原序号,第三行为各项极小值的数值。
S3、确定有效极值:排除步骤S2得到的极大值数列中的异常极大值、得到一组有效极大值数列{YAm},排除步骤S2得到的极小值数列中的异常极小值、得到一组有效极小值数列{YBm},其中有效极大值数列{YAm}中的各项和有效极小值数列{YBm}中的各项都按照其原测量序号的先后顺序排列,m为正整数。由于转子1可能会出现短时间停顿或反向转动的异常,这些异常转动会导致在叶片位移测量值数列{Yn}中产生一些异常极大值和异常极小值,为了防止这些转子1的异常转动对测量数据产生影响,需要找出并排除这些异常极大值和异常极小值。
在本实施例中,异常极大值和异常极小值的判断方法为:算出极大值和极小值的极值平均值,将极大值数列和极小值数列中相同序号的极大值和极小值进行求差、得到极值差值,当该极值差值小于极值平均值的30%时,则该序号的极大值为异常极大值,该序号的极小值为异常极小值。具体地,以表2和表3中的数据为例,根据表2和表3,得到极值差值列表,如下表4所示:
表4极值差值
序号 | 0 | 1 | 2 | 3 | 4 | …… | k |
极大值 | Y<sub>10</sub> | Y<sub>55</sub> | Y<sub>100</sub> | Y<sub>150</sub> | Y<sub>194</sub> | …… | |
极小值 | Y<sub>14</sub> | Y<sub>60</sub> | Y<sub>104</sub> | Y<sub>155</sub> | Y<sub>199</sub> | …… | |
极值差值 | Y<sub>10</sub>-Y<sub>14</sub> | Y<sub>55</sub>-Y<sub>60</sub> | Y<sub>100</sub>-Y<sub>104</sub> | Y<sub>150</sub>-Y<sub>155</sub> | Y<sub>194</sub>-Y<sub>199</sub> | …… |
其中,表4的第一行为极小值数列和极大值数列中各项的序号,第二行为各项极大值的数值,第三行为各项极小值的数值,第四行为极大值减去极小值的极值差值。
将表4中第二行的所有极大值和第三行中的所有极小值相加,再除以极大值和极小值的总数2(k+1),得到极值平均值。将表4中的每一个极值差值与极值平均值相比,当该某个极值差值小于极值平均值的30%时,则该列中的极大值和极小值为异常极大值和异常极小值,例如,若第二列中的极值差值Y55-Y60小于极值平均值的30%时,则Y55为异常极大值、Y60为异常极小值。根据上述方法找出表4中所有的异常极大值和异常极小值,并将其在表4中排除,则表4剩下的即为有效极大值和有效极小值。具体地,以表4中的数据为例,当表4中第二列和第四列的极值差值小于极值平均值的30%时,得到的有效极大值数列{YAm}和有效极小值数列{YBm}分别以表格形式表示,如下表5和表6所示:
表5有效极大值数列
序号 | 1 | 2 | 3 | …… |
原测量序号 | 10 | 100 | 150 | |
有效极大值 | Y<sub>10</sub> | Y<sub>100</sub> | Y<sub>150</sub> | …… |
表6有效极小值数列
序号 | 1 | 2 | 3 | …… |
原测量序号 | 14 | 60 | 155 | |
有效极小值 | Y<sub>14</sub> | Y<sub>60</sub> | Y<sub>155</sub> | …… |
其中,表5和表6中第二行中的原测量序号,即为该列中的有效极大值数和有效极小值的各项在原叶片位移测量值数列{Yn}中的序号。
S4、确定叶片扭转量:取有效极大值数列{YAm}中的第j+1项YA(j+1)、以及有效极小值数列{YBm}中的第j项YBj,根据YBj和YA(j+1)确定所测量的第j+1个被测级叶片2的叶片扭转量S(j+1),其中j为自然数。
在本实施例中,如图2所示,第一传感器3测量的被测级叶片2侧面为平面,根据获得的叶片位移测量值数列{Yn}得到叶片位移测量值的波形图,如图3所示,波形图的横坐标为转子1转动距离X,波形图的纵坐标为叶片位移测量值Y,有效极大值数列{YAm}中各项在波形图中对应的点为Am,有效极小值数列{YBm}中各项在波形图中对应的点为Bm;异常极大值在波形图中对应的点为Ci,异常极小值在波形图中对应的点为Di,i为正整数。采用本实施例中的位移传感器进行测量时,参考图2和图3所示,图2中箭头方向为转子1转动的方向,第1个被测级叶片2上获得的叶片位移测量值Y对应于点A1之前的波段,在第j+1个被测级叶片2上获得的叶片位移测量值Y,对应于图3中点Bj到点A(j+1)段,j为正整数;测量过程中当第一传感器3从第j个被测级叶片2进入第j+1个被测级叶片2,理论上该瞬间获得的叶片位移测量值Y应对应于图3中Aj点到Bj’点,但由于转子1的旋转以及传感器回弹需要时间,实际只测到Bj点处的数值,而无法测到Bj’点的数值,因此产生了一个测量误差。在本实施例中,因为被第一传感器3测量的被测级叶片2侧面为平面,当转子1不出现反向转动或短暂停止时,在一个被测级叶片2上得到的叶片位移测量值Y在波形图中呈线性分布,如图3中的点B2到点A3段。根据叶片位移测量值的波形图,采用相似三角形法进行补偿计算,第j+1个被测级叶片2的叶片扭转量S(j+1)的计算方式为:
其中,XAj和XA(j+1)分别为点Aj和A(j+1)对应的横坐标值;XBj为点Bj对应的横坐标值;e为异常点补偿量。
测量过程中,当转子出现反向转动或短暂停止的异常时,叶片位移测量值Y中会产生异常极大值和异常极小值,参考图3中的点C1和点D1、以及点C2和点D2所示。在本实施例中,异常点补偿量e的确定方式为:
(a)当点A(j+1)和点Bj之间无点Ci和点Di时,e=0;
(b)当点A(j+1)和点Bj之间有点Ci和点Di、且点Ci和点Di的纵坐标值Y不同时,根据点Ci和点Di在叶片位移测量值数列{Yn}中对应的序号确定补偿点Ei,e为补偿点Ei的横坐标值和点Ci横坐标值的差值;参考图3和图4所示,点A2和点B1之间有点C1和点D1,假设点C1在叶片位移测量值数列{Yn}中对应的原测量序号为c、D1在叶片位移测量值数列{Yn}中对应的原测量序号为d,则确定补偿点E1在叶片位移测量值数列{Yn}中对应的原测量序号为2d-c。再通过点D1和补偿点E1在叶片位移测量值数列{Yn}的对应的原测量序号,以此在叶片位移测量值的波形图中来确定点C1横坐标值和补偿点E1的横坐标值,然后两者横坐标的取差值,得到异常点补偿量e。
(c)当点A(j+1)和点Bj之间有点Ci和点Di、且点Ci和点Di的纵坐标值Y相同时,e为点Ci的横坐标值和点Di的横坐标值的差值。参考图3和图5所示,点A4和点B3之间有点C2和点D2,通过点C2和点D2在叶片位移测量值数列{Yn}的对应的原测量序号,以此在叶片位移测量值的波形图中来确定点C2横坐标值和点D2的横坐标值,然后两者横坐标的取差值,得到异常点补偿量e。
进一步地,本发明涉及的测量方法可以用于测量叶片侧面的相对扭转量,被测级叶片2具有依次相连的叶根21、叶身22和围带23;将第一传感器3设置在被测级叶片2的叶根21侧面,通过上述步骤S1至S4获得一组被测级叶片2的叶根扭转量;将第一传感器3设置在被测级叶片2的围带23侧面,通过上述步骤S1至S4获得一组被测级叶片2的围带扭转量。当被测级叶片2的叶根扭转量和围带扭转量两组数据的测量起点相同时,则这两组数据一一对应,相对应的叶根扭转量和围带扭转量属于同一个被测级叶片2,将两组数据一一相减即可以得到每个被测级叶片2的相对扭转量。
由于本发明中测量数据较多,作为优化设计,本发明可设计专用的计算程序进行计算,优选地,采用labview软件进行数据采集和程序设计。
采用本发明涉及的叶片扭转量测量方法,利用传感器自动连续的对一圈被测级叶片2的侧面进行测量,通过测量得到的叶片实际位移量U和转子1的转子轴向窜动量V,经过计算即得到一圈被测级叶片2中多个被测级叶片2的叶片扭转量,操作方便,极大地提高了叶片扭转量的测量效率,并有效地提高了测量的准确性。特别地,采用本实施例的叶片扭转量测量方法,用相似三角形法进行补偿计算,消除了传感器回弹时间和转子1的转动对测量结果的影响,并同时在三角形法补偿计算过程中,有效地排除转子1的反向转动或短暂停止等异常转动对叶片扭转量测量的影响,测量结果更加精确。同时可采用labview软件进行数据采集和程序设计,实现计算机辅助测量和计算,进一步地提高叶片扭转量的测量效率。
综上所述,发明有效克服了现有技术中的种种缺点而具有高度产业利用价值。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (10)
1.一种叶片扭转量的测量方法,其特征在于:包括以下步骤:
S1、获得叶片位移测量值:使转子转动,所述转子上安装有至少一圈被测级叶片,利用传感器同时对被测级叶片的侧面和转子进行测量,获得被测级叶片的叶片实际位移量U和转子的转子轴向窜动量V;根据同时刻测得的叶片实际位移量U和转子轴向窜动量V获得被测级叶片的叶片位移测量值Y:当叶片实际位移量U和转子轴向窜动量V的方向相同时,叶片位移测量值Y等于同时刻的叶片实际位移量U和转子轴向窜动量V的差值;当叶片实际位移量U和转子轴向窜动量V的方向相反时,叶片位移测量值Y等于同时刻的叶片实际位移量U和转子轴向窜动量V的和值;所述叶片实际位移量U、转子轴向窜动量V和叶片位移测量值Y都是一系列值;
S2、确定极值:获取步骤S1得到的叶片位移测量值Y中的所有极大值、得到一组极大值数列,获取步骤S1得到的叶片位移测量值Y中的所有极小值、得到一组极小值数列;
S3、确定有效极值:排除步骤S2得到的极大值数列中的异常极大值、得到一组有效极大值数列{YAm},排除步骤S2得到的极小值数列中的异常极小值、得到一组有效极小值数列{YBm},其中有效极大值数列{YAm}中的各项和有效极小值数列{YBm}中的各项都按照测量的先后顺序排列,m为正整数;
S4、确定叶片扭转量:取所述有效极大值数列{YAm}中的第j+1项YA(j+1)、以及所述有效极小值数列{YBm}中的第j项YBj,根据YBj和YA(j+1)确定所测量的第j+1个被测级叶片的叶片扭转量S(j+1),j为自然数,叶片扭转量S(j+1)的确认方式为:根据所述步骤S1中所获得的叶片位移测量值数列{Yn}得到叶片位移测量值的波形图,所述波形图的横坐标为转子转动距离X,所述波形图的纵坐标为叶片位移测量值Y,所述有效极大值数列
{YAm}中各项在波形图中对应的点为Am,所述有效极小值数列{YBm}中各项在波形图中对应的点为Bm;所述异常极大值在波形图中对应的点为Ci,所述异常极小值在波形图中对应的点为Di;所述第j+1个被测级叶片的叶片扭转量S(j+1)的计算方式为:
其中,XAj和XA(j+1)分别为点Aj和A(j+1)对应的横坐标值;XBj为点Bj对应的横坐标值;e为异常点补偿量。
2.根据权利要求1所述的叶片扭转量的测量方法,其特征在于:所述步骤S1中,所述传感器为位移传感器,所述传感器包括垂直于被测级叶片侧面的第一传感器、以及垂直于转子端面的第二传感器,所述第一传感器用于测量获取叶片实际位移量U,所述第二传感器用于测量获取转子轴向窜动量V;所测得的叶片实际位移量U和转子轴向窜动量V为离散的数值点;根据叶片实际位移量U和转子轴向窜动量V得到离散的叶片位移测量值Y、构成叶片位移测量值数列{Yn},叶片位移测量值数列{Yn}中的各项按照其测量的序号排列,n为自然数。
3.根据权利要求2所述的叶片扭转量的测量方法,其特征在于:所述步骤S1中,当所述第一传感器和第二传感器安装在被测级叶片的同侧时,叶片位移测量值Y等于同时刻的叶片实际位移量U和转子轴向窜动量V的差值;当所述第一传感器和第二传感器分别安装在被测级叶片的两侧时,叶片位移测量值Y等于同时刻的叶片实际位移量U和转子轴向窜动量V的和值。
4.根据权利要求2所述的叶片扭转量的测量方法,其特征在于:所述步骤S2中,所述叶片位移测量值数列{Yn}中极大值和极小值的判断标准为:所述极大值大于其前一个值、且大于或等于其后一个值,所述极小值小于其后一个值、且小于或等于其前一个值。
5.根据权利要求2所述的叶片扭转量的测量方法,其特征在于:所述步骤S2中得到极大值和极小值数量相同,得到所述极大值数列中的各项和极小值数列中的各项都按照测量的先后顺序排列。
6.根据权利要求1或5所述的叶片扭转量的测量方法,其特征在于:所述步骤S3中,所述异常极大值和异常极小值的判断方法为:算出极大值和极小值的极值平均值,将极大值数列和极小值数列中相同序号的极大值和极小值进行求差、得到极值差值,当该极值差值小于极值平均值的30%时,则该序号的极大值为异常极大值、该序号的极小值为异常极小值。
7.根据权利要求2所述的叶片扭转量的测量方法,其特征在于:所述步骤S1中,所述第一传感器测量的被测级叶片侧面为平面。
8.根据权利要求1所述的叶片扭转量的测量方法,其特征在于:所述异常点补偿量e的确定方式为:
(a)当点A(j+1)和点Bj之间无点Ci和点Di时,e=0;
(b)当点A(j+1)和点Bj之间有点Ci和点Di、且点Ci和点Di的纵坐标值Y不同时,根据点Ci和点Di在叶片位移测量值数列{Yn}中对应的序号确定补偿点Ei,e为补偿点Ei的横坐标值和点Ci横坐标值的差值;
(c)当点A(j+1)和点Bj之间有点Ci和点Di、且点Ci和点Di的纵坐标值Y相同时,e为点Ci的横坐标值和点Di的横坐标值的差值。
9.根据权利要求2所述的叶片扭转量的测量方法,其特征在于:所述步骤S1中,所述一圈被测级叶片中末叶片的围带顶部安装有感应块,所述传感器还包括第三传感器,当感应块随末叶片转动到第三传感器处时,所述第三传感器能够感应到感应块,并开始记录第一传感器和第二传感器的测量数据。
10.根据权利要求2所述的叶片扭转量的测量方法,其特征在于:所述被测级叶片具有依次相连的叶根、叶身和围带;将所述第一传感器设置在被测级叶片的叶根侧面,通过所述步骤S1至S4获得被测级叶片的叶根扭转量;将所述第一传感器设置在被测级叶片的围带侧面,通过所述步骤S1至S4获得被测级叶片的围带扭转量;将被测级叶片的叶根扭转量与围带扭转量相减,得到被测级叶片的相对扭转量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811450604.0A CN111256636B (zh) | 2018-11-30 | 2018-11-30 | 一种叶片扭转量的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811450604.0A CN111256636B (zh) | 2018-11-30 | 2018-11-30 | 一种叶片扭转量的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111256636A CN111256636A (zh) | 2020-06-09 |
CN111256636B true CN111256636B (zh) | 2021-11-19 |
Family
ID=70951809
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811450604.0A Active CN111256636B (zh) | 2018-11-30 | 2018-11-30 | 一种叶片扭转量的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111256636B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11248407A (ja) * | 1998-03-03 | 1999-09-17 | Mitsubishi Heavy Ind Ltd | 回転体の変形量計測装置 |
EP1295084A1 (en) * | 2000-06-28 | 2003-03-26 | Snap-on Technologies, Inc. | Method and apparatus for measuring vehicle wheel scrub radius |
CN104634288A (zh) * | 2013-11-06 | 2015-05-20 | 中国科学院沈阳计算技术研究所有限公司 | 压气机转子检测设备 |
CN105180860A (zh) * | 2015-08-31 | 2015-12-23 | 哈尔滨汽轮机厂有限责任公司 | 一种在线补偿式转子跳动测量系统及方法 |
CN205120017U (zh) * | 2015-11-27 | 2016-03-30 | 哈尔滨汽轮机厂有限责任公司 | 用于超临界汽轮机预扭型叶片装配的扭转量测量装置 |
CN106370142A (zh) * | 2016-11-15 | 2017-02-01 | 江苏方天电力技术有限公司 | 一种大量程位移传感器校准装置和校准方法 |
CA3027689A1 (en) * | 2016-07-11 | 2018-01-18 | Wobben Properties Gmbh | Measuring a torsion angle of a rotor blade |
CN107655429A (zh) * | 2016-08-04 | 2018-02-02 | 厦门乃尔电子有限公司 | 叶片振动位移的测量方法和装置 |
US10126119B2 (en) * | 2017-01-17 | 2018-11-13 | General Electric Company | Methods of forming a passive strain indicator on a preexisting component |
WO2018208485A1 (en) * | 2017-05-12 | 2018-11-15 | Siemens Energy, Inc. | Contactless, blade-tip clearance measurement for turbines |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103148764B (zh) * | 2013-03-11 | 2015-07-08 | 上海电气电站设备有限公司 | 汽轮机转子动叶围带间隙的测量方法及夹紧装置 |
FR3037394B1 (fr) * | 2015-06-09 | 2018-11-02 | Safran Aircraft Engines | Procede et dispositif de determination de la vibration d'aubes de rotor |
-
2018
- 2018-11-30 CN CN201811450604.0A patent/CN111256636B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11248407A (ja) * | 1998-03-03 | 1999-09-17 | Mitsubishi Heavy Ind Ltd | 回転体の変形量計測装置 |
EP1295084A1 (en) * | 2000-06-28 | 2003-03-26 | Snap-on Technologies, Inc. | Method and apparatus for measuring vehicle wheel scrub radius |
CN104634288A (zh) * | 2013-11-06 | 2015-05-20 | 中国科学院沈阳计算技术研究所有限公司 | 压气机转子检测设备 |
CN105180860A (zh) * | 2015-08-31 | 2015-12-23 | 哈尔滨汽轮机厂有限责任公司 | 一种在线补偿式转子跳动测量系统及方法 |
CN205120017U (zh) * | 2015-11-27 | 2016-03-30 | 哈尔滨汽轮机厂有限责任公司 | 用于超临界汽轮机预扭型叶片装配的扭转量测量装置 |
CA3027689A1 (en) * | 2016-07-11 | 2018-01-18 | Wobben Properties Gmbh | Measuring a torsion angle of a rotor blade |
CN107655429A (zh) * | 2016-08-04 | 2018-02-02 | 厦门乃尔电子有限公司 | 叶片振动位移的测量方法和装置 |
CN106370142A (zh) * | 2016-11-15 | 2017-02-01 | 江苏方天电力技术有限公司 | 一种大量程位移传感器校准装置和校准方法 |
US10126119B2 (en) * | 2017-01-17 | 2018-11-13 | General Electric Company | Methods of forming a passive strain indicator on a preexisting component |
WO2018208485A1 (en) * | 2017-05-12 | 2018-11-15 | Siemens Energy, Inc. | Contactless, blade-tip clearance measurement for turbines |
Non-Patent Citations (2)
Title |
---|
Blades condition monitoring using shaft torsional vibration signals;Al-Bedoor 等;《Journal of Quality in Maintenance Engineering》;20061231;第12卷(第3期);第275-93页 * |
基于LabVIEW的汽轮机转子动叶扭转量测量系统设计;吴春冬 等;《工具技术》;20150320(第3期);第82-86页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111256636A (zh) | 2020-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105333855B (zh) | 一种测量、计算涡轮导向叶片和导向器喉道面积的方法 | |
CN112067304B (zh) | 发动机整机试验中压气机进口流量的测量方法 | |
US10281297B2 (en) | Blade tip timing | |
CN104075677A (zh) | 航空发动机转子和静子同心度及叶尖间隙测量方法及系统 | |
KR20100038426A (ko) | 회전체의 심 어긋남 산출 방법 및 심 어긋남 산출 시스템 | |
CN103049596B (zh) | 一种提高叶轮零件逆向工程精度的方法 | |
CN108629092B (zh) | 一种基于尺寸效应修正的涡轮盘分区可靠性分析方法 | |
EP2808497A1 (en) | Blade tip timing | |
CN103292760A (zh) | 一种薄壁叶片误差分析方法 | |
CN108052717B (zh) | 一种基于局部应力-应变法的疲劳寿命校准方法 | |
CN105045973A (zh) | 一种变弧长自适应采样方法 | |
EP2818647B1 (en) | Blade tip timing system and method of analysing blade tip displacement | |
CN111256636B (zh) | 一种叶片扭转量的测量方法 | |
Key et al. | An experimental study of vane clocking effects on embedded compressor stage performance | |
CN112033681B (zh) | 一种加力燃烧室出口温度误差修正的方法 | |
CN104567748B (zh) | 狭小缝隙直线度与平面度的测量方法 | |
CN102680012B (zh) | 量测信息自动输出系统及方法 | |
CN108345573B (zh) | 用于汽轮机高压缸胀差测量探头的胀差确定函数计算方法 | |
CN113483720A (zh) | 一种基于蒙特卡洛的最小实体要求下直线度测量方法 | |
CN106525180B (zh) | 天然气压缩因子计算方法 | |
CN103994744B (zh) | 一种齿形测量方法 | |
CN117592313B (zh) | 一种形位误差测量不确定度仿真优化方法 | |
CN110906901A (zh) | 一种正六棱截面轴类偏心测量方法 | |
CN201653329U (zh) | 汽封间隙测量尺 | |
Lou et al. | Accounting for Circumferential Flow Nonuniformity in a Multi-Stage Axial Compressor |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |