A blended approach for evolving phase fields using peridynamics: Cyclic loading in quasi-brittle fracture

Este artículo presenta una teoría de campo de fase de dos puntos dependiente de la historia que combina la peridinámica con la deformación plástica irreversible y el ablandamiento por daño para predecir cuantitativa y cualitativamente la fractura en materiales cuasi-frágiles bajo carga cíclica, cumpliendo con las leyes de la termodinámica y demostrando concordancia con experimentos de flexión en hormigón.

Lo esencial

Imagina que el concreto, el ladrillo o la piedra no son bloques sólidos e indestructibles, sino una enorme red de millones de pequeños resortes conectados entre sí. Cuando cargas un puente o un edificio, estás estirando o comprimiendo esos resortes.

Este artículo presenta una nueva forma de entender y predecir cómo se rompen estos materiales (como el concreto) cuando sufren daños, especialmente cuando se les aplica fuerza repetidamente (como el tráfico en un puente o los sismos).

Aquí tienes la explicación de su investigación usando analogías sencillas:

1. El Problema: ¿Cómo se rompe algo "frágil pero resistente"?

Los materiales como el concreto son "cuasi-frágiles". No se rompen de golpe como un vaso de vidrio (frágil puro), ni se doblan como el metal (plástico puro). Primero se agrietan un poco, se deforman internamente y luego se rompen.

Los métodos antiguos para simular esto eran como intentar arreglar un coche con dos mecánicos que no se hablan: uno calculaba cómo se mueve el coche (física) y el otro calculaba dónde aparecen las grietas (daño), y tenían que coordinarse manualmente. A veces fallaban, especialmente cuando el material se deformaba permanentemente (como cuando doblas un clip y no vuelve a su forma original).

2. La Solución: Una "Red Inteligente" que recuerda todo

Los autores proponen una nueva teoría llamada "Enfoque Mezclado" (Blended Approach). Imagina que en lugar de dos mecánicos, tienes un sistema nervioso único para todo el material.

• La Peridinámica (La Red): En lugar de usar puntos fijos en una cuadrícula (como en los videojuegos antiguos), imagina que cada punto del material tiene una "burbuja de visión" (un radio de acción). Dentro de esa burbuja, el punto "siente" y se conecta con sus vecinos. Si un vecino se aleja demasiado, el resorte que los une se rompe.
• El Campo de Fase (La Memoria): Aquí está la magia. Cada resorte no solo sabe si está estirado o comprimido, sino que tiene memoria.
  • Si lo estiras mucho, el resorte se daña y se vuelve más débil (su elasticidad baja).
  • Si lo sueltas, no vuelve a su estado original; se queda un poco estirado (deformación plástica).
  • El modelo usa dos "campos de fase" (como dos termómetros): uno para medir el daño por estiramiento (tracción) y otro para el daño por aplastamiento (compresión).

3. La Analogía del "Resorte con Memoria"

Imagina un resorte hecho de goma vieja:

1. Carga: Lo estiras. Al principio es duro, pero luego empieza a ceder.
2. Daño: Si lo estiras demasiado, el resorte se "quema" internamente. Ya no es tan fuerte.
3. Descarga: Si lo sueltas, no vuelve a su longitud original; se queda un poco largo.
4. Recarga: Si lo vuelves a estirar, es más fácil porque ya está dañado.

El modelo de los autores hace exactamente esto, pero a escala microscópica para millones de "resortes" a la vez. Lo genial es que todo esto sucede automáticamente siguiendo las leyes de la física (la segunda ley de Newton, la misma que dice que la fuerza es igual a masa por aceleración). No necesitan inventar reglas extra para que aparezca una grieta; la grieta aparece naturalmente cuando los resortes se rompen.

4. ¿Por qué es importante? (El "Efecto Tamaño")

En la vida real, un bloque de concreto pequeño se rompe con menos fuerza relativa que un bloque gigante. Esto se llama "efecto de tamaño". Los métodos antiguos a menudo fallaban al predecir esto o necesitaban trucos matemáticos complicados.

Este nuevo modelo logra predecir automáticamente por qué un edificio grande es más "débil" proporcionalmente que un ladrillo pequeño, simplemente porque la red de resortes y sus memorias interactúan de forma natural.

5. La Prueba: Simulando el mundo real

Los autores probaron su teoría simulando:

