Cours : Simulation numérique par Abaqus Iset Sousse
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
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
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)
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
TPN° 3 :
Etude de la flexion d’une
poutre Encastrée-libre (
modélisation 2D)
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
TP N° 3 : Etude de la flexion d’une poutre encastrée-libre
(modélisation 2D)
3.1 Introduction
Le but principal de ce TP est l’étude de la convergence d’un maillage, pour cela
nous allons considérer la flexion d’une poutre encastrée-libre. Ce cas déjà été
examiné, possède une solution analytique issue de la théorie de la résistance
des matériaux (RDM). Cette solution servira comme solution de référence et
sera comparée par la suite aux résultats numériques obtenus par des simulations
avec différentes densités de maillages.
3.2 Position du problème
Pour la modélisation 2D, on considère les mêmes données du problème déjà
traité précédemment. Soit une poutre de section rectangulaire S =5x2 mm² et
de longueur L=150 mm, encastrée par l’une de ses extrémités, sur l’autre on
applique une charge concentrée F= 15 N. La figure 3.1 illustre les conditions
aux limites du problème. La poutre est supposée élastique conçue en acier
ordinaire avec un module d’Young (E=210 GPa) et un de coefficient de poisson
(υ=0.3).
3.3 Solution analytique
En fonction des caractéristiques du matériau et des paramètres géométriques de la
poutre, la solution analytique de ce problème permettant de calculer la flèche est
donnée par la relation suivante :
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
Les résultats de calcul de la flèche pour différentes valeurs de l’abscisse x sur la
longueur de la poutre pour un pas de 10 mm et qui sont portés sur le tableau 3.1
permettent de la courbe analytique de la flèche (Fig.3.2).
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
3.4 Solution numérique
3.4.1 Réalisation du modèle sous ABAQUS
• Hypothèse de modélisation La poutre est modélisée en 2D (Fig.3.3).
• Modélisation
8. Module PART
Cliquer deux fois sur Parts
Dans la boite de dialogue Create Part
Renommer la pièce Poutre_2D
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
Sélectionner 2D Planar/Déformable /Shell
Choisir Approximative size = 200
Cliquer Continue (Fig.3.4)
Créer un rectangle de 5 x 150 mm
Cliquer Done (Fig.3.5)
Enregistrer le fichier de travail : (POUTRE_2D.cae)
9. Module PROPRETY
• Définition du matériau
Double cliquer sur Materials
Renommer le matériau Acier
Donner une description (au choix)
Cliquer sur l’onglet: Mechanical/Elasticity/Elastic
Introduire le module d’Young (E=210 000 MPa)
Introduire le coefficient de Poisson (ν=0.3)
Cliquer OK (Fig.3.6)
• Affectation de la section
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
Dans l'arbre modèle double-cliquer sur Sections (Fig.3.7)
Dans la boite de dialogue Create Section (Fig.3.8)
Dans la liste Category, sélectionner Solid
Dans la liste Type, sélectionner Homogeneous Fig.3.6
Cliquer sur Continue
La boîte de dialogue Edit Section apparaît (Fig.3.9)
Accepter la sélection par défaut du matériau Acier
Cocher la case Plane stress/strain
Dans la case thikness introduire la valeur 2 qui corresponds à la largeur de la
poutre (Fig.3.9).
• Affectation de la section à la structure
Double-cliquer sur Section Assignements dans la liste des attributs de la pièce
qui apparaît.
Sélectionner la pièce entière comme région sur laquelle la section sera appliquée
ABAQUS/CAE mets la structure en surbrillance (Fig.3.10)
Dans la boite de dialogue Edit Section Assignement accepter la sélection par
défaut de la section Section-1 puis cliquer OK (Fig.3.11)
ABAQUS/CAE affecte la section à la poutre et colore la structure entière en
turquoise (Fig.3.12)
Cliquer sur Done
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
10. Module ASSEMBLY
• Assemblage
Dans l'arbre modèle développer la base Assembly
Double cliquer sur Instances dans la liste qui apparait (Fig.3.13)
Dans la boite de dialogue Create Instance sélectionner Poutre_2D et cliquer OK
ABAQUS/CAE colore la structure entière en bleu et affiche le référentiel
associé (Fig.3.14)
11. Module STEP
• Création du pas d'analyse
On crée un pas d'analyse statique qui suit le pas initial de l'analyse
Double cliquer sur Steps pour créer le pas d’analyse
ABAQUS/CAE passe au module Step et la boîte de dialogue Create Step
apparaît (Fig.3.15)
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 (Fig.3.15)
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
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
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 (Fig.3.16)
Cliquer sur OK pour créer le pas et pour quitter la boîte de dialogue Edit Step
12. Module LOAD
• Application des conditions aux limites à la poutre
Dans l'arbre modèle double cliquer sur la base BCs. ABAQUS/CAE passe au
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
module LOAD et la boîte de dialogue Create Boundry Condition s'affiche
(Fig.3.17). Dans la boîte de dialogue Create Boundry Condition. Renommer la
condition aux limites Encastrement. Dans la liste des pas Steps. Sélectionner
Initial, le pas dans le quel la condition limite sera activée. Toutes les conditions
aux limites mécaniques spécifiées au pas Initial doit avoir grandeurs nulles.
Cette condition est appliquée automatiquement par ABAQUS/CAE. Dans la
liste Category, accepter Mechanical comme catégorie par défaut. Dans la liste
Types for Selected Step, sélectionner Symmetry/Antisymmetry/Encastre et
cliquer sur Continue. Sélectionner l’arrête qui correspond à l’extrémité de la
poutre et cliquer sur Done. Dans la fenêtre Edit Boundary Condition qui apparait
choisir ENCASTRE et cliquer sur OK (Fig.3.18). ABAQUS/CAE affiche la
structure avec encastrement sur le coté fixe(Fig.3.19).
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
• Application du chargement à la poutre
Dans le module LOAD
Cliquer sur Create Load, la fenêtre Create Load apparait (Fig.3.20)
Sélectionner Mechanical/Concentred force et cliquer sur Continue
Sélectionner le point qui se trouve sur le coin haut de l’extrémité libre de la
poutre,
cliquer sur Done
La fenêtre Edit Load apparait. Remplir le champ réservé aux composante de la
force (0,-15,0) et cliquer sur OK (Fig.3.21)
ABAQUS/CAE affiche la structure avec chargement sur le coté libre (Fig.3.22)
13. Module MESH
• Attribution du type d'élément Dans ce module, nous allons attribuer le type
d'élément (quad structured CPS4R) (Fig.3.23), sa taille (5mm) ainsi que le degré
(linéaire) et le type d’interpolation choisi (intégration réduite) (Fig.3.24). Dans
cette première simulation nous utilisons les paramètres du maillage suivants :
Taille de la maille : 5 mm
Degré d’interpolation : Linéaire
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
Type d’intégration : Intégration réduite
Type d’éléments : Quadrilatère
• Démarche
Dans le module MESH cliquer sur Mesh. La poutre est coloriée en vert (Fig.3.25).
Cliquer sur Seed Edges, la poutre est coloriée en rose (Fig.3.26). Sélectionner tout
d’abord les arrêtes qui se trouvent sur la longueur de la poutre et cliquer sur OK
La fenêtre Local Seeds apparait (Fig.3.27), sélectionner By number et choisir 30
comme nombre d’éléments. Ensuite, cliquer de nouveau sur Seed Edges et
sélectionner les arrêtes se trouvant sur les faces encastrée et chargée de la poutre
et cliquer sur OK. La fenêtre Local Seeds apparait de nouveau, choisir 1 comme
nombre d’éléments (Fig.3.29). A la fin cliquer sur Mesh instance et confirmer le
maillage (Fig.3.30).
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
14. Module JOB
• Créer et soumettre le job
Double cliquez sur Jobs pour ouvrir la boîte de dialogue Create Job et créer un
nouveau travail. Nommez le Poutre_2D_Lin_5mm. Ensuite Continue pour créer
le travail. L’Edit Job apparaît (Fig.3.31). Dans description, notez Flexion simple
de la poutre modélisée en 2D. Cliquez OK pour valider. Pour soumettre le travail,
Cliquer sur Submit (Fig.3.32). A la fin de l’analyse du modèle (ABAQUS marque
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
Completed Successfully. Pour visualiser les résultats, cliquer sur Results.
ABAQUS rentre alors dans le module de VISUALISATION.
3.4.2 Problématique
Pour différentes tailles des éléments : 5, 3 et 1 mm, on se propose de tracer la
déformée de la fibre neutre de la poutre suivant l’axe y et de discuter les résultats
trouvés. A la fin nous allons effectuer une comparaison des résultats numériques
avec les résultats théoriques obtenus par le calcul RDM (solution analytique).
3.4.3 Résultats
Dans cette partie nous présentons les résultats obtenus suite à plusieurs
simulations avec différents paramètres du maillage. • Simulation 1 Taille de la
maille : 5 mm Degré d’interpolation : linéaire Type d’intégration : intégration
réduite Type d’éléments : quadrilatère
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
La valeur maximale de la flèche de la fibre neutre est de -1337 mm (Fig.3.33).
L’erreur relative est de presque -345,64 % par rapport à la solution analytique
(yan = -3,857 mm). Les paramètres choisis pour le maillage de cette première
simulation ne sont pas adéquats et ne conduisent pas à une solution précise. On
doit raffiner le maillage pour ce rapproché de la solution analytique.
• Simulation 2
Dans cette deuxième simulation nous allons diminuer la taille de l’élément tout
en conservant les autres paramètres du maillage. Taille de la maille : 3 mm Degré
d’interpolation : linéaire Type d’intégration : intégration réduite Type d’éléments
: quadrilatère
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
La valeur maximale de la flèche de la fibre neutre est de -5.139 mm (Fig.3.33),
l’erreur relative dans ce deuxième cas est de 32,82 % par rapport à la solution
analytique (yan = -3,857 mm). Pour évaluer la différence entre les valeurs
numérique et analytique, nous allons tracer les courbes théorique et fournie par
ABAQUS. • Traçage de la courbe de la flèche Pour tracer la courbe de 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.3.35) La fenêtre Create
Path s’affiche, choisir Node list et cliquer sur Continue (Fig.3.36) La fenêtre Edit
Edge list Path s’affiche, cliquer sur Add Before (Fig.3.37) Revenir sur le
ViewPort, A partir de l’origine, piquer des nœuds sur toute la longueur et de la
ligne moyenne (Fig.3.38) Cliquer sur Donne puis OK Cliquer sur l’onglet XY
Data Manager, la fenêtre XY Data Manager s’affiche (Fig.3.39) Cliquer sur
Create, choisir Path et cliquer sur Continue La fenêtre XY Data from Path apparait
(Fig.3.40) Pour afficher la courbe cliquer sur Plot (Fig.3.41) Pour afficher les
données Data de la courbe (Fig.3.42), cliquer sur Edit On peut aussi transférer ses
données à l’environnement Excel pour pouvoir comparer les résultats numériques
et analytiques (Fig.3.43).
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
La valeur maximale de la flèche de la fibre neutre est de -5.139 mm (Fig.3.38),
l’erreur relative vaut 0,33 % par rapport à la solution analytique (yan = -3,857
mm). Comme on peut le constater, les courbes numérique (ABAQUS) et
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
analytique (RDM) se rapproche mais la convergence n’est pas atteinte. On doit
raffiner encore le maillage pour ce rapproché de la solution analytique.
• Simulation 3
Dans cette troisième simulation nous avons diminué la taille de l’élément jusqu'à
la valeur 1 mm, tout en conservant les autres paramètres du maillage choisi. Taille
de la maille : 1 mm Degré d’interpolation : linéaire Type d’intégration :
intégration réduite Type d’éléments : quadrilatère
La valeur maximale de la flèche de la fibre neutre est de -4.019 mm (Fig.3.44),
l’erreur relative vaut 0,04 % par rapport à la solution analytique (yan = -3,857
mm). On trace et on superpose les courbes théorique et numérique.
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
• Simulation 4
Au lieu de diminuer la taille du maillage, nous augmentons l’ordre de
d’interpolation, de linéaire à quadratique. Lors du calcul précédant, la
convergence n’a pas été atteinte pour une taille de maille de 5 mm. C’est pour cela
nous allons directement passer à une taille de 3 mm de l’élément. Taille de la
maille : 3 mm Degré d’interpolation : quadratique Type d’intégration : intégration
réduite Type d’éléments : quadrilatère
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
La valeur de la flèche maximale est de -3,857 mm (Fig.3.46), l’erreur relative est
égale à 0 % par rapport à la solution analytique. La solution numérique correspond
parfaitement et coïncide avec la solution analytique (yan = -3,857 mm). Le graphe
de la flèche de la fibre neutre (Fig.3.47) nous montre une parfaite concordance
entre les réponses analytique et numérique, ce qui justifie le bon choix des
paramètres du maillage.
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
• Simulation 5
Quoi que la convergence soit atteinte avec un maillage de 3 mm, nous allons
raffiner encore le maillage. On choisit un maillage de 1 mm et nous obtenons les
résultats suivants : Taille de la maille : 1 mm Degré d’interpolation : quadratique
Type d’intégration : intégration réduite Type d’éléments : quadrilatère
La valeur de la flèche maximale est de -3,858 mm (Fig.3.47), l’erreur relative vaut
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
0,0002 %par rapport à la solution analytique. La solution numérique est en
parfaite concordance avec la solution analytique (yan=-3,857 mm).
• Simulation 6
On se propose ici d’étudié l’influence du type de l’élément et du degré
d’interpolation sur le résultat. On remplace les éléments quadrilatère quadratiques
par des éléments triangulaires linéaires, pour une taille de 3 mm. Taille de la
maille : 3 mm Degré d’interpolation : linéaire Type d’intégration : intégration
réduite Type d’éléments : triangulaires
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
La flèche maximale est de -1,714 mm (Fig.3.50), l’erreur relative vaut 0,55 % par
rapport à la solution analytique. On peut constater la divergence entre la solution
numérique et la solution analytique (yan=-3,857 mm) (Fig.3.51).
• Simulation 7
L’influence du degré d’interpolation sur la précision des résultats est examinée
dans cette simulation. On remplace les éléments triangulaires linéaires par des
éléments triangulaires quadratiques. Taille de la maille : 3 mm Degré
d’interpolation : quadratique Type d’intégration : intégration réduite Type
d’éléments : triangulaires
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
La valeur de la flèche maximale est de -3,39 mm (Fig.3.52), l’erreur relative est
égale à 0,12 % par rapport à la solution analytique. Dans ce cas aussi, on constate
une divergence entre la solution numérique et la solution analytique (yan = -3,857
mm) (Fig.3.53).
• Simulation 8
Pour converger le maillage vers la solution précise on raffine le maillage. On
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
considère un maillage avec les paramètres suivants : Taille de la maille : 1 mm
Degré d’interpolation : linéaire Type d’intégration : intégration réduite Type
d’éléments : triangulaires
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
Pour ce cas, la flèche maximale est de -3,389 mm (Fig.3.54), l’erreur relative est
égale à 0,12 % par rapport à la solution analytique. Dans ce cas aussi, on constate
que la divergence entre la solution numérique et la solution analytique persiste
(yan = -3,857 mm) (Fig.3.55).
• Simulation 9
Pour atteindre la convergence vers une solution plus précise on change le degré
d’interpolation de linéaire à quadratique. Taille de la maille : 1 mm Degré
d’interpolation : quadratique Type d’intégration : intégration réduite Type
d’éléments : triangulaires
Année universitaire 2023-2024 Narjes Adalet
Cours : Simulation numérique par Abaqus Iset Sousse
La flèche maximale est de -3, 838 mm (Fig.3.56). L’erreur relative vaut 0,004 %
par rapport à la solution analytique (yan=-3,857 mm). On constate que la solution
numérique est bonne concordance avec la solution analytique (Fig.3.57).
3.4.5 Conclusion Dans le tableau récapitulatif ci-dessous, nous présentons les
résultats obtenus pour les différents types de maillages employés dans les
simulations. En se référant à la solution analytique, nous avons calculé l’erreur
relative entre la solution numérique (yABQ) et la solution analytique (yan = -
3,857 mm).
Année universitaire 2023-2024 Narjes Adalet
