Analyse des Matériaux Gradués Fonctionnellement (FGM) Poreux dans les Applications à Haute Température : Modélisation et Performance

42 sources13 sectionsTexte pas encore vérifiéThis text could not be checked. No percentage was generated.

Introduction aux Matériaux Gradués Fonctionnellement (FGM) et leurs Applications à Haute Température

Les matériaux à gradient fonctionnel, qu'on appelle souvent FGM (Functionally Graded Materials), forment une catégorie de composites hétérogènes assez particulière. Leur truc, c'est une variation continue de la composition et de la structure à travers leur masse . Cette transition progressive permet de leur tailler des propriétés sur mesure, idéales pour des usages industriels où les exigences sont souvent contradictoires au sein d'une même pièce . Comparés aux matériaux homogènes classiques, les FGM se distinguent par de meilleures performances mécaniques, une résistance accrue, une ductilité plus élevée, et des propriétés de surface au top . Leur microstructure graduée a l'avantage de gommer les interfaces nettes et abruptes qu'on trouve dans les assemblages traditionnels, les remplaçant par une transition douce d'un matériau à l'autre .

L'idée des FGM a germé au Japon en 1983. Le but ? Répondre aux contraintes thermomécaniques extrêmes subies par des composants comme les moteurs de fusées et les turboréacteurs . Les composites stratifiés classiques ont souvent un talon d'Achille : le délaminage. C'est la séparation des fibres de la matrice, causée par des taux d'expansion thermique différents quand il fait chaud . En introduisant un gradient de composition continu, sans coupure brutale d'interface, les FGM évitent ce défaut critique. Ils maintiennent ainsi les performances des matériaux, même dans des environnements ultra-agressifs .

Ces caractéristiques font des FGM des solutions privilégiées pour la régulation des transferts thermiques et l'atténuation des contraintes thermiques induites . Les fortes pentes de température – c'est-à-dire des différences de température importantes au sein d'un matériau, générant des contraintes thermiques majeures, un transfert thermique non linéaire et des propriétés variant avec la température – sont courantes dans l'aérospatial, les systèmes énergétiques et l'électronique de haute puissance . Dès leurs débuts, les FGM ont servi de barrières thermiques essentielles pour des applications critiques : pense aux véhicules de rentrée atmosphérique ou aux réacteurs à fusion nucléaire .

Dans les environnements thermiques difficiles, l'intégration de composés spécifiques montre bien l'efficacité des gradients fonctionnels. Prenez par exemple l'intermétallique Ni3Al. Par exemple, l'intermétallique Ni3Al offre une conductivité thermique élevée conjuguée à une résistance remarquable à la traction et à la compression à des températures situées entre 650 et 1100 °C . En plus, le Ni3Al résiste bien au fluage, à la fatigue, à l'usure à haute température, et aux atmosphères riches en oxygène et en carbone. C'est donc un matériau tout indiqué pour les conditions corrosives et chaudes .

Les sources consultées ne fournissent pas de données empiriques sur l'impact des micro-porosités internes sur les performances mécaniques ou thermiques de ces structures. De même, les éléments analysés n'abordent pas l'évolution des théories cinématiques structurales : ni la théorie classique des plaques (CPT), ni le modèle de déformation par cisaillement de premier ordre (FSDT), ni les théories de déformation par cisaillement d'ordre supérieur (HSDT), avec ou sans facteurs de correction. Les fondations élastiques de type Winkler, Pasternak ou Kerr, tout comme les modèles non linéaires et les approches computationnelles sophistiquées de 2023 et 2024, ne sont pas non plus mentionnés dans les documents recensés.

Avantages des FGM par Rapport aux Composites Stratifiés Traditionnels

L'évolution des sciences des matériaux nous a vus passer tranquillement des métaux classiques aux composites perfectionnés et aux matériaux "intelligents". Pourquoi ? On cherchait de meilleures performances structurelles . Mais voilà, si les composites traditionnels sont globalement excellents, ils montrent vite leurs limites dans les environnements difficiles . Leur gros point faible, c'est le délaminage . Ça arrive quand les fibres et la matrice se séparent sous l'effet de la chaleur, tout simplement parce qu'elles ne se dilatent pas de la même façon .

Pour régler le problème des défaillances des composites stratifiés classiques, on a vu apparaître une nouvelle génération : les matériaux à gradient fonctionnel, ou FGM (Functionally Graded Materials). Ils sont nés au Japon dans les années 1980 . Ces matériaux d'ingénierie avancés sont conçus pour tenir le coup dans des conditions très dures, sans que leurs propriétés ne se dégradent. Résultat : on réduit drastiquement les risques de rupture pendant leur durée de vie . Pas étonnant que la communauté scientifique et les industries de pointe s'y intéressent de plus en plus .

