Erste Aktualisierung zur Analyse von Vibrationsschäden in Zentrifugalkompressorenanlagen
Thread Content
Rotationsunausgeglichenheit – Der drehende Teil einer rotierenden Maschine wird als Rotor bezeichnet. Der Rotor weist aufgrund von Fertigungsfehlern, Montagefehlern, ungleichmäßiger Materialstruktur sowie während des Betriebs aufgrund von Korrosion, Abnutzung, Ablagerungen und abfallenden Bauteilen eine exzentrische Masse auf. Abbildung 2-1: Kraftmodell des rotierenden Rotors. Fig.2-1 Force model of rotating. Sei die Masse des Systems M, die exzentrische Masse des Rotors m und der Exzentrizitätsabstand e (siehe Abbildung 2-1 [14]). Wenn sich der Rotor mit einer bestimmten Winkelgeschwindigkeit ω dreht, entsteht durch die exzentrische Masse m eine Zentrifugalkraft F0 = meω². Die Bewegungsgleichungen lauten:Mx″ + Cx′ + Kx = meω²cos(ωt)
My″ + Cy′ + Ky = meω²sin(ωt) (2-1)
Die Amplituden X und Y sind gegeben durch:
X = Y = meω²/(K – Mω²) + (Cω)²
Die Lösungen für x und y sind:
x = X cos(ωt – Φ)
y = Y sin(ωt – Φ) (2-2)
Sei p = …, λ = ω/p, ξ = C/(2Mp). Dann gelten:
X = Y = … (2-3)
Der Phasenunterschied ist gegeben durch:
tanΦ = 2ξλ/(1 – λ²) (2-4)
Die Gleichungen (2-3) und (2-4) können durch Amplitudenzahl-Frequenzkurven sowie Phasenzahl-Frequenzkurven dargestellt werden (Abbildung 2-2 [14]).
Abbildung 2-2: Amplitudenzahl-Frequenz- und Phasenzahl-Frequenzkurven
Wenn sich der Rotor mit konstanter Winkelgeschwindigkeit ω dreht und dabei eine Biegung a aufweist, dann dreht sich der Rotor um seine Achse – gleichzeitig dreht sich auch die Achse um die Verbindungslinie zum Mittelpunkt der Lager. Dieses Phänomen wird als Präzession der Achse bezeichnet. Die Gleichungen für die Bewegung des Schwerpunkts des Rotors lauten laut [14]:
Mx” + Cx’ + Kx = Meω²cosωt
My” + Cy’ + Ky = Meω²sinωt (2-5)
Die Biegung des Mittelpunkts der Drehachse beträgt laut [14]: (2-6)
Aus diesen Gleichungen geht hervor, dass bei λ << 1 die Biegung a sehr klein ist und der Winkel Φ ebenfalls sehr gering ist – er nähert sich null an ; Wenn λ=1 ist und die Dämpfung sehr gering ist, kann die Biegung a sehr groß sein. In diesem Fall beträgt der Phasenwinkel Φ=90° ; Wenn λ >> 1 ist, nähert sich der Wert der Biegung a dem Wert e an. Der Winkel Φ nähert sich 180° an. Der Schwerpunkt G liegt dann in der Nähe des festen Punktes O – wie in Abbildung 2-3 [14] dargestellt. Abbildung 2-3: Phasenänderung des Schwerpunkts des Rotors. In Abbildung 2-3 ist die Bahn der Präzession der Rotorexzentrizität durch eine gedachte Kreislinie dargestellt. Die Drehrichtung der Präzession des Rotors kann entweder der Drehrichtung des Rotors selbst entsprechen oder ihr entgegengesetzt sein. Die Winkelgeschwindigkeit der Präzession kann außerdem der Drehgeschwindigkeit der Achse entsprechen oder auch nicht. Wie aus der Abbildung zu erkennen ist, folgt die Bahn des Rotorendrehpunkts einer Kreisbahn. Doch beim tatsächlichen Rotor unterscheiden sich aufgrund der unterschiedlichen Biegesteifigkeiten der Achse – insbesondere wegen der unterschiedlichen Stützsteifigkeiten – die Reaktionen des Rotors auf ungleichmäßige Massen. In den x- und y-Richtungen sind daher nicht nur die Amplituden unterschiedlich, sondern auch der Phasenunterschied beträgt nicht 90°. Daher ist die Bahn des Rotorendrehpunkts keine Kreisbahn, sondern eine Ellipse. Wie in Abbildung 2-4 [14] dargestellt: Abbildung 2-4: Bahn des Rotorkerns bei Unausgeglichenheit des Rotors. Fig.2-4: Die Bahn von A* bei Unausgeglichenheit des Rotors. Die Merkmale der Vibrationen aufgrund von Unausgeglichenheit sind hohe Spitzen im Frequenzspektrum bei der Nenndrehfrequenz; manchmal treten auch hohe Harmonischen wie die doppelte Frequenz auf. Bei konstanter Drehzahl bleibt die Phasenlage stabil, während die Amplitude proportional zur Unausgeglichenheit ist und mit der Änderung der Drehzahl zunimmt oder abnimmt – mit Ausnahme der Fall, wenn die Drehzahl die kritische Geschwindigkeit überschreitet. Die Bahn des Rotorkerns ist eine elliptische, synchron rotierende Bahn.
Fig.2-8: Radial displacement vibration of rotor
Wenn die Achsen der beiden Rotoren einen Winkelabstand zueinander aufweisen (wie in Abbildung 2-9 [18] dargestellt), dann ist die Drehzahl des Antriebsrotors laut [18] wie folgt gegeben:
ω2 = ω1·cosα / (1 – sin²α·cos²φ1) (2-13)
Das Übersetzungsverhältnis beträgt:
i21 = cosα / (1 – sin²α·cos²φ1) (2-14)
Dabei ist:
ω1 – Drehzahl des Antriebsrotors
α – Winkelabstand zwischen den Achsen der Rotoren
φ1 – Winkel des Antriebsrotors
Abbildung 2-9: Winkelverschiebungsvibration des Rotors
Fig.2-9: Declination displacement vibration of rotor
Wenn die Drehzahl des Antriebsrotors ω1 konstant ist, dann ist die Drehzahl des Antriebsrotors ω2 nicht konstant. Sie hängt vielmehr vom Winkelabstand α sowie vom Winkel φ1 des Antriebsrotors ab. Wenn φ1=0° oder 180° ist, ist i21= 1/cosα am größten. Wenn φ1=90° oder 270° ist, ist i21= cosα am kleinsten – wie in Abbildung 2-10 dargestellt [18]. Abbildung 2-10: Verlaufskurve des Übersetzungsverhältnisses i21. Aus Abbildung 2-10 geht hervor, dass sich das Übersetzungsverhältnis nicht nur bei jeder Umdrehung des Rotors um zwei Werte ändert, sondern dass die Amplitude dieser Änderungen mit zunehmendem Winkel α größer wird. Dadurch wird der Antriebsrotor durch die durch das ändernde Übersetzungsverhältnis entstehende Winkelbeschleunigung in Schwingung versetzt. Die frequenz der radialen Schwingungen beträgt außerdem das Doppelte der Drehfrequenz des Rotors. In der praktischen Produktion ist die Ausrichtung der Verbindungen zwischen den Rotoren einer Anlage in der Regel eine Kombination aus radialer Verschiebung und Winkelverschiebung. Daher ist das Phänomen der radialen Schwingungen der Rotoren das Ergebnis dieser beiden Effekte zusammen. Außerdem führen die zusätzlichen axialen Kräfte, die durch den Zahnkupplungsmechanismus entstehen, sowie die Neigung des Rotors dazu, dass der angetriebene Rotor in jedem Umlauf einmal axial hin- und herbewegt wird. Dadurch ist die Frequenz der axialen Schwingungen des Rotors identisch mit der Drehfrequenz. 2. Schwingungsmechanismus bei Fehlstellung von starren und flexiblen Kupplungen Bei der Verbindung von Rotorblättern mit einer starren Kupplung, die nicht in der richtigen Position sind, biegen sich die Rotorblätter aufgrund der erzwungenen Verbindung. Die Kontaktflächen der Kupplung bewegen sich dabei bei jeder Umdrehung des Rotors um zwei Umdrehungen. Daher ist ihre Schwingungsfrequenz doppelt so hoch wie die Drehfrequenz in Betrieb. Die axiale Schwingungsfrequenz ist mit der Betriebsfrequenz synchronisiert. Das Schwingungsverhalten von flexiblen Kupplungen ist ähnlich wie das von starren Kupplungen. Allerdings wird bei jeder Umdrehung jeder Bolzen der Kupplung einmal gestreckt und wieder zusammengedrückt. Wenn die Anzahl der Bolzen eine gerade Zahl N ist, dann ist die Schwingfrequenz N-fach so hoch. Wenn hingegen die Anzahl der Bolzen eine ungerade Zahl N ist, dann ist die Schwingfrequenz 2N-fach so hoch [1]. In der Praxis werden die durch die Biegeverformungen entstehenden Kräfte größtenteils von den flexiblen Verbindungsbauteilen (z. B. elastischen Membranen) aufgenommen. Dadurch sind die entstehenden Schwingungen relativ gering. Zudem sind flexible Kupplungen weniger anfällig für Schwingungen aufgrund von Ausrichtungsfehlern – daher ist die Intensität der entstehenden Schwingungen ebenfalls relativ gering. Abbildung 2-11: Verbindung mit einem starren Kupplungselement bei Fehlausrichtung des Rotors [14] Abbildung 2-11: Verbindung eines starren Kupplungselements bei Fehlausrichtung des Rotors. Die Merkmale der Schwingungen bei Fehlausrichtung des Rotors sind folgende [19]: Im Spektrum bestehen die radialen Schwingungen aus der Grundfrequenz, der doppelten Frequenz sowie den modulierten Harmonischen. Die Amplitude der harmonischen Wellen mit der doppelten Frequenz ist besonders groß und stellt somit die charakteristische Frequenz dar. Das axiale Schwingungsspektrum setzt sich aus der Grundfrequenz und ihren Harmonischen zusammen, wobei die Grundfrequenz einen Peak aufweist. Bei konstanter Drehzahl bleibt die Phasenlage unverändert. In der Nähe des Kupplungspunkts sind die Vibrationen der Lager groß. Die Bahn des Achsenzentrums besteht aus einem großen Kreis um einen kleinen Kreis. Die Präzession erfolgt synchron in die gleiche Richtung. Die Lager sind empfindlich gegenüber Veränderungen der Last – die Vibrationen nehmen mit zunehmender Last zu.