Prédiction du vieillissement hygrothermique et de la dégradation mécanique des matériaux composites

Résumé

Les matériaux composites à matrice polymère sont largement utilisés dans les secteurs aéronautique, automobile, énergétique et maritime grâce à leur excellent rapport résistance mécanique/masse. Cependant, leurs performances peuvent être fortement affectées par les conditions environnementales sévères, notamment l’exposition simultanée à l’humidité et aux variations de température. Le vieillissement hygrothermique entraîne des phénomènes complexes tels que la diffusion d’eau dans la matrice polymère, la plastification du matériau, la dégradation de l’interface fibre/matrice et l’apparition de contraintes internes.

Cette étude présente une approche numérique visant à analyser l’influence des effets thermiques, hygroscopiques et hygrothermiques sur le comportement mécanique des composites à matrice polymère renforcée par fibres. Un modèle multiphysique couplant diffusion d’humidité, transfert thermique et comportement mécanique est développé par une approche éléments finis. L’évolution des propriétés mécaniques, notamment le module d’élasticité, la résistance à la traction et les contraintes internes, est étudiée sous différentes conditions environnementales.

Les résultats montrent que l’absorption d’humidité et les cycles thermiques provoquent une modification significative du comportement mécanique du composite. Le couplage température-humidité accélère les mécanismes de dégradation en augmentant la mobilité moléculaire de la matrice et en générant des contraintes différentielles entre les constituants du matériau. Cette étude fournit une meilleure compréhension des mécanismes de vieillissement et permet d’améliorer la prédiction de durée de vie des structures composites soumises à des environnements variables.

Mots-clés : Matériaux composites, vieillissement hygrothermique, diffusion d’humidité, éléments finis, matrice polymère, dégradation mécanique

Introduction

Les matériaux composites à matrice polymère représentent une alternative stratégique aux matériaux métalliques traditionnels grâce à leurs propriétés spécifiques élevées, leur faible masse volumique et leur possibilité d’adapter leurs caractéristiques mécaniques par l’orientation des renforts.

Les composites fibres/polymère sont aujourd’hui employés dans des applications nécessitant une grande fiabilité structurelle, notamment les pales d’éoliennes, les pièces aéronautiques, les composants automobiles et les structures marines.

Cependant, contrairement aux matériaux métalliques, les composites polymères présentent une sensibilité importante aux conditions environnementales. L’exposition prolongée à l’humidité et aux variations thermiques peut provoquer une diminution progressive de leurs performances.

L’eau absorbée par la matrice polymère diffuse progressivement dans le matériau selon un mécanisme généralement décrit par la loi de Fick. Cette absorption peut provoquer:

• Une augmentation du volume de la matrice ;
• Une diminution de la température de transition vitreuse ;
• Une réduction de la rigidité ;
• Une dégradation de l’interface fibre/matrice.

Lorsque l’humidité est combinée avec des variations de température, des phénomènes supplémentaires apparaissent. Les différences entre les coefficients de dilatation thermique des fibres et de la matrice génèrent des contraintes internes pouvant favoriser l’apparition de microfissures.

La compréhension et la prédiction du vieillissement hygrothermique constituent donc un enjeu majeur pour garantir la durée de vie des structures composites.

L’objectif de cette étude est de développer un modèle numérique permettant d’évaluer l’évolution du comportement mécanique d’un matériau composite soumis à des conditions thermiques, hygroscopiques et hygrothermiques.

Matériaux et modèle étudié

Le matériau étudié est un composite constitué d’une matrice polymère renforcée par fibres continues.

La structure composite est définie par :

• Une matrice époxy ;
• Un renforcement fibreux unidirectionnel ;
• Une interface fibre/matrice considérée comme zone critique.

Les propriétés initiales considérées sont :

• Module d’élasticité longitudinal ;
• Module transversal ;
• Coefficient de dilatation thermique ;
• Coefficient de diffusion d’humidité ;
• Résistance mécanique initiale.

Le comportement du matériau est étudié avant et après exposition environnementale afin d’évaluer les mécanismes de dégradation.

Modélisation de la diffusion d’humidité

L’absorption d’eau dans le composite est modélisée par l’équation de diffusion de Fick :

∂C/∂t = D∇²C

où :

• C représente la concentration en humidité ;
• D représente le coefficient de diffusion ;
• t représente le temps.

Le coefficient de diffusion dépend de la température selon une relation d’Arrhenius :

D(T)=D₀ exp(-Ea/RT)

où :

