Caveats on formulating finite elasto-plasticity in curvilinear coordinates

Este artículo presenta una metodología práctica para la implementación numérica de la elastoplasticidad finita en coordenadas curvilíneas mediante cambios de base explícitos, aclarando el tratamiento de términos cinemáticos y tensionales específicos que surgen en deformaciones grandes y efectos anelásticos bajo simetría axial.

Lo esencial

Imagina que eres un arquitecto o un ingeniero que debe diseñar un puente, un túnel o incluso simular cómo se comporta un globo cuando lo inflas. Para hacer esto en una computadora, usamos un método llamado "Elementos Finitos". Básicamente, dividimos el objeto en millones de pequeños bloques (como un rompecabezas gigante) y calculamos cómo se mueve cada uno.

Normalmente, estos cálculos se hacen usando una cuadrícula simple y recta, como las líneas de un cuaderno de notas (coordenadas cartesianas). Pero, ¿qué pasa si el objeto que estudias es un cilindro, un tanque de agua o un agujero en el suelo? Usar una cuadrícula recta para algo redondo es como intentar medir la circunferencia de una pizza usando solo una regla recta: es posible, pero muy incómodo y propenso a errores.

Aquí es donde entra este artículo. Los autores nos dicen: "¡Oye! Si usamos coordenadas curvas (como las de un cilindro) para simular cosas que se deforman mucho (como un plástico estirado o un suelo que se hunde), hay trampas ocultas que pueden arruinar tus cálculos si no tienes cuidado."

Aquí te explico los conceptos clave con analogías sencillas:

1. El problema de la "Regla Curva" (Coordenadas Curvilíneas)

Imagina que tienes un mapa del mundo. Si dibujas una cuadrícula sobre un globo terráqueo, las líneas no son rectas; se curvan.

• En el mundo plano (Cartesiano): Si estiras un trozo de goma, la matemática es sencilla.
• En el mundo curvo (Cilíndrico): Si estiras un tubo, la matemática se vuelve extraña. Las distancias cambian no solo porque el material se estira, sino porque la "regla" que usas para medir también cambia de forma.

El papel advierte que muchos programas de computadora simplemente copian las fórmulas del mundo plano y las pegan en el mundo curvo. Esto es un error. Es como intentar usar las instrucciones de un coche de carreras para pilotar un barco; aunque ambos tienen ruedas (o hélices), las reglas de la física son diferentes.

2. El "Deslizador" (El Shifter)

Este es el concepto más importante y curioso del artículo.
Imagina que tienes un mapa antiguo (el estado original de tu objeto) y un mapa actualizado (el estado deformado).

• En un mundo plano, puedes simplemente tomar un punto del mapa antiguo y moverlo al nuevo.
• En un mundo curvo, los "puntos de referencia" (las líneas de latitud y longitud) giran y cambian de tamaño.

Los autores introducen algo llamado el "Shifter" (Deslizador).

• La analogía: Imagina que tienes un grupo de bailarines en un escenario circular. Si todos dan un paso hacia el centro, no solo se mueven; también cambian su orientación relativa. El "Shifter" es como un traductor de baile que le dice a la computadora: "Oye, ese bailarín no solo se movió 1 metro, sino que también giró porque el escenario es redondo. Ajusta el cálculo".
• Si ignoras al "Shifter", la computadora piensa que el material se ha deformado de una manera que en realidad no ha ocurrido, y los resultados serán basura.

3. El "Volumen Fantasma" (El Jacobiano)

Cuando estiras un chicle, su volumen cambia. En matemáticas, usamos algo llamado "Jacobiano" para medir cuánto cambia el volumen.

• En coordenadas rectas, es fácil: solo multiplicas los estiramientos.
• En coordenadas curvas, hay un "factor de corrección" oculto. Es como si el espacio mismo se encogiera o se expandiera solo por el hecho de estar en una curva.
• El artículo explica cómo calcular este volumen real, separando lo que es elasticidad (el material vuelve a su forma, como un resorte) de lo que es plasticidad (el material se queda deformado, como arcilla). Si mezclas estos dos conceptos sin cuidado en un entorno curvo, la computadora puede pensar que el material se ha creado o destruido materia de la nada.

4. ¿Por qué nos importa esto? (La analogía del Túnel)

Piensa en la construcción de un túnel o la extracción de petróleo.

