RÉPUBLIQUE TUNISIENNE
MINISTÈRE DE L’ENSEIGNEMENT SUPÉRIEUR ET DE LA RECHERCHE
SCIENTIFIQUE
INSTITUT SUPÉRIEUR DES ÉTUDES TECHNOLOGIQUES DE SOUSSE
Polycopié de Travaux Pratiques
Logiciels de simulation numérique en mécanique
ABAQUS
Master
Génie Mécanique
Présenté par
Mme : Adalet Narjess
Plan
TPN° 1 : Prise en main d’ABAQUS
TPN° 2 : Etude de la flexion d’une poutre Encastrée-libre ( modélisation 1D)
TPN° 3 : Etude de la flexion d’une poutre Encastrée-libre ( modélisation 2D)
TPN° 4 : Etude de la flexion d’une poutre Encastrée-libre ( modélisation 3D)
TP N° 5 : Déformation plastique (Simulation d’un essai de traction)
TP N° 6 : Flexion 3 points (modélisation 1D)
TPN° 2 :
Etude de la flexion d’une poutre
Encastrée-libre ( modélisation
1D)
TP N° 2 : Etude de la flexion d’une poutre encastrée-libre
(modélisation 1D)
2.1 Introduction
Dans ce travail pratique, nous allons simuler la flexion d’une poutre encastrée par
l’une de ces extrémités et libre de l’autre. Ce problème possède une solution
analytique issue de la théorie des poutres. Ceci nous permettra de comparer les
résultats analytiques et ceux obtenus par les simulations numériques effectuées.
2.2 Position du problème
On considère une poutre de section rectangulaire S = 5x2 mm² et de longueur L=150
mm, encastrée par l’une de ses extrémités et libre de l’autre. Nous appliquons une
charge concentrée F sur le coté libre de la poutre F = 15 N.
La figure (Fig.2.1) présente les conditions aux limites du problème. La poutre est
supposée élastique (E=210 GPa et υ =0.3).
2.3 Solution analytique
Pour déterminer la valeur de la flèche maximale, on se propose d’étudier analytiquement la flexion de
cette poutre, encastrée par l’une de ses extrémités et libre de l’autre.
1. Les réactions au niveau de l’encastrement :
2. Effort tranchant et moment fléchissant :
2.4 Solution numérique
2.4.1Etapes de la modélisation
Pour modéliser sous ABAQUS CAE la flexion de cette poutre, on suit les
étapes suivantes :
1. Créer la géométrie de la poutre (module Part)
2. Créer le matériau (module Proprety)
3. Définir et affecter les propriétés de section (module Proprety)
4. Assembler le modèle (module Assembly)
5. Définir le pas d’analyse (module Step)
6. Appliquer des conditions aux limites et des chargements sur le modèle
(module Load)
7. Mailler le modèle (module Mesh)
8. Créer et soumettre un travail d’analyse (module Job)
9. Visualiser les résultats (module Visualizaion)
2.4.1 Réalisation du modèle sous ABAQUS
Système d’unités
Avant de commencer de définir n'importe quel modèle sous ABAQUS, on
doit décider quel système d’unités, on emploie. ABAQUS n'a aucun
système d'unités intégré. Toutes les données doivent être spécifiées dans
un système d'unités conforme. Nous adoptons dans ce TP le système
d'unités SI (mm).
Longueur Mm
Force Newton
Pression/contrainte MPa
Densité Tonnes/mm3
Masse Tonnes
température Kelvin
Hypothèse de modélisation
La poutre est modélisée sous ABAQUS selon l’hypothèse 1D (Fig.2.3)
. 1D
Fig.2.3
Modélisation
1. Module PART
Dans l’arbre du modèle cliquer deux fois sur Parts
Dans la boite de dialogue Create Part Renommer la pièce
Poutre_1D Sélectionner 3D/Déformable /Solid/Wire Mettre
Approximative size = 200
Cliquer Continue (Fig.2.4)
Créer une ligne (Line) de 150 mm
Cliquer Done (Fig.2.5)
Enregistrer le fichier de travail ([Link])
Fig.2.5 Fig.2.4
2. Module PROPRETY
Définition du matériau et ses propriétés
Nous allons créer un seul matériau élastique linéaire avec les propriétés
suivantes : module de Young E= 210 000 MPa et coefficient de Poisson =0.3.
Dans l’arbre du modèle double-cliquer sur Materials
Renommer le matériau Acier
Donner une description (au choix).
Cliquer sur l’onglet:
Mechanical/Elasticity/Elastic Définir le
module d’Young et le coefficient de Poisson
Cliquer OK (Fig.2.6)
Affectation des propriétés de la section:
Dans l'arbre modèle double-cliquer sur Sections (Fig.2.7)
Dans la boite de dialogue Create Section (Fig.2.8)
Fig.2.7
Renommer la section Section_Poutre-1D Dans la liste Category, sélectionner Beam
Dans la liste Type, sélectionner Beam
Cliquer sur Dans la boite de dialogue Edit Beam Section qui s’affiche
(Fig.2.9)
Cliquer sur Create Beam Profile
La fenêtre Create Profile s’affiche (Fig.2.10)
Choisir Rectangular et cliquer sur Continue
La fenêtre Edit Create Profile apparait (Fig.2.11)
Remplir le champ a=2 et b=5 Cliquer sur OK
Fig.2.8
Fig.2.6
Fig.2.9
Assignement de la section
Pour affecter la section Section_Poutre-1D à la structure, double-cliquer
sur Section Assignements dans la liste des attributs de la pièce qui apparaît.
ABAQUS/CAE affiche des messages-guide dans la zone de message-guide
pour vous aider (Fig. 2.12). Sélectionner la pièce entière comme la région
à laquelle la section sera appliquée. ABAQUS/CAE mets la structure en
surbrillance (Fig.2.13). Dans la boite de dialogue Edit Section
Assignement accepter la sélection par défaut de la section Section_Poutre-
1D puis cliquer OK. ABAQUS/CAE affecte la section à la poutre et colore
la structure entière en turquoise (Fig.2.14). Cliquer sur Done.
Affectation de l’orientation de la
poutre Dans barre de menu double-
cliquer sur Assign Cliquer sur Beam
Section Orientation
Choisir le système par défaut d’ABAQUS (Fig.2.15)
Cliquer sur OK (Fig.2.16)
Fig.2.15 Fig.2.16
3. Module ASSEMBLY
Assemblage
Le module ASSEMBLY permet d’assemblé les différentes pièces créées
précédemment. Cette étape est nécessaire même si un projet n’est composé
que d’une seule pièce. Pour ce problème, on va créer une seule instance de
la poutre objet de cette étude.
Dans l'arbre modèle double-cliquer sur Assembly
Double cliquer sur Instances dans la liste qui apparait
Dans la boite de dialogue Create Instance sélectionner Poutre_1D (Fig.2.17)
Cliquer OK (Fig.2.18).
Fig.2.18
Fig.2.17
4. Module STEP
Ce module permet d’effectuer les étapes de calcul. Dans notre cas il y a
une seule étape de calcul que l’on définit de la manière suivante :
Configuration du pas de l'analyse
Nous nous intéressons à la réponse statique de la poutre sollicitée en
flexion par une charge concentrée F=15 N appliquée à l’extrémité libre de
la poutre, un seul pas d'analyse est nécessaire dans cette simulation.
L'analyse se compose au total de deux pas : Un premier pas, dans lequel
on va appliquer des conditions aux limites de déplacement imposée sur
l’une des extrémités de la poutre.
Un second pas d'analyse, dans lequel on applique la charge
concentrée sur l’autre extrémité libre de la poutre.
Création du pas d'analyse
On crée un pas d'analyse statique qui suit le pas initial d'analyse
généré par ABAQUS/CAE.
Double cliquer sur Steps pour créer un pas.
La boîte de dialogue Create Step apparaît (Fig.2.19). Une liste de
toutes les procédures générales et un nom du pas Step-1 est fourni par
défaut.
Renommer le pas : Chargement
Sélectionner Static General comme Procedure type
Cliquer sur Continue.
La boîte de dialogue Edit Step apparaît avec les paramètres par défaut pour
un pas d'analyse Static General (Fig.2.20)
L'onglet Basic est sélectionné par défaut.
Dans le champ Description (au choix) : taper force concentrée 15 N
Accepter les valeurs par défaut fournies pour le pas.
Cliquer sur OK pour créer le pas et pour quitter la boîte de dialogue Edit Step.
Fig.2.19 Fig.2.20
5. Module LOAD
Le module LOAD permet de définir le chargement subit par la poutre et
d’introduire les conditions aux limites qui sont appliquées à des régions du
modèle constituées de déplacements et/ou des rotations. Les Conditions
imposées, telles que les chargements et les conditions aux limites
dépendent du pas d'analyse. On doit donc spécifier le ou les pas dans
lesquelles ils deviennent actifs.
Application des conditions aux limites à la poutre
Double cliquer sur la base BCs
La boîte de dialogue Create Boundry Condition s'affiche (Fig.2.21)
Renommer la condition aux limites Encastrement
Dans la liste des pas Steps, sélectionner Initial
Dans la liste Category, accepter Mechanical comme catégorie par défaut.
Sélectionner Symmetry/Antisymmetry/Encastre dans la liste Types
for Selected Step (Fig.2.21)
Choisir Encastre (Fig.2.22)
Cliquer sur Continue (Fig.2.23)
Fig.2.21 Fig.2.22
Fig. 2.23
Application du chargement à la
poutre Dans le module LOAD, cliquer
sur Create Load La fenêtre Create
Load apparait
Renommer la charge : Force concentrée
Choisir Mechanical/Concentred force et cliquer sur Continue
(Fig.2.24) Sélectionner le point de l’application de la charge et
cliquer sur Done (Fig.2.25) La fenêtre Edit Load apparait (Fig.2.26)
Remplir le champ réservé aux composante de la force (0,-15,0) et cliquer
sur OK (Fig.2.27)
Fig.2.25
Fig.2.24
Fig.2.27
Fig.2.26
6. Module MESH
Maillage du modèle
Dans cette étape, nous allons générer le maillage. ABAQUS donne la
possibilité de choisir la technique de maillage, le type des éléments ainsi
que leur forme.
Attribution du type d'élément
Dans cette section, on va attribuer le type d'élément pour le modèle. Dans
un premier temps nous avons utilisé un élément de type B31 (Fig.2.28)
avec les paramètres de maillage suivants :
- Taille de la maille : 4 mm
- Degré d’interpolation : Linéaire
- Type: Beam
Démarche
Revenir au module Assembly et cliquer sur Instances
Activer Make Independent
(Fig.2.29) Dans le module MESH
cliquer sur Mesh L’instance est
coloriée en rose (Fig.2.30)
Cliquer sur Seed Edges et sélectionner la poutre cliquer sur OK (Fig.2.31)
La fenêtre Local Speeds apparait (Fig.2.32)
Choisir 38 comme nombre d’éléments et cliquer sur OK (Fig.2.32)-(Fig.2.33)
Cliquer sur Mesh instance et confirmer le maillage (Fig.2.34)
Fig.2.28 Fig.2.29
Fig.2.30
Fig.2.31
Fig.2.32 Fig.2.33
Fig.2.34
7. Module JOB
Créer et soumettre un travail d’analyse
Double cliquez sur Jobs pour ouvrir la boîte de dialogue Create Job
Créer un nouveau travail (Fig.2.35)
Nommez le Poutre_1D_4mm et cliquer OK
Pour soumettre le travail, Cliquer sur Submit (Fig.2.36)
A la fin de l’analyse (ABAQUS affiche Completed Successfully) (Fig.2.37)
Pour visualiser les résultats, cliquer sur Results (Fig.2.38)
Fig.2.35 Fig.2.36 Fig.2.38
Fig.2.37
2.4.2 Problématique
Pour ce problème de modélisation 1D, on se propose d’effectuer le calcul
avec une taille de maillage de 4 mm, de tracer la déformée de la fibre neutre
suivant l’axe y et de discuter les résultats trouvés. Une Comparaison des
résultats numériques avec les résultats théoriques obtenus par le calcul
RDM (solution analytique) sera faite
2.4.3 Résultats
Pour le maillage nous avons choisis les paramètres
suivants : Taille de maille : 4 mm
Degré d’interpolation : linéaire
Type d’éléments : beam
Fig.2.39
La courbe de la flèche de la fibre neutre présentée sur la figure .2.39,
nous montre une valeur maximale de la flèche de -3.86 mm, soit une
erreur relative de presque
0.07 % par rapport à la solution analytique (yan = -3,857 mm).
Traçage de la courbe
Pour tracer la flèche sur toute la longueur de la poutre on suit les étapes
suivantes : Dans le module VISUALISATION, cliquer sur
Tools/Path/Create (Fig.2.40)
La fenêtre Create Path s’affiche (Fig.2.41)
Choisir Edge list et cliquer sur Continue
La fenêtre Edit Edge list Path s’affiche, cliquer sur Add Before (Fig.2.42)
Sur la barre de menu qui s’affiche en bas, choisir By shortest
distance Revenir dans le ViewPort et piquer l’origine et
l’extrémité de la poutre (Fig.2.43) Cliquer sur OK
Cliquer sur l’onglet XY Data Manager (Fig.2.44), la fenêtre XY Data
Manager
s’affiche (Fig.2.45)
Cliquer sur Create, choisir Path et cliquer sur Continue (Fig.2.45)
La fenêtre XY Data from Path apparait (Fig.2.46)
Pour afficher la courbe cliquer sur Plot (Fig.2.47)
Pour afficher les données Data de la courbe : cliquer sur Edit (Fig.2.48)
Fig.2.40 Fig.2.41 Fig.2.42
Fig.2.43
Fig.2.44
Fig.2.45
Fig.2.46
Fig.2.47
Cours : Simulation numérique par Abaqus Iset Sousse
Fig.2.47
20 40 60 80 100 120 140 160
-0,5
-1 Flèche y(x)
Flèche y(x) ABAQUS
-1,5
Déplacement (mm)
-2
Flèche de la fibre neutre (F=15 N)
-3
-3,5
-4 Longeur (mm)
Fig.2.48
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
2.4.4 Conclusion
Les courbes numérique (ABAQUS) et analytique (RDM) sont en
parfaite concordance (une erreur relative de presque 0.07
%).L’élément choisi pour cette simulation conduit à des résultats très
précis et très proches de la solution analytique.
Année universitaire 2023-2024 Narjes Adalet
