Initiation au calcul des

structures
dans le domaine plastique

Elasto - plasticité
en petite transformation
Passer à la
Cours première page
 Objectifs

• Introduction à la plasticité « classique »

• Aborder la résolution de problème non linéaire

• Utiliser un code de calcul en non linéaire

Passer à la
première page
 Plan
• Aspect physique - modèles analogiques
• Elasto-plasticité des barres (1D)
étude des treillis
• Elasto-plasticité des poutres (3D ==> 1D)
étude des portiques (rotule plastique)
• Critères 3D & règles d ’écoulement plastique

Passer à la
première page
Supports :
• Polycopié de cours
• Sur le site intranet de l ’école
• Cours,
• Exercices, projets
• Doc d’utilisation et tutorials des logiciels
• ALGOR
•CASTEM
Evaluation :
• 2 notes :
• TD - TA sur 10
• Projet sur 10
Passer à la
première page
 Aspects physiques
Elastiques (instantanées - réversibles)

 E 

Déformations Visqueuses (fct du temps)


  

Plastiques (irréversibles - non linéaire)


 S 

Passer à la
première page
 Ecrouissage
 A imposé  A

B
B
mesuré
t 

les Essais Fluage - recouvrance

 imposé 
A
A

mesuré
t t

Relaxation


A
imposé A

mesurét
B
t
B
Passer à la
première page
 Modèles rhéologiques
Modèles « linéaires » ==> solides visco - élastiques


Ressort :     E
E

Amortisseur :      

Modèle de Maxwell :   E1   1 

Modèle de Kelvin-Voigt :   E  

Le ressort ou de l ’amortisseur peuvent être non linéaire

Passer à la
première page
Modèles non linéaires
Analogie mécanique Modèles de comportement Essai d’écrouissage

Rigide Plastique Parfait S
S
RPP 
plastique

Élasto-Plastique Parfait S
EPP 
p e

Rigide Plastique avec Écrouissage S
RPE 
E1 
Élasto-Plastique avec Écrouissage
S E2
E2

EPE E1 + E2 

Modèles de base de la plasticité

Passer à la
première page
 Essai de traction
o
F So F

observation Impossible à mesurer ==> Norme

 Ecrouissage monotone
A
o
Domaine élastique

Cycles charges- décharges  //


p(A) e (A)
Passer à la
première page
Ce n ’est pas si simple !!
2ème passage
 C

A
1er passage
o Il faut connaître
l ’historique
du chargement
O O' 
Caractère incrémental des
lois de comportement en
B plasticité
On utilise un temps cinématique =/ temps réel

Passer à la
première page
 Ecrouissage Ecrouissage ISOTROPE

Ecrouissage monotone

Même augmentation en
traction et compression

2o Wdef élastique


Ecrouissage CINEMATIQUE
Effet Baushinger
durcissement dans un sens
adoucissement dans l ’autre

Passer à la
première page
 Modélisation
 ET
Etat présumé
Etat réel
Critère
s(E)
D    ; f ( , E ) :    s ( E )  0 
Loi d ’écoulement plastique
Etat actuel (,E)
1
D p D e 
1
D D p  (  D   S ( E ))
E H
D  D p  D e ET
H  ET / (1  )
E E


Passer à la
première page
 Evolution élasto-plastique
des Treillis
Résolution Explicite
• Analytique (RDM)
• par la MEF
• Modèles EPP et EPE
• Chargements cycliques

Résolution Numérique (MEF)

• Algorithmes de résolution des Pb non linéaire
• Newton-Raphson et N-R modifié
• Algorithme de projection sur le critère
• Calcul du résidus {Fint}-{Fext}
Passer à la
première page
 Application
Structure à étudier Modèle EPP (pas d ’écrouissage)

h o

yo O 

xo
Effort donné

Passer à la
première page
 N 3  N3  N1
N2
N1 
 N 2  F  2 N1
3 5Décharge
Phase élasto-plastique
élastique en compression
N2 21Phase
Phaseélasto-plastique
élastique
En E ruine de la structure
F NN//22à12 oNS12  2 N1
A' oS 7 A 2 B En A
B ruine de la structure
C (F:élastique
6EnPhase = 0)
NF2 (1
o S 2 ): o S
Contraintes
2 résiduelles
1 2 1
7 Phase élasto-plastique
Domaine d ’élasticité  F1  N1   o S o S
1 21  2
On ne peut pas calculer
C' N  (1  2 ) S
directement les2 déformations
o