L'avantage fondamental des FGM sur les composites classiques, c'est leur résistance exceptionnelle aux ruptures causées par les interfaces nettes . Dans les composites stratifiés, la transition brutale entre les couches crée des concentrations de contraintes qui peuvent déclencher le délaminage . Les FGM, eux, suppriment carrément ces interfaces. Comment ? En introduisant progressivement un second matériau dans le premier . Les interfaces nettes font place à une variation continue, une transition tout en douceur entre deux matériaux différents . Ce gradient de composition élimine les discontinuités qui sabotaient jusqu'ici les performances des mélanges homogènes dans les composites classiques .

Cette gradation compositionnelle donne une microstructure non homogène, mais contrôlée. Les FGM affichent ainsi des propriétés mécaniques et thermiques qui varient de manière continue sur leur épaisseur . Grâce à cette souplesse de conception, on peut adapter leurs propriétés aux contraintes spécifiques de nombreuses applications . Au début, les FGM ont surtout servi de barrières thermiques pour les rentrées atmosphériques dans l'aérospatiale, et pour les réacteurs à fusion dans le nucléaire. Il faut dire que ces composants doivent encaisser des températures extrêmes . Mais aujourd'hui, les FGM intéressent aussi les domaines biomédical, automobile, de la défense, les applications pare-balles et les retardateurs de flamme .

Impact de la Micro-porosité sur les Propriétés Mécaniques des FGM

La présence de défauts poreux dans les matériaux à gradient fonctionnel (FGM) modifie profondément leur comportement mécanique global. Ça affecte surtout la distribution des contraintes et déplace le plan neutre de la structure . Ces micro-porosités internes apparaissent souvent à cause des processus de fabrication, particulièrement lors du frittage sans pression de poudres métalliques. La formation de ces micro-porosités internes découle fréquemment des processus de fabrication, notamment lors du frittage sans pression de poudres métalliques, une méthode couramment privilégiée pour sa simplicité relative et son faible coût . Le processus lui-même crée des matériaux poreux lors de la solidification et de la cristallisation . De plus, la température maximale de frittage est limitée par le point de fusion de chaque composant. Ça peut résulter en un frittage imparfait et donc des taux de porosité plus élevés .

Ces porosités ont un impact direct sur la masse volumique et les propriétés mécaniques du composite final. Ces porosités influencent directement la masse volumique ainsi que les propriétés mécaniques du composite final, rendant indispensable la prédiction exacte des variations de son module d'élasticité . Pour comprendre cette dépendance, plusieurs approches de modélisation théorique ont été développées. Il y a un premier modèle basé sur la fraction volumique des pores dans chaque composant du composite. Un deuxième modèle suppose que la porosité n'existe que dans le matériau de la matrice. Et un troisième modèle utilise un ajustement de courbe empirique des modules d'élasticité des matériaux purs en fonction de la porosité . Comparés aux valeurs mesurées sur des structures FGM de Ni/Al2O3, les résultats du troisième modèle correspondent le mieux .

Côté résistance structurelle, les défauts de porosité diminuent la résistance à la flexion des FGM tout en augmentant leur résistance à la compression. Ils peuvent aussi causer des effets transversaux importants . Pour analyser ces effets transversaux, il faut des modèles cinématiques adéquats. L'analyse de ces effets transversaux nécessite la prise en compte de modèles cinématiques adaptés, capables de dissocier deux effets de Poisson distincts et de mettre en évidence le rôle critique des variations du coefficient de Poisson à travers l'épaisseur de la plaque . La modélisation analytique de ces structures, dont les propriétés varient avec l'épaisseur selon des lois exponentielles ou de puissance, permet de garantir que les contraintes de cisaillement sont nulles sur les surfaces supérieure et inférieure, sans avoir besoin d'un facteur de correction de cisaillement . Dans ces cas, le réseau de pores peut être considéré comme vide ou rempli d'air à basse pression. Les équations différentielles d'équilibre sont résolues par la méthode de Navier .

Enfin, l'emplacement de la porosité est déterminant pour gérer les champs de contraintes mécaniques . Si la porosité se concentre au centre de la plaque ou dans des couches spécifiques, cela réduit les contraintes normales de traction le long de l'axe longitudinal . Cela suggère qu'on peut utiliser un agencement stratégique des pores pour atténuer les concentrations de contraintes. L'objectif est de préserver l'intégrité structurelle et d'optimiser les performances des matériaux à gradient fonctionnel soumis à des contraintes mécaniques .

Modélisation de la Porosité et Prédiction des Propriétés des FGM

Les matériaux à gradient fonctionnel (FGM) forment une catégorie de composites avancés. Leur particularité ? Les matériaux à gradient fonctionnel (FGM) constituent une classe de matériaux composites avancés caractérisés par une variation spatiale continue de leur composition et de leurs propriétés structurales . Selon l'usage prévu, on peut les classer : soit par un gradient de matière, soit par un gradient de porosité, soit encore par un gradient microstructural . Les structures poreuses à gradient fonctionnel affichent une répartition graduée des pores. Les structures poreuses à gradient fonctionnel présentent une distribution graduée de la porosité dans leur volume, ce qui leur confère des caractéristiques recherchées telles que la légèreté et une excellente capacité d'absorption d'énergie .

