Peridynamic modeling of the crack velocity dependence via an incubation time fracture criterion

Este estudio utiliza un enfoque peridinámico con un criterio de fractura basado en tiempo de incubación para modelar los experimentos de Ravi-Chandar y Knauss sobre Homalite-100, revelando que las variaciones en el factor de intensidad de tensión del Modo I a velocidades de grieta constantes y el inicio de la ramificación microscópica a velocidades más altas proporcionan nuevas perspectivas sobre la naturaleza de la dependencia de la velocidad de la grieta en la fractura dinámica.

Lo esencial

Imagina que estás observando una grieta correr a través de un trozo de plástico frágil, como una lámina de Homalite-100. En los viejos tiempos de la física, los científicos pensaban que si conocías la velocidad a la que se movía la grieta, podías calcular exactamente cuánta "tensión" (o presión) la empujaba hacia adelante. Era como pensar que si un coche va a 60 millas por hora, el motor debe estar produciendo exactamente 100 caballos de fuerza. ¿Simple, verdad?

Pero los experimentos de la década de 1980 demostraron que esto no era cierto. A veces, la grieta se movía a exactamente la misma velocidad, pero la presión que la empujaba era enormemente diferente. Era como si dos coches fueran ambos a 60 millas por hora, pero uno tuviera un motor diminuto y el otro un propulsor de cohete. Los científicos estaban perplejos: ¿Por qué la misma velocidad tiene diferentes "empujones"?

Este artículo es una historia de detectives donde los autores utilizan un nuevo tipo de simulación por computadora para resolver este misterio.

La herramienta del detective: Peridinámica

La mayoría de los modelos informáticos de grietas son como una cadena de fichas de dominó. Si cae una ficha, empuja a la siguiente. Pero si falta una ficha (una grieta), la cadena se rompe y las matemáticas se atascan.

Los autores utilizaron un método llamado Peridinámica. Piensa en esto no como una cadena, sino como un enjambre de abejas. Cada abeja puede hablar con todas las demás abejas dentro de cierta distancia, incluso si hay un hueco en medio. Si una abeja se aleja (se forma una grieta), las otras abejas simplemente dejan de hablarle, pero el resto del enjambre sigue moviéndose perfectamente bien. Esto permite que la computadora maneje la rotura y la grieta sin confundirse.

El ingrediente secreto: El "tiempo de incubación"

El verdadero avance en este artículo es cómo decidieron cuándo una grieta debería romperse realmente.

A la antigua, si la presión se volvía lo suficientemente alta, el material se rompía instantáneamente. Pero los autores utilizaron una regla llamada Criterio de Tiempo de Incubación.

Imagina que estás intentando romper una ramita seca. No solo tiras y se rompe instantáneamente. Tiras, la mantienes ahí por una fracción de segundo mientras las fibras se estiran y se debilitan, y entonces se rompe. Esa fracción de segundo es el "tiempo de incubación".

Los autores programaron su enjambre de computadoras para recordar los últimos microsegundos de presión. El material solo se rompe si la presión promedio durante ese corto período de "incubación" es lo suficientemente alta. Esto tiene en cuenta el hecho de que los materiales necesitan un pequeño lapso de tiempo para "decidir" romperse.

Lo que descubrieron

Realizaron simulaciones de las placas de plástico siendo separadas, tal como en los experimentos reales. Esto es lo que descubrieron:

