Calcul thermomécanique : quand la température devient un effort

Ressource AVNIR Engineering

Publié le 4 août 2026, à jour des référentiels et bonnes pratiques en vigueur.

Actualité — Le 29 avril 2026, l’ESA a annoncé que le système de protection thermique de Space Rider, ventre et volets de pilotage, avait été porté à 1 600 °C dans la soufflerie à plasma du CIRA à Capoue, sous un jet d’arc à dix fois la vitesse du son. Ces volets de 10 kg, 90 × 70 cm, doivent piloter un module de 3 000 kg en rentrée : tenir la chaleur et transmettre l’effort en même temps, c’est la définition du calcul thermomécanique (ESA, avril 2026).

Une pièce chauffée qui peut se dilater librement ne développe aucune contrainte. C’est la dilatation empêchée qui casse, et cette seule phrase réorganise toute la compréhension du calcul thermomécanique. En aérospatial, la température n’est presque jamais un simple paramètre de confort : elle est un chargement à part entière, capable de plastifier une pièce sans qu’aucun effort extérieur ne s’applique. Cet article reprend les trois modes de transfert et ce qui domine selon le cas, la différence entre régime établi et transitoire, le mécanisme exact par lequel un champ de température devient une contrainte, et les trois situations où le sujet se pose vraiment : équipement embarqué, protection thermique, assemblage de matériaux qui ne se dilatent pas au même rythme.

Qu’est-ce qu’un calcul thermomécanique

Un calcul thermomécanique est un calcul de structure dont le chargement principal n’est pas une force extérieure, mais un champ de température. On résout d’abord la thermique, ce qui donne la température en chaque point et à chaque instant, puis on injecte ce champ dans le modèle mécanique, qui rend des déplacements et des contraintes. Deux physiques, un maillage commun le plus souvent, et une question qui n’a rien d’académique : la pièce tient-elle alors qu’aucun effort extérieur ne lui est appliqué ?

La réponse est fréquemment non. Un support d’équipement peut être largement dimensionné vis-à-vis des efforts mécaniques et se fissurer sur un simple cycle de chauffage. La cause n’est pas la température, mais la manière dont la structure environnante empêche la matière de bouger.

Conduction, convection, rayonnement : ce qui domine selon le cas

Le transfert thermique repose sur trois mécanismes, qui coexistent presque toujours mais dont un seul, en pratique, pilote le résultat. L’identifier avant de lancer le moindre calcul évite de raffiner un maillage là où l’incertitude n’est pas.

• La conduction transporte la chaleur dans la matière, proportionnellement au gradient de température et à la conductivité. L’écart entre matériaux est considérable : 120 W/m·K pour un alliage d’aluminium, 15 W/m·K pour un acier inoxydable austénitique, 7 W/m·K pour un alliage de titane, 0,05 W/m·K pour un feutre en fibres céramiques.
• La convection échange entre une paroi et un fluide en mouvement, à travers un coefficient exprimé en W/m²·K : quelques dizaines en convection naturelle dans l’air, plusieurs centaines en convection forcée. C’est le terme le plus incertain du modèle, et la charnière avec la mécanique des fluides : dès que l’écoulement conditionne le résultat, il faut le calculer, ce qui relève de la simulation d’écoulement.
• Le rayonnement thermique ne demande aucun support matériel et varie comme la puissance quatrième de la température absolue. À 1 600 °C, un corps noir émet de l’ordre de 700 kilowatts par mètre carré : au-delà de quelques centaines de degrés, il écrase les deux autres termes.

Le vide spatial, lui, supprime la convection. Un équipement embarqué sur satellite n’évacue sa puissance que par conduction vers sa structure porteuse et par rayonnement vers l’espace, architecture entièrement différente de celle du même boîtier en soute d’avion.

Régime établi et régime transitoire : pourquoi le transitoire dimensionne

Le régime établi correspond à l’état d’équilibre atteint quand les flux entrants et sortants se compensent : la température ne varie plus dans le temps. Le régime transitoire décrit le chemin pour y parvenir, ou l’évolution permanente quand les conditions changent plus vite que la structure ne peut suivre.