3 plastiques de la barre 2
oS Les déformations plastiques
N1
ne vérifient pas les conditions
6 de compatibilité
Le domaine d'élasticité
 p 2  0
n'est pas modifié aucours p 2  2 p1
C 
 p1    0
du chargement p3
car le
matériau est EPP
4 4 En D (F = - oS)
Il y a plastification de la barre 2 en
E compression
5 D D Passer à la
première page
N2
Ecrouissage ISOTROPE

N1

Passer à la
première page
N2
Ecrouissage CINEMATIQUE

N1

Passer à la
première page
 Résolution numérique
Principe
Incrément de charge donné DF 
Calcul élastique ==> DU   K 1DF
Pour chaque élément D e  Be Due
Élastique ou non ?
ALGORITHME DE PROJECTION ==> D 
 DN i 
==> DN i i e j DN j DFint e  
 DN j 

Assemblage et calcul du résidu {R} = {DFext}-{DFint}

Si R   Il faut itérer DU   K 1R

Passer à la
première page
Algorithme de résolution
•Maillage éléments finis

F •Définition de l ’historique du chargement (incréments DF)

•Définition des lois de comportement o ,E ,ET
•Calcul de [K]
Pour chaque incrément {DF}
F = K(u) U
Initialisation du résidu {R} <== {DF}
DF Tant que || {R} || > 
Calcul de [K(u)]
Solution cherchée Calcul de {DU} = [K]-1 {R}
Pour chaque élément
{R1} Convergence la plus rapide
déformation ==> projection {DFint}e
Convergence lente
fin pour
Nouveau résidu {R} = {R} - {DFint}
fin tant que
Impression des résultats de l ’incrément
fin pour

U Passer à la
première page
D  (u j  ui ) / 
D e  E D
PROJECTION sur le critère de plasticité

si +De > s(E) élément en cours de plastification

Etat présumé
   s(E)
Projection sur R  1
le critère d e
ET Cas particulier  = s(E) incrément plastique
De
Etat réel R 1
s(E)

si non +De < s(E) incrément élastique

Etat actuel (,E) D R0
 Dp De
E
ET
D p  R(1  )D
E
 D  D p  D e
D  RET  (1  R) E D

Passer à la
première page
A vous de jouer ...

Passer à la
première page
 Plasticité des poutres
Rappels bo
 y v,x 
Hypothèses de Bernoulli  
u ( M , t )  v 
  xx   y v,x2 1  xx  E xx 2
 0 
 

Sv  T,x  f L.C. Elastique

Equations « PFD » 3
T   M f ,x
Intégrée sur
la section

y
 +
y
 xx
 M f  EI v
x  xx ,x 2 4
Diagramme
Section des cont raint es

4 équations pour 4 champs ==> résolution

~  M f
En plasticité ce qui change c ’est la LC intégrée on pose ~  v,xx
Passer à la
première page
 ~
Objectif : obtenir la LC élasto-plastique   f (  )
~

Essai de flexion pure M M

==> Mf = cte

 y ~ ~  I
 
Début de plastification des fibres pour  xx       o   o o
I  y  h h

 o
Pour ~  ~o y c
E~
Zone
h
c
~ ~1   o Z ( h )
élast ique
Zones o
 xx
~1
plast iques - o
-c
-h ~0
Section Diagramme des cont raint es

~
LC élasto-plastique ~o
 I( c ) 
~   o   Z ( h )  Z ( c )
 c  Passer à la
première page
 ~1 h Z ( h )
Facteur de forme plastique d ’une section ~ 
o I

b b
S' S h
R R
2h 2h
e << R

