🔬 condensed matter

Reversible viscoelasticity and irreversible elastoplasticity in the power law creep and yielding of gels and fibre network materials under stress

Cette étude computationnelle révèle que les réseaux de fibres athermiques présentent des régimes de fluage viscoélastique réversibles et d'élasto-plasticité irréversible distincts sous contrainte constante, où la connectivité marginale et la rupture des filaments pilotent des dynamiques de déformation suivant des lois de puissance qui déterminent finalement la capacité du matériau à la récupération de forme ou à la défaillance catastrophique.

Auteurs originaux : Michael J. Hertaeg, Suzanne M. Fielding

Publié 2026-01-23

📖 6 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Michael J. Hertaeg, Suzanne M. Fielding

Article original sous licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Imaginez un matériau comme un gel ou un morceau de feutre, non pas comme un bloc solide, mais comme une immense toile d'araignée désordonnée faite de minuscules fils extensibles. Ce document est une simulation informatique qui pose la question suivante : Que se passe-t-il lorsque l'on tire sur cette toile avec une force constante et que l'on maintient cette force ?

Les chercheurs voulaient comprendre deux choses principales :

Comment la toile s'étire-t-elle au fil du temps ?
Quand finit-elle par se rompre et casser ?

Voici le détail de leurs découvertes en utilisant des analogies simples.

La configuration : Une toile de fils

Les scientifiques ont construit un modèle numérique d'un monde en 2D rempli de ces fils.

