Modélisation Quasi-3D Thermomécanique des Matériaux à Gradient de Fonction (FGM) : Une Analyse des Théories de Déformation au Cisaillement et de l'Impact de la Porosité

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

Introduction aux Matériaux à Gradient de Fonction (FGM)

Les matériaux à gradient de fonction, qu'on appelle aussi FGM dans la littérature scientifique anglo-saxonne, sont une classe de composites avancés. Leurs propriétés changent en douceur, dans une ou plusieurs directions, à l'intérieur même de la structure . C'est ce changement progressif qui les distingue vraiment des composites stratifiés classiques. Dans ces derniers, le passage brutal entre les matériaux crée des points de stress et des discontinuités d'interface qui nuisent à la résistance mécanique . C'est pourquoi les FGM sont considérés comme des composites de seconde génération. Ils peuvent marier des propriétés apparemment contradictoires dans une seule pièce, améliorant ainsi leurs performances dans des environnements difficiles, sans défaillance ni dégradation notable de leurs composants .

Tout a commencé dans les années 1980. Des chercheurs japonais ont introduit cette nouvelle famille de composites pour régler un problème majeur des composites traditionnels : le décollement des couches, ou délaminage . Ce phénomène, où les fibres se séparent de la matrice à cause de différences de dilatation thermique lors d'expositions à de hautes températures, limitait beaucoup leur utilisation dans ces conditions . Si l'on veut être plus précis, on attribue le concept à Kawasaki et Watanabe au milieu des années 80. Leur but initial était de réduire les contraintes thermiques dans des situations de haute température. Ils voulaient notamment résoudre le problème de la couche de protection thermique des navettes spatiales qui se décollait du substrat métallique, à cause de différences de température non uniformes entre les deux . Partie d'une motivation aérospatiale, cette idée a lancé un domaine de recherche qui a depuis beaucoup grandi.

Le principe même de conception des FGM est d'éliminer les interfaces nettes entre des matériaux différents. Elles sont remplacées par une transition continue d'une phase à une autre, suivant un gradient bien défini . Il en résulte une microstructure non homogène, permettant une variation progressive des propriétés des matériaux, sans interface de liaison spécifique . Cette architecture aide à réduire les concentrations de contraintes, à atténuer les gradients thermiques et à limiter les phénomènes de décollement ou de délaminage . La variation continue des propriétés sur l'épaisseur, en particulier, rend la structure plus résistante à la propagation des fissures et à la fatigue. Elle diminue aussi le risque de rupture dans les composants critiques .

Les applications des FGM n'ont cessé de se multiplier ces dix dernières années. C'est logique, vu la diversité des environnements sévères où ces matériaux peuvent être utilisés . L'aérospatiale, la défense, l'énergie nucléaire, l'électronique et le biomédical bénéficient directement de ces avancées . Dans l'aérospatiale et les transports, les FGM permettent d'alléger les structures tout en maintenant leur résistance mécanique. Résultat : on économise du carburant et on réduit les coûts d'exploitation . En biomédical, la capacité des FGM à imiter les caractéristiques mécaniques et thermiques des tissus naturels améliore la compatibilité et la fonctionnalité des dispositifs implantables, comme pour les implants dentaires . Cependant, ce survol introductif ne détaille pas toutes les méthodes numériques ni les équations qui régissent l'analyse thermomécanique quasi-tridimensionnelle ; ces aspects seront abordés dans les sections suivantes de ce travail.

Applications et Avantages des FGM

Les matériaux à gradient fonctionnel (FGM) fascinent de plus en plus la science et l'industrie. Pourquoi ? Les matériaux à gradient fonctionnel (FGM) suscitent un intérêt croissant dans les domaines scientifiques et industriels en raison de leurs profils de propriétés configurables pour des utilisations spécifiques, incluant la protection thermique, la résistance à l'usure ainsi que la résistance mécanique . Leur grand avantage ? Ils modulent continuellement composition et microstructure. Fini les interfaces nettes des composites classiques ! Ça réduit considérablement le risque de délaminage, ce décollement si problématique . En adoucissant les concentrations de contraintes et les variations de température, les FGM rendent les structures mécaniques plus efficaces et fiables . Moins de gradients de chaleur, moins de pics de contrainte : voilà qui améliore la résistance à la fatigue et à la propagation des fissures. Les composants durent plus longtemps, tout simplement .

