Nonlinear Incompressible Shear Wave Models in Hyperelasticity and Viscoelasticity Frameworks, with Applications to Love Waves

Este artículo presenta ecuaciones generales para ondas de corte no lineales en materiales hiper-viscoelásticos incompresibles, aplicándolas al análisis de ondas de Love en interfaces y demostrando mediante simulaciones numéricas que, aunque la velocidad de onda variable cumple la condición de existencia lineal, tiende a acercarse a la velocidad de onda del material más rápido a largo plazo.

Lo esencial

¡Hola! Imagina que el mundo está hecho de materiales elásticos, como una mezcla gigante de goma, gelatina y roca. Cuando algo golpea estos materiales (como un terremoto o una explosión), se generan ondas que viajan a través de ellos.

Este artículo científico es como un manual de instrucciones avanzado para entender cómo se comportan esas ondas cuando los materiales no son simples y perfectos, sino que tienen "personalidad": se estiran mucho, se deforman y tienen una memoria viscosa (como la miel).

Aquí tienes la explicación de la investigación, traducida a un lenguaje sencillo y con analogías creativas:

1. El Problema: Ondas en "Goma Elástica" vs. "Miel"

En la física clásica, imaginamos que los materiales se comportan como resortes perfectos (si los estiras un poco, vuelven igual). Pero en la realidad, materiales como la piel, los tejidos biológicos o incluso las capas de la Tierra se comportan de forma más compleja:

• Hiperelasticidad: Imagina una banda elástica muy gruesa. Si la estiras un poco, se comporta normal. Pero si la estiras mucho, se pone dura y cuesta más estirarla. Eso es lo que estudian aquí: ondas que viajan por materiales que cambian de "dureza" según cuánto se deforman.
• Viscoelasticidad: Ahora imagina que esa banda elástica está hecha de miel. Cuando la mueves, no solo rebota, sino que se queda un poco pegajosa y pierde energía en forma de calor. El modelo incluye esta "pegajosidad" para que las ondas no viajen para siempre, sino que se apaguen con el tiempo, como en la vida real.

2. El Experimento: La "Explosión de Galleta"

Para probar sus teorías, los autores hicieron simulaciones por computadora. Imagina que tienes una capa de gelatina (la corteza terrestre) flotando sobre un océano de roca más dura (el manto).

• La Prueba: En el centro de esta gelatina, hacen una pequeña "explosión" (como lanzar una piedra a un estanque, pero en 3D).
• El Objetivo: Ver cómo se mueve la onda. ¿Se queda atrapada en la capa superior? ¿Viaja por la frontera entre las dos capas? ¿Cuánto tarda en llegar a la superficie?

3. Los Descubrimientos Clave

A. La Regla de la "Carrera de Velocidad"

En la física lineal (la versión simple), existe una regla estricta para que las ondas de Love (un tipo de onda que viaja por la superficie) existan: deben ir más rápido que el material de arriba pero más lento que el de abajo.

• La Analogía: Imagina un coche de carreras (la onda) en una pista. Si la pista de arriba es de barro (lento) y la de abajo es de hielo (rápido), el coche puede "surfear" en la frontera entre ambos.
• El Hallazgo: Los autores descubrieron que incluso cuando el material se deforma mucho (comportamiento no lineal), la onda sigue respetando esta regla. Aunque la velocidad de la onda cambia mientras viaja, siempre se mantiene dentro de ese rango permitido. Es como si la onda tuviera un "GPS" que le impide ir demasiado rápido o demasiado lento.

B. El Efecto de la "Memoria" (Viscosidad)

Cuando añadieron la parte viscosa (la miel), las ondas se comportaron de forma más realista:

• Se dispersaron más rápido.
• Perdieron energía.
• Analogía: Si gritas en una habitación vacía (sin viscosidad), el eco dura mucho. Si gritas en una habitación llena de almohadas y miel (con viscosidad), el sonido se ahoga rápido. El modelo muestra cómo la Tierra "absorbe" el sonido de un terremoto, evitando que las ondas viajen eternamente.

C. La Simetría y las Soluciones Mágicas

Los autores usaron matemáticas avanzadas (simetrías) para encontrar soluciones exactas en una dimensión (como una onda en una cuerda).

• Encontraron fórmulas que describen ondas que viajan con una forma muy específica, aunque matemáticamente son un poco "locas" (se hacen infinitas en un punto).
• Sin embargo, cuando compararon estas fórmulas con sus simulaciones por computadora, vieron que las ondas reales se parecen mucho a esas fórmulas mágicas, especialmente al principio. Es como si la naturaleza intentara seguir una receta matemática perfecta, pero la "suciedad" de la realidad (la viscosidad) la suaviza un poco.

4. ¿Por qué importa esto?

Este trabajo es como actualizar el mapa de navegación de los sismólogos.