Les fils : Ils agissent comme des élastiques. Si on les étire un peu, ils tirent en arrière. Si on les étire trop (au-delà d'un certain « point de rupture »), ils cassent de façon permanente.
Les connexions : Les fils sont attachés ensemble par des nœuds. Le nombre de fils attachés à chaque nœud est appelé la « connectivité » ( $Z$ $Z$ ).
- Toile serrée ( $Z > 4$ ) : Beaucoup de fils par nœud. Cette toile est rigide et forte.
- Toile lâche ( $Z < 4$ ) : Moins de fils par nœud. Cette toile est « molle » et vacillante.
- Toile « juste ce qu'il faut » ( $Z = 4$ ) : C'est le point « marginal » où la toile tient à peine debout.

L'expérience : La traction constante

Ils ont appliqué une traction constante et inchangée (contrainte de cisaillement) sur la toile et ont observé ce qui se passait au fil du temps. Ils ont mesuré l'étirement de la toile (déformation) et la vitesse à laquelle elle s'étirait (vitesse de déformation).

Découverte n°1 : La toile « juste ce qu'il faut » et l'étirement magique

Lorsque la toile était rigide (beaucoup de connexions) ou très lâche (peu de connexions), l'étirement se comportait de manière prévisible :

Toile rigide : Elle s'étirait un peu, puis s'arrêtait, se stabilisant dans une forme constante.
Tole lâche : Elle oscillait, s'étirait beaucoup, puis se stabilisait également.

Mais la toile « juste ce qu'il faut » ( $Z=4$ ) était spéciale.
Lorsqu'ils tiraient sur cette toile aux connexions marginales, elle ne se contentait pas de s'arrêter. Elle continuait de s'étirer selon un rythme mathématique très spécifique. La vitesse d'étirement ralentissait avec le temps, suivant une « loi de puissance ».

Analogie : Imaginez une personne marchant dans un couloir. Une personne normale peut marcher vite, puis ralentir, puis s'arrêter. Cette toile « juste ce qu'il faut » est comme une personne qui continue de marcher éternellement, mais dont les pas deviennent de plus en plus lents selon un motif parfait et prévisible. Les chercheurs ont découvert que cela se produit parce que la toile réorganise ses nœuds internes de manière complexe et non uniforme, et non parce que quelque chose se brise encore.

Découverte n°2 : La fuite lente et la rupture soudaine

Les chercheurs ont ensuite permis aux fils de se rompre s'ils étaient trop étirés. Cela a tout changé.

La fuite lente (Fluage élastoplastique) : Même si la toile était « lâche » (ce qui, normalement, ne présenterait pas cet étirement spécial en loi de puissance), le fait de permettre aux fils de casser a créé une longue période d'étirement lent et régulier. À mesure que la toile s'étirait, certains fils se cassaient. Cela endommageait la toile, la rendant plus faible, ce qui la faisait s'étirer davantage, provoquant la rupture de plus de fils.
- Analogie : Pensez à un trombone que l'on plie d'avant en arrière. Au début, il se plie facilement. Puis, après quelques pliages, une petite fissure se forme. Cette fissure rend le pliage suivant plus facile, ce qui agrandit la fissure. Finalement, l'ensemble casse. Dans la simulation, cette phase de « fissuration » pouvait durer très longtemps, créant ce même comportement d'étirement en « loi de puissance » observé pour la toile « juste ce qu'il faut », mais cette fois-ci causé par des dommages, et non par un simple réarrangement.
La rupture soudaine (Seuil d'écoulement) : Si la traction était assez forte, la toile atteignait finalement un « point de bascule ». Après une longue période d'étirement lent et d'accumulation de dommages, la toile subissait une défaillance catastrophique soudaine.
- Le délai : Plus la traction était forte, plus la rupture était rapide. Mais si la traction était juste légèrement inférieure au point de rupture, la toile pouvait tenir pendant un temps étonnamment long avant de finir par céder. Les chercheurs ont cartographié précisément comment ce « temps d'attente » variait en fonction de la force de la traction et de la résistance des fils individuels.

Découverte n°3 : La toile peut-elle guérir ? (Récupération)

Enfin, ils ont demandé : Si nous arrêtons de tirer, la toile reprend-elle sa forme d'origine ?

La partie élastique : Si la toile était rigide et que les fils ne s'étaient pas encore cassés, elle se comportait comme un élastique. Quand on lâchait la pression, la toile reprenait parfaitement sa forme d'origine.
La partie plastique : Si la toile était lâche, ou si certains fils avaient déjà cassé, elle se comportait comme un morceau d'argile. Quand on lâchait la pression, la toile restait étirée. Elle ne pouvait pas récupérer totalement car les dommages (fils cassés) étaient permanents.
La surprise : Même dans une toile rigide, si quelques fils se cassaient, la toile ne pouvait pas récupérer totalement. Cependant, si la toile était si bien connectée que même après la rupture de certains fils, elle conservait suffisamment de connexions pour rester « rigide », elle pouvait encore rebondir presque parfaitement. C'est comme un pont : si vous perdez quelques boulons, il peut encore garder sa forme ; si vous en perdez trop, il s'affaisse de façon permanente.

Résumé

Ce document explique pourquoi certains matériaux mous (comme les gels ou les tissus biologiques) s'étirent lentement et de manière prévisible pendant longtemps avant de se rompre.

Mécanisme A : Si le matériau est parfaitement équilibré (juste assez de connexions), il s'étire lentement en raison d'un réarrangement interne.
Mécanisme B : Si le matériau est plus faible ou si la traction est forte, il s'étire lentement parce qu'il se brise lentement lui-même.

Les deux mécanismes mènent au même comportement d'étirement en « loi de puissance », mais l'un est réversible (élastique) et l'autre est un dommage permanent (plastique). L'étude aide à expliquer pourquoi ces matériaux tiennent parfois bon pendant longtemps sous une charge constante avant de subir une défaillance soudaine.

Résumé technique : Viscoélasticité réversible et élastoplasticité irréversible dans le fluage en loi de puissance et l'écoulement des gels et des matériaux à réseaux de fibres

Énoncé du problème
Les matériaux à réseaux, allant des cytosquelettes biologiques et des hydrogels aux papiers et gels colloïdaux, présentent fréquemment un fluage en loi de puissance sous contrainte constante, où la vitesse de déformation décroît selon $\dot{\gamma} \sim t^n$ ( $n < 0$ ) et la déformation croît de manière sous-linéaire. Bien que ce comportement soit ubiquitaire dans la matière molle, les mécanismes physiques sous-jacents font l'objet de débats. Plus précisément, il n'est pas clair si le fluage en loi de puissance provient de déformations viscoélastiques purement réversibles (réarrangements non affines) ou d'une dégradation élastoplastique irréversible (rupture de liaisons). De plus, les conditions sous lesquelles ces matériaux passent d'un fluage lent à un écoulement (yielding) et une rupture catastrophiques, ainsi que la mesure dans laquelle ils recouvrent leur forme originale lors du déchargement, nécessitent une caractérisation systématique. Ce travail aborde ces questions en utilisant un modèle computationnel minimal de réseaux de fibres athermiques.

Méthodologie
Les auteurs emploient des simulations numériques d'un modèle de réseau désordonné athermique minimal en $d=2$ dimensions spatiales.

Construction du modèle : Les réseaux sont générés par un empilement aléatoire de disques bidispersés, reliant les centres chevauchants par des ressorts de type Hookean, et par élagage itératif pour atteindre un nombre de coordination moyen $Z$ cible. Chaque ressort est élastique linéaire jusqu'à une déformation de rupture $b$ , au-delà de laquelle il se rompt de manière irréversible.
Protocole rhéologique : L'étude utilise un protocole de contrainte de cisaillement imposée constante $\Sigma$ , moins courant dans les simulations que les protocoles à vitesse de déformation constante.
Algorithme numérique : Pour maintenir une contrainte constante, les auteurs ont développé et validé deux algorithmes. L'algorithme I ajuste itérativement la vitesse de déformation via une méthode de Newton-Raphson pour correspondre à la contrainte cible. L'algorithme II impose une vitesse de déformation $\dot{\gamma} = (\Sigma - \sigma)/\eta$ en utilisant un solvant Newtonien fictif de viscosité $\eta$ . Les auteurs démontrent que les deux algorithmes convergent vers les mêmes résultats dans la limite $\eta \to 0$ , utilisant par la suite l'algorithme II, plus efficace, avec $\eta = 0,16$ .
Paramètres : Les simulations sont réalisées avec $N=16\,000$ nœuds. Les variables clés incluent la connectivité du réseau $Z$ , le seuil de rupture des filaments $b$ , et la contrainte appliquée $\Sigma$ .

Contributions clés et résultats

Régime purement viscoélastique (sans rupture de filament) :
- Réseaux rigides ( $Z > 4$ ) : Pour les réseaux situés au-dessus du seuil de rigidité de Maxwell ( $Z > 2d = 4$ ), le matériau se comporte comme un solide de Hooke à long terme, atteignant une déformation stationnaire $\gamma_{ss} \propto \Sigma$ . Le fluage transitoire est purement viscoélastique, piloté par des déformations internes non affines.
- Réseaux mous/« floppy » ( $Z < 4$ ) : Les réseaux situés sous le seuil de rigidité présentent une réponse « molle » où une contrainte infime induit une déformation finie. La vitesse de déformation décroît selon $\dot{\gamma} \sim t^{-2}$ , et la déformation sature vers une constante.
- Connectivité marginale ( $Z \approx 4$ ) : Un régime unique émerge pour les réseaux à connectivité marginale ( $Z \gtrsim 4$ ). Ici, le système présente un fluage en loi de puissance soutenu avec $\dot{\gamma} \sim t^{-0,5}$ et $\gamma \sim t^{0,5}$ sur une durée prolongée avant de saturer. Ce comportement est purement viscoélastique et provient de déformations non affines sans aucune rupture de liaison.
Régime élastoplastique (avec rupture de filament) :
- Fluage en loi de puissance dans les gels sous-rigides : Lorsque la rupture de filament est autorisée ( $b < \infty$ ), la plage de conditions supportant le fluage en loi de puissance s'élargit. Spécifiquement, les réseaux initialement mous ( $Z < 4$ ) qui présenteraient autrement une saturation rapide dans le cas purement viscoélastique présentent désormais un régime prolongé de fluage en loi de puissance ( $\dot{\gamma} \sim t^n$ avec $-1 < n < 0$ ) piloté par l'accumulation de dommages plastiques.
- Écoulement (Yielding) et rupture : Pour des contraintes dépassant une contrainte d'écoulement critique $\Sigma_y$ , le réseau subit une rupture catastrophique après un temps de retard $\tau_y$ . La contrainte d'écoulement varie approximativement de manière linéaire avec le seuil de rupture $b$ . Le temps de retard $\tau_y$ augmente de façon spectaculaire lorsque $\Sigma$ approche $\Sigma_y$ par le haut.
- Validation : La contrainte d'écoulement mesurée dans les simulations imposées par la contrainte concorde avec la limite de vitesse de déformation nulle dans les simulations de démarrage imposées par la vitesse de déformation, confirmant la robustesse du point d'écoulement mesuré.
Récupération de déformation :
- L'étude quantifie la déformation recouvrable versus non recouvrable lors du déchargement.
- Réseaux rigides ( $Z > 4$ ) : Une récupération complète de la forme se produit si la contrainte est retirée avant qu'un dommage plastique significatif ne s'accumule. Remarquablement, même si une petite fraction de ressorts se brise, le réseau peut recouvrer entièrement sa forme, à condition que la connectivité restante reste au-dessus du seuil de rigidité ( $Z > 4$ ).
- Réseaux mous ( $Z < 4$ ) : Ces réseaux présentent une récupération partielle ou nulle, même à de faibles déformations, car ils manquent de la structure élastique nécessaire pour le rappel. Une déformation non recouvrable s'installe immédiatement dès le chargement, même avec une rupture de liaison négligeable.
Localisation des dommages :
- L'analyse visuelle des modes de rupture révèle des motifs distincts. Dans les réseaux initialement rigides, les dommages plastiques se localisent en bandes macroscopiques traversant le système. En revanche, les dommages dans les réseaux initialement mous sont distribués de manière plus diffusive à travers le matériau.

Signification et affirmations
L'article affirme identifier deux mécanismes physiques distincts capables de générer un fluage en loi de puissance dans les matériaux à réseaux :

Mécanisme viscoélastique : Émanant des réseaux à connectivité marginale ( $Z \gtrsim 4$ ) dû aux déformations non affines sans dommage.
Mécanisme élastoplastique : Émanant des réseaux initialement mous ( $Z < 4$ ) dû à l'accumulation progressive de ruptures de liaisons.

Les auteurs concluent que, bien que ces deux mécanismes expliquent le fluage en loi de puissance dans leur modèle minimal, ils ne prétendent pas que ce sont les seuls mécanismes responsables du phénomène observé dans la grande variété de matériaux mous expérimentaux. Ils suggèrent modestement que d'autres facteurs, tels que le mouvement des chaînes de type Rouse, les structures de réseaux fractales ou la dynamique thermique des lieurs, peuvent également contribuer à l'ubiquité du fluage en loi de puissance dans les systèmes réels. Ce travail fournit un cadre computationnel pour distinguer la viscoélasticité réversible de l'élastoplasticité irréversible dans le contexte du fluage et de l'écoulement.