Modélisation et simulation numérique

TP2 : Simulation thermomécanique

sous ANSYS

Objectif du TP

L’objectif de ce TP est de se familiariser avec les outils de simulation thermomécanique sous

ANSYS. Les modules de calculs thermiques et élastiques en statique seront abordés, un exemple
d’étude transitoire thermique sera introduit et la mise en place d’études couplées thermique et
élastique sera réalisée. Le TP est divisé en 4 parties distinctes. La première partie consiste en
l’étude mécanique en statique d’une poutre cantilever en flexion. La seconde partie est dédiée
à l’étude thermique statique et transitoire de la pièce d’entrée d’un catalyseur de véhicule
automobile. La troisième partie est dédiée à une étude mécanique d’un composite fibre/polymère
sous condition thermique. La dernière partie introduit le calcul couplé thermoélastique via un
exemple simple de pièce circulaire.

1 Etude mécanique d’une poutre cantilever

Nous considérons l’étude mécanique d’une poutre en alliage d’aluminium de section carrée
(40×40mm2 ) et de longueur 500mm. La poutre de fibre neutre droite orientée selon z est
encastrée en une extrémité et soumise à une force verticale (direction y) descendante de 500N
à son extrémité opposée.

Figure 1 : Poutre cantilever

Leclerc Willy 1 PUSQ

Modélisation et simulation numérique

1.1 Schéma du projet

Ouvrir un nouveau projet et consid-

érer le schéma de projet ci-contre.

Figure 2 : Schéma du projet d’étude statique

1.2 CAO

Dessiner la géométrie de la poutre de section carrée (40×40mm2 ) et de longueur 500mm dans

le modeleur de votre choix.

1.3 Etude mécanique

Vérifier que la géométrie est bien chargée et mise à jour. Pour l’affectation du matériau,
sélectionner l’Alliage d’aluminium (module de Young 71 GPa, coefficient de Poisson 0,33) dans
l’onglet matériaux de l’arborescence du module Structure Statique. Celui-ci doit normale-
ment s’ajouter à l’acier dans la liste des matériaux du projet. Afin d’affecter ce matériau au
solide étudié, aller vers l’onglet géométrie de l’arborescence du module Structure Statique
et remplacer le matériau d’affectation au besoin.

Figure 3 : Chargement d’un matériau existant dans l’arborescence : Figure 4 : Affectation

l’Alliage d’aluminium du matériau au solide

Leclerc Willy 2 PUSQ

Modélisation et simulation numérique

Pour l’étape de maillage, considérer un maillage par défaut (grossier, résolution 2) puis définir
les chargements et les supports suivants. La poutre est encastrée en une extrémité et soumise
à une force orthogonale transverse à la poutre de module 500N à son extrémité opposée. Nous
nous intéressons à l’étude du déplacement vertical de la poutre (direction y) et de la contrainte
normale (direction z, longitudinale à la fibre neutre) au sein de la poutre. Evaluer les champs
correspondants, sonder le déplacement à l’extrémité chargée de la poutre et faire une compara-
ison avec la valeur théorique donnée par la théorie des poutres d’Euler-Bernoulli:

P L3
uy (z = L) = −
3EI

où E est le module de Young du matériau de la poutre, I et L, le moment quadratique et la

longueur de la poutre, et P la force appliquée. De même, évaluer la contrainte à l’extrémité
encastrée et faire une comparaison avec la valeur théorique donnée par la formule suivante :
P Lh
σz (z = 0) =
I
où σz est la contrainte normale selon l’axe de la fibre neutre de la poutre, h est la demi hauteur
de section de la poutre (soit 20mm) et I est le moment quadratique de celle-ci.

Figure 5 : Champ de déplacement directionnel Figure 6 : Champ de contrainte normale

(dir. y) (dir. z)

1.4 Etude en convergence de maillage

Réaliser une étude en convergence de maillage de la solution en contraintes déterminée

précédemment. Pour cela, utiliser l’outil de convergence disponible en cliquant droit sur l’onglet
Contrainte Normale. Nous considérons les paramètres de convergence donnés sur les Figures
7a et 7b et les paramètres de raffinement par défaut (disponibles directement dans Solution).

Figure 7a : Définition de la modification Figure 7b : Définition des paramètres de

maximale autorisée raffinement adaptatif
Etudier la courbe obtenue. Qu’en déduit-on sur la convergence de la contrainte normale ?
Expliquer le résultat.

Leclerc Willy 3 PUSQ

Modélisation et simulation numérique

2 Etude thermique d’une pièce de catalyseur

