TP-ABAQUS 1/5
Premiers calculs avec ABAQUS
Objectif : prise en main du logiciel sur un cas simple
Analyse des résultats avec évolution du maillage et du modèle
Prise en main du logiciel
Problème : statique linéaire en contrainte plane
Géométrie de la plaque. 30
p
épaisseur de la plaque est de 5 mm, e=5
pression de résultante 600 N sur 1cm2.
100 r
Premier modèle r = 0 sans congé, vous pourrez
ensuite faire le calcul avec un congé.
20
yo 300
xo
Démarche : Vous devez passer successivement dans les modules (à partir de la liste déroulante):
Part définition et création des géométries
Property Propriétés des matériaux
Assembly Création d’un ensemble de part
Step Type de calcul et sorties souhaitées
Interaction liaisons (contacts) entre les parts
Load Conditions aux limites du modèle
Mesh Réalisation du maillage
Job Paramètres et lancement du calcul
PART
pour créer la géométrie 2D (ouverture de la fenêtre et des outils de dessin)
Utiliser la barre d’outils pour
Dessiner la géométrie (ne pas oublier les éléments utiles pour les CL)
Donner les dimensions
La validation se fait dans le menu du bas
PROPERTY
Création du matériau et de ses propriétés, entrez les caractéristiques d’un acier
Définissez une section solide homogène en acier d’épaisseur 5mm
Associez la section à votre géométrie (domaine)
ASSEMBLY
Création d’instance (assemblage des entités précédentes)
STEP
Vous allez effectuer un calcul statique (pour l’instant on ne précise pas de grandeurs à conserver)
INTERACTION : pour ce modèle il n’y en a pas
LOAD
TP-ABAQUS 2/5
Pour définir les chargements,
ici une pression de 6MPa appliqué sur la surface définie par un segment de 2cm
il faut que le segment ai été défini dans PART
Pour définir vos conditions aux limites en déplacement
Rotule en A et appui simple en B
MESH
Vous devez préalablement sélectionner Part pour mailler votre géométrie (message d’erreur)
Pour définir la taille moyenne des éléments (vous pouvez utiliser 5mm pour un premier calcul)
Pour définir le type d’élément (utilisez une interpolation de degré 1 et l’intégration réduite)
Pour réaliser le maillage de la pièce
JOB : après avoir créé votre Job, il faut le lancer, vous pouvez le monitorer, puis voir les résultats à partir
de cette fenêtre :
Résultats
C’est maintenant que commence votre travail
TP-ABAQUS 3/5
Analyse des résultats
Avant toute analyse :
Il faut savoir ce que vous voulez observez (quelle grandeur ? Comment est-elle calculée ?)
Et comprendre comment est construite l’image que vous observez
Dans Abaqus les possibilités de visualisation sont extrêmement nombreuses, vous apprendrez à les
utiliser au fur et à mesure et cela en fonction des problèmes que vous aurez à traiter.
Pour débuter :
Dans le menu : Result / Field output vous permettra de sélectionner le champ à visualiser
Regardez les déplacements, les contraintes, les réactions, les déformations
Testez différents type de visualisation
Outils utiles
Permet de définir les options de tracé (Nombre d’intervalle et type de contour)
Permet de choisir l’information à visualiser (valeur moyenne ou par élément)
Permet de lister les valeurs avec l’option probe values
En jouant sur les options et ces outils essayez de comprendre ce que vous observez
Outils pratiques
Viewport : Vous pouvez modifier, entre autre, la taille des caractères avec "annotation option".
Tool : Vous pouvez créer un chemin (path par exemple une suite de noeuds) puis visualiser une
grandeur le long de ce chemin (c'est bien pratique).
Premières analyses ( Début du TP : Étude de cas)
Que pensez-vous des résultats
En déplacement ?
Quelle est votre analyse globale ?
Pouvez-vous donner les valeurs maximales de u ? de v ? (Relevez les valeurs)
En contrainte ?
Quelle est votre analyse globale ? (Von Mises ? Tenseur ?)
Quelle est la précision sur la contrainte de Von Mises ? (Relevez les valeurs)
En conclusion la précision sur la contrainte maximale peut-elle être jugée suffisante pour ce modèle?
TP : Étude de cas
L'objectif de ce TP est de comparer les résultats obtenus pour différentes modélisations d'un même
problème.
On se propose de réaliser un modèle 2D en élasticité plane de la pièce sans modéliser le congé.
Les grandeurs d'intérêt. Le comportement général de la pièce
Le déplacement maximal en x et y
La contrainte de Von Mises et les contraintes principales dans la pièce.
Les facteurs d'influence à étudier sont : La prise en compte d'un congé
La modélisation des conditions d'appui
La modélisation du chargement.
TP-ABAQUS 4/5
Savoir analyser un problème est indispensable pour mener une étude numérique dans de bonnes
conditions, c'est la phase la plus complexe de l'étude nous y reviendrons un peu plus tard, il faut savoir
choisir le modèle, définir les objectifs de l'étude et les grandeurs d'intérêt, être capable de prévoir les
facteurs qui influencerons les résultats.
Premier modèle : Nous avons utilisé un maillage grossier avec des éléments de Type Q4 & T3.
Un maillage de ce type donne de mauvais résultats.
Pouvez-vous justifier cette assertion, en regardant et analysant différents résultats ?
Étude de convergence
Modifiez le type d'intégration en utilisant 4 point de gauss. (Dans Mesh modifier le type d'élément)
Cela modifie t'il les résultats ? (Expliquer pourquoi)
Que pensez-vous de la convergence: Globale ? Locale ?
Améliorer ce maillage et étudiez la convergence des résultats.
Dans un premier temps en affinant le maillage de degré 1.
Dans un second temps en utilisant des éléments à base quadratique (degré 2).
Pour étudier cette convergence vous pouvez utiliser les grandeurs d'intérêt.
Pour analyser les résultats chercher les bonnes options de visualisation.
Le fait d'améliorer les résultats numériques du modèle 2D, pose de nouveaux problèmes numériques.
Identifier ces problèmes.
Pouvez-vous donner la valeur de la contrainte maximale dans la pièce ? (avec quelle précision)
Pouvez-vous déterminer la valeur des contraintes au voisinage des appuis ? (à justifier)
Comment peut-on améliorer ce modèle?
Modèle avec congé
La contrainte maximale est située dans l'angle, pour pouvoir la calculer nous devons modéliser un
congé.
Prise en compte d'un congé de 1mm dans l'angle
Revenir dans le module PART, dans l'arbre allez à part puis Features et enfin éditer votre sketch,
Utilisez l'outil Fillet pour créer le congé, donnez lui sa dimension.
Validez le sketch et remailler la structure avant de relancer le calcul.
Que pensez-vous des résultats ?
Affinez votre maillage pouvez-vous conclure sur la convergence du modèle ?
Donnez les valeurs maximales (déplacements et contraintes) avec leur précision
Modifiez le rayon du congé pour introduire un congé numérique (congé < à 1mm)
Donnez les valeurs maximales avec leur précision
Qu'en pensez-vous?
Quel rayon de congé proposez-vous pour un acier de limite élastique Re = 250 MPa
Validez numériquement votre choix de congé.
Faut-il s'y intéresser aux contraintes au voisinage des conditions aux limites ?
Travail à rendre
Une note de calcul synthétisant les principaux résultats et les éléments d'analyse que vous avez
effectuée.
Synthèse
Le cours "Processus de modélisation" s'appuie sur une synthèse des résultats de ce TP.
TP-ABAQUS 5/5
Pour aller un peu plus loin
Influence des conditions aux limites
Vérifier que l'on peut introduire le chargement par une force ponctuelle équivalente.
Qu'elle est l'influence du frottement au niveau des appuis sur les résultats ?
Proposer des modèles permettant d'évaluer cette influence.
Prise en compte de la plasticité du matériau
Pour un congé numérique, et un chargement de 800N une zone plastique existe localement dans le
congé, l’objectif de l’étude est de simuler le comportement élasto-plastique de la pièce et de
déterminer les contraintes et déformations résiduelles après décharge.
Vous utiliserez un acier dont la loi de comportement est définie par la courbe contrainte – déformation
plastique suivante :
σ MPa εp 300
250
207 0 200
210 0,001028 150 Série1
230 0,001763 100
50
250 0,002718 0
270 0,003925 0 0,0005 0,001 0,0015 0,002 0,0025 0,003 0,0035 0,004 0,0045
Modifiez en conséquence le jeu de données
Création des propriétés plastiques (courbe d’écrouissage isotrope)
Création du « Step » de décharge
Analyser les étapes du calcul (processus itératif associé aux non linéarités)
Analyser les résultats.
Refaire le calcul avec une loi de comportement élasto-plastique parfaite (pas d'écrouissage) et un rayon
de congé nul, que pensez-vous de ces résultats ?