⚡ electrical engineering

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

Este artículo utiliza el Método de Cuadratura Armónica y el análisis de estabilidad de Floquet para revelar cómo las ganancias proporcionales y derivativas, la rigidez cúbica del eje y la resistencia aerodinámica no lineal interactúan para impulsar cascadas de bifurcación en forma de S complejas, una proliferación extensiva de bucles con más de 36 bifurcaciones de silla-nodo y regímenes de amortiguamiento distintos en rotores flexibles soportados por cojinetes magnéticos activos.

Autores originales: Farouk Thaljaoui, Thabet Guesmi, Thaljaoui Wathek

Publicado 2026-08-05

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Farouk Thaljaoui, Thabet Guesmi, Thaljaoui Wathek

Artículo original bajo licencia CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Imagina un mundo donde las máquinas pesadas y giratorias no tocan el suelo en absoluto. En lugar de rodar sobre ruedas o deslizarse sobre aceite, flotan en el aire, suspendidas por manos magnéticas invisibles. Esta es la magia de los Rodamientos Magnéticos Activos (AMB, por sus siglas en inglés), una tecnología utilizada en todo, desde motores de reacción de alta velocidad hasta volantes de inercia para el almacenamiento de energía. Pero aquí está el truco: los imanes son complicados. La fuerza con la que tiran no crece en línea recta; se vuelve increíblemente fuerte cuanto más cerca están las cosas, como una banda de goma que retrocede con más fuerza cuanto más se estira. Cuando combinas este "latigazo" magnético con el bamboleo natural de un eje giratorio y la resistencia del aire de un objeto que se mueve rápido, el sistema puede empezar a comportarse como un compañero de baile caótico. Podría dar un salto repentino a un temblor enorme y peligroso, o quedarse atrapado en un bucle donde no sabe hacia qué lado ir. Los ingenieros necesitan entender estos saltos y bucles para evitar que sus máquinas se despedacen por el propio temblor.

Este artículo profundiza en esa danza caótica, utilizando una poderosa simulación computacional para mapear exactamente cómo se comporta un rotor flotante cuando se ajustan sus controles. Los investigadores trataron al rotor como una vara flexible con un disco pesado en el medio, sostenida por rodamientos magnéticos que utilizan un "controlador PD" (una forma elegante de decir un cerebro inteligente que ajusta la fuerza del imán basándose en la posición y la velocidad). Añadieron dos ingredientes extra: el hecho de que el eje se vuelve más rígido a medida que se dobla (como un resorte que se vuelve más difícil de presionar cuanto más lo presionas) y la resistencia del aire que cambia dependiendo de qué tan rápido se mueva el rotor. Al ejecutar miles de simulaciones, descubrieron que el "cerebro" de la máquina tiene dos perillas: una Perilla Proporcional (P) que controla con qué fuerza tiran los imanes, y una Perilla Derivativa (D) que controla cuánta resistencia ofrece al movimiento.

El equipo descubrió que aumentar la Perilla Proporcional (P) hace que el sistema sea cada vez más salvaje. Con ajustes bajos, el rotor se comporta bien, como un péndulo tranquilo. Pero a medida que subían P, la curva suave del movimiento del rotor comenzó a plegarse sobre sí misma, creando "formas en S". Es como una pista de montaña rusa que de repente vuelve sobre sí misma en un bucle. Al principio, había un bucle (forma de 1S), luego dos (forma de 2S) y finalmente tres (forma de 3S). En estas zonas con forma de S, el rotor tiene múltiples opciones para vibrar a la misma velocidad. Podría estar vibrando suavemente o sacudiéndose violentamente, y un pequeño empujón podría hacerlo saltar de lo seguro a lo peligroso instantáneamente. Los investigadores calcularon que, al aumentar P de 2 a 4, la amplitud máxima de la sacudida creció de 0.2328 a 0.3839, y el número de estos "puntos de salto" (llamados bifurcaciones de nodo de silla) se multiplicó, creando una amplia zona de confusión donde el comportamiento de la máquina depende enteramente de su historia.

Sin embargo, la Perilla Derivativa (D) actúa como un superhéroe opuesto. Mientras que P crea el caos, D lo elimina. El artículo muestra que aumentar D actúa como echar agua sobre un fuego, suavizando esos salvajes bucles en forma de S. Cuando D es bajo, el sistema es un desastre de bucles. Pero a medida que los investigadores subían D, los bucles se encogían. En un ajuste específico de D = 1.23, el caos desapareció por completo y el sistema regresó a un estado tranquilo y predecible. Es como si la perilla P construyera un laberinto, y la perilla D fuera la llave que abre la puerta hacia un camino recto.