~o 1  2 S
b h2 R3 b h2  e R2 h( S '  )
o 6 4 3 6
~1 1 2 4 3 S
bh R b h2 4 e R2 h( S '  )
o 3 3 4
~1 16 4 12  3( S S ')
2  1,7 1,5  1,27
~o 3  12  2( S S ')

Représente la réserve vis à vis de la plastification totale : 1 pas de réserve

Passer à la
première page
 M f  ~
F F
4
Essai de flexion simple
x

2

4~o
Fo 
Négligeons T et utilisons la LC ~  f ( ~ ) ==>

4~1
F1 


Phase élastique F  Fo Phase élasto-plastique Fin de plastificationF  F1

Fo  F  F1

~ ~ ~
~
~o
1
~o ~
o
x x

a 
x
a
2 2 
2

  
2 2 2

Passer à la
première page
 Etude de ~ ~
F  F1

Fo  F  F1

F  Fo
x

3
a 
2

Comportement asymptotique ==> notion de rotule plastique

Modèle simplifié - rotule plastique ~

~1
Modèle simplifié
Application ~0 LC élasto-plastique
Etude des portiques
charges limites ~
~o
Passer à la
première page
 Effort donné F

/2 /2

yo A C B


xo
Phase élastique
RA M A RB  RA  RB  F  0
F   11
M  F   R   0 
 A A
R
 A 16 F
A

B 2 ==>  3
W 5 M A  F
 0  RB  F  16
RB 16
~
critère 16 ~1
M f  ~1 ==> Rotule plastique en A pour F1 
MA

3 
/2 x


MC

Phase élasto-plastique
~1  ~1 F critère
RA F RB RA    2 ~1
==>  M C   1  F2  6
~
~1 F 
A
B
 RB   
   2 F2 charge de ruine
Passer à la
première page
A vous de jouer ...

Passer à la
première page
 Plasticité 3D
Critère (fonction de charge)
f (  , E )  0 État actuel est intérieur au domaine d ’élasticité
f (  , E )  0 État actuel se situe sur la frontière du domaine (état plastique)
Un état extérieur au domaine d’élasticité est physiquement impossible à obtenir
Le domaine d’élasticité représente l’ensemble des états  admissibles.
Ecrouissage « E »
2 Modèle EPE 2 Modèle EPP
B Domaine init ial
f (E) 
A A Domaine act uel
Domaine actuel
B

O 1 O 1
fo Domaine init ial fo  f ( E )

f (  , E ) Doit respecter les symétries matérielles Passer à la

première page
Critère de Von Mises (1910)
Pour les métaux
P 
P    P   De
 ==> utilisation du déviateur des 
 P
1
 d    tr(  ) 1
3
Fonction de charge Partie sphérique

fo (  )   d :  d  k 2
2
==> essai de traction : k o
3
 3 
 ; d :d -s  0 
 2 

3
1,1,1 VM Contrainte de Von Mises
1
 VM   1   2 )2   2   3 )2   3   1 )2 
fo() 2  
2
Passer à la
1 première page
Critère de Tresca (1870)
 ; i, j  1,3 i  j -s  0 
Contraintes principales

T Contrainte de Tresca  T  max  i   j

ij

Critère dit de cisaillement maximal

si 1   2   3   max  ( 1   3 ) 2 2 max -  s  0

1,1,1
3
3
Von Mises

fo() T resca

2 1 2

Tresca plus sévère de 13%

1
Passer à la
première page
Loi d ’écoulement plastique
Principe de Hill   *  Domaine d'élasticté Dp   * : d  p

D p  d : d  p Dissipation plastique
d  p associée à d 
==> De est convexe
dp est normal à la surface de charge

d p f (  , E )
d p     grad  ( f (  , E ))

Loi de normalité
n

d p d p  
i 1
i grad  ( fi )

Point singulier,
f ( , E ) la direction de la déformation
plastique appartient à un cône

Von Mises
fo (  )   d :  d  k 2 ==> grad  ( f )  2  d ==> d  p    d
Passer à la
première page
Principe de la résolution numérique
Soit dU  la solution de K( u ) dU   dF 
Pour chaque élément d e   DBdU  Aux points de Gauss

La sol. Supposée élastique      d e 

Si f (  , E )  0 OK

Si non il faut projeter sur De

Solution supposée  Solution supposée
élastique
n élastique
Etat actuel, avant d e
incrém ent de charge B
( , E )
( , E )
d
Solution cherchée Solution numérique

f ( , E ) f ( , E )
Schéma idéal schéma numérique
Passer à la
première page