🔬 materials science

Strength-degradation phase-field regularization of cohesive fracture: the antiplane case

Este artículo introduce un modelo de campo de fase de degradación de resistencia para la fractura antiplana que unifica el análisis de límite, la plasticidad y diversos regímenes de fractura al desacoplar la resistencia, la rigidez y la tenacidad, permitiendo así que la nucleación y propagación de grietas estén gobernadas por una superficie de resistencia convexa arbitraria y una escala de longitud elasto-cohesiva en lugar de estar restringidas por la energía elástica.

Autores originales: Blaise Bourdin, Corrado Maurini

Publicado 2026-08-03

📖 8 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Blaise Bourdin, Corrado Maurini

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 estás tratando de predecir cómo se romperá un trozo de vidrio, un bloque de arcilla o incluso una viga de metal. Durante mucho tiempo, los científicos han tenido dos formas muy diferentes de abordar este problema. Por un lado, está la visión "frágil": los materiales se quiebran repentinamente, como una rama seca, y lo único que importa es cuánta energía se necesita para crear una nueva superficie de grieta. Por otro lado, está la visión "dúctil": los materiales se aplastan y se estiran antes de romperse, como un caramelo masticable, y la clave es cuánto cede o se deforma el material bajo presión.

La parte difícil es que los materiales del mundo real suelen hacer ambas cosas. Un trozo de metal puede estirarse un poco (plasticidad) antes de romperse finalmente (fractura). Durante décadas, los modelos computacionales tuvieron dificultades para manejar esta mezcla. Los populares modelos de "campo de fase" (phase-field), que son como simulaciones digitales que permiten que las grietas crezcan de forma suave en lugar de saltar de un píxel a otro, tenían un fallo importante: vinculaban la resistencia del material directamente a qué tan "difusa" se veía la grieta en la pantalla de la computadora. Si querías simular un material fuerte, tenías que hacer que la "difusidad" de la computadora fuera diminuta, lo que requería supercomputadoras. Si hacías la difusidad más grande para ahorrar tiempo, el material se volvía mágicamente más débil. Era un juego frustrante de "elige dos": podías tener resistencia, podías tener velocidad o podías tener una forma de grieta realista, pero no podías tener las tres al mismo tiempo.

Este artículo presenta una nueva forma de modelar la ruptura de materiales que finalmente desenreda este nudo. Los autores, Blaise Bourdin y Corrado Maurini, proponen un modelo donde la resistencia del material es un ingrediente separado e independiente, tal como lo es su rigidez o su tenacidad. Probaron esta idea en un entorno simplificado llamado "cizalladura antiplana" (piensa en deslizar un mazo de cartas lateralmente) y descubrieron que su nuevo modelo funciona de maravilla. Puede simular desde un metal que se estira y se deforma (plasticidad) hasta un quiebre limpio y repentino (fractura frágil) sin necesidad de cambiar las reglas del juego. Lo más importante es que demostraron que la "difusidad" de la simulación por computadora es solo una herramienta numérica para ayudar a que las matemáticas funcionen, no una propiedad física que cambie la fuerza del material. Esto significa que los científicos finalmente pueden simular comportamientos de ruptura complejos en computadoras estándar sin que los resultados cambien solo por ajustar el tamaño de la rejilla.

El nuevo enfoque de "Resistencia Primero"

El núcleo de este artículo es una nueva receta matemática para simular cómo se rompen las cosas. Los autores lo llaman "regularización de campo de fase por degradación de la resistencia". Es un nombre complicado, así que vamos a desglosarlo con una analogía sencilla.

Imagina que un material es como una multitud de personas tomadas de la mano. En los modelos antiguos, cuando la multitud empezaba a sufrir estrés, los brazos de las personas (la rigidez) se debilitaban cada vez más hasta que soltaban la mano. El problema era que el punto en el que soltaban dependía enteramente de qué tan de cerca las estuvieras observando (la "longitud de regularización"). Si las observabas de cerca, aguantaban más; si las observabas desde lejos, soltaban antes.

En este nuevo modelo, las personas no debilitan sus brazos. En su lugar, tienen una fuerza de agarre máxima. Mientras la tracción esté por debajo de esa fuerza, se mantienen firmes. Una vez que la tracción alcanza ese límite, el agarre comienza a degradarse y eventualmente sueltan. Crucialmente, esta "fuerza de agarre máxima" es una propiedad fija de las personas, no un efecto secundario de cómo las estás observando. Esto permite al modelo manejar materiales que son fuertes pero no muy elásticos, o materiales que se estiran mucho antes de romperse, todo dentro del mismo marco de trabajo.

