⚡ electrical engineering

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

This paper utilizes the Harmonic Quadrature Method and Floquet stability analysis to reveal how proportional and derivative gains, cubic shaft stiffness, and nonlinear aerodynamic drag interact to drive complex S-shaped bifurcation cascades, extensive loop proliferation with over 36 saddle-node bifurcations, and distinct damping regimes in active magnetic bearing-supported flexible rotors.

Original authors: Farouk Thaljaoui, Thabet Guesmi, Thaljaoui Wathek

Published 2026-08-05

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Farouk Thaljaoui, Thabet Guesmi, Thaljaoui Wathek

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a world where heavy, spinning machines don't touch the ground at all. Instead of rolling on wheels or sliding on oil, they float in mid-air, suspended by invisible magnetic hands. This is the magic of Active Magnetic Bearings (AMBs), a technology used in everything from high-speed jet engines to energy storage flywheels. But here's the catch: magnets are tricky. The force they pull with doesn't grow in a straight line; it gets crazy strong the closer things get, like a rubber band that snaps back harder the more you stretch it. When you combine this magnetic "snap" with the natural wobble of a spinning shaft and the air resistance of a fast-moving object, the system can start behaving like a chaotic dance partner. It might suddenly jump to a huge, dangerous shake, or get stuck in a loop where it doesn't know which way to go. Engineers need to understand these jumps and loops to keep their machines from shaking themselves apart.

This paper dives deep into that chaotic dance, using a powerful computer simulation to map out exactly how a floating rotor behaves when you tweak its controls. The researchers treated the rotor like a flexible stick with a heavy disk in the middle, held up by magnetic bearings that use a "PD controller" (a fancy way of saying a smart brain that adjusts the magnet strength based on position and speed). They added in two extra ingredients: the shaft getting stiffer as it bends (like a spring that gets harder to push the more you push it) and air resistance that changes depending on how fast the rotor is shaking. By running thousands of simulations, they discovered that the "brain" of the machine has two knobs: a Proportional knob (P) that controls how hard the magnets pull, and a Derivative knob (D) that controls how much they resist motion.

The team found that turning up the Proportional knob (P) makes the system wilder and wilder. At low settings, the rotor behaves nicely, like a calm pendulum. But as they cranked P up, the smooth curve of the rotor's movement started to fold over itself, creating "S-shapes." It's like a rollercoaster track that suddenly loops back on itself. At first, there was one loop (1S-shaped), then two (2S-shaped), and finally three (3S-shaped). In these S-shaped zones, the rotor has multiple choices for how to vibrate at the same speed. It could be vibrating gently or shaking violently, and a tiny nudge could make it jump from safe to dangerous instantly. The researchers calculated that as they increased P from 2 to 4, the maximum shaking amplitude grew from 0.2328 to 0.3839, and the number of these "jump points" (called saddle-node bifurcations) multiplied, creating a wide zone of confusion where the machine's behavior depends entirely on its history.

However, the Derivative knob (D) acts like a superhero opposite. While P creates the chaos, D kills it. The paper shows that increasing D acts like pouring water on a fire, smoothing out those wild S-shaped loops. When D is low, the system is a mess of loops. But as the researchers turned D up, the loops shrank. At a specific setting of D = 1.23, the chaos vanished completely, and the system returned to a calm, predictable state. It's as if the P knob builds a maze, and the D knob is the key that unlocks the door to a straight path.

The most mind-bending discovery happened when they turned on the "stiffening" of the shaft (the cubic stiffness). When they combined the magnetic "softening" with the shaft's "hardening," something bizarre happened: the S-shaped loops didn't just fold; they proliferated into a giant, connected closed loop. Imagine a tangled ball of yarn that somehow connects back to itself. In this specific zone, the researchers found a single branch of the solution that twisted and turned with more than 36 different "fold points" (places where the path doubles back). This happened when the stiffness of the shaft grew strong enough to compete with the magnets, reaching a ratio where the shaft's force was about 8.6% of the magnetic force. In this zone, the rotor has dozens of different stable ways to vibrate at the same speed, making its behavior incredibly hard to predict without a perfect map.

Finally, the team looked at air resistance. They found that the air fights the rotor differently depending on how hard it's shaking. When the rotor is shaking violently (the main resonance), the air acts like thick, sticky mud (nonlinear damping), which is great at stopping the big jumps. But when the rotor is doing tiny, high-frequency wiggles (superharmonics), the air acts like thin, watery syrup (linear damping). If you only modeled the air as one or the other, you'd get the wrong answer. The paper proves that you need both models to get the full picture, as the "mud" dominates the big shakes while the "syrup" handles the tiny ones.

In short, this paper provides a detailed map of the dangerous and fascinating terrain of magnetic bearings. It shows engineers exactly where the "S-shaped" traps are, how to use the derivative knob to escape them, and warns them about the hidden "loop proliferation" that can happen when shaft stiffness and magnetic forces fight each other. It's a guide to keeping high-speed machines from turning into unpredictable, shaking nightmares.

Technical Summary: S-Shaped Bifurcation Cascades and Loop Proliferation in an Active Magnetic Bearing Rotor

Problem Statement
Active Magnetic Bearings (AMBs) are critical for high-speed rotating machinery, offering contactless suspension via proportional-derivative (PD) feedback control. However, the electromagnetic force in AMBs is inherently nonlinear (inversely proportional to the square of the air gap), which can induce complex dynamic behaviors such as multi-valued frequency responses, amplitude jumps, and hysteresis. Existing literature has largely examined electromagnetic nonlinearity, cubic shaft stiffness, and nonlinear aerodynamic drag in isolation or in simple combinations. This study addresses the gap in understanding how these competing nonlinear mechanisms—specifically the softening nature of magnetic forces versus the hardening nature of cubic shaft stiffness, combined with amplitude-dependent aerodynamic drag—interact within a single flexible rotor model. The primary objective is to characterize the resulting bifurcation topology and identify the conditions under which complex phenomena, such as closed-loop branches with numerous saddle-node bifurcations, emerge.

