A 2D discrete model with a phase-field regularisation applied to dynamic crack propagation

Cuando un material frágil como el vidrio o el acero endurecido se rompe, no simplemente se parte en dos. En su lugar, una grieta recorre el sólido, ramificándose a menudo en una compleja red de fracturas que hace añicos el objeto en muchas piezas. Predecir exactamente por dónde viajarán estas grietas y con qué rapidez se moverán es uno de los desafíos más difíciles en la ingeniería. Si un modelo utilizado para simular esta destrucción es defectuoso, podría sugerir que un puente fallará en una dirección segura cuando en realidad colapsará catastróficamente. Durante décadas, los científicos han intentado construir modelos informáticos que imiten la estructura atómica de los materiales, tratándolos como una cuadrícula de diminutos puntos conectados por resortes. Aunque estos modelos son potentes, sufren un fallo peculiar: el camino que sigue la grieta a menudo depende de cómo se rote la cuadrícula de la computadora, en lugar de depender de la física del material. Es como si una grieta en un trozo de vidrio real cambiara repentinamente de dirección simplemente porque el ingeniero giró la pantalla de la computadora.

Para resolver este problema, Matías Braun, un investigador de la Universidad Nebrija en España, ha desarrollado una nueva forma de simular la fractura dinámica. Su enfoque combina una cuadrícula tradicional de puntos con un método de "campo de fase" (phase-field), una técnica que trata la grieta no como una línea afilada que salta de punto en punto, sino como una zona de daño suave y difusa que se propaga a través del material. Al combinar estas dos ideas, Braun creó un modelo que permite que las grietas viajen en cualquier dirección, independientemente de cómo esté orientada la cuadrícula informática subyacente. Probó este nuevo método en tres escenarios distintos: una placa de acero golpeada por un impacto de alta velocidad, una placa con una muesca estirada hasta que se rompe y un cilindro metálico grueso que explota desde el interior hacia afuera. En cada caso, el nuevo modelo produjo resultados que coincidieron con experimentos del mundo real y eliminaron el sesgo artificial causado por la cuadrícula de la computadora, ofreciendo una visión más clara y fiable de cómo fallan realmente los materiales.

El núcleo del problema que abordó Braun reside en cómo las computadoras manejan el mundo invisible de los átomos. En su trabajo anterior, utilizó un modelo discreto donde el material estaba representado por una red triangular, un patrón de puntos conectados por líneas. Cuando el material era sometido a tensión, las conexiones entre los puntos se rompían, simulando una griza. Sin embargo, debido a que los puntos estaban dispuestos en un patrón fijo, la grieta se veía obligada a seguir las líneas de la red. Si la red se rotaba, la grieta cambiaba su trayectoria, aunque el material en sí no hubiera cambiado. Este es un problema significativo porque, en la realidad, una grieta debería seguir el camino de menor resistencia determinado por las propiedades del material, no el ángulo arbitrario de una malla informática.

La solución de Braun fue introducir una regularización de campo de fase. En lugar de forzar la ruptura de conexiones específicas entre los puntos, asignó un valor a cada punto que representa qué tan dañado está. Este valor oscila entre cero, lo que significa que el material está perfectamente intacto, y uno, lo que significa que está completamente roto. La grieta se convierte en una región donde este valor aumenta gradualmente, creando una transición suave en lugar de una línea dentada. Esto permite que la grieta se curve y se ramifique en cualquier dirección, ignorando eficazmente la estructura rígida de la red. Crucialmente, Braun calibró este sistema utilizando dos propiedades fundamentales del material: su resistencia a la tracción, que es la fuerza necesaria para iniciar una grieta, y su energía de fractura, que es la energía necesaria para mantener el crecimiento de la grieta. Al anclar la simulación a estas constantes físicas reales, el modelo asegura que los resultados dependan del material, no de las matemáticas.