1. El acertijo de la velocidad frente a la presión: Al igual que en los experimentos reales, su computadora mostró que para la misma velocidad de grieta, la presión (Factor de Intensidad de Tensión) no era un solo número. Era un rango. A veces era baja, a veces alta.
2. El efecto de "micro-ramificación": Cuando la grieta se movía lentamente, iba recta. Pero cuando aceleraba (por encima de 400 metros por segundo), empezaba a volverse nerviosa. Comenzó a brotar pequeñas grietas laterales microscópicas, como una rama de árbol que se divide en ramitas.
   • La analogía: Imagina a un corredor haciendo un sprint. A un trote constante, corre en línea recta. Pero cuando sprinta a máxima velocidad, empieza a tambalearse y zigzaguear ligeramente para mantener el equilibrio.
   • El resultado: Estos pequeños "tambaleos" (micro-ramificaciones) hicieron que la lectura de la presión subiera y bajara salvajemente. Esto explicó por qué la presión no era única para una velocidad dada; la grieta estaba cambiando físicamente su forma ligeramente mientras corría.

La conclusión

El artículo concluye que la razón por la que vemos diferentes valores de presión para la misma velocidad de grieta es porque la grieta no es una línea suave y perfecta. Es una cosa caótica y viva que fluctúa.

• A velocidades más bajas: La grieta es estable y la presión es relativamente constante.
• A velocidades más altas: La grieta comienza a "micro-ramificarse" (brotar pequeñas grietas laterales). Este caos hace que la presión rebote, creando la dispersión observada en los experimentos.

Al utilizar este "enjambre de abejas" (Peridinámica) combinado con el "período de espera" (Tiempo de Incubación), los autores recrearon con éxito la relación desordenada y no única entre la velocidad de la grieta y la presión que los experimentos del mundo real habían mostrado durante décadas. Demostraron que el "ruido" en los datos no es un error; es una característica física real de cómo se comportan las grietas que se mueven rápidamente.

Resumen técnico

Resumen Técnico: Modelado Peridinámico de la Dependencia de la Velocidad de la Grieta mediante un Criterio de Fractura por Tiempo de Incubación

Planteamiento del Problema
Este estudio aborda un problema central no resuelto en la mecánica de fractura dinámica: la dependencia del factor de intensidad de tensiones (FIT) instantáneo en modo I ($K_I$) respecto a la velocidad de propagación de la grieta ($\dot{a}$). La Mecánica de Fractura Elástica Lineal (MFEL) clásica y los criterios dinámicos basados en Freund y Rosakis a menudo sugieren una relación única entre $K_I$ y $\dot{a}$. Sin embargo, los datos experimentales de Ravi-Chandar y Knauss sobre placas de polímero Homalite-100 frágil demuestran una dependencia no única, donde una sola velocidad de grieta corresponde a un rango de valores de FIT. Además, a velocidades más altas, se produce ramificación microscópica, lo que genera una dispersión significativa en las mediciones de FIT. El desafío radica en reproducir numéricamente esta no unicidad y las fluctuaciones asociadas sin depender de propiedades materiales dependientes de la tasa ajustadas empíricamente.

Metodología
Los autores emplean un marco peridinámico no local para simular la propagación dinámica de grietas en Homalite-100, replicando las condiciones experimentales de Ravi-Chandar y Knauss. La innovación metodológica central es la integración del Criterio de Fractura por Tiempo de Incubación (ITFC) en la formulación peridinámica.

• Marco Peridinámico: Las simulaciones utilizan un modelo peridinámico de estado basado en estados ordinarios implementado en el software de código abierto Peridigm. El material se discretiza en puntos que interactúan dentro de un horizonte esférico ($\delta$), permitiendo el manejo natural de discontinuidades como grietas sin operadores diferenciales espaciales.
• Criterio de Fractura por Tiempo de Incubación (ITFC): A diferencia de los criterios clásicos que asumen una falla instantánea, el ITFC postula que la fractura requiere un tiempo finito ($\tau$) para que ocurran procesos preparatorios (por ejemplo, microfisuración). El criterio evalúa el esfuerzo promedio en el tiempo durante el período de incubación $\tau$ y una longitud característica $d$.
  • La longitud característica $d$ se deriva del criterio de Irwin cuasiestático para garantizar la consistencia con el FIT estático crítico ($K_{Ic}$).
  • La implementación peridinámica adapta el criterio de fractura por Esfuerzo Remoto (RS). La falla de los enlaces se determina integrando el esfuerzo normal del enlace sobre el tiempo de incubación $\tau$ (Ecuación 11).
  • Para gestionar la complejidad computacional y centrarse en la ramificación microscópica en lugar de la ramificación macroscópica completa, el criterio se evalúa únicamente en planos paralelos ($0^\circ$) y perpendiculares ($90^\circ$) al eje principal de la grieta.
