A multi-phase-field model for fiber-reinforced composite laminates based on puck failure theory

Cet article propose un modèle multi-champs de phase bidimensionnel basé sur la théorie de la rupture de Puck et utilisant une méthode de superposition de maillages pour prédire avec précision les modes de défaillance distincts (fibres et interfibres) et leurs interactions dans les laminés composites renforcés de fibres.

L'essentiel

Voici une explication simple et imagée de ce papier scientifique, conçue pour être comprise par tous, sans jargon technique compliqué.

🛩️ Le Défi : Pourquoi les avions ne tombent pas (et comment on les teste)

Imaginez que vous construisez un avion avec des couches de tissu très résistant, comme un mille-feuille géant fait de fibres de carbone et de résine. C'est ce qu'on appelle un composite. C'est léger et solide, parfait pour les avions comme le Boeing 787.

Mais il y a un problème : si vous tirez ou poussez trop fort sur ce mille-feuille, il se casse. Et ce n'est pas une simple cassure comme du verre. C'est compliqué :

1. Les fibres (les "spaghetti" solides) peuvent se rompre.
2. La colle entre les fibres (la "matrice") peut se fissurer.
3. Les couches peuvent se décoller les unes des autres (comme les feuillets d'un livre).

Les ingénieurs doivent prédire exactement où et comment ça va casser avant de construire l'avion. Jusqu'à présent, ils devaient faire des centaines de tests physiques (casser des vrais morceaux d'avion), ce qui coûte cher et prend du temps.

🧠 La Solution : Un "Super-Détective" Numérique

Les auteurs de ce papier ont créé un nouveau logiciel (un modèle mathématique) qui agit comme un super-détective. Au lieu de casser des vrais matériaux, il simule la casse sur un ordinateur avec une précision incroyable.

Voici comment ils ont fait, avec des analogies simples :

1. La Théorie "Puck" : Le Code Secret de la Cassure

Imaginez que chaque couche de votre matériau a un "code secret" qui dit : "Si je suis tiré ici, je casse. Si je suis poussé là, je plie."
Les chercheurs ont utilisé une règle célèbre appelée la théorie de Puck. C'est comme un manuel d'instructions qui dit exactement quand les fibres cassent et quand la colle entre elles se fissure.

• L'analogie : C'est comme si vous aviez deux alarmes différentes. L'une sonne si les "spaghetti" (fibres) cassent, l'autre si la "colle" (matrice) se fissure. Ce modèle utilise les deux alarmes en même temps.

2. La Méthode "Calque" (Mesh Overlay) : Le Mille-Feuille Virtuel

Le plus dur avec les composites, c'est qu'ils sont faits de plusieurs couches avec des orientations différentes (une couche horizontale, une verticale, une en diagonale...).

• L'ancienne méthode : C'était comme essayer de dessiner un mille-feuille en 3D, ce qui est très lourd pour l'ordinateur.
• La nouvelle méthode (Calque) : Imaginez que vous prenez une feuille de papier transparente (une couche). Vous la posez sur votre écran. Ensuite, vous prenez une deuxième feuille transparente, vous la tournez de 45 degrés, et vous la posez par-dessus la première, exactement au même endroit. Vous faites ça pour toutes les couches.
  • Le génie : Toutes les couches partagent les mêmes points (les nœuds), mais chacune garde sa propre "personnalité" (sa direction). Cela permet de simuler un matériau épais et complexe sans alourdir l'ordinateur. C'est comme superposer des calques de dessin : on voit tout, mais c'est léger à manipuler.

3. Le Champ de Phase : La "Tache d'Encre" qui se Propage

Comment le logiciel voit-il la fissure ? Il n'utilise pas une ligne noire nette. Il utilise une sorte de tache d'encre floue.

• L'analogie : Imaginez une goutte d'encre qui commence à se répandre sur du papier.
  • Là où l'encre est blanche (0), le matériau est intact.
  • Là où l'encre est noire (1), le matériau est cassé.
  • Entre les deux, c'est la zone de transition où le matériau est en train de se fissurer.
    Ce système permet de voir la fissure grandir, bifurquer et se propager de manière très naturelle, comme une vraie fissure dans la glace.

🧪 Les Tests : A-t-il raison ?

Pour vérifier si leur "Super-Détective" fonctionne, ils l'ont mis à l'épreuve avec quatre défis classiques, comme des examens de conduite :

1. Le Test de Traction/Compression (Le "Coupon") : Ils ont tiré et poussé sur des petits morceaux de composite. Le logiciel a prédit la force de rupture presque exactement comme dans la réalité.
2. Le Trou dans le Tissu (Open Hole) : Imaginez un tissu avec un trou au milieu. Où va se fissurer le tissu autour du trou ? Le logiciel a prédit que la fissure partait du trou et suivait les fibres, exactement comme sur les photos de vrais tests.
3. Le Test de Tension Compacte : Un test pour voir comment le matériau résiste à la déchirure. Là encore, le logiciel a suivi la courbe de force et de déformation de très près.
4. Le Double Entaille (Double Notched) : Deux trous dans le matériau. C'est très complexe car les fissures peuvent aller dans toutes les directions. Le logiciel a réussi à montrer comment les fissures de la "colle" et des "fibres" interagissaient entre elles.

🏆 Le Résultat Final

Ce papier nous dit essentiellement : "Nous avons créé un outil numérique rapide et précis qui peut prédire comment les matériaux composites cassent, couche par couche, sans avoir besoin de casser des milliers d'échantillons réels."

• Pourquoi c'est important ? Cela permet aux ingénieurs de concevoir des avions, des voitures de course ou des éoliennes plus sûrs et plus légers, en testant des milliers de scénarios virtuels avant de fabriquer le moindre objet physique.
• La limite : Pour l'instant, le modèle ne simule pas le décollage des couches (délaminage) entre elles, car il travaille en 2D (comme un dessin), mais il est déjà très performant pour tout ce qui se passe dans les couches.

En résumé, c'est comme avoir une machine à voyager dans le temps pour voir la casse avant qu'elle n'arrive, en utilisant des règles mathématiques intelligentes et une astuce de superposition de calques.

Résumé technique

Voici un résumé technique détaillé de l'article « A Multi-Phase-Field Model for Fiber-Reinforced Composite Laminates Based on Puck Failure Theory », rédigé en français.

1. Problématique et Contexte

Les composites à renfort de fibres (FRC) sont largement utilisés dans l'aérospatiale et l'industrie automobile en raison de leur rapport poids/résistance exceptionnel. Cependant, la prédiction de leur rupture reste un défi majeur en raison de la complexité de leurs mécanismes de défaillance, qui incluent :

• Le fissurage matriciel intralaminar (transversal et longitudinal).
• La rupture des fibres.
• Le délaminage (fissuration interlaminar).
• L'interaction complexe entre ces modes de rupture selon la séquence d'empilement des plis.

Les méthodes numériques existantes (éléments cohésifs, modèles de fissures lissées, XFEM, peridynamique) présentent souvent des limitations, notamment pour capturer la nucléation des fissures et la propagation dans des séquences d'empilement complexes sans recourir à des maillages 3D extrêmement coûteux. De plus, les modèles de champ de phase (Phase-Field - PF) classiques peinent souvent à gérer les charges de compression et à distinguer correctement les modes de rupture des fibres et de la matrice dans des laminés multicouches.

2. Méthodologie Proposée

Les auteurs proposent un modèle multi-champ de phase (Multi-Phase-Field) bidimensionnel couplé à une méthode de superposition de maillages (mesh overlay) et basé sur la théorie de la rupture de Puck.

A. Formulation Variationnelle et Champs de Phase

Le modèle utilise deux variables de champ de phase indépendantes :

• $d_f$ : Pour la rupture dominée par les fibres.
• $d_m$ : Pour la rupture dominée par l'interface fibre/matrice (inter-fibre).

L'énergie totale du système est décomposée en énergie élastique stockée (dégradée par les fonctions d'endommagement) et énergie de surface (liée à la création de nouvelles surfaces de fissure). Deux échelles de longueur caractéristiques distinctes ($\ell_f$ et $\ell_m$) et deux tenseurs structuraux sont utilisés pour décrire précisément les modes de dommage anisotropes.

B. Intégration de la Théorie de Puck

Contrairement aux modèles PF classiques qui reposent souvent sur une séparation énergétique simple (tension/compression), ce modèle intègre les critères de rupture de Puck pour déclencher l'initiation de l'endommagement.

• La théorie de Puck distingue les modes de rupture : rupture de fibre en traction/compression, et trois modes de rupture inter-fibre (Mode A, B, C).
• Les facteurs d'exposition ($F$) calculés à partir des contraintes locales déterminent si l'énergie motrice est suffisante pour activer les variables de champ de phase.
• Cela permet de capturer correctement la nucléation des fissures, un point faible des modèles variationnels classiques.

C. Méthode de Superposition de Maillages (Mesh Overlay)

Pour simuler des laminés avec différentes orientations de plis sans utiliser un maillage 3D complexe :

• Chaque pli est représenté par un maillage séparé mais superposé.
• Les nœuds sont partagés entre tous les plis (couplage des déplacements), assurant la continuité cinématique.
• Les contraintes et les propriétés matérielles sont découplées et dépendent de l'orientation des fibres ($\theta_i$) de chaque pli.
• Limitation : Cette approche 2D suppose des interfaces parfaites et ne modélise pas explicitement le délaminage (rupture interlaminar), se concentrant sur la propagation des fissures dans le plan.

D. Implémentation Numérique

Le modèle est implémenté via un élément utilisateur (UEL) dans le code Abaqus, résolvant le problème variationnel via une méthode de Newton-Raphson.

3. Contributions Clés

1. Modélisation Multi-échelles et Multi-modes : Introduction de deux variables de champ de phase distinctes et de deux échelles de longueur pour séparer physiquement la rupture des fibres de celle de la matrice.
2. Couplage avec Puck : Utilisation de la théorie de Puck pour définir les forces motrices d'endommagement, permettant une prédiction réaliste de la nucléation des fissures sous compression et traction, là où les modèles PF classiques échouent.
3. Efficacité Computationnelle : La méthode de superposition de maillages permet de simuler des laminés complexes (multi-plis, orientations variées) en 2D, évitant le coût prohibitif des modèles 3D tout en capturant les interactions de dommages dans le plan.
4. Validation Étendue : Le modèle est validé sur une large gamme de géométries et de séquences d'empilement, comparant des résultats qualitatifs (morphologie des fissures) et quantitatifs (courbes force-déplacement, charges maximales).

4. Résultats et Validation

Quatre cas de référence ont été simulés et comparés à des données expérimentales issues de la littérature :

• Essais de coupons (Traction/Compression) :
  • Simulation de laminés à plis croisés ($\pm 45^\circ, \pm 60^\circ, \pm 75^\circ$).
  • Le modèle prédit correctement les modes de rupture (fissuration matricielle suivie de rupture de fibre) et les charges maximales, bien qu'il tende à légèrement surestimer la résistance en l'absence de plasticité.
• Traction avec trou (Open Hole Tension - OHT) :
  • Cas de laminés $[0/90]_s$ et $[45/-45]_s$.
  • Pour $[0/90]_s$, le modèle reproduit la propagation de fissures inter-fibres parallèles à la charge, suivie de la rupture des fibres perpendiculairement.
  • Pour $[45/-45]_s$, le modèle capture l'asymétrie de la propagation des fissures (effet d'ombrage ou crack shielding) et la charge de rupture avec une grande précision.