Leur atout majeur ? Ils excellent dans les environnements chauds . Le Japon, dès les années 80, a imaginé ces matériaux pour contrer les contraintes thermiques extrêmes et éviter que les couches de protection ne se décrochent aux interfaces . Avec leur gestion fine des températures et leur pouvoir isolant, ils empêchent la dégradation et la rupture prématurée des matériaux . Cette régulation thermique progressive évite les catastrophes tout en préservant l'intégrité structurelle. Parfait pour les composants critiques des dispositifs de sécurité . En plus, la transition douce des propriétés allège les structures sans sacrifier la résistance mécanique. Économie de carburant et de coûts dans le transport : un vrai plus . Et ce n'est pas tout ! Les FGM réduisent les déchets de matière. Ils servent de pont entre des matériaux différents, sans recourir aux assemblages compliqués .

Concrètement, l'aérospatiale est un gros client. Pensez aux véhicules de rentrée atmosphérique, aux moteurs de turbines à gaz . L'industrie nucléaire s'en sert aussi, dans les réacteurs à fusion ou les centrales, notamment pour les tuyauteries soumises à de fortes chaleurs et à des contraintes thermomécaniques sévères . L'automobile les adopte pour des pièces comme les soupapes. Les secteurs chimique et minier y trouvent aussi leur compte . Et dans le domaine biomédical ? Les implants dentaires en FGM imitent les tissus naturels : meilleure compatibilité et fonctionnalité accrue . Les sources ne détaillent pas les formulations thermomécaniques quasi-3D précises pour chaque cas, mais elles confirment bien que les FGM sont utilisés partout dans ces branches de l'ingénierie .

Impact de la Porosité sur les Propriétés des FGM

Les structures poreuses à gradient fonctionnel, ou FGPSs, forment une catégorie d'éléments légers. Leurs propriétés ? Des gradients continus qui jouent sur la porosité et les caractéristiques matérielles, le tout de manière spatiale . L'introduction ou la présence inhérente de défauts de porosité modifie de manière substantielle le comportement mécanique global des matériaux à gradient fonctionnel (FGM) . En fait, la porosité a tendance à affaiblir la résistance globale du matériau . Concrètement, dans un composant, elle peut réduire la résistance totale, la rigidité, ou les deux simultanément . Ces changements sont si importants qu'on ne peut les ignorer. En raison de ces altérations notables, les effets de la porosité ne peuvent pas être négligés lors de la modélisation numérique et doivent impérativement être pris en compte dans l'étude du comportement mécanique des structures en FGM .

L'effet des défauts poreux sur les performances mécaniques va bien au-delà d'une simple dégradation. Ils provoquent des redistributions complexes des contraintes et des modifications dans les mouvements . Les défauts de porosité influencent fortement le déplacement du plan neutre de la structure et la manière dont les contraintes se répartissent globalement . Sur le plan des sollicitations de base, la porosité dans les FGM peut diminuer la résistance à la flexion, mais paradoxalement, augmenter la résistance à la compression . De surcroît, la présence de ces cavités peut générer des effets transversaux significatifs qu'il convient de caractériser avec précision .

La façon dont les pores se répartissent dans l'épaisseur ou le volume du matériau dicte largement les réponses structurales . Que la distribution des pores soit uniforme ou non, cela modifie considérablement les comportements mécaniques des structures FGM . Cette variabilité ouvre la porte à l'optimisation structurale : placer les pores stratégiquement peut améliorer les champs de contraintes . Par exemple, concentrer la porosité au centre d'une plaque, ou dans certaines couches, réduit les contraintes de traction normales le long de l'axe longitudinal . Ce mécanisme peut ainsi atténuer les points de forte contrainte, aidant l'intégrité de la structure sous charge .

Défauts de Fabrication et Modélisation de la Porosité

La porosité, l'un des défauts de fabrication les plus courants dans les Matériaux à Gradient Fonctionnel (FGM), influence grandement leur comportement thermomécanique. Lors de la création de plaques sandwiches à gradient fonctionnel (FGSP), des défauts de porosité peuvent apparaître, affaiblissant la rigidité de la structure et affectant directement ses réactions de vibration et de flambement . Ces imperfections sont généralement modélisées comme des critères de réduction de rigidité ; il est donc crucial de les prendre en compte pour prédire avec fiabilité le comportement structural . Le procédé de frittage, méthode fréquente pour fabriquer les FGM, favorise particulièrement la formation de micro-vides ou de porosités, ce qui en fait un défaut récurrent dans ces matériaux .