• Configuración de la Simulación: El modelo simula una probeta de 300 mm $\times$ 500 mm con una entalla previa de 50 mm. La carga dinámica se aplica mediante un pulso de fuerza volumétrica trapezoidal con un tiempo de subida de 25 $\mu$s. El FIT dinámico se calcula utilizando un contorno rectangular móvil para la integral J dinámica, asegurando que la trayectoria de integración encierre la punta principal de la grieta mientras excluye las ramas microscópicas en desarrollo.

Resultados Clave
Las simulaciones numéricas reproducen con éxito los fenómenos observados experimentalmente respecto a la relación $K_I$-$\dot{a}$:

1. Dependencia No Única $K_I$-$\dot{a}$: Los resultados confirman que, para una velocidad de propagación de grieta casi constante, el FIT instantáneo exhibe fluctuaciones significativas. En consecuencia, el modelo predice un rango de valores de FIT para una velocidad dada, reflejando la no unicidad observada en los experimentos.
2. Estabilización de la Velocidad y Transitorios: Mientras que los experimentos muestran que la velocidad de la grieta se estabiliza casi instantáneamente, las simulaciones exhiben una fase inicial de aceleración que dura aproximadamente 30–40 $\mu$s antes de alcanzar un estado estacionario. Durante esta fase transitoria, la historia del FIT muestra variaciones pronunciadas.
3. Ramificación Microscópica y Dispersión: A tasas de carga más altas (resultando en velocidades > 400 m/s), las simulaciones exhiben el inicio de la ramificación microscópica. Este fenómeno se correlaciona con un aumento marcado en la amplitud de las oscilaciones del FIT. La dispersión numérica en los valores de FIT a altas velocidades (0.53–1.25 MPa$\sqrt{m}$) se alinea razonablemente bien con los rangos experimentales (0.45–1.5 MPa$\sqrt{m}$).
4. Validación del Mecanismo: El estudio demuestra que la dispersión observada en el FIT está físicamente justificada por la dinámica de la punta de la grieta, específicamente la fase de aceleración y el inicio de la ramificación microscópica, en lugar de ser un artefacto de errores de medición.

Significado y Afirmaciones
El artículo afirma que la implementación peridinámica del ITFC proporciona un enfoque valioso para evaluar la tenacidad a la fractura dinámica y modelar la propagación de grietas en un amplio rango de tasas de carga. El significado principal radica en la capacidad del modelo para capturar la sensibilidad a la tasa de la fractura frágil y la no unicidad de la relación $K_I$-$\dot{a}$ sin introducir parámetros materiales dependientes de la tasa arbitrarios.

Los autores notan modestamente que, aunque el modelo captura cualitativamente las tendencias experimentales, persisten discrepancias cuantitativas, particularmente en cuanto al momento de la estabilización de la velocidad y la magnitud precisa de las fluctuaciones del FIT. Atribuyen la predicción exitosa de la dependencia no única a la inclusión del tiempo de incubación, que cuenta con la historia temporal del esfuerzo, y a la naturaleza no local de la peridinámica, que maneja naturalmente los campos singulares y las discontinuidades asociadas con la fractura dinámica. El estudio concluye que la interacción entre la fase de aceleración de la grieta, las fluctuaciones del FIT y el inicio de la ramificación microscópica son los mecanismos gobernantes detrás del comportamiento $K_I$-$\dot{a}$ observado experimentalmente.