⚡ electrical engineering

Time-Domain Two-Way Fluid–Structure Interaction Analysis of the HIRENASD Transonic Wing: Static Aeroelastic Response and Flutter Prediction

Este estudio valida un marco de Interacción Fluido-Estructura en el dominio del tiempo y de acoplamiento débil en ANSYS Workbench frente a los bancos de pruebas de la NASA para predecir con precisión la respuesta aeroelástica estática y los límites de flameo del ala transónica HIRENASD, demostrando su viabilidad como una herramienta reproducible y accesible para la industria en la evaluación preliminar de alas de aeronaves.

Autores originales: Siddalingappa P Kodigaddi

Publicado 2026-07-13

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Siddalingappa P Kodigaddi

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 ala gigante y flexible planeando por el cielo, bailando con el viento. Ahora, imagina que ese viento no solo pasa a su lado; sino que en realidad está empujando y tirando del ala, mientras el ala se dobla y se retuerce de vuelta, cambiando la forma en que fluye el viento. Esta danza invisible se llama Interacción Fluido–Estructura (FSI), y cuando se sale de control, el ala puede empezar a sacudirse violentamente hasta romperse. Este temible temblor se llama flameo (o flutter).

En este estudio, un investigador llamado Siddalingappa P. Kodigaddi construyó un túnel de viento virtual súper preciso dentro de una computadora para observar cómo ocurre esta danza en un ala específica llamada HIRENASD. Este no es un ala cualquiera; es un ala "supercrítica" de alta tecnología diseñada para jets transónicos rápidos (volando justo por debajo de la velocidad del sonido). El objetivo era ver si un programa de computadora estándar y comercial (ANSYS) podía predecir exactamente cuándo este ala comenzaría a flamear, sin necesidad de un código personalizado y sumamente costoso.

La prueba del túnel de viento virtual

Primero, el investigador tuvo que demostrar que su túnel de viento virtual era real. Configuró una simulación a una velocidad de Mach 0.8 (unas 614 mph) y una densidad de aire muy alta (número de Reynolds de 7×10⁶). Comparó los resultados de su computadora con un famoso punto de referencia de la NASA.

El Resultado: La computadora predijo la sustentación, la resistencia y las fuerzas de torsión con una precisión increíble, fallando por menos del 4% respecto a los números de la NASA. Fue como adivinar el peso de un gato y fallar por solo unos pocos gramos.
El Ritmo Natural del Ala: Antes de que soplara el viento, el investigador comprobó cómo vibraba el ala por sí sola (como al pulsar la cuerda de una guitarra). La frecuencia principal de flexión del ala se midió en la vida real en 26.0 Hz. La computadora predijo 26.504 Hz. Esa es una diferencia de solo un 1.94%, lo que demuestra que el ala virtual era rígida y pesada tal como la real.

El estiramiento estático (Antes del temblor)

Después, empujaron el ala a diferentes ángulos (desde -1.5° hasta +4.5°) para ver cuánto se doblaba bajo una presión constante.

La Flexión: En el ángulo más alto (4.5°) y velocidad (Mach 0.8), la punta del ala se dobló hacia abajo 12 mm.
La Torsión: Debido a que el ala está retrasada (angulada como un bumerán), esta flexión en realidad torció ligeramente el ala, lo que ayudó a reducir la carga en la raíz. Es como la rama flexible de un árbol que se dobla con el viento para evitar romperse. La computadora coincidió con estos números de flexión casi perfectamente con otras simulaciones de alto nivel.

El Gran Momento: Encontrando la Velocidad de Flameo

Este es el evento principal. El investigador desactivó el "amortiguamiento" (la fricción interna que normalmente detiene el temblor) y aumentó lentamente la velocidad del viento en la simulación para ver cuándo el ala se volvería loca. Probaron dos "atmósferas" diferentes:

Gran Altitud (Aire Delgado): Velocidad Mach 0.8 con una densidad de 0.2525 kg/m³.
Baja Altitud (Aire Espeso): Velocidad Mach 0.7 con una densidad de 0.50 kg/m³.

