💻 computer science

A new approach for the determination of through-thickness and free-edge stresses in composite laminates based on structural elements

Cet article propose une nouvelle méthode utilisant des éléments de cohésion structurels et une rigidité de pénalité précise dérivée des couches riches en résine pour récupérer les composantes de contraintes hors plan dans les éléments de coquilles minces basés sur la théorie de Kirchhoff-Love, permettant ainsi une prédiction précise de l'initiation des dommages dans les stratifiés composites.

Auteurs originaux : Xiaopeng Ai, Christos Kassapoglou, Boyang Chen

Publié 2026-07-31

📖 7 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Xiaopeng Ai, Christos Kassapoglou, Boyang Chen

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 que vous construisez un gratte-ciel avec des milliers de petites briques Lego colorées. Dans le monde réel, les ingénieurs n'utilisent pas de plastique ; ils utilisent des « composites stratifiés », qui sont comme des sandwichs ultra-résistants composés de couches de tissu imprégnées de colle (résine). Ces matériaux sont le ingrédient secret derrière les avions légers et les voitures rapides car ils sont incroyablement solides tout en étant très légers. Cependant, il existe un problème délicat : lorsque l'on plie ou étire ces sandwichs, les couches ont tendance à glisser les unes par rapport aux autres ou à se décoller, surtout près des bords. C'est ce qu'on appelle la « délamination », et c'est un tueur silencieux pour les structures.

