CZ306833B6 - Způsob detekce a lokalizace částečného kontaktu rotor-stator při provozu turbíny - Google Patents
Způsob detekce a lokalizace částečného kontaktu rotor-stator při provozu turbíny Download PDFInfo
- Publication number
- CZ306833B6 CZ306833B6 CZ2012-947A CZ2012947A CZ306833B6 CZ 306833 B6 CZ306833 B6 CZ 306833B6 CZ 2012947 A CZ2012947 A CZ 2012947A CZ 306833 B6 CZ306833 B6 CZ 306833B6
- Authority
- CZ
- Czechia
- Prior art keywords
- rotor
- frequency
- signals
- turbine
- stator
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
- G01M13/045—Acoustic or vibration analysis
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D11/00—Preventing or minimising internal leakage of working-fluid, e.g. between stages
- F01D11/08—Preventing or minimising internal leakage of working-fluid, e.g. between stages for sealing space between rotor blade tips and stator
- F01D11/14—Adjusting or regulating tip-clearance, i.e. distance between rotor-blade tips and stator casing
- F01D11/20—Actively adjusting tip-clearance
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D19/00—Starting of machines or engines; Regulating, controlling, or safety means in connection therewith
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D21/00—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
- F01D21/003—Arrangements for testing or measuring
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D21/00—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
- F01D21/04—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for responsive to undesired position of rotor relative to stator or to breaking-off of a part of the rotor, e.g. indicating such position
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D25/00—Component parts, details, or accessories, not provided for in, or of interest apart from, other groups
- F01D25/04—Antivibration arrangements
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D25/00—Component parts, details, or accessories, not provided for in, or of interest apart from, other groups
- F01D25/16—Arrangement of bearings; Supporting or mounting bearings in casings
- F01D25/162—Bearing supports
- F01D25/164—Flexible supports; Vibration damping means associated with the bearing
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02C—GAS-TURBINE PLANTS; AIR INTAKES FOR JET-PROPULSION PLANTS; CONTROLLING FUEL SUPPLY IN AIR-BREATHING JET-PROPULSION PLANTS
- F02C9/00—Controlling gas-turbine plants; Controlling fuel supply in air- breathing jet-propulsion plants
- F02C9/26—Control of fuel supply
- F02C9/28—Regulating systems responsive to plant or ambient parameters, e.g. temperature, pressure, rotor speed
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
- G01H1/003—Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines
- G01H1/006—Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines of the rotor of turbo machines
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/14—Testing gas-turbine engines or jet-propulsion engines
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
Vynález se týká způsobu detekce a lokalizace částečného kontaktu rotor-stator při provozu turbíny, u kterého se rotor (1) turbíny snímá alespoň jedním párem snímačů (3, 30, 4, 40) rotorových vibrací umístěným v jedné měřicí rovině (330, 440) kolmé k ose rotoru (1) a alespoň jedním párem snímačů absolutních vibrací ložisek (2, 20). U tohoto způsobu se nejprve u synchronně digitalizovaných signálů snímačů (3, 30, 4, 40) rotorových vibrací vypočte aritmetický průměr jednotlivých harmonických pásem v úplném spektru ve frekvenční oblasti, čímž se převedou informace o subharmonických složkách a jejich násobcích v úplném spektru těchto signálů do frekvenčního pásma -1X až 1X (kde X označuje otáčkovou frekvenci rotoru (1)). Zároveň se výpočtem mediánu odhadne úroveň šumu, kterým jsou zatížena jednotlivá frekvenční pásma úplného spektra, a která se od amplitud v jednotlivých frekvenčních pásmech úplného spektra signálů snímačů rotorových vibrací odečte. Tím se vytvoří normované akumulované úplné spektrum signálů snímačů rotorových vibrací, ve kterém se následně porovnává amplituda předem stanovené subharmonické složky/složek signálů rotorových vibrací s amplitudami frekvenčních složek v jejím/jejich okolí stanovenými Fourierovou transformací nebo aproximací Fourierovou řadou. Přitom se usuzuje na přítomnost sledované subharmonické složky signálů rotorových vibrací ve spektru těchto signálů, na základě čehož se následně usuzuje na existenci kontaktu rotor-stator. V případě existence tohoto kontaktu se dále lokalizuje místo/místa turbíny, kde k němu dochází.
Description
Oblast techniky
Vynález se týká způsobu detekce a lokalizace částečného kontaktu rotor-stator při provozu turbíny, u kterého se rotor turbíny snímá alespoň jedním párem snímačů rotorových vibrací umístěným v jedné měřicí rovině kolmé k ose rotoru.
Dosavadní stav techniky
Při provozu parních i plynových turbín může v některých situacích, např. během najíždění turbíny, kdy při překonávání vlastních frekvencí rotoru jsou jeho vibrace většinou nejvyšší, dojít k nežádoucímu kontaktu mezi rotorem a statorem („stator-rotor rub“ nebo „rubbing“). Během tohoto kontaktu se v první fázi především obrušují ucpávky uložené mezi rotujícími a statickými částmi turbíny, v důsledku čehož se zvyšuje množství unikajícího média a snižuje účinnost turbíny. V případech, kdy není kontakt rotor-stator včas detekován a vhodným zásahem do provozu turbíny, např. změnou otáček při jejím nájezdu nebo doběhu, nebo změnou parametrů zvedacího oleje při jejím provozu na natáčecím zařízení, odstraněn, může dojít k těžkému poškození rotoru i statoru, resp. celkové havárii turbíny se značnými ekonomickými škodami.
Kontakt rotor-stator může být buď částečný, kdy se jedná o velmi krátký, avšak alespoň několikrát opakovaný kontakt mezi rotorem a statorem, nebo úplný, kdy se jedná o kontinuální nebo téměř kontinuální kontakt mezi rotorem a statorem. Úplný kontakt rotor-stator je přitom vždy předcházen, alespoň na krátký okamžik, částečným kontaktem rotor-stator.
V současné době je detekce částečného kontaktu rotor-stator založena zejména na offline analýze vibračních signálů, kdy obsluha zařízení sleduje v průběhu měření celkovou úroveň vibrací a také fázor první harmonické složky otáčkové frekvence ve vibračních signálech. Pokud dojde ke skokové změně celkových vibrací nebo rotaci fázoru první harmonické složky s proměnnou nebo periodicky se měnící amplitudou, je detekován kontakt rotor-stator, a po ukončení měření je provedena detailní analýza dat pro vyloučení možnosti falešně pozitivní detekce. Tento přístup tak není vhodný pro detekci kontaktu rotor-stator během reálného provozu turbíny, a je využitelný v podstatě pouze pro laboratorní či experimentální účely.
Další přístup k detekci částečného kontaktu rotor-stator je založen na tom, že kromě změny fázoru první harmonické složky vibračních signálů je tento kontakt doprovázen také vznikem subharmonických složek spektra, jejichž frekvence odpovídá četnosti nárazů rotoru na stator. Nevýhodou detekce kontaktu rotor-stator na základě sledování těchto subharmonických složek je však to, že frekvence některých z nich jsou velmi blízké frekvencím odpovídajícím projevům jiných závad - např. nestabilitě olejového filmu kluzného ložiska, apod., v důsledku čehož může dojít buď k falešně pozitivní detekci kontaktu rotor-stator, nebo naopak k chybnému vyhodnocení projevů tohoto kontaktu jako závady jiného typu.
Kromě toho jsou z patentové i nepatentové literatury dále známé také jiné způsoby pro detekci a/nebo lokalizaci částečného kontaktu rotor-stator založené na zcela odlišných principech.
Např. z EP 1 533 479 je známý způsob detekce kontaktu rotor-stator, u kterého se sledují provozní veličiny turbíny, jako např. teplota pláště, excentricita rotoru, tlak v kondenzátoru a zatížení generátoru. Příslušný algoritmus přitom vyhodnocuje samostatně změnu jednotlivých veličin, a pokud u některé z nich dojde k abnormální, resp. skokové, změně, je tato změna vyhodnocena jako průvodní jev kontaktu rotor-stator. Nevýhodou tohoto způsobuje, že skoková změna sledované veličiny může souviset i s jinými jevy nebo závadami, ke kterým během provozu turbíny dochází, takže existuje velmi vysoká pravděpodobnost falešně pozitivní detekce kontaktu rotor-1 CZ 306833 B6 stator. Vzhledem k dynamice sledovaných veličin a předpokládané citlivosti použitého algoritmu navíc nelze vyhodnocování změny těchto veličin a následnou detekci kontaktu rotor-stator považovat za dostatečně rychlý způsob reálně použitelný během provozu turbíny.
Z US 2008/240902 je pak známý způsob detekce kontaktu rotor-stator, který spočívá ve sledování teploty statoru turbíny alespoň na části jeho obvodu, přičemž případné lokální zvýšení teploty je přisouzeno tření mezi rotorem a statorem, a je vyhodnoceno jako přítomnost kontaktu rotor-stator. Dle umístění teplotního snímače, který zvýšení teploty zaznamenal, je navíc možné místo kontaktu na statoru hrubě lokalizovat. Nevýhodou tohoto postupu je, že ve chvíli, kdy je třením mezi rotorem a statorem vytvořeno dostatečně velké množství tepla, aby bylo tepelným snímačem zaznamenáno, může od jejich prvního kontaktu uplynout doba až v řádu minut, během které už může být rotor, stator a případně i další součásti turbíny značně poškozeny, a ohýbání rotoru již nelze zamezit jinak než odstavením stroje. Další nevýhodou je, že v případě závažnějšího kontaktu rotor-stator, ke kterému dochází po celém obvodu statoru nebo jeho větší části, je zvýšení teploty statoru zaznamenáno u většiny nebo u všech tepelných snímačů umístěných ve statoru v podstatě rovnoměrně, díky čemuž nemusí být vždy vyhodnoceno jako kontakt rotor-stator. Podstatnou nevýhodou je i to, že dodatečná instalace většího množství teplotních snímačů do tělesa statoru turbíny vyžaduje odstávku turbíny a poměrně značné investice.
Kromě toho je z US 2009/0 003 991 dále známý způsob detekce kontaktu rotor-stator, který spočívá v kontinuálním měření vzdálenosti lopatek rotoru od statoru a/nebo výšky lopatek rotoru. Nevýhodou tohoto způsobuje stejně jako v předcházejícím případě nutnost dodatečné instalace speciálních senzorů do statoru turbíny a s tím související zvýšení celkových nákladů. Měřením změny velikosti prostoru mezi lopatkou rotoru a lopatkovou ucpávkou ve statoru pak sice lze detekovat kontakt mezi lopatkou a statorem, avšak vzhledem k tomu, že vůle mezi lopatkou a statorem jsou větší než vůle mezi hřídelí a statorem, a ke kontaktu rotor-stator tak dochází primárně na hřídeli turbíny, reálná pravděpodobnost úspěšné a včasné detekce kontaktu rotor-stator tímto způsobem je jen velmi malá. Další nevýhodou je možnost falešně pozitivní detekce kontaktu rotor-stator při změně senzory snímaného rozměru, ke které však může dojít i zjiného důvodu, např. v důsledku excentricity rotoru, změnou rychlosti jeho rotace, apod.
US 4 478 082 popisuje zařízení pro detekci kontaktu rotor-stator na základě vyhodnocení akustických signálů snímaných alespoň jedním akustickým snímačem umístěným na rotoru nebo statoru rotačního stroje, např. turbíny. Nevýhodou tohoto řešení je nejen to, že rotační stroj musí být dodatečně vybaven tímto novým, jednoúčelovým snímačem a jeho příslušenstvím, ale také to, že akustický signál, který se následně vyhodnotí jako kontakt rotor-stator, může být ve skutečnosti způsoben jiným zdrojem nebo jevem, ke kterému dochází při provozu turbíny.
Stejná je situace také v případě US 4 377 947, který popisuje podobné zařízení, které je opatřeno alespoň dvěma akustickými snímači.
WO 2008/060 164 pak popisuje způsob pro detekci kontaktu rotor-stator u turbosoustrojí, který je založený na použití akcelerometru, snímače vibrací vysokých frekvencí nebo mikrofonu, přičemž se signál získaný některým z těchto snímačů následně srovnává s předem stanovenými hodnotami nebo vzory, a na základě tohoto srovnání se usuzuje na přítomnost nebo nepřítomnost kontaktu rotor-stator. Tento způsob je jedním z výše popsaných způsobů pro offline analýzu vibračních signálů.
Článek PATEL T H ET AL: “Coupled bending-torsional vibration analysis if rotor with rub and crack”, JOURNAL OF SOUND & VIBRATION, LONDON, GB, vol. 326, no. 3-5, 9 October 2009, pages 740-752, XP026422078, ISSN: 0022-460X popisuje způsob modelování a analýzy vibračního vzorce rotoru při kontaktu rotor-stator a případně i existenci příčné praskliny hřídele. Použitá Hilbert-Huangova transformace pak v kombinaci s výpočtem okamžitých frekvencí umožňuje zobrazovat tyto nelinearity v časo-frekvenčním spektrogramu.
-2CZ 306833 B6
Cílem vynálezu je odstranit nevýhody stavu techniky a navrhnout nový způsob detekce částečného kontaktu rotor-stator během provozu turbíny, který by odstranil nevýhody stavu techniky, a který by umožnil detekci kontaktu rotor-stator v reálném čase, co nejdříve po jeho vzniku, a současně i jeho co nejpřesnější lokalizaci.
Podstata vynálezu
Cíle vynálezu se dosáhne způsobem detekce a lokalizace částečného kontaktu rotor-stator při provozu turbíny, u kterého se rotor turbíny snímá alespoň jedním párem snímačů rotorových vibrací umístěným v jedné měřicí rovině kolmé kose rotoru a alespoň jedním párem snímačů absolutních vibrací ložisek. Podstata tohoto vynálezu přitom spočívá v tom, že u synchronně digitalizovaných signálů snímačů rotorových vibrací se vypočte aritmetický průměr jejich harmonických pásem vymezených otáčkovou frekvencí X rotoru turbíny a jejími vyššími harmonickými násobky v úplném spektru ve frekvenční oblasti, čímž se převedou informace o subharmonických složkách a jejich násobcích v úplném spektru těchto signálů do frekvenčního pásma -IX až IX (kde X označuje otáčkovou frekvenci rotoru), zároveň se výpočtem mediánu odhadne úroveň šumu, kterým jsou zatížena jednotlivá frekvenční pásma úplného spektra, a který se od amplitud v jednotlivých frekvenčních pásmech úplného spektra signálů snímačů rotorových vibrací odečte, čímž se vytvoří normované akumulované úplné spektrum signálů snímačů rotorových vibrací, ve kterém se porovnává amplituda předem stanovené subharmonické složky/složek signálů rotorových vibrací s amplitudami frekvenčních složek v jejím/jejich okolí stanovenými Fourierovou transformací nebo aproximací Fourierovou řadou. Vizuálním hodnocením ze strany obsluhy a/nebo automatickým vyhodnocením se přitom určuje přítomnost sledované subharmonické složky signálů rotorových vibrací ve spektru těchto signálů, která je znakem existence kontaktu rotorstator.
V případě existence kontaktu rotor-stator se dále lokalizuje místo/místa turbíny, kde k tomuto kontaktu dochází. Při axiální lokalizaci se zpracují úseky signálů snímačů absolutních vibrací ložisek turbíny v okamžiku vzniku, resp. průběhu kontaktu rotor-stator, v časo-frekvenční oblasti, ve které se identifikují širokopásmová amplitudová navýšení těchto signálů a následně se určí časy nárůstu amplitud signálů nad úroveň šumu, a ze znalosti těchto časů a jejich časových diferencí, rozmístění snímačů absolutních vibrací ložisek, rychlosti šíření rázu v materiálu rotoru/statoru se axiální lokalizací určí axiální poloha místa/míst kontaktu rotor-stator, a srovnáním histogramu absolutních četností těchto míst s výkresem nebo osou turbíny se určí axiální poloha místa/míst, kde ke kontaktu rotor-stator skutečně dochází. Pro obvodovou lokalizaci kontaktu rotor-stator se amplitudová navýšení signálů snímačů absolutních vibrací ložisek zobrazí do orbity 1. harmonické složky otáčkové frekvence rotoru vytvořené ze signálů snímačů rotorových vibrací, čímž se určí úhlová poloha maxima amplitudových navýšení způsobených kontaktem rotorstator, která odpovídá obvodové poloze místa/míst, kde dochází ke kontaktu rotor-stator.
Výhodné je přitom použití snímačů relativních rotorových vibrací, které jsou standardní součástí polní instrumentace turbín. Kromě nich nebo v kombinaci s nimi však lze použít i snímače absolutních rotorových vibrací.
Při vizuálním vyhodnocování sleduje obsluha amplitudu předem stanovené subharmonické složky/složek signálů rotorových vibrací v intervalu 1/9X až 1/2X, kterou porovnává s amplitudami frekvenčních složek v jejím/jejich okolí, přičemž pokud amplituda této subharmonické složky/složek převýší amplitudy okolních frekvenčních složek, znamená to, že sledovaná subharmonická složka/složky je/jsou přítomná/přítomné ve spektru signálů snímačů rotorových vibrací, a že dochází ke kontaktu rotor-stator.
Při automatickém vyhodnocování se s výhodou stanoví charakteristická veličina subharmonických složek s frekvencemi rovnými celočíselným podílům otáčkové frekvence rotoru, přičemž pokud ke kontaktu rotor-stator nedochází, nabývá hodnota této charakteristické veličiny pro
-3 CZ 306833 B6 danou subharmonickou složku hodnoty blízké 1, a pokud ke kontaktu rotor-stator dochází, nabývá hodnota charakteristické veličiny pro danou subharmonickou složku hodnotu převyšující 1.
Charakteristická veličina je definována vztahem fsf (0 = max
kde Acaj(t, ±sf) označuje amplitudu normovaného akumulovaného spektra pro subharmonickou frekvenci Sf v čase t a Acqf označuje amplitudu frekvenční linie v blízkém okolí této frekvence, M je celkový počet frekvenčních linií, které se použijí k průměrování, Aje index součtu a fk je hodnota frekvenční linie, která se nachází v okolí subharmonické frekvence Sf.
V případě, že ke kontaktu rotor-stator nedochází, lze charakteristickou veličinu chápat jako veličinu náhodnou a aproximovat ji součtem její střední hodnoty μ a směrodatné odchylky σ vypočtených rekurzivním algoritmem. Na existenci kontaktu rotor-stator se pak usuzuje ze změny střední hodnoty μ a směrodatné odchylky σ.
Při axiální lokalizaci místa/míst kontaktu rotor-stator se úseky signálů snímačů absolutních vibrací ložisek turbíny v okamžiku vzniku, resp. průběhu kontaktu rotor-stator, s výhodou zpracují vínkovou transformací, ve které se jako bázová funkce použije Morletova vlnka.
Objasnění výkresů
Na obr. 1 je schematicky bez lopatek a dalších konstrukčních prvků znázorněn hřídel turbíny, na kterém je aplikován způsob detekce částečného kontaktu rotor-stator při provozu turbíny podle vynálezu, na obr. 2 příkladný úplný spektrogram signálů snímačů relativních rotorových vibrací uložených ve dvou různých měřicích rovinách, na obr. 3 příkladný průběh úplného spektra vypočteného ze signálů snímačů relativních rotorových vibrací uložených v jedné z měřicích rovin, na obr. 4 příkladný spektrogram akumulovaného úplného spektra vypočteného ze signálů snímačů relativních rotorových vibrací v jedné z měřicích rovin, na obr. 5 průběh úplného spektra signálů snímačů relativních rotorových vibrací podle levé části obr. 2 v řádové oblasti, na obr. 6 srovnání průběhu úplného spektra signálů snímačů relativních rotorových vibrací podle obr. 5 s průběhem normovaného úplného spektra signálů snímačů relativních rotorových vibrací v řádové oblasti, na obr. 7 příkladný spektrogram normovaného akumulovaného úplného spektra vypočteného ze signálů snímačů relativních rotorových vibrací v jedné z měřicích rovin, na obr. 8 příkladný průběh charakteristických veličin pro subharmonické frekvence 1/2 a 1/3 otáčkové frekvence rotoru, na obr. 9 příkladný spektrogram spektra signálů snímačů absolutních vibrací ložiska turbíny v okamžiku kontaktu rotor-stator, na obr. 10 příkladný výkres turbíny s přeneseným histogramem absolutních četností kontaktu rotor-stator a na obr. 11 příkladná úhlová poloha maxima amplitudových změn způsobených částečným kontaktem rotor-stator zobrazená do orbity 1. harmonické složky otáčkové frekvence rotoru turbíny.
Příklady uskutečnění vynálezu
Způsob detekce částečného kontaktu rotor-stator při provozu turbíny a jeho následné lokalizace podle vynálezu bude popsán na příkladu turbíny schematicky znázorněné na obr. 1. Rotor 1 turbíny je otočně uložen ve dvou radiálních ložiscích 2 a 20, přičemž v blízkosti každého z nich je snímán dvojicí snímačů 3, 30 a 4, 40 relativních rotorových vibrací. Každá dvojice je přitom uložena ve společné měřicí rovině 330 a 440 kolmé k ose 10 rotoru 1, ve stejné vzdáleností od ní, a snímače 3, 30 a 4, 40 jsou v příslušné měřicí rovině 330, resp. 440 vzájemně posunuly, s výhodou o úhel 90°. Kromě toho jsou radiální ložiska 2 a 20 turbíny snímána neznázoměnými
-4CZ 306833 B6 snímači absolutních vibrací těchto ložisek. Všechny tyto snímače jsou běžnou součástí polní instrumentace turbín, takže způsob detekce částečného kontaktu rotor-stator a jeho lokalizace při provozu turbíny podle vynálezu nevyžaduje dodatečnou instalaci žádných speciálních součástí. Další jeho výhodou je, že umožňuje detekci kontaktu rotor-stator online v reálném čase, v podstatě okamžitě po jeho vzniku.
Signály snímačů 3, 30, 4, 40 relativních rotorových vibrací lze v dalších variantách doplnit a/nebo nahradit signály neznázoměných snímačů absolutních rotorových vibrací, které však u stávajících instalovaných turbín nejsou standardně používány ve dvojicích, případně nejsou používány vůbec, takže v takovém případě je nutné turbínu, resp. její stator těmito snímači osadit.
Pro další zpracování se využijí vždy signály alespoň jedné z dvojic snímačů 3, 30 a/nebo 4, 40 rotorových vibrací umístěných v jedné měřicí rovině 330, 440. Signály těchto snímačů 3, 30 a/nebo 4, 40 se přitom synchronně digitalizují neznázorněným známým A/D převodníkem/převodníky.
U takto vytvořených digitálních signálů snímačů 3, 30, 4, 40 rotorových vibrací se následně vyšetřuje přítomnost, resp. nepřítomnost subharmonických složek s frekvencí rovnou celočíselným podílům otáčkové frekvence rotoru 1 v jejich spektru, neboť tyto složky jsou zřejmým projevem částečného kontaktu rotor-stator. V případě přítomnosti takových složek se s výhodou dále sleduje také jejich frekvence, která vypovídá o intenzitě kontaktu rotor-stator.
Základním způsobem detekce přítomnosti subharmonických složek signálů snímačů 3, 30, 4, 40 rotorových vibrací v jejich spektru je přitom sledování změn amplitud těchto složek v úplném spektru. Tyto amplitudy se přitom s výhodou stanoví Fourierovou transformací signálů snímačů 3, 30, 4, 40 rotorových vibrací nebo na základě aproximace těchto signálů Fourierovou řadou. Takto získané amplitudy se následně srovnávají s předem stanoveným prahem (viz charakteristická veličina dále), přičemž v případě, že ho překročí, usuzuje se na přítomnost částečného kontaktu rotor-stator.
Vzhledem ktomu, že částečný kontakt rotor-stator má obvykle proměnnou intenzitu, což se projevuje změnou frekvence subharmonických složek signálů snímačů 3, 30, 4, 40 rotorových vibrací např. z 1/3 a 2/3 otáčkové frekvence rotoru 1 (v daném příkladě je otáčková frekvence 8 Hz) na 1/2 otáčkové frekvence rotoru 1, což je zřejmé např. z pravé části obr. 2 (pro měřicí rovinu 330), zatímco jiné subharmonické složky signálů snímačů 3, 30, 4, 40 rotorových vibrací mají frekvenci např. v oblasti -3 až -4 násobku otáčkové frekvence rotoru 1 - viz levá strana obr. 2 (pro měřicí rovinu 440), je vhodné přenést informaci o subharmonických složkách v úplném spektru signálů snímačů 3, 30, 4, 40 rotorových vibrací z tohoto širokého frekvenčního pásma do užšího intervalu a sledovat jejich přítomnost, resp. chování v něm. Výhodný je například interval -IX až IX (kde X označuje otáčkovou frekvenci rotoru 1), který umožňuje určit o jakou subharmonickou složku signálů snímačů 3, 30, 4, 40 rotorových vibrací vzhledem k otáčkové frekvenci rotoru 1 se jedná. Klíčovými složkami pro detekci kontaktu rotor-stator jsou přitom složky s frekvencí rovnající se celočíselným podílům otáčkové frekvence rotoru f. Na obr. 3 je zobrazen příkladný průběh úplného spektra vypočteného ze signálů snímačů relativních rotorových vibrací v měřicí rovině 440, přičemž hodnoty amplitudového spektra jsou zobrazeny v decibelech. Svislými přerušovanými čarami je pak toto spektrum rozděleno do frekvenčních intervalů podle násobků otáčkové frekvence rotoru 1, přičemž, jak je patrné, subharmonické složky spektra a jejich násobky jsou obsaženy ve všech těchto frekvenčních intervalech. Konkrétně ve znázorněném příkladu provedení se přitom jedná o komponenty signálů s frekvencí rovnou 1/2X a jejím celočíselným násobkům (na obr. 3 označené hodnotami 1/2X až 7/2X). K přenesení informací o subharmonických složkách signálů do užšího frekvenčního intervalu se s výhodou použije metoda akumulovaného úplného spektra signálů.
Akumulované úplné spektrum signálů je přitom možno vyjádřit následujícím vztahem:
-5CZ 306833 B6
Va-1
Λ/ (t, ±a>) = — X Af (z, ±a>t), 7Va k=0 = (ka>r(ť), (k + 1)ωΓ(0), kde Aaft,±o>) je akumulované úplné spektrum v čase t, Α/ί,±ω) je úplné spektrum v čase t, a»r(t) je otáčková frekvence rotoru 1 v čase t, a Na počet intervalů, do kterých byla levá a pravá část úplného spektra rozdělena (obr. 3). Na se přitom na základě dlouhodobých zkušeností s používanými snímači zejména relativních rotorových vibrací a experimentálního ověření s výhodou volí v intervalu 5 až 10.
Výpočet akumulovaného úplného spektra tak spočívá ve výpočtu aritmetického průměru harmonických frekvenčních pásem signálů snímačů 3, 30, 4, 40 rotorových vibrací z úplného spektra ve frekvenční oblasti, a umožňuje rozhodnout o vztahu subharmonických složek těchto signálů k otáčkové frekvenci rotoru f Jeho podstatou je odhalení násobků základních subharmonických složek signálů, které jsou sice identifikovatelné ve vyšších pásmech, ale v základním pásmu 0X až ±1X je jejich identifikace ztížena nebo dokonce znemožněna přítomností provozního šumu.
Jak je zřejmé ze znázornění spektrogramu akumulovaného úplného spektra, který je pro měřicí rovinu 330 znázorněn na obr. 4, může během jeho výpočtu dojít k téměř úplné ztrátě informace o přítomnosti subharmonických složek signálů snímačů 3, 30, 4, 40 rotorových vibrací na vyšších frekvencích.
Na obr. 5 je pak vykreslen průběh úplného spektra vypočteného ze signálů snímačů 3, 30 nebo 4, 40 rotorových vibrací, z jehož průběhu je patrné, že šum dosahuje v intervalu -3X až 3X téměř totožné úrovně jako amplitudy násobků subharmonických složek signálů, takže v tomto pásmu nejsou vidět složky s frekvencí 1/2X, 3/2X a 5/2X a viditelné jsou až složky s frekvencí od 7/2X, resp. -7/2X, které již nejsou šumem překryty.
Pro snížení úrovně šumu v jednotlivých frekvenčních pásmech se úplné spektrum dále normuje. Přitom se výpočtem mediánu šumu (znázorněn přerušovanou čárou na obr. 5) pro jednotlivá frekvenční pásma úplného spektra odhadne úroveň šumu. Tento medián se následně od amplitud v jednotlivých frekvenčních pásmech úplného spektra odečte. Z porovnání úplného spektra (světlejší dolní křivka na obr. 6) a normovaného úplného spektra, tj. úplného spektra s potlačeným šumem (horní křivka na obr. 6), je zřejmé, že úplné spektrum se normováním „narovnalo“ a jeho subharmonické složky více vynikly.
U způsobu detekce částečného kontaktu rotor-stator při provozu turbíny podle vynálezu se tedy nejprve úplné spektrum signálů snímačů 3, 30, 4, 40 rotorových vibrací v každém časovém okamžiku akumuluje, resp. průměruje ve frekvenční oblasti, a zároveň se normuje, resp. se od něj odečítá odhadnutá úroveň šumu v jednotlivých frekvenčních pásmech, čímž se vytvoří normované akumulované úplné spektrum signálů snímačů 3, 30, 4, 40 rotorových vibrací, jehož příkladný spektrogram pro měřicí rovinu 330 a pro snímače 3 a 30 rotorových vibrací je znázorněn na obr.
7. Z jeho průběhu je patrné, že při provozu turbíny docházelo k částečnému kontaktu rotor-stator - viz subharmonické složky s frekvencí ±1/2X, ±1/3X a ±2/3 X přítomné ve spektrogramu.
Spektrogram normovaného akumulovaného úplného spektra se přitom v reálném čase vykresluje na monitor či displej, případně tiskne, přičemž obsluha zařízení sleduje amplitudu určité subharmonické složky/složek signálů rotorových vibrací (např. 1/9X až 1/2X). Pokud amplituda této subharmonické složky/složek naroste, obsluha usuzuje na existenci kontaktu rotor-stator. V takovém případě obsluha provede kroky k potlačení tohoto kontaktu, nebo odstaví turbínu.
-6CZ 306833 B6
Tento postup diagnostiky lze dále automatizovat sledováním průběhu charakteristické veličiny subharmonických složek signálů snímačů 3, 30, 4, 40 rotorových vibrací s frekvencemi rovnými celočíselným podílům otáčkové frekvence rotoru 1 definované následovně:
= max
Acaf(t> Sf)
kde Acaf(t, ±Sf) označuje amplitudu normovaného akumulovaného úplného spektra pro subharmonickou frekvenci sf v čase t a Acaf (t, ±/fe) označuje amplitudu frekvenční linie v blízkém okolí této frekvence v čase t, M je celkový počet frekvenčních linií, které se použijí v průměrování, k je index součtu a fk je frekvenční linie, která se nachází v okolí subharmonické frekvence Sf.
Charakteristická veličina se přitom srovnává s předem zvolenou prahovou hodnotou, jejíž překročení je známkou přítomnosti částečného kontaktu rotor-stator. Tato prahová hodnota se přitom vždy volí vyšší než 1 (typicky 2) a udává tak faktor, s jakým překračuje amplituda subharmonické složky signálů hladinu amplitud frekvenčních linií ve svém blízkém okolí. Pokud ke kontaktu rotor-stator nedochází a tedy subharmonická složka signálů rotorových vibrací a ani její násobky nejsou v signálech obsaženy, nabývá hodnota charakteristické veličiny pro danou subharmonickou složku hodnoty blízké 1. Přítomnost konkrétní subharmonické složky v signálech pak způsobbí, že hodnota charakteristické veličiny bude hodnotu 1 a předem zvolenou prahovou hodnotu převyšovat.
Z průběhu charakteristické veličiny subharmonických složek signálu s frekvencemi rovnými celočíselným podílům otáčkové frekvence rotoru 1, na obr. 8, ve kterém je znázorněn průběh charakteristické veličiny pro frekvence 1/2X a 1/3X, je zřejmá jejich velká fluktuace. V případě, že ke kontaktu rotor-stator nedochází, lze charakteristickou veličinu chápat jako veličinu náhodnou popsanou svojí střední hodnotou a směrodatnou odchylkou. Průběh charakteristické veličiny pak lze aproximovat součtem její střední hodnoty a směrodatné odchylky, čímž se získá její homí obálka ve formě:
= M(t) + ^(t), kde p(t) je střední hodnota charakteristické veličiny a o(t) je směrodatná odchylka charakteristické veličiny.
Střední hodnota a směrodatná odchylka však musejí být vzhledem k nestacionaritě charakteristické veličiny, určeny rekurzivním algoritmem:
M(/c) = Z./z(/c-l) + (l-A).es/(/c) a2(k) = λ.σ2(/ί - 1) + (1 -A).[fs/(k) -μ(/τ)]2 , kde λ představuje faktor zapomínání a s ohledem na vzorkovací frekvenci signálu a zvolenou časovou konstantu t se volí nej častěji v intervalu 0,9 až 0,9999, k je krok a s každým vzorkem signálu se tedy zvyšuje o 1.
Faktor zapomínání λ je možné vyjádřit vztahem kde fs je vzorkovací frekvence signálu a t je zvolená časová konstanta.
Tento algoritmus přitom současně slouží jako filtr, neboť v případě, kdy dochází ke kontaktu rotor-stator, se mění střední hodnota i směrodatná odchylka, což je patrné z obr. 8, resp. jeho spodní části, kde například v časovém rozmezí 1400 až 1500 s došlo ke kontaktu rotor-stator.
Po pozitivní detekci kontaktu rotor-stator následuje jeho lokalizace, resp. určení místa/míst na stroji, kde k tomuto kontaktu došlo, resp. dochází. Pro axiální lokalizaci tohoto místa/míst se použijí signály neznázorněných snímačů absolutních vibrací ložisek 2, 20 turbíny (které jsou standardní součástí polní instrumentace turbín), přičemž úseky těchto signálů v okamžiku vzniku, resp. průběhu kontaktu rotor-stator, se zpracují v časo-frekvenční oblasti. Přitom se identifikují širokopásmová amplitudová navýšení (rázy) ve spektru těchto signálů a následně se podle jejich charakteru určí frekvenční pásmo, ve kterém jsou tyto rázy dostatečně zřejmé, a v tomto pásmu se dále zpracovávají (např. libovolné frekvenční pásmo mezi frekvencemi 0 až 150 Hz či 250 až 300 Hz ve spektrogramu na obr. 9). Pro zpracování signálů je možné použít libovolné časofrekvenční zpracování, např. krátkodobou Fourierovu transformaci, avšak vzhledem kjejich nestacionaritě je výhodnější použít vínkovou transformaci, při níž se jako bázová funkce použije Morletova vlnka.
Po identifikaci širokopásmových amplitudových změn (rázů) následuje určení času příchodu těchto rázů k jednotlivým snímačům absolutních vibrací ložisek 2, 20 turbíny - tzv. paty rázu, což je určení okamžiků nárůstu amplitud signálu nad úroveň šumu. Znalost těchto pat rázů a jejich vzájemných časových diferencí pak umožňuje v kombinaci se znalostí rozmístění snímačů absolutních vibrací ložisek 2, 20 v turbíně a rychlosti šíření rázu v materiálu rotoru/statoru axiální lokalizaci místa/míst kontaktu rotor-stator. Při této lokalizaci se získá několik desítek až stovek výsledků ukazujících místa vzniku kontaktu v turbíně, přičemž tyto výsledky se zobrazí pomocí histogramu absolutních četností, u kterého je na ose x zavedena axiální vzdálenost snímačů absolutních vibrací ložisek 2, 20 v turbíně. Po jeho spojení s výkresem dané turbíny nebo její osou je pak zřejmá axiální poloha míst, ve kterých došlo ke kontaktu rotor-stator po délce rotoru 1 turbíny, přičemž jsou jasně rozlišena místa s častým kontaktem od míst s pouze občasným případně falešným kontaktem - viz obr. 10.
Pro upřesnění lokalizace místa/míst, kde dochází ke kontaktu rotor-stator, lze dále ze získaných dat určit i jeho/jejich obvodovou polohu. Přitom se s využitím synchronně vzorkovaných signálů snímačů 3, 30 a 4, 40 rotorových vibrací a signálů neznázorněných snímačů absolutních vibrací ložisek, které jsou standardní součástí polní instrumentace turbín, identifikuje úhel natočení rotoru 1, kdy dochází ke kontaktu rotor-stator. Amplitudové změny (rázy) způsobené kontaktem rotor-stator se přitom zobrazí do orbity 1. harmonické složky otáčkové frekvence rotoru 1, čímž se určí úhlová poloha maxima těchto amplitudových změn, která současně představuje obvodovou polohu místa kontaktu rotor-stator - viz úsečka a na obr. 11.
Claims (8)
1. Způsob detekce a lokalizace částečného kontaktu rotor-stator při provozu turbíny, u kterého se rotor (1) turbíny snímá alespoň jedním párem snímačů (3, 30, 4, 40) rotorových vibrací umístěným v jedné měřicí rovině (330, 440) kolmé kose rotoru (1) a alespoň jedním párem snímačů absolutních vibrací ložisek (2, 20), vyznačující se tím, žeu synchronně digitalizovaných signálů snímačů (3, 30, 4, 40) rotorových vibrací se vypočte aritmetický průměr jejich harmonických pásem vymezených otáčkovou frekvencí X rotoru (1) turbíny a jejími vyššími harmonickými násobky v úplném spektru ve frekvenční oblasti, čímž se převedou informace o subharmonických složkách a jejich násobcích v úplném spektru těchto signálů do frekvenčního pásma -IX až IX (kde X označuje otáčkovou frekvenci rotoru (1)), zároveň se výpočtem mediánu odhadne úroveň šumu, kterým jsou zatížena jednotlivá frekvenční pásma
-8CZ 306833 B6 úplného spektra, a která se od amplitud v jednotlivých frekvenčních pásmech úplného spektra signálů snímačů rotorových vibrací odečte, přičemž se v takto vytvořeném normovaném akumulovaném úplném spektru signálů snímačů rotorových vibrací porovnává amplituda předem stanovené subharmonické složky/složek signálů rotorových vibrací s amplitudami frekvenčních složek v jejím/jejich okolí stanovenými Fourierovou transformací nebo aproximací Fourierovou řadou, přičemž se vizuálním vyhodnocením ze strany obsluhy a/nebo automatickým vyhodnocením určuje přítomnost sledované subharmonické složky signálů rotorových vibrací ve spektru těchto signálů, která je znakem existence kontaktu rotor-stator, a v případě jeho existence se dále lokalizuje místo/místa turbíny, kde k tomuto kontaktu dochází, přičemž při axiální lokalizaci se zpracují úseky signálů snímačů absolutních vibrací ložisek (2, 20) turbíny v okamžiku vzniku, resp. průběhu kontaktu rotor-stator, v časo-frekvenční oblasti, ve které se identifikují širokopásmová amplitudová navýšení těchto signálů a následně se určí časy nárůstu amplitud signálů nad úroveň šumu, a ze znalosti těchto časů a jejich časových diferencí, rozmístění snímačů absolutních vibrací ložisek (2, 20), rychlosti šíření rázu v materiálu rotoru/statoru se axiální lokalizací určí axiální poloha místa/míst kontaktu rotor-stator, a srovnáním histogramu absolutních četností těchto míst s výkresem nebo osou turbíny se určí axiální poloha místa/míst, kde ke kontaktu rotor-stator skutečně dochází, a pro obvodovou lokalizaci se amplitudová navýšení signálů snímačů absolutních vibrací ložisek zobrazí do orbity 1. harmonické složky otáčkové frekvence rotoru (1) vytvořené ze signálů snímačů (3, 30, 4, 40) rotorových vibrací, čímž se určí úhlová poloha maxima amplitudových navýšení způsobených kontaktem rotor-stator, která odpovídá obvodové poloze místa/míst, kde dochází ke kontaktu rotor-stator.
2. Způsob podle nároku 1, vyznačující se tím, že snímače (3,30, 4, 40) rotorových vibrací jsou snímače relativních rotorových vibrací.
3. Způsob podle nároku 1 nebo 2, vyznačující se tím, že snímače (3, 30, 4, 40) rotorových vibrací jsou snímače absolutních rotorových vibrací.
4. Způsob podle libovolného z předcházejících nároků, vyznačující se tím, že při vizuálním vyhodnocování sleduje obsluha amplitudu předem stanovené subharmonické složky/složek signálů rotorových vibrací v intervalu 1/9X až 1/2X, kterou porovnává s amplitudami frekvenčních složek v jejím/jejich okolí, přičemž pokud amplituda této subharmonické složky/složek převýší amplitudy okolních frekvenčních složek, znamená to, že sledovaná subharmonická složka/složky je/jsou přítomná/přítomné ve spektru signálů snímačů rotorových vibrací, a že dochází ke kontaktu rotor-stator.
5. Způsob podle libovolného z nároků 1 až 3, vyznačující se t í m , že při automatickém vyhodnocování se stanoví charakteristická veličina subharmonických složek s frekvencemi rovnými celočíselným podílům otáčkové frekvence rotoru 1 definovaná vztahem = w-ax kde Acaf(t, ±sý označuje amplitudu normovaného akumulovaného spektra pro subharmonickou frekvenci Sf v čase t a Acaj ít.fjJ označuje amplitudu frekvenční linie v blízkém okolí této frekvence, M je celkový počet frekvenčních linií, které se použijí k průměrování, k ']C index součtu a fk je hodnota frekvenční linie, která se nachází v okolí subharmonické frekvence Sf, přičemž pokud ke kontaktu rotor-stator nedochází a tedy subharmonická složka ani její násobky nejsou v signálech rotorových vibrací obsaženy, nabývá hodnota charakteristické veličiny pro danou subharmonickou složku hodnoty blízké 1, a pokud ke kontaktu rotor-stator dochází a tedy subharmonická složka signálu a/nebo její násobky jsou v signálech rotorových vibrací obsaženy,
-9CZ 306833 B6 nabývá hodnota charakteristické veličiny pro danou subharmonickou složku hodnotu převyšující 1.
6. Způsob podle nároku 5, vyznačující se tím, že průběh charakteristické veličiny se aproximuje součtem střední hodnoty μ a směrodatné odchylky σ vypočtených rekurzivním algoritmem ve formě:
μΜ = λ. μ(]ί -1) + (1- Λ). fs/(k) σ2(/0 = λ. a\k -1) + (1- λ). (fc) - μΜ]2 , kde λ představuje faktor zapomínání a s ohledem na vzorkovací frekvenci signálu a zvolenou časovou konstantu t se volí nejčastěji v intervalu 0,9 až 0,9999, k je krok a s každým vzorkem signálu se tedy zvyšuje o 1, přičemž ze změny střední hodnoty μ a směrodatné odchylky σ se usuzuje na existenci kontaktu rotor-stator.
7. Způsob podle libovolného z předcházejících nároků, vyznačující se tím, že při axiální lokalizaci se zpracují úseky signálů snímačů absolutních vibrací ložisek (2, 20) turbíny v okamžiku vzniku, resp. průběhu kontaktu rotor-stator, vínkovou transformací.
8. Způsob podle nároku 7, vyznačující se tím, že jako bázová funkce vínkové transformace se použije Morletova vlnka.
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CZ2012-947A CZ306833B6 (cs) | 2012-12-20 | 2012-12-20 | Způsob detekce a lokalizace částečného kontaktu rotor-stator při provozu turbíny |
| PL13196683T PL2746541T3 (pl) | 2012-12-20 | 2013-12-11 | Sposób wykrywania i lokalizacji częściowego ścierania w układzie wirnik-stojan podczas eksploatacji turbiny |
| US14/103,514 US9903787B2 (en) | 2012-12-20 | 2013-12-11 | Method of detecting and localizing partial rotor-stator rubbing during the operation of a turbine |
| EP13196683.0A EP2746541B1 (en) | 2012-12-20 | 2013-12-11 | A method of detecting and localizing partial rotor-stator rubbing during the operation of a turbine |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CZ2012-947A CZ306833B6 (cs) | 2012-12-20 | 2012-12-20 | Způsob detekce a lokalizace částečného kontaktu rotor-stator při provozu turbíny |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CZ2012947A3 CZ2012947A3 (cs) | 2014-07-02 |
| CZ306833B6 true CZ306833B6 (cs) | 2017-08-02 |
Family
ID=49911149
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CZ2012-947A CZ306833B6 (cs) | 2012-12-20 | 2012-12-20 | Způsob detekce a lokalizace částečného kontaktu rotor-stator při provozu turbíny |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US9903787B2 (cs) |
| EP (1) | EP2746541B1 (cs) |
| CZ (1) | CZ306833B6 (cs) |
| PL (1) | PL2746541T3 (cs) |
Families Citing this family (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6251658B2 (ja) * | 2014-09-12 | 2017-12-20 | 株式会社神戸製鋼所 | 回転機の異常検知装置、該方法および該システムならびに回転機 |
| CN104634526B (zh) * | 2015-01-27 | 2017-06-30 | 西安交通大学 | 一种基于非线性压缩变换的转子碰摩故障检测方法及系统 |
| JP6533141B2 (ja) * | 2015-10-08 | 2019-06-19 | 株式会社メガチップス | 振動解析システム、ユーザーインターフェースシステムおよび振動解析データ生成方法 |
| US10263553B2 (en) * | 2015-12-14 | 2019-04-16 | Rolls-Royce North American Technologies Inc. | Synchronous electrical power distribution system |
| CN106248379B (zh) * | 2016-08-19 | 2018-08-21 | 北京航空航天大学 | 一种固体润滑轴承加速寿命试验载荷谱设计方法 |
| CN106285799B (zh) * | 2016-08-23 | 2017-08-08 | 中国大唐集团科学技术研究院有限公司华东分公司 | 一种基于谐分量分解的转子振动保护方法 |
| DE102017200964A1 (de) * | 2017-01-20 | 2018-07-26 | Rolls-Royce Deutschland Ltd & Co Kg | Messvorrichtung und Messverfahren zur Erfassung von Mischreibungsereignissen und / oder Stick-Slip-Ereignissen |
| US10801359B2 (en) * | 2017-03-14 | 2020-10-13 | General Electric Company | Method and system for identifying rub events |
| RU2668358C1 (ru) * | 2017-11-17 | 2018-09-28 | Публичное акционерное общество "ОДК-Уфимское моторостроительное производственное объединение" (ПАО "ОДК-УМПО") | Способ обнаружения резонансных колебаний ротора газотурбинного двигателя |
| CN108952837B (zh) * | 2018-06-21 | 2021-02-26 | 东南大学 | 监测汽轮机启动至热态过程中猫爪负荷变化的装置及方法 |
| CN108917916A (zh) * | 2018-07-25 | 2018-11-30 | 北京振测智控科技有限公司 | 一种转轴绝对振动的测量方法 |
| CN109141893B (zh) * | 2018-09-29 | 2020-03-13 | 天津大学 | 气动双转子-机匣多功能试验台 |
| DE102019212369A1 (de) * | 2019-08-19 | 2021-01-21 | Zf Friedrichshafen Ag | Verifikationsverfahren für Sensordaten |
| CN112746875B (zh) * | 2019-10-31 | 2022-08-19 | 中国航发商用航空发动机有限责任公司 | 航空发动机转子轴系复杂振动的主动控制系统、方法 |
| US11639670B2 (en) | 2019-11-14 | 2023-05-02 | General Electric Company | Core rub diagnostics in engine fleet |
| CN111780857B (zh) * | 2020-06-05 | 2022-02-15 | 南京曦光信息科技有限公司 | 一种基于谐波累加的p-otdr系统的多点扰动定位检测方法 |
| CN114070158B (zh) * | 2020-08-04 | 2023-12-29 | 美的威灵电机技术（上海）有限公司 | 基于转矩信息的电机的控制方法、电机和存储介质 |
| US11698287B2 (en) * | 2020-08-31 | 2023-07-11 | Rolls-Royce Deutschland Ltd & Co Kg | System and method for detecting vibrations in rotating machinery |
| CN114061729A (zh) * | 2021-10-22 | 2022-02-18 | 中国大唐集团科学技术研究总院有限公司华东电力试验研究院 | 转子突变振动故障定位方法及系统 |
| DE102022113888A1 (de) | 2022-06-01 | 2023-12-07 | Rheinisch-Westfälische Technische Hochschule (RWTH) Aachen, Körperschaft des öffentlichen Rechts | Verfahren und System zur Überwachung eines Gleitlagers |
| EP4308797B1 (en) * | 2022-06-02 | 2024-02-28 | Západoceská Univerzita V Plzni | Method of automatic detection of synchronous rubbing in turbine |
| CN115015756A (zh) * | 2022-07-14 | 2022-09-06 | 深圳市文浩科技有限公司 | 一种细碎融合的新材料混合机故障诊断方法 |
| CN117053920B (zh) * | 2023-08-21 | 2024-09-03 | 哈尔滨东安汽车发动机制造有限公司 | 一种变速器摩擦片振动噪声故障测试装置及其测试方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4377947A (en) * | 1979-11-16 | 1983-03-29 | Hitachi, Ltd. | Device for identifying rubbing locations in rotating machines |
| US4478082A (en) * | 1981-08-26 | 1984-10-23 | Hitachi, Ltd. | Method and apparatus for detecting rubbing in a rotary machine |
| US20040120809A1 (en) * | 2002-12-20 | 2004-06-24 | Peter Loftus | Rotor system |
| WO2008060164A1 (en) * | 2006-11-17 | 2008-05-22 | Dynatrend As | A method and a device for detection of rubbing between a stator and a rotor of a turbo-machine |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS54149696A (en) * | 1978-05-16 | 1979-11-24 | Hitachi Ltd | Rubbing detection between rotary body and stationary body |
| US7409319B2 (en) | 2003-11-24 | 2008-08-05 | General Electric Company | Method and apparatus for detecting rub in a turbomachine |
| US20080240902A1 (en) | 2007-03-28 | 2008-10-02 | General Electric Company | Method and system for rub detection in a steam turbine |
| US8177474B2 (en) | 2007-06-26 | 2012-05-15 | General Electric Company | System and method for turbine engine clearance control with rub detection |
| US7866213B2 (en) * | 2008-06-18 | 2011-01-11 | Siemens Energy, Inc. | Method of analyzing non-synchronous vibrations using a dispersed array multi-probe machine |
-
2012
- 2012-12-20 CZ CZ2012-947A patent/CZ306833B6/cs unknown
-
2013
- 2013-12-11 PL PL13196683T patent/PL2746541T3/pl unknown
- 2013-12-11 US US14/103,514 patent/US9903787B2/en active Active
- 2013-12-11 EP EP13196683.0A patent/EP2746541B1/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4377947A (en) * | 1979-11-16 | 1983-03-29 | Hitachi, Ltd. | Device for identifying rubbing locations in rotating machines |
| US4478082A (en) * | 1981-08-26 | 1984-10-23 | Hitachi, Ltd. | Method and apparatus for detecting rubbing in a rotary machine |
| US20040120809A1 (en) * | 2002-12-20 | 2004-06-24 | Peter Loftus | Rotor system |
| WO2008060164A1 (en) * | 2006-11-17 | 2008-05-22 | Dynatrend As | A method and a device for detection of rubbing between a stator and a rotor of a turbo-machine |
Non-Patent Citations (1)
| Title |
|---|
| XP026422078 PATEL T. H., et al.: "Coupled bending-torsional vibration analysis of rotor with rub and crack", JOURNAL OF SOUND & VIBRATION, LONDON, GB, Vol. 326, No. 3-5, ISSN 0022-460X, (20091009), pages 740 - 752, 07.06.2009 * |
Also Published As
| Publication number | Publication date |
|---|---|
| EP2746541A2 (en) | 2014-06-25 |
| EP2746541B1 (en) | 2017-07-19 |
| US20140230555A1 (en) | 2014-08-21 |
| CZ2012947A3 (cs) | 2014-07-02 |
| EP2746541A3 (en) | 2016-04-20 |
| PL2746541T3 (pl) | 2018-05-30 |
| US9903787B2 (en) | 2018-02-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CZ306833B6 (cs) | Způsob detekce a lokalizace částečného kontaktu rotor-stator při provozu turbíny | |
| US10359473B2 (en) | Detecting faults in turbine generators | |
| Saruhan et al. | Vibration analysis of rolling element bearings defects | |
| Liang et al. | An energy operator approach to joint application of amplitude and frequency-demodulations for bearing fault detection | |
| US8994359B2 (en) | Fault detection based on current signature analysis for a generator | |
| Klausen et al. | Cross-correlation of whitened vibration signals for low-speed bearing diagnostics | |
| CN109883703B (zh) | 一种基于振动信号相干倒谱分析的风机轴承健康监测诊断方法 | |
| Klausen et al. | Multi-band identification for enhancing bearing fault detection in variable speed conditions | |
| JP2011259624A (ja) | 転がり軸受部振動データの高周波電磁振動成分除去方法および高周波電磁振動成分除去装置、回転機械の転がりの軸受診断方法および軸受診断装置 | |
| CN111397877B (zh) | 一种旋转机械拍振故障检测与诊断方法 | |
| CN109029999B (zh) | 基于增强调制双谱分析的滚动轴承故障诊断方法 | |
| Corne et al. | Comparing MCSA with vibration analysis in order to detect bearing faults—A case study | |
| KR101857355B1 (ko) | 웨이브렛과 순환대수 엔빌로프 스펙트럼을 이용한 회전기기의 고장을 검출하는 방법 | |
| Liang et al. | Intelligent bearing fault detection by enhanced energy operator | |
| Luo et al. | Effective and accurate approaches for wind turbine gearbox condition monitoring | |
| Kumar et al. | A combined approach for weak fault signature extraction of rolling element bearing using Hilbert envelop and zero frequency resonator | |
| Leite et al. | Bearing fault detection in induction machine using squared envelope analysis of stator current | |
| JP7098399B2 (ja) | 状態監視装置および状態監視方法 | |
| JP5958932B2 (ja) | 動的設備の状態監視システムとその方法とそのプログラム | |
| JP7394031B2 (ja) | 転がり軸受の異常検出装置、及び異常検出方法 | |
| Matic et al. | Vibration Based Broken Bar Detection in Induction Machine for Low Load Conditions. | |
| Zhu et al. | Fault diagnosis method for rolling element bearings under variable speed based on TKEO and Fast-SC | |
| Lim et al. | Improved blade fault diagnosis using discrete Blade Passing Energy Packet and rotor dynamics wavelet analysis | |
| RU2670771C9 (ru) | Способ определения характера касания лопатки вращающегося колеса о корпус турбомашины | |
| CN111594893B (zh) | 一种吸油烟机抖动故障诊断方法 |