Nous considérons l’étude thermique d’une pièce divergente située à l’entrée d’un catalyseur
automobile. La pièce est axisymétrique, c’est à dire qu’elle est invariante par rotation autour
d’un axe central. Nous nous intéressons à la déperdition thermique au sein de l’épaisseur du
solide. Deux études seront à réaliser dans ce but : une stationnaire et une transitoire.

Figure 8 : Pièce de catalyseur

2.1 CAO

L’esquisse à réaliser est donnée en Figure 9. Dessiner cette esquisse et générer la pièce par
révolution autour de l’axe central (l’axe x sur la Figure 9). Les dimensions à considérer sont les
suivantes:

Figure 9 : Esquisse de la pièce de catalyseur

• V1=40mm • H4=45mm

• V3=10mm • H5=47,5mm

• V2=75mm • H6=20mm

Leclerc Willy 4 PUSQ

Modélisation et simulation numérique

2.2 Propriétés matérielles

La pièce de catalyseur est constituée d’alliage d’aluminium. Définir ce matériau dans les pro-
priétés matérielles, vérifier qu’il est bien affecté à notre solide. Noter sa conductivité thermique.

2.3 Etude thermique stationnaire

L’étude thermique stationnaire est ef-

fectuée à l’aide du module Ther-
mique stationnaire. Le plan du pro-
jet à réaliser est donnée ci-contre.

Figure 10 : Schéma du projet d’étude thermique

stationnaire
A l’aide de l’arborescence, définir une température initiale de 300◦ C. Nous supposons les con-
ditions aux limites suivantes :

• Température fixée à 500 ◦ C sur la surface intérieure (3 surfaces à sélectionner dans le

modèle)
• Convection sur la surface extérieure (3 surfaces à sélectionner également) avec un coeffi-
cient de convection thermique de 3,5e-5 W/(mm2 K) et une température ambiante de 300
◦ C.

Figure 11 : Carte de température au sein du solide (cadre stationnaire)

Nous nous intéressons à la température au sein du solide. Définir un résultat de température et

lancer le calcul. Evaluer le différentiel de température entre les surfaces interne et externe de la
pièce. Comparer le résultat avec la théorie. Commenter.

Leclerc Willy 5 PUSQ

Modélisation et simulation numérique

2.4 Etude thermique transitoire

Nous nous intéressons maintenant à une étude thermique transitoire. Le schéma du projet à
définir est donné sur la Figure 12.

Figure 12 : Schéma du projet d’étude thermique transitoire

Suivre les étapes de l’assistant de

calcul thermique transitoire. Vérifier
le matériau (Alliage d’aluminium)
et contrôler l’incrément avec les
paramètres suivants :

• Nombre d’incréments : 1

• Incrément courant : 1

• Temps final pour cet incrément

: 10s

• Incrément de temps initial 0,1s

• Incrément de temps minimal 1e-

2s

• Incrément de temps maximal 1s

Figure 13 : Etapes de l’assistant de calcul thermique
transitoire
Lancer le calcul et visualiser la montée en température sur la surface extérieure de la pièce.
Comparer le résultat final en régime permanent avec le résultat obtenu dans le calcul station-
naire.

3 Etude élastique sous condition thermique d’un composite fi-

bre/polymère
Nous nous intéressons à l’étude élastique sous condition thermique d’un composite fi-
bre/polymère. Le modèle considéré est une plaque carrée avec un disque au centre représentant

Leclerc Willy 6 PUSQ

Modélisation et simulation numérique

Figure 14 : Carte de température au sein du solide (cadre transitoire)

la fibre dans un plan transverse. Nous nous intéressons à une étude mécanique (module Struc-
ture statique) sous condition de température de 22◦ C. Le modèle à réaliser est visible sur la
Figure 15.

3.1 CAO

Dessiner l’esquisse de la plaqué carré de côté 40mm×40mm (garder le disque de côté pour
l’instant). Générer d’abord la plaque carré seule par extrusion, puis sélectionner la surface
supérieure et définir une nouvelle esquisse sur laquelle vous dessinerez le disque en vous référant
à la Figure 15.

Figure 15 : Modèle de composite fibre/polymère

Leclerc Willy 7 PUSQ

Modélisation et simulation numérique

Les cotations à respecter sont les suivantes :

• R1 (rayon du disque du carré) = 10mm

• V2 (demi-longueur verticale du carré) = 20mm

• H3 (demi-longueur horizontale du carré) = 20mm

• L’épaisseur de la plaque est de 4mm