Las tres etapas de la ruptura

Los autores realizaron simulaciones para ver cómo se comporta su modelo bajo diferentes condiciones. Descubrieron que su único modelo puede transicionar naturalmente a través de tres etapas distintas de falla, dependiendo del tamaño del objeto y las propiedades del material:

La fase "elástica/maleable" (Yielding a pequeña escala): Cuando una grieta comienza a formarse en un área pequeña, el material alrededor de la punta no solo se quiebra. Se estira y se deforma, creando una "zona de proceso" donde el material está cediendo. El modelo captura esto perfectamente, mostrando una zona de deformación que coincide con las teorías clásicas de la plasticidad.
La fase "cohesiva": A medida que la grieta se abre, el material no se rompe instantáneamente. En su lugar, mantiene una fuerza "cohesiva", como una cinta adhesiva que se va despegando lentamente. La fuerza necesaria para separarlo disminuye a medida que el espacio se ensancha. Este es el comportamiento de la "grieta cohesiva", que es esencial para entender cómo fallan los materiales reales.
La fase "frágil": Finalmente, si el espacio es lo suficientemente ancho, el material se rinde por completo y se quiebra. La grieta se convierte en una rotura limpia y afilada, tal como en los modelos más simples de fractura frágil.

La belleza de este trabajo es que no necesita cambiar entre diferentes fórmulas matemáticas para describir estas tres etapas. Todo es un proceso continuo. El modelo simplemente evoluciona de un estado al siguiente basándose en la física, no en la elección de un programador humano.

El experimento de "Surfear" y la estela plástica

Una de las partes más fascinantes del artículo es una simulación que los autores llaman "surfear". Imagina una grieta que ya existe, y tú la arrastras a través del material a una velocidad constante, como un surfista sobre una ola. Los autores realizaron esta simulación de dos maneras: una donde la deformación del material era reversible (como una banda elástica) y otra donde era irreversible (como la arcilla que permanece deformada).

En el caso reversible, la energía necesaria para romper el material coincidió exactamente con lo que la teoría predecía. Pero en el caso irreversible, algo interesante sucedió. A medida que la grieta se movía, dejaba atrás una "estela plástica": un rastro de material permanentemente deformado, como la estela de un bote. Debido a que el material tenía que ser deformado para permitir el paso de la grieta, se requería más energía para romper el material de lo que la teoría para una simple grieta frágil sugeriría. Los autores midieron esta "tenacidad efectiva" y descubrieron que era aproximadamente un 30% mayor que la tenacidad base del material. Esto explica por qué algunos materiales parecen más tenaces en la vida real de lo que predicen sus fórmulas básicas: la deformación plástica delante de la punta de la grieta está realizando un trabajo adicional.

Cerrando la brecha con una muesca en V

Para demostrar que su modelo funciona para formas complejas, los autores simularon una "muesca en V" (un corte afilado en una pieza de material, como la esquina de una forma de Pac-Man). Aplicaron fuerza a esta muesca y observaron qué sucedía.

Descubrieron que, a medida que aumentaba la carga, el material pasaba por las mismas tres etapas: primero, se formó una pequeña zona plástica en la punta (la fase "maleable"); luego, una grieta cohesiva comenzó a crecer (la fase de "cinta adhesiva"); y finalmente, una grieta frágil de repente nucleó y atravesó el material.

El resultado más impresionante aquí es que la fuerza requerida para iniciar la grieta siguió una "ley universal". Independientemente de si la muesca era muy afilada (como una grieta) o muy roma (como un borde recto), las predicciones del modelo colapsaron en una sola curva. Esta curva conecta suavemente la resistencia del material (qué tan fuerte hay que tirar para romperlo) con su tenacidad (cuánta energía se necesita para romperlo). Esto sugiere que el modelo puede predecir la falla para cualquier forma sin necesidad de adivinar qué fórmula usar.

Por qué esto es importante