• Traction Compacte (Compact Tension - CT) :
  • Comparaison sur des laminés à plis groupés (blocked ply).
  • Le modèle reproduit bien le comportement linéaire initial et la propagation des fissures, bien que l'absence de délaminage dans le modèle entraîne une surestimation de la charge finale par rapport aux essais où le délaminage dissipe de l'énergie.
• Traction Double Entaillée (Double-Edged Notched Tension) :
  • Étude de laminés croisés, quasi-isotropes et isotropes.
  • Le modèle capture efficacement l'interaction compétitive entre la rupture des fibres et de la matrice, ainsi que la propagation de fissures subcritiques.

5. Signification et Conclusion

Cet article présente une avancée significative dans la modélisation numérique des composites stratifiés. En combinant la robustesse de la théorie de Puck pour l'initiation de la rupture avec la flexibilité des méthodes de champ de phase pour la propagation, l'approche proposée surmonte les limitations des modèles existants concernant la prédiction de la rupture en compression et la complexité des séquences d'empilement.

Bien que le modèle soit actuellement limité aux phénomènes 2D (excluant le délaminage explicite), il offre un outil computationnel efficace et précis pour l'évaluation qualitative et quantitative de la résistance des laminés. Les auteurs mettent à disposition les codes et les données, favorisant la transparence et la reproductibilité pour la communauté scientifique. Ce travail constitue une base solide pour le développement futur de modèles de dommages plus complets intégrant potentiellement le délaminage et le comportement plastique.