Modéliser des structures poreuses soulève plusieurs difficultés dues à leur complexité microstructurale. Les paramètres de fabrication qui contrôlent le taux de croissance et la taille des pores ne sont pas encore totalement compris, ni optimisés . Cette méconnaissance partielle des conditions de fabrication entraîne une grande variabilité dans les caractéristiques effectives des matériaux poreux. De plus, l'homogénéisation est essentielle pour modéliser ces structures, étant donné la complexité microstructurale élevée créée par les variations de taille des pores, d'épaisseur des parois, de forme et de distribution au sein d'un même matériau . Une étape d'homogénéisation, qu'elle soit locale ou globale, est donc indispensable avant toute analyse à l'échelle macroscopique. Cela permet de tenir compte de la différence d'ordre de grandeur entre la longueur caractéristique de la microstructure et les dimensions réelles de la pièce .

Les propriétés mécaniques effectives des matériaux poreux dépendent fortement de plusieurs facteurs, dont le procédé de fabrication, ainsi que la taille, la forme et la distribution des pores . Le module d'Young et le module de cisaillement sont particulièrement sensibles à ces paramètres. La prédiction analytique des propriétés des matériaux poreux devient alors complexe, en raison du caractère aléatoire de leurs structures . Il devient donc nécessaire d'avoir une connaissance aussi précise que possible de la microstructure pour obtenir une prédiction numérique significative . Dans ce contexte, plusieurs modèles micromécaniques ont été développés au fil des ans pour obtenir les propriétés effectives de matériaux composites ou poreux macroscopiquement homogènes. Cependant, les informations précises sur la taille, la forme et la distribution des particules ou des pores ne sont pas toujours disponibles .

Les travaux de Nguyen et al., cités dans la littérature, ont examiné le comportement mécanique de plaques poreuses à gradient fonctionnel. Ils ont considéré deux distributions de porosité distinctes variant dans l'épaisseur : les distributions dites paire et impaire . Par ailleurs, Zhang et Wang ont créé huit structures différentes de matériaux poreux avec diverses distributions, y compris des distributions à gradient. Ils les ont soumises à des essais mécaniques pour évaluer des propriétés importantes comme le module d'Young . Ces résultats expérimentaux ont révélé une relation évidente entre cette propriété et la porosité, bien que cette relation ne soit pas forcément linéaire. Cela offre une base solide pour valider les propriétés obtenues par des modèles théoriques . Les sources consultées ne précisent cependant pas toutes les distributions de porosité envisageables, ni les seuils critiques au-delà desquels la dégradation des propriétés mécaniques devient irréversible pour les FGM.

Théories Classiques de Déformation des Plaques

L'analyse des structures en plaques s'est toujours appuyée sur des hypothèses cinématiques simplifiées. Le but ? Éviter la complexité des modèles d'élasticité tridimensionnelle. La plus basique, c'est la théorie classique des plaques, aussi appelée théorie de Kirchhoff-Love . Elle ignore complètement les déformations dues au cisaillement transversal . Résultat : elle ne s'applique qu'aux plaques minces .

Très vite, on a cherché des modèles plus fidèles. Les théories affinées sont apparues dès le dix-neuvième siècle . Pour la dynamique, Bresse et Lord Rayleigh ont ouvert la voie . En parallèle, pour la statique, des pionniers comme Poncelet, Rankine, Grashof, Engesser et Föppl ont développé des théories plus poussées . Ces travaux ont permis des descriptions mécaniques plus précises, corrigeant les lacunes des hypothèses classiques .

Pourtant, malgré ces avancées, les théories classiques sont encore couramment utilisées dans les logiciels de calcul par éléments finis . Le problème ? Elles ne modélisent pas toujours bien les propriétés mécaniques transversales . Ni les conditions aux limites sur les surfaces externes . Ces limites apparaissent surtout avec les plaques épaisses, mais aussi dès qu'on analyse des effets de bord, des déformations locales ou des interactions complexes comme les problèmes thermiques .

Théories de Déformation au Cisaillement d'Ordre Supérieur (FSDT)

La théorie de la déformation au cisaillement du premier ordre (FSDT) est née pour pallier les lacunes de la théorie classique des plaques (CPT), aussi appelée théorie de Kirchhoff-Love. Elle intègre en effet les déformations dues au cisaillement transversal. Or, la CPT part du principe que les sections planes restent planes et perpendiculaires à la surface moyenne. Dans le cadre de la CPT, les sections planes restent planes et normales à la surface moyenne, ce qui revient à négliger les déformations de cisaillement transverse et limite par conséquent son applicabilité aux plaques minces . La FSDT, elle, relâche cette hypothèse : les sections planes restent planes après déformation, mais plus forcément perpendiculaires à la surface moyenne. Elle permet ainsi de prendre en compte les effets du cisaillement sur toute l'épaisseur de la plaque . Cette théorie constitue ainsi une extension de la CPT et représente l'approche la plus simple pour analyser les stratifiés composites minces à modérément épais .