• Antes: Usaban mapas simples que asumiían que la Tierra era un resorte perfecto.
• Ahora: Tienen un mapa más detallado que sabe que la Tierra es como una goma elástica que se estira y una miel que se calienta.
• Aplicación: Esto ayuda a entender mejor cómo se mueven las ondas sísmicas cerca del epicentro de un terremoto (donde los movimientos son grandes) y cómo afectan a los edificios. También es útil para estudiar tejidos biológicos o materiales compuestos.

En Resumen

Los autores crearon una fórmula maestra que describe cómo se mueven las ondas en materiales que se deforman mucho y tienen "pegajosidad". Demostraron que, incluso en el caos de una deformación grande, las ondas siguen reglas de velocidad muy estrictas y que la viscosidad es clave para entender por qué las ondas se apagan. Es una mezcla de matemáticas puras, simulaciones por computadora y una visión más realista de cómo vibra nuestro planeta.

Resumen técnico

A continuación presento un resumen técnico detallado del artículo de investigación "Modelos de Ondas de Corte No Lineales Incompresibles en Marcos de Hiperelasticidad y Viscoelasticidad, con Aplicaciones a Ondas de Love", escrito por Shawn Samuel Carl McAdam, Samuel Opoku Agyemang y Alexei Cheviakov.

1. Problema y Contexto

El estudio de las ondas de corte (ondas S) que se propagan a través de materiales elásticos es fundamental en sismología, geología y ciencia de materiales. Tradicionalmente, los modelos se basan en la teoría de la elasticidad lineal (Ley de Hooke), que asume deformaciones pequeñas y relaciones lineales entre tensión y deformación. Sin embargo, muchos materiales reales (polímeros, tejidos biológicos, formaciones geológicas) exhiben comportamientos no lineales bajo grandes deformaciones.

El problema central abordado en este trabajo es la derivación y análisis de ecuaciones gobernantes para ondas de corte de amplitud finita (no lineales) que se propagan en materiales hiperelásticos incompresibles y, posteriormente, en materiales hiper-viscoelásticos. Específicamente, el estudio se centra en la generalización de las ondas de Love (ondas superficiales polarizadas horizontalmente) en un modelo de "Tierra plana" de dos capas, donde se investigan los efectos de la no linealidad y la disipación de energía en la velocidad de fase y la estabilidad de la onda.

2. Metodología

El enfoque metodológico se divide en tres etapas principales: derivación teórica, análisis de simetría y simulación numérica.

A. Marco Teórico y Derivación de Ecuaciones

• Hiperelasticidad: Se utiliza el marco de mecánica de medios continuos en coordenadas materiales (Lagrangianas). Se asume un material incompresible, isotrópico y homogéneo.
• Energía de Deformación: Se define una función de densidad de energía de deformación $W^H$. El trabajo se especializa en materiales de tipo neo-Hookeano generalizado, donde $W^H$ depende solo del primer invariante ($I_1$) del tensor de deformación de Cauchy-Green, ignorando el segundo invariante ($I_2$) para simplificar el sistema y garantizar la compatibilidad de la presión hidrostática.
• Modelo Constitutivo: Se aplica específicamente el modelo cúbico de Yeoh, que permite representar curvas tensión-deformación no lineales mediante un polinomio cúbico en términos de $(I_1 - 3)$.
• Viscoelasticidad: Se introduce un potencial de disipación $W^V$ (función pseudo-energía) dependiente de la tasa de cambio del tensor de deformación, modelando efectos viscosos mediante un coeficiente $\eta$.
• Derivación: A partir de las leyes de conservación (masa, momento lineal y angular), se derivan ecuaciones de onda no lineales parciales (PDEs) para el desplazamiento transversal $v$.
  • Para el caso hiperelástico, se obtiene una ecuación de onda con términos polinómicos cúbicos y quínticos.
  • Para el caso hiper-viscoelástico, se añade un término de viscosidad que involucra derivadas mixtas (ej. $u_{xxt}$).

B. Análisis de Simetría (Lie)

Para el caso unidimensional (1D) de las ecuaciones viscoelásticas, se emplea el método de simetrías de Lie. Esto permite:

• Identificar generadores de simetría puntual.
• Derivar soluciones exactas invariantes (soluciones de onda viajera).
• Analizar la estructura de las soluciones y su comportamiento asintótico.

C. Simulación Numérica

Se utilizan simulaciones numéricas en 2+1 dimensiones (dos espaciales, una temporal) mediante el método de líneas:

• Discretización: Se utiliza diferencias finitas de segundo orden en una malla no uniforme (más densa cerca de la fuente y la interfaz).
• Integración Temporal: Se emplea el solucionador ode23 de MATLAB (Runge-Kutta de orden 2 y 3) para resolver el sistema de EDOs resultante.
• Condiciones Iniciales: Se simula una "explosión gaussiana" súbita en diferentes ubicaciones (capa superior o semiespacio inferior) para generar ondas.
• Escenarios: Se comparan tres casos:
  1. Ondas de Love lineales (referencia).
  2. Ondas de Love hiperelásticas (modelo Yeoh cúbico, sin viscosidad).
  3. Ondas de Love hiper-viscoelásticas (modelo Yeoh + viscosidad).