• D₀ est le facteur pré-exponentiel ;
• Ea est l’énergie d’activation ;
• R est la constante des gaz parfaits ;
• T est la température absolue.

Cette relation permet de représenter l’accélération de la diffusion lorsque la température augmente.

Couplage thermo-hygro-mécanique

Afin de reproduire les conditions réelles d’utilisation, un modèle couplé est développé.

Trois phénomènes sont considérés :

Effet thermique

Les variations de température provoquent des déformations thermiques :

εth = αΔT

où α représente le coefficient de dilatation thermique.

Effet hygroscopique

L’absorption d’eau entraîne une déformation hygroscopique :

εh = βC

où β représente le coefficient de gonflement hygroscopique.

Interaction thermo-hygromécanique

Les contraintes totales sont obtenues par :

σ = C(ε – εth – εh)

Cette formulation permet d’évaluer les contraintes internes générées par les interactions entre humidité, température et chargement mécanique.

Méthodologie numérique

L’étude numérique a été réalisée à l’aide d’une approche par éléments finis afin d’analyser le comportement hygrothermique et mécanique d’un matériau composite à matrice polymère soumis à des conditions environnementales sévères. La modélisation a été effectuée principalement sous ANSYS Workbench, en utilisant le module ACP Pre pour la définition du stratifié composite, puis une analyse modale et de stabilité a été conduite afin d’évaluer la réponse dynamique et le risque de flambement thermique.

La démarche numérique suivie peut être résumée comme suit :

Définition de la géométrie du modèle composite

Une plaque composite stratifiée a été modélisée en tenant compte de ses dimensions géométriques, de son épaisseur et du rapport géométrique (a/h), paramètre important dans l’étude de la stabilité thermique.

Définition de l’architecture du stratifié

Le module ACP Pre a permis de définir l’empilement des couches composites, l’orientation des fibres, l’épaisseur de chaque pli ainsi que les propriétés orthotropes associées au matériau.

Attribution des propriétés mécaniques et thermiques

Les propriétés élastiques, thermiques et hygrothermiques du composite ont été introduites dans le modèle. Ces propriétés tiennent compte de la différence de comportement entre la matrice polymère et les fibres de renforcement.

Maillage par éléments finis

La structure a été discrétisée à l’aide d’un maillage régulier permettant une bonne représentation des champs de déplacement et des zones de concentration de contraintes. Le raffinement du maillage permet d’améliorer la précision des résultats, notamment pour l’analyse des modes propres et des zones critiques.

Application des conditions aux limites et chargements environnementaux 

Des conditions aux limites mécaniques ont été appliquées aux bords de la plaque afin de représenter les contraintes de fixation. Les effets thermiques et hygroscopiques ont été introduits sous forme de variations de température et de prise d’humidité, représentant les conditions de vieillissement hygrothermique.

Analyse modale et analyse de stabilité 

Une analyse modale a été réalisée pour déterminer les premières formes propres de vibration de la plaque composite. Ensuite, une étude de flambement thermique a permis d’estimer la température critique à partir de laquelle la structure perd sa stabilité.

Validation numérique par comparaison MATLAB/ANSYS 

Les résultats obtenus sous ANSYS ont été comparés avec ceux calculés par MATLAB. Cette comparaison permet de vérifier la cohérence du modèle numérique et de confirmer la fiabilité des tendances observées.

Trois scénarios principaux ont été considérés :

• exposition à une humidité constante ;
• cycles thermiques alternés ;
• couplage simultané humidité-température.

Cette méthodologie permet d’étudier séparément et conjointement les effets de l’humidité et de la température sur la réponse mécanique globale du composite.

Résultats et discussion

Les résultats numériques obtenus mettent en évidence l’influence significative du vieillissement hygrothermique sur le comportement mécanique et la stabilité des structures composites à matrice polymère.

Analyse du modèle numérique

Le schéma de simulation sous ANSYS Workbench montre le couplage entre le module ACP Pre et l’analyse Modal. Cette organisation permet de transférer correctement les informations relatives au stratifié composite — orientation des fibres, propriétés des plis et empilement — vers le module de calcul mécanique.

Cette étape confirme que le modèle composite a été correctement préparé pour l’analyse par éléments finis. Les étapes de géométrie, maillage, définition des matériaux, conditions aux limites et résolution sont validées dans l’environnement ANSYS.

Analyse des formes propres

Les figures représentant les premiers modes propres montrent différentes formes de déformation de la plaque composite. Le premier mode présente généralement une déformation globale centrée, caractéristique d’un comportement de flexion dominant. Les modes suivants montrent des formes plus complexes avec plusieurs zones de déplacement positif et négatif.