Historiquement, les équations dynamiques des poutres, qui intègrent la déformation au cisaillement et l'inertie de rotation, portent la signature de Timoshenko et Ehrenfest. On les appelle les équations de Timoshenko-Ehrenfest, ou encore la théorie de la déformation au cisaillement du premier ordre (FSDT) . Pour les plaques, des théories plus sophistiquées ont vu le jour grâce à Basset, Reissner, Hencky, Bolle et Hildebrand et al. Quant aux équations dynamiques des plaques déformables en cisaillement avec inertie de rotation, elles ont d'abord été établies par Uflyand en utilisant la méthode de l'équilibre dynamique. Puis, Mindlin les a redérivées en s'appuyant sur le principe de Hamilton . Ces équations, aujourd'hui connues sous le nom d'équations de Uflyand-Mindlin, sont aussi qualifiées de FSDT .

La FSDT a cependant un défaut majeur. Son hypothèse d'une variation linéaire des déplacements sur l'épaisseur conduit à des contraintes et des déformations de cisaillement transversales constantes sur toute cette épaisseur . Cette hypothèse va à l'encontre de deux conditions fondamentales : d'une part, les contraintes nulles sur les surfaces supérieure et inférieure de la plaque, et d'autre part, la distribution parabolique des contraintes de cisaillement transversales sur l'épaisseur . L'intégrale des contraintes de cisaillement transversales calculée avec la FSDT est donc supérieure aux valeurs réelles. Il faut donc introduire un facteur de correction de cisaillement, généralement inférieur à un, pour prédire correctement le comportement mécanique réel des plaques et obtenir des résultats analytiques acceptables pour les plaques moyennement épaisses et minces . Choisir ce facteur de correction n'est pas simple, même si K = 5/6 est une valeur constante souvent utilisée dans les simulations numériques. Timoshenko, lui, avait proposé un coefficient dépendant du coefficient de Poisson .

En fin de compte, la FSDT s'avère assez précise pour les plaques minces à moyennement épaisses. Mais ses limites ont poussé au développement de la théorie de la déformation au cisaillement d'ordre supérieur (HSDT). Cette dernière offre un outil plus puissant et plus précis pour analyser les composites stratifiés, capable de résoudre les problèmes que la FSDT rencontre .

Théories de Déformation au Cisaillement d'Ordre Supérieur (HSDT)

Dans le cadre de l'analyse structurelle des plaques et des stratifiés composites, les théories cinématiques ont progressivement évolué afin de surmonter les limites inhérentes aux approches conventionnelles . La théorie classique des plaques, désignée par CPT ou CLPT, convient uniquement aux plaques composites minces car elle néglige les effets du cisaillement transversal . Pour pallier ce manque, la théorie de déformation au cisaillement du premier ordre, ou FSDT, propose une formulation simple adaptée aux stratifiés composites minces à modérément épais . Néanmoins, l'hypothèse d'une variation linéaire des déplacements dans l'épaisseur adoptée par la FSDT produit des déformations et des contraintes de cisaillement transversal constantes . Cette cinématique viole les conditions aux limites de contraintes nulles sur les surfaces supérieure et inférieure, imposant ainsi le recours à un facteur de correction de cisaillement pour ajuster les résultats .

Afin de s'affranchir de cette contrainte, les théories de déformation au cisaillement d'ordre supérieur, communément abrégées en HSDT, ont été introduites . Ces modèles formulent les composantes de déplacement à l'aide de séries entières suivant les coordonnées de l'épaisseur . En enrichissant la cinématique et en augmentant le nombre de variables du champ de déplacement, les HSDT permettent au cisaillement transversal de varier de façon non linéaire à travers l'épaisseur de la structure . Cette modélisation permet de satisfaire rigoureusement les conditions de nullité des contraintes sur les surfaces supérieure et inférieure de la plaque, éliminant par conséquent le besoin d'introduire un facteur de correction de cisaillement artificiel . Les théories HSDT s'imposent ainsi comme des outils particulièrement puissants et précis pour l'analyse des stratifiés composites .

Parmi les diverses variantes développées au sein de cette catégorie, le modèle proposé par Reddy en 1984, fondé sur la théorie de déformation au cisaillement du troisième ordre, ou TSDT, représente une approche largement exploitée . Ce modèle particulier est couramment retenu en raison de sa simplicité et de son niveau de précision élevé dans l'évaluation du comportement mécanique . L'approximation cinématique de la théorie d'ordre trois permet une représentation satisfaisante des déformations sans introduire de complexité excessive dans les formulations régissantes .