Toutefois, l'apparition de micropores internes résulte aussi des procédés de fabrication . En effet, quand on fabrique des FGM par frittage de poudres métalliques sans pression, on obtient inévitablement des matériaux poreux. C'est une conséquence directe des mécanismes de solidification et de cristallisation . De plus, la température de frittage est limitée par le point de fusion de chaque composant. Ce qui peut mener à un frittage incomplet et donc à une porosité plus importante . Cette porosité résiduelle pèse directement sur la densité et les performances mécaniques globales du composite fini .

Pour évaluer et prédire avec justesse les propriétés réelles des structures FGM poreuses, il faut des formulations mathématiques bien adaptées . Dans la littérature scientifique, on estime souvent les propriétés matérielles des coques poreuses en FGM à l'aide de la règle des mélanges. On y ajoute un terme supplémentaire qui représente la porosité dans l'épaisseur . La règle des mélanges modifiée, ainsi que la loi généralisée des mélanges, servent à prédire les propriétés de ces composites. La loi généralisée prend notamment en compte des facteurs comme la forme des grains et la répartition des phases .

Pour estimer le module d'élasticité en fonction de la porosité, plusieurs approches existent. Il y a un premier modèle basé sur la fraction volumique des pores dans chaque matériau du composite. Un deuxième modèle suppose que la porosité se trouve uniquement dans la matrice. Et un troisième modèle intègre un ajustement empirique des modules d'élasticité des matériaux purs en fonction de la porosité . Des comparaisons expérimentales menées sur des FGM en Ni/Al2O3 ont montré que ce troisième modèle correspondait le mieux aux valeurs mesurées .

Les propriétés matérielles d'une plaque FGM poreuse évoluent progressivement selon l'épaisseur, suivant une loi de puissance . Pour aller au-delà des formulations classiques qui supposent une dépendance exclusive à la règle des mélanges, des modèles analytiques établissent un lien direct entre la porosité et la fraction volumique du matériau . Deux schémas distincts de porosité sont ainsi examinés : la porosité dépendant de la fraction volumique (VFD) et la porosité dépendant de la règle des mélanges (RMD) . Ces schémas intègrent différents profils de distribution dans la géométrie. On trouve ainsi quatre types distincts : Even, Uneven, Linear (1) et Linear (2) . Ou encore cinq types de fonctions de distribution de porosité appliqués à l'analyse des plaques . Dans le cas particulier des nanomatériaux composites poreux à gradient fonctionnel, il a été prouvé que les nanotubes de carbone et les distributions de porosité ont des effets opposés et distincts sur le comportement mécanique des structures .

Le contexte opérationnel des structures en FGM fait souvent intervenir des fondations élastiques ou des gradients thermiques non linéaires. Or, les sources consultées ne décrivent pas explicitement la prise en compte des fondations de Kerr ou des gradients thermiques non linéaires dans ces schémas de porosité. De même, les documents disponibles n'expliquent pas en détail les formulations cinématiques raffinées à quatre variables d'ordre logarithmique ou hyperbolique pour ces lois de mélange spécifiques.

Évolution des Théories de Modélisation des Plaques : De CPT à HSDT

L'étude du comportement mécanique des structures en plaques s'appuie sur une hiérarchie de modèles cinématiques issus de la réduction 2D de l'élasticité 3D. On y retrouve la théorie classique des plaques (CPT), la théorie du premier ordre et les théories d'ordre supérieur . La CPT, aussi appelée théorie de Kirchhoff-Love, est le cadre le plus élémentaire. Le cadre le plus élémentaire est la théorie classique des plaques (CPT), également désignée sous le nom de théorie de Kirchhoff-Love, dont l'extension aux plaques composites stratifiées constitue la théorie classique des plaques stratifiées (CLPT) . Ce formalisme suppose que les segments rectilignes initiaux restent droits et perpendiculaires au plan moyen après déformation. Les déformations transversales de cisaillement et normales s'annulent donc . La CPT donne des résultats satisfaisants pour les plaques fines, là où ces effets sont négligeables. Mais elle perd en précision quand l'épaisseur augmente. L'écart avec les théories plus poussées devient alors notable .