Los autores advierten cuidadosamente que estos resultados provienen de simulaciones y pruebas matemáticas, no de experimentos físicos con materiales reales. Sin embargo, las simulaciones son robustas y coinciden perfectamente con los límites teóricos conocidos. Han demostrado que, al cambiar la forma en que su modelo degrada la resistencia en lugar de la rigidez, han creado una herramienta que es tanto flexible como precisa.

Anteriormente, si querías simular un material que fuera fuerte pero frágil, estabas atrapado. O tenías que usar un modelo que era computacionalmente imposible (requiriendo una malla tan fina que colapsaría tu computadora) o un modelo que daba la respuesta incorrecta. Este nuevo enfoque permite que la "longitud de regularización" (la difusidad de la computadora) sea solo una herramienta numérica. Puedes hacerla pequeña para obtener una grieta afilada, o más grande para ahorrar potencia de cómputo, y la resistencia y la tenacidad del material permanecerán exactamente iguales.

En resumen, este artículo ofrece un lenguaje unificado para la fractura. Reúne los mundos de la plasticidad (deformación), la fractura cohesiva (despegado) y la fractura frágil (quiebre) en un único marco consistente. Sugiere que no tenemos que elegir entre estos diferentes tipos de falla; son solo capítulos diferentes en la misma historia de cómo se rompen los materiales. Para ingenieros y científicos, esto significa que pronto podríamos ser capaces de simular fallas complejas en puentes, aviones o tejidos biológicos con un nivel de detalle y precisión que antes estaba fuera de nuestro alcance.

Resumen Técnico: Regularización de campo de fase por degradación de resistencia de la cohesión: el caso de plano antiplano

Planteamiento del Problema
Los modelos tradicionales de fractura de campo de fase, aunque exitosos en la regularización de la fractura frágil, sufren una limitación fundamental: la superficie de resistencia del material no es un parámetro constitutivo independiente. En su lugar, emerge de la interacción entre la rigidez elástica, la tenacidad a la fractura y la longitud de regularización ( $\ell$ ). Este acoplamiento obliga a que la longitud de regularización actúe como una propiedad del material para ajustar las cargas de nucleación, lo que hace que este enfoque sea poco práctico para materiales rígidos o estructuras a gran escala donde $\ell$ debe ser resuelto por la malla. Además, la forma de la superficie de resistencia resultante está dictada por la forma de la energía elástica, fallando al intentar reproducir datos de resistencia multiaxial arbitrarios. Por el contrario, mientras que la teoría de la plasticidad permite dominios de resistencia convexos arbitrarios, los modelos de plasticidad de ablandamiento estándar sufren patologías de localización y dependencia de la malla, y los modelos cohesivos de interfaz nítida presentan dificultades en la implementación numérica y la falta de un marco unificado para la nucleación de grietas a lo largo de trayectorias desconocidas.

Metodología
Los autores analizan un modelo de campo de fase variacional introducido en un trabajo complementario (Bourdin, Marigo, et al., 2025), diseñado específicamente para el entorno de cizalladura de plano antiplano. La innovación central de este modelo es la degradación de la resistencia en lugar de la rigidez.

Funcional de Energía: El modelo utiliza un funcional de energía donde la densidad de energía elástica se define mediante una minimización sobre un campo de deformación no lineal $p$ . Esta formulación incorpora un dominio de resistencia convexo $K_0$ (definido por la resistencia al corte $\tau_c$ ) y una variable de daño $\alpha$ que degrada el dominio de resistencia ( $K(\alpha) = k(\alpha)K_0$ ) mientras deja inalterada la rigidez elástica ( $\mu$ ).
Estructura Variacional: El modelo se interpreta como una regularización de campo de fase de la plasticidad de ablandamiento (estilo Hencky) sin una condición de irreversibilidad sobre la deformación plástica $p$ en la formulación primaria. La variable de daño $\alpha$ dirige la evolución, gobernada por un término de regularización de tipo gradiente.
Esquema Numérico: Para resolver el problema de minimización de energía no suave y separadamente convexo, los autores proponen un esquema de minimización alternante. Los subproblemas se reestructuran como Programas de Cono de Segundo Orden (SOCP) y se resuelven mediante optimización cónica (MOSEK), evitando el suavizado o la penalización del término de resistencia.
Marco Analítico: El estudio deriva soluciones de forma cerrada para un problema de cizalladura simple, estableciendo la ley cohesiva equivalente y analizando las condiciones para el snap-back tanto en regímenes homogéneos como localizados.

