Ministère de l’Agriculture et Ministère de l’éducation
des Ressources hydrauliques National de l’Enseignement
Supérieur et de la Recherche
Institution de la Recherche et de
l‘Enseignement Supérieur
Université de Jendouba
Agricoles
Ecole Supérieure des Ingénieurs
Medjez El Bab
2019/2020
Travaux pratiques
ABAQUS
Réalisé par : AZOUNA Iheb
KRAYNI Bilel
BEN HLIMA Med Atef
SOMMAIRE
I. TP numéro 1 : Simulation d’un essai de traction d’une éprouvette plate ..........................3
1. Introduction .................................................................................................................3
2. Etape 1 : Définition de la géométrie .............................................................................3
3. Etape 2 : Définition du matériau et affectation du matériau à la géométrie ...................4
a. Définition du matériau ..............................................................................................4
b. Création d’une section : ........................................................................................4
c. Affectation des propriétés aux différentes zones de l’éprouvette ...............................4
4. Etape 3 : appel à la géométrie.......................................................................................5
5. Etape 4 : création des steps ..........................................................................................6
6. Etape 5 : Définition des conditions aux limites et du chargement .................................7
7. Etape 6 : Maillage de la structure .................................................................................7
8. Etape 7 : Exécution des calculs ....................................................................................8
9. Visualisation : ..............................................................................................................9
II. TP numéro 2 :Etude thermique Refroidissement d’un processeur ..................................10
1. Position du problème : ............................................................................................... 10
2. Etape 1 : Définition de la géométrie ........................................................................... 11
3. Etape 2 : Définition du matériau et affectation du matériau à la géométrie ................. 12
4. Etape 3 : appel à la géométrie..................................................................................... 12
5. Etape 4 : création des steps ........................................................................................ 13
6. Etape 5 : Définition des conditions aux limites et du chargement ............................... 13
7. Etape 6 : Maillage de la structure ............................................................................... 14
8. Etape 7 : Exécution des calculs .................................................................................. 14
Page | 2
I. TP numéro 1 : Simulation d’un essai de traction d’une éprouvette
plate
1. Introduction
Pour réaliser des calculs par éléments finis, il faut passer par les différentes étapes suivantes
(il existe d’autres étapes mais qui sont nécessaires dans les problèmes de contact ou bien si on
a plusieurs pièces qui sont modélisées) Etape 1 : Définition de la géométrie Etape 2 :
Définition des matériaux et affectation de ses matériaux à la géométrie Etape 3 : Appel aux
géométries qui seront utilisées (assembly) Etape 4 : Définition de la nature de l’analyse, des
étapes de calcul et des résultats recherchés (outputs) Etape 5 : Définition des conditions aux
limites et du chargement Etape 6 : Choix de la stratégie de maillage, de la densité de maillage,
de la nature des éléments et finalement maillage de la structure Etape 7 : Résolution du
problème (exécution du job) Etape 8 : Post-traitement et analyse des résultats
2. Etape 1 : Définition de la géométrie
Dans ABAQUS et dans la majorité des logiciels industriels de calculs par éléments finis, il
existe plusieurs façons pour définir la géométrie de la structure à calculer. Les deux méthodes
les plus utilisées sont : La première consiste à construire soit même la géométrie à partir des
outils disponibles dans le préprocesseur du logiciel et la deuxième consiste à importer cette
géométrie à partir d’un fichier issu d’un logiciel de CAO (format neutre). Pour cette première
séance, vous allez construire la géométrie de l’éprouvette en utilisant le préprocesseur
d’ABAQUS qui est ABAQUS CAE, afin de maîtriser ses différentes fonctionnalités et
options.
Page | 3
3. Etape 2 : Définition du matériau et affectation du matériau à la
géométrie
a. Définition du matériau
Pour définir un ou plusieurs matériaux, il faut se placer dans le module property Puis cliquer
sur create material, type mechanical, elasticity, elastic Ensuite rentrer le module de Young et
le coefficient de Poisson E=210 GPa et nu=0.3
b. Création d’une section :
Category : solid (car c’est un problème plan et non pas une plaque) Homogenious Puis choisir
le matériau relatif à cette section et l’épaisseur de la tranche (2 mm)
c. Affectation des propriétés aux différentes zones de l’éprouvette
Cliquer sur assign section Sélectionner les 5 zones de l’éprouvette et affecter les la section
qu’on vient de définir.
Page | 4
4. Etape 3 : appel à la géométrie
Dans le module assembly on fait appel aux géométries qui seront utilisées lors de l’analyse et
de les positionner les unes par rapport aux autres. Choisir l’option make independant
Page | 5
5. Etape 4 : création des steps
Un step est une étape de calcul. Dans un calcul on peut avoir plusieurs steps. Dans ce cas on
aura besoin que d’un seul step et qui sera de nature static general puisqu’on va simuler un
essai de traction quasi-statique sans effets dynamiques.
1/ aller dans module step
2/ cliquer sur l’icône create step
3/ choisir le type static, general
4/ mettre les autres options par défauts NLGEOM : non (on reste dans le cadre linéaire
géométrique HPP)
Définition des sorties (output)
Il existe deux types de sorties : les field output (sortie de type champ, on obtient l’évolution
spatiale de certaines grandeurs). Les history output (sortie de certaines grandeurs en des
points ou éléments particuliers). On s’intéresse ici à l’évolution des champs de contraintes et
de déformations dans l’éprouvette (pour les field output) et aux déplacements des deux lignes
de repère tracées sur l’éprouvette et de la réaction de l’éprouvette dans les history output.
Pour les field output on les laisse par défauts Pour les history output on va définir trois
ensembles (set) Le premier correspond à la marque du premier repère Le deuxième
correspond à la marque du deuxième repère Le dernier correspond à la frontière de
l’éprouvette qui sera encastrée.
Pour les sets relatifs aux marques des repères, on demande comme sortie le déplacement
suivant l’axe de l’éprouvette Pour le set de la frontière encastrée, on demande comme sortie la
réaction suivant l’axe de l’éprouvette.
Page | 6
6. Etape 5 : Définition des conditions aux limites et du chargement
Les conditions aux limites choisies pour cet essai sont blocage d’une partie de l’éprouvette en
contact avec le mors fixe et un déplacement axial de la partie de l’éprouvette en contact avec
le mors mobile.
Aller dans module load Cliquer sur
create Boundary condition
L’encastrement sera appliquer dans
le step initial Category : mechanical
Types : displacement Mettre U1=0 et
U2=0
Choisir les 3 arêtes définissant la
frontière de l’éprouvette en contact
avec le mors fixe
Refaire la même démarche pour la partie en contact avec le mors mobile sauf que le step sera
le step-1 et on impose comme déplacement U1=-0.1
7. Etape 6 : Maillage de la structure
Page | 7
Aller dans le module mesh
1/ definir la technique adoptée pour
le maillage ainsi que la forme des
éléments (quadrangulaires,
triangulaires ou un mélange de deux)
Cliquer sur assign mesh control
Choisir les zones désirées (les parties
droites de l’éprouvette) Choisir la
forme des éléments (quad) Définir la
technique de maillage : structured
Refaire la même démarche pour les deux zones de l’éprouvette en contact avec les mors et
choisir une technique free avec une option advancing front
2/ définir la taille ou le nombre d’élément par arête (il existe plusieurs techniques, on utilise la
technique edge by nember)
3/ choix des types des éléments (triangulaire, quadrangulaire, linéaire ou quadratique ainsi
que la nature de l’intégration réduite ou complète) et choix de la famille de l’élément qui est
ici plane stress (contraintes planes) On choisit ici des éléments quadrangulaires linéaires à
intégration complète (CPS4) CPS pour dire c’est plane stress et 4 c’est quadrangulaire).
4/ maillage de la géométrie (part entier en un seul coup ou par région)
8. Etape 7 : Exécution des calculs
Module job
Create job : nom du job test traction Exécution du job : job manager submit
Page | 8
9. Visualisation :
Page | 9
II. TP numéro 2 :Etude thermique Refroidissement d’un processeur
1. Position du problème :
La majorité des circuits électroniques sont refroidis par un système comportant un ventilateur
et radiateur. On se propose d'effectuer une analyse de la dissipation de la chaleur dans
l'élément radiateur. La géométrie du dispositif (domaine d'étude) est donnée. Elle est décrite
par la longueur L et la section transversale S, la température T0 générée par le processeur et la
température ambiante Te. Pour déterminer : la distribution de la température T dans l'ailette,
nous supposons que la largeur W est grande devant les autres dimensions : W >> H et W >> L
de sorte que la distribution va être la même dans toute section transversale (x,y). Où H est
l’épaisseur de l’ailette et L sa longueur.
Considérons le domaine de calcul défini par la section transversale de l'ailette, représenté par
la figure suivante.
L’équation de la chaleur est :
∇ .(k ∇ T(x , y))=ρCp ∂T (x , y)/ ∂ t
Les pertes latérales par convection thermique sont prises en compte au niveau des conditions
aux limites qui s'écrivent :
Page | 10
Sur Γ1( y= 0) la température est imposée : T = T0.
Sur Γ3(y = L) nous prenons la condition de flux nul.
Sur Γ2(x = H/2) et Γ4 (x=-H/2) nous prenons une condition de convection avec hc coefficient
de convection et Te température à l’infini.
On remarque que le problème est symétrique par rapport à l’axe des abscisses (y=0). On
pourra prendre uniquement la moitié de l’ailette avec une condition de symétrie pour la
frontière (y=0). L’ailette est en alliage d’aluminium dont les propriétés mécaniques et
thermiques sont :
• Conductivité thermique k=200 W.m-1.K-1.
• Capacité calorifique Cp=900 J.Kg-1.K-1.
• Masse volumique ρ=2700 Kg/m3.
Les dimensions des ailettes sont : W=60 mm ; H=1 mm et L=16 mm. L’épaisseur de la plaque
qui supporte les ailettes E=4 mm et la distance entre deux ailettes est p=5 mm hc=1000W/m2
K, T0=80°C et Te=20°C.
2. Etape 1 : Définition de la géométrie
Page | 11
3. Etape 2 : Définition du matériau et affectation du matériau à la
géométrie
4. Etape 3 : appel à la géométrie
Page | 12
5. Etape 4 : création des steps
6. Etape 5 : Définition des conditions aux limites et du chargement
Page | 13
7. Etape 6 : Maillage de la structure
8. Etape 7 : Exécution des calculs
Page | 14
Page | 15