L'accroissement du nombre de variables cinématiques inhérent aux approches d'ordre supérieur améliore la précision globale mais induit simultanément une augmentation du temps de calcul . Il a toutefois été observé que la précision tend à converger à mesure que l'ordre des développements polynomiaux augmente . En ce qui concerne l'évaluation du déplacement transversal, les théories de quatrième ordre ne démontrent d'ailleurs qu'une amélioration marginale par rapport aux théories de troisième ordre .

Théories Quasi-3D et Modèles Raffinés

L'analyse thermomécanique des structures en matériaux à gradient fonctionnel (FGM) exige des représentations cinématiques qui rendent compte des déformations de cisaillement et des variations dimensionnelles en épaisseur. Les théories dites quasi-tridimensionnelles (quasi-3D) incorporent explicitement l'effet d'étirement dû à l'épaisseur . Contrairement aux approches classiques qui ignorent ce facteur, les modèles quasi-3D et les théories affinées de cisaillement d'ordre élevé (HSDT) garantissent une distribution parabolique des contraintes et déformations de cisaillement transversal à travers l'épaisseur . Cette modélisation permet de satisfaire rigoureusement les conditions aux limites de traction nulle sur les surfaces libres supérieure et inférieure des plaques sans nécessiter l'introduction d'un facteur de correction de cisaillement .

Une caractéristique majeure de ces formulations réside dans la structuration du champ de déplacement au moyen de fonctions de forme spécifiques et de variables intégrales indéterminées . L'utilisation d'intégrales indéterminées et de fonctions de forme trigonométriques, hyperboliques ou du troisième ordre permet de décrire fidèlement le gauchissement de la section transversale . Les théories affinées réussissent ainsi à réduire significativement le nombre de degrés de liberté requis, tout en maintenant une précision comparable à celle des HSDT conventionnelles qui emploient plus de variables . Certains modèles limitent même le champ cinématique à quatre inconnues seulement, ce qui est inférieur aux théories de cisaillement d'ordre premier (FSDT) . Dans le cadre d'une formulation à quatre variables pour l'analyse thermomécanique de plaques sandwichs FGM sur une fondation élastique, le système complet ne comprend que quatre équations directrices . D'autres méthodes, comme la théorie de cisaillement trigonométrique à cinq inconnues (FTSD), sont aussi proposées pour intégrer l'étirement transversal dans l'étude du flambement et de la flexion de plaques poreuses .

L'établissement des équations d'équilibre ou de mouvement et des conditions aux limites associées s'appuie généralement sur le principe des déplacements virtuels . Ces relations variationnelles permettent d'étudier des géométries variées, allant des plaques sandwichs FGM jusqu'aux poutres fissurées à gradient de propriétés fonctionnelles . De nombreux modèles étendent cette analyse aux structures FGM reposant sur des fondations élastiques modélisées par des paramètres de type Winkler, Pasternak ou selon le modèle à trois paramètres de Kerr .

Pour obtenir les solutions sous chargement, des procédures analytiques telles que la méthode de Navier sont appliquées pour des plaques FGM simplement appuyées présentant différentes distributions de porosité . Slimani, Menasria, Rachedi, Mourad, Refrafi, Nimer, Bouhadra et Mamen ont notamment formulé une théorie quasi-3D d'ordre élevé qui réduit le nombre d'inconnues pour la flexion statique de plaques poreuses FGM, résolue par la méthode de Navier . En parallèle des formulations analytiques, des approches numériques telles que la méthode sans maillage (*meshless method*) sont utilisées pour traiter la théorie quasi-3D de cisaillement d'ordre élevé appliquée aux plaques FGM .

Méthodes Numériques et Analytiques pour l'Analyse Thermomécanique

L'analyse thermomécanique quasi-tridimensionnelle des plaques et poutres en FGM mobilise un éventail de méthodes numériques et analytiques, dont le choix conditionne la précision et le coût de calcul des résultats. Parmi les techniques sans maillage, la méthode Galerkin sans élément (element-free Galerkin) a été appliquée à l'analyse thermomécanique de plaques FGM . Pour les plaques sandwich à cœur FGM, la méthode de quadrature différentielle généralisée (generalized differential quadrature method) a été employée afin d'obtenir la solution thermoélastique . Dans le cadre de la stabilité de microplaques FGM minces soumises à des chargements mécaniques et thermiques, la méthode des bandes finies par splines (spline finite strip method), combinée à la théorie du couple de contrainte modifié, a été utilisée .

