Computation of stresses in jammed packings modeled with Tresca friction

Resumen Técnico: Computación de Esfuerzos en Empaquetamientos Atascados Modelados con Fricción de Tresca

Planteamiento del Problema
Este artículo aborda la computación de campos de esfuerzos dentro de empaquetamientos atascados (jammed) de celdas poligonales rígidas, tales como estructuras de mampostería o materiales granulares densos. A diferencia de los Métodos de Elementos Discretos (DEM) tradicionales que a menudo dependen de deformaciones elásticas para penalizar la interpenetración, este trabajo se centra en sistemas donde las celdas son completamente rígidas y no se deforman elásticamente. En tales configuraciones "atascadas", el equilibrio mecánico se mantiene exclusivamente mediante fuerzas de contacto y fricción, sin una deformación elástica significativa. El desafío específico radica en determinar la distribución de esfuerzos internos en estos sistemas estáticamente indeterminados (hiperestáticos), donde el número de variables de fuerza excede el número de restricciones de equilibrio, lo que conduce a un espacio de soluciones no único. Los autores modelan la física de contacto utilizando una ley de fricción de Tresca, la cual postula un umbral fijo para el deslizamiento tangencial, distinto de la fricción de Coulomb dependiente del esfuerzo.

Metodología
La metodología propuesta sigue un marco computacional de tres etapas:

1. Minimización de Energía Primal:
   Los autores formulan un problema de minimización restringido donde se minimiza la energía friccional total. El funcional de energía consiste en la suma del trabajo de fricción de Tresca sobre todas las interfaces internas, definida como $s_T |[u]_e \cdot t_e|$, donde $s_T$ es el coeficiente de fricción, $[u]_e$ es el salto de desplazamiento, y $t_e$ es el vector tangente. Esta minimización está sujeta a:

   • Restricciones de no interpenetración: El componente normal del salto de desplazamiento debe ser no negativo ($[u]_e \cdot n_e \geq 0$).
   • Confinamiento: Se aplica una condición de contorno de Neumann ($g$) para evitar que las celdas se deslicen hacia el infinito, ya que la energía carece de fuerzas de restauración elásticas.
   • Invarianza de traslación: Se añade una restricción global para fijar el centro de masa.
     Las celdas se modelan como cuerpos rígidos con grados de libertad de traslación únicamente; los grados de libertad rotacionales se descartan explícitamente. Esta simplificación está motivada por la observación de que la resistencia a la rotación (que requiere dilatancia) es órdenes de magnitud superior a la resistencia al deslizamiento.

2. Maximización Dual:
   Dado que el problema primal busca desplazamientos (que son efectivamente cero en un estado atascado estable), los autores deducen un problema de maximización dual. Las variables duales corresponden a las fuerzas que actúan en los bordes internos. La resolución de este problema dual proporciona los esfuerzos normales y tangenciales en las interfaces. La existencia de minimizadores para la energía primal y soluciones para el problema dual se establece bajo hipótesis apropiadas respecto al término de confinamiento. La formulación dual selecciona efectivamente una solución específica del espacio de la variedad hiperestática de posibles estados de esfuerzo.

3. Reconstrucción de Esfuerzos:
   El resultado del problema dual proporciona valores de esfuerzo discretos en los bordes. Para reconstruir un campo de esfuerzos continuo $\sigma$ dentro de cada celda, los autores emplean elementos finitos de Raviart–Thomas de bajo orden ($RT_0$). Estos elementos son conformes a $H(\text{div})$, lo que significa que sus grados de libertad corresponden a los esfuerzos normales en los baricentros de los bordes. La reconstrucción implica resolver un problema mixto local en cada cela para encontrar un tensor de esfuerzo $\sigma_c$ tal que:

   • $\text{div}(\sigma_c) = 0$ (equilibrio en ausencia de cargas volumétricas).
   • El esfuerzo normal en el contorno coincide con los valores obtenidos de la solución dual.
   • La simetría del tensor de esfuerzo ($\sigma = \sigma^T$) se impone débilmente mediante un multiplicador de Lagrange. Esto es necesario porque el modelo discreto subyacente descuida las ecuaciones de balance de momentos, lo que significa que el esfuerzo reconstruido puede poseer un componente antisimétrico que representa torques residuales.

