Analytic Nonlinear Theory of Shear Banding in Amorphous Solids

Cet article présente une théorie analytique non linéaire de la localisation de cisaillement dans les solides amorphes cisailés de manière athermale, en déduisant des équations qui prennent en compte l'écrantage des dipôles induit par la plasticité pour expliquer le mécanisme d'instabilité, prédire la largeur de la bande de cisaillement et déterminer le seuil de contrainte critique pour la rupture.

L'essentiel

Imaginez un bloc de verre ou un tas de sable. Dans le monde de la physique, on les appelle des « solides amorphes ». Contrairement à un cristal (comme un diamant), où les atomes sont alignés en rangées parfaites, les atomes de ces matériaux sont emmêlés de façon aléatoire, comme une foule de gens à un concert sans places assignées.

Pendant longtemps, les scientifiques ont tenté de prédire comment ces matériaux se briseraient ou se déformeraient en utilisant les mêmes règles que pour les cristaux parfaits. Mais ces règles ont échoué. Lorsque vous poussez sur du verre ou du sable, il ne se plie pas simplement ; il se brise soudainement ou forme une ligne étroite et tranchante de dommages appelée bande de cisaillement. Pensez-y comme à une fissure se formant dans un pare-brise, mais au lieu d'une simple ligne, c'est une zone où le matériau glisse sur lui-même.

Ce papier d'Avanish Kumar et Itamar Procaccia offre une nouvelle « recette » mathématique pour prédire exactement comment et pourquoi ces bandes de cisaillement se forment, et à quoi elles ressemblent. Voici l'explication en termes simples :

1. Le Problème : Le Désordre « Caché »

Lorsque vous poussez sur un cristal parfait, il s'étire de manière fluide. Mais lorsque vous poussez sur des solides amorphes, de minuscules réarrangements chaotiques se produisent à l'intérieur. Les auteurs appellent ces « événements plastiques ».

• L'Analogie : Imaginez une pièce bondée. Si vous poussez la foule, les gens ne se déplacent pas simplement en ligne droite ; ils se cognent les uns contre les autres, se décalent sur le côté et créent de petits tourbillons de mouvement. Dans l'article, ces tourbillons sont appelés « quadrupôles » (formes de mouvement à quatre points).
• L'Ancienne Théorie : Les théories précédentes traitaient ces tourbillons comme s'ils étaient répartis uniformément, comme du sucre dissous dans du thé. Cela fonctionnait pour de petites poussées, mais échouait à expliquer la formation soudaine et violente des bandes de cisaillement.
• La Nouvelle Insight : Les auteurs ont réalisé que lorsque le matériau est soumis à une contrainte, ces tourbillons cessent d'être répartis uniformément. Ils commencent à se regrouper, créant des « dipôles » (forces à deux points) qui agissent comme des charges d'écran.
  • Métaphore : Imaginez ces dipôles comme une foule de gens tenant des parapluies. S'ils sont répartis uniformément, la pluie (la contrainte) frappe tout le monde également. Mais s'ils se regroupent, ils créent un « bouclier » ou un « écran » qui bloque la pluie à certains endroits et la laisse s'abattre à d'autres. Cet écran crée une « échelle de longueur » spécifique — une largeur naturelle pour la zone de dommages.

2. La Grande Avancée : Mathématiques Non Linéaires

L'article soutient que pour comprendre les bandes de cisaillement, on ne peut pas utiliser des mathématiques simples et linéaires (équations linéaires). Il faut des mathématiques non linéaires.

• L'Analogie : Imaginez conduire une voiture. À basse vitesse, si vous tournez légèrement le volant, la voiture tourne légèrement (linéaire). Mais à haute vitesse, un tout petit tour de volant peut envoyer la voiture dans une embardée (non linéaire).
• Les auteurs ont dérivé un nouvel ensemble d'équations qui tiennent compte de ce comportement « haute vitesse » du matériau. Ils ont inclus deux effets non linéaires principaux :
  1. Comment la forme du matériau change au fur et à mesure qu'il se déforme (la relation contrainte-déplacement).
  2. Comment les « tourbillons » de mouvement interagissent entre eux lorsqu'ils sont entassés (les interactions dipolaires).

3. Le Résultat : Prédire la « Fissure »