Les résultats montrent que :

• le premier mode propre est dominé par une déformation globale de la plaque ;
• les modes supérieurs présentent plusieurs lobes de déformation ;
• l’augmentation du numéro de mode entraîne une distribution plus complexe des déplacements ;
• les zones de déplacement maximal apparaissent principalement au centre ou dans des régions localisées de la plaque.

Ces formes propres sont importantes, car elles permettent d’identifier les zones les plus sensibles à l’instabilité mécanique. Dans le contexte du vieillissement hygrothermique, ces zones peuvent devenir critiques lorsque les propriétés de la matrice diminuent sous l’effet de l’humidité et de la température.

Distribution des déplacements

Les cartes de déplacement obtenues sous ANSYS montrent une répartition non uniforme des déformations dans la plaque composite. Les zones en rouge correspondent aux déplacements maximaux, tandis que les zones en bleu représentent les déplacements minimaux.

Les résultats indiquent que :

• les déplacements maximaux apparaissent dans les zones centrales ou dans les régions modales localisées ;
• les bords de la plaque présentent généralement des déplacements plus faibles en raison des conditions aux limites ;
• les modes supérieurs génèrent plusieurs zones de déplacement maximal ;
• la distribution des déplacements dépend fortement de l’orientation des fibres et de l’empilement du stratifié.

Ces observations confirment le caractère anisotrope du matériau composite. Contrairement aux matériaux isotropes, le comportement mécanique dépend fortement de la direction des fibres, de la séquence d’empilement et des conditions environnementales.

Influence de l’humidité sur les propriétés mécaniques

L’absorption d’humidité provoque une dégradation progressive de la matrice polymère. L’eau pénètre principalement dans la matrice et peut modifier les interactions entre les chaînes polymères. Cela entraîne une diminution de la rigidité globale du composite.

L’effet de l’humidité se traduit par :

• une réduction du module d’élasticité ;
• une augmentation des déformations ;
• une diminution de la résistance mécanique ;
• une baisse de l’efficacité du transfert de charge entre fibre et matrice ;
• une fragilisation progressive de l’interface fibre/matrice.

La matrice étant plus sensible à l’humidité que les fibres, le vieillissement hygrothermique provoque une hétérogénéité locale du comportement mécanique. Cette différence de sensibilité peut entraîner des concentrations de contraintes à l’interface fibre/matrice.

Influence de la température et des cycles thermiques

La température joue un rôle accélérateur dans le processus de vieillissement. Lorsque la température augmente, la mobilité moléculaire dans la matrice devient plus importante, ce qui facilite la diffusion de l’eau à l’intérieur du matériau.

Les cycles thermiques provoquent également des dilatations et contractions successives. Comme les fibres et la matrice ne possèdent pas les mêmes coefficients de dilatation thermique, des contraintes internes apparaissent dans le composite.

Ces contraintes peuvent provoquer :

• l’apparition de microfissures ;
• une décohésion progressive à l’interface fibre/matrice ;
• une perte de rigidité ;
• une réduction de la stabilité thermique ;
• une diminution de la durée de vie de la structure.

Ainsi, l’effet combiné de l’humidité et de la température est plus sévère que l’effet de chaque facteur pris séparément.

Température critique de flambement

La courbe de température critique de flambement montre l’évolution de (Delta T) en fonction du rapport géométrique (a/h). Les résultats obtenus par MATLAB et ANSYS présentent une tendance similaire, ce qui confirme la cohérence entre les deux approches numériques.

On observe que :

• la température critique augmente avec le rapport (a/h) ;
• l’augmentation est rapide pour les faibles valeurs de (a/h) ;
• au-delà d’une certaine valeur, la courbe tend vers un palier ;
• les résultats MATLAB sont légèrement supérieurs à ceux obtenus sous ANSYS ;
• l’écart entre les deux méthodes reste faible, ce qui valide le modèle éléments finis.

Cette évolution indique que la géométrie de la plaque influence fortement sa stabilité thermique. Les plaques plus minces ou présentant un rapport géométrique élevé peuvent montrer une sensibilité particulière au flambement sous chargement thermique.

La légère différence entre MATLAB et ANSYS peut être expliquée par :

• les hypothèses analytiques utilisées dans le calcul MATLAB ;
• le type d’éléments finis utilisé dans ANSYS ;
• les conditions aux limites numériques ;
• la discrétisation du maillage ;
• les approximations liées au modèle de stratifié composite.