Contribuciones Clave y Resultados

Parámetros de Material Independientes: El análisis demuestra que, en el límite $\ell \to 0$ , el modelo recupera un comportamiento de fractura cohesiva donde la resistencia ( $\tau_c$ ) y la tenacidad ( $G_c$ ) son datos de material independientes, desacoplados de la longitud de regularización $\ell$ y de la rigidez elástica $\mu$ . La longitud de regularización actúa puramente como un parámetro numérico.
Marco Variacional Unificado: El modelo unifica el análisis de límite, la plasticidad perfecta, la fractura cohesiva y la fractura frágil en un solo marco. La transición entre estos regímenes está gobernada por la relación entre el tamaño estructural y la longitud elasto-cohesiva ( $\ell_{ch} = \mu G_c / \tau_c^2$ ) y la relación de regularización $\ell/\ell_{ch}$ .
Soluciones de Forma Cerrada: Para el problema de cizalladura de plano antiplano, los autores derivan expresiones explícitas para:
- La ley cohesiva equivalente, que es independiente de $\ell$ y $\mu$ .
- Las condiciones para la ausencia de snap-back, que dependen de la relación de fragilidad $L/\ell_{ch}$ para la respuesta estructural y $\ell/\ell_{ch}$ para la respuesta del material.
- La densidad de energía superficial, mostrando que es una energía cohesiva de tipo Barenblatt con una pendiente inicial $\tau_c$ y una meseta $G_c$ .
Verificación Numérica:
- Cizalladura Simple: Las simulaciones verifican la recuperación de la ley cohesiva analítica y la independencia de la resistencia máxima respecto a $\ell$ . El modelo captura correctamente el efecto de tamaño, transitando de una falla progresiva (barras cortas) a un snap-back frágil (barras largas) a medida que aumenta $L/\ell_{ch}$ .
- Simulación de Surfeo (Surfing): Al simular una grieta propagándose bajo un factor de intensidad de esfuerzo prescrito, los autores distinguen entre deformación no lineal reversible e irreversible. En el caso irreversible, se forma un rastro plástico residual, aumentando la tenacidad efectiva ( $G_{eff}^c \approx 1.3 G_c$ ) en comparación con la $G_c$ prescrita.
- Entalla en V: El modelo une la plasticidad de zona de fluencia pequeña (SSY), la fractura cohesiva y la fractura frágil. Reproduce la transición de una zona plástica confinada (Hult-McClintock) a una grieta cohesiva (Dugdale/Barenblatt) y, finalmente, a una grieta frágil. Las cargas de nucleación para varios ángulos de entalla colapsan en una "ley universal" que interpola entre el límite de resistencia y el límite de Griffith.

Significancia y Reivindicaciones
El artículo afirma que esta formulación resuelve el problema de larga data de la regularización de la resistencia en los modelos de campo de fase. Al degradar la resistencia en lugar de la rigidez, el modelo permite:

Superficies de Resistencia Arbitrarias: La capacidad de prescribir dominios de resistencia convexos arbitrarios, una característica previamente exclusiva de la teoría de la plasticidad pero difícil de implementar en la fractura.
Independencia de la Malla: La ley cohesiva equivalente y la resistencia del material son independientes de la longitud de regularización, permitiendo que $\ell$ se elija únicamente para la resolución numérica.
Fenomenología Unificada: El modelo proporciona una descripción variacional única que abarca naturalmente todo el espectro de modos de falla —desde la localización plástica de tipo dúctil hasta la fractura frágil— sin requerir hipótesis a priori o el cambio entre diferentes modelos constitutivos.

Los autores señalan que, si bien el entorno de plano antiplano simplifica el análisis al asegurar la compatibilidad de saltos, la extensión a la elasticidad vectorial 3D implica condiciones de compatibilidad de saltos tensoriales que son objeto de trabajo continuo. Además, aunque el modelo maneja el comportamiento de tipo dúctil mediante la irreversibilidad plástica, una teoría completa de la fractura dúctil requeriría un mayor desarrollo respecto al endurecimiento y las deformaciones plásticas permanentes.

El nuevo enfoque de "Resistencia Primero"

Las tres etapas de la ruptura

El experimento de "Surfear" y la estela plástica

Cerrando la brecha con una muesca en V

Por qué esto es importante

Más como este