Du côté des formulations par éléments finis, Sahoo, Verma et Singh ont proposé en 2024 un modèle éléments finis C0 à formulation pénalisée pour une théorie de déformation par cisaillement d'ordre supérieur (HSDT) à cinq variables, dédié à l'analyse de flexion linéaire de plaques FGM sous charges mécaniques et thermiques, en s'appuyant sur la règle des mélanges et sur sept degrés de liberté par nœud, et en comparant la théorie de déformation par cisaillement d'ordre trois (TSDT) à une théorie non polynomiale (NPSDT) . L'analyse isogéométrique (IGA) fondée sur la TSDT et utilisant les fonctions de base NURBS a été étendue en 2023 à l'analyse vibratoire de plaques FGP en environnement thermique reposant sur une fondation de Kerr, le champ de température à travers l'épaisseur satisfaisant l'équation de transfert thermique en régime permanent et les équations gouvernantes étant dérivées du principe de Hamilton . Plus récemment, l'approche éléments finis de continuation d'ordre élevé (FE-HOCA), reposant sur la HOSDT et la méthode asymptotique numérique, a été mobilisée pour étudier le flambement et le post-flambement non linéaires de plaques FGM poreuses sur fondations élastiques; cette méthode utilise un pas adaptatif pour s'ajuster aux non-linéarités locales et reformule les équations d'équilibre par un développement en série de Taylor en une séquence récursive de sous-problèmes linéaires, avec un élément quadrilatéral à huit nœuds comportant neuf degrés de liberté par nœud .

Sur le plan analytique, les équations gouvernantes des plaques FGM soumises à des chargements mécaniques sont dérivées par l'application du principe des travaux virtuels, en substituant l'énergie de déformation virtuelle, l'énergie potentielle virtuelle des charges appliquées et l'énergie de déformation virtuelle de la fondation élastique dans ce principe . Les équations de mouvement sont quant à elles établies à partir du principe de Hamilton, la variation de l'énergie de déformation totale étant représentée par une intégrale sur le volume . Pour les poutres FGM, les modèles de fondation Winkler, Pasternak et Kerr sont incorporés dans le Hamiltonien afin de dériver les équations du mouvement . Enfin, la solution de Navier pour une plaque rectangulaire simplement appuyée sur ses bords est utilisée pour résoudre analytiquement les équations gouvernantes dérivées du principe de Hamilton, comme dans le cas de l'analyse en vibration libre de plaques sandwich à cœur auxétique et faces FGM poreuses .

Les sources consultées ne précisent pas de comparaison systématique de coût de calcul ni de critères quantitatifs de convergence entre ces méthodes; elles ne discutent pas non plus de benchmarks unifiés permettant de hiérarchiser leur précision respective. Ces aspects restent donc à documenter dans la littérature examinée.

Modélisation des Propriétés Matérielles et des Chargements

La modélisation thermomécanique quasi-3D des matériaux à gradient fonctionnel (FGM) exige une description précise de la façon dont leurs propriétés varient dans l'espace et des contraintes qu'ils subissent. Pour les plaques poreuses en FGM, on calcule souvent leurs propriétés effectives, en tenant compte de la porosité, grâce à une loi de puissance modifiée . Cette approche permet de bien décrire la transition progressive entre les composants, typiquement une céramique et un métal, dont le volume fractionnaire change sur l'épaisseur selon un exposant spécifique . On combine souvent des céramiques, peu conductrices de chaleur, avec des métaux, qui apportent ductilité et résistance à la fracture, en se basant sur des distributions suivant une loi de puissance ou exponentielle sur l'épaisseur .

En environnement thermomécanique, les équations des plaques FGM doivent intégrer la dépendance des propriétés du matériau par rapport à la température . Cette dépendance s'exprime généralement en considérant chaque paramètre effectif – comme le module d'Young, le coefficient de Poisson, la conductivité thermique ou le coefficient de dilatation thermique – comme une somme pondérée des propriétés du métal et de la céramique, elles-mêmes fonctions de la température . Une telle formulation capture mieux le comportement des plaques FGM lorsqu'elles sont exposées à des montées en température importantes .

Pour modéliser les charges thermiques, on suppose souvent une conduction thermique unidimensionnelle, stable, à travers l'épaisseur. La température varie alors de façon non linéaire entre la surface supérieure et inférieure . On résout cette distribution par des approximations polynomiales et on impose des conditions de continuité thermique aux interfaces . Les contraintes mécaniques incluent des compressions uniaxiales ou biaxiales, des combinaisons de compression et traction, ainsi que des chargements thermiques critiques. Leur analyse de stabilité nécessite de prendre en compte les imperfections géométriques et les non-linéarités géométriques .

