Autoparametric Resonance and Fractional-Order Viscoelastic Damping in Cable-Stayed Bridges: A Multiple-Scales Analysis with Finite-Element Verification

Resumen Técnico: Resonancia Autoparamétrica y Amortiguamiento Viscoelástico de Orden Fraccionario en Puentes Atirantados

Planteamiento del Problema
La dinámica de los puentes atirantados se caracteriza por dos fenómenos distintos y bien estudiados que, históricamente, se han analizado por separado:

1. Resonancia Autoparamétrica (Interna): Un mecanismo no lineal donde la energía se transfiere de un modo del tablero (viga) excitado globalmente a un modo local del tirante cuando sus frecuencias naturales están sintonizadas cerca de una relación 2:1 ($\omega_p \approx 2\omega_q$). Esto puede conducir a oscilaciones de gran amplitud en los tirantes incluso bajo una excitación moderada del tablero.
2. Amortiguamiento Viscoelástico de Orden Fraccionario: La evidencia experimental muestra que los amortiguadores de los tirantes exhiben una disipación dependiente de la frecuencia que los modelos clásicos de viscosidad de orden entero no pueden reproducir con precisión. Los modelos de derivada fraccionaria (específicamente las derivadas de Caputo) se han aplicado con éxito al amortiguamiento estructural, pero aún no se han integrado en el marco específico de la resonancia autoparamétrica entre tablero y tirante.

Aunque existen metodologías de perturbación que combinan el cálculo fraccionario con la resonancia interna (mediante el Método de las Escalas Múltiples) para puentes colgantes y osciladores genéricos, y existen modelos de resonancia entre tablero y tirante con amortiguamiento clásico, ningún estudio ha combinado estos elementos para derivar un umbral de inestabilidad de forma cerrada para puentes atirantados que incorpore el amortiguamiento de orden fraccionario.

Metodología
El artículo formula un modelo no lineal de dos modos que acopla un modo de tablero forzado ($p$) y un modo de tirante ($q$). Las ecuaciones gobernantes incluyen términos de acoplamiento geométrico cuadrático y un término de amortiguamiento de orden fraccionario de Caputo ($D^\alpha_t$) aplicado específicamente al modo del tirante.

1. Enfoque Analítico: Se aplica el Método de las Escalas Múltiples (MMS) al sistema débilmente no lineal. Un paso clave consiste en la descomposición armónica de la derivada fraccionaria que actúa sobre el modo del tirante. Esta descomposición divide el término fraccionario en dos contribuciones físicas distintas:
   • Un término de amortiguamiento efectivo ($\mu_{q,eff}$), que escala con $\sin(\alpha\pi/2)$.
   • Un desplazamiento de frecuencia efectivo ($\Delta\omega_{q,eff}$), que escala con $\cos(\alpha\pi/2)$.
2. Derivación: El análisis produce ecuaciones de modulación de amplitud-fase en forma cerrada. A partir de estas, el autor deriva la rama de respuesta "trivial" (solo del tablero) y, crucialmente, un umbral de inestabilidad autoparamétrica en forma cerrada. Este umbral define la amplitud del tablero por encima de la cual el modo del tirante se activa espontáneamente.
3. Verificación Numérica:
   • Integración Directa: Las ecuaciones fraccionarias no lineales originales se integran directamente utilizando una discretización de Grünwald-Letnikov junto con un esquema de marcha temporal implícito-explícito (IMEX) para verificar las tasas de crecimiento analíticas y el comportamiento de saturación.
   • Extensión de Elementos Finitos (FE): Se ensambla un modelo de puente de cables múltiples plano (29 nodos, 81 grados de libertad) para verificar que la formulación del amortiguamiento fraccionario sigue siendo numéricamente tratable a una escala de múltiples grados de libertad, independientemente del acoplamiento no lineal en esta etapa.

Contribuciones Clave

• Síntesis: El artículo integra con éxito el amortiguamiento de orden fraccionario en el marco establecido de la resonancia autoparamétrica entre tablero y tirante, un vacío previamente no abordado en la literatura.
• Umbral de Inestabilidad en Forma Cerrada: El principal resultado analítico es la Ecuación (12), que proporciona una expresión en forma cerrada para la amplitud crítica del tablero ($a_{0,crit}$) requerida para activar la inestabilidad del tirante. Esta expresión depende explícitamente del orden fraccionario $\alpha$ a través tanto del amortiguamiento efectivo como del desplazamiento de frecuencia efectivo.
• Implicación de Diseño No Monotónica: El análisis revela que aumentar el orden fraccionario $\alpha$ (que generalmente aumenta la disipación aparente) no necesariamente eleva el umbral de inestabilidad. Debido a los efectos competitivos del término de desplazamiento de frecuencia, un $\alpha$ más alto puede, de hecho, reducir el umbral de inestabilidad en ciertas configuraciones. Esto desafía la suposición de que maximizar el amortiguamiento por sí solo maximiza la estabilidad contra la resonancia autoparamétrica.
• Validación de la Saturación: El estudio confirma numéricamente la firma clásica de la "saturación de amplitud" de la resonancia autoparamétrica: una vez que el tirante se activa, la amplitud del tablero se fija en el valor del umbral mientras el exceso de energía se desvía hacia el tirante.

Resultados

• Predicciones Analíticas: El umbral de inestabilidad derivado muestra que el orden fraccionario $\alpha$ gobierna el inicio de la respuesta de gran amplitud del tirante. Para los parámetros probados, el umbral disminuye monótonamente a medida que $\alpha$ aumenta de 0.2 a 0.8.
• Acuerdo Numérico: La integración numérica directa de las ecuaciones no lineales confirma las predicciones analíticas. La tasa de crecimiento lineal predicha para el modo del tirante coincide con los resultados de la simulación con un error inferior al 1% para casos representativos.
• Escalabilidad de FE: La extensión de elementos finitos demuestra que la matriz de amortiguamiento fraccionario puede ensamblarse y resolverse eficientemente para un sistema de cables múltiples (81 DOFs), completándose en menos de un minuto en una estación de trabajo estándar. Sin embargo, el modelo de FE utilizado para esta verificación fue lineal y aún no incluía los términos de acoplamiento autoparamétrico cuadráticos.

Significado y Reclamaciones
El artículo se posiciona como una contribución metodológica y analítica en la intersección de dos campos maduros. No pretende haber descubierto una nueva teoría de perturbación ni haber realizado una validación experimental. Su importancia radica en:

1. Derivar un umbral específico orientado al diseño: Proporcionar una fórmula que vincula el orden fraccionario de un amortiguador directamente con los límites de estabilidad de un puente atirantado, resaltando que la selección del amortiguador debe tener en cuenta tanto el efecto de amortiguamiento como el de desplazamiento de frecuencia.
2. Demostrar la viabilidad: Probar que el cálculo fraccionario puede aplicarse rigurosamente a la dinámica compleja de los puentes atirantados y que la maquinaria resultante escala a modelos estructurales realistas.
3. Identificar el trabajo futuro: El autor establece explícitamente que los siguientes pasos consisten en mapear los límites de estabilidad secundaria, incorporar el acoplamiento no lineal en el modelo de FE y, lo más crítico, obtener la validación experimental física utilizando un banco de pruebas equipado con un amortiguador de orden fraccionario caracterizado.

El trabajo mantiene un alcance modesto en su alcance, reconociendo que los resultados actuales son aproximaciones de primer orden válidas para una no linealidad débil y que la validación experimental de la extensión fraccionaria en sí misma es una prioridad para la investigación futura.