• Pruebas de flexión: Como doblar una viga de concreto hasta que se rompa.
• Cargas cíclicas: Como subir y bajar una carga repetidamente (simulando un puente con tráfico).
• Grietas en esquinas: Cómo se propagan las grietas en formas extrañas (como una pieza en forma de L).

El resultado: Sus simulaciones computacionales coincidieron casi perfectamente con experimentos reales de laboratorio. Pudieron ver cómo se formaban las grietas, cómo el material se deformaba permanentemente y cómo absorbía energía, todo sin usar una malla de puntos rígidos (es un método "sin malla", más flexible).

En resumen

Esta investigación es como darles consciencia y memoria a los átomos de un material de construcción. En lugar de tratar el concreto como un bloque estático, lo tratan como una red viva de conexiones que recuerdan cuánto han sufrido, cuándo se han dañado y cómo deben comportarse al siguiente golpe.

Esto permite a los ingenieros diseñar edificios y puentes más seguros, sabiendo exactamente cómo se comportarán bajo estrés repetido, terremotos o cargas pesadas, sin necesidad de construir y destruir cientos de prototipos físicos.

Resumen técnico

Aquí presento un resumen técnico detallado del artículo "A blended approach for evolving phase fields using peridynamics: Cyclic loading in quasi-brittle fracture" (Un enfoque híbrido para la evolución de campos de fase utilizando peridinámica: Carga cíclica en fractura cuasi-frágil), escrito por Hayden Bromley y Robert Lipton.

1. Planteamiento del Problema

La identificación de una teoría de campo adecuada para formular el problema de valor inicial y de frontera (IBVP) en la mecánica de fractura ha sido un desafío histórico. Los métodos tradicionales de campo de fase a menudo requieren ecuaciones de evolución separadas para el desplazamiento (basadas en la segunda ley de Newton) y para el campo de fase (basadas en ecuaciones tipo Ginzburg-Landau o restricciones variacionales). Esto puede complicar la simulación de fenómenos dinámicos, la histéresis en cargas cíclicas y la localización de daños en materiales cuasi-frágiles (como el concreto), donde la disipación de energía plástica y el ablandamiento elástico juegan un papel crucial.

El problema central abordado es cómo modelar la degradación del material, la plasticidad irreversible y la fractura dinámica bajo cargas cíclicas utilizando un marco unificado que respete estrictamente las leyes de la termodinámica y la mecánica clásica, sin depender de ecuaciones de evolución ad hoc para el daño.

2. Metodología

Los autores proponen un enfoque híbrido (blended) que combina la peridinámica (una teoría no local) con un campo de fase de dos puntos dependiente de la historia.

• Fundamento Teórico:

  • Se utiliza exclusivamente la segunda ley de Newton para la evolución del desplazamiento. No existe una ecuación separada para la evolución del campo de fase; este emerge como parte de una ley constitutiva dependiente de la historia.
  • El campo de fase se define en pares de puntos $(y, x)$ dentro del horizonte de interacción no local.
  • La deformación se descompone aditivamente en deformación elástica y deformación irreversible (plástica).

• Leyes Constitutivas con Memoria:

  • Se introducen dos campos de fase independientes: $\gamma^+$ para tracción y $\gamma^-$ para compresión. Esto permite modelar comportamientos asimétricos (el material puede fallar en tracción pero resistir elásticamente en compresión).
  • Se define un parámetro de descarga ($\beta$) que relaciona la deformación plástica con el ablandamiento elástico. Este parámetro controla la transición entre comportamiento frágil (sin plasticidad, ablandamiento puro) y cuasi-frágil (con plasticidad y ablandamiento).
  • La ley constitutiva describe un "envolvente de resistencia" (strength envelope) y una familia de leyes de descarga que folian el dominio de resistencia, garantizando que los pares tensión-deformación sean admisibles.

• Formulación Matemática:

  • El problema se formula como un IBVP bien planteado en un espacio de Banach.
  • Se demuestra la existencia y unicidad de la solución utilizando el teorema del punto fijo de Banach, basándose en la continuidad Lipschitz del campo de tensiones respecto al desplazamiento.
  • Se establece un balance de energía riguroso. La tasa de disipación de energía es no negativa, cumpliendo con la segunda ley de la termodinámica. La energía de fractura se recupera directamente del modelo mediante teoría de medida geométrica, sin depender de la longitud de escala de no-localidad ($\epsilon$).

