The fracture resistance of elastic networks increases with the density of defects like a random walk

Résumé technique : La résistance à la fracture des réseaux élastiques augmente avec la densité de défauts comme une marche aléatoire

Énoncé du problème
Les réseaux de ressorts désordonnés servent de modèle fondamental pour étudier la fracture dans divers matériaux, allant des céramiques et des réseaux polymères aux métamatériaux mécaniques. Une grandeur centrale en mécanique de la fracture est l'énergie de fracture apparente, $G_c$, qui mesure la résistance à la propagation d'une fissure préexistante dominante. Bien qu'il soit établi que le désordre peut augmenter $G_c$ par l'arrêt de fissure dû aux hétérogénéités locales, la dépendance quantitative de ce durcissement par rapport au degré de désordre (spécifiquement la densité de défauts) reste mal comprise. Des travaux antérieurs ont montré qu'en l'absence de fissure préexistante, les défauts favorisent généralement la nucléation et réduisent la résistance à la fracture. Cependant, dans les systèmes pré-fissurés, la relation est plus complexe : à faible densité de défauts, l'augmentation du désordre peut améliorer la résistance à la propagation de la fissure, tandis qu'à haute densité, le matériau s'affaiblit. Les mécanismes spécifiques reliant les statistiques microscopiques des défauts à l'énergie de fracture macroscopique dans le régime des fissures préexistantes nécessitent une clarification théorique et quantitative supplémentaire.

Méthodologie
Les auteurs emploient un réseau triangulaire bidimensionnel de ressorts linéaires sous chargement quasi-statique pour étudier l'effet de l'absence aléatoire de liaisons sur la propagation de la fissure.

• Système modèle : Le réseau est composé de ressorts de longueur au repos $\ell_r$ et de rigidité $k$. Un ressort se brise brusquement lorsque sa déformation dépasse un seuil $\varepsilon_{max} = 1$.
• Introduction de défauts : Une fraction $\nu$ de liaisons est retirée aléatoirement d'un réseau par ailleurs parfait.
• Protocole de simulation : Le système est chargé par des incréments de déplacement vertical uniformes. Le réseau est relâché vers l'équilibre statique en utilisant l'algorithme de minimisation FIRE. Les ressorts dépassant le seuil de déformation sont retirés séquentiellement, en commençant par les plus étirés, suivant un protocole spécifique pour simuler la propagation de la fissure.
• Mesures clés : L'énergie de fracture apparente $G_c(a)$ est définie comme le taux de libération d'énergie élastique critique requis pour propager la fissure d'une longueur initiale $L_c$ à $L_c + a$. L'étude se concentre sur l'énergie d'initiation $G_c(0)$ et l'énergie de rupture $G_c(a_f)$ à un avancement de fissure spécifique $a_f$.
• Cadre théorique : Les auteurs cartographient l'hétérogénéité élastique causée par les liaisons manquantes sur un paysage d'énergie de fracture locale équivalent, $\Gamma_{loc}(a)$. Ils appliquent des théories établies de l'arrêt de fissure planaire, en supposant que l'énergie de fracture apparente en tout point est déterminée par l'énergie de fracture locale maximale rencontrée le long du trajet de la fissure jusqu'à ce point.

Résultats clés

1. Durcissement par les défauts : Pour une concentration fixe de liaisons manquantes $\nu$, l'énergie de fracture apparente $G_c(a)$ augmente avec l'avancement de la fissure $a$. Contrairement aux réseaux parfaits qui échouent instantanément dès la rupture de la première liaison, les réseaux désordonnés présentent une propagation de fissure intermittente caractérisée par une alternance de croissance instable et d'arrêt.
2. Échelle avec la densité de défauts : L'écart type des fluctuations de l'énergie de fracture locale équivalente $\Gamma_{loc}$ augmente avec la racine carrée de la fraction de liaisons manquantes, $\sqrt{\nu}$. Par conséquent, l'énergie moyenne de fracture à la rupture $\langle G_c(a_f) \rangle$ suit également une échelle en $\sqrt{\nu}$.
3. Mécanisme de marche aléatoire : Cette échelle en $\sqrt{\nu}$ s'explique en modélisant les perturbations causées par les liaisons manquantes individuelles comme des variables aléatoires indépendantes. La perturbation totale de l'énergie de fracture locale est la superposition de ces effets individuels, analogue à une marche aléatoire. La variance de la perturbation totale croît linéairement avec le nombre de liaisons manquantes, conduisant à un écart type proportionnel à $\sqrt{\nu}$.
4. Queue exponentielle et croissance logarithmique : La fonction de densité de probabilité (PDF) de l'énergie de fracture locale $\Gamma_{loc}$ présente une queue exponentielle plutôt qu'une distribution purement gaussienne. Cette queue exponentielle conduit à une distribution de valeurs extrêmes de Gumbel pour l'énergie de fracture apparente. En conséquence, la médiane de $G_c(a)$ augmente logarithmiquement avec l'avancement de la fissure $a$ ($G_c \propto \ln a$).
5. Domination de l'arrêt de fissure : L'étude démontre que l'arrêt de fissure par les fluctuations locales de $\Gamma_{loc}$ est le mécanisme de durcissement dominant dans le régime étudié. Bien que les trajets de fissure puissent dévier légèrement des lignes droites (rugosité), le modèle d'arrêt de fissure planaire prédit avec précision le comportement observé, suggérant que la rugosité du trajet joue un rôle mineur par rapport à l'arrêt induit par l'hétérogénéité locale.

Signification et revendications
L'article revendique l'établissement d'un lien quantitatif entre le désordre microstructural et l'énergie de fracture macroscopique. En démontrant que les fluctuations de l'énergie de fracture locale suivent une échelle en $\sqrt{\nu}$ due à la superposition de type marche aléatoire des perturbations provenant des défauts individuels, les auteurs fournissent un cadre théorique pour prédire le durcissement des matériaux désordonnés. L'identification de la queue exponentielle dans le paysage d'énergie local explique l'augmentation logarithmique de l'énergie de fracture avec l'avancement de la fissure, un comportement distinct des systèmes sans fissures préexistantes.

Les auteurs notent que ces résultats devraient s'étendre à d'autres structures de réseau et systèmes continus où les liaisons manquantes correspondent à des microvides ou des microfissures, à condition que la densité de défauts reste faible et que la fracture se poursuive le long d'un trajet presque droit. Ils énoncent explicitement des limitations, notamment l'hypothèse d'une propagation de fissure planaire et l'indépendance des perturbations de liaison, reconnaissant que ces hypothèses peuvent être invalidées à des densités de défauts très élevées ou dans des systèmes tridimensionnels où se produit un moyennage du front de fissure. Ce travail vise à ouvrir la voie à des prédictions quantitatives de l'énergie de fracture dans une grande variété de matériaux désordonnés en reliant les statistiques microscopiques des défauts à la ténacité macroscopique.