• Si usas la fórmula "plana" incorrecta en un entorno "curvo", tu simulación podría decirte que el túnel es seguro cuando en realidad se va a derrumbar.
• O podría decirte que necesitas más acero del necesario, desperdiciando millones de dólares.

¿Qué logran los autores?

Ellos han escrito un manual de instrucciones paso a paso para los programadores de simulaciones. Han creado una "caja de herramientas" matemática que:

1. Usa el "Shifter" para traducir correctamente entre el mapa antiguo y el nuevo.
2. Calcula el volumen real considerando la curvatura.
3. Separa la parte elástica de la plástica para que la computadora no se confunda.

En resumen:
Este artículo es como un aviso de seguridad para los ingenieros que usan simulaciones por computadora. Les dice: "No asumas que lo que funciona en una línea recta funciona en un círculo. Si quieres simular cosas reales (como túneles, pozos o crecimiento de tejidos) que sufren grandes deformaciones, necesitas usar estas nuevas reglas de traducción matemática, o tus resultados serán tan falsos como un mapa plano de la Tierra".

Gracias a este trabajo, los ingenieros pueden hacer simulaciones más rápidas (usando solo una "rebanada" del objeto en 2D en lugar de todo el objeto en 3D) y, lo más importante, más precisas y seguras.

Resumen técnico

Aquí tienes un resumen técnico detallado del artículo en español, estructurado según los puntos solicitados:

Título: Advertencias sobre la formulación de la elastoplasticidad finita en coordenadas curvilíneas

1. Planteamiento del Problema

El análisis de mecánica de sólidos a menudo utiliza aproximaciones bidimensionales (como la simetría axial) para reducir los costos computacionales de métodos numéricos como el Método de Elementos Finitos (MEF). Sin embargo, existe una brecha crítica entre la formulación teórica de la mecánica de medios continuos, que se basa en tensores invariantes de marco, y su implementación numérica, que depende de representaciones matriciales en bases seleccionadas por el usuario.

El problema central abordado en este trabajo es que, al extender formulaciones cartesianas (ortogonales y ortonormales) a sistemas de coordenadas curvilíneas (como cilíndricas o esféricas) para problemas de grandes deformaciones (estrain finito) y elastoplasticidad, surgen errores fundamentales si no se tratan cuidadosamente ciertos términos geométricos. Específicamente:

• La definición del gradiente de deformación ($F$) y su determinante (el Jacobiano, $J$) no es trivial en coordenadas curvilíneas.
• La inclusión de efectos inelásticos (plasticidad) complica la distinción entre componentes elásticos y anelásticos del gradiente de deformación y del Jacobiano.
• La linealización consistente necesaria para la convergencia del método de Newton-Raphson en MEF se vuelve más compleja debido a la dependencia de las métricas y los símbolos de Christoffel.

El artículo utiliza el problema de la expansión de una cavidad en un cilindro hueco como caso de estudio inicial para exponer estas dificultades.

2. Metodología

Los autores adoptan un enfoque práctico en lugar de puramente geométrico (evitando el uso de variedades riemannianas complejas). En su lugar, trabajan dentro de representaciones cartesianas estándar y utilizan cambios de base explícitos para obtener las formas curvilíneas requeridas.

Los pasos metodológicos clave incluyen:

• Definición de Bases: Diferenciación clara entre bases cartesianas ($e_J, e_j$) y bases curvilíneas covariantes/contravariantes ($g_A, g_a$) en las configuraciones de referencia y actual.
• Gradiente de Deformación y Desplazadores: Se define rigurosamente el gradiente de deformación ($F$) como un tensor de dos puntos. Se introduce el concepto del tensor desplazador (shifter), $S$, que mapea las bases de la configuración original a la actual. En coordenadas cartesianas, este tensor es trivial (identidad), pero en curvilíneas tiene componentes no nulas que afectan los cálculos.
• Descomposición Multiplicativa: Se aplica la descomposición multiplicativa del gradiente de deformación ($F = F_e \cdot F_p$) para separar la respuesta elástica y plástica. Se define un Jacobian elástico ($J_e$) y un Jacobian plástico ($J_p$) para rastrear cambios de volumen reversibles e irreversibles.
• Medidas de Tensión y Equilibrio: Se derivan las formas fuertes y débiles del equilibrio de la cantidad de movimiento lineal en coordenadas curvilíneas, utilizando medidas de tensión adecuadas (como el tensor $\zeta = J_e \sigma$) que garantizan la satisfacción de la desigualdad de Clausius-Duhem y preservan una configuración intermedia libre de tensiones.
• Linealización Consistente: Se proporciona un procedimiento paso a paso para la linealización de las ecuaciones de equilibrio, incorporando derivadas covariantes y símbolos de Christoffel, esenciales para la convergencia cuadrática en el algoritmo de Newton-Raphson.
• Implementación en MEF: Se particularizan las formulaciones para problemas axisimétricos, reduciendo el problema 3D a 2D (plano R-Z) pero manteniendo la integridad de los componentes fuera del plano a través de los factores métricos.

