WO2016091233A1 - Method of monitoring residual life of turbomachinery rotor blades - Google Patents

Method of monitoring residual life of turbomachinery rotor blades Download PDF

Info

Publication number
WO2016091233A1
WO2016091233A1 PCT/CZ2015/000149 CZ2015000149W WO2016091233A1 WO 2016091233 A1 WO2016091233 A1 WO 2016091233A1 CZ 2015000149 W CZ2015000149 W CZ 2015000149W WO 2016091233 A1 WO2016091233 A1 WO 2016091233A1
Authority
WO
WIPO (PCT)
Prior art keywords
damage
reconstructed
residual life
blade
rotor blades
Prior art date
Application number
PCT/CZ2015/000149
Other languages
English (en)
French (fr)
Inventor
Miroslav BALDA
Olga ČERVENÁ
Original Assignee
Výzkumný A Zkušební Ústav Plzeň S.R.O.
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Výzkumný A Zkušební Ústav Plzeň S.R.O. filed Critical Výzkumný A Zkušební Ústav Plzeň S.R.O.
Publication of WO2016091233A1 publication Critical patent/WO2016091233A1/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H1/00Measuring characteristics of vibrations in solids by using direct conduction to the detector
    • G01H1/003Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines
    • G01H1/006Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines of the rotor of turbo machines

