A Projection Method for an Elasto-plasticity Model with Linear Kinematic Hardening

1. Planteamiento del Problema

El artículo aborda un sistema dinámico de elastoplasticidad que incorpora dos características físicas clave:

1. Viscosidad de Kelvin-Voigt: Introduce un término viscoso en la ecuación de movimiento para regularizar el problema y permitir el control espacial de la velocidad.
2. Endurecimiento Cinemático Lineal (Tipo Melan-Prager): Modela el fenómeno de endurecimiento por deformación donde el conjunto de restricciones (superficie de fluencia) se traslada en el espacio de tensiones debido a una variable interna llamada tensión de retroceso (backstress, $\alpha$).

El Modelo Matemático:
El sistema acopla la ecuación de movimiento con una ley de flujo plástico y una ecuación de evolución para la variable interna.

• Variables: Velocidad $v$, tensión $\sigma$, y variable interna $\alpha$.
• Ecuación de Movimiento: Incluye inercia, viscosidad ($\eta$), y divergencia de la tensión.
• Ley de Flujo: Se formula mediante una inclusión variacional. La tensión $\sigma(t)$ debe pertenecer a un conjunto convexo cerrado $K(t; \alpha)$ que depende del tiempo y de la variable interna $\alpha$.
• Conjunto de Restricciones: $K(t; \alpha) = \tilde{K}(t) + \alpha(t)$, donde $\tilde{K}(t)$ está definido por un criterio de von Mises con un límite de fluencia $g(t, x)$ que puede depender del tiempo y del espacio.
• Estructa Cuasi-Variacional: A diferencia de los modelos donde el conjunto de restricciones es fijo o dado externamente, aquí el conjunto $K$ depende de la solución desconocida $\alpha$, lo que introduce una estructura cuasi-variacional compleja.

El objetivo es probar la existencia y unicidad de soluciones débiles para este sistema bajo condiciones generales sobre el límite de fluencia $g$ (sin asumir continuidad espacial estricta ni cotas inferiores estrictamente positivas).

2. Metodología Propuesta

Los autores emplean el Método de Rothe (discretización temporal) combinado con un esquema de proyección para construir soluciones.

Esquema de Discretización (P):
Para cada paso de tiempo $n$, el problema se divide en tres sub-pasos secuenciales:

1. Paso de Ensayo (Trial Step): Se resuelve un subproblema lineal visco-elástico para obtener una tensión de ensayo $\sigma^*_n$ y una velocidad $v_n$, ignorando temporalmente la restricción de plasticidad.
2. Actualización de la Tensión de Retroceso: Se actualiza la variable interna $\alpha_n$ utilizando la tasa de deformación plástica aproximada $\xi_n$.
3. Paso de Proyección: Se proyecta la tensión de ensayo desplazada ($\sigma^*_n - \alpha_n$) sobre el conjunto de restricciones trasladado $\tilde{K}(t_n)$. Esto se realiza mediante un operador de proyección puntual $P_{g_n}$ que corrige la tensión para satisfacer el criterio de fluencia.

Análisis de Solvabilidad:
El artículo demuestra que, a pesar de la acoplamiento no lineal entre $\sigma_n$ y $\alpha_n$ en el paso de proyección, el sistema puede resolverse explícitamente punto a punto (casi en todo $\Omega$) mediante un análisis geométrico de las partes desviadoras de las tensiones en el espacio de matrices simétricas.

3. Resultados Principales

Los autores establecen los siguientes teoremas fundamentales:

• Estabilidad del Esquema Discreto (Teorema 4.1):
  Se prueban estimaciones de estabilidad uniformes en el paso de tiempo $\Delta t$ para las soluciones discretas. Estas estimaciones controlan normas en espacios de Sobolev y $L^2$ de la velocidad, la tensión, la variable interna y sus derivadas temporales.

  • Se demuestra que la diferencia entre la tensión de ensayo y la tensión proyectada tiende a cero cuando $\Delta t \to 0$.
  • Las cotas dependen de los datos iniciales y de las fuerzas externas, pero no del tamaño del paso de tiempo.

• Existencia de Solución Débil (Teorema 4.2):
  Utilizando las estimaciones de estabilidad, se aplican teoremas de compacidad (incluyendo el teorema de Aubin-Lions y argumentos de convergencia débil) para pasar al límite cuando $\Delta t \to 0$.

  • Se demuestra que el límite de las soluciones discretas satisface la formulación variacional del problema continuo (Problema 1.1).
  • Se verifica que la condición de pertenencia al conjunto $K(t; \alpha)$ se mantiene en el límite, incluso con un límite de fluencia $g$ que carece de regularidad espacial estricta.

• Unicidad de la Solución:
  Se prueba la unicidad mediante un argumento de energía. Considerando dos soluciones, se analiza la diferencia entre ellas y se demuestra que la energía de la diferencia es no positiva, lo que implica que las soluciones deben coincidir.

4. Contribuciones Clave

1. Generalización de las Condiciones sobre $g$: A diferencia de trabajos anteriores (como [2]), este artículo no asume continuidad espacial ni una cota inferior estrictamente positiva ($0 < c \le g$) para el límite de fluencia. Esto permite modelar materiales con límites de fluencia que pueden anularse o ser discontinuos.
2. Enfoque Constructivo: El método de proyección no solo prueba la existencia teórica, sino que proporciona un algoritmo directo y constructivo para la aproximación numérica. La estructura del esquema es compatible con la discretización espacial (por ejemplo, mediante Elementos Finitos), facilitando la implementación computacional.
3. Tratamiento de la Estructura Cuasi-Variacional: El artículo maneja rigurosamente la dependencia del conjunto de restricciones respecto a la variable interna $\alpha$, resolviendo la complejidad inherente a los problemas cuasi-variacionales en plasticidad dinámica.

5. Significado y Conclusión

Este trabajo es significativo porque cierra la brecha entre la teoría analítica de la elastoplasticidad dinámica con endurecimiento cinemático y la implementación numérica práctica.

• Rigor Matemático: Proporciona una base sólida para la existencia y unicidad bajo supuestos más débiles que los tradicionales.
• Aplicabilidad Numérica: El esquema propuesto es una base directa para desarrollar algoritmos de Elementos Finitos (FEM) para problemas complejos de plasticidad.
• Futuro: Los autores señalan que el siguiente paso natural es extender este análisis a la discretización completa (espacio-tiempo), investigar otros criterios de fluencia (como Tresca) y abordar el caso no viscoso ($\eta = 0$), que presenta desafíos analíticos mayores debido a la pérdida de control espacial de la velocidad.

En resumen, el artículo ofrece un marco robusto y constructivo para el análisis y la simulación de materiales que exhiben comportamiento elasto-plástico con endurecimiento cinemático lineal y viscosidad.