3. Contribuciones Clave

• Clarificación de Términos Ocultos: El artículo demuestra que términos que son triviales o implícitos en coordenadas cartesianas (como el determinante del Jacobiano y los componentes fuera del plano del gradiente de deformación) requieren tratamientos explícitos y correctos en coordenadas curvilíneas para evitar errores fundamentales en el MEF.
• Tratamiento del Tensor Desplazador (Shifter): Se destaca el papel crucial del tensor desplazador en la definición correcta de los desplazamientos y el gradiente de deformación en sistemas no ortogonales, mostrando cómo su omisión o mal uso lleva a resultados incorrectos.
• Formulación de Jacobianos Elástico y Plástico: Se propone una metodología robusta para calcular y linealizar los Jacobianos elásticos y plásticos por separado en el contexto de coordenadas curvilíneas, utilizando la relación con la deformación logarítmica elástica para simplificar la derivación.
• Algoritmos de Implementación: Se presentan flujos de trabajo algorítmicos detallados (Algoritmos 1 y 2) para formulaciones Lagrangianas Actualizadas (UL) y Totales (TL) en problemas axisimétricos, especificando cómo actualizar métricas, símbolos de Christoffel y gradientes de funciones de forma.
• Validación Numérica: Se valida la formulación propuesta comparándola con un análisis 3D completo de un cilindro de pared gruesa sometido a grandes deformaciones y plasticidad.

4. Resultados

• Validación de la Formulación Axisimétrica: Los resultados numéricos muestran una coincidencia casi perfecta entre la formulación axisimétrica propuesta y el análisis 3D completo. Se compararon invariantes de tensión (presión de Cauchy, tensión desviadora) y Jacobianos (elástico y plástico) en puntos críticos del cilindro.
• Eficiencia Computacional: La formulación axisimétrica logra resultados tridimensionales precisos con una discretización bidimensional, reduciendo significativamente el costo computacional.
• Convergencia del Algoritmo: El análisis de la convergencia del error en el método de Newton-Raphson demuestra que la linealización consistente propuesta es robusta. Se observó que el número de iteraciones necesarias para la convergencia es muy similar entre las formulaciones 3D y 2D axisimétricas (promedio de ~3.5 iteraciones).
• Comportamiento del Material: El ejemplo de simulación (cilindro de arcilla) mostró correctamente la transición a la plasticidad y la evolución de las tensiones, confirmando que la formulación maneja adecuadamente la dependencia de la tensión de Cauchy respecto a la parte plástica del Jacobiano.

5. Significado e Impacto

Este trabajo es fundamental para el desarrollo de códigos de elementos finitos robustos para problemas de grandes deformaciones en geometrías axisimétricas o curvilíneas.

• Prevención de Errores: Advierte a los investigadores y desarrolladores de software que la simple extensión de fórmulas cartesianas a coordenadas curvilíneas es insuficiente y puede llevar a inconsistencias graves (por ejemplo, entre formulaciones Lagrangianas Actualizadas y Totales).
• Marco Práctico: Proporciona una guía paso a paso y explícita para implementar estas correcciones geométricas sin necesidad de recurrir a marcos matemáticos abstractos complejos, haciendo la teoría accesible para la ingeniería computacional.
• Aplicabilidad Amplia: Las conclusiones son relevantes para diversas áreas de la ingeniería, incluyendo geotecnia (expansión de cavidades, pruebas de penetración), metalurgia (forjado, extrusión) y biomecánica, donde la simetría axial y las grandes deformaciones son comunes.

En resumen, el artículo cierra la brecha entre la teoría tensorial abstracta y la implementación práctica en MEF, asegurando que las simulaciones de elastoplasticidad finita en coordenadas curvilíneas sean matemáticamente consistentes y numéricamente precisas.