Pfe Master Mécanique Blida
Pfe Master Mécanique Blida
Pfe Master Mécanique Blida
Faculté de technologie
Thème
La détermination numérique de ces deux facteurs (Kt et KI) par Abaqus nos a aidé à
distinguer l’effet de la différente forme d’entaille telle que l’entaille en V, l’entaille circulaire
et l’entaille carrée sur le comportement mécanique de la tôle en acier doux.
2
Abstract
1/ The simulation of notched specimens without cracks for the determination of the
stress concentration factor Kt for different specimen shapes with v-notch, circular notch
and square notch.
2/ Simulation of notched specimens with crack for the determination of stress intensity
factor KI for different specimen shapes with v-notch, circular notch and square notch.
Numerical determination of these two factors (Kt and KI) by Abaqus helped us to
distinguish the effect of different notch shape such as V-notch, circular notch and square
notch on the mechanical behavior of mild steel sheet.
Keywords: Abaqus. Stress intensity factor, KI. Stress concentration factor, Kt. Tensile.
Fracture mechanics. Crack.
3
ملخص
تم تخصيص هذا العمل لمحاكاة عينات الشد (وفقًا للمعيار األوروبي )EN 10002-1: 1990مع الشق
،Vالدائري والمربع ،من ألواح الفوالذ الطري بسمك 2مم ،واستخدام المحاكاة له أهمية متزايدة ،ألن االختبارات
باهظة الثمن .محور هذا العمل موجه نحو المحاكاة باستخدام برنامج Abaqusفي جزئين متميزين:
/1محاكاة العينات المحززة بدون تشققات لتحديد عامل تركيز اإلجهاد Ktألشكال عينات مختلفة ذات
الشق Vوالشق الدائري والشق المربع.
/ 2محاكاة العينات المحززة مع التشققات لتحديد عامل شدة اإلجهاد KIألشكال العينات المختلفة ذات الشق
Vوالشق الدائري والشق المربع.
ساعدنا التحديد العددي لهذين العاملين ( Ktو )KIبواسطة Abaqusعلى التمييز بين تأثير أشكال الشق
المختلفة مثل الشق Vوالشق الدائري والشق المربع على السلوك الميكانيكي لصفائح الفوالذ الطري.
كلمات مفتاحية .Abaqus :عامل شدة اإلجهاد .KIعامل تركيز االجهاد .Ktميكانيكا الكسر .كسر.
4
DEDICACES
Dédicaces
A nos collèges.
5
REMERCIEMENTS
Selon la tradition, et à travers cette page de remerciements, nous tenons à saluer toutes
les personnes qui, de près ou de loin, ont contribué à la concrétisation de ce travail.
Nous remercions Monsieur DAOUD Lotfi et Madame RAMDOUM Sara pour avoir
accepté de diriger ce travail, et leurs conseils, leurs observations et leurs aides pour la
réalisation de ce mémoire.
Nous voudrions remercier les membres du jury pour leur précieux temps accordé à
la lecture et à la critique de ce travail.
Nos très cordiaux remerciements vont aussi à tous nos enseignants pour tous les
efforts conjugués tout au long de notre formation.
6
LISTE DES FIGURES
7
Figure 40 : Résultats de la simulation. ..................................................................................... 53
Figure 41 : Géométrie avec l’emplacement de la fissure. ........................................................ 54
Figure 42 : section de fissure.................................................................................................... 54
Figure 43 : Création de la fissure. ............................................................................................ 55
Figure 44 : Pointe de la fissure. ................................................................................................ 55
Figure 45 : Direction de propagation. ...................................................................................... 56
Figure 46 : Singularité. ............................................................................................................. 56
Figure 47 : Paramètres de sortie(1). ......................................................................................... 57
Figure 48 : Paramètres de sortie (2). ........................................................................................ 58
Figure 49 : Création de job (2). ................................................................................................ 58
Figure 50 : Analyse (2). ............................................................................................................ 59
Figure 51 : Résultat de propagation de la fissure. .................................................................... 59
Figure 52 : Kt en fonction de b/a pour l’entaille Circulaire. .................................................... 64
Figure 53 : Kt en fonction de b/a pour l’entaille en V. ............................................................. 65
Figure 54 : Kt en fonction de b/a pour l’entaille carrée. ........................................................... 66
Figure 55 : Comparaison entre les facteurs de concentration de contraintes Kt des éprouvettes
de différents forme d’entaille (en V, circulaire et carrée). ....................................................... 66
Figure 56 : Variation de facteur d’intensité de contrainte (KI) en fonction ............................. 68
Figure 57 : Variation de facteur d’intensité de contrainte (KI) en fonction ............................. 69
Figure 58 : Variation de facteur d’intensité de contrainte (KI) en fonction ............................. 70
Figure 59 : Entaille en V (b=3mm) .......................................................................................... 77
Figure 60 : Entaille V (b=6mm) ............................................................................................... 78
Figure 61 : Entaille en V (b=9mm) .......................................................................................... 79
Figure 62 : Entaille V (b=12mm) ............................................................................................. 80
Figure 63 : Entaille en V (b=15mm) ........................................................................................ 81
Figure 64 : Entaille Circulaire (R = 3mm) ............................................................................... 82
Figure 65 : Entaille Circulaire (R = 6mm) ............................................................................... 83
Figure 66 : Entaille Circulaire (R = 9 mm) .............................................................................. 84
Figure 67 : Entaille Circulaire (R = 12 mm) ............................................................................ 85
Figure 68 : Entaille Carrée (b = 3 mm) .................................................................................... 86
Figure 69 : Entaille Carrée (b = 6 mm) .................................................................................... 87
Figure 70 : Entaille Carrée (b = 9 mm) .................................................................................... 88
Figure 71 : Entaille Carré (b = 12 mm) .................................................................................... 89
Figure 72 : Entaille Carrée (b = 15 mm) .................................................................................. 90
8
LISTE DES TABLEAUX
9
NOMENCLATURE
a : profondeur de l’entaille.
b : largeur de l’entaille.
L : longueur de fissure.
G : paramètre qui peut s’exprimer en J/m2 ou en N/m, correspondant à une énergie libérée par
unité de surface, parfois appelée force d’extension de la fissure.
Re : limite élastique.
M : moment fléchissant.
T : couple de torsion.
𝝉 : contrainte tangentielle.
10
TABLE DES MATIERES
Résumé ................................................................................................................................................... 2
Abstract .................................................................................................................................................. 3
ملخص........................................................................................................................................................ 4
DEDICACES ......................................................................................................................................... 5
REMERCIEMENTS ............................................................................................................................. 6
LISTE DES FIGURES .......................................................................................................................... 7
LISTE DES TABLEAUX ..................................................................................................................... 9
NOMENCLATURE ............................................................................................................................ 10
INTRODUCTION ............................................................................................................................... 13
Chapitre I ............................................................................................................................................. 15
Chapitre I : Notions de base sur La Mécanique de la rupture. ....................................................... 16
I. 1 Introduction............................................................................................................................... 16
I. 2 Aperçu historique sur la mécanique de la rupture ................................................................ 17
I. 3 Objectifs de la mécanique de la rupture ................................................................................. 17
I. 4 Principes de la mécanique de la rupture................................................................................. 18
I. 5 Concept d’intensité des contraintes......................................................................................... 18
I. 6 Modes de rupture ...................................................................................................................... 19
I. 6.1 Champs des contraintes et des déplacements au voisinage de l’extrémité de fissure .. 20
I. 7 Principe de superposition en mécanique de la rupture ......................................................... 22
I. 8 Facteur critique d'intensité de contrainte............................................................................... 23
I. 9 Énergie de rupture .................................................................................................................... 23
I. 9 .1 Théorie de Griffith (paramètre G) .................................................................................. 23
I. 9. 2 Application à un cas simple .............................................................................................. 24
I. 10 Relations entre G et K ............................................................................................................ 25
Conclusion ............................................................................................................................................ 26
Chapitre II ........................................................................................................................................... 27
Chapitre II : Concentration de contraintes....................................................................................... 28
II. 1 Introduction ............................................................................................................................. 28
II. 2 Facteur de concentration de contraintes ............................................................................... 29
II. 3 Méthodes utilisées pour connaître les facteurs de concentration de contraintes ............... 31
II. 3. 1 Théorie de l’élasticité ...................................................................................................... 31
11
II. 3. 2 Méthodes expérimentales................................................................................................ 31
II. 4 Moyens de diminuer les effets de concentration de contraintes .......................................... 32
II. 5 Facteurs de concentration de contraintes en phase de déformations plastiques ............... 32
Conclusion ............................................................................................................................................ 35
Chapitre III .......................................................................................................................................... 36
Chapitre III : Simulation avec Abaqus ............................................................................................. 37
III. 1 Introduction............................................................................................................................ 37
III. 2 Le logiciel Abaqus : ............................................................................................................... 37
III. 2. 1 L’interface d’ABAQUS ................................................................................................. 39
III. 3 Modèle géométrique et mécanique de la simulation ........................................................... 40
III. 4 Modélisation avec Abaqus : .................................................................................................. 43
III. 4. 1 Définition de la géométrie : ........................................................................................... 43
III. 4. 2 Définition du matériau : ................................................................................................ 45
III. 4. 3 Création et assignation d’une section : ........................................................................ 46
III. 4. 4 Création d’un assemblage : ........................................................................................... 47
III. 4. 5 Définition des Steps : ..................................................................................................... 48
III. 4. 6 Définir les conditions aux limites et le chargement : .................................................. 49
III. 4. 7 Création du maillage : ................................................................................................... 51
III. 4. 8 Définir et lancer un job : ............................................................................................... 52
III. 4. 9 Création de la fissure :................................................................................................... 54
III. 5 Conclusion : ............................................................................................................................ 60
Chapitre IV .......................................................................................................................................... 61
Chapitre IV : Principaux résultats de simulation avec Abaqus ...................................................... 62
IV. 1 Introduction ............................................................................................................................ 62
IV. 2 Facteur de concentration de contraintes (Kt) : .................................................................... 63
IV.2 Facteur d’intensité de contraintes (KI) : ............................................................................... 67
CONCLUSION .................................................................................................................................... 71
REFERENCES .................................................................................................................................... 74
ANNEXES ............................................................................................................................................ 76
Quelques géométrie et iso-valeurs des contraintes obtenues par Abaqus .................................. 77
12
INTRODUCTION
L'approche de résistance des matériaux, dans laquelle les sollicitations appliquées sont
seulement comparées à la limite d'élasticité ou à la résistance à la traction du matériau sous
contrainte, mais l’approche mécanique de la rupture, où on s'intéresse non seulement à la
contrainte appliquée, mais également à l'occurrence et à la taille de fissures présentes dans le
matériau ainsi qu'à la ténacité du matériau, ce travail met en évidence une approche qui se
repose les résultats de la simulation sur le logiciel des éléments finis Abaqus que nous allons
développer dans notre mémoire intitulé : effet de l’entaille sur le comportement mécanique
d’une plaque en traction.
Nous allons étudier cet effet par le calcul de deux facteurs qui sont très utilisé pour la
construction d’une structure mécanique, ce sont le facteur de concentration de contrainte (Kt)
et le facteur d’intensité de contraintes en mode I (KI), sur des éprouvettes de traction (selon la
norme européenne EN 10002-1 :1990) à entaille en V, circulaire et carré, d’une tôle en acier
doux d’épaisseur 2mm.
-Le premier chapitre est dédié aux quelques notions de base sur la mécanique de la rupture, le
premier volet donne : aperçu, objectifs, principes de la mécanique de la rupture et catégorie et
mode de rupture ; le deuxième volet est consacré au facteur d’intensité de contraintes et énergie
de rupture.
- Le deuxième chapitre est dédié aux notions de base sur la concentration de contraintes, le
premier volet donne : une définition au facteur de concentration de contraintes et ces méthodes
de calcul ; le deuxième volet est consacré aux moyens de diminuer les effets de concentration
de contraintes et aux facteurs de concentration de contraintes en phase de déformations
plastiques.
- Le troisième chapitre est dédié à la simulation avec le logiciel Abaqus, le premier volet sur
Abaqus, le deuxième volet sur le modèle géométrique et mécanique à simuler, le troisième volet
donne les étapes à suivre faire faire la simulation sur notre modèle.
13
- Le quatrième chapitre est dédié aux principaux résultats de simulation avec Abaqus, le premier
volet les résultats de facteur de concentration de contraintes ; et le deuxième volet les résultats
de facteur d’intensité de contraintes.
Enfin la conclusion réitère les points les plus importants de l’étude ci présente suivie d'une liste
de bibliographie et les annexes clôturent ce mémoire.
14
Chapitre I :
Notions de base sur
la mécanique de la rupture
15
Chapitre I : Notions de base sur La Mécanique de la rupture.
I. 1 Introduction
Ce chapitre est une initiation à la mécanique de la rupture, et on se limite seulement à des
notions de base pour connaitre le phénomène de la rupture des matériaux en générale ; et nous
soulignerons les titres suivants :
•Catégories de rupture.
•Modes de rupture.
-Champs des contraintes et des déplacements au voisinage de l’extrémité de fissure.
•Énergie de rupture.
-Théorie de Griffith (paramètre G).
-Relations entre G et K.
16
I. 2 Aperçu historique sur la mécanique de la rupture
Théoriquement, la contrainte à la pointe de fissure devient infinie et ne peut pas être utilisée
pour décrire son état. La mécanique de la rupture est utilisée pour caractériser les charges sur
une fissure, en utilisant généralement un seul paramètre pour le décrire. Un certain nombre de
paramètres différents ont été développés. Lorsque la zone plastique à l'extrémité de la fissure
est petite par rapport à la longueur de la fissure, l'état de contrainte à l'extrémité de la fissure est
le résultat de forces élastiques à l'intérieur du matériau et est appelé mécanique de rupture
élastique linéaire (MREL) et peut être caractérisé à l'aide de facteur d'intensité de contrainte K.
Bien que la charge sur une fissure puisse être arbitraire, en 1957, G. Irwin a découvert que
n'importe quel état pouvait être réduit à une combinaison de trois facteurs d'intensité de
contrainte indépendants :
• Quelques dates :
17
La mécanique de la rupture devrait tenter de fournir des réponses quantitatives aux questions
suivantes :
2/Combien de temps faut-il pour qu'une fissure se développe d'une certaine taille initiale, par
exemple la taille de fissure minimale détectable, à la taille de fissure maximale autorisée ?
3/Quelle est la durée de vie d'une structure lorsqu'une certaine taille de défaut préexistant (par
exemple un défaut de fabrication) est supposée exister ?
4/Pendant la période disponible pour la détection des fissures, à quelle fréquence la structure
doit-elle être inspectée pour les fissures ? [3]
18
𝑲
𝝈𝒙 = 𝒇𝒙 (𝜽)…………………………………………………………...…………..……(1)
√𝟐 𝝅 𝒓
𝑲
𝝈𝒚 = 𝒇𝒚 (𝜽) ………………………………………………………..……..……………(2)
√𝟐 𝝅 𝒓
𝑲
𝝉𝒙 = 𝒇𝒙𝒚 (𝜽) …………………………………………………………………...………(3)
√𝟐 𝝅 𝒓
𝑲 = 𝜶 𝝈 √ 𝝅 𝒂 …………………………………………..…………………………......……(4)
•Mode I :
Figure 4 : Mode I.
𝜽 𝟑𝜽
𝟏 − 𝒔𝒊𝒏 (𝟐) 𝒔𝒊𝒏 ( 𝟐 )
𝝈𝒙𝒙
𝑲 𝜽 𝜽 𝟑𝜽
[ 𝝉𝒙𝒚 ] = 𝝅𝑰 𝒓 𝒄𝒐𝒔 (𝟐) 𝒔𝒊𝒏 (𝟐) 𝒄𝒐𝒔 ( 𝟐 ) ………………….…………………………(6)
√𝟐
𝝈𝒚𝒚 𝜽 𝟑𝜽
[𝟏 + 𝒔𝒊𝒏 (𝟐) 𝒔𝒊𝒏 ( 𝟐 )]
𝜽 𝜽
𝒖𝒙 𝑲𝑰 𝒓
𝒄𝒐𝒔 (𝟐) (𝝌 − 𝟏 + 𝟐 𝒔𝒊𝒏𝟐 (𝟐))
[𝒖 ] = 𝟐𝑮 √𝟐 𝝅 [ 𝜽 𝜽
]……………………….…………………(7)
𝒚
𝒔𝒊𝒏 (𝟐) (𝝌 + 𝟏 − 𝟐 𝒄𝒐𝒔𝟐 (𝟐))
20
•Mode II :
𝜽 𝜽 𝟑𝜽
− 𝒔𝒊𝒏 (𝟐) (𝟐 + 𝒄𝒐𝒔 (𝟐) 𝒄𝒐𝒔 ( 𝟐 ))
𝝈𝒙𝒙
𝑲 𝜽 𝜽 𝜽 𝟑𝜽
[ 𝝉𝒙𝒚 ] = 𝑰𝑰𝝅 𝒓 𝒄𝒐𝒔 (𝟐) 𝒄𝒐𝒔 (𝟐) (𝟏 − 𝒔𝒊𝒏 (𝟐) 𝒔𝒊𝒏 ( 𝟐 )) ……………………………… (8)
√𝟐
𝝈𝒚𝒚 𝜽 𝜽 𝟑𝜽
[ 𝒔𝒊𝒏 ( ) 𝒄𝒐𝒔 ( ) 𝒄𝒐𝒔 ( ) ]
𝟐 𝟐 𝟐
𝜽 𝜽
𝒖𝒙 𝑲𝑰𝑰 𝒓
𝒔𝒊𝒏 (𝟐) (𝝌 + 𝟏 + 𝟐 𝒄𝒐𝒔𝟐 (𝟐))
[𝒖 ] = 𝟐𝑮 √𝟐 𝝅 [ 𝜽 𝜽
]…………………….…………………(9)
𝒚
−𝒄𝒐𝒔 (𝟐) (𝝌 − 𝟏 − 𝟐 𝒔𝒊𝒏𝟐 (𝟐))
Avec :
E : Module d’élasticité longitudinal du matériau (Module d’Young).
G : Module d’élasticité transversal ou de cisaillement du matériau (Module de Coulomb).
𝐸
Où : 𝐺 = et υ : coefficient de Poisson du matériau.
2(1−𝜐)
21
•Mode III :
𝜽
𝝉𝒙𝒛 𝑲
−𝒔𝒊𝒏 (𝟐)
[𝝉 ] = 𝑰𝑰𝑰 [ 𝜽
]…………………………………………………..………………(10)
𝒚𝒛 √𝟐 𝝅𝒓
𝒄𝒐𝒔 (𝟐)
𝒖𝒛 =
𝑲𝑰𝑰𝑰 𝒓 𝜽
𝟐 √𝟐 𝝅 𝒔𝒊𝒏 (𝟐)………………………………………...………………………...………(11)
𝑮
Avec :
σzz = 0 ;en contrainte plane.
σzz = υ(σxx – σyy) ;en déformation plane.
Et avec aussi :
E : Module d’élasticité longitudinal du matériau (Module d’Young).
G : Module d’élasticité transversal ou de cisaillement du matériau (Module de Coulomb).
𝐸
Où : 𝐺 = et υ : coefficient de Poisson du matériau.
2(1−𝜐)
22
- On ne peut additionner des FIC (Facteur d’intensité de contrainte) que s’ils concernent le
même mode de sollicitation (mode I, II ou III).
Exemple :
Si une fissure est sollicitée en mode I par trois chargements différents (σ1 ,σ2 ,σ3) alors on a:
KItot = KI1 + KI2 + KI3 ……………………………………………...……………...………(12)
• Si la même fissure est sollicitée simultanément en mode I, en mode II et en mode III, dans
ce cas : Ktot ≠ KI + KII+ KIII ………………………….……………………………………(13)
I. 8 Facteur critique d'intensité de contrainte
Une rupture se produit lorsque la contrainte appliquée est supérieure à une certaine valeur
critique σc , alors on a :
𝑲𝑪 = 𝜶 𝝈𝒄 √𝝅 𝒂 ………………………………….………….…………….………………(14)
Cette valeur critique de Kc , appelée parfois ténacité, porte le nom de facteur critique d'intensité
de contrainte.
• Irwin introduit la notion de ténacité d’un matériau en postulant qu’il existe une valeur critique
KIC du facteur d’intensité des contraintes KI du mode d’ouverture qui ne peut jamais être
dépassée et en deçà de laquelle la fissure ne peut pas se propager. [1]
Matériau KIC(MPa.m1/2)
- Alliage d’Aluminium 30
- Alliage de Titane 100
- Acier trempé 120
- Polymère 3
- Bois 2
- Béton 1
Tableau 1 : Ténacité de quelques matériaux courants.
I. 9 Énergie de rupture
I. 9 .1 Théorie de Griffith (paramètre G)
Si nous supposons que les forces extérieures Fe dérivent d’un potentiel V, un accroissement
virtuel δA d’aire fissurée libère une énergie G δA telle que :
23
𝜹𝑷 = − 𝑮 𝜹𝑨 ………………………………….………….………………..………………(15)
Et : P=W+V =W-Te……………………………..………….…………….……………..…(16)
Avec :
G : paramètre qui peut s’exprimer en J/m2 ou en N/m, correspondant à une énergie libérée par
unité de surface, parfois appelée force d’extension de la fissure.
Il existe un seuil critique Gc au-delà duquel une extension de la fissure libère plus d’énergie
qu’elle n’en absorbe, et la fissure est instable pour G tel que : G δA > Gc δA.
La figure suivante représente une éprouvette d’épaisseur constante e, contenant une fissure de
longueur a, soumise à une force constante F dont le point d’application se déplace suivant une
loi linéaire : z = C(a) F avec : C (a) compliance.
À un accroissement virtuel δa de la fissure, correspond :
1
δW = δ (2 𝐹𝑧) : Variation de l′énergie de déformation élastique.
24
𝐹2 𝑑𝐶
Dans ce cas, l’énergie libérée par unité de surface G vaut : 𝐺 = 2 𝑒 𝑑𝑎
C’est une relation très générale, valable en théorie linéaire pour tous les chargements et
déplacements.
I. 10 Relations entre G et K
Le facteur d’intensité de contraintes K est relié au facteur G (et ce quel que soit le mode I, II,III)
par les relations :
𝐾2
𝐺= (1 − 𝜐 2 ) en déformations planes (DP) ………………………….……………..…(17)
𝐸
𝐾2
𝐺= en contraintes planes (CP) ………………………………………....…………..…(18)
𝐸
Avec :
25
Conclusion
26
Chapitre II :
Concentration de contraintes
27
Chapitre II : Concentration de contraintes
II. 1 Introduction
Un état de contrainte uniforme est généralement illustré par une colonne droite soumise
à une charge uniaxiale à condition que cette charge soit également appliquée de manière
uniforme aux extrémités de la colonne. Un cas comme celui-ci est rarement, voire jamais trouvé
dans la pratique. En général, l'état de contrainte change d'un point à l'autre et si l'une des
composantes de contrainte (ou les résultantes) est représentée par une famille de lieux, des pics
et des creux de valeurs apparaissent. La contrainte est donc répartie de manière non uniforme
et la notion pas toujours précisément définie de concentration de contrainte est associée à des
pics, et des gradients autour des pics, des composantes de contrainte.[5]
Le calcul des efforts internes par la résistance des matériaux pour des exemples usuelles :
1/ Contrainte normale induite par traction ou compression pure :
𝑵
𝝈= ………………………………….……………………….……….…………………(19)
𝑨
28
1/ Contrainte normale induite par flexion pure :
𝑴𝒛
𝝈(𝒛) = ………………………………….…..…………….….…………………………(20)
𝑰
Où I est le moment d’inertie par rapport à un axe orthogonale au plan moyen de la poutre, z la
côte comptée à partir de la fibre neutre et M le moment fléchissant.
𝝉
𝑲𝒕𝒔 = en torsion………………….….…..………...…..…….…………………………(23)
𝝉𝟎
Figure 12 Barre de tension avec encoches : (a) éprouvette ; (b) photographie de frange
photoélastique (Doz. Dr.-Ing. habil. K. Fethke, Universität Rostock); (c) solution par éléments
finis (Guy Nerad, Université de Virginie) [6]
30
D’après les formulaires qui donnent le facteur de concentration de contraintes, on a :
𝑟
= 0,15
𝑑
𝐷 } ⇒ 𝐾𝑡𝑠 = 1,35
= 1,25
𝑑
16 𝑇
D’où : 𝜏 = 𝐾𝑡𝑠 𝜏0 = 𝐾𝑡𝑠 𝜋 𝑑3
La solution analytique des équations de l’élasticité dans le cas d’une plaque à grande
dimension avec un trou elliptique chargé en traction nous donne le facteur de concentration de
contraintes. (voir la figure suivante)
Remarque : si les contraintes seront compliquées il est impossible de résoudre les équations
de l’élasticité analytiquement et trouver le facteur de concentration de contraintes
Remarquons que dès que la géométrie ou bien les charges deviennent compliquées, il est
impossible d’exhiber analytiquement les facteurs de concentration
II. 3. 2 Méthodes expérimentales
31
- la technique de vernis craquelant;
- la technique de la grille incluse;
- la technique des franges de Moiré;
- les techniques analogiques.
Remarque : La technique de l’extensomètrie par les jauges de déformation est la plus utilisé.
Les moyens les plus utilisés pour diminuer les effets de concentration des contraintes sont les
suivantes :
- il faut que les sections ne varient pas brusquement et que les rayons de raccordement soient
les plus grands possibles.
-les contraintes appliquées soient dans le sens de l’axe du grand trou elliptique.
- utilisant des arrondis et en enlèvement plus du matériau pour réduire l’effet de concentration
de contraintes.
II. 5 Facteurs de concentration de contraintes en phase de déformations plastiques
32
1/ En élasticité pur (σ < σy): avec σy c’est la limite élastique du matériau.
-Le facteur de concentration de contraintes est : K1= σ / σ0
-pour d/w= 0.5, on obtient K1 =2.8.
-Les contraintes sont représentées sur la figure suivante.
33
Remarques :
1/ Matériau ductile :
Si on décharge la section en domaine plastique, alors on a la naissance des contraintes
résiduelles dans cette section, et si on l'a chargé un autrefois, on a les contraintes
représentées sur la figure suivante.
34
Conclusion :
Dans ce chapitre on a donné les notions de base (telle quelle sont dans la bibliographie)
sur la concentration de contraintes, le premier volet donne : une définition au facteur de
concentration de contraintes et ces méthodes de calcul ; le deuxième volet est consacré aux
moyens de diminuer les effets de concentration de contraintes et aux facteurs de concentration
de contraintes en phase de déformations plastiques.
de la tôle en traction.
35
Chapitre III :
Simulation avec Abaqus
36
Chapitre III : Simulation avec Abaqus
III. 1 Introduction
Dans ce chapitre nous avons présenté la procédure de modélisation avec Abaqus d’une
plaque en forme d’éprouvette avec entaille en v, entaille circulaire et entaille carré, ces entailles
ont la même profondeur (a = 3mm) mais des largeurs des entailles variables (3mm – 6mm –
9mm –12mm –15mm), et le modèle géométrique et mécanique utilisé avec le mode de
chargement et les conditions aux limites dans ce cas.
Abaqus est un progiciel de calcul d'éléments finis qui a été développé par Hibbit,
Karlsson & Sorensen (HKS) devenu ABAQUS, Inc depuis 30 ans en début de la théorie des
éléments finis et a été amélioré au fur et à mesure pour intégrer toutes les nouveautés de la
théorie et des besoins de l’entreprise, jusqu’à ce que l’entreprise soit rachetée par Dassault
industries en Octobre de l’année 2005. ABAQUS est avant tout un logiciel de simulation par
éléments finis de problèmes très variés en mécanique. Il est connu et répandu, en particulier
pour ses traitements performants de problèmes non-linéaires.
ABAQUS est donc un progiciel de calcul éléments finis qui se compose de trois produits :
37
• ABAQUS/CAE : constitue une interface intégrée de visualisation et de modélisation pour les
dits solveurs.
Chacun de ces produits est complété par des modules additionnels et/ou optionnels, spécifiques
à certaines applications. Les produits ABAQUS, ABAQUS/CAE notamment, sont écrits
intégralement avec les langages C++, Fortran pour les parties calcul et Python pour les scripts
et les paramétrisations. La gestion de l'interface graphique est assurée par fox-toolkit. [8]
ABAQUS est très largement utilisé dans les industries automobiles et aéronautiques. En raison
du large spectre de ses capacités d'analyse et de sa bonne ergonomie, il est également très
populaire dans les milieux universitaires, pour la recherche et l'éducation.
ABAQUS fut d'abord conçu pour analyser les comportements non-linéaires. Il possède en
conséquence une vaste gamme de modèles de matériau. Ses modélisations d'élastomères, en
particulier, méritent d'être reconnues.
En 2005, ABAQUS, Inc a été acquise par Dassault Systèmes et est devenu SIMULIA Le coeur
du logiciel ABAQUS est donc ce qu’on pourrait appeler son "moteur de calcul". À partir d’un
fichier de données (caractérisé par le suffixe .inp), qui décrit l’ensemble du problème
mécanique, le logiciel analyse les données, effectue les simulations demandées et fournit les
résultats dans un fichier .odb.
Deux tâches restent à accomplir : générer le fichier de données (cela s’appelle aussi effectuer
le prétraitement), et exploiter les résultats contenus dans le fichier .odb (ou post traitement).
La structure du fichier de données peut se révéler rapidement complexe : elle doit contenir
toutes les définitions géométriques, les descriptions les maillages, les matériaux, les
chargements, etc…, suivant une syntaxe précise. [9]
Il faut savoir que le prétraitement et le post traitement peuvent être effectués par d’autres
logiciels.
ABAQUS propose le module ABAQUS CAE, interface graphique qui permet de gérer
l’ensemble des opérations liées à la modélisation :
38
- L’exploitation des résultats. [10], [11].
- une première barre de menus en haut : ces menus sont relatifs au module dans lequel on se
trouve.
- une deuxième barre horizontale d’outils : y sont disposées les icônes correspondant aux
commandes les plus courantes, c’est-à-dire les commandes d’ouverture et de sauvegarde de
fichiers, et les commandes réglant l’affichage des vues (perspective, ombrage, zoom, rotations,
etc.)
- une troisième barre de menus déroulant permettant d’accéder aux autres modules, ou de
sélectionner le modèle ou la pièce sur lesquels on souhaite travailler.
Plus important :
- on dispose à gauche d’une colonne d’icônes permettant d’accéder aux outils disponibles à
l’intérieur du module dans lequel on se trouve. ATTENTION : les icônes comportant une petite
flèche en bas à droite ( ) cachent d’autres icônes ! Il faut maintenir enfoncée la touche M1 de la
souris sur l’cône pour voir apparaître les icônes cachées et pouvoir sélectionner les outils
correspondants. Présentation du logiciel. L’espace juste sous la fenêtre d’affichage graphique
est celui dans lequel ABAQUS vous parle : les messages affichés à cet endroit sont faits pour
vous guider dans l’action que vous avez entreprise. Lisez-les, suivez-les attentivement, ils vous
seront d’une grande aide.
- Enfin, une fenêtre en-dessous sur fond blanc est utilisée pour afficher des messages informatifs
: réponses à des questions posées, avertissements liés à des problèmes rencontrés, etc. [12] .
39
Figure 19 : L'interface Abaqus.
III. 3 Modèle géométrique et mécanique de la simulation
Notre étude dans ce mémoire consiste à simuler une tôle d’acier chargé en traction avec
trois types d’entaille, et pour faciliter cette tache nous allons prendre une éprouvette de traction
plate normalisé comme modèle de base, car dans l’essai de traction on utilise pour les tôles des
éprouvettes de section rectangulaire normalisés selon la norme européenne EN 10002-1 :1990.
[13] comme suit:
e = 2 mm
b0 = 40 mm
40
B = 50 mm
R = 12 mm
𝑳𝟎 = 𝒌 √𝑺 = 𝟓. 𝟔𝟓 √𝒃𝟎 . 𝒆 ……………………………………………………………….(25)
A.N :
On prend : L0 = 50 mm
L0 + b0 ≤ Lc ≤ L0 + 2 b0 ……………………………………………………………...….(26)
A.N :
On prend : Lc = 60 mm
L Lc + 2 à 4 b0 ……………………………………………………………………......….(27)
A.N :
L 60 + 2 à (4 x 40) ⇒ L 60 + 2 à 160
On prend : L = 200 mm
On prend aussi : h = 65 mm
Les dimensions de notre modèle géométrique est illustré dans le tableau suivant :
41
Les entailles utilisées sont comme suit :
• Entaille en V.
• Entaille en cercle.
• Entaille en carré.
1/ Géométrie de l’entaille en V :
42
Figure 23: Dimensions de l’entaille carré avec fissure(L).
En résume dans le tableau suivant les Dimensions de l’entaille en V, en cercle et carré.
*La longueur de fissure (L) prend les valeurs : 4mm – 7mm – 10mm –15mm –20mm –25mm.
III. 4 Modélisation avec Abaqus :
F
Figure 26 : Création de la géométrie 3D.
44
Figure 27: La géométrie 3D dans Abaqus.
45
Figure 28 : Définition du matériau.
46
Figure 29 : Définition de la section.
47
Figure 31 : Création de l’assemblage.
48
Figure 32 : Création de l’étape « step-1 ».
49
Figure 33 : Définition de la charge.
50
III. 4. 7 Création du maillage :
2. Sélectionnez Mesh part instance , vous allez choisir la taille approximative des
éléments ici 3.3 puis OK.
51
4. Selection Mech-Part-Instance et Cliquer yes.
Figure 39 : Analyse(1).
3. Quand Job est marqué Completed, c’est que le travail est fini est qu’il est alors
possible de visualiser les résultats. Cliquez sur results.
1. Retour à la section Module, sélectionner et cliquer sur Part puis Cliquer sur
3. Cliquer sur l’icône Create Circle : Center and Perimeter pour créer un
cercle autour de la pointe de fissure nécessaire pour la création des contours et
Cliquer sur Done.
55
Figure 45 : Direction de propagation.
8. Dans la nouvelle fenêtre Edit Crack.et dans l’onglet singularity, on va définir les
paramètres de la singularité on met 0.5 dans Midside node parameter (0<t<1), on
choisit collapsed element side, single node et on termine par ok.
Figure 46 : Singularité.
56
9. Entrer dans le module Mesh ,Sélectionnez Seed – Edge By Number et spécifiez
que 8 éléments seront créés au nouveau de cercle et cliquez sur ok. La même chose
pour le rayons de cercle, entrer 4 et cliquez sur ok.
Maintenant :
3. Dans la nouvelle fenêtre, Dans cette boîte, choisissez Crack dans Domain. Et last
increment dans Frequency pour le nombre de contours on met 5.
4. Ensuite, dans type : choisissez Maximum tangential stress.
5. Et enfin validez avec OK.
57
Figure 48 : Paramètres de sortie (2).
6. Entrer dans le module JOB et créez un nouveau job appelé « fissure-v3»,la meme
methode que nous avons utilisee dans le cas sans fissure precedent.
58
7. Sauvegardez votre modèle et lancez (submit) le job pour une analyse.
Affichez le Monitor pour voir la progression de votre analyse.
59
III. 5 Conclusion :
60
Chapitre IV :
Principaux résultats de
simulation avec Abaqus
61
Chapitre IV : Principaux résultats de simulation avec Abaqus
IV. 1 Introduction
Dans ce chapitre on va présenter les différents résultats obtenus par la simulation avec
Abaqus qui sont partagés en deux parties :
62
IV. 2 Facteur de concentration de contraintes (Kt) :
Nous présentons dans ce qui suit la variation des facteurs de concentration de
contraintes Kt en fonction de la variation de la géométrie de l’entaille (b/a) pour différent forme
d’éprouvette avec entaille en v, circulaire et carré.
Avec :
σ0 : Contrainte maximale appliqué au bord de l’éprouvette lors de l’essai de traction (on prend
la contrainte maximale de la courbe de traction de l’éprouvette sans entaille égale à 300Mpa).
1/Entaille Circulaire :
63
3,5
Kt = f(b/a) pour entaille circulaire
3
2,5
2
Kt
1,5
0
0 0,5 1 1,5 2 2,5 3 3,5 4 4,5 5
b/a
Kt = f(b/a) entaille V
2,65
2,6
Kt
2,55
2,5
Kt = -0.0013(b/a) 4 + 0.0108(b/a) 3 - 0.0337(b/a) 2 + 0.0142(b/a) + 2.69
2,45 R² = 1
2,4
0 0,5 1 1,5 2 2,5 3 3,5 4 4,5 5
b/a
2/ Entaille carré :
65
3,4 Kt = f(b/a) entaille carré
3,2
3
2,8
Kt
2,6
2,4
2,2 Kt = 0.0108(b/a)4 - 0.14(b/a)3 + 0.6642(b/a)2 - 1.605(b/a) + 4.42
R² = 1
2
0 1 2 3 4 5
b/a
Series2 Series3
3,25
2,75
Kt
2,5
2,25
1,75
1,5
0 1 2 3 4 5
b/a
66
IV.2 Facteur d’intensité de contraintes (KI) :
Nous présentons dans ce qui suit la variation des facteurs d’intensité de contraintes (KI)
en fonction de la variation de la propagation de longueur de fissure pour différent forme
d’éprouvette avec entaille en v, circulaire et carrée.
1/ Entaille Circulaire :
Largeur
Facteur d’intensité de contraintes (Mpa.mm0.5)
de l’entaille (b)
↓ ↓ ↓ ↓ ↓ ↓ ↓
3mm 2350 3215 5012 10312 15033 16908
67
18000
Entaille 3 mm
16000 Entaille 6 mm
Entaille 9 mm
14000 Entaille 12 mm
Entaille 15 mm
)
0,5
12000
K I (mpa , mm
10000
8000
6000
4000
2000
0 5 10 15 20 25
Taille de la fissure a (mm)
2/ Entaille en V :
Largeur
Facteurs d’intensité de contraintes (Mpa.mm0.5)
de l’entaille
↓ ↓ ↓ ↓ ↓ ↓ ↓
3mm 1702.28 2334 2946 4931 6695 7534
68
La variation de facteur d’intensité de contrainte (KI) en fonction de longueur de fissure (L)
pour l’entaille en V est représentée sur la figure suivante :
65000
60000 Entaille 3 mm
55000 Entaille 6 mm
50000 Entaille 9 mm
45000 Entaille 12 mm
Entaille 15 mm
KI (Mpa,mm )
40000
0,5
35000
30000
25000
20000
15000
10000
5000
0
0 5 10 15 20 25
Largeur
Facteur d’intensité de contraintes (Mpa.mm0.5)
de l’entaille
↓ ↓ ↓ ↓ ↓ ↓ ↓
3mm 2655 3152 5372 5326 5598 25011
69
25000
Entaille 3 mm
Entaille 6 mm
20000 Entaille 9 mm
Entaille 12 mm
)
Entaille 15 mm
0,5
15000
K I (mpa , mm
10000
5000
0 5 10 15 20 25
Taille de la fissure a (mm)
70
CONCLUSION
71
Conclusion
La simulation des éprouvettes de traction (selon la norme européenne EN 10002-
1 :1990) à entaille en V, circulaire et carré, d’une tôle en acier doux d’épaisseur 2mm, est faite
avec succès, nous donnons en dessous les principales conclusions des deux parties de l’étude :
1/ première partie :
Kt = -0.0046 (b/a) 4 - 0.0475 (b/a) 3 + 0.6546 (b/a) 2 - 2.1425 (b/a) + 4.68 (Entaille circulaire)
Kt = -0.0013 (b/a) 4 + 0.0108 (b/a) 3 - 0.0337 (b/a) 2 + 0.0142 (b/a) + 2.69 (Entaille en V)
Kt = 0.0108 (b/a) 4 - 0.14 (b/a) 3 + 0.6642 (b/a) 2 - 1.605 (b/a) + 4.42(Entaille carrée)
-On remarque que les contraintes au fond des entailles sont maximales, et ils diminuent qu’on
s’éloigne de celui-ci. [14]
- Le facteur de concentration de contrainte l’entaille circulaire varis selon trois phases avec
l’accroissement de l’ouverture de l’entaille. (voir figure)
1/ deuxième partie :
La simulation des éprouvettes entaillées avec fissure pour la détermination de facteur d’intensité
de contraintes KI pour différent forme d’éprouvette avec entaille en v, entaille circulaire et
entaille carrée.
- Les courbes de facteur d’intensité de contraintes pour l’entaille circulaire. (voir figures)
- Les courbes de facteur d’intensité de contraintes pour l’entaille en V. (voir figures)
- Les courbes de facteur d’intensité de contraintes pour l’entaille carrée. (voir figures)
72
-La figure de propagation de fissure suivante montre bien les iso-valeurs des contraintes de Von
Mises sur l’éprouvette de traction à entaille en V (profondeur de l’entaille a=3mm et largeur
b=3mm)
Donc le recours à l’approche de simulation nous à aider à gagner le temps et l‘argent, car
l’expérimental et le control destructif est plus couteux.
73
REFERENCES
74
Références
[2] DAVID B. (1982). Elementary engineering fracture mechanics. USA. Martinus nijhoff
publishers.
[3] DAVID B. (1982). The Practical Use of Fracture Mechanics. London. Kluwer Academic
Publishers.
[6] PILKEY W. D. PILKEY D. F. (2008). Peterson’s stress concentration factors. Third Edition.
John Wiley & Sons, Inc., Hoboken, New Jersey. p38
[10] Hamid R. K. Abdul K. M. (April 2012). Abaqus For Reinforced Concrete Structures.
Faculty of Civil Engineering – UTM. http://www.win2pdf.com.
[11] Dassault Systèmes Simulia Corp. (2014). Abaqus 6.14 example problems guide volume 1:
static and dynamic analyses. Providence, RI, USA.
[12] Benkrira S. (2017/2018). Calcul du facteur d'intensité de contrainte pour des sections
fissurées circulaire et elliptique avec Abaqus. Mémoire de master. Université Saad Dahlab de
Blida 1. Faculté de technologie. Département de génie mécanique.
[13] European standard. (1990). English version. Metallic materials - Tensile testing - Part 1:
Method of test at ambient temperature. EN 10002-1: 1990.
[14] Brick chaouche A. Tala- ighil N. (2012). Effet du type et de la taille d’un défaut sur le
comportement mécanique d’une plaque en alliage de titane soumise à un chargement axial. 3ème
Conférence Internationale sur le Soudage, le CND et l’Industrie des Matériaux et Alliages (IC-
WNDT-MI’12) Oran du 26 au 28 Novembre 2012. (pp. 74-8). http://www.csc.dz/ic-wndt-
mi12/index.php.
75
ANNEXES
76
Quelques géométrie et iso-valeurs des
contraintes obtenues par Abaqus
- Entaille en V (b=3mm) :
77
- Entaille en V (b=6mm) :
78
- Entaille en V (b=9mm) :
79
- Entaille en V (b=12mm) :
80
- Entaille V (b=15mm) :
81
- Entaille Circulaire (R = 3mm) :
82
- Entaille Circulaire (R = 6mm) :
83
- Entaille Circulaire (R = 9 mm) :
84
- Entaille Circulaire (R = 12 mm) :
85
- Entaille Carrée (b = 3 mm) :
86
- Entaille Carrée (b = 6 mm) :
87
- Entaille Carrée (b = 9 mm) :
88
- Entaille Carrée (b = 12 mm) :
89
- Entaille Carrée (b = 15 mm) :
90