Pour prédire si une aile d'avion tiendra le coup, les ingénieurs utilisent des simulations informatiques. Pour les structures minces, ils utilisent généralement un raccourci appelé « éléments de coque » (shell elements). Considérez un élément de coque comme une feuille de papier dans un programme informatique. Il est excellent pour calculer comment le papier s'étire ou se plie sur sa surface, mais il est incapable de comprendre ce qui se passe à l'intérieur de l'épaisseur du papier. C'est comme essayer de comprendre la garniture d'un sandwich en ne regardant que la tranche de pain supérieure ; vous manquez entièrement le fromage et la viande. Parce que ces modèles informatiques standards ne peuvent pas voir les contraintes internes (les forces qui poussent les couches à s'écarter), ils passent souvent à côté des signes avant-coureurs avant qu'une structure ne cède. Ce document intervient pour corriger cet angle mort, permettant aux ingénieurs de voir les forces cachées à l'intérieur du sandwich sans avoir à construire un modèle 3D massif, lent et complexe de chaque grain de colle.

La grande idée du papier : Voir les couches invisibles

Les auteurs, Xiaopeng Ai, Christos Kassapoglou et Boyang Chen, ont concocté une nouvelle façon de regarder à l'intérieur de ces sandwichs composites grâce à une astuce ingénieuse. Ils s'appuient sur une méthode qui traite chaque couche du stratifié comme un « élément structurel » (une version sophistiquée de notre feuille de papier) et les connecte avec des « éléments cohésifs ». Vous pouvez considérer ces éléments cohésifs comme des bandes de colle invisibles et extensibles entre les couches.

Par le passé, ces bandes de colle étaient un peu mystérieuses. Les ingénieurs devaient deviner leur rigidité. S'ils se trompaient, l'ordinateur pensait soit que les couches étaient collées avec une époxy ultra-résistante (trop rigide), soit avec un ruban adhésif faible (trop mou), menant à des réponses erronées sur la résistance de la structure. Le papier soutient que les anciennes méthodes de devinette de cette rigidité étaient erronées car elles traitaient la colle comme si elle n'avait aucune épaisseur, ce qui n'est pas vrai dans le monde réel où il existe toujours un minuscule espace riche en résine entre les couches.

Le secret de la « couche riche en résine »

La principale découverte de l'équipe est une nouvelle formule pour calculer exactement la rigidité de ces « bandes de colle ». Au lieu de deviner, ils ont réalisé que la rigidité dépend des propriétés réelles de la couche riche en résine — cette zone minuscule, chargée de colle, entre les couches de tissu.

Imaginez que vous étirez une pile de notes autocollantes. La force que vous ressentez dépend de l'épaisseur de la partie collante et de la force de la colle. Les auteurs ont dérivé une règle mathématique qui stipule : « La rigidité de la bande de colle est simplement la résistance de la colle divisée par l'épaisseur de cette minuscule couche de résine. »

Ils ont testé cette idée en effectuant des simulations sur différents types de sandwichs composites :

Laminés à plis croisés : Couches empilées comme un damier (0° et 90°).
Laminés à plis angulaires : Couches empilées selon des angles comme un zigzag (+45° et -45°).
Laminés quasi-isotropes : Empilements complexes qui agissent de la même manière dans toutes les directions.

Dans ces simulations, ils ont appliqué un étirement uniforme aux extrémités du stratifié, comme si l'on tirait sur un élastique. Ils ont ensuite comparé leur nouvelle méthode à deux choses :

Solutions analytiques : Ce sont des réponses mathématiques parfaites dérivées par d'autres scientifiques (comme Pagano et Mittelstedt) qui servent de « référence absolue ».
Modèles à éléments solides : Ce sont les modèles informatiques lourds et lents qui décomposent le sandwich en millions de minuscules blocs 3D pour voir chaque détail.

Ce qu'ils ont trouvé

Les résultats ont été impressionnants. Lorsque les auteurs ont utilisé leur nouvelle formule de rigidité « riche en résine », leurs modèles de coques simples ont prédit les contraintes internes de manière presque parfaite.

Précision : Leurs résultats correspondaient aux solutions mathématiques de « référence absolue » ainsi qu'aux modèles 3D de blocs lourds. Ils pouvaient prédire avec précision la contrainte normale interlaminaire (la force qui tente de décoller les couches, notée $\sigma_z$ ) et la contrainte de cisaillement (la force qui tente de faire glisser les couches les unes par rapport aux autres, notées $\tau_{xz}$ et $\tau_{yz}$ ).
Effets de bord libre : L'une des choses les plus difficiles à prédire est ce qui se passe au bord même du matériau, là où les couches sont exposées à l'air. C'est là que se produisent les « singularités de contrainte » (des pics de force extrêmes). La nouvelle méthode a capturé ces pics avec précision, alors que les anciennes méthodes les manquaient souvent ou nécessitaient des maillages incroyablement fins et lents pour les voir.
Robustesse : L'équipe a également vérifié si sa méthode était sensible à de petits changements. Ils ont demandé : « Et si la couche de résine est légèrement plus épaisse ou plus mince ? Et si la colle est légèrement plus forte ou plus faible ? » Ils ont constaté que leur méthode était très stable. Même si l'épaisseur de la couche de résine variait entre 0,6 × 10⁻³ pouce et 0,9 × 10⁻³ pouce, ou si la rigidité de la colle variait entre 3 GPa et 5 GPa, les résultats ne changeaient pas beaucoup. Cela signifie que les ingénieurs n'ont pas besoin de connaître les détails microscopiques exacts de la colle pour obtenir une bonne réponse.

Le gain de vitesse

La partie la plus excitante est peut-être la vitesse. Les modèles de blocs 3D lourds (les modèles à « éléments solides ») sont précis mais lents. Pour obtenir une bonne réponse, ils devaient diviser le stratifié en huit petits blocs pour chaque couche. La nouvelle méthode, cependant, n'avait besoin que d'un seul élément de coque par couche.

Dans leurs tests, la nouvelle méthode a réduit le temps de traitement informatique d'environ 45 %. Pour un cas de test spécifique, l'ancienne méthode a pris 9,52 secondes par étape de calcul, tandis que la nouvelle méthode n'a pris que 5,21 secondes. C'est une économie énorme lorsque l'on conçoit une aile d'avion entière composée de milliers de pièces.

Ce que cela signifie pour l'avenir

Le papier ne prétend pas avoir résolu tous les problèmes de l'univers. Ils admettent que pour des motifs de contraintes très complexes à l'intérieur d'une seule couche, leur méthode utilise une « interpolation linéaire » (une ligne droite) pour deviner les valeurs entre les couches, ce qui pourrait manquer de petits détails sinueux. Cependant, pour l'aspect global de la façon dont les couches se décollent et glissent, cela fonctionne magnifiquement.

En prouvant que l'on peut obtenir des données de contrainte 3D précises à l'aide d'éléments de coque simples et rapides, ce travail ouvre la voie aux ingénieurs pour simuler les dommages et la délamination de manière beaucoup plus efficace. C'est comme passer d'un télescope lent et flou à un télescope net et rapide, nous permettant de voir les dangers cachés dans nos structures légères avant même qu'elles ne prennent leur envol. Les auteurs suggèrent que ce cadre pourrait être la base d'outils futurs capables de simuler simultanément les dommages dans le plan (fissures dans le tissu) et la délamination (décollement des couches), rendant nos avions et nos voitures plus sûrs et plus légers.

Résumé technique : Une nouvelle approche pour les contraintes à travers l'épaisseur et de bord libre dans les stratifiés composites

Énoncé du problème
Les composites renforcés de fibres sont essentiels pour les structures aérospatiales légères, pourtant la prédiction de leur résistance et de leur tolérance aux dommages reste un défi. Bien que les éléments de coque minces basés sur l'hypothèse de Kirchhoff–Love soient largement utilisés pour leur efficacité de calcul, ils ne prennent intrinsèquement en compte que trois composantes de contrainte dans le plan ( $\sigma_x$ , $\sigma_y$ et $\tau_{xy}$ ). Par conséquent, les composantes de contrainte hors plan nécessaires aux critères de dommage tridimensionnels (3D) conventionnels sont indisponibles. Cette limitation empêche la prédiction précise de l'amorçage des dommages, particulièrement le délaminage interlaminaire et les effets de bord libre, qui sont pilotés par les contraintes normales et de cisaillement hors plan.

Les méthodes existantes pour récupérer ces contraintes reposent souvent sur des théories de coques de Mindlin–Reissner ou des techniques de post-traitement qui négligent soit les éléments cohésifs (limitant la simulation du délaminage), soit font face à des goulots d'étranglements de densité de maillage imposés par les modèles de zone cohésive. De plus, les formulations précédentes pour la rigidité de pénalité des éléments cohésifs supposent souvent une épaisseur nulle ou reposent sur des hypothèses cinématiques incompatibles avec l'hypothèse de Kirchhoff–Love, ce qui conduit à surestimer les niveaux de contrainte.

Méthodologie
Les auteurs proposent un cadre unifié basé sur des éléments structurels qui capture les effets de bord libre dans les stratifiés composites multicouches avec des empilements arbitraires. La méthodologie se compose de trois composantes fondamentales :

Formulation de l'élément structurel : Le stratifié est modélisé à l'aide d'éléments de coque cubiques pour chaque pli, connectés par des éléments cohésifs structurels. Cette approche conserve les hypothèses cinématiques de Kirchhoff–Love (où les normales restent droites et perpendiculaires au plan moyen) tout en éliminant la limitation de densité de maillage de la zone cohésive trouvée dans la modélisation de zone cohésive traditionnelle.
Mécanisme de récupération des contraintes : Les composantes de contrainte hors plan ( $\sigma_z$ , $\tau_{xz}$ , $\tau_{yz}$ ) sont récupérées directement à partir des tractions transmises par les éléments cohésifs structurels. En imposant l'équilibre à l'interface, les tractions cohésives sont équivalentes aux contraintes hors plan dans les plis adjacents. Ces contraintes interfaciales sont ensuite interpolées à travers l'épaisseur pour reconstruire le champ de contraintes 3D complet.
Rigidité de pénalité physiquement représentative : Une innovation critique est la dérivation d'une nouvelle formulation de la rigidité de pénalité pour les éléments cohésifs. Contrairement aux approches précédentes qui supposent une épaisseur nulle ou des déformations de cisaillement simples, ce travail définit la rigidité de pénalité comme une fonction des propriétés mécaniques et de l'épaisseur de la couche riche en résine entre les plis adjacents.
- La rigidité de pénalité normale ( $K_n$ ) est dérivée de $E_3 / h_{rr}$ , où $E_3$ est le module d'élasticité hors plan et $h_{rr}$ est l'épaisseur de la couche riche en résine.
- La rigidité de pénalité de cisaillement ( $K_s, K_t$ ) est dérivée de $G / h_{rr}$ , où $G$ est le module de cisaillement de la couche riche en résine.
- Cette formulation reconnaît que l'épaisseur effective de l'élément de coque cohésif est la moitié de la somme des épaisseurs des coques adjacentes, s'alignant ainsi sur les hypothèses cinématiques de la théorie de Kirchhoff–Love.

Résultats clés
La méthode proposée a été validée par rapport à des solutions analytiques (ex. Pagano, Mittelstedt) et des modèles d'éléments solides 3D haute fidélité (C3D8 dans Abaqus) à travers diverses configurations de stratifiés :

Stratifiés à plis croisés ([0°/90°]s) : La méthode a prédit avec précision les contraintes normales interlaminaires ( $\sigma_z$ ) et de cisaillement ( $\tau_{yz}$ ) à l'interface et au plan médian, montrant un excellent accord avec les solutions analytiques.
Stratifiés à plis angulaires ([+θ/−θ]s) : La méthode a capturé avec succès la contrainte de cisaillement interlaminaire ( $\tau_{xz}$ ) significative près des bords libres, correspondant aux résultats analytiques et des éléments solides.
Laminés quasi-isotropes et empilements généraux : L'approche a été étendue à des séquences d'empilement complexes (ex. [45/0/−45/90]s et empilements arbitraires) où aucune solution analytique n'existe. Dans ces cas, les résultats des éléments structurels ont démontré un accord étroit avec les modèles d'éléments solides, capturant avec précision les singularités de contrainte et les valeurs de pointe.
Analyse de sensibilité : Les résultats ont montré que les distributions de contraintes prédites sont relativement peu sensibles aux variations du module d'élasticité et du coefficient de Poisson de la couche riche en résine dans des plages réalistes. La méthode reste robuste même avec des propriétés approximatives de la couche riche en résine.
Efficacité de calcul : En utilisant un seul élément de coque par pli plutôt que plusieurs éléments solides à travers l'épaisseur, la méthode proposée a réduit le temps CPU par incrément d'environ 45 % tout en maintenant une précision comparable dans la prédiction des contraintes.

Signification et revendications
L'article affirme fournir une approche dédiée pour la récupération du champ de contraintes 3D complet au sein des éléments de coque de Kirchhoff–Love, une capacité auparavant indisponible pour les formulations de coques d'ordre élevé. En introduisant une formulation de rigidité de pénalité basée sur la couche riche en résine, les auteurs répondent à la divergence entre les hypothèses théoriques d'épaisseur nulle et l'épaisseur effective finie des éléments cohésifs structurels.

La signification de ce travail réside dans sa capacité à :

Permettre la prédiction précise des contraintes à travers l'épaisseur et de bord libre sans recourir à des modèles d'éléments solides 3D coûteux en calcul.
Faciliter la simulation combinée de dommages dans le plan et de délaminage sur des maillages grossiers en supprimant les limitations de densité de maillage des modèles de zone cohésive traditionnels.
Fournir un fondement physique pour les développements futurs des modèles d'éléments cohésifs structurels pour les composites.

Les auteurs concluent que cette méthode offre une alternative efficace et précise pour l'analyse de stratifiés composites complexes, particulièrement pour les applications nécessitant l'évaluation de l'amorçage du délaminage et des effets de bord libre.

La grande idée du papier : Voir les couches invisibles

Le secret de la « couche riche en résine »

Ce qu'ils ont trouvé

Le gain de vitesse

Ce que cela signifie pour l'avenir

Articles similaires