El descubrimiento más asombroso ocurrió cuando activaron el "endurecimiento" del eje (la rigidez cúbica). Cuando combinaron el "ablandamiento" magnético con el "endurecimiento" del eje, algo extraño sucedió: los bucles en forma de S no solo se plegaron; proliferaron en un gigantesco bucle cerrado y conectado. Imagina una bola de estambre enredada que, de alguna manera, se conecta consigo misma. En esta zona específica, los investigadores encontraron una única rama de la solución que giraba y se retorcía con más de 36 diferentes "puntos de pliegue" (lugares donde el camino dobla hacia atrás). Esto sucedió cuando la rigidez del eje fue lo suficientemente fuerte como para competir con los imanes, alcanzando una proporción donde la fuerza del eje era aproximadamente el 8.6% de la fuerza magnética. En esta zona, el rotor tiene docenas de formas diferentes de vibrar de manera estable a la misma velocidad, lo que hace que su comportamiento sea increíblemente difícil de predecir sin un mapa perfecto.

Finalmente, el equipo analizó la resistencia del aire. Descubrieron que el aire lucha contra el rotor de manera diferente dependiendo de qué tan fuerte esté sacudiéndose. Cuando el rotor se sacude violentamente (la resonancia principal), el aire actúa como un lodo espeso y pegajoso (amortiguamiento no lineal), lo cual es excelente para detener los grandes saltos. Pero cuando el rotor realiza pequeños temblores de alta frecuencia (superarmónicos), el aire actúa como un jarabe fino y acuoso (amortiguamiento lineal). Si solo modelaran el aire como uno u otro, obtendrían la respuesta incorrecta. El artículo demuestra que se necesitan ambos modelos para obtener la imagen completa, ya que el "lodo" domina las grandes sacudidas mientras que el "jarabe" se encarga de las pequeñas.

En resumen, este artículo proporciona un mapa detallado del terreno peligroso y fascinante de los rodamientos magnéticos. Muestra a los ingenieros exactamente dónde están las trampas en "forma de S", cómo usar la perilla derivativa para escapar de ellas, y les advierte sobre la "proliferación de bucles" oculta que puede ocurrir cuando la rigidez del eje y las fuerzas magnéticas luchan entre sí. Es una guía para evitar que las máquinas de alta velocidad se conviertan en pesadillas impredecibles de sacudidas constantes.

Resumen Técnico: Cascadas de Bifurcación en Forma de S y Proliferación de Bucles en un Rotor con Cojinetes Magnéticos Activos

Planteamiento del Problema
Los Cojinetes Magnéticos Activos (AMB, por sus siglas en inglés) son críticos para la maquinaria rotativa de alta velocidad, ofreciendo una suspensión sin contacto mediante control de retroalimentación proporcional-derivativo (PD). Sin embargo, la fuerza electromagnética en los AMB es inherentemente no lineal (inversamente proporcional al cuadrado del entrehierro), lo que puede inducir comportos dinámicos complejos tales como respuestas de frecuencia multivaluadas, saltos de amplitud e histéresis. La literatura existente ha examinado mayoritariamente la no linealidad electromagnética, la rigidez cúbica del eje y el arrastre aerodinámico no lineal de forma aislada o en combinaciones simples. Este estudio aborda la brecha en la comprensión de cómo estos mecanismos de no linealidad competitivos —específicamente la naturaleza de ablandamiento de las fuerzas magnéticas frente a la naturaleza de endurecimiento de la rigidez cúbica del eje, combinadas con un arrastre aerodinámico dependiente de la amplitud— interactúan dentro de un único modelo de rotor flexible. El objetivo principal es caracterizar la topología de bifurcación resultante e identificar las condiciones bajo las cuales emergen fenómenos complejos, tales como ramas de bucle cerrado con numerosos puntos de bifurcación de silla-nodo.

Metodología
Los autores modelan un rotor flexible soportado por AMB como un sistema de cuatro grados de libertad, incorporando:

Fuerzas Electromagnéticas No Lineales: Gobernadas por el control de retroalimentación PD, donde la corriente es una función del desplazamiento y la velocidad del jurnal.
Rigidez Cúbica del Eje: Que representa la no linealidad geométrica para deflexiones moderadas a grandes.
Arrastre Aerodinámico No Lineal: Modelado como una combinación de amortiguamiento viscoso lineal (baja velocidad/pequeña amplitud) y amortiguamiento turbulento cuadrático (alta velocidad/gran amplitud).

Las ecuaciones de movimiento gobernantes se derivan mediante el formalismo de Lagrange y se desdimensionalizan. Para resolver el sistema, el estudio emplea el Método de Cuadratura Armónica (HQM), una técnica de discretización en el dominio del tiempo que evita las costosas integrales de Galerkin manteniendo una alta precisión de orden superior. Este método se combina con la continuación de pseudo-longitud de arco para rastrear tanto las ramas de solución periódica estables como las inestables a través de los espacios de parámetros, así como el análisis de estabilidad de Floquet para verificar la estabilidad y detectar bifurcaciones de silla-nodo. El marco numérico se valida mediante estudios de convergencia respecto al orden armónico y al tamaño del paso de longitud de arco.

Contribuciones Clave y Resultados
El estudio investiga sistemáticamente la influencia de cuatro parámetros clave: la ganancia proporcional ( $P$ ), la ganancia derivativa ( $D$ ), el coeficiente de rigidez no lineal ( $\eta_{NL}$ ) y los coeficientes de amortiguamiento aerodinámico.

