S-Shaped Bifurcation Cascades and Loop Proliferation in an Active Magnetic Bearing Rotor with Cubic Shaft Stiffness and Nonlinear Aerodynamic Drag

Technische Samenvatting: S-vormige Bifurcatiecascades en Loop-proliferatie in een Actieve Magnetische Lager Rotor

Probleemstelling
Actieve magnetische lagers (AMB's) zijn cruciaal voor snel draaiende roterende machines en bieden contactloze ophanging via proportioneel-derivatie (PD) feedbackcontrole. De elektromagnetische kracht in AMB's is echter inherent nietlineair (omgekeerd evenredig met het kwadraat van de luchtspleet), wat complexe dynamische gedragingen kan induceren zoals meerwaardige frequentieresponsies, amplitude-sprongen en hysteresis. Bestaande literatuur heeft de elektromagnetische nietlineariteit, kubische asstijfheid en nietlineaire aerodynamische weerstand grotendeels afzonderlijk of in eenvoudige combinaties onderzocht. Deze studie adresseert de kloof in het begrip van hoe deze concurrerende nietlineaire mechanismen — specifiek de verzachtende aard van magnetische krachten versus de verhardende aard van kubische asstijfheid, gecombineerd met amplitude-afhankelijke aerodynamische weerstand — interageren binnen een enkel flexibel rotormodel. Het primaire doel is om de resulterende bifurcatietopologie te karakteriseren en de condities te identificeren waaronder complexe fenomenen, zoals gesloten lus-takken met talrijke saddle-node bifurcaties, ontstaan.

Methodologie
De auteurs modelleren een flexibele rotor ondersteund door AMB's als een vier vrijheidsgraden systeem, waarbij rekening wordt gehouden met:

• Nietlineaire Elektromagnetische Krachten: Gereguleerd door PD-feedbackcontrole, waarbij de stroom een functie is van de journal-verplaatsing en snelheid.
• Kubische Asstijfheid: Representatief voor geometrische nietlineariteit bij matige tot grote deflecties.
• Nietlineaire Aerodynamische Weerstand: Gemodelleerd als een combinatie van lineaire viskeuze demping (lage snelheid/kleine amplitude) en kwadratische turbulente demping (hoge snelheid/grote amplitude).

De leidende bewegingsvergelijkingen worden afgeleid met behulp van de Lagrange-formalisme en genon-dimensionaliseerd. Om het systeem op te lossen, gebruikt de studie de Harmonische Kwadratuur Methode (HQM), een techniek voor tijdsdomein-discretisatie die Galerkin-integralen vermijdt terwijl een hoge orde nauwkeurigheid behouden blijft. Deze methode wordt gecombineerd met pseudo-arc-length continuatie om zowel stabiele als onstabiele periodieke oplossings-takken te volgen over parameterruimten, en Floquet-stabiliteitsanalyse om stabiliteit te verifiëren en saddle-node bifurcaties te detecteren. Het numerieke kader wordt gevalideerd via convergentiestudies met betrekking tot harmonische orde en arc-length stapgrootte.

Belangrijkste Bijdragen en Resultaten
De studie onderzoekt systematisch de invloed van vier sleutelparameters: proportionele winst ($P$), derivatieve winst ($D$), nietlineaire stijfheidscoëfficiënt ($\eta_{NL}$) en aerodynamische dempingscoëfficiënten.

1. Proportionele Winst ($P$) en S-vormige Cascades:
   Het verhogen van $P$ drijft een progressieve bifurcatiecascade aan door vier onderscheidende regimes:

• Lineair Regime ($P < 2,68$): Enkelwaardige, stabiele respons.
• 1S-vormig Regime ($2,69 \le P \le 2,75$): Emergentie van een enkele S-vormige vouw met twee saddle-node bifurcaties.
• 2S-vormig Regime ($2,76 \le P \le 3,0$): Twee S-vormige takken, wat leidt tot toenemende multi-stabiliteit.
• 3S-vormig Regime ($3,1 \le P \le 4,1$): Drie S-vormige takken met acht vouwpunten.
  Naarmate $P$ stijgt van 2 naar 4, stijgt de primaire resonantie-amplitude van 0,2328 naar 0,3839, en de multi-stabiele frequentieband wordt aanzienlijk breder.

2. Derivatieve Winst ($D$) en Inverse Cascade:
   De derivatieve winst fungeert als een stabiliserend tegenmaatregel. Het verhogen van $D$ produceert een exacte inverse cascade aan die van $P$:

• Bij lage $D$ ($0,5 \le D \le 1,16$), blijft het systeem in het 3S-vormige regime.
• Bij intermediaire $D$ ($1,16 \le D \le 1,23$), transiteert het systeem naar een 1S-vormig regime.
• Voor $D \ge 1,23$, wordt de multi-stabiliteit vernietigd en wordt het quasi-lineaire gedrag hersteld.
  Dit onthult een symmetrische, complementaire relatie waarbij $P$ complexiteit induceert en $D$ deze systematisch onderdrukt.

3. Nietlineaire Stijfheid ($\eta_{NL}$) en Loop-proliferatie:
   Wanneer het systeem vaststaat in het 3S-vormige regime en $\eta_{NL}$ wordt verhoogd, concurreert een verhardend mechanisch nietlineair effect met de verzachtende magnetische nietlineariteit. Bij een kritische drempelwaarde ($\eta_{NL} \approx 0,55$) genereert het systeem een verbonden gesloten lus-tak die bovenop de hoofdresponscurve ligt.

• Deze lus bevat meer dan 36 door Floquet geverifieerde saddle-node bifurcatiepunten verspreid over het frequentie-interval $\Omega \in [0,2466, 0,2922]$.
• Een krachtvergelijking bevestigt dat dit een fysiek fenomeen is: de ratio van de kubische naar de magnetische kracht stijgt van 0,002% naar 8,6% bij het begin van de lus, wat aangeeft dat de verhardende stijfheid sterk genoeg wordt om de verzachtende magnetische respons te vervormen.

4. Regimeovergang van Aerodynamische Demping:
   De studie identificeert een amplitude-afhankelijke transitie in de dominantie van demping:

• Nietlineaire turbulente demping reguleert de primaire resonantie, waardoor grote-amplitude S-vormige lussen effectief worden onderdrukt.
• Lineaire viskeuze demping controleert superharmonische pieken, waar de amplitudes kleiner zijn.
  Een crossover vindt plaats nabij $\Omega \approx 0,27$. Modellen die slechts één dempingsmechanisme gebruiken, zouden de responsamplitudes systematisch foutief voorspellen (bijv. overschat een model met enkel lineaire demping de primaire resonantie-amplitude met ~57% vergeleken met een nietlineair model).

Betekenis en Claims
Het artikel claimt de eerste uitgebreide karakterisering te bieden van AMB-rotor dynamica onder de gelijktijdige werking van verzachtende elektromagnetische krachten, verhardende kubische stijfheid en amplitude-afhankelijke weerstand. De auteurs stellen dat:

• De geïdentificeerde S-vormige bifuratiecascades een kwantitatieve ontwerptabel bieden voor het afstemmen van PD-controllers, waarbij specifieke drempels (bijv. $P < 2,68$) worden gedefieerd om enkelwaardige operatie te garanderen.
• De inverse cascade door $D$ een kwantitatieve stabiliteitsgrens vaststelt ($D \ge 1,23$) voor het elimineren van multi-stabiliteit, wat een praktische controlestrategie biedt voor systemen die in hoog-gain regimes werken.
• De ontdekking van de verbonden gesloten lus-tak met >36 vouwpunten een complexe topologische transformatie vertegenwoordigt die verschillend is van simpelere S-vormige vouwen, wat de onvoorspelbaarheid van de rotorrespons in specifieke parameterrange door geschiedenis-afhankelijkheid benadrukt.
• De duale dempingsformulering wordt gepresenteerd als een modelleringsnoodzaak voor snelle AMB-rotoren om zowel primaire als superharmonische gedragingen accuraat te vangen.

De studie concludeert dat deze bevindingen de nietlineaire bifuratietopologie van AMB-systemen karakteriseren, wat directe implicaties heeft voor de veilige afstemming van PD-controllers en de voorspelling van rotordynamica in aanwezigheid van gecombineerde nietlineariteiten. Toekomstig werk wordt gesuggereerd om passieve stabilisatie (squeeze film dampers) en actieve controlestrategieën te verkennen om de geïdentificeerde loop-proliferatie te mitigeren.