Dans le cas des micro-plaques poreuses axisymétriques, la formulation quasi-3D de Reddy, d'ordre trois, a été utilisée avec une distribution en loi de puissance à deux composants et trois distributions de porosité définies par des fonctions cosinus . La non-linéarité géométrique et la théorie modifiée du couple de contrainte, qui introduit un paramètre unique de longueur caractéristique, complètent cette modélisation . Pour les plaques poreuses saturées, la théorie de Reddy d'ordre supérieur s'associe à la relation cinématique non linéaire de von Kármán et à la théorie élastique de Biot pour matériaux poreux saturés. La résolution s'effectue analytiquement par la méthode de Galerkin et la fonction de contrainte d'Airy pour des appuis simples .

Il faut cependant noter plusieurs limites à ces approches. L'hypothèse d'une loi de puissance pour décrire la composition spatiale des FGM peut être trop restrictive pour certaines variations de propriétés . De plus, les capacités de calcul limitées par le passé restreignaient le nombre de variables de conception qu'on pouvait considérer lors de la définition de la forme fonctionnelle de la composition matérielle . Prédire analytiquement les propriétés des matériaux poreux reste complexe, car leurs structures sont aléatoires et nécessitent une connaissance précise de la microstructure . Enfin, la relation entre les propriétés du matériau et la porosité n'est pas forcément linéaire, ce qui rend l'élaboration de modèles prédictifs fiables plus difficile . Les modèles d'ordre supérieur devraient mieux prédire le comportement des plaques épaisses avec divers modèles de porosité . Les sources consultées ne précisent pas toutes les distributions de chargement thermique et mécanique utilisées récemment, ni comment les conditions aux limites affectent la modélisation des propriétés effectives.

Analyse des Comportements Spécifiques (Vibrations, Flambage, Stabilité)

Les structures en FGM se comportent de manière mécanique très particulière. Pour les étudier, il faut des modèles qui tiennent compte de la nature hétérogène de ce matériau. Le délaminage, c'est-à-dire la séparation entre les couches, est une cause d'endommagement fréquente dans les composites stratifiés. Il peut mener à une rupture prématurée, réduire la capacité portante et dégrader les performances des composants à base de FGM. Il est donc crucial d'en examiner les mécanismes en détail . L'analyse de ces comportements fait appel à des théories de déformation en cisaillement d'ordre supérieur et à des méthodes numériques capables de saisir les non-linéarités géométriques et matérielles.

Pour ce qui est des vibrations non linéaires, une étude a porté sur des plaques FGM poreuses, saturées et soumises à un environnement thermique. Elle a utilisé la théorie de déformation en cisaillement d'ordre supérieur de Reddy, couplée à la relation de déformation non linéaire géométrique de von Kármán et à la théorie élastique de Biot pour les matériaux poreux saturés. Une solution analytique de Galerkin associée à une fonction de contrainte d'Airy a été employée pour les plaques simplement appuyées, en considérant différentes distributions de porosité : uniforme, symétriquement irrégulière et asymétriquement irrégulière . Cette même recherche a révélé plusieurs points importants : l'augmentation du coefficient de porosité fait baisser la fréquence propre et monte l'amplitude de déflexion ainsi que le rapport de fréquence de la plaque poreuse saturée. Le coefficient de Skempton, bien que mineur, a un effet positif sur la rigidité. Les fondations élastiques augmentent la fréquence propre. Enfin, l'accroissement de la variation de température réduit significativement la fréquence propre tout en augmentant l'amplitude de déflexion .

L'analyse du flambage et du post-flambage des plaques FGM poreuses s'appuie sur des approches numériques avancées. Une étude s'est intéressée aux caractéristiques de stabilité et de post-flambage de telles plaques reposant sur une fondation élastique non linéaire géométrique de type Winkler/Pasternak. Pour ce faire, elle a utilisé une approche par éléments finis de continuation d'ordre supérieur (FE-HOCA) intégrant la théorie de déformation en cisaillement d'ordre supérieur (HOSDT) . Cette méthode se distingue par son pas adaptatif, qui ajuste la résolution aux non-linéarités locales de la solution. Elle facilite ainsi l'identification des points de bifurcation. De plus, elle reformule les équations d'équilibre en une série récursive de sous-problèmes linéaires grâce à un développement en série de Taylor .