En résolvant ces équations complexes, les auteurs ont trouvé un moyen de prédire le profil de la bande de cisaillement.

• Le Cas « Ductile » (Souple) : Dans les matériaux un peu plus flexibles, la bande de cisaillement est large et lisse.
  • Métaphore : Comme une pente lente et douce. Le matériau glisse progressivement sur une large zone. Les mathématiques prédisent que cette forme ressemble à une courbe tangente hyperbolique (tanh) — une forme de S lisse.
• Le Cas « Brittle » (Dur) : Dans les matériaux très rigides, la bande de cisaillement est incroyablement tranchante et étroite.
  • Métaphore : Comme le bord d'une falaise. Le matériau reste immobile d'un côté et glisse instantanément de l'autre. Les mathématiques montrent que dans ce cas, le « cœur » de la bande se comporte différemment des bords, créant une transition très nette.

4. Le Commutateur « d'Instabilité »

L'article explique aussi quand cela se produit.

• L'Analogie : Imaginez équilibrer un crayon sur sa pointe. Tant que le vent est léger, il reste debout. Mais à une vitesse de vent critique spécifique, il devient instable et tombe.
• Les auteurs ont calculé la « contrainte critique » exacte (la vitesse du vent) où le matériau perd sa stabilité. Ils ont découvert que cela se produit lorsqu'une valeur mathématique spécifique (une valeur propre du « Hessien », qui est simplement une façon élégante de dire un calculateur de stabilité) tombe à zéro.
• Une fois ce point atteint, le matériau ne peut plus maintenir sa forme uniformément, et la bande de cisaillement « se déclenche » pour exister.

