A blended approach for evolving phase fields using peridynamics: Cyclic loading in quasi-brittle fracture

Cet article présente une théorie de champ de phase hybride couplée à la dynamique non locale pour modéliser la fracture quasi-fragile sous chargement cyclique, en intégrant la déformation plastique et l'adoucissement des modules élastiques tout en respectant les lois de la thermodynamique et en reproduisant avec succès les résultats expérimentaux sur le béton.

L'essentiel

🏗️ Le "Super-Simulateur" de la Rupture : Quand le Béton se Fend

Imaginez que vous tenez un morceau de béton ou de céramique. Si vous le pliez trop, il casse. Mais avant de casser, il commence à se fissurer de l'intérieur, il s'abîme petit à petit. Les ingénieurs veulent pouvoir prédire exactement où et quand cela va arriver, même si le matériau a déjà été plié plusieurs fois (comme un trombone qu'on plie et déplie jusqu'à ce qu'il casse).

Ce papier présente une nouvelle méthode mathématique, un peu comme un jeu vidéo ultra-réaliste pour les matériaux, capable de simuler ces cassures complexes.

1. L'Analogie du "Filet de Pêche" (La Méthode Péri-dynamique)

Traditionnellement, pour simuler la rupture, on divisait le matériau en petits blocs (comme des briques Lego) reliés entre eux. Si une brique cassait, le lien se brisait. Mais si la fissure passait entre les briques, ça posait problème.

Les auteurs de ce papier utilisent une approche différente, appelée Péri-dynamique.

• L'image : Imaginez que le matériau n'est pas fait de briques, mais d'un filet de pêche où chaque nœud est relié à ses voisins par des élastiques invisibles.
• Le fonctionnement : Chaque point du matériau "parle" directement à ses voisins immédiats. Si un élastique entre deux points s'étire trop, il se détend ou se rompt. Pas besoin de "briques" rigides. C'est plus fluide, comme de l'eau ou du tissu.

2. Le "Mémorisateur" (Le Champ de Phase)

Le vrai génie de ce papier, c'est qu'ils ont ajouté une mémoire à ce filet.

• Le problème : Si vous pliez une règle en plastique, elle ne revient pas exactement à sa forme initiale. Elle garde une trace de la déformation (c'est la déformation plastique). Si vous la pliez encore, elle s'affaiblit.
• La solution : Le modèle utilise ce qu'ils appellent un "champ de phase". C'est comme un thermomètre de santé pour chaque petit lien du filet.
  • Si le lien est neuf : le thermomètre est à 100% (vert).
  • Si le lien a été étiré et a gardé une trace : le thermomètre descend (jaune/orange).
  • Si le lien est cassé : le thermomètre est à 0% (rouge).
• L'innovation : Ce modèle a deux mémoires séparées : une pour quand on tire (traction) et une pour quand on pousse (compression). Le béton réagit différemment selon qu'on l'écrase ou qu'on l'étire.

3. Le "Cycle de Vie" (Chargement Cyclique)

La plupart des modèles anciens étaient bons pour dire "ça casse quand on tire fort une seule fois". Mais dans la vraie vie, les ponts et les bâtiments subissent des vibrations, des vents, des tremblements de terre (des charges qui vont et viennent).

• L'analogie : Imaginez un vieux jean. Si vous le tirez une fois, il résiste. Si vous le tirez, relâchez, tirez encore, il finit par s'user et se déchirer à un endroit précis.
• Ce nouveau modèle est capable de simuler cette fatigue. Il sait que si vous avez déjà étiré un lien, il sera plus faible la prochaine fois. Il calcule l'usure cumulée, comme un compteur de kilomètres sur une voiture.

4. La "Loi de la Conservation de l'Énergie" (Pas de Magie)

En physique, l'énergie ne disparaît pas. Quand un matériau casse, l'énergie utilisée pour le casser doit venir de quelque part (votre main qui pousse, le vent, etc.).

• Les auteurs ont prouvé mathématiquement que leur modèle respecte cette règle sacrée. Toute l'énergie dépensée pour créer une fissure est comptabilisée. C'est comme un compte en banque : on ne peut pas dépenser plus que ce qu'on a gagné. Cela rend les prédictions très fiables.

5. Les Résultats : Une Prédiction Parfaite

Pour vérifier leur théorie, ils ont comparé leur simulation à de vraies expériences sur du béton :

• Test 1 : Un bloc de béton qu'on plie jusqu'à la cassure.
• Test 2 : Un bloc qu'on plie, qu'on relâche, qu'on replie (chargement cyclique).
• Test 3 : Des blocs de différentes tailles.

Le résultat ? Leurs simulations ont donné des courbes de force et de fissuration qui correspondent presque parfaitement à la réalité physique. Même chose pour les effets de taille : un petit morceau de béton ne se comporte pas exactement comme un gros, et leur modèle a réussi à prédire cette différence sans avoir besoin de changer les règles du jeu.

En Résumé

Ce papier propose une nouvelle façon de voir la rupture des matériaux. Au lieu de les voir comme des blocs rigides qui se brisent, ils les voient comme un réseau de liens élastiques avec une mémoire.

C'est comme si on donnait à chaque atome du béton un petit carnet de notes où il écrit : "J'ai déjà été étiré à ce niveau, je suis un peu fatigué". Grâce à cela, les ingénieurs pourront mieux concevoir des bâtiments et des ponts qui résistent mieux aux séismes et aux usures du temps, en prédisant exactement où la fissure va apparaître avant qu'elle n'arrive.

C'est une avancée majeure qui mélange la physique des matériaux, les mathématiques pures et l'informatique pour mieux comprendre comment notre monde se brise.

Résumé technique

1. Problématique

La modélisation de la fissuration et de l'endommagement dans les matériaux quasi-fragiles (comme le béton, les maçonneries ou les géomatériaux) sous chargement cyclique reste un défi majeur en mécanique des milieux continus. Les approches traditionnelles de champ de phase couplent souvent une équation d'évolution pour le champ de phase (basée sur une équation de Ginzburg-Landau ou une contrainte variationnelle) avec les équations de mouvement de Newton. Ces méthodes peuvent souffrir de complexité numérique et de difficultés à garantir la cohérence thermodynamique, en particulier lors de cycles de chargement-déchargement impliquant une déformation plastique irréversible et un adoucissement élastique.

L'objectif de cet article est de formuler une théorie de champ unifiée capable de prédire l'endommagement, la plasticité et la propagation de fissures dans des matériaux quasi-fragiles sous chargement cyclique, en n'utilisant que la deuxième loi de Newton, sans équation d'évolution séparée pour le champ de phase.

2. Méthodologie : L'Approche Hybride (Blended Approach)

Les auteurs proposent une formulation basée sur la péridynamique (une théorie non locale) combinée à un concept de champ de phase à deux points dépendant de l'histoire.

• Loi de comportement non locale avec mémoire :

  • La déformation est définie par une décomposition additive d'une déformation élastique et d'une déformation irréversible (plastique).
  • Le modèle introduit deux champs de phase distincts, $\gamma^+$ (pour la traction) et $\gamma^-$ (pour la compression), qui évoluent en fonction de l'histoire de la déformation entre deux points $x$ et $y$.
  • Un paramètre clé, le rapport de déchargement ($\beta$), est introduit. Il représente le rapport entre la déformation plastique associée à l'adoucissement élastique et la déformation irréversible totale. Ce paramètre permet de moduler le comportement entre une rupture purement fragile ($\beta=0$) et un comportement quasi-fragile avec plasticité ($\beta > 0$).
  • La loi de contrainte-déformation est définie par une enveloppe de résistance et une famille de lois de déchargement élastique-endommagement-plastique.

• Équations d'évolution :

  • Le déplacement du matériau $u(t, x)$ est gouverné exclusivement par la deuxième loi de Newton ($\rho \ddot{u} + L_\epsilon[u] = b$), où $L_\epsilon$ est un opérateur non local représentant les forces internes.
  • Le champ de phase n'a pas d'équation d'évolution propre ; il est déterminé implicitement par le champ de déplacement via la loi constitutive dépendante de l'histoire.
  • La fissuration correspond à la rupture des "liens" (bonds) entre les points lorsque la force devient nulle en traction, tandis que les liens peuvent continuer à résister en compression (modélisant le frottement des faces de fissure).

• Existence et unicité :

  • Les auteurs établissent l'existence et l'unicité de la solution du problème aux limites initial (IBVP) dans un espace de Banach, en utilisant le théorème du point fixe de Banach, grâce à la continuité Lipschitz de la loi de contrainte par rapport au déplacement.

3. Contributions Clés

1. Unification Thermodynamique : Le modèle satisfait automatiquement le bilan d'énergie. Le taux de dissipation d'énergie est strictement positif, conforme aux lois de la thermodynamique, et découle directement des équations du mouvement sans hypothèses ad hoc.
2. Récupération de l'Énergie de Fracture : L'énergie de fissuration pour des fissures planes est récupérée directement du modèle via la théorie de la mesure géométrique. Elle correspond au produit du taux de libération d'énergie critique ($G_c$) et de la surface de la fissure, indépendamment de l'échelle de non-localité.
3. Modélisation de l'Hystérésis et de la Plasticité : Contrairement aux modèles élastiques-fragiles, cette approche capture explicitement l'hystérésis (boucles de chargement-déchargement) et l'accumulation de déformation plastique irréversible, essentiels pour les matériaux quasi-fragiles comme le béton.
4. Indépendance de l'Échelle de Non-localité ($\epsilon$) : L'échelle de non-localité $\epsilon$ n'est pas une propriété matérielle intrinsèque mais un paramètre numérique de régularisation. La calibration des propriétés matérielles (module d'Young, résistance, énergie de fracture) ne dépend pas de $\epsilon$, ce qui évite les problèmes de dépendance à la maille.
5. Méthode Sans Maillage (Mesh-free) : La discrétisation numérique est basée sur une intégration quadratique sans éléments ni connexions géométriques entre les nœuds, permettant une prédiction naturelle de la propagation de fissures complexes.

4. Résultats Numériques et Validation

Les auteurs ont validé le modèle sur plusieurs problèmes benchmarks en comparant les simulations avec des données expérimentales réelles :

• Convergence de maillage : Des études de convergence ($\epsilon$-convergence et $m$-convergence) montrent que le modèle converge vers la solution de la théorie des champs lorsque le rapport entre l'horizon et la taille de la grille augmente.
• Fissuration Mode I (Béton haute résistance) :
  • Simulation de tests de flexion trois points avec chargement monotone et cyclique.
  • Le modèle reproduit avec précision les courbes Force vs Ouverture de fissure (CMOD), y compris les boucles d'hystérésis et la dégradation de la rigidité, en accord avec les données de Chen et Liu (2023).
  • Visualisation de la localisation de la déformation plastique (zone de processus) devant l'extrémité de la fissure.
• Fissuration Mixte (Mode I + Mode II) :
  • Simulation de poutres en béton avec des entailles décalées (chargement mixte).
  • Le modèle prédit correctement les trajectoires de fissure et les courbes Force-CMOD pour différents rapports de décalage, en accord avec Jenq et Shah (1988).
• Échantillon en L (Sans entaille préexistante) :
  • Reproduction de la propagation de fissure complexe dans un panneau en forme de L (Winkler, 2001), montrant une transition du mode I au mode mixte.
• Effet de taille (Size Effect) :
  • Le modèle capture l'effet de taille déterministe dans les poutres en béton de différentes dimensions (García-Álvarez et al., 2012).
  • La contrainte nominale à la rupture diminue avec l'augmentation de la taille de la structure, un phénomène que les théories classiques de la plasticité ou de la contrainte maximale ne peuvent pas prédire sans règles ad hoc. Le modèle le reproduit naturellement grâce à la formulation non locale.

5. Signification et Impact

Ce travail établit la validité théorique et computationnelle d'une approche hybride péridynamique-champ de phase pour les problèmes physiques complexes impliquant endommagement et plasticité.

• Simplicité et Robustesse : En éliminant l'équation d'évolution séparée pour le champ de phase, le modèle réduit la complexité numérique tout en maintenant une précision élevée.
• Prédiction Réaliste : La capacité à modéliser l'hystérésis, la plasticité et l'effet de taille avec un seul jeu de constantes matérielles (module, résistance, énergie de fracture, rapport de déchargement) rend cette méthode très attractive pour la simulation de structures en béton et autres matériaux quasi-fragiles.
• Perspectives Futures : Le cadre théorique ouvre la voie à l'extension vers des milieux anisotropes, des matériaux composites et des formulations péridynamiques basées sur l'état complet (state-based), utilisant des champs de phase vectoriels ou tensoriels.

En résumé, cette approche offre une alternative puissante aux méthodes de champ de phase traditionnelles, garantissant la conservation de l'énergie et la cohérence thermodynamique tout en capturant fidèlement le comportement cyclique et la localisation de l'endommagement dans les matériaux quasi-fragiles.