💻 computer science

Phase-Field Peridynamics

Este trabajo introduce una nueva formulación de peridinámica de campo de fase que evita las inestabilidades numéricas asociadas a la eliminación progresiva de enlaces mediante la degradación continua de su contribución energética y el uso de una función de degradación cinemática separada, garantizando así la precisión y consistencia termodinámica en la simulación de fracturas.

Autores originales: Kai Partmann, Christian Wieners, Michael Ortiz, Kerstin Weinberg

Publicado 2026-04-13

📖 4 min de lectura☕ Lectura para el café

CC BY 4.0

Autores originales: Kai Partmann, Christian Wieners, Michael Ortiz, Kerstin Weinberg

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 tienes una pieza de goma muy elástica y quieres predecir exactamente cómo se romperá si la estiras o la golpeas. En el mundo de la ingeniería y la física, esto es un gran desafío: los modelos tradicionales a menudo se "atascan" o fallan justo cuando la goma empieza a romperse, porque no saben cómo manejar el hecho de que el material deja de ser continuo.

Este artículo presenta una nueva y brillante forma de simular estas roturas, llamada Peridinámica de Campo de Fase. Aquí te lo explico con analogías sencillas:

1. El problema: La "tijera" digital

Imagina que la materia está hecha de millones de puntos conectados por pequeños "hilos" o cuerdas (llamados bonds en el texto).

El método antiguo: Cuando un hilo se estira demasiado, el modelo lo corta de golpe (como usar unas tijeras). El problema es que si cortas muchos hilos de golpe en una zona, la red se vuelve inestable, como una telaraña que pierde su forma y colapsa. Esto hace que las simulaciones se vuelvan locas y den resultados erróneos, especialmente en materiales complejos.
La analogía: Es como intentar arreglar un puente viejo quitando vigas enteras de golpe. El puente se tambalea y puede caerse antes de tiempo.

2. La solución: El "atenuador de volumen"

Los autores proponen no cortar los hilos, sino debilitarlos suavemente.

En lugar de usar tijeras, usan un control de volumen (un dimmer de luz). Cuando un hilo se estira demasiado, en lugar de desaparecer, su "fuerza" se va bajando poco a poco hasta que es casi cero.
La analogía: Imagina que tienes una red de luces de Navidad. Si un cable se daña, no lo cortas y dejas un hueco negro. En su lugar, haces que esa bombilla se vaya apagando suavemente hasta que se vuelve invisible. La red sigue manteniendo su forma y estructura, pero esa parte ya no aporta fuerza.

3. El truco de ingeniería: Dos tipos de debilitamiento

El gran secreto de este nuevo método es que usa dos tipos de control para los hilos dañados:

Control de Energía (El "Atenuador"): Reduce la fuerza que el hilo puede ejercer. Si el hilo está muy dañado, ya no tira con fuerza. Esto simula que el material se está rompiendo.
Control de Movimiento (El "Estabilizador"): Aquí está la magia. Aunque el hilo esté casi "apagado" (sin fuerza), el modelo sigue usándolo para calcular cómo se mueve el resto de la red.
- Por qué es importante: Si quitas los hilos rotos de los cálculos de movimiento, la red pierde su forma y se vuelve inestable (como un edificio sin vigas). Al mantenerlos "en la ecuación" pero con fuerza cero, la red mantiene su geometría perfecta y la simulación no se rompe.

4. La "Regla de Oro" (La constante de normalización)

Para que todo esto funcione y sea realista, los autores tuvieron que calcular una fórmula matemática muy precisa (llamada constante de normalización).

La analogía: Es como si fueras a cocinar un pastel y tuvieras que saber exactamente cuánta harina usar para que, al hornearse, el pastel tenga el tamaño exacto que promete la receta, sin importar si usas un molde grande o pequeño.
Ellos demostraron matemáticamente cómo ajustar esta "cantidad de harina" para cualquier tipo de material y cualquier tamaño de simulación, asegurando que la energía gastada en romper el material sea la misma que predice la teoría clásica de la física (la teoría de Griffith).