5. Pourquoi Cela Compte (Selon l'Article)

Les théories précédentes pouvaient dire « une bande de cisaillement se formera », mais elles ne pouvaient pas vous dire quelle serait sa largeur ou à quoi ressemblerait sa forme sans simplement deviner ou utiliser des simulations informatiques.

• Cet article fournit une théorie analytique, ce qui signifie qu'il donne une formule directe.
• Il explique que la largeur de la bande de cisaillement est déterminée par une compétition entre la rigidité du matériau et l'effet d'« écran » de ces tourbillons internes.
• Il distingue les matériaux fragiles (ruptures nettes et soudaines) des matériaux ductiles (glissements lents et larges) sur la base des mathématiques de ces équations.

Résumé

En bref, les auteurs ont construit un nouveau modèle mathématique qui traite les solides amorphes (comme le verre ou le sable) non pas comme de simples ressorts, mais comme des foules complexes de particules en mouvement. En tenant compte de la façon dont ces particules « s'écrant » mutuellement leurs mouvements et de leur comportement non linéaire sous contrainte, ils ont dérivé une formule qui prédit exactement quand un matériau se brisera et à quoi ressemblera la « fissure » résultante (bande de cisaillement), depuis un glissement lisse jusqu'à une rupture nette.

Résumé technique

Résumé technique : Théorie analytique non linéaire du cisaillement localisé dans les solides amorphes

Énoncé du problème
La réponse mécanique des solides amorphes (par exemple, les verres, les milieux granulaires) sous cisaillement quasi-statique apériodique diffère fondamentalement de celle des solides cristallins. Alors que les cristaux présentent des réponses purement élastiques jusqu'à la rupture, les solides amorphes subissent des déformations plastiques à tout niveau de déformation, se manifestant généralement par des « inclusions d'Eshelby » quadrupolaires. Une instabilité critique dans ces matériaux est le cisaillement localisé (shear banding), où la déformation se concentre dans des bandes étroites. Ce phénomène varie depuis des sauts brusques, de type fragile, du déplacement à travers un plan, jusqu'à une localisation plus douce, de type ductile.

Les théories existantes présentent des limitations significatives :

1. Théories linéaires : Les approches analytiques antérieures reposaient sur des approximations linéaires entre la déformation et le déplacement, ainsi que sur des développements quadratiques des lagrangiens. Bien qu'elles aient réussi à décrire les effets d'écran et les modules élastiques renormalisés, elles ne pouvaient pas prédire le profil spécifique d'une bande de cisaillement ni la contrainte critique pour le début de l'instabilité.
2. Modèles élastoplastiques : Les modèles mésoscopiques standards capturent les statistiques d'avalanches et la redistribution des contraintes, mais sont souvent locaux dans l'espace. Ils manquent d'une échelle de longueur continue intrinsèque, ce qui signifie que la largeur d'une bande de cisaillement est déterminée par une régularisation numérique ou la taille des blocs plutôt que par des propriétés physiques du matériau. Par conséquent, ils ne peuvent pas déterminer de manière unique le profil de déformation interne.

Les auteurs visent à développer une théorie analytique non linéaire qui dérive les équations du champ de déplacement, prédit la contrainte accumulée critique pour l'instabilité, fournit une solution analytique pour le profil de la bande de cisaillement, et distingue les caractéristiques ductiles et fragiles en fonction des propriétés du matériau.

Méthodologie
Les auteurs construisent une théorie variationnelle basée sur une formulation lagrangienne qui intègre deux non-linéarités essentielles précédemment négligées :

1. Relation non linéaire déformation-déplacement : Au-delà de la relation linéaire $u_{\alpha\beta} = \frac{1}{2}(\partial_\alpha d_\beta + \partial_\beta d_\alpha)$, la théorie utilise la définition géométrique complète :
   $$u_{\alpha\beta} = \frac{1}{2} [\partial_\alpha d_\beta + \partial_\beta d_\alpha] + \frac{1}{2}\partial_\beta d_\gamma \partial_\alpha d_\gamma$$
2. Développement dipolaire d'ordre supérieur : Bien que les événements plastiques induisent des quadrupôles, les gradients dans la densité de quadrupôles créent des dipôles effectifs. Les auteurs développent le lagrangien pour inclure des termes quartiques dans le champ de dipôle ($P^\alpha$), nécessaires pour capturer la saturation non linéaire de l'instabilité.

La dérivation procède par les étapes suivantes :

• Dérivation des équations non linéaires : Partant d'un fonctionnel d'énergie incluant la contrainte de fond ($\Sigma_{\alpha\beta}$), l'énergie élastique, et les interactions quadrupôle/dipôle, les auteurs minimisent l'énergie par rapport au champ de déplacement ($d$) et au champ de quadrupôle ($Q$). Cela aboutit à une équation aux dérivées partielles non linéaire (Éq. 32) pour le champ de déplacement.
• Projection vers une équation d'amplitude scalaire : Pour résoudre l'équation vectorielle complexe, les auteurs emploient une stratégie basée sur le mode mou (soft mode) de la hessienne du système. Ils soutiennent que près de l'instabilité, le champ de déplacement peut être projeté sur une fonction scalaire $f(\xi)$ variant le long de la normale à la bande $\xi = \mathbf{m} \cdot \mathbf{x}$. Cela réduit le problème à une équation d'amplitude unidimensionnelle (Éq. 61).
• Fonctionnel d'énergie et analyse de la hessienne : Un fonctionnel d'énergie est construit dont l'équation d'Euler-Lagrange correspond à l'équation d'amplitude dérivée. La hessienne de ce fonctionnel est analysée pour identifier le point critique où la plus petite valeur propre s'annule, signalant le début de l'instabilité.
• Solutions analytiques : L'équation d'amplitude non linéaire est résolue dans deux limites :
  • Limite ductile : Où les pentes sont faibles, en négligeant les termes de gradient d'ordre supérieur.
  • Limite fragile : Où le cœur de la bande est dominé par des termes de gradient quartiques, nécessitant un développement asymptotique apparié (solutions intérieure et extérieure).

Contributions et résultats clés

1. Dérivation de l'équation d'amplitude non linéaire :
   Le résultat central est l'Éq. (61), une équation différentielle non linéaire pour le profil de la bande de cisaillement $f(\xi)$ :
   $$(\mu + \Sigma(m)) f''(\xi) + \frac{3}{2}(\lambda + 2\mu) (f'(\xi))^2 f''(\xi) = -A f(\xi) + B f^3(\xi)$$
   Ici, $\Sigma(m)$ représente la projection de la contrainte de fond, et $A$ et $B$ sont des paramètres dérivés des tenseurs de couplage et des paramètres d'écran.

2. Prédiction de la contrainte critique :
   En analysant la hessienne linéarisée du fonctionnel d'énergie, les auteurs dérivent un critère pour le début de l'instabilité. La condition critique se produit lorsque la plus petite valeur propre de la hessienne s'annule :
   $$(\mu + \Sigma(m)) \left(\frac{2\pi}{L}\right)^2 = A$$
   Cela fournit une prédiction pour la contrainte accumulée critique nécessaire pour déclencher le cisaillement localisé, la reliant directement au paramètre d'écran $A$ (lié à $\tilde{\kappa}_e^2$).

3. Profils analytiques pour les bandes ductiles et fragiles :

   • Solution ductile : Dans la limite des faibles pentes, la solution est un profil tanh :
     $$f(\xi) = f_0 \tanh\left(\frac{\xi - \xi_0}{\ell}\right)$$
     La largeur $\ell$ est contrôlée par le paramètre d'écran et les modules élastiques : $\ell \sim \sqrt{(\mu + \Sigma)/A}$.
   • Solution fragile : Dans la limite des fortes pentes, le cœur de la bande est gouverné par le terme quartique, conduisant à un cœur sinusoïdal apparié à une queue tanh. La largeur du cœur $L_c$ est déterminée par les coefficients du terme cubique de dipôle ($G$), distincte de la longueur d'écran. Cela explique la nature plus aiguë et plus localisée des bandes de cisaillement fragiles.

4. Échelles de longueur intrinsèques :
   La théorie démontre que le cisaillement localisé nécessite une échelle de longueur intrinsèque générée par les tenseurs de couplage dans le lagrangien. Cette échelle résulte de la compétition entre les termes de contrainte déstabilisants et les termes de gradient régularisants (écran). Contrairement aux modèles élastoplastiques où la largeur est arbitraire, ici la largeur est une propriété physique déterminée par les paramètres du matériau ($\mu, \lambda, A, B, G$).

5. Mode nul de la hessienne non linéaire :
   Les auteurs prouvent que la dérivée de la solution de la bande de cisaillement, $f'(\xi)$, est le mode nul (fonction propre avec valeur propre nulle) de la hessienne non linéaire. Cela relie l'invariance par translation de la solution de type « kink » directement à l'analyse de stabilité.

Portée et affirmations
L'article affirme offrir la première théorie analytique capable de prédire à la fois la contrainte critique pour le cisaillement localisé et le profil spatial détaillé de la bande, en distinguant les comportements ductiles et fragiles.

• Résolution des échecs précédents : Les auteurs soutiennent que les théories antérieures (par exemple, Rudnicki et Rice) ont échoué à fournir une solution car elles manquaient des non-linéarités nécessaires pour saturer l'instabilité. Le terme cubique ($B f^3$) dans l'équation d'amplitude dompte l'instabilité linéaire, tandis que le terme de gradient non linéaire contrôle la forme du cœur.
• Sélection physique de la localisation : La théorie postule que la localisation du cisaillement est un « objet continu véritablement écranté ». L'écran des interactions élastiques à longue portée (champs d'Eshelby) par les dipôles plastiques est essentiel ; sans cet écran, la déformation se propagerait globalement. La théorie génère naturellement une bande de largeur finie sans régularisation ad hoc.
• Distinction par rapport aux modèles élastoplastiques : L'article souligne que les modèles élastoplastiques standards sont incomplets pour prédire les profils de bandes de cisaillement car ils manquent d'échelles de longueur intrinsèques. En revanche, cette théorie intègre directement l'échelle de longueur dans le lagrangien via les paramètres de couplage, rendant la largeur de la bande une quantité prédictive dépendante du matériau.
• Modestie concernant les paramètres : Les auteurs reconnaissent que, bien que la forme fonctionnelle de la solution soit dérivée analytiquement, les valeurs numériques spécifiques pour des paramètres comme $B$ (qui détermine l'échelle de longueur interne dans la limite fragile) ne sont pas encore mesurées. Ils proposent que des simulations numériques futures sur des modèles classiques de formation de verres sont nécessaires pour déterminer ces paramètres et valider pleinement la théorie.

En résumé, l'article fournit un cadre variationnel rigoureux qui unifie les concepts d'écran plastique, d'élasticité non linéaire et de théorie de l'instabilité pour expliquer l'émergence, la criticité et la morphologie des bandes de cisaillement dans les solides amorphes.