Gm2 MFDP Cours
Gm2 MFDP Cours
de Monastir Mécanique
UE 8 : Fabrication Mécanique
2ème année Génie Mécanique
Séance Sujet
- Présentation du cours.
S01
- Introduction.
- Rappel : Elastoplasticité isotrope , Loi de Prandtl-Reuss : formulation,
S02
réponse en traction et en cisaillement ; identification.
S03 - Méthodes analytiques de calcul : les modèles ; méthode des tranches.
Figures de “Fundamentals
of modern manufacturing”
ISBN 978-0470-467002
Emboutissage Pliage Découpage
• Selon la Température
- Déformation à froid
→ Plasticité
Emboutissage, Forgeage, Laminage, Etirage
- Déformation à chaud
T ≥ 0.5 température de fusion
→ Viscoplasticité
Forgeage, Laminage, Filage
- Procédé
- Machine et outillages
• Choix: - Conditions d’essai (vitesse, température)
- Pièce initiale : matériau et dimensions
- Lubrification
UE 8 : Fabrication Mécanique
2ème année Génie Mécanique
Cadre général
Elastoplasticité 1D
Elastoplasticité 3D
Dépendance
Réversibilité Linéarité
/vitesse
Métaux :
- Elasticité linéaire
- Plasticité incompressible (à froid)
- Viscoplasticité (à chaud)
• Plasticité
ε
- Coulomb : σ = Ro s Ro : seuil de plasticité
Plasticité parfaite (sans écrouissage)
ε
ε
- Prager : σ -β ε = Ro Ro
Ro : seuil initial de plasticité
Ecrouissage cinématique linéaire ε β : Coef. Écrouissage ciném.
β
• Viscoélasticité
- Newton : σ = µ ε µ : coefficient de viscosité
Viscoélasticité linéaire
- Norton-Hoff :
Viscoélasticité non linéaire
(viscoplasticité)
• Loi élastique :
Linéaire isotrope, Loi de Hooke
E : Module d’Young
ν : Coefficient de Poisson
• Critère de plasticité :
Plasticité incompressible, Critère de Von Mises
: Fonction d’écrouissage
: Contrainte équivalente de Von Mises
σ 0 0 ε 0 0 ε = εe + εp ; εt = εet + εpt
σ : �0 0 0� ⇒ ε : �0 εt 0� Comportement élastique ⇒ ε t = -ν ε ;
e e
εe = σ/Ε
0 0 0 0 0 εt Incompressibilité plastique ⇒ εpt = - ½ εp ; ε� = εp
σ� = σ
Ecoulement plastique ⇒ pour σ < σο= σs(0) : εp = 0
pour σ > σο : σ = σs(εp)
→ Identification complète à partir d’un essai de traction
1000
Contrainte (Mpa)
800
600
⇒ E = 203 GPa
400
200
0
0 0,1 0,2 0,3 0,4 0,5
Déformation
: Fonction d’écrouissage
UE 8 : Fabrication Mécanique
2ème année Génie Mécanique
Abdelwaheb DOGUI
Février 2015
Modèles et Méthodes
• Comportement plastique:
Matériau rigide parfaitement plastique avec critère de Von Mises
⇒ ⇒
• Modèles de frottement:
- Modèle mixte : τ = τs
Si µ|p| > τs
- Modèle de Tresca : τ = m τs
- Tresca : τs = σs / 2
0≤m≤1
- V. Mises :
h x • Hypothèses :
- Matériau rigide parfaitement plastique
- Critère de Von Mises
- Les directions principales des contraintes sont constantes
2a / temps et espace
Largeur L >> a et h - Contraintes indépendantes de y : σx(x) ; σy(x)
Déformation plane - Le frottement de contact n’influence pas l’état de
contrainte au sein du matériau
p
- Le coefficient de frottement est constant
τ
σx σx + dσx
• Méthode : Equilibre de tranches fines
τ
p
p - Equilibre (x>0) ⇒ ; σy = p
- Déformations planes ⇒ σz = ½ (σx + σ y )
- Frottement de Tresca ⇒
τ = m τs
- Frottement de Coulomb ⇒
τ = µ |p|
mt = ; mc =
P : Force d’écrasement
- Equilibre ⇒
- Incompressibilité ⇒ σr = σθ
- Frottement de Tresca ⇒
τ = m τs
- Frottement de Coulomb ⇒
τ = µ |p|
P : Force d’écrasement
Fin S03
Mise en forme par déformation plastique 29
Ecole Nationale d’Ingénieurs Département de Génie
de Monastir Mécanique
UE 8 : Fabrication Mécanique
2ème année Génie Mécanique
Abdelwaheb DOGUI
Février 2015
Méthodes des bornes
La puissance dissipée par les efforts extérieurs est maximale pour la solution exacte
CSPA CCPA
Borne inférieure
Borne supérieure
P
x - CSA
⇒ σy2 =0 ; σx2 = σx1
2 1 2 (CL statiques)
Poinçon rigide
Longueur L
Pas de frottement - CPA σy1 = p = - 4τs
⇒
Déformation plane
(Critère) σx1 = - 2τs
v P
2. Détermination des vitesses
2a x Mouvement par glissement relatif entre les frontières
h (1) (3)
(2)
(4)
Poinçon rigide
Longueur L ⇑
Pas de frottement
Déformation plane
⇒
Puissance plastique
totale
4. Théorème Cinématique
5. Optimisation
Pinf = 8La τs ≤ P ≤
Fin S04
Mise en forme par déformation plastique 39
Ecole Nationale d’Ingénieurs Département de Génie
de Monastir Mécanique
UE 8 : Fabrication Mécanique
2ème année Génie Mécanique
Abdelwaheb DOGUI
Février 2015
Simulation numérique
Généralités
σ1= F/ S σ1 / σ2 (ε1)
σ2 / σ2 (ε1)
σ2= F/ So
Ιncompressibilité ⇒ S = S o e -ε1
Hauteur h = 20mm
- Plateaux Bas et Haut
Diamètre D = 20mm Rigides
• Eléments : CAX8
Axisymétrique, quadrilatère, quadratique
• Finesse de maillage
Maillage homogène, 3200 éléments
• Méthode de résolution
Implicite, Quasi statique
I. Hmila – A. Dogui
PFE-ENIM - 2006
1 17,5
8,8
R4 23°,28
118
Objectifs :
– Forme Dimensions et formes
– Répartition la plus homogène de l’épaisseur des outillages
– Maximiser l’épaisseur minimale
1000
contrainte (MPa)
800
600
400
200
0
0 0,1 0,2 0,3 0,4 0,5
déformation
Poinçon U1=UR3=0 et U2= -0,001m U1=UR3=0 et U2= -0.0228m U1=UR3=0 et U2= 0 propagé
20 -10
Déplacement serre flanc (µ) 0
(jeu) (écrasement)
• Objectifs :
– Répartition la plus homogène de l’épaisseur
– Maximiser l’épaisseur minimale
0,850 0,850
µ = 0,2
µ = 0,3
0,800 0,800 µ = 0,4
0,750 0,750
Epmin
Ep min
0,700
0,700
OPTIMISATION
0,650
0,650
0,600
0,600
0,130 0,150 0,170 0,190 0,210 0,230 0,250
0,020 0,040 0,060 0,080 0,100 0,120 0,140 Epmax - Epmin
sup\mi-mj\
0,850 Fs1<Fs2<Fs3
R2 0,850
Fs1
R3
0,800 Fs2
R4 0,800 Fs3
0,750
Epmin
Epmin 0,750
0,700 0,700
0,650 0,650
0,600 0,600
0,130 0,150 0,170 0,190 0,210 0,230 0,250 0,130 0,150 0,170 0,190 0,210 0,230 0,250
Epmax - Epmin Epmax - Epmin
2
1
Rmatrice=2 mm
Force d’emboutissage
µ = 0.4
Fs = 25 kN
µ = 0.3
µ = 0.2
µ = 0.2
Fs = 25 kN
Force d’emboutissage
R = 2 mm
R = 3 mm
R = 4 mm
Fs = 275 kN
Fs = 150 kN
Rmatrice = 4mm
Fs = 25 kN
µ = 0.3
0,85
0,84
0,83
Epmin 0,82
0,81
0,8
0,79
0,78
OPTIMISATION
0,77
0,76
0,1 0,12 0,14 0,16 0,18 0,2 0,22 0,24
EpMax - Epmin
0,86
0,86
0,85 µ = 0,21
§=0
µ = 0,25 0,85
0,84 § = 2,00E-05
µ = 0,3 0,84
§ = -1,00E-05
0,83 0,83
0,82 0,82
Epmin
Epmin
0,81 0,81
0,8 0,8
0,79 0,79
0,78 0,78
0,77 0,77
0,76 0,76
0,1 0,15 0,2 0,25 0,1 0,15 0,2 0,25
Epmax - Epmin Epmax - Epmin
Rmatrice = 4, µ = 0.21 et δ = 0
Rmatrice = 4, µ = 0.21 et δ = +2E-5m
Effort sur serre flanc
250000 R4 µ=0,3
R4 µ=0,25
R4 µ=0,21
R3 µ=0,3
R3 µ=0,25
200000
R3 µ=0,21
R2 µ=0,3
Effort sur serre flanc
R2 µ=0,25
R2 µ=0,21
150000
Fs(N)
100000
50000
(δ = 0)
0
0,00E+00 5,00E-03 1,00E-02 1,50E-02 2,00E-02
U2(m)
δ = 0, R=4 (R = 3, µ = 0.21)
Fin S05
UE 8 : Fabrication Mécanique
2ème année Génie Mécanique
Plasticité orthotrope
Abdelwaheb DOGUI
Janvier 2018
Isotropie / Orthotropie
Directions principales d’orthotropie
𝑚𝑚
��⃗3
𝒎𝒎𝟏𝟏 = 𝑚𝑚
��⃗1 ⨂𝑚𝑚
��⃗1
𝑚𝑚
��⃗2 𝒎𝒎𝟐𝟐 = 𝑚𝑚
��⃗2 ⨂𝑚𝑚
��⃗2 m1 + m2 + m3 = 1
𝑚𝑚
��⃗1
𝒎𝒎𝟑𝟑 = 𝑚𝑚
��⃗3 ⨂𝑚𝑚
��⃗3
• Isotrope : � 𝟐𝟐𝒗𝒗𝒗𝒗 = ½(σ22 – σ33)2 + ½(σ33 – σ11)2 + ½(σ11 – σ22)2 + 3σ232 +3σ312 + 3σ122
𝝈𝝈
Von Mises ∀ le repère
• Orthotrope : � 𝟐𝟐𝒉𝒉𝒉𝒉𝒉𝒉𝒉𝒉 = F(σ22 – σ33)2 + G(σ33 – σ11)2 + H(σ11 – σ22)2 + 2Lσ232 + 2Mσ312 + 2Nσ122
𝝈𝝈
Hill
���⃗𝒊𝒊
Dans le repère 𝒎𝒎
𝑚𝑚
��⃗3
σ σ 𝑚𝑚
��⃗2
11 11 𝑚𝑚
��⃗1
σ σ
22 22
σ σ
33 33
� 𝟐𝟐𝒉𝒉𝒉𝒉𝒉𝒉𝒉𝒉 = 2σ
𝝈𝝈 2σ = σ.H.σ = σ : H : σ
23 23
2σ 2σ
31 31
2σ 2σ
12 12
p 𝐻𝐻 p 𝐺𝐺 εp
��⃗1 ⇒
𝑒𝑒⃗1 = 𝑚𝑚 εt = − 𝐻𝐻+𝐺𝐺 εp εn = − 𝐻𝐻+𝐺𝐺 εp 𝜀𝜀̅ = 𝜎𝜎� = 𝜎𝜎√𝐻𝐻 + 𝐺𝐺
√𝐻𝐻+𝐺𝐺
p 𝐹𝐹 p 𝐻𝐻 εp
��⃗2 ⇒
𝑒𝑒⃗1 = 𝑚𝑚 εt = − εp εn = − 𝐻𝐻+𝐹𝐹 εp 𝜀𝜀̅ = 𝜎𝜎� = 𝜎𝜎√𝐻𝐻 + 𝐹𝐹
𝐻𝐻+𝐹𝐹 √𝐻𝐻+𝐹𝐹
p 𝐺𝐺 p 𝐹𝐹 εp
��⃗3 ⇒
𝑒𝑒⃗1 = 𝑚𝑚 εt = − 𝐹𝐹+𝐺𝐺 εp εn = − 𝐹𝐹+𝐺𝐺 εp 𝜀𝜀̅ = 𝜎𝜎� = 𝜎𝜎√𝐹𝐹 + 𝐺𝐺
√𝐹𝐹+𝐺𝐺
εtp
Coefficient de Lankford : r = r1 = H/G r2 = F/H r3 = G/F
εnp
𝑚𝑚
��⃗1 𝒑𝒑
ψ
σ 0 0 𝛆𝛆𝐩𝐩 𝛆𝛆𝟏𝟏𝟏𝟏 0 𝑒𝑒⃗𝑖𝑖
𝑒𝑒⃗1 𝒑𝒑 𝒑𝒑
σ : �0 0 0� ⇒ ε : �𝛆𝛆𝟏𝟏𝟏𝟏
p
𝛆𝛆𝒕𝒕 0�
0 0 0 0 0
𝒑𝒑
𝛆𝛆𝒏𝒏
c11 = (G + H )(−a1 sin 2 (2ψ ) + 2a2 sin 2 (ψ )) c12 = − H + (G + H )a1 sin 2 (2ψ )
c13 = (G + H )(−a1 sin( 4ψ ) + a2 sin( 2ψ )) / 2
(F + G + 4H − 2N )
a1 =
1 − 2a3 − 2a1 sin (2ψ )
2 4 (G + H )
r (ψ ) =
1 + 2a3 + 4a2 sin 2 (ψ ) F −G
a2 =
2 (G + H )
ro = H/G r90 = F/H r45 = (2N-F-G)/2(F+G)
G−H
a3 =
2 (G + H )
50 100
0° 0°
(daN/mm²)
40 80 15°
(daN/mm²)
15°
30° 30°
30 45° 60 45°
60° 60°
75°
σ
20
σ
40 75°
90°
90°
10 20
0 0
0 0,05 0,1 0,15 0,2 0,25 0,3 0,35
0 0,1 0,2 0,3 0,4 0,5 0,6
ε
3 ε
2,5
2 Fe P04 S2 DT
coefficient "r"
Fe P01
1,5 Inox18-10
S1
Inox19-11
1 alu A5 ψ
DL
0,5
0
0 20 40 60 80 100
Ψ°