Definitions

  • This invention relates to measurement and evaluation of vibration in solids.
  • BTT Band Tip Timing
  • the equipment is relatively simple in principle, despite the state- of-the-art devices used. It comprises a high-precision clock counter of very high frequency and two or more stator-mounted probes. The first of these is a reference probe which determines the start of each revolution, while the others monitor the passage of blade tips in front of them. As the machine rotates, the probes trigger pulses which govern the logging of instantaneous states of the counter to a memory device whose content is finally written to data files.
  • the records from a single measurement run take a binary form where the data from all probes is coded and stored in either a single file or in as many files as the probes.
  • the time histories in these files are ordinarily used for evaluating differences in blade passing times at the probes from which blade tip deflections are calculated, whose maxima are an indication of the vibration severity.
  • BTT systems are typically employed to monitor radial clearance between the rotor and the stator and to finding natural frequencies of blades in service.
  • the resulting signals are typically undersampled which is why the discrete Fourier transform of the time series of blade tip position samples does not yield a true Fourier amplitude spectrum but a spectrum with aliasing distortion in the frequency domain. This may result in some complications to the signal analysis, as described below.
  • Signal undersampling with all its consequences is known in this field from professional literature, e.g.: M. Balda: Lfvod do statisticki mechaniky. Zapadoceska univerzita v Plzni, 2001, ISBN 80-7082-820-X.
  • the ordinary BTT system includes an additional accelerometer which is attached to the machine's stator and transmits a composite signal that comprises rotational-frequency imbalance-induced signals and signals from the vibrating blades. Using appropriate filtering, rotational components are removed from this acceleration signal, leaving a signal which should contain blade vibration signals and, possibly, other still unidentified signals. Integer multiples of the sampling (rotational) frequency are added to or subtracted from the resonance frequencies of the aliased Fourier blade spectra, until a match is found of one of resonance frequencies from the accelerometer.
  • the system for generating damage estimates requires two simultaneous and independent measuring processes and an analysis of the discrete signal from the BTT system and the continuous accelerometer signal.
  • a signal filtering system is necessary for preparing the accelerometer signal for the subsequent analysis of stator vibrations.
  • the amplitude of the measured blade signal from the aliased Fourier spectrum can only be determined accurately if it is a line spectrum. However, this can only be carried out in exceptional cases where the resonance frequency lies in the middle of the frequency increment of df width in Fourier analysis. If that is not the case, the resonance frequency decomposes across several frequency bands and its value cannot be determined accurately.
  • the invention relates to a method of monitoring residual life of turbomachinery rotor blades.
  • the method uses records of sampled movement of blade tips obtained from clock counter data in a BTT system.
  • the records are decoded and clock counter overflows, mean value nonstationarities and signal trends are removed from the data.
  • the discrete Fourier transform is then applied to the data and its results for each rotating blade are translated into aliased amplitude spectra.
  • the aliased amplitude spectra are linked to the natural frequencies of blade vibrations which are known from calculations or found by testing on test equipment.
  • the discrete Fourier transform produces aliased amplitude spectra in frequency bands (-/N, +fii), where f N is Nyquist frequency, which is half the rotational frequency.
  • Resonance bands of the aliased amplitude spectra are then translated, while respecting the polarity, into reconstructed spectra of individual blades within the (-f & interval, where f R is half the sampling frequency of the reconstructed spectrum and its value is an order of magnitude larger than fy.
  • Damage increments caused by sequences of extrerna in time series of damaging stresses decomposed by the rainflow counting technique are then determined from the effective damaging stress values with the aid of a chosen damage rule. These damage increments are accumulated separately for each blade.
  • the accumulated damage increments are stored in processing logs for further use. Where required, the accumulated damage increments are used for plotting diagrams which provide information on the damage growth in and the residual life of each rotating blade.
  • the relative residual life of each blade is a complement of its accumulated relative damage to 1 , which corresponds to complete destruction of the blade.
  • the records are evaluated independently for each probe, up to the stage of damage increment calculation. From these results, the data from the probe which represents the largest damage increment will be used for evaluating the accumulated damage. This makes the estimate of accumulated damage a conservative result.
  • the damage rule for calculating the damage increments is either Palmgren-Miner rule or another rule derived from it.
  • the data can be displayed and analyzed during processing. This means that the operator may suspend the automatic processing sequence of measured BTT data, e.g. in order to transform the data in detail to a graphic form, and then return to the processing sequence without losing its results.
  • the turbomachine operator receives a warning at the time when a pre-defined blade damage extent is reached.
  • the initial damage level for the new replaced blades is set as zero.
  • the value of half the sampling frequency of the reconstructed spectrum is 25 times higher than fy.
  • the width of the frequency band of the reconstructed signal can be chosen by the user. The broader this band is, the better is the reconstruction of the extrema in the reconstructed signal, although it comes at the cost of higher computational demands.
  • the extrema in the reconstructed signal govern the intensity of the fatigue damage in the rotating blades.
  • Fig. 1 - shows moduli of the aliased Fourier vibration spectrum of a single blade sparsely sampled at the operating speed, as obtained from the discrete Fourier transformation of the data file generated by the BTT system.
  • the numbers at the top of the plot show the positions of aliased natural frequencies calculated using a formula from professional literature: M. Balda: Uvod do statisticke mechaniky. Zdpaduceskd univerzita v Plzni, 2001, ISBN 80-7082- 820-X;
  • Fig. 2 - shows moduli of reconstructed vibration spectrum of the same blade, in which natural frequency bands from the aliased spectrum shown in Fig. 1 were transferred to the correct frequencies within the reconstructed spectrum
  • Fig. 3 - shows a selected interval of the reconstructed blade tip movement signal derived from the reconstructed spectrum in Fig. 2 by means of inverse Fourier transform, and the points which represent the original sparse sampling of the original signal in the BTT file
  • Fig. 4 - illustrates how relative damage increased with time in individual blades over the course of a long single BTT measuring run.
  • the x axis indicates blade numbers
  • the y axis shows numerical codes of measurement intervals for a single probe, which also represent the time scale
  • the z axis represents relative damage increments.
  • the x - z plane contains relative damage increments for individual blades accumulated from all measurement intervals from the relevant probe. The figure clearly shows that substantial damage only occurs in some blades. These are the blades whose vibration is extensive, causing higher dynamic loads which, in turn, lead to fatigue crack initiation in critical locations.
  • Fig. 5 - is a graphical representation of accumulated relative damage in blades of the entire rotor. This figure shows, as does Fig. 4, that relative damage does not increase equally in all blades.
  • the relative residual life of a blade is a complement of the radius of the damage (dark) sector for the given blade to the radius of 1 of the circle. For a blade with very little damage, the circle with the radius of 1 which represents total relative damage can be outside the area captured in the graph.
  • the BTT system type is entered and, the input procedure is selected accordingly.
  • Blade parameters are defined, including the blade material, methods to be used for file processing, determination of frequency transfer functions relating the time differences and the loads in critical locations of the blades, and output types.
  • the subsequent processing of measurement data files takes place automatically without operator's intervention.
  • the measurement processing cycle is as follows: The output file(s) from the BTT system is/are read and, if required, decoded; upon long measurement runs, the files are split into shorter intervals, if needed.
  • results of the step described under point 2 are transformed using the discrete Fourier transform and an aliased Fourier spectrum is obtained for each probe and for each blade or blade measurement interval.
  • the aliased resonance frequencies are converted into real frequencies, as known in this field (see M, Balda: Uvod do statisticke mechaniky. ZdpadocesM univerzita v Plzni, 2001, ISBN 80-7082-820-X).
  • a concrete embodiment of the invention offers the processing operator a value of fa which is equal to 25 times the /N. The operator may accept it or enter their own value. Currently, the said value appears to be optimal, as it leads to a sufficient amount of data without being extremely computationally demanding.
  • the resonance bands are transferred from the aliased Fourier spectra to the correct frequencies in the reconstructed spectra determined according to point 4, while respecting the polarity.
  • the remainder of the aliased Fourier spectrum can be distributed across the reconstructed spectrum as desired by the user.
  • the currently implemented variant distributes frequency components uniformly across all bands of the reconstructed spectrum up to and including the highest resonance frequency.
  • the spectra reconstructed according to point 5 are converted using the inverse Fourier transformation to obtain reconstructed profiles which are close to the real movement of blade tips in units of time. If sparsely sampled signal data and densely sampled reconstructed signal data is plotted in diagrams, it becomes apparent that the original data is very close to the reconstructed signals and one can therefore expect a very good agreement between the evaluated and real damage levels.
  • the relative damage increments are accumulated separately for each blade, including the periods during which no measuring took place.
  • the accumulated relative damage levels are written into a file and, simultaneously, backups are created.
  • the PvFLB system waits for the BTT system to supply further measurement data.
  • the operator may interrupt the automatic processing cycle in order to deal with immediate tasks (e.g. in order to generate images from any portion of the processing sequence completed up to that point, to assess operating parameters with respect to the growth of relative damage, etc.) and, having completed these additional activities, return to the automatic operation cycle without losing any information.
  • immediate tasks e.g. in order to generate images from any portion of the processing sequence completed up to that point, to assess operating parameters with respect to the growth of relative damage, etc.
  • the present method allows sampled variable-amplitude signals to be processed as well, which represents a considerable strength because measuring runs involving constant amplitudes are very rare.
  • One of the extraordinary aspects of the method according to the present invention is that its processing time is shorter than the actual BTT measuring time, thanks to which earlier data can be processed without the new data building up.
  • the system enables the accumulation of relative damage to be continued for the remaining blades and to be started anew for the new blades.
  • the subject of the present invention will take various forms, differing in certain special aspects, most notably in the measurement data input. For this reason, the subject of the invention may take the following forms:
  • Figs. 1 through 5 The exemplary embodiment is shown in Figs. 1 through 5.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
PCT/CZ2015/000149 2014-12-11 2015-12-10 Method of monitoring residual life of turbomachinery rotor blades WO2016091233A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CZPV2014-893 2014-12-11
CZ2014-893A CZ305615B6 (cs) 2014-12-11 2014-12-11 Způsob monitorování zbytkové životnosti oběžných lopatek turbostrojů

Publications (1)

Publication Number Publication Date
WO2016091233A1 true WO2016091233A1 (en) 2016-06-16

Family

ID=55080294

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CZ2015/000149 WO2016091233A1 (en) 2014-12-11 2015-12-10 Method of monitoring residual life of turbomachinery rotor blades

Country Status (2)

Country Link
CZ (1) CZ305615B6 (cs)
WO (1) WO2016091233A1 (cs)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102019075B1 (ko) * 2018-04-19 2019-09-06 한국해양대학교 산학협력단 스펙트럼의 피로손상도 계산 최적화 시스템 및 그 방법
CN110567574A (zh) * 2019-08-02 2019-12-13 西安交通大学 一种旋转叶片叶端定时振动参数辨识方法与系统
CN111175033A (zh) * 2020-01-10 2020-05-19 中国人民解放军国防科技大学 一种变转速下的同步振动叶端定时信号重构方法
JP2021071354A (ja) * 2019-10-30 2021-05-06 株式会社日立製作所 軸受診断システム、および、軸受診断方法
CN113239613A (zh) * 2021-04-09 2021-08-10 国网新源控股有限公司 一种水轮机组摆度报警判断方法
CN113533530A (zh) * 2021-05-18 2021-10-22 西安交通大学 单个叶端定时传感器的叶片固有频率检测方法
CN119885792A (zh) * 2025-03-28 2025-04-25 肯佐控制设备(上海)有限公司 基于动态负载监测的电动执行机构寿命评估方法及系统
CN120257897A (zh) * 2025-06-06 2025-07-04 西北工业大学 基于积累功测量和疲劳寿命的非线性叶片颤振评估方法

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CZ307571B6 (cs) * 2017-06-27 2018-12-12 Ăšstav termomechaniky AV ÄŚR, v.v.i. Způsob identifikace poškozené lopatky lopatkového stroje za provozu
CN112084583B (zh) * 2020-07-24 2024-04-05 西安交通大学 基于数字孪生的转子叶片寿命预测方法和预测系统
CN113533529B (zh) * 2021-05-18 2022-10-28 西安交通大学 单个或均布叶端定时传感器提取叶片间固有频率差值方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2466817A (en) 2009-01-09 2010-07-14 Rolls Royce Plc Determining rotor blade vibration frequencies and/or amplitudes by combining blade tip timing measurements and rotor support structure vibration measurements

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7245786B2 (en) * 2002-05-10 2007-07-17 976076 Alberta Inc. Filtering artifact from fMRI data using the stockwell transform
CZ303476B6 (cs) * 2009-10-15 2012-10-10 Štarman@Stanislav Zpusob bezkontaktního monitorování turbín, zejména jednotlivých lopatek parní nebo plynové turbíny v elektrárnách, a systém k provádení tohoto zpusobu

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2466817A (en) 2009-01-09 2010-07-14 Rolls Royce Plc Determining rotor blade vibration frequencies and/or amplitudes by combining blade tip timing measurements and rotor support structure vibration measurements
US20100179775A1 (en) * 2009-01-09 2010-07-15 Rolls-Royce Plc Determination of blade vibration frequencies and/or amplitudes

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
ANDREW HALFPENNY: "A Frequency Domain Approach for Fatigue Life Estimation from Finite Element Analysis", KEY ENGINEERING MATERIALS, vol. 167-168, 30 June 1999 (1999-06-30), pages 401 - 410, XP055257789, DOI: 10.4028/www.scientific.net/KEM.167-168.401 *
M BALDA, CLVOD DO STATISTICKE MECHANIKY, 2001
M BALDA, ÍVOD DO STATISTICKE MECHANIKY, 2001
M BALDA; J. SVOBODA; 1! FRÖHLICH: "Proc. Conf Engineering Mechanics 2007", 2007, article "An estimation of fatigue life under general stress"
M. BALDA, ÚVOD DO STATISTICKE MECHANIKY, 2001
M. BALDA,; J. SVOBODA; V FROHLICH: "Proc. Conf Engineering /'vlechanics 2007", 2007, IT CAS, article "An estimation of fatigue life under general stress."
MIROSLAV BALDA: "Prediction of Damage Cumulation in Vibrating Rotors", FIFTH INTERNATIONAL CONFERENCE ON ROTOR DYNAMICS DARMSTADT, 10 September 1998 (1998-09-10), XP055257794 *
MIROSLAV BALDA: "Processing of Blade Monitoring System Data", PROCEEDINGS OF THE 7TH INTERNATIONAL CONFERENCE ON VIBRATION IN ROTATING MACHINES, 14 September 2000 (2000-09-14), XP055257660, Retrieved from the Internet <URL:http://www.cdm.cas.cz/publications/balda/Nottingham00.pdf> [retrieved on 20160311] *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102019075B1 (ko) * 2018-04-19 2019-09-06 한국해양대학교 산학협력단 스펙트럼의 피로손상도 계산 최적화 시스템 및 그 방법
CN110567574A (zh) * 2019-08-02 2019-12-13 西安交通大学 一种旋转叶片叶端定时振动参数辨识方法与系统
CN110567574B (zh) * 2019-08-02 2020-11-06 西安交通大学 一种旋转叶片叶端定时振动参数辨识方法与系统
JP2021071354A (ja) * 2019-10-30 2021-05-06 株式会社日立製作所 軸受診断システム、および、軸受診断方法
CN111175033A (zh) * 2020-01-10 2020-05-19 中国人民解放军国防科技大学 一种变转速下的同步振动叶端定时信号重构方法
CN113239613A (zh) * 2021-04-09 2021-08-10 国网新源控股有限公司 一种水轮机组摆度报警判断方法
CN113239613B (zh) * 2021-04-09 2024-05-31 国网新源控股有限公司 一种水轮机组摆度报警判断方法
CN113533530A (zh) * 2021-05-18 2021-10-22 西安交通大学 单个叶端定时传感器的叶片固有频率检测方法
CN119885792A (zh) * 2025-03-28 2025-04-25 肯佐控制设备(上海)有限公司 基于动态负载监测的电动执行机构寿命评估方法及系统
CN120257897A (zh) * 2025-06-06 2025-07-04 西北工业大学 基于积累功测量和疲劳寿命的非线性叶片颤振评估方法

Also Published As

Publication number Publication date
CZ2014893A3 (cs) 2016-01-06
CZ305615B6 (cs) 2016-01-06

Similar Documents

Publication Publication Date Title
WO2016091233A1 (en) Method of monitoring residual life of turbomachinery rotor blades
KR101718251B1 (ko) 회전 블레이드 강건성 모니터링 방법 및 시스템
EP2136189B1 (en) Method for analysing vibration in rotor blades
EP2072975A1 (en) Method and apparatus for vibration-based automatic condition monitoring of a wind turbine
US9016132B2 (en) Rotating blade analysis
US9657588B2 (en) Methods and systems to monitor health of rotor blades
EP2345894B1 (en) Trending of vibration data taking into account torque effect
CN108709635B (zh) 确定旋转机械振动信号主要频率成分信息的方法及装置
CN111397877B (zh) 一种旋转机械拍振故障检测与诊断方法
WO2018002819A1 (en) A method and system for measuring rotor blade tip deflection using blade tip timing (btt)
US20110291850A1 (en) Machine vibration monitoring
JP2015505058A5 (cs)
WO2021049068A1 (ja) 振動検出装置、振動検出方法および異常判定システム
Neri et al. Non-harmonic fourier analysis for bladed wheels damage detection
EP4139694A1 (en) Rotating machine speed estimation
JP2015125147A (ja) ロータブレードの健全性をモニターする方法およびシステム
CN118936894A (zh) 透平机械转子振动全息监测及动平衡方法、设备和介质
RU2670771C9 (ru) Способ определения характера касания лопатки вращающегося колеса о корпус турбомашины
RU2584723C1 (ru) Способ определения параметров колебаний лопаток вращающегося колеса турбомашины и устройство для его осуществления
CN115683644A (zh) 航空发动机双源拍振特征识别方法
Oberholster et al. Eulerian laser Doppler vibrometry: Online blade damage identification on a multi-blade test rotor
RU2659428C1 (ru) Устройство для анализа динамических процессов в рабочих колесах турбомашин
Alekseev et al. Data measurement system of compressor units defect diagnosis by vibration value
Rajendran et al. Wavelet Transform‐Based Damage Identification in Bladed Disks and Rotating Blades
Procházka et al. Sensors and methods for blade damage operational assessment in low-pressure steam turbine stages

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15823490

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 15823490

Country of ref document: EP

Kind code of ref document: A1