Trouer la plaque à l’aide de l’outil d’extrusion et nommer le nouveau corps obtenu "matrice".
Ensuite, générer la fibre circulaire toujours à l’aide de l’outil d’extrusion (option ajout d’un
corps bloqué). Nommer ce second corps "fibre".

3.2 Propriétés matérielles

Afin de tenir compte du caractère hétérogène du matériau composite, il nous faut affecter à
chaque solide un matériau distinct. La fibre est constituée d’alliage d’aluminium, le reste de la
plaque est constitué d’un polymère aux propriétés suivantes :

• Coefficient de dilatation : 2,34e-6 K−1

• Module de Young : 3 GPa

• Coefficient de Poisson : 0,3

Figure 16 : Création du matériau "polymer" dans une nouvelle bibliothèque personnalisée

"Perso"

Le matériau polymère n’existe pas dans la base de données par défaut de ANSYS. Il nous faut
donc le créer, ce qui peut se faire via les données techniques (case B2 du module Structure
statique). Dans ce but, il nous faut tout d’abord créer une nouvelle bibliothèque personnalisée
(nommée "Perso" sur la Figure 16) puis éditer cette bibliothèque pour créer le nouveau matériau.
Afficher les propriétés matérielles et définir celles du polymère comme indiqué sur la Figure 17.

Leclerc Willy 8 PUSQ

Modélisation et simulation numérique

Figure 17 : Définition des propriétés du matériau "polymer"

3.3 Etude statique

Réaliser une étude statique du composite. Nous supposons une condition thermique de 22◦ C.
Vérifier les matériaux. Le matériau polymère doit être associé au corps "matrice" et l’alliage
d’aluminium au corps circulaire "fibre". Vérifier également la surface de contact. Les conditions
aux limites à considérer sont les suivantes :

• Bord surfacique gauche fixé

• Bord surfacique droit soumis à une force de 50 N dirigée selon la normale extérieure à la
surface

Nous nous intéressons au calcul de la contrainte équivalente de Von-Mises et à l’intensité de

déformation élastique. Effectuer les calculs et expliquer les résultats.

Figure 18 : Carte des contraintes au sein du Figure 19 : Carte des déplacements au sein du
composite composite

4 Etude thermoélastique d’une pièce circulaire

Nous nous intéressons à l’étude couplée thermique et élastique (thermoélastique donc) d’une
pièce simple. La pièce à étudier est soumise à la fois à un gradient de température et à un
chargement mécanique.

Leclerc Willy 9 PUSQ

Modélisation et simulation numérique

Figure 20 : Pièce circulaire à étudier

4.1 Organigramme du projet

L’organigramme du projet à générer est celui visible sur la Figure 21. Il fait intervenir à la fois
un calcul thermique stationnaire et un calcul statique mécanique. Afin de le générer rapidement,
sélectionner l’onglet Thermomécanique dans Systèmes Personnalisés.

Figure 21 : Organigramme du projet

4.2 Propriétés matérielles

Nous supposerons que l’anneau circulaire est constitué d’acier (matériau par défaut).

4.3 CAO

La pièce à générer a un rayon intérieur de 5 pouces et un rayon extérieur de 10 pouces.

L’épaisseur de la pièce est de 1 pouce.

4.4 Etude thermomécanique

Les conditions aux limites thermiques à définir sont les suivantes :

• Température imposée de 250◦ C sur la surface intérieure

• Température imposée de 25 ◦ C sur la surface extérieure.

Leclerc Willy 10 PUSQ

Modélisation et simulation numérique

Figure 22 : Conditions aux limites considérées (étude thermomécanique)

La surface intérieure est soumise à une condition aux limites mécanique également. Il s’agit d’une
pression de 2 MPa dirigée vers la normale intérieure à cette surface intérieure. Dans un premier
temps, nous considérons le problème thermique stationnaire. Considérer une température ini-
tiale de 22◦ C (onglet Température Initiale). Nous nous intéressons à la température au sein
du solide. Lancer la résolution du problème thermique. Dans un second temps, nous considérons
l’étude statique élastique. Les résultats de l’étude thermique déjà réalisée sont automatiquement
reliés à l’étude élastique. Définir une température d’environnement de 22◦ C (dans les options de
l’onglet du module Structure statique). Nous nous intéressons au contraintes de Von-Mises
au sein de la pièce. Lancer la résolution du problème statique.

Figure 23 : Carte des températures au sein du Figure 24 : Carte des contraintes de Von-Mises
composite au sein du composite

Leclerc Willy 11 PUSQ