Computation of stresses in jammed packings modeled with Tresca friction

Résumé Technique : Calcul des contraintes dans les empilements de blocage modélisés par la friction de Tresca

Énoncé du problème
Cet article traite du calcul des champs de contraintes au sein d'empilements de cellules polygonales rigides en état de blocage (jammed), tels que les structures de maçonnerie ou les matériaux granulaires denses. Contrairement aux méthodes classiques d'éléments discrets (DEM) qui reposent souvent sur des déformations élastiques pour pénaliser l'interpénétration, ce travail se concentre sur des systèmes où les cellules sont entièrement rigides et ne se déforment pas élastiquement. Dans de telles configurations de « blocage », l'équilibre mécanique est maintenu exclusivement par les forces de contact et la friction, sans déformation élastique significative. Le défi spécifique réside dans la détermination de la distribution de contrainte interne dans ces systèmes hyperstatiques, où le nombre de variables de force excède le nombre de contraintes d'équilibre, conduisant à un espace de solutions non unique. Les auteurs modélisent la physique du contact en utilisant une loi de friction de Tresca, qui postule un seuil fixe pour le glissement tangentiel, distinct de la friction de Coulomb dépendante de la contrainte.

Méthodologie
La méthodologie proposée suit un cadre de calcul en trois étapes :

1. Minimisation de l'énergie primale :
   Les auteurs formulent un problème de minimisation contraint où l'énergie frictionnelle totale est minimisée. Le fonctionnel d'énergie consiste en la somme du travail de la friction de Tresca sur toutes les interfaces internes, définie par $s_T |[u]_e \cdot t_e|$, où $s_T$ est le coefficient de friction, $[u]_e$ est le saut de déplacement, et $t_e$ est le vecteur tangent. Cette minimisation est soumise à :

• Contraintes de non-interpénétration : La composante normale du saut de déplacement doit être non négative ($[u]_e \cdot n_e \geq 0$).
• Confinement : Une condition aux limites de Neumann ($g$) est appliquée pour empêcher les cellules de glisser vers l'infini, car l'énergie ne possède pas de forces de rappel élastiques.
• Invariance par translation : Une contrainte globale est ajoutée pour fixer le centre de masse.
  Les cellules sont modélisées comme des corps rigides possédant uniquement des degrés de liberté de translation ; les degrés de liberté de rotation sont explicitement négligés. Cette simplification est motivée par l'observation que, dans les configurations de blocage, la résistance à la rotation (nécessitant de la dilatance) est de plusieurs ordres de grandeur supérieure à la résistance au glissement.

2. Maximisation duale :
   Puisque le problème primal recherche les déplacements (qui sont effectivement nuls dans un état de blocage stable), les auteurs en déduisent un problème de maximisation dual. Les variables duales correspondent aux forces agissant sur les arêtes internes. La résolution de ce problème dual permet d'obtenir les contraintes normales et tangentielles aux interfaces. L'existence de minimiseurs pour l'énergie primale et de solutions pour le problème dual est établie sous des hypothèses appropriées concernant le terme de confinement. La formulation duale sélectionne efficacement une solution spécifique parmi le manifold de contraintes possibles de l'état hyperstatique.

3. Reconstruction de la contrainte :
   Le résultat du problème dual fournit des valeurs de contrainte discrètes aux arêtes. Pour reconstruire un champ de contrainte continu $\sigma$ au sein de chaque cellule, les auteurs utilisent des éléments finis de type Raviart–Thomas de bas ordre ($RT_0$). Ces éléments sont conformes à $H(\text{div})$, ce qui signifie que leurs degrés de liberté correspondent aux contraintes normales aux barycentres des arêtes. La reconstruction implique la résolution d'un problème mixte local dans chaque cellule pour trouver un tenseur de contrainte $\sigma_c$ tel que :