Los Hallazgos:

A Mach 0.8 (Aire Delgado): El ala se mantuvo tranquila hasta que la velocidad alcanzó los 266.56 m/s. En esta velocidad exacta, el ala comenzó a sacudirse con un ritmo constante y estable (sin crecer ni morir). Esta es la velocidad de flameo. El temblor ocurrió a una frecuencia de aproximadamente 22 Hz.
A Mach 0.7 (Aire Espeso): El ala se mantuvo tranquila hasta los 230.46 m/s. Aquí, la frecuencia de temblor fue de unos 25.5 Hz.

El Momento de Revelación ("Aha!"):
Podrías pensar: "Espera, el ala flameó a una velocidad mayor en el aire delgado (266.56 m/s) que en el aire espeso (230.46 m/s). ¿Significa eso que es más segura a grandes altitudes?".
No exactamente. El artículo explica que, aunque la velocidad fue mayor, la fuerza del viento (presión dinámica) era en realidad mucho más fuerte en el caso del aire espeso.

La "presión de flameo" en el caso de menor velocidad (aire espeso) fue de 13.27 kPa.
La "presión de flameo" en el caso de mayor velocidad (aire delgado) fue de solo 8.97 kPa.
Esto significa que el ala está en realidad 47.9% más estresada en el escenario de aire espeso, a pesar de que se movía más lento. El artículo confirma una regla general: el peligro escala con la densidad multiplicada por la velocidad al cuadrado. Así que, un ala que es segura a baja altitud puede volverse peligrosa al ascender, incluso si mantiene la misma velocidad, porque el aire se vuelve más delgado y la "rigidez" aerodinámica disminuye.

¿Qué tipo de temblor fue?

El investigador observó de cerca el temblor para ver cómo sucedía.

Oscilación de un Solo Modo: El ala no se retorció y se dobló junta en un nudo complejo. En cambio, fue un movimiento de flexión de arriba hacia abajo simple (como un trampolín).
¿Por qué? La frecuencia de torsión del ala es mucho más alta (279.13 Hz) que su frecuencia de flexión (26.504 Hz). Están tan alejadas que nunca se "fusionan" (coalescen) para crear un flameo complejo y mortal. Es una flexión simple de un solo modo.
El Desplazamiento: A la velocidad más alta (Mach 0.8), la presión del aire en realidad hizo que el ala se sintiera más "suave", bajando su frecuencia de temblor de 26.504 Hz a 22 Hz. Al Mach 0.7 (aire menos "suave"), la frecuencia solo bajó a 25.5 Hz.

Qué significa esto para el futuro

El artículo no pretende haber "resuelto" el flameo para todos los aviones para siempre. En cambio, muestra que el software comercial disponible en el mercado (ANSYS Fluent y ANSYS Mechanical) puede utilizarse para predecir estas velocidades peligrosas con alta confianza, sin necesidad de códigos secretos de uso exclusivo gubernamental.

La Confianza: Los resultados se basan en simulaciones que han sido rigurosamente verificadas contra datos del mundo real y otros puntos de referencia de supercomputadoras. Los autores estiman que su método es preciso dentro de un margen del 5–8%.
El Límite: La simulación no tuvo en cuenta el frío extremo del túnel de viento (que hace que el metal sea ligeramente más blando) o cambios de temperatura complejos, pero el modo de flexión principal fue lo suficientemente preciso como para confiar en la predicción del flameo.

En resumen, este estudio construyó un gemelo digital de un ala de alta tecnología, demostró que se comporta como la real y lo utilizó para encontrar la velocidad exacta donde el ala comienza a bailar sin control. Mostró que, aunque el ala pueda volar más rápido antes de temblar en aire delgado, la fuerza del viento es el verdadero peligro, y este nuevo método computacional puede ayudar a los ingenieros a diseñar alas más seguras para el futuro.

Resumen Técnico: Análisis de Interacción Fluido-Estructura Bidireccional en el Dominio del Tiempo del Ala Transónica HIRENASD