Pour corriger l'omission du cisaillement transversal, on a introduit la théorie de déformation par cisaillement de premier ordre (FSDT), issue des travaux de Reissner et Mindlin . Ce modèle prend en compte les déformations de cisaillement transversal. Il formule une variation linéaire des déplacements dans l'épaisseur de la plaque . Cependant, l'hypothèse cinématique linéaire de la FSDT crée une déformation de cisaillement uniforme. Elle viole les conditions aux limites de contraintes nulles sur les surfaces libres de la plaque . La FSDT requiert par conséquent l'intervention d'un facteur de correction de cisaillement afin d'ajuster l'évaluation des contraintes sur ces surfaces libres . Déterminer ce facteur précisément s'avère très ardu et complexe en pratique. Sa valeur dépend de multiples paramètres : géométrie, type de chargement, conditions aux limites .

Face à ces limites, les théories de déformation par cisaillement d'ordre supérieur (HSDTs) ont vu le jour. Elles visent à dépasser les contraintes de la FSDT et à supprimer le besoin d'un facteur de correction de cisaillement . Les formulations HSDT intègrent des termes d'ordre supérieur dans le champ des déplacements, selon la coordonnée de l'épaisseur. Elles utilisent des séries de puissances ou des développements en série de Taylor jusqu'au troisième ordre . Cette modélisation capte avec précision les variations non linéaires des déplacements dans le plan à travers l'épaisseur. Cette modélisation permet de capturer avec précision les variations non linéaires des déplacements dans le plan à travers l'épaisseur et de réconcilier les écarts entre les distributions réelles et supposées des déformations de cisaillement . Parmi les HSDTs, le modèle de Reddy, basé sur la théorie de déformation par cisaillement au troisième ordre (TSDT), est souvent utilisé. Son équilibre entre simplicité et précision le rend populaire . Enfin, des théories de plaques raffinées (RPTs) ont été conçues pour simplifier les HSDTs. Elles réduisent le nombre de fonctions inconnues tout en gardant une exactitude suffisante pour l'analyse des structures .

Nécessité et Avantages des Théories de Déformation de Cisaillement d'Ordre Supérieur (HSDT)

Le développement des théories structurelles raffinées, incluant les théories de déformation de cisaillement d'ordre supérieur (HSDT), découle directement des exigences de la pratique de l'ingénierie, qui cherche à modéliser fidèlement les comportements mécaniques tout en améliorant l'efficacité de résolution par rapport aux modèles tridimensionnels complets . Au sein des approches de type monocouche équivalente (ESL), les modèles bidimensionnels sont dérivés de la théorie de l'élasticité 3D en introduisant des hypothèses cinématiques ou d'état de contrainte selon l'épaisseur de la structure . Alors que la théorie classique des plaques (CPT) néglige totalement les déformations de cisaillement transversal, conduisant à des résultats satisfaisants uniquement pour les plaques minces mais insuffisants pour les plaques épaisses, la théorie de déformation de cisaillement du premier ordre (FSDT) intègre une variation linéaire des déplacements plans . Néanmoins, la FSDT suppose une déformation de cisaillement transversal uniforme à travers l'épaisseur et nécessite l'introduction d'un facteur de correction de cisaillement pour approximer la distribution parabolique réelle des contraintes de cisaillement transversal et satisfaire les conditions sur les surfaces libres .

La détermination de ce facteur de correction de cisaillement pose des difficultés majeures en raison de sa dépendance vis-à-vis des conditions aux limites, des charges appliquées et de la géométrie de la structure . Les théories HSDT s'avèrent dès lors indispensables pour modéliser avec précision les plaques en matériaux composites et en matériaux à gradient fonctionnel (FGM) . Elles sont formulées en introduisant des composantes de déplacement sous forme de séries de puissances dans la coordonnée de l'épaisseur, ce qui permet de concilier les écarts observés entre les distributions effectives et supposées du cisaillement . L'utilisation de polynômes d'ordre supérieur pour les déplacements au travers de l'épaisseur génère des variations quadratiques des déformations et des contraintes de cisaillement transversal . Par conséquent, les modèles HSDT représentent un champ de cisaillement transversal non uniforme et éliminent totalement le recours à un facteur empirique de correction de cisaillement .

Par rapport aux formulations de niveau inférieur, les théories HSDT capturent plus efficacement les effets de la déformation par cisaillement, fournissant une meilleure prédiction de la déformation de cisaillement transversal ainsi que des déformations standard dans les structures FGM . Elles sont considérées comme les plus précises de la famille des modèles monocouches équivalents pour les structures fonctionnellement graduées . De plus, certaines variantes des HSDT peuvent prendre en compte les effets d'étirement selon l'épaisseur, qui ne sont alors pas négligés dans le calcul de la déformation normale transversale . Bien qu'elles engendrent des coûts de calcul supérieurs à ceux des théories d'ordre inférieur, les HSDT présentent une très bonne concordance avec les solutions exactes issues de la théorie de l'élasticité 3D, tout en constituant une approche analytique et numérique plus simple et moins exigeante en ressources que les formulations tridimensionnelles pures .

Avancées Récentes : Théories de Déformation de Cisaillement Non-Polynomiales et Raffinées pour les FGM Poreux