• $\text{div}(\sigma_c) = 0$ (équilibre en l'absence de charges volumiques).
• La contrainte normale sur la frontière correspond aux valeurs obtenues de la solution duale.
• La symétrie du tenseur de contrainte ($\sigma = \sigma^T$) est imposée faiblement via un multiplicateur de Lagrange. Ceci est nécessaire car le modèle discret sous-jacent néglige les équations d'équilibre des moments, ce qui signifie que la contrainte reconstruite peut posséder une composante antisymétrique représentant des couples résiduels.

Contributions clés

• Modèle de friction rigide : L'article introduit un modèle pour les empilements de blocage qui repose uniquement sur la friction rigide (critère de Tresca) sans paramètres élastiques, se distinguant ainsi des approches DEM standards.
• Gestion de l'hyperstaticité : Le travail traite explicitement de l'indétermination statique des empilements de blocage. Il démontre que l'approche d'optimisation sélectionne une solution unique dans l'espace des états d'auto-contrainte en minimisant le fonctionnel d'énergie de friction.
• Reconstruction de contrainte cohérente : Une nouvelle technique de reconstruction est présentée en utilisant des éléments $RT_0$ pour générer un champ de contrainte cohérent à partir des forces d'interface discrètes, garantissant l'équilibre des forces locales au sein de chaque cellule.
• Négligence de la rotation : Les auteurs proposent un modèle simplifié ignorant les degrés de liberté de rotation, arguant que cela est physiquement justifié pour les empilements compressés où le glissement de translation domine, tout en reconnaissant que cela conduit à des tenseurs de contrainte non symétriques.

Résultats numériques
La méthodologie a été testée sur diverses mailles polygonales 2D (spécifiquement des tessellations de Voronoi) avec les conclusions suivantes :

• Vérification : Dans des tests avec des conditions aux limites uniformes, le champ de contrainte reconstruit correspondait à l'état de contrainte uniforme attendu, confirmant la cohérence de la méthode.
• Test de cisaillement : Dans un scénario de chargement par cisaillement, le tenseur de contrainte reconstruit présentait une asymétrie ($\sigma_{12} \neq \sigma_{21}$). Les auteurs attribuent cela à l'absence de bilan de moment imposé dans le modèle discret, notant que la partie antisymétrique rend compte des couples résiduels.
• Compression homogène : Malgré des conditions aux limites de compression uniforme, le solveur a produit un champ de contrainte non homogène. Les auteurs expliquent cela par une conséquence de l'hyperstaticité du système ; l'optimisation sélectionne un état spécifique minimisant l'énergie au sein du noyau (null space), ce qui conduit naturellement à la formation de « chaînes de force » (réseaux épars de contraintes élevées) et de régions de faible ou nulle contrainte.
• Tension en compression : Les résultats ont montré des cellules en tension (contraintes principales positives) même sous compression globale. Les auteurs identifient cela comme un phénomène physique découlant de l'imbrication frictionnelle et des effets de pontage, plutôt qu'un artefact numérique.

Signification et affirmations
L'article affirme que la méthodologie proposée offre une manière robuste et cohérente de calculer les contraintes dans les empilements de blocage en exploitant la physique de la friction rigide et la minimisation de l'énergie. Les auteurs soutiennent que leur approche capture mieux la physique sous-jacente des systèmes de blocage par rapport aux DEM basées sur l'élasticité, particulièrement dans les régimes où la déformation élastique est négligeable.

La signification des résultats réside dans la démonstration que la minimisation de l'énergie de friction (un objectif de type $L_1$) conduit naturellement à des réseaux de force épars, ce qui concorde avec les observations expérimentales des chaînes de force dans les matériaux granulaires. L'article reconnaît modestement ses limites : la négligence de l'équilibre de rotation produit des tenseurs de contrainte non symétriques, et l'utilisation de la friction de Tresca (seuil fixe) simplifie le problème par rapport aux inégalités variationnelles quasi-convexes plus complexes requises pour la friction de Coulomb. Les auteurs suggèrent que des travaux futurs pourraient étendre le modèle pour inclure les degrés de liberté de rotation et des lois de friction plus réalistes afin d'atteindre une pleine cohérence mécanique.