L'impact de la porosité sur le comportement mécanique des plaques FGM a également été documenté. L'introduction d'un paramètre de porosité augmente la déflexion et les contraintes de cisaillement transverse. Elle diminue par contre la contrainte axiale dans les plaques FGM poreuses . On observe des déplacements plus importants pour les plaques FGM poreuses que pour les plaques non poreuses. La distribution uniforme de porosité accentue ces déplacements pour une même fraction volumique de porosité . Par ailleurs, une distribution de porosité non uniforme crée un module de Young minimal au plan médian de la plaque et maximal aux surfaces externes, car sa valeur dépend du signe de la coordonnée z .

L'analyse des vibrations libres de structures sandwich, combinant des noyaux auxétiques et des faces FGM poreuses, a été traitée via une solution analytique. Celle-ci repose sur la théorie de déformation en cisaillement d'ordre supérieur à quatre variables (HSDT-4). Les équations gouvernantes ont été dérivées du principe de Hamilton et résolues par la solution de Navier pour une plaque rectangulaire simplement appuyée. Cette plaque reposait sur des fondations élastiques de type Winkler, Pasternak ou Kerr . Les résultats ont été validés par comparaison avec des études antérieures. Les effets des propriétés matérielles, des paramètres géométriques du noyau auxétique et de la fondation élastique sur les caractéristiques vibratoires ont été analysés . Les sources consultées n'indiquent pas tous les paramètres quantitatifs qui régissent ces comportements spécifiques, ni l'étendue complète des configurations de chargement thermique et mécanique étudiées dans la littérature dite quasi-3D.

Limitations et Perspectives Futures

Les matériaux à gradient fonctionnel (FGM) ont des atouts structurels indéniables, comme cette réduction des contraintes concentrées et un meilleur contrôle des gradients thermiques . Ça améliore l'efficacité, la fiabilité, la résistance à la fatigue et même la propagation des fissures, bref, ça fait durer les composants et les structures plus longtemps . Mais voilà, analyser et modéliser thermomécaniquement ces composites, c'est un vrai casse-tête, avec des blocages méthodologiques et physiques qui limitent les approches actuelles.

Un problème majeur ? La délamination. C'est quand les couches se séparent, un mode de rupture classique pour les composites stratifiés . Ça peut causer une rupture rapide, une perte de capacité de charge, et globalement, de mauvaises performances pour les FGM. Il faut donc vraiment creuser ces mécanismes de délamination . Sur le plan de la modélisation cinématique, les théories classiques, celles qu'on trouve souvent dans les logiciels d'éléments finis, montrent leurs limites . Elles peinent à modéliser les propriétés mécaniques transverses et à respecter scrupuleusement les conditions aux limites sur les surfaces externes . C'est encore plus vrai pour les plaques épaisses, les effets de bord, les distorsions locales et les interactions complexes, comme celles dues à la chaleur .

À ces soucis de formulation s'ajoutent des contraintes historiques et actuelles liées à la puissance de calcul et à la conception géométrique . Avant, avec des ordinateurs moins performants, on ne pouvait pas prendre en compte beaucoup de variables pour définir la composition du matériau . Même si l'optimisation numérique progresse, la fabrication reste un défi de taille . En particulier, on ne maîtrise pas encore parfaitement les paramètres du procédé pour contrôler la croissance et la taille des pores .

La présence de pores dans ces structures complique sacrément la modélisation micromécanique . L'étape d'homogénéisation est primordiale pour les structures poreuses, vu la complexité microstructurale due aux variations de taille, d'épaisseur de paroi, de forme et de distribution des pores . Et puis, il faut bien se servir de modèles théoriques pour concevoir les FGM, car il y a une diversité folle de compositions et de distributions possibles pour arriver aux propriétés spécifiques qu'on vise .

Pour éviter les coûts de calcul énormes des formulations 3D complètes ou des théories d'ordre élevé, développer des formulations raffinées avec moins de variables est une piste prometteuse . On a montré qu'une théorie raffinée à quatre variables demande moins d'efforts de calcul que d'autres méthodes analytiques ou numériques basées sur des théories de plaques d'ordre trois ou des théories raffinées à cinq variables . Les recherches futures devront donc combiner ces réductions cinématiques efficaces avec une description rigoureuse de l'homogénéisation microstructurale et des conditions aux limites thermomécaniques complexes.

Questions ouvertes

— Les paramètres de traitement pour contrôler la vitesse et la taille de croissance des pores sont-ils entièrement compris et optimisés ?
— La relation entre la propriété du matériau et la porosité est-elle toujours non linéaire ?
— Les modèles de plus haut niveau se comportent-ils mieux pour prédire le comportement des structures de plaques plus épaisses avec différents modèles de porosité ?