L'évolution récente de la modélisation structurelle des plaques en matériaux à gradient fonctionnel (FGM) intégrant des porosités s'est concentrée sur le développement de théories d'ordre supérieur (HSDT) et non-polynomiales (NPSDT) afin de surmonter les limitations des approches classiques . Les recherches menées sur la période 2023-2024 appliquent ces théories raffinées aux plaques FGM poreuses exposées à des environnements thermiques . Contrairement aux formulations du premier ordre qui imposent un facteur de correction de cisaillement, ces théories cinématiques enrichies satisfont rigoureusement les conditions aux limites de contraintes nulles sur les surfaces libres et reproduisent une distribution parabolique des déformations et des contraintes de cisaillement transversal dans l'épaisseur de la structure .

Pour réduire les coûts de calcul associés aux modèles cinématiques standards à cinq variables ou plus, plusieurs approches raffinées à quatre variables inconnues ont été formulées . Dans ces cadres cinématiques, le déplacement transversal est scindé en composantes de flexion et de cisaillement . Ces théories à quatre variables capturent les effets de déformation par cisaillement sans recours aux facteurs de correction transversaux, comme démontré lors de l'analyse vibratoire de plaques piézoélectriques à porosités ou de poutres FGM sous chargements thermiques uniforme, linéaire, non linéaire ou sinusoïdal . Pour des plaques sandwichs FGM avec divers profils de porosité, un modèle à quatre variables donne des résultats très proches des formulations d'ordre supérieur du type TSDT tout en réduisant le nombre d'inconnues . Parallèlement, des approches quasi-tridimensionnelles (quasi-3D) intègrent l'effet d'étirement suivant l'épaisseur tout en diminuant le nombre d'équations gouvernantes pour la flexion statique de plaques FGM poreuses .

L'introduction de fonctions de forme non-polynomiales constitue une autre avancée majeure dans l'évaluation thermo-mécanique précise des structures poreuses. Une théorie HSDT novatrice formulée au moyen d'une fonction logarithmique-hyperbolique a été établie pour l'analyse des vibrations libres thermo-mécaniques de plaques solaires FGM poreuses soumises à un environnement thermique . Cette formulation, régie par quatre inconnues, s'affranchit de tout facteur de correction de cisaillement et fournit une précision supérieure à plusieurs HSDT par rapport aux solutions exactes tridimensionnelles . En outre, des comparaisons directes entre les théories de cisaillement du troisième ordre (TSDT) et les modèles non-polynomiaux (NPSDT) ont été établies au moyen de modèles éléments finis C0 à pénalisation pour la flexion de plaques FGM sous sollicitations mécaniques et thermiques .

La prise en compte des interactions avec les systèmes de support élastique constitue également un axe central des travaux récents. Les plaques FGM poreuses soumises à des effets thermiques et dotées de propriétés dépendant de la température sont fréquemment modélisées sur des fondations élastiques de type Kerr à trois paramètres . Ces études montrent que les fondations de Kerr permettent d'ajuster le comportement dynamique des structures soumises à un champ thermique satisfaisant l'équation de transfert de chaleur en régime stationnaire . Les comportements statiques de plaques saturées poreuses FGM (FGSP) sont également investigués en association avec la fondation élastique de Pasternak au travers de théories quasi-3D .

Enfin, diverses techniques d'approximation numérique ont été adaptées à ces théories de déformation avancées. L'analyse isogéométrique (IGA) reposant sur les fonctions de base B-splines rationnelles non uniformes (NURBS) a été couplée à la théorie TSDT pour modéliser exactement la géométrie et respecter la continuité C1 requise sans facteur de correction sous fondation de Kerr . Des formulations non linéaires par éléments finis fondées sur la théorie TSDT ont également été appliquées aux micro-plaques circulaires et annulaires sous gradients thermiques et mécaniques en incorporant la théorie des contraintes de couple modifiées et divers profils de distribution de pores .

Modélisation de l'Interaction avec les Systèmes de Support Élastique

Pour bien cerner le comportement mécanique des structures en matériaux à gradient fonctionnel (FGM) poreux, il faut impérativement tenir compte de leurs conditions d'appui, surtout quand elles reposent sur un milieu continu déformable . Dans la littérature, ce sont les modèles de fondation élastique de Winkler et de Pasternak qu'on retrouve le plus souvent pour décrire cette interaction . Le modèle classique de Winkler, c'est un seul paramètre qui représente un lit de ressorts indépendants. La formulation de Pasternak, elle, utilise deux paramètres : la rigidité de Winkler (k1) et la rigidité de cisaillement de Pasternak (k2) . Cette dernière, tout comme des modèles plus sophistiqués, offre une prédiction plus proche de la réalité du comportement structural que le modèle de Winkler, car elle prend en compte le cisaillement dans le milieu support . On a même développé des extensions du modèle de Winkler-Pasternak, qui ajoutent un paramètre de masse à la fondation, sous le nom de modèle EFTIM .