Methodology
The authors model a flexible rotor supported by AMBs as a four-degree-of-freedom system, incorporating:

Nonlinear Electromagnetic Forces: Governed by PD feedback control, where current is a function of journal displacement and velocity.
Cubic Shaft Stiffness: Representing geometric nonlinearity for moderate to large deflections.
Nonlinear Aerodynamic Drag: Modeled as a combination of linear viscous damping (low speed/small amplitude) and quadratic turbulent damping (high speed/large amplitude).

The governing equations of motion are derived using Lagrange formalism and non-dimensionalized. To solve the system, the study employs the Harmonic Quadrature Method (HQM), a time-domain discretization technique that avoids costly Galerkin integrals while maintaining high-order accuracy. This method is combined with pseudo-arc-length continuation to track both stable and unstable periodic solution branches across parameter spaces and Floquet stability analysis to verify stability and detect saddle-node bifurcations. The numerical framework is validated against convergence studies regarding harmonic order and arc-length step size.

Key Contributions and Results
The study systematically investigates the influence of four key parameters: proportional gain ( $P$ ), derivative gain ( $D$ ), nonlinear stiffness coefficient ( $\eta_{NL}$ ), and aerodynamic damping coefficients.

Proportional Gain ( $P$ ) and S-Shaped Cascades:
Increasing $P$ drives a progressive bifurcation cascade through four distinct regimes:
- Linear Regime ( $P < 2.68$ ): Single-valued, stable response.
- 1S-Shaped Regime ( $2.69 \le P \le 2.75$ ): Emergence of a single S-shaped fold with two saddle-node bifurcations.
- 2S-Shaped Regime ( $2.76 \le P \le 3.0$ ): Two S-shaped branches, increasing multi-stability.
- 3S-Shaped Regime ( $3.1 \le P \le 4.1$ ): Three S-shaped branches with eight fold points.
  As $P$ increases from 2 to 4, the primary resonance amplitude rises from 0.2328 to 0.3839, and the multi-stable frequency band widens significantly.
Derivative Gain ( $D$ ) and Inverse Cascade:
The derivative gain acts as a stabilizing countermeasure. Increasing $D$ produces an exact inverse cascade to that of $P$ :
- At low $D$ ( $0.5 \le D \le 1.16$ ), the system remains in the 3S-shaped regime.
- At intermediate $D$ ( $1.16 \le D \le 1.23$ ), the system transitions to a 1S-shaped regime.
- For $D \ge 1.23$ , the multi-stability is destroyed, restoring quasi-linear behavior.
  This reveals a symmetric, complementary relationship where $P$ induces complexity and $D$ systematically suppresses it.
Nonlinear Stiffness ( $\eta_{NL}$ ) and Loop Proliferation:
When the system is fixed in the 3S-shaped regime and $\eta_{NL}$ is increased, a hardening mechanical nonlinearity competes with the softening magnetic nonlinearity. At a critical threshold ( $\eta_{NL} \approx 0.55$ ), the system generates a connected closed-loop branch superimposed on the main response curve.
- This loop contains more than 36 Floquet-verified saddle-node bifurcation points spanning the frequency interval $\Omega \in [0.2466, 0.2922]$ .
- Force comparison confirms this is a physical phenomenon: the ratio of cubic-to-magnetic force increases from 0.002% to 8.6% at the loop onset, indicating the hardening stiffness becomes strong enough to distort the softening magnetic response.
Aerodynamic Damping Regime Transition:
The study identifies an amplitude-dependent transition in damping dominance:
- Nonlinear turbulent damping governs the primary resonance, effectively suppressing large-amplitude S-shaped loops.
- Linear viscous damping controls superharmonic peaks, where amplitudes are smaller.
  A crossover occurs near $\Omega \approx 0.27$ . Models using only one damping mechanism would systematically mispredict response amplitudes (e.g., a linear-only model overestimates primary resonance amplitude by ~57% compared to a nonlinear model).

Significance and Claims
The paper claims to provide the first comprehensive characterization of AMB rotor dynamics under the simultaneous action of softening electromagnetic forces, hardening cubic stiffness, and amplitude-dependent drag. The authors assert that:

The identified S-shaped bifurcation cascades offer a quantitative design chart for tuning PD controllers, defining specific thresholds (e.g., $P < 2.68$ ) to guarantee single-valued operation.
The inverse cascade by $D$ establishes a quantitative stability boundary ( $D \ge 1.23$ ) for eliminating multi-stability, offering a practical control strategy for systems operating in high-gain regimes.
The discovery of the connected closed-loop branch with >36 fold points represents a complex topological transformation distinct from simpler S-shaped folds, highlighting the unpredictability of rotor response in specific parameter ranges due to history dependence.
The dual-damping formulation is presented as a modeling necessity for high-speed AMB rotors to accurately capture both primary and superharmonic behaviors.

The study concludes that these findings characterize the nonlinear bifurcation topology of AMB systems, providing direct implications for the safe tuning of PD controllers and the prediction of rotor dynamics in the presence of combined nonlinearities. Future work is suggested to explore passive stabilization (squeeze film dampers) and active control strategies to mitigate the identified loop proliferation.

More like this