Le réflexe naturel consiste à calculer l’état final, plus rapide à poser. C’est presque toujours le mauvais choix. Pendant le transitoire, la peau chauffe avant le cœur, le gradient dans l’épaisseur passe par un maximum, et c’est ce gradient qui crée la contrainte. Une paroi de 10 mm dont la face exposée monte de 500 °C en quelques secondes voit ses deux peaux se dilater très différemment alors que rien n’a bougé au centre.

La diffusivité thermique, rapport de la conductivité à la capacité calorifique volumique, donne le premier repère de cinétique : 5 × 10⁻⁵ m²/s pour un alliage d’aluminium contre 4 × 10⁻⁶ m²/s pour un acier inoxydable. Le nombre de Biot dit ensuite si la pièce peut être traitée comme isotherme ou si son gradient interne doit être modélisé pour lui-même.

Ce qui crée vraiment la contrainte : la dilatation empêchée

Une contrainte d’origine thermique est le produit d’une dilatation empêchée, jamais d’une dilatation en soi. Le point est décisif et souvent manqué. Une barre chauffée uniformément, libre à ses deux extrémités, s’allonge selon son coefficient de dilatation thermique et ne développe strictement aucune contrainte : elle change de dimension, c’est tout. Placez la même barre entre deux appuis rigides, ne changez rien d’autre, et elle développe une contrainte de compression égale au produit du module d’Young, du coefficient de dilatation et de l’écart de température.

Le calcul est éloquent. Pour un alliage d’aluminium, module d’Young voisin de 70 GPa et coefficient de dilatation de 23 µm/m·K, un échauffement uniforme de 100 °C sur une pièce parfaitement bridée produit environ 161 MPa, soit près de la moitié de la limite d’élasticité d’un 2024-T3, de l’ordre de 345 MPa. Sans le moindre effort extérieur. Pour un acier, 210 GPa et 12 µm/m·K, le même échauffement donne 252 MPa : l’acier se dilate deux fois moins que l’aluminium, mais sa raideur trois fois supérieure fait plus que compenser.

Le bridage prend trois formes, qu’il faut savoir reconnaître sur une pièce réelle :

• un bridage externe, imposé par les conditions aux limites : une pièce boulonnée sur une structure restée froide, une bride encastrée, un montage hyperstatique ;
• un bridage interne, dû au gradient : les fibres chaudes voudraient s’allonger davantage que les fibres froides du même solide, s’en empêchent mutuellement, et la pièce s’auto-contraint sans liaison extérieure ;
• un bridage d’assemblage, quand deux matériaux de coefficients différents sont liés et se contraignent mutuellement.

D’où une conséquence de conception immédiate : la parade n’est pas toujours de refroidir, elle consiste très souvent à libérer. Un trou oblong, un appui glissant ou un point fixe unique suppriment la contrainte en autorisant le déplacement.

Assembler deux matériaux à coefficients de dilatation différents

L’assemblage de deux matériaux à coefficients de dilatation différents est le cas thermomécanique le plus fréquent et le plus sous-estimé, parce qu’il produit des contraintes sans qu’aucune source de chaleur spectaculaire ne signale le problème. Les ordres de grandeur usuels suffisent à voir venir la difficulté :

Une ferrure en alliage d’aluminium boulonnée sur un panneau en composite carbone illustre le mécanisme. L’écart de coefficient atteint environ 22 µm/m·K : sur une interface de 500 mm et un écart de 100 °C, l’allongement différentiel vaut 1,1 mm. Cette différence ne disparaît pas, elle se retrouve intégralement en cisaillement dans les fixations et en pelage au bord d’un collage. Sur une structure exposée alternativement au soleil et à l’ombre, l’amplitude encaissée dépasse la centaine de degrés, et elle se répète.

Les parades relèvent du choix matériau autant que de l’architecture : intercaler un alliage de titane à coefficient intermédiaire, scinder une grande interface en zones courtes, adapter la séquence d’empilement du stratifié. Le choix des alliages métalliques employés en aéronautique ne se joue donc jamais sur la seule tenue mécanique.

Fatigue thermique et cyclage : l’usure par la température