Contribuciones Clave

• Modelo de Fricción Rígida: El artículo introduce un modelo para empaquetamientos atascados que depende únicamente de la fricción rígida (criterio de Tresca) sin parámetros elásticos, distinguiéndose de los enfoques DEM estándar.
• Manejo de la Hiperestaticidad: El trabajo aborda explícitamente la indeterminación estática de los empaquetamientos atascados. Demuestra que el enfoque de optimización selecciona una solución única del espacio de estados de auto-esfuerzo al minimizar el funcional de energía friccional.
• Reconstrucción de Esfuerzos Consistente: Se presenta una técnica de reconstrucción novedosa utilizando elementos $RT_0$ para generar un campo de esfuerzos consistente a partir de las fuerzas de interfaz discretas, asegurando el equilibrio de fuerzas local dentro de cada celda.
• Omisión de la Rotación: Los autores proponen un modelo simplificado que ignora los grados de libertad rotacionales, argumentando que esto está físicamente justificado para empaquetamientos comprimidos donde el deslizamiento traslacional domina, reconociendo al mismo tiempo que esto conduce a tensores de esfuerzo no simétricos.

Resultados Numéricos
La metodología fue probada en varias mallas poligonales 2D (específicamente teselaciones de Voronoi) con los siguientes hallazgos:

• Verificación: En pruebas con condiciones de contorno uniformes, el campo de esfuerzos reconstruido coincidió con el estado de esfuerzo uniforme esperado, confirmando la consistencia del método.
• Prueba de Cizalla: En un escenario de carga de cizalla, el tensor de esfuerzo reconstruido exhibió asimetría ($\sigma_{12} \neq \sigma_{21}$). Los autores atribuyen esto a la falta de balance de momentos impuesto en el modelo discreto, señalando que la parte antisimétrica da cuenta de los torques residuales.
• Compresión Homogénea: A pesar de las condiciones de contorno de compresión uniforme, el resolvedor produjo un campo de esfuerzos no homogéneo. Los autores explican esto como una consecuencia de la hiperestaticidad del sistema; la optimización selecciona un estado específico de minimización de energía del espacio nulo, lo que naturalmente conduce a la formación de "cadenas de fuerza" (redes dispersas de alto esfuerzo) y regiones de bajo o cero esfuerzo.
• Tensión en Compresión: Los resultados mostraron celdas en tensión (esfuerzos principales positivos) incluso bajo compresión global. Los autores identifican esto como un fenómeno físico derivado del trabazón friccional (frictional interlocking) y efectos de puenteo, más que un artefacto numérico.

Significancia y Reivindicaciones
El artículo afirma que su metodología propuesta ofrece una forma robusta y consistente de computar esfuerzos en empaquetamientos atascados aprovechando la física de la fricción rígida y la minimización de la energía. Los autores argumentan que su enfoque captura mejor la física subyacente de los sistemas atascados en comparación con los DEM basados en elasticidad, particularmente en regímenes donde la deformación elástica es despreciable.

La significancia de los resultados radica en la demostración de que la minimización de la energía friccional (un objetivo de tipo $L_1$) conduce naturalmente a redes de fuerza dispersas, alineándose con las observaciones experimentales de cadenas de fuerza en materiales granulares. El artículo reconoce modestamente sus limitaciones: la omisión del equilibrio rotacional resulta en tensores de esfuerzo no simétricos, y el uso de la fricción de Tresca (umbral fijo) simplifica el problema en comparación con las desigualdades de cuasi-variación no convexas requeridas para la fricción de Coulomb. Los autores sugieren que trabajos futuros podrían extender el modelo para incluir grados de libertad rotacionales y leyes de fricción más realistas para lograr una plena consistencia mecánica.