Planteamiento del Problema
La integridad estructural de las alas de aeronaves de transporte en el régimen transónico está gobernada por acoplamientos aeroelásticos complejos, particularmente el flutter (flameo), que puede conducir a fallos catastróficos. Los métodos lineales en el dominio de la frecuencia, como el Método de la Doble Red (DLM), a menudo carecen de la fidelidad necesaria para capturar fenómenos no lineales críticos inherentes al flujo transónico, incluyendo ondas de choque móviles, separación de la capa límite inducida por choque y fuertes interacciones fluido-estructura. Si bien el proyecto High Reynolds Number Aero-Structural Dynamics (HIRENASD) proporciona una base de datos experimental exhaustiva para un ala supercrítica (perfil BAC3-11), persisten brechas significativas en la literatura abierta con respecto a la validación de estas configuraciones utilizando solvers CFD-CSD comerciales de propósito general. Además, la identificación sistemática en el dominio del tiempo de los límites de flutter a través de múltiples números de Mach y el análisis cuantitativo de cómo la densidad de flujo libre y el número de Mach influyen conjuntamente en la presión dinámica de flutter no han sido abordados completamente para este estándar de referencia.

Metodología
El estudio emplea un marco de Interacción Fluido-Estructura (FSI) bidireccional, débilmente acoplado, en el dominio del tiempo implementado dentro de ANSYS Workbench. El enfoque particiona la solución entre:

CFD: ANSYS Fluent (solver basado en densidad) utilizando las ecuaciones de Navier–Stokes promediadas por Reynolds (RANS) con el modelo de turbulencia SST $k-\omega$ . El dominio fluido cuenta con una malla tetraédrica no estructurada (~600,000 celdas) con capas prismáticas para lograr un $y^+ \approx 1$ .
CSD: ANSYS Mechanical utilizando un modelo de elementos finitos hexaédrico (41,923 elementos) de la estructura del ala HIRENASD (acero maraging de níquel). El análisis modal utiliza el algoritmo Block-Lanczos, y la integración transitoria emplea el método Newmark- $\beta$ .
Acoplamiento: ANSYS System Coupling gestiona el intercambio de datos mediante la Interfaz de Rejilla General (GGI) para la conservación de fuerzas y la interpolación Smart Bucket para el mapeo de desplazamientos. La deformación de la malla se gestiona mediante suavizado por analogía de resortes con remallado local.

El estudio valida el marco de trabajo contra los benchmarks aerodinámicos de NASA FUN3D y las frecuencias estructurales experimentales de HIRENASD. Se realizan análisis aeroelásticos estáticos a números de Mach ( $M$ ) de 0.7 y 0.8 a través de ángulos de ataque desde $-1.5^\circ$ hasta $+4.5^\circ$ . El análisis dinámico de flutter identifica el límite de flutter mediante la variación sistemática de la velocidad de flujo libre para caracterizar los regímenes de respuesta subcrítico (amortiguado), crítico (neutro) y supercrítico (divergente) en dos condiciones de vuelo distintas: $M=0.8$ ( $\rho = 0.2525$ kg/m³) y $M=0.7$ ( $\rho = 0.50$ kg/m³).

Contribuciones Clave

Validación del Marco de Trabajo: El estudio establece con éxito un flujo de trabajo FSI reproducible y accesible industrialmente utilizando exclusivamente herramientas comerciales (ANSYS Fluent, Mechanical, System Coupling) para la aeroelasticidad transónica.
Benchmarking Aeroelástico Estático: El marco reproduce los coeficientes aerodinámicos en estado estacionario ( $C_L, C_D, C_M$ ) dentro de un 1.87% a 3.52% de los datos de referencia de NASA FUN3D y predice la frecuencia de flexión fundamental (26.504 Hz) con un error del 1.94% respecto a los resultados experimentales de la prueba de vibración en tierra.
Caracterización Dinámica del Flutter: El estudio proporciona una determinación sistemática en el dominio del tiempo de los límites de flutter, caracterizando explícitamente la transición de oscilaciones amortiguadas a divergentes. Identifica que el flutter está gobernado por un mecanismo de flexión de un solo modo en ambos regímenes de Mach probados.
Perspectiva Física sobre el Escalamiento: La investigación cuantifica la relación entre el número de Mach, la densidad de flujo libre y la presión dinámica de flutter, demostrando que, si bien la velocidad de flutter puede ser mayor a densidades más bajas (mayor altitud), la presión dinámica de flutter es significativamente mayor a densidades más altas debido a la ley de escalamiento $q_f \propto \rho V_f^2$ .