Para demostrar que este nuevo método funcionaba, Braun realizó tres pruebas rigurosas. La primera fue la prueba de Kalthoff–Winkler, un experimento clásico donde una placa de acero con una muesca precortada es golpeada por un proyectil a alta velocidad. En este escenario, una grieta sale disparada desde la muesca en un ángulo específico. Los modelos anteriores a menudo predecían el ángulo incorrecto dependiendo de cómo se rotara la red. En las simulaciones de Braun, la grieta viajó consistentemente a un ángulo de aproximadamente 70 grados, coincidiendo perfectamente con las observaciones experimentales. Este resultado se mantuvo constante independientemente de si la red estaba alineada con la placa o rotada 15 o 30 grados. La velocidad de la grieta también se mantuvo estable, viajando a aproximadamente el 70 por ciento de la velocidad del sonido en el material, de forma independiente a la orientación de la red.

La segunda prueba involucró una placa con una muesca de un solo borde bajo tensión, una configuración utilizada para estudiar cómo crecen las grietas lentamente bajo un estiramiento constante. Aquí, el objetivo era ver si la grieta viajaría en línea recta desde la muesca hasta el borde opuesto. Los modelos de red tradicionales a menudo forzaban a la grieta a zigzaguear o desviarse según la alineación de la red. El modelo de Braun, sin embargo, produjo una grieta recta para cada orientación de la red. La fuerza requerida para romper la placa también fue consistente, variando en menos de un dos por ciento independientemente de cómo se rotara la red. Esto demostró que el modelo podía predecir con precisión tanto el camino como la resistencia del material sin la influencia artificial de la malla informática.

La tercera prueba y la más compleja fue la fragmentación de un cilindro grueso sometido a una explosión interna. En este escenario, ninguna grieta única es guiada por una muesca; en su lugar, muchas grietas deben nuclearse espontáneamente y lanzarse hacia afuera, destrozando el cilindro en docenas de piezas. Este es un evento caótico donde la simetría de la red podría dictar fácilmente el patrón de los fragmentos. Cuando Braun ejecutó la simulación con la red centrada en el eje del cilindro, las grietas formaron un patrón simétrico. Sin embargo, cuando desplazó la red ligeramente para que ya no estuviera perfectamente centrada, las grietas se distribuyeron de forma irregular y el número de fragmentos cambió. Esto confirmó que el modelo era sensible a la configuración física más que a la simetría de la red. La simulación predijo entre 12 y 19 fragmentos, un número que se encuentra dentro del rango de otros métodos avanzados, pero que es inferior a algunos, probablemente porque la zona de daño suave en el modelo evita la formación de grietas que están demasiado cerca unas de otras.

A lo largo de estas pruebas, los investigadores rastrearon cuidadosamente la energía del sistema. Compararon la energía perdida a medida que el material se rompía contra la energía requerida para crear las nuevas superficies de la grieta. En algunos casos, hubo un pequeño desajuste entre la energía calculada por el modelo de daño y la energía realmente liberada por las conexiones rotas. Esta es una limitación conocida del enfoque híbrido, donde la fuerza impulsora de la grieta se calcula de forma diferente a la energía de los enlaces rotos. Sin embargo, el equilibrio de la energía mecánica se mantuvo preciso dentro de un medio por ciento, confirmando que la simulación era físicamente sólida. Los resultados sugieren que, si bien el modelo no es perfecto, representa un paso significativo adelante en la comprensión de la fractura dinámica.

Las implicaciones de este trabajo se extienden más allá de simplemente mejores gráficos por computadora. Al eliminar la dependencia de la orientación de la red, los ingenieros pueden confiar en que sus simulaciones de pruebas de choque, fallos estructurales o pruebas de materiales están reflejando la verdadera física del material. El modelo logra cerrar la brecha entre la naturaleza discreta de los átomos y el comportamiento continuo de los objetos sólidos. Demuestra que, con las herramientas matemáticas adecuadas, es posible simular el proceso caótico y violento de la ruptura sin dejar que la cuadrícula de la computadora dicte el resultado. El estudio concluye que este enfoque ofrece una forma robusta de predecir cómo fallan los materiales bajo condiciones extremas, proporcionando un camino más claro para diseñar estructuras más seguras y resilientes.