Malgré ces différences, la bonne concordance des courbes confirme la validité de l’approche numérique.

Effet couplé hygrothermique

Le couplage humidité-température représente le cas le plus critique pour le composite. L’humidité réduit les propriétés mécaniques de la matrice, tandis que la température accélère la diffusion de l’eau et augmente les contraintes thermiques internes.

L’effet combiné provoque :

• une diminution progressive de la rigidité ;
• une augmentation des déplacements ;
• une réduction de la température critique de stabilité ;
• une concentration des contraintes à l’interface fibre/matrice ;
• une augmentation du risque de microfissuration ;
• une perte progressive de l’intégrité mécanique du composite.

Les résultats numériques permettent donc d’identifier les zones critiques avant l’apparition d’une défaillance macroscopique. Cette capacité prédictive est essentielle pour les structures composites utilisées dans des environnements sévères.

Applications industrielles

La méthodologie numérique développée présente un intérêt important pour plusieurs secteurs industriels utilisant des matériaux composites à matrice polymère.

Elle peut être appliquée notamment dans :

• l’aéronautique, pour prédire la durée de vie des panneaux composites soumis à l’humidité et aux variations de température ;
• le secteur naval et offshore, où les structures composites sont exposées à des milieux humides et salins ;
• l’énergie éolienne, pour étudier le vieillissement des pales soumises aux cycles thermiques et à l’humidité atmosphérique ;
• l’automobile, pour optimiser les pièces composites exposées à des conditions climatiques variables ;
• le génie civil, pour évaluer la durabilité des renforts composites utilisés dans les structures exposées à l’environnement.

L’utilisation de la simulation numérique permet de réduire le nombre d’essais expérimentaux longs et coûteux. Elle permet également d’optimiser la conception des structures composites dès les premières étapes du développement.

Cette approche peut ainsi contribuer à :

• améliorer la fiabilité des structures composites ;
• augmenter leur durée de vie ;
• réduire les coûts de maintenance ;
• identifier les zones sensibles à la dégradation ;
• proposer des stratifiés plus résistants aux environnements hygrothermiques.

Conclusion

Cette étude numérique a permis d’analyser le vieillissement hygrothermique des matériaux composites à matrice polymère en tenant compte de l’absorption d’humidité, des variations thermiques et de leurs effets sur les propriétés mécaniques.

La modélisation par éléments finis réalisée sous ANSYS, associée au module ACP pour la définition du stratifié composite, a permis d’étudier la réponse modale, la distribution des déplacements et la stabilité thermique de la structure. Les résultats montrent que l’humidité entraîne une diminution progressive de la rigidité mécanique, tandis que la température accélère la diffusion de l’eau dans la matrice et favorise l’apparition de contraintes internes.

L’analyse des formes propres montre que les modes de déformation deviennent plus complexes lorsque le numéro de mode augmente. Les cartes de déplacement mettent en évidence des zones critiques où les déformations sont maximales. Ces zones peuvent constituer des sites privilégiés pour l’initiation de microfissures ou de décohésion interfaciale.

La comparaison entre MATLAB et ANSYS pour la température critique de flambement montre une bonne concordance des résultats. Cette comparaison confirme la fiabilité du modèle numérique proposé. Les faibles écarts observés peuvent être attribués aux différences entre les hypothèses analytiques et la discrétisation par éléments finis.

Les résultats confirment que le couplage humidité-température constitue un facteur majeur de dégradation des matériaux composites. L’interface fibre/matrice apparaît comme une zone particulièrement sensible, car elle assure le transfert de charge entre les constituants du composite. Sa dégradation peut entraîner une perte progressive des performances mécaniques.

En conclusion, l’approche numérique développée constitue un outil efficace pour prédire le comportement et la durée de vie des structures composites exposées à des environnements hygrothermiques sévères. Elle permet d’identifier les zones critiques, d’évaluer la stabilité thermique et d’optimiser la conception des stratifiés composites.

Les travaux futurs pourront intégrer des modèles de dommage plus avancés, prenant en compte :

• l’initiation et la propagation des fissures ;
• la dégradation progressive de l’interface fibre/matrice ;
• les effets irréversibles de l’absorption d’humidité ;
• le vieillissement à long terme sous cycles thermiques répétés ;
• la validation expérimentale du modèle numérique.

Ainsi, cette étude ouvre la voie à une meilleure compréhension du comportement durable des matériaux composites à matrice polymère dans les conditions réelles de service.