Resultos

Respuesta Estática: A $M=0.8$ , la pendiente de la curva de sustentación es aproximadamente 0.107 por grado. La deflexión máxima de la punta del ala alcanza los 12 mm (deflexión adimensional $\delta/b = 0.0093$ ) a un ángulo de ataque de $4.5^\circ$ . La polar de resistencia a $M=0.8$ exhibe un aumento agudo de resistencia inducido por choque, mientras que la polar a $M=0.7$ es más parabólica y dominada por la resistencia inducida.
Análisis Modal: El primer modo de flexión fuera del plano se predice en 26.504 Hz (error del 1.94% vs. experimento). Los modos de orden superior muestran una sobrepredicción sistemática (9–13%) atribuida a simplificaciones estructurales (por ejemplo, falta de cumplimiento de los sujetadores) en el modelo de elementos finitos.
Límites de Flutter:
- A $M=0.8$ ( $\rho = 0.2525$ kg/m³): La velocidad de flutter ( $V_f$ ) es de 266.56 m/s con una presión dinámica de flutter ( $q_f$ ) de 8.97 kPa. La frecuencia de flutter es de ~22 Hz, lo que representa una reducción del 16.9% respecto a la frecuencia estructural en el vacío debido a la modificación de la rigidez aerodinámica transónica.
- A $M=0.7$ ( $\rho = 0.50$ kg/m³): La velocidad de flutter es de 230.46 m/s con una presión dinámica de flutter ( $q_f$ ) de 13.27 kPa. La frecuencia de flutter es de ~25.5 Hz, mostrando una reducción menor (3.8%) de la frecuencia estructural, lo que indica una no linealidad transónica más débil.
Ley de Escalamiento: A pesar de la menor velocidad de flutter a $M=0.7$ , la presión dinámica de flutter es un 47.9% mayor que a $M=0.8$ debido a la densidad de flujo libre un 98% superior. Esto confirma el escalamiento dimensional $q_f \propto \rho V_f^2$ .

Significancia y Reivindicaciones
El artículo sostiene que el marco de trabajo validado de ANSYS Workbench ofrece una vía práctica y relevante para la industria para la evaluación preliminar de los límites de flutter de las alas de aeronaves de transporte sin requerir código propio o solvers propietarios. El estudio enfatiza que el uso exclusivo de herramientas comerciales permite la transferencia inmediata a entornos de diseño industrial.

Los autores afirman que la confirmación cuantitativa del escalamiento de la presión dinámica de flutter dependiente de la densidad proporciona una base científicamente rigurosa para extrapolar los resultados de modelos de túnel de viento a las condiciones de vuelo a escala real, particularmente con respecto al riesgo de flutter dependiente de la altitud. El trabajo destaca que un ala puede estar libre de flutter a bajas altitudes (alta densidad) pero acercarse a su límite de flutter en altitud de crucero (baja densidad), incluso a un número de Mach constante, debido a la reducción en la rigidez aerodinámica efectiva. El estudio concluye que las velocidades de flutter identificadas conllevan una incertidumbre metodológica estimada de 5–8%, basada en comparaciones con marcos validados similares para el ala AGARD 445.6, y que el mecanismo de flexión de un solo modo observado es consistente con la gran separación de frecuencias entre los modos de flexión y torsión en esta configuración específica.

La prueba del túnel de viento virtual

El estiramiento estático (Antes del temblor)

El Gran Momento: Encontrando la Velocidad de Flameo

¿Qué tipo de temblor fue?

Qué significa esto para el futuro

Más como este