Analytic theory of shear localization in amorphous solids confined by Couette geometry

Resumen Técnico: Teoría Analítica de la Localización de Cizalla en Sólidos Amorfos

Planteamiento del Problema
La localización de cizalla y las bandas de cizalla son fenómenos prevalentes en sólidos amorfos bajo deformación de cizalla, lo que a menudo conduce a la falla catastrófica en materiales como los vidrios metálicos. Si bien está establecido que estos fenómenos surgen de eventos plásticos (avalanchas) en lugar de una elasticidad pura, los marcos teóricos actuales carecen de la capacidad para predecir o describir los campos de desplazamiento detallados asociados con las caídas de tensión durante estos eventos. Específicamente, no existe una teoría actual que pueda proporcionar predicciones a priori del campo de desplazamiento (movimiento de partículas) durante las avalanchas plásticas en un sólido amorfo cargado cuasiestáticamente confinado por una geometría de Couette. Además, aunque la banda de cizalla (eventos que atraviesan el sistema en el punto de fluencia) es distinta, los mecanismos que gobiernan la localización de cizalla a menor escala, que ocurre tanto antes como después de la fluencia, requieren un enfoque teórico dedicado.

Metodología
Los autores combinan simulaciones numéricas con una teoría analítica recién derivada para abordar esta brecha.

• Simulaciones Numéricas: El estudio utiliza un modelo bidimensional de un sólido amorfo compuesto por $N=3534$ discos repulsivos bidispersos (radios 1 y 1.4) confinados entre dos anillos concéntricos ($R_{in}=28$, $R_{out}=80.8$). El sistema se prepara a una fracción de área $\phi=0.91$ y se lleva al equilibrio mecánico. Se aplica un protocolo de carga cuasiestática donde el anillo interior se rota en sentido antihorario en pasos pequeños ($\delta\theta = 0.024^\circ$), y el sistema se reequilibra mediante dinámica molecular amortiguada tras cada paso. El límite exterior permanece fijo. El estudio se centra en los campos de desplazamiento asociados con caídas bruscas de tensión (eventos plásticos) en lugar de la respuesta elástica.
• Teoría Analítica: Los autores extienden marcos teóricos previos que tratan los eventos plásticos como cuadrupolos y dipolos topológicos en el campo de desplazamiento.
  • Formulan un Lagrangiano ($U$) que incluye el costo de los cuadrupolos, las interacciones cuadrupolo-cuadrupolo y el acoplamiento entre el campo de deformación y los campos cuadrupolares.
  • Crucialmente, la teoría tiene en cuenta un estrés de fondo sustancial ($\Sigma_{r\theta}$) acumulado antes del evento plástico, un factor que no fue abordado completamente en estudios previos de carga de un solo paso.
  • Al minimizar el Lagrangiano con respecto al campo cuadrupolar ($Q$) y al campo de desplazamiento ($d$), los autores derivan un tensor elástico renormalizado.
  • La inclusión de contribuciones dipolares (derivadas de los gradientes en el campo cuadrupolar) conduce a una ecuación diferencial gobernante para el campo de desplazamiento: $\Delta d + (1 + \tilde{\lambda})\nabla(\nabla \cdot d) + \Gamma d = 0$.
  • Esta ecuación se resuelve analíticamente utilizando expansiones de series de Fourier para las componentes radial ($d_r$) y tangencial ($d_\theta$), sujetas a las condiciones de contorno específicas de la geometría de Couette.

Contribuciones Clave y Resultados

1. Predicción Analítica de Campos de Desplazamiento: La principal contribución es la derivación de una solución analítica (Ec. 14) para el campo de desplazamiento tangencial promediado angularmente $d_\theta(r)$. Esta solución, expresada en términos de funciones de Bessel ($J_1$ y $Y_1$) y un parámetro de apantallamiento $\kappa_e$, reproduce con éxito los complejos perfiles de desplazamiento observados en las simulaciones, incluyendo la reversión de la dirección de rotación de las partículas (sentido horario vs. antihorario) que la teoría de elasticidad clásica no logra predecir.
2. Mecanismo de Localización de Cizalla: La teoría demuestra que la localización de cizalla está determinada por valores discretos del parámetro de apantallamiento $\kappa_e$. Estos valores están dictados únicamente por los parámetros geométricos ($R_{in}$ y $R_{out}$) del sistema, no por el estado específico del material.
   • Cuando $\kappa_e \approx 0.067$, la localización de cizalla ocurre adyacente al límite interno.
   • Cuando $\kappa_e \approx 0.123$, la localización de cizalla ocurre en el seno del material (bulk).
   • El sistema selecciona preferentemente valores de $\kappa_e$ donde el denominador de la solución analítica tiende a cero (singularidades), lo que corresponde a una respuesta máxima del sistema.
3. Validación: Las predicciones analíticas muestran una excelente concordancia con las simulaciones numéricas. La teoría predice correctamente la ubicación y la forma de la localización de cizalla para la mayoría de las 200 caídas de tensión más grandes analizadas. El artículo señala que las caídas de tensión grandes tienden a localizarse cerca del límite, mientras que las más pequeñas se localizan en el seno del material, lo cual es consistente con la selección de diferentes modos de $\kappa_e$.
4. Distinción de la Banda de Cizalla: Los autores aclaran que su teoría lineal (expandida a orden cuadrático en el Lagrangiano) se aplica a pasos de cizalla cuasiestáticos donde el estrés de fondo es divergencia-nula. Se establece explícitamente que esta teoría no describe la banda de cizalla asociada con la fluencia, la cual requiere una extensión no lineal (por ejemplo, añadiendo términos de dipolo de cuarto orden) que queda fuera del alcance de este trabajo.

Significancia
El artículo afirma proporcionar la primera teoría analítica capaz de predecir el campo de desplazamiento promediado angularmente asociado con las caídas plásticas en sólidos amorfos sometidos a cizalla. Al incorporar cargas topológicas (cuadrupolos y dipolos) y tener en cuenta el estrés de fondo, la teoría cierra la brecha entre la elasticidad pura (que falla al describir estos fenómenos) y las complejas simulaciones numéricas. La capacidad de predecir a priori dónde ocurrirá la localización de cizalla basándose en restricciones geométricas ofrece una comprensión fundamental de la respuesta dúctil en materiales amorfos. Los autores enfatizan que, si bien la teoría explica la localización en sólidos dúctiles, el mecanismo para la banda de cizalla frágil sigue siendo un tema para futuras extensiones no lineales.