Pour dépasser les limites des approches à un ou deux paramètres, le modèle de fondation de Kerr est aussi utilisé pour analyser les poutres, les plaques et les nanostructures FGM . Intégrer une fondation élastique de type Kerr permet d'ajuster finement la réponse dynamique et le comportement vibratoire libre des plaques FGM poreuses, particulièrement face aux effets thermiques et à des propriétés matérielles réparties de façon complexe . Concernant le flambement, on étudie les nanoplates sandwich FGM poreuses sur fondation de Kerr, en considérant des propriétés dépendantes de la température et suivant une loi de puissance . De manière générale, augmenter les paramètres d'une fondation élastique, comme les constantes de rappel élastique et de cisaillement, renforce la rigidité de la structure. Cela fait généralement monter les fréquences propres et les charges critiques de flambement des structures FGM poreuses .

La modélisation cinématique des plaques et poutres FGM poreuses interagissant avec ces fondations s'appuie sur diverses théories mécaniques. Beaucoup de recherches emploient la théorie de déformation par cisaillement de premier ordre (FSDT) . D'autres études se tournent vers des théories de déformation par cisaillement d'ordre élevé (HSDT) ou des théories quasi-tridimensionnelles (quasi-3D) pour examiner la propagation des ondes et les vibrations libres. Souvent, les propriétés matérielles sont décrites par une loi de puissance de la fraction volumique sur l'épaisseur . Plus récemment, une nouvelle théorie de plaque à quatre variables a vu le jour pour analyser les vibrations libres de plaques FGM avec porosités sur fondations de Winkler, Pasternak ou Kerr .

Sur le plan des calculs et du numérique, plusieurs méthodes ont été mises au point pour gérer de manière fiable l'interaction avec le milieu élastique. La méthode des éléments finis lissés basée sur les arêtes (ES-FEM), associée à la technique d'interpolation mixte des composantes tensorielles (MITC), a été utilisée pour analyser des plaques FGM poreuses sur fondation élastique. Le but ? Éviter le verrouillage en cisaillement et améliorer la précision numérique . Par ailleurs, la méthode des fonctions R, couplée à la méthode variationnelle de Ritz, a permis de résoudre avec succès le problème des vibrations libres de plaques FGM poreuses à épaisseur variable sur des fondations élastiques biparamétriques .

Distribution de la Porosité et ses Effets sur la Performance

La présence de micro-porosités au sein des structures en matériaux à gradient fonctionnel (FGM) change radicalement leur comportement mécanique et dynamique. En gros, introduire des porosités réduit la rigidité globale de la structure, ce qui fait baisser les fréquences propres . La littérature distingue plusieurs manières dont ces porosités peuvent être réparties, que ce soit à travers l'épaisseur ou le volume des composants. Elles peuvent être uniformes (« even ») ou non uniformes (« uneven »), symétriques ou asymétriques . Des formulations mathématiques spécifiques ont aussi été mises au point pour analyser des distributions linéaires, non linéaires dites « Nonlinear (1) » et « Nonlinear (2) », ou encore trigonométriques .

L'impact de ces défauts dépend vraiment de la façon dont on modélise le gradient de vide dans l'espace. L'étude de différents schémas de répartition révèle des sensibilités variables : par exemple, le schéma « Porosity I » influence fortement les fréquences, tandis que « Porosity II » a l'effet le plus faible . Dans le cadre d'études vibratoires sur des plaques FGM posées sur des fondations élastiques, la distribution trigonométrique de la porosité montre de bonnes performances dynamiques. C'est parce qu'elle conserve la rigidité dans les zones de surface des plaques . En plus, les distributions non uniformes des propriétés des matériaux ou de la porosité peuvent provoquer un comportement de torsion inattendu sous contrainte .

Pour évaluer la performance mécanique, il faut des cadres cinématiques précis qui simulent les variations de déformation sans artifice. L'utilisation d'une théorie de déformation par cisaillement d'ordre supérieur, quasi-3D, permet d'intégrer des coefficients de porosité différents pour le module de Young et la masse volumique. Ces coefficients sont adaptés à chaque type de distribution. Ça autorise des comparaisons rigoureuses entre différentes répartitions, tout en maintenant une masse totale constante . Par ailleurs, appliquer un cadre théorique de déformation par cisaillement à quatre variables permet de comparer systématiquement l'effet de trois motifs distincts de répartition de porosité sur la flèche adimensionnelle . Les fonctions de déformation par cisaillement utilisées dans ces modèles prédisent avec précision comment la porosité affecte le comportement en déflexion des panneaux sandwichs .

Complexité Computationnelle et Simplifications dans la Modélisation des FGM Poreux