5. ¿Qué lograron probar?

Pusieron a prueba su nuevo modelo en tres situaciones difíciles:

Romper tirando (Modo I): Como estirar un chicle hasta que se parte. Funcionó perfecto.
Romper cortando (Modo II): Como intentar cortar una goma con un movimiento de tijera. Aquí es donde los modelos antiguos fallaban (se volvían inestables), pero el nuevo modelo lo hizo con estabilidad.
El experimento Kalthoff-Winkler: Un choque a alta velocidad contra una placa de acero. El modelo predijo que las grietas se abrirían en un ángulo de 70 grados, exactamente igual que en los experimentos reales de laboratorio.

En resumen

Este trabajo es como pasar de construir con bloques de Lego que se caen si quitas uno a construir con una masa de arcilla inteligente. Si la arcilla se rompe, no desaparece de la mesa; se vuelve suave y flexible, manteniendo la forma del objeto mientras se rompe.

Esto permite a los ingenieros simular cómo fallan los puentes, los aviones o los materiales de construcción en situaciones extremas con una precisión y seguridad mucho mayores, sin que la computadora se "confunda" cuando las cosas empiezan a romperse.

Resumen Técnico: Peridinámica de Campo de Fase (PFPD)

1. Planteamiento del Problema

La modelización computacional de la fractura en estructuras mecánicas es un desafío central. Los métodos tradicionales basados en la mecánica del continuo (como FEM) requieren derivadas espaciales que se vuelven indefinidas en las discontinuidades (grietas), obligando a técnicas costosas de remallado o enriquecimiento (XFEM).

La Peridinámica (PD), introducida por Silling, reformula el equilibrio de momento lineal como una ecuación integro-diferencial, eliminando la necesidad de derivadas espaciales y siendo naturalmente adecuada para modelar discontinuidades. Sin embargo, las formulaciones existentes presentan limitaciones críticas:

Formulaciones de Correspondencia: Permiten usar leyes constitutivas clásicas y evitan el "efecto de superficie", pero sufren de inestabilidades numéricas (modos de energía cero).
Enfoque de Borrado de Enlaces (Bond Deletion): La forma tradicional de modelar fractura en PD consiste en eliminar irreversiblemente los enlaces cuando superan un umbral crítico. En las formulaciones de correspondencia, la eliminación progresiva de enlaces degrada el tensor de forma necesario para aproximar el gradiente de deformación. Esto puede llevar a que el tensor se vuelva mal condicionado o singular, causando inestabilidades numéricas y fallos en la simulación, especialmente en modos de carga complejos (como el Modo II).

2. Metodología Propuesta

Los autores introducen un nuevo enfoque llamado Peridinámica de Campo de Fase (PFPD) que integra la teoría de campo de fase dentro del marco de correspondencia asociada a enlaces (Bond-Associated Quadrature Point - BAQP).

Innovaciones Clave:

Degradación Continua en lugar de Borrado: En lugar de eliminar instantáneamente los enlaces al alcanzar un umbral crítico, la contribución energética de cada enlace se degrada continuamente mediante un parámetro de campo de fase de enlace ( $s$ ).
Funciones de Degradación Separadas:
1. Degradación Energética ( $g(s)$ ): Reduce la respuesta mecánica (tensión) del enlace suavemente a medida que el daño avanza ( $s \to 1$ ).
2. Degradación Cinemática ( $h(s)$ ): Una función separada que degrada la contribución cinemática de los enlaces dañados al tensor de forma. Esto es crucial: mantiene la precisión de la aproximación del gradiente de deformación y la estabilidad del tensor de forma incluso cuando el material está severamente dañado, evitando la singularidad numérica.
Derivación Analítica de la Constante de Normalización ( $c_0$ ):
- Para garantizar la consistencia termodinámica con la teoría de fractura de Griffith, se deriva analíticamente una constante de normalización para funciones de núcleo esféricas generales.
- La constante $c_0$ se expresa como el cociente de dos integrales unidimensionales sobre el intervalo unitario. Esto permite calcular $c_0$ exactamente para cualquier función de núcleo (cúbica, lineal, constante, etc.) sin depender de simulaciones numéricas empíricas.
Criterio de Evolución del Daño: Se utiliza un principio variacional basado en la energía para derivar la evolución del parámetro $s$ . Se introduce una variable de historia para asegurar la irreversibilidad del daño. El criterio se basa en la energía de deformación elástica (o tensión principal máxima) para activar el crecimiento de la grieta.

