🔢 mathematics

A 2.5D NURBS-Trace Infinite-Element Method for Moving-Load Wave Propagation and Soil--Structure Interaction in Semi-Infinite Ground

Este artículo presenta un método de elementos infinitos de traza NURBS 2.5D (NBIEM) que combina el análisis isogeométrico para el campo cercano con NURBS de producto tensorial y funciones radiales exponenciales para el campo lejano para modelar de manera eficiente y precisa la propagación de ondas por carga móvil y la interacción suelo-estructura en medios geotécnicos heterogéneos semiinfinitos sin requerir cuadratura radial o límites artificiales.

Autores originales: Yanhui Zhong, Hao Hong, Bei Zhang, Quansheng Zang, Hussein Rappel, Stéphane P. A. Bordas

Publicado 2026-07-15

📖 7 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Yanhui Zhong, Hao Hong, Bei Zhang, Quansheng Zang, Hussein Rappel, Stéphane P. A. Bordas

Artículo original bajo licencia CC BY 4.0 (http://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 que intentas escuchar un tren pasando velozmente por un campo, pero el campo es en realidad un océano infinito de tierra que se extiende eternamente en todas las direcciones. Si intentaras construir un modelo computacional de todo este océano infinito para ver cómo tiembla el suelo, es probable que tu computadora se derritiera. Es como intentar contar cada grano de arena en una playa solo para ver cómo un paso hace que la arena ondule.

Este artículo presenta un nuevo y astuto truco llamado el Método de Elementos Infinitos con Traza NURBS en 2.5D (o NBIEM por sus siglas en inglés) para resolver este problema del "océano infinito" sin necesidad de una supercomputadora del tamaño de una ciudad.

El truco de magia: 2.5D

Primero, hablemos de la parte del "2.5D". Normalmente, para modelar un tren moviéndose por una vía, necesitas un modelo 3D completo. Pero si el suelo y la vía se ven iguales a lo largo de toda la línea (como un túnel largo y recto o una carretera recta), no necesitas modelar todo el mundo 3D. Solo necesitas modelar una sola sección transversal del suelo (la sección transversal) y luego usar un "hechizo matemático" (una transformada de Fourier) para imaginar cómo esa sección se repite mientras el tren pasa a toda velocidad. Esto convierte un enorme rompecabezas 3D en una serie de rompecabezas 2D más pequeños y manejables.

El problema: El borde infinito

Incluso con esta sección 2D, hay un inconveniente. El suelo no se detiene; continúa para siempre. En los modelos computacionales antiguos, los científicos simplemente cortaban el suelo a cierta distancia y fingían que terminaba ahí. Pero eso es como cortar el extremo de un trampolín; las ondas rebotarían en ese falso borde cortado y arruinarían los resultados.

Para solucionar esto, utilizaron "elementos infinitos". Piensa en ellos como extensiones especiales y mágicas conectadas al borde de tu modelo que saben cómo dejar que las ondas escapen hacia el vacío sin rebotar.

La nueva solución: NURBS-Trace

La gran innovación de los autores es cómo conectan el "campo cercano" (el suelo justo debajo del tren) con el "campo infinito" (el resto del mundo).

La forma antigua: Imagina que el campo cercano está hecho de piezas de Lego y el campo infinito está hecho de plastilina. Para conectarlos, tienes que aplastar la plastilina para que encaje con el Lego, lo que a menudo crea huecos o conexiones inestables.
La nueva forma (NBIEM): Los autores utilizan una técnica llamada Análisis Isogeométrico. En lugar de Lego y plastilina, utilizan el mismo material suave y flexible (curvas matemáticas llamadas NURBS) tanto para el campo cercano como para el campo infinito.

Lo llaman un método de "NURBS-trace". Es como tener una única cinta continua y sin costuras que fluye desde las vías del tren hacia el infinito. Debido a que la cinta es del mismo material en ambos lados, no necesitan forzar la unión. La conexión es perfecta, suave y matemáticamente exacta.

El atajo de "Sin Cuadratura"

Aquí está la parte realmente genial. Usualmente, cuando se calcula cómo se comportan estas ondas infinitas, las computadoras tienen que realizar millones de pasos matemáticos diminutos y repetitivos (llamados cuadratura numérica) para adivinar la respuesta. Es lento y puede volverse complicado.

Los autores encontraron una manera de escribir la respuesta para la parte infinita utilizando una fórmula de forma cerrada. Piensa en esto como: en lugar de contar cada uno de los pasos que da un corredor para llegar a la meta, encontraron una única ecuación que dice exactamente cuánto dura la carrera. Esto significa que su computadora no tiene que dividir el suelo infinito en trozos diminutos ni hacer millones de conjeturas. Simplemente introduce los números en la fórmula y, boom, la respuesta está ahí. Esto hace que el método sea increíblemente rápido.

Lo que probaron (y lo que no)

El equipo no solo imaginó esto; lo sometió a prueba en una serie de simulaciones:

La "Prueba de Velocidad": Simularon trenes moviéndose a diferentes velocidades:
- Sub-Rayleigh: Más lento que la velocidad de vibración natural del suelo (como un coche conduciendo por una carretera con baches).
- Super-shear: Más rápido que las ondas de cizalla pero más lento que las ondas de compresión (un poco como el inicio de un estallido sónico).
- Super-compresional: Más rápido que todas las ondas del suelo (un verdadero estallido sónico en el suelo).
- El Resultado: En todos estos casos, su nuevo método coincidió casi exactamente con las soluciones matemáticas "perfectas" conocidas, incluso para el esfuerzo y la deformación (cuánto se está comprimiendo o estirando el suelo).
La "Comparación": Compararon su nuevo método con otras dos formas populares de hacer esto:
- FEM/IEM: La forma estándar y más antigua.
- IGA/SBIGA: Una forma más compleja y nueva.
- El Resultado: Su nuevo método fue mucho más rápido que la forma estándar y mucho más rápido que la forma nueva y compleja (¡que tardó 168 segundos en resolver un problema que su método resolvió en 0.023 segundos!). También fue más preciso que la forma estándar, especialmente en frecuencias más altas.
Escenarios del Mundo Real: Simularon:
- Suelo Estratificado: Como un pastel con diferentes sabores de tierra.
- Túneles Enterrados: Como la Línea 1 del Metro de Zhengzhou. Comprobaron qué tan profundo debe estar el túnel para evitar que las vibraciones lleguen a la superficie. Encontraron que una vez que el túnel tiene unos 12 a 14 metros de profundidad, las vibraciones por encima de él disminuyen significativamente, y profundizar más ofrece beneficios cada vez menores.
- Capas Blandas: Probaron qué sucede si hay una capa de suelo suave y blanda bajo la carretera. Encontraron que tratar el suelo hasta una profundidad de aproximadamente 1.2 a 1.5 metros es lo que más ayuda; ir más profundo no ayuda mucho más.

Lo que explícitamente descartan

El artículo es muy claro sobre lo que este método no es:

No es un método que requiera conocer la respuesta "perfecta" de antemano para ajustar sus configuraciones. Crearon un flujo de trabajo "libre de referencia" que funciona sin necesidad de una solución analítica para contrastar.
No es un método que dependa de "Capas Perfectamente Acopladas" (PML), que son capas absorbentes artificiales gruesas que requieren una configuración cuidadosa. Su método utiliza los elementos infinitos conectados directamente al límite.
No es un método que sugiera que la versión "calibrada" (que utiliza matemáticas adicionales complejas para ajustar la atenuación) sea siempre necesaria. Encontraron que la versión más simple "all-S" (que utiliza solo un tipo de información de onda) funciona de maravilla para la mayoría de las necesidades de ingeniería, y la calibración compleja es solo un extra para casos específicos de baja frecuencia.

¿Qué tan seguros están?

Los autores están muy seguros de sus simulaciones. Han demostrado a través de sus propios modelos computacionales que:

El método funciona para diferentes velocidades de tren.
Es más rápido y preciso en comparación con los métodos existentes.
Maneja bien las formas complejas (como túneles curvos y suelos estratificados).

Sin embargo, son cuidadosos al decir que estas son simulaciones. Aún no han comparado sus resultados con mediciones del mundo real de una vía de tren real (aunque sugieren que ese sería un excelente siguiente paso). También señalan que sus configuraciones "calibradas" se derivaron de un benchmark específico y podrían necesitar ajustes si las condiciones del suelo cambian drásticamente.

La Conclusión

Este artículo presenta una forma nueva, súper rápida y fluida de modelar cómo tiembla el suelo cuando un tren o un camión pesado pasa por su lado. Al utilizar una cinta matemática continua y sin costuras para conectar el suelo cerca de la carretera con el suelo infinito lejos de ella, y al utilizar una "fórmula de atajo" para calcular la parte infinita, pueden resolver estos problemas mucho más rápido que antes. Es una herramienta que los ingenieros pueden usar para diseñar mejores vías y túneles, asegurando que las vibraciones no molesten a las personas que viven cerca, todo sin necesidad de una computadora que cueste más que el PIB de un pequeño país.

Resumen Técnico: Un Método de Elementos Infinitos de Traza NURBS 2.5D para la Propagación de Ondas por Cargas Móviles

Planteamiento del Problema
La vibración del suelo inducida por cargas móviles (p. ej., trenes de alta velocidad, tráfico pesado) es una preocupación crítica en la geotecnia de transporte. Predecir estas vibraciones requiere modelar la propagación de ondas en dominios terrestres semi-infinitos y no acotados. Aunque la respuesta es inherentemente tridimensional, los sistemas con geometría y propiedades de materiales invariantes a lo largo de la dirección de desplazamiento (eje longitudinal $z$ ) pueden analizarse eficientemente mediante una formulación 2.5D. Este enfoque reduce el problema a una secuencia de problemas de sección transversal 2D en el dominio de la frecuencia-número de onda mediante la transformación de Fourier. Sin embargo, modelar con precisión el dominio exterior no acotado, satisfaciendo las condiciones de radiación y manteniendo la consistencia geométrica con interfaces complejas (túneles curvos, suelos estratificados), sigue siendo un desafío. Los métodos existentes, como el acoplamiento de Elementos Finitos/Infinitos (FEM/IEM), los Métodos de Elementos de Contorno (BEM) y los Métodos de Elementos Finitos de Frontera Escalada (SBFEM), a menudo enfrentan compensaciones entre el costo computacional, la flexibilidad geométrica y la necesidad de soluciones fundamentales o capas absorbentes artificiales.

Metodología
El artículo propone un Método de Elementos Infinitos de Traza NURBS (NBIEM) 2.5D No Uniforme. Esta formulación acopla un dominio de campo cercano acotado con un dominio exterior semi-infinito utilizando un marco isogeométrico unificado.

Reducción 2.5D: Las ecuaciones de movimiento de los medios viscoelásticos lineales se transforman al dominio de la frecuencia-número de onda. La carga móvil se trata como una excitación armónica donde la frecuencia temporal $\omega$ y el número de onda longitudinal $k$ están acoplados mediante la velocidad de la carga $c_{load}$ ( $k = (\omega - \omega_0)/c_{load}$ ).
Campo Cercano Acotado (IGA): La sección transversal del campo cercano se discretiza mediante Análisis Isogeométrico (IGA). Tanto la geometría como el campo de desplazamientos se representan mediante las mismas funciones base NURBS, lo que permite la representación exacta de fronteras curvas y continuidad de alto orden.
Exterior Semi-Infinito (Elementos Infinitos de Traza NURBS): El dominio exterior se modela utilizando elementos infinitos construidos como productos tensoriales de:
- La traza NURBS de la frontera (compartida con la discretización del campo cercano).
- Funciones radiales admisibles ( $P_q$ ) que representan ondas salientes o evanescentes (p. ej., decaimiento exponencial).
- Crucialmente, el campo cercano y el campo lejano comparten la misma traza NURBS discreta y los grados de libertad de los puntos de control de la frontera. Esto impone la continuidad del desplazamiento directamente sin necesidad de proyecciones, variables de mortero o operaciones de transferencia de interfaz.
Integración en Forma Cerrada: Para funciones radiales exponenciales seleccionadas, los momentos radiales requeridos para las matrices de rigidez y masa del campo lejano se evalúan en forma cerrada. Esto elimina la necesidad de cortes radiales finitos o integración numérica radial, asegurando la satisfacción exacta de la condición de radiación dentro del espacio funcional elegido.
Manejo de la Condición de Radiación: El método emplea una realización "all-S" donde se asignan funciones radiales informadas por ondas S a todos los segmentos de la frontera. Los parámetros de decaimiento se determinan mediante las propiedades del material exterior y las escalas geométricas. Para garantizar la robustez sin referencias analíticas, los autores introducen dos indicadores de sensibilidad: $I_\beta$ (sensibilidad a los parámetros de decaimiento radial) e $I_D$ (sensibilidad a la ubicación de la frontera).

Contribuciones Clave
El artículo realiza tres contribuciones principales:

Formulación: Desarrollo de una formulación de elementos infinitos de traza NURBS 2.5D de tres componentes. Acopla de manera única el dominio IGA acotado y el exterior semi-infinito a través de una traza NURBS discreta idéntica, preservando el acoplamiento completo de deformación e inercia dependiente de la frecuencia y el número de onda longitudinal entre los tres componentes de desplazamiento.
Eficiencia Analítica: Evaluación de las contribuciones de rigidez y masa del campo lejano mediante momentos radiales en forma cerrada para funciones radiales exponenciales seleccionadas. Esto evita la truncación radial finita y la integración numérica en la dirección radial.
Caracterización Sistemática: Verificación y evaluación exhaustiva del método, incluyendo:
- Validación contra soluciones de semiespacio en forma cerrada a través de regímenes sub-Rayleigh, super-corte (super-shear) y super-compresional.
- Desarrollo de indicadores de sensibilidad sin referencia ( $I_\beta$ e $I_D$ ) para evaluar la adecuación de la frontera.
- Comparación directa con FEM/IEM clásico e IGA/SBFEM (IGA/SBIGA) en cuanto a precisión, fidelidad de fase y costo computacional.

Resultos

Precisión: El NBIEM reproduce las funciones de respuesta en frecuencia de desplazamiento y esfuerzo con alta precisión, coincidiendo con las soluciones analíticas tanto para cargas móviles sub-Rayleigh como super-compresionales. Los resultados de esfuerzo (campos dependientes de la derivada) confirman la capacidad del método para capturar gradientes, beneficiándose de la mayor continuidad de las bases NURBS.
Comparación con Alternativas:
- vs. FEM/IEM: El NBIEM logra una precisión de desplazamiento y fase comparable o superior con significativamente menos grados de libertad. Es aproximadamente un orden de magnitud más rápido en tiempo de ejecución para mallas moderadas a grandes debido a la eliminación de la construcción de la rigidez dinámica de frontera escalada.
- vs. IGA/SBIGA: El NBIEM logra una precisión similar al IGA/SBIGA en el rango de frecuencias intermedias a altas, pero con un costo computacional sustancialmente menor. A diferencia del IGA/SBIDA, no requiere resolver un sistema de rigidez dinámica de frontera escalada, evitando así los altos costos de ensamblaje asociados.
Aplicaciones de Ingeniería: El método modela con éxito:
- Suelo Estratificado: Manejo de interfaces de materiales y contrastes de rigidez.
- Sistemas de Vía-Subrasante: Análisis de configuraciones de terraplén frente a corte con geometrías de múltiples parches, revelando la amplificación dependiente de la velocidad en interfaces profundas.
- Estructuras Enterradas: Simulación de túneles de metro con variaciones en la profundidad de cobertura, demostrando cómo el aumento de la profundidad de enterramiento atenúa la transmisión de vibraciones al suelo sobre la clave del túnel.
Sensibilidad: Se encontró que la realización por defecto "all-S" es robusta. Los indicadores de sensibilidad identificaron eficazmente los casos que requerían la ampliación de la frontera (p. ej., bajas frecuencias o velocidades super-shear), asegurando la confiabilidad de la solución sin necesidad de comparaciones analíticas.

Significado y Reivindicaciones
El artículo sostiene que el NBIEM proporciona un marco geométricamente consistente y computacionalmente eficiente para la propagación de ondas 2.5D y la interacción suelo-estructura en dominios semi-infinitos. Su importancia radica en:

Consistencia Geométrica: Al utilizar trazas NURBS tanto para el campo cercano como para el lejano, maneja naturalmente interfaces curvas complejas y configuraciones de múltiples parches, evitando los errores de aproximación geométrica comunes en el FEM de bajo orden.
Eficiencia Computacional: Ofrece un equilibrio favorable entre precisión y costo, particularmente para análisis repetidos en el dominio de la frecuencia, al evitar la costosa construcción de la rigidez dinámica de SBFEM y las dependencias de la solución fundamental de BEM.
Robustez: Los indicadores de sensibilidad propuestos permiten un cribado automatizado y sin necesidad de referencias de la configuración de la frontera, haciendo que el método sea práctico para aplicaciones de ingeniería donde no existen soluciones analíticas.

Los autores posicionan al NBIEM no como un reemplazo para todos los métodos existentes, sino como una alternativa distinta y eficiente que cierra la brecha entre la flexibilidad geométrica del análisis isogeométrico y el tratamiento de dominios no acotados de los elementos infinitos, diseñado específicamente para problemas de carga móvil en la ingeniería geotécnica.