La modélisation numérique et analytique des structures à gradient fonctionnel (FGM) poreuses pose de sérieux défis, surtout à cause des coûts de calcul et de la difficulté à représenter fidèlement la microstructure . Dans la recherche actuelle, la complexité des calculs nécessaires est une vraie limite . Pourquoi ? Parce que ces plaques et coques FGM poreuses ont un comportement mécanique intrinsèque qui exige une analyse en 3D pour être bien décrit, alourdissant considérablement les besoins en calcul .

Quand on évalue les théories de déformation par cisaillement, on voit bien le compromis entre précision et charge de calcul. La théorie de premier ordre (FSDT), elle, est souvent choisie pour sa simplicité et ses faibles besoins en calcul . Néanmoins, la formulation classique de la FSDT requiert des facteurs de correction de cisaillement, dont la détermination et l'application s'avèrent impraticables pour les configurations de FGM poreux . À l'inverse, les théories d'ordre supérieur (HSDTs) évitent ces facteurs empiriques en représentant un cisaillement transversal non uniforme . Par contre, utiliser les HSDTs rend l'étude des plaques et coques FGM poreuses bien plus complexe que la FSDT, ce qui alourdit sérieusement la résolution numérique .

Ces exigences analytiques et de calcul s'accentuent encore pour des configurations complexes. Pensez à la prédiction précise des vibrations et de la propagation des ondes dans des poutres FGM poreuses bidirectionnelles . Là, la variation multidirectionnelle des propriétés matérielles dans la section droite interagit avec les effets combinés du cisaillement et de l'inertie de rotation. Ça modifie la réponse modale, rendant les modèles cinématiques plus compliqués . De plus, les études sur les formulations 3D indépendantes par éléments finis montrent l'importance de la discrétisation spatiale. On voit des analyses de convergence de maillages avec une variation de 0,290 % entre les grilles les plus fines .

À côté de ces contraintes cinématiques, simplifier la distribution de la porosité pose des problèmes critiques dans les approches habituelles. L'aléa géométrique local des mousses métalliques rend la prédiction des performances des structures poreuses très compliquée . La densité relative est souvent le seul facteur utilisé dans les simplifications analytiques courantes. Pourtant, le caractère aléatoire de la taille et de la forme des cellules internes influence clairement le comportement structural des matériaux poreux . Mesurer expérimentalement cette variabilité cellulaire reste d'une difficulté extrême .

Pour contourner ces limitations géométriques sans modéliser explicitement chaque cellule, d'autres voies ont vu le jour. On combine des éléments de volume représentatifs (RVE) avec des morphologies cellulaires aléatoires, qui imitent les mousses d'aluminium à cellules fermées, et des réseaux neuronaux convolutifs (CNN) . Ces outils d'apprentissage profond établissent une corrélation directe entre les caractéristiques mésoscopiques et les modules effectifs. L'erreur moyenne est de 5,92 % pour les modèles CNN sélectionnés . Toutefois, l'intervalle d'erreur de ces modèles d'estimation introduit une incertitude mécanique. Il faut la répercuter dans l'analyse globale, par exemple sous forme de nombres flous triangulaires biparamétriques . Les sources consultées ne donnent aucune information sur les limitations spécifiques des lois de distribution sigmoïde (S-FGM) ou exponentielle (P-FGM), ni sur des formulations particulières de gradients thermiques non linéaires dans ce contexte.

Cadres Cinématiques et Modélisation des Matériaux/Porosités Novateurs

L'élaboration de modèles cinématiques avancés est une étape clé pour décrire précisément le comportement structurel des composants à gradient fonctionnel (FGM) qui présentent des porosités. Les théories affinées à quatre variables marquent un progrès considérable par rapport aux formulations classiques. Elles permettent de saisir fidèlement les effets de la déformation par cisaillement sans avoir recours à des facteurs de correction . Ces cadres cinématiques dérivent rigoureusement les équations différentielles du comportement mécanique et les conditions aux limites via le principe de Hamilton . L'élimination du facteur de correction transversal simplifie la structure computationnelle tout en offrant une représentation réaliste de la distribution des contraintes sur l'épaisseur de l'élément structurel .

Concernant les poutres à gradient fonctionnel et poreuses, une théorie affinée de déformation par cisaillement à quatre variables considère un déplacement transversal constant. Elle y associe une variation d'ordre supérieur du déplacement axial selon l'épaisseur . Ces formulations permettent d'étudier la réponse vibratoire thermomécanique. C'est utile quand les propriétés des matériaux dépendent de la température et qu'ils sont soumis à diverses élévations thermiques — répartitions de température uniformes, linéaires, non linéaires ou sinusoïdales le long de l'axe vertical . Une approche analytique aide à déterminer les fréquences propres. Elle permet aussi d'examiner l'influence de l'exposant de la loi en puissance, de la fraction volumique de porosité, des différentes distributions de pores et des sollicitations thermiques .