La fatigue thermique désigne l’endommagement progressif d’une pièce soumise à des variations de température répétées, chaque cycle thermique se comportant comme un cycle de chargement mécanique. Deux particularités la distinguent : les amplitudes de déformation sont souvent élevées, ce qui place le problème dans le domaine oligocyclique, et les propriétés du matériau évoluent au cours du cycle.

Les points de concentration sont prévisibles : bords de collage, pieds de cordon de soudure, brasures d’un composant sur son circuit imprimé, congés entre zone épaisse et zone mince. Les profils de qualification d’un équipement embarqué balayent couramment une plage de -55 °C à +85 °C, répétée sur des dizaines de cycles. Ces essais d’environnement, encadrés par les référentiels DO-160, MIL-STD-810 et ECSS, s’articulent avec les campagnes de qualification vibratoire : une même pièce doit tenir les deux.

Le CIRA l’a illustré en février 2026 sur la protection thermique de Space Rider : un assemblage de quatre panneaux dépassant le mètre carré, isolant jusqu’à 1 600 °C, a d’abord dû franchir une qualification vibratoire sur un pot vibrant de 200 kN. Une structure thermique doit survivre au lancement avant d’affronter la rentrée.

Trois cas d’usage du calcul thermomécanique en aérospatial

Le calcul thermomécanique sert dans trois familles de situations, dont les difficultés diffèrent nettement.

L’équipement embarqué et sa dissipation. Un boîtier d’électronique dissipant quelques dizaines de watts doit évacuer cette puissance sans dépasser la température de jonction admissible de ses composants, souvent 125 °C, tout en tenant les contraintes induites dans ses cartes et ses fixations. La résistance thermique de contact aux interfaces pèse fréquemment plus lourd dans le bilan que la conductivité des matériaux eux-mêmes.

La protection thermique. Il s’agit d’empêcher un flux considérable d’atteindre une structure qui ne le supporterait pas, sur une durée limitée. Le dimensionnement se joue sur le transitoire : épaisseur d’isolant, tenue de l’accrochage, fixations traversantes qui font pont thermique. La NASA a publié en avril 2026 ses premières évaluations après le vol d’Artemis II. Le bouclier ablatif d’Orion, entré dans l’atmosphère à près de Mach 35, s’est comporté comme prévu, avec une perte de matériau carbonisé nettement réduite par rapport au vol précédent et une performance cohérente avec les essais au sol en soufflerie à arc. Le calcul se valide contre l’essai au sol, l’essai au sol se valide contre le vol.

L’assemblage hétérogène. Traité plus haut, c’est le cas le plus banal et celui qui produit le plus de retours d’expérience désagréables : aucun événement thermique remarquable ne le signale, et personne ne pense à ouvrir une analyse thermique pour une ferrure.

Chaîner le calcul thermique et le calcul de structure

Le chaînage thermomécanique consiste à transférer le champ de température issu du solveur thermique vers le modèle mécanique, qui l’utilise comme déformation imposée. Le couplage faible résout les deux physiques l’une après l’autre : il suffit tant que la déformation mécanique ne modifie pas le problème thermique. Le couplage fort les résout simultanément, et devient nécessaire quand la déformation change la thermique elle-même, par exemple lorsqu’un jeu s’ouvre à une interface et fait varier la résistance de contact.

Trois précautions gouvernent la qualité du résultat. Le transfert de champ entre deux maillages différents doit conserver l’énergie et ne pas lisser les gradients, faute de quoi le pic de contrainte disparaît en route. Les propriétés introduites doivent être celles de la température atteinte. Enfin l’état initial doit correspondre à un état sans contrainte identifié, température de montage ou de cuisson d’un collage, sans quoi tout le champ est décalé d’un bloc. La démarche, les hypothèses à consigner et le format du livrable rejoignent ceux de tout calcul de structure justifié, note de calcul comprise.

Ce qu’un calcul thermique ne dit pas

Un calcul thermique repose sur des conditions aux limites dont l’incertitude dépasse presque toujours celle du modèle lui-même : c’est sa première limite, et elle est structurelle. Un coefficient d’échange convectif est rarement connu à mieux qu’un facteur deux. Une émissivité dépend de l’état de surface, du revêtement et de son vieillissement : elle passe de 0,1 pour un métal poli à plus de 0,9 pour une peinture noire mate. Une résistance de contact entre deux pièces boulonnées dépend de la pression de serrage et de la rugosité.