3. Contribuciones Principales

Estabilidad Numérica en Formulaciones de Correspondencia: El método resuelve el problema de la inestabilidad en la formulación de correspondencia asociada a enlaces (BAQP) durante la propagación de grietas, un problema que limita la aplicación de esta formulación más precisa en problemas de fractura.
Consistencia Termodinámica Rigurosa: Proporciona una derivación analítica exacta para la constante de normalización $c_0$ para núcleos esféricos generales, asegurando que la disipación de energía total coincida con la tasa de liberación de energía de Griffith.
Marco Unificado: Integra exitosamente la regularización de campo de fase (que evita el seguimiento explícito de grietas) con la peridinámica de alta precisión (BAQP), combinando las ventajas de ambos mundos.
Algoritmo Eficiente: Se presenta una implementación algorítmica que mantiene todos los enlaces en la aproximación cinemática (evitando la reconstrucción de matrices por eliminación binaria), degradando solo sus contribuciones de manera continua.

4. Resultados Numéricos

Los autores validan el modelo mediante varios ejemplos numéricos en 3D:

Propagación de Grietas en Modo I y Modo II:
- El modelo predice correctamente la trayectoria de la grieta en placas agrietadas.
- Hallazgo Crítico: En el caso de Modo II (cizalladura), el modelo estándar BAQP con borrado de enlaces falla debido a la degradación del tensor de forma. El modelo PFPD permanece estable y produce resultados precisos, demostrando su superioridad en cargas complejas.
Prueba de Tensión de Borde (Boundary Tension Test - BTT):
- Se utilizó para verificar la constante de normalización $c_0$ .
- Se probaron diferentes relaciones de horizonte ( $m_\delta$ ) y funciones de núcleo (cúbico B-spline y lineal).
- Resultado: Los patrones de ramificación de grietas, los ángulos de ramificación y los puntos de inicio fueron consistentes e independientes de la elección del núcleo o la relación de horizonte, confirmando la corrección de la derivación analítica de $c_0$ .
- La velocidad de la punta de la grieta se mantuvo por debajo de la mitad de la velocidad de la onda de Rayleigh, como se espera teóricamente.
Experimento de Kalthoff-Winkler:
- Simulación de impacto en una placa de acero con entallas predefinidas.
- El modelo PFPD reprodujo con precisión el ángulo de propagación de la grieta (~70°), característico de la fractura dominada por Modo II observada experimentalmente.
- Comparado con modelos OSB (Ordinary State-Based) y BAQP estándar, el PFPD mostró perfiles de velocidad de punta de grieta más suaves y una mejor concordancia con los datos experimentales.

5. Significado e Impacto

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

Robustez: Permite utilizar la formulación de correspondencia más precisa (BAQP) en problemas de fractura dinámica sin sufrir las inestabilidades numéricas que históricamente la han limitado.
Versatilidad: Al no depender de la eliminación de enlaces, el modelo es más robusto frente a la malla y a la geometría de la grieta, facilitando la simulación de fenómenos complejos como la ramificación y la coalescencia de grietas.
Fundamento Teórico: La derivación analítica de la constante de normalización cierra una brecha teórica importante, permitiendo que los modelos de peridinámica de campo de fase sean predictivos y consistentes con la mecánica de fractura clásica de Griffith para una amplia gama de funciones de influencia.

En conclusión, el enfoque PFPD ofrece una herramienta robusta, precisa y termodinámicamente consistente para la simulación de fractura dinámica, superando las limitaciones de los métodos de borrado de enlaces tradicionales.