3. Contribuciones Clave

1. Derivación de Nuevas Ecuaciones de Onda: Se presentan ecuaciones de onda no lineales generales para desplazamientos de corte en materiales hiperelásticos incompresibles con una función de energía de deformación arbitraria. Específicamente, se derivan las ecuaciones explícitas para el modelo cúbico de Yeoh, que incluyen términos de orden cúbico y quíntico.
2. Inclusión de Viscoelasticidad No Lineal: Se generaliza el modelo para incluir efectos viscoelásticos mediante un potencial de disipación, resultando en ecuaciones que contienen términos de dispersión con derivadas mixtas ($u_{xxt}$), lo cual es crucial para modelar la atenuación realista en materiales geológicos y biológicos.
3. Validación de la Condición de Existencia No Lineal: Se demuestra numéricamente que, incluso en el régimen totalmente no lineal, la velocidad de fase de las ondas de Love interfaciales y superficiales satisface la condición clásica de existencia lineal: $c_1 < |v| < c_2$ (donde $c_1$ y $c_2$ son las velocidades de onda de los materiales adyacentes).
4. Análisis de Comportamiento Temporal: Se observa que, aunque las velocidades de fase pueden variar inicialmente, tienden asintóticamente hacia la mayor de las dos velocidades de onda del material ($c_1$ o $c_2$) a largo plazo.
5. Corrección de Errores y Aclaraciones: El artículo incluye un "Addendum" que corrige una suposición implícita en la versión original sobre la dependencia de $I_2$ en la energía de deformación, restringiendo explícitamente el modelo a materiales neo-Hookeanos generalizados ($W^H_2 = 0$) para garantizar la consistencia matemática de la presión.

4. Resultados Principales

• Comportamiento de Ondas Lineales vs. No Lineales: Las simulaciones muestran que el comportamiento cualitativo de las ondas no lineales es similar al lineal en términos de propagación y reflexión en la interfaz. Sin embargo, la no linealidad introduce oscilaciones de alta frecuencia y modula la velocidad de fase localmente.
• Efecto de la Viscosidad: La inclusión de términos viscoelásticos provoca una dispersión más rápida de la energía y una atenuación significativa de la amplitud de la onda en comparación con los casos puramente elásticos. Las ondas se vuelven más predecibles y estables a medida que aumenta el coeficiente de viscosidad $\eta$.
• Velocidad de Fase: En todos los casos simulados (lineal, hiperelástico, hiper-viscoelástico), la velocidad de la onda a lo largo de la interfaz y la superficie tiende a converger hacia el máximo de las velocidades de onda de los dos medios ($\max(c_1, c_2)$) a medida que avanza el tiempo.
• Soluciones Exactas (1D): El análisis de simetría en 1D proporcionó soluciones de onda viajera exactas. Aunque estas soluciones analíticas son no acotadas, las simulaciones numéricas con condiciones iniciales físicas (explosión gaussiana) muestran un comportamiento cualitativo muy cercano a estas soluciones, confirmando la validez del modelo.
• Ruptura Numérica: En el apéndice, se reporta el descubrimiento de soluciones que experimentan "ruptura numérica" (formación de oscilaciones espurias y cambios bruscos en la pendiente) en configuraciones radiales simétricas, lo que sugiere la necesidad de esquemas numéricos de alta resolución para capturar fenómenos de choque en ondas de gran amplitud.

5. Significado e Impacto

Este trabajo es significativo por varias razones:

• Avance en Sismología Teórica: Proporciona un marco matemático más robusto para modelar ondas sísmicas superficiales (Love) que ocurren cerca del foco de un terremoto, donde las deformaciones pueden ser grandes y los efectos no lineales no despreciables, a diferencia de los modelos lineales estándar que asumen deformaciones infinitesimales.
• Aplicabilidad Multidisciplinaria: El modelo es aplicable no solo a la geofísica, sino también a la investigación de tejidos biológicos y materiales compuestos, donde la viscoelasticidad y la no linealidad son inherentes.
• Fundamento para Futuras Investigaciones: Establece las ecuaciones base para estudiar ondas no lineales en medios anisotrópicos y para otras ondas superficiales (Rayleigh, Scholte, Stoneley), abriendo la puerta a modelos más complejos que incluyan compresibilidad y anisotropía.
• Validación Numérica: Demuestra que las condiciones de existencia de las ondas de Love derivadas de la teoría lineal se mantienen válidas incluso en regímenes no lineales fuertes, lo que refuerza la utilidad de los modelos lineales como aproximaciones de primer orden, mientras que los modelos no lineales son necesarios para capturar la dinámica transitoria y la atenuación precisa.

En resumen, el artículo conecta la teoría constitutiva avanzada de materiales con la dinámica de ondas no lineales, ofreciendo tanto herramientas analíticas (simetrías) como numéricas para comprender mejor la propagación de ondas de corte en medios complejos.