Ganancia Proporcional ( $P$ ) y Cascadas en Forma de S:
Aumentar $P$ impulsa una cascada de bifurcación progresiva a través de cuatro regímenes distintos:

Régimen Lineal ( $P < 2.68$ ): Respuesta de valor único y estable.
Régimen con 1S en forma de S ( $2.69 \le P \le 2.75$ ): Emergencia de un solo pliegue en forma de S con dos bifurcaciones de silla-nodo.
Régimen con 2S en forma de S ( $2.76 \le P \le 3.0$ ): Dos ramas en forma de S, aumentando la multiestabilidad.
Régimen con 3S en forma de S ( $3.1 \le P \le 4.1$ ): Tres ramas en forma de S con ocho puntos de pliegue.
A medida que $P$ aumenta de 2 a 4, la amplitud de la resonancia primaria sube de 0.2328 a 0.3839, y la banda de multiestabilidad se ensancha significamente.

Ganancia Derivativa ( $D$ ) y Cascada Inversa:
La ganancia derivativa actúa como una contramedida estabilizadora. Aumentar $D$ produce una cascada inversa exacta a la de $P$ :

Con un $D$ bajo ( $0.5 \le D \le 1.16$ ), el sistema permanece en el régimen de 3S en forma de S.
Con un $D$ intermedio ( $1.16 \le D \le 1.23$ ), el sistema transiciona a un régimen de 1S en forma de S.
Para $D \ge 1.23$ , la multiestabilidad se destruye, restaurando un comportamiento cuasi-lineal.
Esto revela una relación simétrica y complementaria donde $P$ induce complejidad y $D$ la suprime sistemáticamente.

Rigidez No Lineal ( $\eta_{NL}$ ) y Proliferación de Bucles:
Cuando el sistema se fija en el régimen de 3S en forma de S y se incrementa $\eta_{NL}$ , una no linealidad mecánica de endurecimiento compite con la no linealidad magnética de ablandamiento. En un umbral crítico ( $\eta_{NL} \approx 0.55$ ), el sistema genera una rama de bucle cerrado conectado superpuesta en la curva de respuesta principal.

Este bucle contiene más de 36 puntos de bifurcación de silla-nodo verificados por Floquet que abarcan el intervalo de frecuencia $\Omega \in [0.2466, 0.2922]$ .
La comparación de fuerzas confirma que este es un fenómeno físico: la relación de la fuerza cúbica sobre la magnética aumenta de 0.002% a 8.6% en el inicio del bucle, indicando que la rigidez de endurecimiento se vuelve lo suficientemente fuerte como para distorsionar la respuesta de ablandamiento magnético.

Transición de Régimen de Amortiguamiento Aerodinámico:
El estudio identifica una transición dependiente de la amplitud en la dominancia del amortiguamiento:

El amortiguamiento turbulento no lineal gobierna la resonancia primaria, suprimiendo eficazmente los bucles en forma de S de gran amplitud.
El amortiguamiento viscoso lineal controla los picos superarmónicos, donde las amplitudes son menores.
Ocurre un cruce cerca de $\Omega \approx 0.27$ . Los modelos que utilizan un solo mecanismo de amortiguamiento subestiman o sobreestiman sistemáticamente las amplitudes de respuesta (por ejemplo, un modelo puramente lineal sobreestima la amplitud de la resonancia primaria en un ~57% comparado con un modelo no lineal).

Significancia y Reivindicaciones
El artículo afirma proporcionar la primera caracterización exhaustiva de la dinámica de rotores AMB bajo la acción simultánea de fuerzas electromagnéticas de ablandamiento, rigidez cúbica de endurecimiento y arrastre dependiente de la amplitud. Los autores sostienen que:

Las cascadas de bifurcación en forma de S identificadas ofrecen un cuadro de diseño cuantitativo para el ajuste de controladores PD, definiendo umbrales específicos (por ejemplo, $P < 2.68$ ) para garantizar una operación de valor único.
La cascada inversa mediante $D$ establece un límite de estabilidad cuantitativo ( $D \ge 1.23$ ) para eliminar la multiestabilidad, ofreciendo una estrategia de control práctica para sistemas que operan en regímenes de alta ganancia.
El descubrimiento de la rama de bucle cerrado conectado con más de 36 puntos de pliegue representa una transformación topológica compleja distinta de los pliegues simples en forma de S, resaltando la imprevisibilidad de la respuesta del rotor en rangos de parámetros específicos debido a la dependencia de la historia.
La formulación de doble amortiguamiento se presenta como una necesidad de modelado para capturar con precisión tanto los comportamientos de la resonancia primaria como los superarmónicos en rotores AMB de alta velocidad.

El estudio concluye que estos hallazgos caracterizan la topología de bifurcación no lineal de los sistemas AMB, proporcionando implicaciones directas para el ajuste seguro de los controladores PD y la predicción de la dinámica del rotor en presencia de no linealidades combinadas. Se sugiere como trabajo futuro explorar la estabilización pasiva (amortiguadores de película de presión/squeeze film dampers) y estrategias de control activo para mitigar la proliferación de bucles identificada.

Más como este