Deuxième limite, la dépendance des propriétés à la température. Conductivité, coefficient de dilatation, module d’Young et limite d’élasticité évoluent tous avec la température : le module d’un alliage d’aluminium chute sensiblement au-delà de 150 °C, sa limite d’élasticité davantage encore. Entrer des propriétés relevées à 20 °C dans un calcul à 300 °C produit un résultat propre et faux. Cette matière vient de la caractérisation expérimentale des matériaux, à la température d’emploi.

Un calcul thermomécanique ne dit rien, enfin, de plusieurs mécanismes qu’il faut traiter à part. Il ignore le fluage, dimensionnant au-delà d’environ 40 % de la température de fusion absolue du matériau, qui fait dépendre la déformation du temps de maintien autant que de la charge. Il ne préjuge ni de l’oxydation ni du vieillissement d’un revêtement. Et il ne remplace pas l’essai instrumenté : des thermocouples posés aux bons endroits, un cycle réel appliqué, puis la corrélation entre températures mesurées et calculées restent le seul moyen de recaler les conditions aux limites. L’angle mort est toujours le même : un modèle thermique non recalé rend des contraintes justes à la méthode près, et fausses à l’hypothèse près.

FAQ : calcul thermomécanique

Qu’est-ce qu’un calcul thermomécanique ?

C’est un calcul de structure dont le chargement principal est un champ de température plutôt qu’une force extérieure. On résout d’abord la thermique, en conduction, convection et rayonnement, pour obtenir la température en chaque point et à chaque instant. Ce champ est ensuite transféré au modèle mécanique, qui le traite comme une déformation imposée et rend des déplacements et des contraintes. La sanction porte sur la tenue de la pièce alors qu’aucun effort extérieur ne lui est appliqué.

Pourquoi une pièce libre de se dilater ne développe-t-elle aucune contrainte ?

Parce que la contrainte naît d’une déformation empêchée, pas d’un changement de dimension. Une barre chauffée uniformément et libre à ses deux extrémités s’allonge simplement de son allongement thermique : elle change de longueur sans que la matière soit sollicitée. Dès qu’un appui, un gradient interne ou un matériau voisin s’oppose à cet allongement, la part bloquée se convertit en déformation élastique, donc en contrainte. C’est pourquoi une bonne parade consiste souvent à libérer un degré de liberté plutôt qu’à refroidir.

Comment estimer la contrainte d’origine thermique dans une pièce bridée ?

Dans le cas simple d’un bridage complet et d’un échauffement uniforme, la contrainte vaut le produit du module d’Young, du coefficient de dilatation et de l’écart de température. Pour un alliage d’aluminium, 70 GPa et 23 µm/m·K, un écart de 100 °C donne environ 161 MPa. Ce calcul de coin de table sert de premier tri : s’il approche la limite d’élasticité du matériau, une analyse par éléments finis devient nécessaire, car le bridage réel n’est jamais total ni uniforme.

Faut-il toujours mener un calcul en transitoire ?

Non, mais il faut toujours vérifier que le régime établi est bien le cas le plus sévère, ce qui est rarement vrai. Pendant une montée rapide en température, la face exposée se dilate avant le cœur de la paroi et le gradient dans l’épaisseur passe par un maximum, souvent bien au-delà de la valeur d’équilibre. Le nombre de Biot et la diffusivité thermique du matériau donnent le premier repère. Un calcul en régime établi seul suffit quand les sollicitations évoluent lentement devant le temps de réponse de la pièce.

Quel écart de coefficient de dilatation devient problématique dans un assemblage ?

Il n’existe pas de seuil universel : ce qui compte est le produit de l’écart de coefficient, de la longueur d’interface et de l’amplitude de température. Un écart de 22 µm/m·K entre aluminium et composite carbone, sur 500 mm et 100 °C, produit 1,1 mm d’allongement différentiel, qui part intégralement en cisaillement dans les fixations. Le même écart sur 20 mm reste anodin. Réduire la longueur des zones liées, ou intercaler un matériau à coefficient intermédiaire, agit sur les deux termes.