• Implementación Numérica:

  • Se utiliza un método sin malla (mesh-free) basado en peridinámica, donde la discretización se realiza mediante cuadratura sobre nodos.
  • Se emplea un método de relajación dinámica para simulaciones cuasi-estáticas y un esquema de integración temporal de diferencias centrales para dinámicas.

3. Contribuciones Clave

1. Unificación de Dinámica y Daño: Eliminación de la necesidad de ecuaciones de evolución separadas para el campo de fase. El daño y la fractura son consecuencias directas de la evolución del desplazamiento bajo una ley constitutiva no local con memoria.
2. Modelado de Histéresis Cíclica: Capacidad de capturar la histéresis en cargas de tracción y compresión, incluyendo la acumulación de deformación plástica irreversible y el ablandamiento elástico, mediante el uso de dos campos de fase y el parámetro $\beta$.
3. Independencia de la Escala de No-Localidad: Se demuestra que la energía de fractura de grietas planas es independiente de la longitud de horizonte $\epsilon$, recuperando directamente la energía de Griffith.
4. Efecto de Tamaño Determinista: El modelo predice naturalmente el efecto de tamaño en materiales cuasi-frágiles (la resistencia aparente disminuye al aumentar el tamaño de la muestra) sin necesidad de reglas empíricas adicionales, emergiendo directamente del problema de valor inicial.
5. Calibración Simplificada: El modelo se calibra utilizando solo propiedades materiales esenciales: módulo elástico, resistencia a la tracción, tasa de liberación de energía crítica y la relación de descarga ($\beta$).

4. Resultados

Los autores validaron el modelo mediante simulaciones numéricas comparadas con datos experimentales en varios casos de prueba:

• Fractura Modo I (Tracción) en Concreto: Simulación de pruebas de flexión en tres puntos con carga monótona y cíclica. El modelo reprodujo con precisión la curva Fuerza vs. Apertura de la grieta (CMOD), incluyendo la histéresis observada en experimentos reales.
• Fractura Mixta (Cortante + Tracción): Se simuló una viga de concreto con una entalla desplazada. El modelo predijo correctamente las trayectorias de grieta y las curvas de fuerza para diferentes ángulos de carga, mostrando concordancia cualitativa y cuantitativa con los datos de Jenq y Shah (1988).
• Specimen en Forma de L: Se modeló la fractura mixta en un panel de concreto sin entalla previa. El modelo capturó la transición de modo I a modo mixto (deslizamiento) a medida que la grieta se propagaba desde la esquina, coincidiendo con los resultados experimentales de Winkler (2001).
• Efecto de Tamaño: Se probaron vigas de concreto de tres tamaños diferentes. El modelo capturó la disminución de la tensión nominal en la falla a medida que aumentaba el tamaño de la viga, validando la capacidad del modelo para predecir el efecto de tamaño sin parámetros ajustables adicionales.
• Convergencia: Se realizó un estudio de convergencia de malla ($m$-convergence y $\epsilon$-convergence), demostrando que la solución numérica converge a la teoría de campo a medida que se refina la discretización.

5. Significancia

Este trabajo representa un avance significativo en la mecánica de la fractura computacional por varias razones:

• Rigor Termodinámico: Proporciona un marco matemáticamente bien planteado donde la disipación de energía y el balance de energía son consecuencias directas de las ecuaciones de movimiento, no postulados.
• Versatilidad: El enfoque híbrido es capaz de manejar tanto la fractura frágil como la cuasi-frágil, y transiciones entre carga estática y dinámica, utilizando un único conjunto de ecuaciones.
• Eficiencia Computacional y Predictiva: Al eliminar ecuaciones de evolución adicionales para el daño, el método simplifica la implementación numérica mientras mantiene la capacidad de predecir patrones de grietas complejos y fenómenos de histéresis.
• Aplicabilidad Práctica: La capacidad de calibrar el modelo con propiedades materiales estándar y predecir el efecto de tamaño lo hace altamente relevante para la ingeniería civil y la evaluación de estructuras de concreto y materiales geológicos.

En conclusión, los autores demuestran que es posible derivar la evolución de daños y fracturas, incluyendo plasticidad y ciclos de carga, directamente de la segunda ley de Newton en un marco no local, ofreciendo una alternativa robusta y teóricamente sólida a los métodos de campo de fase tradicionales.