Pour les structures bidimensionnelles, on utilise aussi un modèle de plaque affiné à quatre inconnues. Il sert notamment à l'analyse des vibrations électromécaniques des plaques piézoélectriques à gradient fonctionnel (FGP) poreuses . Les équations du mouvement, issues du principe de Hamilton et résolues analytiquement, tiennent compte de multiples conditions aux limites. Elles permettent d'évaluer l'impact combiné de la tension appliquée, de la distribution des porosités, de la gradation matérielle et des paramètres géométriques de la plaque .

Plus récemment, une théorie de déformation par cisaillement d'ordre supérieur (HSDT) a été établie. Elle utilise une fonction logarithmique-hyperbolique pour l'analyse des vibrations libres thermomécaniques de plaques solaires en FGM poreux, plongées dans un environnement thermique . Cette théorie, nommée 2D-LHHSDT, limite le champ de déplacement à quatre inconnues seulement. La troisième composante du déplacement est divisée en une partie flexion et une partie cisaillement, éliminant ainsi le besoin de facteurs de correction de cisaillement transversal . Comparée aux solutions tridimensionnelles exactes, la formulation 2D-LHHSDT donne de meilleurs résultats que diverses autres théories HSDT .

La modélisation matérielle s'appuie souvent sur des lois en puissance simples ou modifiées. Cela permet d'intégrer la gradation continue des constituants et l'apparition de micro-vides ou de défauts de fabrication . Les porosités peuvent être réparties uniformément ou non (types even et uneven) le long de l'épaisseur. Elles influencent notablement les caractéristiques mécaniques . De plus, les interactions avec le milieu environnant sont prises en compte via des modèles de fondations élastiques. Un exemple est la fondation de Kerr à trois paramètres. Elle permet d'étudier conjointement les effets de la dépendance thermique des propriétés, du paramètre de loi de puissance et de la fondation sur le comportement dynamique des structures solaires en FGM poreux .

Conditions Environnementales Réalistes et Défis Futurs

L'évaluation du comportement structural des matériaux à gradient fonctionnel (FGM) impose de considérer des conditions environnementales complexes et réalistes, incluant les interactions avec des fondations élastiques et de sévères sollicitations thermiques . Récemment, des études ont ainsi exploré le flambement de nanplaques sandwich poreuses en FGM, simplement appuyées et reposant sur une fondation de type Kerr . Dans ce contexte, les propriétés mécaniques et thermiques du matériau sont liées à la température et évoluent continuellement selon l'épaisseur . Afin de décrire cette variation continue des propriétés de la plaque sandwich FGM, un modèle de loi de puissance (power-law model) est couramment adopté . Les équations analytiques du système peuvent d'ailleurs être formulées via le principe variationnel de Hamilton. On y intègre les énergies de déformation dues aux charges thermiques, ainsi qu'une théorie de plaques à gradient de déformation non local d'ordre élevé, qui saisit les effets de cisaillement sans nécessiter de facteur de correction .

Ces modélisations environnementales sophistiquées incluent une solution exacte pour la variation non linéaire de température à travers l'épaisseur de la nanoplque sandwich. Elles prennent rigoureusement en compte la conductivité thermique, le paramètre d'inhomogénéité et les schémas d'agencement du sandwich . Les résultats numériques issus de ces approches révèlent l'influence combinée de nombreux paramètres physiques et géométriques sur la réponse au flambement. On pense ici à l'indice de fraction volumique, au paramètre non local, au paramètre d'échelle de longueur, au coefficient de porosité, à la fondation de Kerr et à la différence de température appliquée .

Certes, des progrès notables ont été accomplis dans la représentation des champs de cisaillement et des interactions élastiques. Cependant, des défis méthodologiques importants subsistent pour les recherches futures . Nombre d'approches actuelles imposent des limites et des vérifications de faisabilité numérique, mais elles négligent plusieurs contraintes techniques essentielles à l'ingénierie pratique . Notamment, les contraintes admissibles, les aspects économiques et les contraintes de fabricabilité spécifiques aux procédés industriels ne sont pas appliqués dans les modèles actuels; ils sont explicitement réservés à de futures extensions . Une optimisation pratique de ces structures devra donc impérativement intégrer des limites de porosité calibrées, des dimensions minimales fabricables, des critères de rupture et de ténacité, des exigences de stabilité sous axe faible ou de déversement latéral-torsionnel, ainsi que des contraintes strictes de procédé et de coût . Les sources disponibles ne fournissent pas d'autres détails sur l'application industrielle directe de ces structures.

Questions ouvertes

— Comment la distribution stochastique de la porosité affecte-t-elle précisément le comportement des structures poreuses FGM ?
— Quelles sont les contraintes de fabricabilité et de coût qui ne sont pas encore intégrées dans les modèles actuels ?
— Comment optimiser la distribution de la porosité pour maximiser la performance mécanique et thermique dans des conditions environnementales complexes ?