Plasticité 1
Plasticité 1
Plasticite
Théorie et applications
20/08/2024
Table des matières
Préambule .................................................................................................................................. iv
CHAPITRE 1 ............................................................................................................................. 1
GENERALITÉS ......................................................................................................................... 1
1.1. Les Grandes classes des matériaux .................................................................................. 1
1.1.1. Classification selon les propriétés mécaniques ......................................................... 1
1.1.2. Classification selon leur microstructure .................................................................... 1
1.2. Echelle et moyens d’observation ..................................................................................... 2
1.3. Méthodes expérimentales ................................................................................................ 2
1.3.1. Selon le chargement .................................................................................................. 3
1.3.2. Selon que la contrainte ou la déformation soit imposée ........................................... 4
1.4. Les grandes classes de comportement macroscopique .................................................. 11
1.4.1. Sous charges monotones ......................................................................................... 11
1.4.2. Charge-décharge ..................................................................................................... 12
1.4.3. Essai cyclique.......................................................................................................... 14
1.5. Diversité des mécanismes de déformation microscopiques .......................................... 16
CHAPITRE 2 ........................................................................................................................... 17
PLASTICITE UNIAXIALE .................................................................................................... 17
2.1. Généralités ..................................................................................................................... 17
2.2. Modélisation du comportement élasto-plastique ........................................................... 17
2.2.1. Modèles parfaits ...................................................................................................... 18
2.2.2. Modèles avec écrouissage ....................................................................................... 18
2.3. Fonctions de charge ....................................................................................................... 25
2.3.1. Modèle élastique-parfaitement plastique ................................................................ 25
2.3.2. Modèle élastique-plastique avec écrouissage isotrope............................................ 26
2.2.2. Modèle élastique-plastique avec écrouissage cinématique ..................................... 27
2.4. Loi d’écoulement plastique ........................................................................................... 28
2.4.1 Comportement parfaitement plastique .................................................................... 28
2.4.2. Ecrouissage cinématique linéaire ............................................................................ 29
2.4.3. Ecrouissage isotrope ............................................................................................... 29
2.5. Viscoplasticité ............................................................................................................... 30
2.5.1. Influence de la vitesse de chargement..................................................................... 30
2.5.2. Influence de la température .................................................................................... 31
2.5.3. Modèles Rhéologiques ............................................................................................ 31
CHAPITRE 3 ........................................................................................................................... 38
i
PLASTICITE MULTIAXIALE ............................................................................................... 38
3.1. Critères de plasticité ...................................................................................................... 38
3.1.1. Critère de Rankine .................................................................................................. 39
3.1.2. Critère de Tresca ..................................................................................................... 39
3.1.3. Critère de von Mises ............................................................................................... 42
3.1.4. Critère de Mohr-Coulomb....................................................................................... 46
3.1.5. Critère de Drucker-Prager ....................................................................................... 48
3.2. Ecrouissage .................................................................................................................... 50
3.2.1. Ecrouissage isotrope ............................................................................................... 53
3.2.2. Ecrouissage Cinématique ........................................................................................ 54
3.3. Loi d’écoulement plastique ........................................................................................... 55
3.3.1. Règle de partition .................................................................................................... 55
3.3.2. Principe du travail plastique maximal..................................................................... 56
3.3.3. Potentiel d'état et relations d'état ............................................................................. 56
3.3.4. Règle de normalité des matériaux dits standards ................................................... 57
3.3.5. Règle de normalité des matériaux dits non standards ............................................. 68
CHAPITRE 4 ........................................................................................................................... 75
Analyse limite à la flexion ....................................................................................................... 75
4.1. Introduction ................................................................................................................... 75
4.2. Moment plastique .......................................................................................................... 75
4.3. Loi moment-courbure .................................................................................................... 77
4.4. Rotule plastique ............................................................................................................. 78
4.5. Redistribution plastique des efforts et mécanismes de ruine ......................................... 79
4.6. Méthodes de calcul plastique ......................................................................................... 80
4.6.1. Méthode de l’équilibre ultime ................................................................................. 83
4.6.2. Méthode pas à pas ................................................................................................... 87
4.6.3. Méthodes règlementaires ........................................................................................ 89
4.6.4. Capacité de rotation ................................................................................................ 89
4.7. Influence de l’effort tranchant sur le moment plastique ................................................ 91
4.7.1. Section rectangulaire ............................................................................................... 92
4.7.2. Section en I ............................................................................................................. 95
4.7.3. Formules de l’EC3 .................................................................................................. 96
4.8. Influence de l’effort normal sur le moment plastique ................................................... 98
4.8.1. Section rectangulaire ............................................................................................... 98
4.8.2. Section en I ............................................................................................................. 99
4.8.3. Formules de l’EC3 ................................................................................................ 101
CHAPITRE 5 ......................................................................................................................... 102
ii
INTÉGRATION NUMÉRIQUE DES LOIS DE COMPORTEMENT ................................. 102
5.1. Introduction ................................................................................................................. 102
5.2. Schémas d’intégration numérique ............................................................................... 102
5.2.1. Méthode d’Euler avant .......................................................................................... 103
5.2.2. Méthode d’Euler arrière ........................................................................................ 103
5.2.3. Méthode de Runge-Kutta d’ordre 2 ...................................................................... 104
5.2.4. La Thêta méthode ................................................................................................. 106
5.3. Intégration de lois élasto-plastiques ............................................................................ 106
5.3.1. Algorithme Euler explicite ................................................................................... 107
5.3.2. Algorithme Euler implicite ................................................................................... 108
5.3.3. Algorithmes de type Return Mapping................................................................... 109
Exemples d’applications ........................................................................................................ 123
Bibliographie .......................................................................................................................... 199
iii
Préambule
Je serai reconnaissant aux lectrices et aux lecteurs qui voudront bien me faire part par mail de
leurs critiques, commentaires, suggestions ou corrections.
[Link]@[Link]
soltani_mohamed_redha@[Link].
iv
CHAPITRE 1
GENERALITÉS
Monocristal Polycristal
1
SiO2 cristallin SiO2 amorphe
SiO2 : Silicate de sodium
Structure mixte
Fig. 1.1. Microstructures des matériaux
2
Il existe de nombreux essais qui permettent de caractériser les propriétés mécaniques des
matériaux. Certains sont normalisés :
AFNOR - Association Française de NORmalisation ;
ISO - International Standardisation Organisation ;
Allemagne - DIN 50 125, 50 145 ;
ASTM - American Society for Testing and Materials.
3
[Link]. Charge-décharge
On augmente la force F puis on la ramène à 0. La Figure 1.4 représente la réponse à une série
de charges-décharges consécutives dans le cas d’un comportement avec écrouissage.
Essai d’écrouissage : essai pour lequel la déformation est imposée à vitesse constante.
4
Fig. 1.6. Ecrouissage
Fluage + recouvrance : essai pour lequel la contrainte est imposée grâce à une fonction
échelon suivi d’attente à contrainte nulle. On mesure la déformation longitudinale supposée
homogène dans le temps. La courbe de réponse caractérise la déformation différée (viscosité).
La déformation croît en fonction de t de façon continue et monotone jusqu’au point de
décharge A (fluage). Après décharge, la réponse comporte une réponse instantanée suivie
d’une évolution monotone décroissante de en fonction de t (recouvrance). Cette
recouvrance est dite totale si la déformation s’annule pour t suffisamment grand.
Retour instantanée
Fluage
Recouvrance
Relaxation
Retour
instantané
Effacement
ment
5
1.3.2. Selon la direction des contraintes
[Link]. Les essais uniaxiaux classiques
L’essai de traction simple : pour identifier le comportement macroscopique d'un matériau, on
réalise des essais uniaxiaux sur éprouvette, i.e. des essais où l'on exerce des forces de traction
ou de compression dans une direction (soit e1 ). Le tenseur des contraintes sera alors de la
forme :
11 0 0
0 0 0 e1 e1
0 0 0
Machine de traction
6
Fig. 1.9. Essai de traction
L’essai de torsion : réalisé sur éprouvette pleine, cet essai est essentiellement utilisé à haute
température pour connaître l'aptitude à la mise en forme des métaux. L'avantage de ce type
d’essai est d'éviter la striction. Par contre, il est d'interprétation difficile, dans la mesure où
l'état de contrainte et déformation n'est pas uniforme. Il est possible de remédier à ce dernier
inconvénient, en adoptant comme éprouvettes des tubes minces qui peuvent être instrumentés
localement, à l'aide de jauges ou d'extensomètres.
L’essai de torsion effectué sur un tube mince (Figure 1.10) permet d’imposer dans la zone
utile de l’éprouvette un état de contrainte de la forme :
0 0 0
0 0 e e z e z e
0 0
M
2 R 2 e
7
de l'essai de traction.
L’essai de flexion : L’essai de flexion est un test mécanique statique. Cet essai permet de
caractériser des propriétés intrinsèques et extrinsèques du matériau étudié.
Cet essai permet par exemple de connaitre :
- la résistance à la traction du béton ;
- la déformabilité des matériaux (courbes charge-flèche) ;
- le comportement mécanique du matériau.
8
F M
2 Re 2 R 2 e
9
Les fentes dans les bras sont utilisées pour rendre uniforme les distributions de contrainte et
de déformation dans la zone centrale.
10
[Link]. Essais complexes
Outre les essais de traction-torsion sur tube, il existe d’autres moyens de générer des états de
contraintes multiaxiales contrôlés dans des éprouvettes. C’est le cas d’essais de traction-
pression interne sur tube (Figure 1.17).
Fig. 1.18. Loi de comportement comportant une partie linéaire et une partie non linéaire
- Le durcissement du matériau. Le durcissement correspond au cas d’une courbe (,)
croissante (toute augmentation de déformation requiert une augmentation de contrainte).
11
- L’adoucissement du matériau. L’adoucissement correspond à une courbe(,) décroissante
(la contrainte diminue alors que la déformation augmente). Pour détecter correctement
l’adoucissement, il est nécessaire de contrôler l’essai en déformation (le pilotage en contrainte
ne peut se poursuivre au-delà de la contrainte maximale et provoque une instabilité de l’essai).
1.4.2. Charge-décharge
- L’anélasticité : si les courbes de charge-décharge coïncident, le milieu est élastique
(éventuellement non linéaire). Dans le cas contraire il est anélastique. Après décharge
complète, il subsiste une déformation résiduelle. Cette déformation résiduelle peut s’effacer
avec le temps, signe d’une viscosité du matériau. Dans les matériaux insensibles à la vitesse
de chargement (métaux laminés à froid par exemple), cette déformation résiduelle est
permanente, tant que l’application d’une contrainte ne vient pas la perturber.
12
Ecrouissage : L’anélasticité ne se manifeste que lorsque la force atteint un certain seuil, qui
est la limite d’élasticité du matériau. Initialement ce seuil est 0 (Figure 1.21b). Dans la
plupart des métaux, si on effectue une décharge de A à B puis on recharge en B, la charge se
fait à nouveau le long du trajet AB et de façon élastique (le trajet AB est réversible) et selon le
coefficient de proportionnalité E. La limite d’élasticité pour une charge à partir de B est donc
A. Le matériau est dit écrouissable.
13
Fig. 1.23. Endommagement couplé à la plasticité
14
Sous chargement cyclique, les éléments de construction se rompent sous des contraintes
inférieures à la résistance ultime mesurées dans les essais statiques, parfois même inférieures
à la limite d’élasticité (s’il s’agit d’un matériau ductile). On observe également que plus la
variation des contraintes est forte, plus le nombre de cycles de sollicitations nécessaires à
produire la rupture diminue. Ce phénomène de la décroissance de la résistance du matériau
aux efforts variables dans le temps s’appelle fatigue.
La courbe de Wöhler est la seule qui permet de visualiser la tenue de la pièce ou des
matériaux dans le domaine de fatigue. Elle définit une relation entre la contrainte appliquée σ
(parfois notée S) et le nombre de cycles à la rupture NR (en fait nombre de cycles pour lequel
on observe P% de ruptures). La courbe de Wöhler est appelée courbe S-N (Stress - Number of
cycles) dans les pays anglo-saxons. En pratique, la courbe de Wöhler est généralement donnée
pour une probabilité de rupture p = 50%.
Sur la courbe de Wöhler, on peut distinguer trois domaines (Figure 1.26) :
Domaine de fatigue oligocyclique où les ruptures surviennent après un petit nombre de
cycles (N 103….105). Il correspond à des contraintes élevées pour lesquelles se produit une
déformation plastique macroscopique du matériau.
Domaine de l’endurance limitée où les ruptures surviennent après un nombre de cycles
croissant avec la décroissance de la contrainte (103….105 < N 107….108). Dans ce cas, la
rupture survient sans être accompagnée d’une déformation plastique d’ensemble mesurable.
La réponse de l’éprouvette peut être purement élastique (comportement élastique dès les
premiers cycles) ou bien adapté (comportement devenu élastique après stabilisation cyclique).
Domaine d’endurance illimitée ou zone de sécurité où les ruptures ne se produisent pas
avant un nombre de cycles supérieur à la durée de vie envisagée de la pièce (N > 107….108).
En deçà d’une valeur seuil D, dite limite de fatigue, il n’y a jamais de rupture par fatigue quel
que soit le nombre de cycles appliqué. Cette limite peut ne pas exister ou être mal définie pour
certains matériaux (aciers à haute résistance, métaux non ferreux).
103 107
15
1.5. Diversité des mécanismes de déformation microscopiques
Le comportement macroscopique est en fait le résultat de déformations locales à une échelle
microscopique, voire à une échelle inférieure. Cet aspect microscopique est fondamental pour
la compréhension physique des phénomènes et relève du domaine des matériaux.
A l’échelle microscopique, les différents mécanismes de déformation observés sont :
clivage : rupture des liaisons interatomiques dans une direction perpendiculaire au plan de
rupture. Ce mécanisme est dominant dans les matériaux fragiles (Figure 1.27).
a) Glissement en bloc
PLASTICITE UNIAXIALE
2.1. Généralités
Les chargements uniaxiaux correspondent à un tenseur de contrainte comportant une seule
composante non nulle. La courbe conventionnelle de l'essai de traction (Figure 2.1),
indépendante de tout effet de viscoplasticité ou de vieillissement s’appuie sur les concepts
fondamentaux suivants :
- les déformations restent infinitésimales tout au long du trajet de chargement ;
- la partition des déformations, qui décompose les déformations totales en déformations
élastiques e (récupérable par décharge) et déformations plastiques p (déformations
permanentes après décharge) : = e + p ;
- le seuil de plasticité ( = y) en sollicitation uniaxiale défini la frontière entre les domaines
élastique et plastique ;
- la règle d’écoulement plastique, qui définit la façon dont évoluent les déformations
plastiques.
18
domaine élastique (Figure 2.4b). Le modèle cinématique respecte l’effet Bauschinger
couramment observé pour les matériaux métalliques, à savoir un durcissement dans un sens
(sens de l’écoulement plastique) et un adoucissement d’égale amplitude dans le sens contraire
(Figure 2.4b). Dans la réalité, les matériaux ont un écrouissage qui est la combinaison des
deux phénomènes.
a) b)
Fig. 2.4. Illustration des deux principaux types d’écrouissage
Considérons un incrément de contrainte dprit à partir d’un état actuel situé sur la courbe
d’écrouissage (Figure 2.5). Soit ET la pente de la courbe contrainte-déformation supposée
constante dans l’intervalle d. ET est appelé le module tangent. L’accroissement de la
contrainte est défini par :
d ETd
avec :
d
d d p
E
d’où :
d d
d p
ET E
1 1
d d p
ET E
E
d 1 T E T d p
E
On définit le coefficient d’écrouissage h comme étant :
d
h
d p
soit :
ET
h
E
1 T
E
19
h
d
dp ET
h d
de dp
d
20
Fig. 2.7. Association ressort-patin en parallèle (Modèle de Prager)
Si on rajoute un ressort assemblé en série au patin, Le modèle analogique devient élasto-
plastique (Figure 2.8).
21
s’exprime en fonction de la déformation nominale à partir l’expression suivante :
L
dL L L L0 L0
v ln L ln L 0 ln ln ln 1 n
L0
L L0 L0
Les formules précédentes sont valables si l’allongement est uniforme, ce qui n’est pas le cas
lorsque la striction apparait après le dépassement de la contrainte maximale.
Dans le domaine des petites déformations, la courbe conventionnelle et la courbe rationnelle
sont confondus (ln(1 + n) n). C’est le cas des courbes dans le domaine élastique. Les
courbes rationnelles sont moins arrondies que les courbes conventionnelles et se placent au-
dessus de celles-ci après plastification du matériau.
La Figure 2.9 compare les courbes conventionnelle et rationnelle pour un matériau donné. On
constate que les deux courbes divergent d’autant plus que les déformations plastiques sont
importantes.
Courbe rationnelle
Courbe rationnelle corrigée
pour tenir compte de la
striction
Courbe conventionnelle
Striction
F F F F
S S durant la striction
n
v k pv
ln v ln k n ln pv
22
C’est l’équation d’une droite de la forme y = ax+b avec a = n (pente de la droite) et b = log K.
a) Aluminium b) Acier
Fig. 2.11. Comparaison entre les points de la courbe de traction rationnelle et la modélisation
de Hollomon [13]
Cette loi très simple convient surtout pour les aciers doux ou mi-durs et les métaux ou alliages
à faible limite d’élasticité.
En dehors de la loi puissance, différentes relations ont été proposées pour la caractérisation du
comportement plastique des métaux. Les principales lois proposées sont :
n
Loi de Ludwik v Y k pv ou ln v ln Y ln k n ln pv
Cette loi est bien adaptée aux aciers doux et les aciers faiblement alliés
23
n
Loi de Swift ou loi de Krupkowski v k 0 pv ou ln v ln k n ln 0 pv
Y k 0
n
Loi de Voce v Y S Y 1 exp pv
où s est la contrainte de saturation.
Les paramètres S et sont ajustées par la méthode des moindres carrés à partir de courbes
expérimentales.
Dans la littérature, on trouve une autre forme de la loi de Voce :
v Y 1 A exp B pv
où A et B sont les paramètres rhéologiques spécifiques à la loi de Voce qu’on détermine à
partir de la pente et l’ordonnée à l’origine de la droite d’équation :
ln v ln Y ABpv ax b
où a = ln(V) et b = ABpv
La loi de Voce comporte une contrainte seuil contrairement à celle de Ludwik où la contrainte
continue toujours à monter.
c) Loi de voce
Fig. 2.12. Lois d’écrouissage non linéaire
24
1
Loi Ramberg Osgood v Y k Y p M
v
Y ou ln v Y ln k
1
My
ln pv
25
point donné, à partir de l’état initial naturel, il est possible de déterminer l’état dans lequel se
trouve le matériau selon les conditions suivantes :
d
- Si f < 0 : comportement purement élastique, d d e d p 0
E
d
- Si f = 0 et df = d < 0 : décharge élastique d d e d p 0
E
- Si f = 0 et df = d = 0 : écoulement plastique : d d p 0 d e 0
- Si f > 0 est un état de contraintes impossible.
La condition f = 0 désigne le critère de plasticité.
Ee E p
2.3.2. Modèle élastique-plastique avec écrouissage isotrope
L’écrouissage isotrope (Taylor et Quinney, 1931) correspond à l’expansion que peut subir le
domaine d’élasticité. En d’autres termes, le domaine d’élasticité voit sa taille augmenter, mais
qui reste centré sur l’origine (Figure 2.15). La variable d’écrouissage qui intervient dans f est
la dimension du domaine d’élasticité, notée R. C’est un paramètre scalaire mesurant la taille
du domaine. L’évolution de cette variable est la même quel que soit le signe de la vitesse de
déformation plastique. La fonction de charge peut se mettre alors sous la forme :
f , y , R y R p
Dans bien des cas, les utilisateurs ne prennent pas la peine de définir une forme explicite de la
loi de comportement, et décrivent la courbe de traction point par point. Cela revient
implicitement à considérer un écrouissage isotrope.
L’écriture la plus simple pouvant être considérée pour R(p) est une fonction linéaire de la
déformation plastique, de sorte que :
R p hp p
26
h est le module d’écrouissage et p étant la déformation plastique cumulée (voir § 2.4.3).
R(p) peut être aussi défini par une fonction puissance ou une fonction exponentielle comme
par exemple la loi de Hollomon, loi de Ramberg–Osgood (voir § [Link]).
f , X X P y
y
2y
- y
a) Ecrouissage non linéaire b) Ecrouissage linéaire de Prager
Fig. 2.16. Ecrouissage cinématique
27
Dans le cas d’un écrouissage mixte, la fonction de charge s’écrit :
f , X, R X y R
En régime élastique, la vitesse de déformation plastique est bien entendu nulle, en d’autres
termes on peut écrire :
f 0 p 0
Pendant l'écoulement ( = y), la vitesse de déformation élastique devenant à son tour nulle et
la surface de charge et l’incrément de la déformation plastique sont orthogonaux. La vitesse
de déformation plastique ̇ , on peut écrire :
f ( ) 0 et 0
28
En résumé, pour un modèle élasto-plastique parfait on peut écrire :
f ( ) 0, 0 et f ( ) 0
Sous chargement uniaxial, le tenseur de vitesse de déformation plastique est une diagonale, la
variation de la déformation plastique cumulée est égale à la valeur absolue de la vitesse de la
déformation plastique (voir § [Link]) :
p p
29
2.5. Viscoplasticité
Les modèles précédemment décrits correspondent à des réponses dépendants de l'histoire du
chargement, mais instantanée : aucun temps caractéristique n'intervient, contrairement à ce
qui est observé lors des essais de fluage, de relaxation et des essais de traction à vitesse de
déformations constante ou avec des sauts de vitesse ou sous hautes températures. Dans ces
cas, les effets visqueux doivent être pris en compte. A noter que la viscosité n'influe pas sur
les déformations réversibles (élastiques), mais n'entre en jeu qu'à l'apparition des
déformations permanentes (plastiques).
On parle de viscoélasticité lorsqu’on effectue une décharge complète, aucune déformation
résiduelle ne subsiste après un certain temps (phénomène de retard). Le matériau ne retrouve
pas instantanément sa configuration initiale comme dans le cas d’un comportement élastique
parfait (Figure 2.18a). Le comportement d'un matériau viscoélastique linéaire est ainsi
intermédiaire entre celui d'un solide élastique et celui d'un liquide visqueux.
On parle de viscoplasticité lorsqu’après décharge, la déformation résiduelle diminue de façon
continue en fonction du temps mais dont une fraction reste toujours présente (Figure 2.18b).
a) b)
Fig. 2.18. Comportement viscoélastique et viscoplastique
30
Fig. 2.19. Courbe de traction typique avec sauts de vitesse
Une façon de mettre en évidence la dépendance de la courbe de comportement à la vitesse de
déformation est de faire des essais de traction uniaxiaux en pilotant en déformation et en
imposant une vitesse de déformation constante. On constate que, plus la vitesse de
déformation est grande, plus la contrainte est élevée pour le même niveau de déformation
(Figure 2.20). La sensibilité à la vitesse de chargement traduit le caractère visqueux du
matériau.
31
Fig. 2.21. Influence de la température
32
C’est une équation différentielle du premier ordre, dont la solution peut s’écrire sous la forme
suivante :
E
t
t Ae
a) b)
Fig. 2.23. Modèles de Maxwell et Voigt
Pour un essai de fluage, l’équation différentielle du solide de Voigt s’intègre avec la condition
(0) = 0 pour donner la solution suivante :
0
t
(t) 1 e
E
où la constante , homogène à un temps, est le « temps de fluage » ou « de retard ».
E
Sous l’effet d’une contrainte σ0 constante en fonction du temps, la déformation tends vers la
valeur asymptotique σ0/E.
33
La particularité du modèle de Voigt est de ne pas présenter d’élasticité instantanée. L’application
d’un saut de déformation en t = 0 produit une contrainte infinie.
La Figure 2.24 montre la réponse du modèle de Kelvin-Voigt (ou de PT) à un essai de fluage
et de relaxation.
Modèle de Zener
Le modèle de Voigt n’est pas utilisable en relaxation, sauf si la mise en charge est
progressive. Pour cette raison, on rajoute au modèle de Voigt un ressort pour effectuer des
calculs de structure (modèle de Zener ou modèle du solide linéaire standard (SLS)).
=
a) Fluage à contrainte constante b) Relaxation à déformation constante
Fig. 2.24. Réponse du modèle de Voigt
a) b)
Fig. 2.26. Modèles rhéologiques de Zener
34
[Link]. Modèles Viscoplastiques
Modèle de Norton-Hoff
Norton a proposé en 1929 une loi pour représenter le fluage stationnaire et unidimensionnel et
décrire le comportement viscoplastique parfait de nombreux matériaux et de certains aciers à
haute température. Dans ce modèle, la viscosité linéaire ou non linéaire est schématisée par un
simple amortisseur sans considération des déformations élastiques et du seuil initial (Figure
2.27).
Elle a été généralisée par Odqvist en1934 puis par Hoff en 1954 pour le cas des sollicitations
multidimensionnelles. Elle est appelée depuis : loi de Norton-Hoff (Power law).
La loi de Norton-Hoff sans écrouissage s’écrit sous la forme suivante :
1
vp N
où :
vp est la vitesse de déformation viscoplastique ;
N est le paramètre du matériau.
e vp signe pour Y
E
e
pour Y
E
35
avec :
Y si Y 0
f
0 si Y 0
Si n = 1, on obtient le modèle plastique de Bingham à seuil (Bingham plastic) qui est une loi
de Newton à seuil.
En relevant le ressort en série (E→∞), on obtient le modèle rigide viscoplastique.
Modèle élasto-viscoplastique avec écrouissage de Bingham généralisé
Le modèle de Bingham généralisé s'obtient par combinaison d’un patin, d'un amortisseur et
d’un patin assemblé en parallèle auxquels on ajoute un ressort en série (Figure 2.29). La
détermination des équations de ce modèle s’effectue en considérant les équations de
comportement individuelles de chacun des éléments :
X H vp v vp p Y
- vp = 0 p H vp Y
36
La Figure 2.30 représente compare entre les réponses des modèles considérés ci-dessus.
vp
Fig. 2.30. Comparaison entre divers modèles de comportement viscoplastiques
Loi de Johnson-Cook
Johnson et Cook ont proposé une loi de comportement destinée aux grandes déformations et
aux hautes températures. Elle est aussi généralement employée lorsqu’une large gamme des
vitesses de déformation est utilisée. Elle se compose, dans l'ordre d'apparition dans l’équation
donnée ci-dessous, d'un module traduisant l’écrouissage du matériau, d'un module traduisant
la sensibilité à la vitesse de déformation et d’un module traduisant la sensibilité à la
température :
A B n 1 C ln * 1 T *m
avec :
*
r
T*
T Tr
Tm Tr
* est la vitesse de déformation homologue correspondant au rapport entre la vitesse de
déformation imposée et une vitesse de déformation de référence r .
T est la température homologue, Tr est une température de référence, Tm est la température de
fusion, n est le coefficient de la déformation à la contrainte et m le coefficient de la température à la
contrainte. A , B , n , C et m sont les paramètres du matériau à identifier.
Afin de traduire les couplages existants entre déformation, température et vitesse de
déformation, d’autres variantes de la loi de Johnson-Cook ont été proposées.
37
CHAPITRE 3
PLASTICITE MULTIAXIALE
38
est le tenseur des contraintes (ordre 2), de composantes ij (i, j =1, 2 ou 3), où i indique la
direction du vecteur force et j indique la direction de la normale à la face sur laquelle la force
est appliquée ;
, et sont respectivement les tenseurs (d’ordre 2) des vitesses de déformation totale, de
e P
D :
p
3.1.1. Critère de Rankine
Dans le cas des matériaux fragiles sujets au clivage (qu’ils soient polycristallins ou non), le
critère de la contrainte normale maximale exprime que la rupture se produit lorsque la plus
grande contrainte sur la facette principale (contrainte principale majeure) dépasse une valeur
critique 0 indépendante de l’orientation de la facette. Le critère de Rankine s’écrit donc :
f Sup (n) 0 0
Fig. 3.2. Représentation du critère de la contrainte normale maximale dans le plan de Mohr.
L’état de contrainte est représenté par les cercles de Mohr. a) Etat de contrainte physiquement
admissible, le critère n’est pas atteint. b) Etat de contrainte non physiquement admissible. c)
Etat limite : critère atteint.
39
produit lorsque la contrainte de cisaillement maximale sur une facette n atteint une valeur
critique 0. Le critère de Tresca s’écrit donc :
f () Sup (n) 0 0
Si Sup (n ) 0 : il n’y a pas de glissement
Si n tel que (n) 0 il y a un glissement sur le plan de la normale n où le cisaillement
limite est atteint.
Les états de contrainte autorisés par le critère sont représentés sur la Figure 3.3 par des cercles
situés sous la droite (horizontale dans le plan de Mohr) =0. Tout état de contrainte
conduisant à des cercles de Mohr coupant cette droite ne respecte pas le critère (Figure 3.3b).
Fig. 3.3. Représentation du critère du cisaillement maximal dans le plan de Mohr. a) Etat de
contrainte physiquement admissible. b) Etat de contrainte non physiquement admissible. c)
Etat limite : critère atteint.
Le cisaillement critique 0 peut être déterminé par un essai de traction simple où les
contraintes principales intermédiaire 2 et mineure 3 sont nulles. Le cisaillement maximal
étant égal à la moitié du diamètre du grand cercle de Mohr vaut 0/2 dans un essai de traction
simple où la limite élastique vaut 0 (ou y). Le cisaillement critique y est donc lié à la limite
d’élasticité y en traction simple par : y = y/2.
40
ou encore :
f ( ) max i j y 0
i, j
f 1 y si 0 2 1
f 1 2 y si 2 0 1
Les 3 autres expressions sont obtenues par symétrie par rapport à l’axe d’équation 1 = 2.
Fig. 3.5. Représentation du critère de Trésca dans le plan des contraintes principales σ1 et σ2
Dans l’espace des contraintes principales, le critère de Tresca est représenté par un cylindre à
base hexagonale (Figure 3.6).
41
3.1.3. Critère de von Mises
Le critère dit de von Mises fut formulé par Maxwell en 1865. Huber (1904) le développa
partiellement dans un article en polonais. Cependant, sa paternité est généralement attribuée à
von Mises (1913). On parle aussi parfois de la théorie de Maxwell–Huber–Hencky–von
Mises, ou de critère de Prandtl-Reuss.
Fig. 3.6. Représentation du critère de Tresca dans l’espace des contraintes principales σ1, σ2 et
σ3
Le critère de von Mises, dit critère du plus grand travail de distorsion, suppose que la rupture
est due à la part de l’énergie de déformation provoquée par les seules déformations angulaires
(ij). Il peut ainsi s’énoncer comme suit : l’énergie fournie pour augmenter ou diminuer le
volume initial ne joue aucun rôle dans la rupture de l’équilibre élastique, seule l’énergie
fournie pour déformer le volume entre en ligne de compte. L'écoulement plastique est dit
alors incompressible. L'hypothèse d'incompressibilité plastique est admise si à tout instant :
Tr
p p
11 22
p
33
p
0
L’énergie due au changement de forme sans changement de volume se calcule par la formule :
Uf
1
12G
2 2 2
11 22 22 33 11 33 6 122 223 312
1
12G
1 2 2 3 1 3 2
2 2 2
J
2G
où ij (i 3, j 3) est le tenseur des contraintes dans une base quelconque. i (1 i 3) est le
tenseur des contraintes principales, J2 est le second invariant du déviateur du tenseur des
contraintes et G est le module de cisaillement.
J2
1
6
2 2
1
1 2 1 3 2 3 Tr s
2
2
2
1
6
2 2
11 22 22 33 33 11 6 12
2 2
223 31
2
1
s:s
2
où s est le tenseur des contraintes déviatoriques.
42
Tr 1 ij kk ij
1
s sij
3 3
11 p 12 13 1 p 0 0
21 22 p 22 0 2 p 0
31 32 33 p 0 0 3 p
avec :
i 1 j1
oct
1
3
s:s
1
3
Tr s
2 1
3
1 2 1 3 2 3
2 2 2
oct lim
La contrainte octaédrique limite octlim peut être déterminée par la formule suivante obtenue en
traction uniaxiale :
1 1
oct 2 I2 oct lim 2 2y
3 3
Les matériaux de l’ingénieur sont le plus souvent considérés comme isotropes, c’est-à-dire
qu’aucune direction de l’espace n’est privilégiée, la fonction de charge f est donc invariable
quel que soit le repère choisi pour exprimer le tenseur des contraintes.
Le critère de von Mises peut s’écrire sous les formes suivantes :
- Dans la base principale :
1 2 2 3 1 3
2 2 2
22y 0
43
Le critère de von Mises peut aussi s’exprimer par :
f 6 12 y 0
1
11 22 22 33 11 33
2 2 2 2
23
2
31
2
f eq y
3 3
s : s y s y 3J 2 y 0
2 2
eq est la contrainte équivalente de von Mises qui traduit un état de contrainte complexe en
une contrainte de traction qui lui est − vis-à-vis du critère de von Mises – équivalente :
1
2
1 2 2 3
2 2
1 3
2 3
2
Tr s
2
3 24 1 1
f 3J 2 y y y
2 9 9 9
Dans le plan traction-cisaillement, lorsque seules les composantes σ = σ11 et τ = σ12 sont non
nulles; les expressions des critères se réduisent à :
f
- von Mises (forme Huber) : f , 2 32 y et m y
3
- Tresca : f , 2 42 y et t y
2
La comparaison des critères de Tresca et de von Mises dans le plan traction-cisaillement est
représentée sur la Figure 3.7.
von Mises
Tresca
Fig. 3.7. Comparaison des critères de Tresca et de von Mises dans le plan traction-
cisaillement
44
Dans les deux cas, les surfaces prédites par les deux critères sont des ellipses passant chacune
par les deux points de coordonnées ( = y, = 0) et ( = - y, = 0). On remarque que le
critère de Tresca est conservatif par rapport à Von Mises : en effet, les limites admissibles
avec Tresca sont inférieures ou égales à celles admissibles avec Von Mises. L’écart maximal
entre les deux critères est de l’ordre de 13%.
Dans le plan des contraintes principales σ1 et σ2, le critère de von Mises s’écrit :
f 12 22 12
0.5
y
La surface seuil prédite par le critère de Tresca est une surface polygonale présentant des
points anguleux, tandis que la surface prédite par le critère de von Mises est régulière (Figure
3.8).
Fig. 3.8. Comparaison des critères de Tresca et de von Mises dans la plan des contraintes
principales σ1 et σ2
Dans l’espace des contraintes principales, la surface limite est un cylindre à base circulaire,
circonscrit au cylindre à base hexagonale de Tresca, axé sur la trisectrice du repère (1, 2,
2
3). L’intersection du cylindre de Von Mises avec le plan Π est un cercle de rayon 0
3
(Figure 3.9).
45
Le critère de von Mises peut s’écrire sous la forme suivante :
2 2 2 2
s:s y s : s y s y
3 3 3
L’isotropie des deux critères de Tresca et von Mises impose que la permutation de deux
arguments ne change pas la valeur de la fonction :
f (1 , 2 , 3 ) f (2 , 1 , 3 )
Les critères de Tresca et de von Mises sont tous deux insensibles à la pression hydrostatique
(1 = 2 = 3). Les deux critères sont toujours vérifiées quelques soient la valeur des
contraintes 1, 2 et 3.
Les critères de Tresca et de von Mises vérifient l’hypothèse d’équivalence de la traction et de
la compression. Cette hypothèse se traduit, pour tout état de contrainte par :
f () f
Il est à noter que la plupart des métaux ne vérifient pas cette hypothèse simplificatrice (effet
Bauschinger).
46
Fig. 3.10. Représentation du critère de Mohr Coulomb dans un diagramme -
Les valeurs courantes de l’angle de frottement interne sont comprises entre 15 et 45°. Les
valeurs inférieures ou autour de 30° sont typiques des argiles, tandis que des valeurs
supérieures caractérisent les sables (entre 25 et 45°). Les sols cohérents ont une cohésion
comprise entre quelques kilopascals et plusieurs centaines de kilopascals.
Le critère de Mohr-Coulomb peut également s’écrire sous la forme :
1 k p 3 c
où :
1 sin
kp tan 2
1 sin 4 2
c est la résistance à la compression simple :
2c cos
c
1 sin
47
Dans le plan (p,q), le critère de Mohr-Coulomb s’écrit :
sin sin
f () J 2 cos p sin c cos 6 6
3
C’est l’équation d’une droite (Figure 3.12).
a) b)
Fig. 3.13 Représentation de la surface de charge de Mohr-Coulomb a) dans l'espace
tridimensionnelle des contraintes principales, b) dans le plan π des déviateurs des
contraintes
48
deuxième invariant du déviateur et de la trace du tenseur des contraintes I1. C'est toujours un
cercle dans le plan déviateur, ce qui permet d’éviter les singularités associées aux angles, mais
qui dépend de 1’altitude sur la trisectrice des axes σ1, σ2, σ3 de contraintes principales (Figure
3.14). L’expression de la fonction de charge s’écrit :
f ( i ) J 2 I1 k
et k sont deux paramètres qui sont déterminés à partir des résultats d’essais. si = 0, la loi
se réduit à celle de von Mises.
Dans le repère des contraintes principales, elle s’écrit sous la forme :
1
6
2 2 2
1 2 2 3 1 3 K 1 2 3
Si t est la contrainte limite en traction simple du matériau, le critère de Drucker-Prager
implique :
t
K t
3
2 c t 1 t c
k
3 c t 3 c t
2sin 6c cos
k
3 3 sin 3 3 sin
49
Fig. 3.15. Comparaison entre les critères de Drucker-Prager et le critère de Mohr-Coulomb
dans l’espace des contraintes principales
3.2. Ecrouissage
Pour les matériaux écrouissables (aciers, argiles, etc.), lorsqu’on continue la charge au-delà de
la surface seuil initiale, le matériau se déforme plastiquement. Puis, si l’on décharge
l’éprouvette, on revient en régime élastique. Si après recharge la nouvelle surface seuil de
plasticité est identique à la surface seuil initiale, le matériau est élastique parfaitement
plastique. Si la nouvelle surface seuil est différente de la surface initiale, celle-ci traduit
l’écrouissage du matériau. La forme de la nouvelle surface seuil dépend du trajet de
chargement parcouru au-delà du domaine d’élasticité initial et les variables d’écrouissage du
matériau. La Figure 3.16, tirée des travaux de Bui, présente l’évolution de la surface seuil
d’un acier doux, pour deux trajets de chargement différents. Le premier trajet est purement
radial (Figure 3.16a), comme ceux employés pour la détermination de la surface seuil initiale,
mais le chargement a été poursuivi au-delà de l’offset. Le second trajet (Figure 3.16b) n’est
pas radial et consiste en une traction simple (suivant l’axe horizontal) suivi d’une torsion à
force de traction constante (portion verticale). La surface seuil actuelle est différente dans les
deux cas confirmant ainsi le rôle essentiel du trajet de chargement.
a)
50
b)
Fig. 3.16. Evolution de la surface seuil pour deux histoires de chargement.
A la différence du comportement élasto-plastique parfait, pour un comportement avec
écrouissage, il existe un domaine d'élasticité initial (généralement convexe) du matériau dans
l'espace des contraintes de Cauchy et un domaine d'élasticité actuel (subsequent yield surface)
, dépendant de l’histoire de déformation plastique caractérisé par les variables cachées
nommées variables d’écrouissage, comme il est illustré sur la Figure 3.17.
51
a)
b)
c)
Fig. 3.18. Illustration de l’écrouissage isotrope, cinématique et mixte
L’écrouissage des matériaux peut se manifester par une distorsion (changement de forme de
la surface seuil) et une translation, comme il est illustré sur la Figure 3.19.
52
3.2.1. Ecrouissage isotrope
Dans le cas d’un écrouissage isotrope, la fonction de charge s’écrit sous la forme :
f , R f y R(p) 0
R hp
53
3.2.2. Ecrouissage Cinématique
Dans le cas d’un écrouissage cinématique, la fonction de charge s’écrit sous la forme :
f , X f X y 0
X est le tenseur des contraintes internes qui existent à l’échelle microscopique dues aux
défauts et aux hétérogénéités. Le tenseur X décrit le déplacement du centre du domaine
élastique dans l’espace des contraintes. Il est symétrique et de trace nulle.
Le modèle le plus simple consiste à relier le tenseur des contraintes internes X aux tenseurs
des vitesses de déformations plastique par la relation :
X H: P
Le modèle linéaire connue sous le nom de loi de Prager est le plus simple. Il est caractérisé
par un module d’écrouissage dépendant du seul scalaire H :
2 p
X H
3
H est la pente de l’écrouissage cinématique linéaire.
Dans le cas de la traction uniaxiale, le tenseur des contraintes internes X a pour
composantes :
X11 0 0
2 X11
X 0 0
3 2
X
0 0 11
2
Le tenseur (déviateur) des vitesses des contraintes internes a pour composantes :
P
11 0 0
2 P11
X H 0 0
3 2
P11
0 0
2
En raison de l'incompressibilité des déformations plastiques, on doit avoir :
tr(X) tr X tr( ) 0
p
La loi d’écrouissage cinématique non linéaire proposée par Armstrong et Frédéricks conduit à
54
l’équation différentielle suivante :
X C Xp
p
Dans l’hypothèse des petites déformations, la déformation totale d’un matériau élasto-
e P
S:
e
55
où S désigne le tenseur des complaisances élastiques.
Le tenseur de déformation élastique s’écrit :
1
tr 1
e
E E
E et désignant respectivement le module deYoung et le coefficient de Poisson du matériau ;
1 est le tenseur unité.
L’incrément de déformation plastique dp dépend a priori du trajet de chargement et de l'état
d'écrouissage actuel du matériau, de l'état de contrainte et de l'incrément de chargement
appliqué et qu’on peut calculer à l’aide de l’expression :
S :
p e
:
* p
0
La dissipation plastique dans l’état de contrainte réel est donc maximale parmi tous les états
de contrainte physiquement admissibles. Le principe du travail plastique maximal a des
conséquences importantes puisqu'il impose que le domaine de réversibilité soit convexe et que
la vitesse d'écoulement plastique, lorsque les contraintes sont sur le seuil de plasticité, est
orientée suivant la normale extérieure au domaine La vitesse de déformation plastique étant
nulle lorsque les contraintes sont à l’intérieur du domaine d’élasticité.
56
– VK (K ≥ 2) représentent l´écrouissage, si besoin l’endommagement ; pour
l’endommagement l’un des VK sera noté D, variable de dommage.
Une manière de décrire complètement l’état thermodynamique du milieu est de définir un
potentiel thermodynamique, à partir duquel se déduisent les lois d’état.
L’énergie libre est une fonction des variables d’état observables et internes ; elle doit être
concave par rapport à la température, et convexe par rapport aux autres variables d’état pour
vérifier a priori les conditions de stabilité thermodynamiques.
e ,T, p , VK
57
Fig. 3.21. Règle de normalité pour une surface seuil régulière
avec :
f
n ij est le gradient de f par rapport à ij.
ij
Le tenseur vitesse de déformation s’écrit donc :
f f f
11
p
12p
13p
11 12 13
f f f
p21 p22 p23
p
31
p
32p
3"p 21 22 23
f f f
31 32 33
Le comportement élasto-plastique peut être synthétisé par la condition de consistance de
Karush-Kuhn-Tucker :
f (, A k ) 0 limite du comportement élastique
0 écoulement plastique
f (, A k ) 0 synthèse de l'alternative élastique plastique
Les variables internes sont définies par la règle de normalité généralisée dans le cas d‘un
écrouissage mixte :
f , X, R
p
f , X, R
X X
f , X, R
r
R R
58
direction de la vitesse de déformation plastique est dans le cône des normales extérieures aux
sommets (Figure 3.22).
La règle de normalité est alors généralisée au sens suivant : la vitesse de déformation
plastique en un point singulier de la surface seuil de plasticité appartient au cône des normales
extérieures à P en ce point. Elle s’exprime par l’expression suivante :
fi
, Ak
N
i
P
i 1
avec :
i 0 si fi , Ak 0
i 0 si fi , Ak 0
Fig. 3.22. La règle de normalité pour une surface de charge avec un point singulier
f eq y
3
sij : sij y
2
Le gradient de f par rapport à est calculé à partir de la règle de dérivation en
chaîne (ou règle de la chaîne, chain rule selon l'appellation anglaise) :
f f s mn
n
ij s mn ij
où :
mn
s mn mn kk
3
s mn mn mn kk
ij ij 3 ij
avec :
mn
im jn 1seulement si i m et j n
ij
kk 11 22 33 1 si i j
ij ij 0 si i j
59
On en déduit :
s mn 1
im jn mn ij
ij 3
f f 1
im jn mn ij
ij s mn 3
f
1 f 1 f
im jn mn ij im jn mn ij
s mn
3 s mn 3 s mn
f 1 f
ij
sij 3 s mm
or smm =J1 et la fonction de charge f ne dépend que de J2, d’où :
1 1
f f 3 2 1 3 2
n ij sij : sij sij : sij 3sij
ij sij sij 2 22
3 sij sij
2 eq 2
sij
3
La loi d´écoulement plastique déduite de la règle de normalité associée s’écrit alors sous
forme tensorielle :
3 s s s
p
2 eq 2 2
s y
3 3
ou sous forme indicielle :
3 sij sij sij
ijp
2 eq 2 2
sij y
3 3
Pour un matériau régi par le critère de von Mises, la vitesse de déformation plastique est
colinéaire au déviateur des contraintes.
Le critère de von Mises est applicable aux matériaux incompressibles, la déformation
s’effectue donc sans variation de volume, d’où :
tr
p p
11 22
p
33
p
0
La vitesse de déformation plastique cumulée est colinéaire au déviateur des contraintes, d’où :
2 p 2 p
p : :
p
3 3
La déformation plastique cumulée mesure la longueur du trajet représentant l’écoulement
dans l’espace des déformations plastiques, toutes directions confondues, soit, au temps t :
t
p(t) p(s)ds
0
60
Sous chargement uniaxial, on obtient :
1 0 0 1 0 0
2 3 s
s 0 0.5 0 eq n
0 0.5 0 signe()
3 2 eq
0 0 0.5 0 0 0.5
11
p
0 0 1 0 0
0 0.5 0 0 0.5 0 signe
p p
11
0 0 p
0.511 0 0 0.5
(e1 e1 0.5e2 e2 0.5e3 e3 )signe
Fig. 3.23. Règle de normalité dans la plan des contraintes principales σ1 et σ2 pour le critère
de von Mises
On en déduit la vitesse de déformation cumulée :
2 p
2 2
11
p
11
p
2
p 11 11
p
3 2 2
f
1p 1
f
p2
p
p 2
3 f
3
61
La règle de normalité s’exprime à l’aide d’un ou deux multiplicateurs selon les cas :
Point régulier du segment BC : 1 2 3
f 3 1 y
1
0
p
1
Point anguleux E :
f1 1 3 y f 2 1 2 y
p 1 0
3 1
avec :
1p 2p 3p 0
3
3 - 2 = y
3 - 1 = y
1<2< 3
2<1<3
1 - 2 = y 2<3<1
2 - 1 = y
1<3< 2
3<2< 1 3<1< 2
2
1
1 - 3= y 2 - 3 = y
62
module d’Young E et du coefficient de Poisson . Comme il a été mentionné au § 3.1.4, le
critère de rupture peut être défini par l’une des équations suivantes.
- Dans le plan de Mohr (,) :
tan c
- Dans le plan des contraintes principales :
1 k p 3 c
Lorsque le point représentatif de l'état de contraintes est situé sur la surface de charge, les
composantes du vecteur vitesse de déformation plastique ont pour valeurs :
1p
p f
n p2 0
p
3 k p
1 sin
kp 1 3p 1p
1 sin
Fig. 3.25. Critère de Mohr-Coulomb et loi de normalité associée dans le plan (1, 3)
En conséquence, il se produit un raccourcissement dans la direction verticale et une dilatation
dans la direction latérale. Pour = 30°, kp = 3, la dilatation latérale vaut dans ce cas trois fois
le raccourcissement dans la direction verticale, ce qui n’est pas conforme avec les résultats
des essais triaxiaux.
63
Le taux de déformation volumique plastique vaut donc :
pp 1p 2p 3p 1 k p 0 puisque k p 1
D:
p
D : n
En multipliant les deux membres de l’équation par la normale à la surface d’écoulement, on
obtient la relation suivante :
n : n : D : n n n : D : n : D : n 0
On en déduit :
n:D:
n:D:n
On voit bien que le multiplicateur plastique est indéterminé étant donné qu’il n’est pas
possible d’obtenir directement le tenseur vitesse des déformations en fonction du tenseur
vitesse de contrainte. Au même niveau de contrainte peut correspondre plusieurs déformations
plastiques.
Pour un matériau régi par le critère de von Mises, le multiplicateur plastique vaut :
2
n:
3
64
Comportement avec écrouissage mixte
Dans le cas d’un écrouissage mixte, la condition de cohérence s’écrit :
f f R
f
X :
X : X
p
p
n : n : H : hp
p
n : n : H : n h 0
On en déduit :
n: n:D:
n:H:nh n:H:nh
La courbe décrite par R(p) est donc la courbe d’écrouissage en chargement uniaxial
monotone. R(p) peut être défini point par point ou par une fonction telle qu’il a été décrit au §
[Link]. 0 = y.
On en déduit :
n: n:
p
dR h
dp
n : n
2 p 2 2
p :
p
3 3
65
puisque :
3 s 3 9 s:s 9 s:s 3
n eq
2
s:s n :n
2 eq 2 4 eq 4
2
3 2
s:s
2
Le module plastique h est défini comme étant la pente de la courbe eq-p (eq= R+y) :
deqdR
h
dp dp
La loi de Prandtl-Reuss permet de déterminer la direction et l’intensité de l’écoulement
plastique :
f n :
p
n
h
On admet que la courbe (p) obtenue sous chargement uniaxial est similaire à la courbe
eq(p) représentant le comportement sous chargement multiaxial (Figure 3.27) :
d deq
h
d p dp
eq
p
a) b)
Fig. 3.27. Courbes contrainte-déformation plastique et R-déformation plastique cumulée
Loi de Levy-Mises
Dans les cas des matériaux isotrope, il est raisonnable d’admettre que les composantes de la
vitesse de déformation plastique sont raisonnablement proportionnelles aux composantes du
tenseur des contraintes déviatoriques. La loi d’écoulement de Levy-Mises s’écrit :
pxx yy zz pxy pyz pyz
p p
0
s xx s yy s zz s xy s yz s yz
Soit :
ijp sij
f ij eq 0
3
sij : sij y R 0
2
66
Par ailleurs, en utilisant en plus les expressions suivantes :
ijp sij
2
sij : sij eq
3
On trouve :
2 p p 2 2 2 2
dp ij : ij sij : sij eq eq
3 3 3 3 3
Donc :
3 dp 3 deq 3 deq
et ijp sij
2 eq 2 heq 2 heq
avec :
deq
h
dp
Loi de Prager
C’est la loi obtenue en utilisant le critère de von Mises et une règle d’écrouissage cinématique
linéaire. La loi d’écrouissage linéaire cinématique proposée par Prager introduit un tenseur
déviatorique d’écrouissage X associée à la déformation plastique :
2 2
X H: H : n
p
3 3
La fonction de charge s’écrit :
f , X eq X y
avec :
f s X 3 s X
eq X s X : s X et n
3
2 2 eq
67
La condition de cohérence d’écrit :
f f
: : X n : n : X n : n : Hn 0
X
3
n : n , on en déduit :
2
n :
H
Le multiplicateur plastique a la même expression formelle que dans le cas de l’écrouissage
isotrope. Il faut néanmoins noter que la définition de n est modifiée, et que H est constant.
0 si f 0
p
g
p
0 si f 0
p
Fig. 3.29. Directions des déformations plastiques reportées au potentiel plastique dans le plan
(p,q)
La fonction seuil f et le potentiel plastique g sont deux fonctions scalaires différentes (Figure
3.30) et sont déterminées sur une base expérimentale.
Des indications sur le potentiel plastique peuvent être fournies par les résultats expérimentaux
si l’on reporte la direction des incréments de déformations plastiques dpp et dqq construite
au préalable sur la surface de charge dans le plan (p, q) ou à tout le moins sur la frontière du
domaine d’élasticité initial. Généralement, les directions des déformations plastiques ne sont
pas strictement normales à la surface d’état limite (Figure 3.31). Le plus délicat reste ensuite
68
la détermination de la fonction analytique capable de décrire les directions observées
expérimentalement.
q f
i 1 i 1
Fig. 3.31. Directions des déformations plastiques reportées sur la surface de charge de Mohr-
Coulomb obtenues expérimentalement [Graham, Noonan and Law, 1983]
70
La représentation graphique sur la Figure 3.34 de la surface de charge de Mohr-Coulomb dans
l'espace des contraintes principales et dans le plan montre six plans qui s'intersectent deux à
deux suivant une arête anguleuse, et se rejoignant au sommet d’un cône. Ces arêtes, au
nombre de six, ainsi que le sommet du cône, forment des singularités (les fonctions f et g ne
sont pas dérivables en ces points). La surface seuil de plasticité appartient donc au cône des
normales extérieures aux points angleux (voir § [Link]), comme il est illustré sur la Figure
3.34. Dans ce cas, la fonction seuil est définie par deux fonctions f1 et f2. A titre d’exemple,
pour le point anguleux pour lequel 1 > 2 = 3 = 0, les fonctions de charge s’écrivent :
Aux points anguleux, la règle de normalité s’écrit pour les deux mécanismes :
g1 g 2
1 2
p
Dans les directions principales, les gradients g1 et g2 par rapport à ont pour expressions :
g1 g 2
1 sin 1 sin
1 1
g1 g1
g 2 g 2
0 1 sin
2 2
g g
1 1 sin 3 0
3
3
71
On en déduit le vecteur vitesse de déformation plastique :
1p 1 sin
1 sin 1 sin 1 2
p2 1 0 2 1 sin 2 1 sin
p
p 1 sin 0
3 1 1 sin
Le taux de déformation volumique plastique vaut :
pp 2sin 1 2sin 2 2sin 1 2 0
Les taux des déformations totales sont directement mesurables lors des essais.
Une fois les taux des déformations volumiques et des déformations axiales mesurés, l’angle
de dilatance peut être estimé à partir de la formule obtenue par Vermeer en 1984 :
pp
p
arc sin p 1
p
p 2
1
La Figure 3.35 représente la modélisation de l’essai triaxial drainé sur un sable par la loi
élastique - linéaire parfaitement plastique avec critère de rupture de Mohr-Coulomb et la
variation volumique plastique de l’échantillon de sol sur la variation de la déformation
verticale.
A partir du signe de l’angle de dilatance, on peut déduire :
- si > 0 la déformation s’effectue avec une augmentation de volume (dilatance)
- si = 0 la déformation s’effectue à volume constant
- si < 0 la déformation s’effectue avec une diminution de volume (contractance).
L’angle de dilatance est généralement compris entre 0 et 15°. Les sables lâches et les
argiles ont des valeurs très faibles, valant quelques degrés à peine, voire zéro. Si = 0, la
rupture se produit à volume constant. D’une manière générale, l’angle de frottement est
pratiquement toujours supérieur à l’angle de dilatance. Une relation empirique simple et
72
généralement bien vérifiée reliant l’angle de dilatance à l’angle de frottement interne a été
proposée par Vermeer pour les sables : = - 30°. Par ailleurs, dans le cas de matériaux
fortement contractants à la rupture, l’estimation de l’angle de dilatance peut conduire à des
valeurs négatives.
Fig. 3.35. Modélisation de l’essai triaxial non drainé avec critère de rupture de Mohr-Coulomb
L’angle de distance peut aussi être déterminé à partir d’un essai de cisaillement direct par la
relation :
p22
tan
12p
avec :
p22 : déformation axiale plastique dans la direction verticale
p12 : distorsion.
En compression biaxiale en conditions drainées après avoir subi une consolidation isotrope,
les directions axiale et transversale sont confondues avec les directions principales majeure et
mineure respectivement. On peut alors écrire les égalités suivantes :
g
1 p
1 1 sin
1
2 p
g 1 sin
2
2
D’où l’on déduit l’expression obtenue par Roscoe en 1970 :
1p 2p
sin
1p p2
73
Pour l’ensemble des essais sur des sols isotropes décrits ci-dessus, on définit le taux de
dilatance comme étant le rapport des taux de déformation plastique volumique et
déviatorique :
pp
p
q
74
CHAPITRE 4
4.1. Introduction
En étudiant le comportement des structures en béton armé, en acier ou mixte, il a été constaté
qu’en phase élasto-plastique jusqu’à la rupture, la distribution des efforts différent de celle
obtenue à l’aide du calcul élastique. Ce processus de redistribution est dû essentiellement à la
plastification ou à la fissuration des matériaux.
Les méthodes d’analyse limite ont pour objectifs la détermination des chargements ultimes
des structures constituées de matériaux ayant un comportement élastique parfaitement
plastique ( y et illimité) et chargée de façon monotone dans l’hypothèse des petites
perturbations. L’hypothèse des petites perturbations implique que les changements de
géométrie de la structure étudiée sont négligeables au cours de son évolution. On se réfère
alors à la configuration initiale de la structure lorsqu’il s’agit d’exprimer les équations
d’équilibre et éventuellement de compatibilité des déformations.
On peut admettre que les changements de géométrie sont négligeables dans les cas suivants :
- structures rigides ;
- portiques tenues latéralement (contreventés).
dA A2 A1 y 0 A2 A1
A
En d’autres termes, les aires des deux côtés de l’axe neutre plastique doivent être égales (A1=
A2).
Le moment plastique de la section transversale s’obtient alors par la formule :
A
M pl Wpl y y1 y 2 A1 y y1 y 2 y
2
A
où Wpl y1 y 2 est le module de résistance plastique de la section transversale.
2
76
16 4 12 3S / S'
2 1.7 1.5 1.27
3 12 2S / S'
Tableau 4.1. Facteurs de forme plastique pour différentes sections
Plus ce rapport tend vers l’unité, meilleure est la section du point de vue élastique (toutes les
fibres plastifient en même temps). La section transversale en I du Tableau 4.1 est optimale.
77
Fig. 4.4. Loi moment-courbure pour diverses formes de section droite
78
Les méthodes de calcul plastique admettent qu’il se forme dans une section de moment
maximal une rotule plastique qui diffère d’une articulation mécanique puisque elle équilibre
un moment non nul égal à Mpl (Figure 4.6).
79
rapport à ceux obtenus par les méthodes de calcul élastiques (méthode des déplacements,
méthode des forces etc.), comme il est illustré sur la Figure 4.8..
Fig. 4.8. Redistribution des moments fléchissants dans une poutre continue
Si la structure est p fois hyperstatique et s’il y’a création de p rotules plastiques
supplémentaires, elle devient isostatique. La naissance d’une rotule plastique supplémentaire
a pour effet de transformer la structure en un système hypostatique, c'est-à-dire elle se
transforme en un mécanisme qui s’écroule par suite d’une augmentation quelconque de la
charge (Figure 4.9a). Une structure isostatique atteint l’état limite ultime dès l’apparition de la
première rotule plastique. La ruine de la structure peut se produire aussi de deux autres
façons :
- mécanisme partiel qui n’affecte qu’une partie de la structure (Figure 4.9b) ;
- mécanisme surabondant, qui pour des raisons de symétrie par exemple, fait intervenir un
nombre de rotules plus grand que le nombre nécessaire (Figure 4.9c).
80
- le calcul à l’état limite ultime (méthode de l’équilibre ultime) en faisant abstraction de son
comportement antérieur. Pour le calcul à l’état limite ultime, il existe deux approches :
l’approche cinématique et l’approche statique.
Considérons une poutre bi-encastrée soumise à une charge uniformément répartie d’intensité
q que l’on fait croitre de façon monotone. La poutre est deux fois hyperstatique. Trois rotules
plastiques doivent se former dans les sections critiques pour produire un mécanisme. Pour cet
exemple, les rotules se forment aux encastrements et au milieu de la poutre.
Tant que les contraintes produites dans les fibres extrêmes des sections critiques sont
inférieures à la limite élastique y, les moments fléchissants critiques, aux encastrements M1
et au milieu de la poutre M2, sont obtenus à partir de la théorie de la résistance des matériaux
par les formules :
ql2 ql2
M1 M2
12 24
Donc :
M1
2
M2
q
Mp
EI
l
-Mp
q1=12Mp/l2
0.5Mp
q=4Mp/l2
0.5Mp -Mp
q2=16Mp/l 2 q
16Mp/l2
Mp
12Mp/l2
Mécanisme
Déformée + w
Courbe charge- flèche 2 2
5Mp l2/96EI pour q =q 3Mp l /96EI 8Mp l /96EI
Moments résiduels
81
Au-delà de cette valeur de la charge, d’autres fibres au droit des encastrements vont se
plastifier sans avoir une influence importante sur la distribution des moments fléchissants
jusqu’à l’apparition de deux rotules plastiques au droit des encastrements. Les moments
d’encastrement atteignent alors le moment de plastification complète Mpl. Les deux rotules
plastiques apparaissent pour une charge égale à la valeur suivante :
12M p
q1
l2
Le moment au milieu de la poutre vaut :
q1l2
M2 0.5M p
24
Avec deux rotules plastiques aux encastrements, la poutre se comporte comme une poutre sur
deux appuis simples. L’augmentation de la charge q au-delà de la valeur q1 entraîne une
augmentation de la courbure et des moments fléchissants en travée où les sections ne sont au
plus que partiellement plastifiées tandis que les moments d’encastrement restent constants, si
on admet un comportement élastique-parfaitement plastique du matériau. Par suite, le rapport
M1/M2 décroît, mais reste tout de même supérieur à 1. Il y’a redistribution des moments.
L’augmentation du moment fléchissant M2 due à une charge supplémentaire q vaut :
ql2
M 2
8
En poursuivant l’augmentation de la charge, il nait une troisième rotule plastique au milieu de
la poutre qui transforme celle-ci en un mécanisme. L’état d’équilibre ultime de la poutre est
atteint, il n’est plus possible d’augmenter encore plus la charge. La charge supplémentaire à
appliquer pour provoquer la ruine de la poutre est déterminée à partir de l’équation
d’équilibre :
ql2 4M p
0.5M p M p q 2
8 l
Puisque M1 M 2 M p , on en déduit la valeur de la charge ultime :
12M p 4M p 16M p
qu
l2 l2 l2
La charge de ruine peut être déterminée à partir de l’équation d’équilibre ultime des moments,
sans avoir recours à un calcul élasto-plastique (méthode pas-à-pas). La condition de
mécanisme donne :
q u l2
M1 M 2 2M p
8
En en déduit :
16M p
qu
l2
On peut remarquer que le calcul plastique pour cette poutre permet d’avoir un gain de 33% de
la charge ultime par rapport à un calcul élastique. La poutre étudiée possède donc un
coefficient de redistribution égal à 1.33. L’hyperstaticité permet une redistribution des efforts
82
le long de la poutre. Pour une structure isostatique, le coefficient de redistribution est toujours
égal à 1.
Comme indiqué sur la Figure 4.9, chaque phase de comportement est caractérisée par un
segment de droite dont la pente décroît au fur et à mesure qu'apparaissent des rotules
plastiques et s'annule lorsque la charge limite est atteinte. La flèche maximale juste avant la
formation des deux rotules plastiques aux encastrements vaut :
M p l2
w
32EI
L’augmentation de la charge q provoquant la ruine de la poutre induit une déformée
supplémentaire égale à :
5 ql 4
2
5 M pl
w
384 EI 96 EI
L’augmentation de la charge q induit une flèche 60% plus grande que celle provoquée par la
charge q1 produisant la formation des deux rotules plastiques, bien que : q = q1/3.
Lorsque le comportement de la poutre demeure dans le domaine élastique, les rotations au
niveau des deux encastrements sont nulles. A l’état ultime, les rotations angulaires
deviennent :
ql3 M p l
24EI 6EI
Si, après avoir atteint la charge limite qu, on effectue une décharge la poutre, les rotations des
rotules plastiques cessent d’augmenter et la poutre se comporte élastiquement. Le moment
résiduel est provoqué par la rotation qui demeure constante durant la décharge. La charge
appliquée étant nulle. Le moment résiduel dans toutes les sections de la poutre est déterminé à
partir de la théorie d’élasticité par la formule suivante :
2EI 2EI M p l
M(x) 0.33M p
l l 6EI
Par opposition à la méthode pas à pas qui consiste à analyser l’évolution des efforts dans
une structure depuis le début jusqu’à la ruine, la méthode de l’équilibre ultime ne fournit que
les efforts internes à l’état ultime, d’où le nom de la méthode.
Les hypothèses principales de la méthode de l’équilibre ultime sont :
- si dans l’une des sections de la structure le moment Mpl est atteint, il naisse une rotule
plastique parfaite ;
- la rotation de la section où se situe la rotule plastique peut augmenter sans restriction et sans
que le moment fléchissant s’accroisse ;
- on néglige les déformations élastiques, ceci conduit à admettre une loi du type rigide-
plastique pour le diagramme moment-courbure des sections droites comme il est représenté
sur la Figure 4.11 ;
- avant la ruine, il ne se produit aucun phénomène d’instabilité dans la structure ou dans l’une
de ses parties.
83
[Link]. Théorèmes fondamentaux
Le calcul à l’état limite ultime des structures doit satisfaire les conditions suivantes :
1) équilibre : les forces extérieures et les réactions appliquées à la structure doivent
s’équilibrer ;
2) résistance : les moments fléchissants ne doivent pas être supérieurs aux valeurs limites (M
Mpl) ;
3) mécanisme : un nombre suffisant de rotules plastiques doit se former afin que la structure
ou une partie de celle-ci puisse se transformer en un mécanisme.
Notons que dans le calcul élastique, la troisième condition est remplacée par la condition de
compatibilité et de continuité des déformations.
Il n’est, en général, pas possible de satisfaire en une seule fois les trois conditions du calcul
plastique. On est obligé de partir de deux de ces conditions et essayer de satisfaire la troisième
par tâtonnement.
Fig. 4.11. Loi moment-courbure du type rigide-plastique
Q M pl
avec :
est le déplacement virtuel du mécanisme de ruine au droit de la force extérieure Q ;
est la rotation de la rotule plastique.
84
borne inférieure (lower bound theorem). Elle donne des valeurs du côté de la sécurité (safe
theorem).
85
- deux mécanismes de nœuds.
Le nombre total de mécanismes à considérer vaut 26- 1= 63.
a) b) c)
EA
L
12 EI
0 3
l
0 0
6 EI
l
2
EA EA
0 0
L 12 EI 6 EI
l
12 EI
0 0
l
3
l
2
l
3
6 EI 2 EI 6 EI 4 EI
0 l
2
l
0
l
2
l
i j
87
EA
L
3EI
0 3
l
0 0
6 EI
l
2
EA EA
0 0
L 3EI
l
3EI
0 0 0
l
3
l
3
3EI 3EI 3EI
0 l
2
0 0
l
2
l
EA
L
0 0
0 0 0
EA 0 0
EA
L l
0 0 0 0 0
0 0 0 0 0 0
En écrivant, à chaque pas dans un système général ces rigidités pour toutes les barres, on
arrive à l’équation matricielle :
P k D
[K] est la matrice carrée dite de rigidité ;
D est la matrice colonne des déplacements aux nœuds;
{P} est la matrice colonne des actions aux nœuds.
88
4.6.3. Méthodes règlementaires
Dans certaines conditions, les normes de construction autorisent une redistribution forfaitaire
et limitée des moments issus d’une analyse élastique. La redistribution des moments consiste
à réduire les moments dans les sections où le rapport entre moment sollicitant et moment
résistant est le plus élevé (en général, aux appuis intermédiaires), et à augmenter les moments
de signes opposés (en général dans la partie courante des travées).
Les conditions à satisfaire pour les poutres en acier selon l’EC3 sont les suivantes :
- les sollicitations dans l'ossature restent en équilibre avec les forces appliquées ;
- toutes les barres dans lesquelles les moments sont ainsi réduits ont des sections transversales
« compactes » de classe 1 ou de classe 2 (voir § 5.6 de l’EN 1993-1-1);
- le déversement des poutres est empêché.
Le diagramme des moments élastiques peut être modifié en redistribuant jusqu’à 15 % du
moment « aux appuis », à condition que les moments et les efforts tranchants obtenus restent
en équilibre avec les charges extérieures appliquées et que les nouveaux moments « en travées
» ne dépassent pas le moment résistant plastique Mpl (voir § 5.4 de l’EN 1993-1-1).
Pour réduire les moments sur appuis ou pour vérifier la résistance structurale des poutres ou
dalles, le BAEL propose deux méthodes : la méthode des moments forfaitaires (§ B.6.2.21 et
Annexe E.1) et la méthode dite Caquot « minoré » (§ B.6.2.210 et Annexe E.2)
L’EN 1992-1-1 autorise un calcul élastique avec redistribution limitée à l'état-limite ultime
uniquement. Les moments déterminés par l’analyse élastique-linéaire, peuvent être
redistribués (jusqu'à un maximum de 30%), sous réserve que la nouvelle distribution des
efforts internes continue à équilibrer les charges appliquées et soit soumise à certaines limites
et certains critères de ductilité, (limitations de la profondeur de l'axe neutre, pour garantir une
capacité de rotation plastique suffisante, par exemple).
Les détails de l’Analyse linéaire élastique avec redistribution limitée des poutres mixtes des
bâtiments sont exposés dans le § 5.4.4 de l’EN 1994-1-1. Les méthodes applicables aux
poutres continues des ponts sont exposées dans le § 5.4.2 de l’EN 1994-1-2.
89
- Les sections transversales de Classe 3 sont celles où la contrainte calculée dans la fibre
comprimée extrême peut atteindre la limite d’élasticité mais où le voilement local empêche le
développement du moment résistant plastique. En effet, le voilement local apparaît «
immédiatement » après la plastification des fibres extrêmes.
- Les sections transversales de Classe 4 sont celles où le voilement local limite le moment
résistant (ou la résistance à la compression pour les éléments sous charges normales). Ce
voilement apparaît bien avant la plastification des fibres extrêmes.
90
Pour la vérification des éléments en béton armé aux ELU, l’EC2 propose trois méthodes pour
tenir compte de la redistribution des moments :
- une analyse élastique-linéaire avec redistribution limitée des moments sans vérification
explicite de la capacité de rotation des sections (Art 5.5 de l’EN 1992-1-1) ;
- une analyse plastique (cinématique ou statique) sans vérification directe de la capacité de
rotation des sections (§ 5.6.2 de l’EN 1992-1-1) à condition de vérifier les conditions de
l’Art 5.6.2(2);
- une analyse plastique avec vérification directe de la capacité de rotation des sections (Art
5.6.3 de l’EN 1992-1-1).
La vérification de la rotation plastique à l'état-limite ultime est considérée comme satisfaite si
l'on montre que, sous l'action considérée, la rotation calculée, s (Figure 4.18), est inférieure
ou égale à la rotation plastique admissible pl.
Fig. 4.18. Rotation plastique s d'une section droite de béton armé [Figure 5.5 de l’EN 1992-
1-1]
La courbure étant supposée constante sur 1,2.h, la rotation plastique est égale à :
1 c s
pl dx dx
r 1.2h d
s est l’allongement plastique des aciers tendus ;
c est le raccourcissement plastique du béton comprimé.
La rotation plastique admissible peut être déterminée en multipliant la valeur de base de la
rotation admissible, qpl,d par un coefficient de correction k qui dépend de l'élancement vis-à-
vis de l'effort tranchant (Art 5.6.3 de l’EN 1992-1-1). Les valeurs de qpl,d recommandées sont
données par la Figure 5.6N de l’EN 1992-1-1.
91
Isolons une tranche de poutre comprise entre les abscisses x et x + dx (Figure 4.19). La
tranche de longueur dx est soumise en ses extrémités aux contraintes et + (/x)dx. Les
contraintes normales doivent donc être équilibrés pat des contraintes tangentes longitudinales
dont l’intensité vaut :
dx
x
y y
92
a)
y y y
max
s max =
- y - y - y
b) Section A3A3 c) Section A2A2 d) Section A1A1
Fig. 4.20. Répartition élasto-plastique des contraintes normales et tangentielles
Cette répartition des contraintes tangentielles se produit dans les sections situées entre A1A1
et A2A2. Dans les sections situées à droite de la section A2A2, des zones de plastification par
cisaillement apparaissent près du centre de gravité de la section transversale en plus des
zones de plastification par flexion dans les fibres extrêmes. La section A3A3 est entièrement
plastifiée par flexion et par cisaillement (Figure 4.20b).
La Figure 4.21 indique trois cas de répartition de contraintes normales et tangentielles
utilisées par les normes de construction pour déterminer les critères de résistance des
sections rectangulaires soumises à un moment de flexion et effort tranchant.
- Cas 1 (Figure 4.21a) :
Soient :
T l’effort tranchant sollicitant ;
Tpl la résistance plastique au cisaillement de la section transversale.
y
T bs
3 T s
y Tpl h
Tpl bh
3
La résistance plastique du moment compte tenu de l’effort tranchant vaut :
1 2 1 1
M v,pl M pl bs y bh 2 y bs 2 y
4 4 4
Mpl est le moment de résistance plastique en absence d’effort tranchant.
93
Il en résulte les critères de plasticité :
2
M v,pl s
1
M pl h
2
M v,pl T
1
M pl T
pl
y
- y
y y / 3
a)
y y
- y
b)
y
yr y
- yr
- y
c)
Fig. 4.21. Répartition possibles des contraintes normales et tangentielles à l’état ultime
- Cas 2 (Figure 4.21b) :
2 y
T bs
3 3 T 2 s 2
y Tpl 3 h 3
Tpl bh
3
94
s
b
1 2 2 s s 1 2 bs 2
M v,pl bh y y bh y y
4 2 6 6 4 12
Il en résulte les critères de plasticité :
2
M v,pl 1 s
1
M pl 3 h
2
M v,pl 3 T
1
M pl 4 Tpl
- Cas 3 (Modèle de Heyman et Dutton, Figure 4.21c) :
T bh T
Tpl y bh Tpl y
1 2
M v,pl bh yr
4
Il en résulte les critères de plasticité :
M v,pl yr 32
1
M pl y 2y
2 2 2
M v,pl T M T
1 1
M pl T M
pl pl Tpl
4.7.2. Section en I
Pour les sections en I, on admet que l’effort tranchant est entièrement repris par l’âme. Le
moment plastique réduit est calculé en appliquant une limite élastique réduite yr sur l'aire de
l’âme Aweb et en considérant, comme dans le modèle de Heyman et Dutton, les répartitions
des contraintes normales et tangentielles représentées sur la Figure 4.22.
T A web
T 3
y
Tpl A web Tpl y
3
L’âme est complètement plastifiée lorsqu’en tout point, le critère de von Mises est vérifié :
2
32 T
3 yr y 1 2 y 1
T
2 2 2
y
yr y
pl
yr y
y yr
Fig. 4.22. Répartition des contraintes avec une limite d’élasticité réduite dans l’âme
97
4.8. Influence de l’effort normal sur le moment plastique
Le moment de plastification d’une section est réduit du fait de la présence d’un effort normal.
L’interaction moment plastique-effort axial appliqué est déterminée pour les sections
transversales entièrement plastifiées. Cependant, les aires des deux côtés de l’axe neutre
plastique ne sont pas égales comme c’est le cas de la flexion simple.
y y
y
y y
N bd y
bh 2 y bd 2 y
M N,pl
4 4
On en déduit :
bh 2 y bd 2 y
M N,pl 4 4 d2
1
M pl bh 2 y h2
4
Ce qui permet d’obtenir le critère de plasticité pour une section rectangulaire :
2
N M N,pl
1
N pl M pl
98
La Figure 4.26 indique les limites du domaine conduisant à la plastification des fibres
extrêmes de la section rectangulaire et le diagramme d’interaction plastique établies dans un
repère orthonormé réduit (N/Npl, M/Mpl). L’utilisation de paramètres adimensionnels permet
de les utiliser quel que soit les unités des moments fléchissants et des efforts normaux.
4.8.2. Section en I
Dans le cas d’une section en I, deux cas sont envisagés.
1. L’effort normal appliqué est inférieur à la résistance plastique de la section de l’âme, soit :
N y t w h 2t f
Dans ce cas, l’axe neutre plastique tombe dans l’âme. Le diagramme des contraintes normales
est décomposé selon les indications de la Figure 4.27.
La hauteur d1 est déterminée à partir de l’équation d’équilibre :
N
N y t w d1 d1
y t w
La résistance plastique de la section transversale vaut :
N pl y A
d’où :
N t w d1
N pl A
Le moment plastique réduit est déterminé à partir des expressions suivantes :
t w d12 t N2 N2 A2
M N,pl M pl y M pl w2 2 y M pl
4t w y 4t w 2y
y
4 A2
2
N A2
M pl y
N pl 4t
w
99
y y
y
y y
M N,pl bh 2 bd 22
M pl 4Wpl
y y
y y y
Fig. 4.28. Axe neutre plastique tombe dans l’âme
100
La Figure 4.29 représente les diagrammes d’interaction plastique MN,pl/Mpl-N/Npl plastique
des profilés HEA300 et IPE 360. Les valeurs 0.25 et 0.40 du rapport N/Npl sont les limites,
respectivement, pour les deux profilés à partir desquelles l’axe neutre plastique tombe dans
une des semelles. On peut remarquer que le moment de résistance plastique est peu influencé
par l’effort normal lorsque le rapport N/ Npl est de l’ordre de 0.5.
N/Npl
MN,pl/Mpl
Fig. 4.29. Diagramme d’interaction plastique des profilés HEA 300 et IPE 360
N
2
M N,Rd M pl,Rd 1 Ed
N pl,Rd
où :
Mpl,Rd est la valeur de calcul de la résistance d'une section transversale à la flexion ;
NEd est la valeur de calcul de l’effort normal ;
Npl,Rd est la valeur de calcul de la résistance d'une section transversale à l’effort normal.
Pour les sections bi-symétriques en I ou H et autres sections bi-symétriques à semelles, il n'est
pas nécessaire de considérer l'incidence de l'effort normal sur le moment résistant plastique
lorsque :
N Ed 0.5h w t w y
101
CHAPITRE 5
5.1. Introduction
Le problème élasto-plastique se formule comme un problème d’évolution non linéaire pour
lequel une stratégie de résolution incrémentale est nécessaire. Supposons qu’à l’instant tn
toutes les grandeurs mécaniques {un, n, n…….} sont connues en un point d’intégration
de la structure. L’approche incrémentale consiste à déterminer ces grandeurs à l’instant tn+1.
La discrétisation temporelle permet de construire une succession de problèmes non linéaires
et indépendants du temps. Les conditions d’admissibilité statique et cinématique sont
uniquement vérifiées aux extrémités des incréments temporels.
En résumé, chaque itération de calcul est composée de deux étapes :
1. Étape globale lors de laquelle les déplacements sont calculés aux nœuds du maillage de la
structure. À l’itération 0, il est supposé que l’incrément des déformations plastiques p est
nul (calcul purement élastique). Aux itérations suivantes, les déformations à une itération
n joue le rôle d’une déformation initiale pour l’itération n+1. Les équations de la
plasticité ne sont pas utilisées durant cette étape.
2. Étape locale qui consiste à intégrer la loi de comportement en chaque point de Gauss de la
structure. Le principe de l’état local considère que le comportement en un point ne dépend
que des variables définies en ce point, et non de son voisinage. Cependant, la résolution
des équations du modèle de comportement intervient dans le contexte de vérification de
l’équilibre globale de la structure. Connaissant l’incrément des déformations totales , on
calcule l’incrément des déformations plastiques p et des variables internes Vk, ainsi que
les contraintes et les forces thermodynamiques associées Ak en utilisant les équations de
la plasticité.
102
Pour les méthodes à un pas, le calcul de la valeur discrète yn+1 au nœud tn+1 fait intervenir la
valeur yn obtenue à l'abscisse précédente.
Pour les méthodes multi-pas, le calcul de la valeur discrète yn+1 au nœud tn+1 fait intervenir
plusieurs valeurs yn ; yn-1; yn-2… obtenues aux abscisses précédentes.
Il existe principalement deux types de schémas d’intégration :
- méthodes d’intégration explicites (Euler-Avant, loi du point milieu généralisée loi
trapézoïdale généralisée etc.) ;
- méthodes d’intégration implicites (Euler-arrière et Schéma de Vermeer).
D’une manière générale, un schéma d’intégration doit assurer trois conditions :
1. La précision de la solution.
2. Le coût de la solution, lié au temps de calcul.
3. La stabilité, étroitement liée à la nature du schéma d’intégration.
Dans ce cours, nous nous limiterons à présenter brièvement les méthodes à un pas :
- méthode d'Euler avant ;
- méthode d'Euler arrière ;
- méthodes de Runge et Kutta ;
- la thêta méthode.
5.2.1. Méthode d’Euler avant
La méthode d’intégration d’Euler avant (forward Euler scheme) est une méthode explicite qui
consiste à approximer une fonction y(t) en connaissant sa dérivée temporelle
dy
f (y t , y) à partir d’une condition initiale y(0) = 0. En connaissant un point yn = y(tn) à
dt
un temps tn, le point yn+1 = y(tn+1) à un temps tn+1 = tn + t peut s’approximer par :
yn 1 yn tf (t n , yn )
103
Comme yn+1 est encore inconnue et écrite sous sa forme implicite, cette méthode d’intégration
est considérée implicite. Dans ce cas, calculer yn+1 revient à résoudre une équation non
linéaire en yn+1, ce qui demande l'emploi d'un algorithme approprié tel que celui de Newton-
Raphson.
y n 1 y n f (t
tn
n , y(t)) dt
Le principe consiste à approcher cette intégrale de façon plus précise que ne le fait la méthode
d'Euler.
Si on approcher l'intégrale par la méthode du rectangle à gauche (Figure 5.2), on obtient le
schéma de la méthode d’Euler explicite, soit :
t n 1
f (t
tn
n , y(t)) dt t f (t n , y n )
f (t
tn
n , y(t)) dt t f (t n 1 , y n 1 )
yn 1 yn tf (t n 1 , yn 1 )
Comme on le voit sur le graphe de la Figure 5.3, la méthode implicite produit une erreur
opposée à celle de la méthode d'Euler explicite.
104
t
105
Cette méthode présente l'avantage d'être précise et assez simple à programmer.
avec 0 1.
Pour chaque pas d’intégration [tn, tn+1], l´évaluation de la fonction f se fait dans un point y(tn +
tθ) à l’intérieur de l’intervalle. Pour ce point on considère le développement en t jusqu’ à
l’ordre 1 :
y t n t 1 yn yn 1 yn (yn 1 yn )
Pour θ = 0, on retrouve le schéma d’Euler explicite.
Pour θ = 1, on retrouve le schéma d’Euler implicite.
Pour θ = 1/2, on obtient la méthode de type Midpoint dite semi implicite :
t yn yn 1
yn 1 yn tf (t n , )
2 2 2
qui est la plus précise, car c’est une méthode d’ordre 2 en t, alors que toutes les autres sont
d’ordre 1.
f n , A k,n 0
L’étape suivante consiste à calculer à l’incrément n+1 les contraintes n 1 , les déformations
plastiques n 1 et les variables d’écrouissages X n 1 et Rn+1.
p
106
5.3.1. Algorithme Euler explicite
Les étapes de calcul selon le schéma d’intégration d’Euler explicite sont résumées ci-après.
1. Calcul des composantes du tenseur de la déformation totale au point d’intégration au temps
tn+1 :
n 1 n B u
2. On ignore si le passage de tn à tn+1 se fait en régime élastique ou en régime plastique. Pour
cette raison, on commence par admettre que la réponse de la structure entre tn et tn+1 est
purement élastique (prédiction élastique). Le champ de contrainte essai, correspondant à
une réponse purement élastique de la structure entre tn et tn+1, est défini par l’équation
suivante :
essai
D : n 1 n D : n n n D :
p p
3. On vérifie si le champ de contraintes est effectivement purement élastique (critère de
plasticité). Deux cas sont possibles :
- 1er Cas : f ( , Ak,n ) 0 , l’état de contrainte est à l'intérieur du critère ou sur le critère; il
essai
plasticité. On doit apporter une correction plastique à l’état de contrainte essai pour
déterminer un état de contrainte qui soit plastiquement admissible.
Dans ce cas, on commence par évaluer l’état de contrainte c situé sur la surface de charge
dans la direction de l’incrément de prédiction élastique (Figure 5.5) :
c n (1 r)
Ak,n représente les paramètres d’écrouissage obtenus à l’instant tn. Les composantes du
tenseur des contraintes c sont supposées croître proportionnellement au seul paramètre
scalaire r.
L’incrément de contrainte re , dirigé vers l’extérieur de la surface de charge, viole le
critère de plasticité. Le champ de contrainte re doit être projeté sur la surface de charge
actualisée. En général, la contrainte calculée en explicite ne vérifie pas exactement f = 0 à
cause de la linéarisation du critère au voisinage de la solution recherchée. Par conséquent, la
contrainte re est appliquée progressivement en réalisant des sous-pas successifs pour
réduire les erreurs de calcul des codes E.F. Le nombre des sous-pas m doit être suffisant
pour obtenir une solution satisfaisante. Selon la publication [5], le nombre nécessaire pour
obtenir la convergence peut atteindre jusqu’à 200. A chaque incrément r / m, l’état de
contrainte corrigé est obtenu par l’expression suivante :
re
( j)
( j)
D: ( j)
( j)
p
m
( j) D : n
( j)
107
est le multiplicateur de charge et n le gradient de f par rapport à et qui ont été défini au §
3.3.3.
Fig. 5.5. Projection par incréments successifs du champ de contrainte sur la surface de charge
L’état de contrainte et les variables internes au temps tn+1 sont alors déterminés à partir des
expressions suivantes :
m r( j)
n 1 n e
( j) D : n ( j)
j1 m
Xn 1 Xn H R n 1 R n dR( )
p p
n 1 n D : n 1 n 1
p
n 1 n n n 1
p p p p
f n 1 , X n 1 , R n 1
n 1 X n 1
f n 1 , X n 1 , R n 1 0
A la différence des schémas explicites où toutes les variables de fin d’incrément sont déduites
sans satisfaire le critère de plasticité, c.-à-d. sans garantir f n 1 , X n 1 , R n 1 0 , la résolution
implicite tient compte de la vérification du critère de plasticité en fin d’incrément.
En écrouissage mixte, la condition de consistance plastique au temps t+t lorsqu’on néglige
l’effet de la température s’écrit :
f n 1 , X n 1 , R n 1 0
On peut à priori résoudre le système non linéaire. Une fois les inconnues : n 1 et n 1 , sont
déterminées toutes les autres quantités seront connues. A ce niveau, plusieurs choix sont
possibles.
108
On peut chercher à résoudre directement ce système par une méthode de type Newton-
Raphson. Cependant, cela présente l’inconvénient d’avoir à inverser la matrice jacobienne à
chaque itération, matrice qui peut être mal conditionnée. De plus, les équations non linéaires à
résoudre sont inter-reliées entre elles et les dérivées nécessaires des équations en questions
dans le calcul de la matrice jacobienne sont souvent difficilement calculables et trop coûteuses
numériquement. Cette complexité rend l’algorithme de la solution lourde et instable. Pour ces
raisons, on a recours à des algorithmes qui surpassent ces difficultés. Les algorithmes de type
Return Mapping est une forme spéciale de la méthode d’Euler et sont souvent utilisé pour
l’intégration des relations constitutives en élastoplasticité.
La Figure 5.6 compare entre la méthode Euler explicite et les méthodes implicites. Dans le cas
de la méthode Euler explicite, la correction plastique s’effectue à partir du point d’intersection
de la prédiction élastique et de la surface de charge. Dans le cas des méthodes implicites, la
correction plastique, et par suite la direction d’écoulement, s’effectue à partir de la prédiction
élastique vers la surface de charge.
Fig. 5.6. Comparaison entre la méthode Euler explicite et les méthodes implicites
On effectue un test d’essai, fondé sur une réponse élastique du milieu, qui ne dépend que des
conditions initiales à l’instant tn. Cet état d’essai est défini par le tenseur de déformation
élastique et la déformation plastique cumulée :
n 1 n
essai e e
n 1 n
P essai p
L’incrément de contrainte d’essai pour un milieu isotrope est déterminé à partir de la loi de
109
Hooke généralisée :
ij tr 1 2 kk ij 2 ij
e
Si le critère de plastification est vérifié, les variables sont actualisées avec les valeurs
prédites :
n 1 n 1 n
P essai p p
Vk n 1 Vk n 1 Vk n
essai
n 1 n 1 n
essai
f essai
n 1 ,essai A n 1 0
La contrainte élastique prédite doit être projetée, avec un correcteur plastique, sur la surface
de charge actualisée, dans le cas d’un comportement élasto-plastique avec écrouissage, ou sur
la surface de charge initiale, dans le cas d’un comportement élasto-plastique parfait, (Figure
5.7). Cette projection est normale à la surface de charge finale.
La correction plastique passe par la résolution du système d’équations non linéaire suivant :
n 1
e essai
n 1 n n 1 , A n 1 0
e
Vn 1 essai
Vn 1 H n 1 , A n 1 0
f n 1 , A n 1 0
Différents méthodes de projection de type «Return Mapping» ont été proposés en plus des
méthodes standards de Newton-Raphson ou de Quasi-Newton. . Parmi ces méthodes on peut
citer :
- la méthode du retour radial ;
- la méthode des trapèzes généralisée ;
- la méthode de return mapping généralisée ;
110
- la méthode du point milieu généralisée
- la méthode de prédiction élastique retour sécant la correction est réalisée suivant la
contrainte moyenne ;
- la méthode de prédiction tangente retour radial la prévision provient cette fois de
l'intégration des équations d'évolution par la méthode d'Euler explicite.
a) b)
Fig. 5.7. Schéma d’intégration de type «Return Mapping» : a) comportement élasto-plastique
avec écrouissage b) comportement élasto-plastique parfait
2 s n 1
pn 1 pn dp pn
f n 1 , R n 1 3J 2 (s n 1 ) R n 1 (p n 1 ) 0
Après résolution du système non linéaire, le tenseur des déformations plastiques mis à jour
s’écrit :
3 s n 1
n 1 n
p p
2 s n 1
où :
s n 1 2Gd n 1
e
111
Le système d’équation non linéaire peut être réduit à la résolution d’une seule équation non
linéaire, ce qui simplifie le problème comme il sera montré ci-dessous.
Compte tenu de l’incompressibilité des déformations plastiques (tr 0) induite par le
p
critère de von Mises, la partie sphérique et la partie déviatorique du tenseur des déformations
sont :
3
p essai
s n 1
2 essai
eq
n 1
3 s n 1
d n 1 d n 1
e essai e
2 s n 1
La décomposition du tenseur des contraintes donne :
essai
p p
n 1 n 1
3 s n 1
s n 1 s n 1 2G
essai
2 s n 1
Soit :
3 2G
1 s n 1
essai
s n 1
2 s n 1
et que les tenseurs déviatoriques d’essai et le tenseur mis à jour coïncident. On peut ainsi
déterminer le tenseur déviatorique mis à jour en fonction du tenseur déviatorique d’essai :
3 2G essai 3G essai
s n 1 1 s 1 s n 1
2 essai s n 1 essai
n 1 (n 1) eq
3 essai
où : essai (n 1) eq s n 1 est la contrainte équivalente d’essai de von Mises (voir §
2
3.1.3).
3G
La quantité intervient ici comme un correcteur plastique c’est la correction qui
essai
(n 1) eq
doit être apportée au tenseur déviatorique calculé à partir de la prédiction élastique pour
atteindre un état plastiquement admissible telle qu’il est illustré sur la Figure 5.8.
112
a) Espace des contraintes principales b) plan déviatorique
Fig. 5.8. Interprétation géométrique du schéma d’intégration «Retour Radial»
Le multiplicateur plastique est déterminé à partir de la condition de consistance :
n 1
f n 1 , R n 1 eq R n 1 (p) y
3
s : s R n 1 (p) y 0
2 n 1 n 1
En écrivant les relations suivantes :
essai
s n 1
s n 1 essai
(n 1) eq 3G essai
(n 1) eq
essai essai
s n 1 : s n 1 2
essai
2
(n 1) eq 3
3
p essai
s n 1
2 essai
eq
n 1
Les autres variables d’état sont ensuite calculées par les expressions suivantes :
3G essai
s n 1 1 s n 1
essai
(n 1) eq
113
n 1 s n 1
essai
p
n 1
v n 1
essai e
s n 1
n 1 D : n 1
1
n 1
e e
2G 3
3 s n 1
n 1 n
p p
2 s n 1
pn 1 pn
Dans le cas d’un écrouissage isotrope linéaire :
R n 1 y ph
La condition de consistance s’écrit :
essai n 1
eq 3G y p h 0
On en déduit l’incrément du multiplicateur plastique :
essai
(n 1)eq y ph
essai
f essai n 1 , R n 1
3G h 3G h
Dans le cas d’un comportement élasto-plastique parfait (h = 0), l’incrément du multiplicateur
plastique devient :
(i 1)
essai
y
essai
f essai n 1 , y
( n1) ( n1) D (i)
3G 3G
114
La vitesse de déformation plastique cumulée pour la plasticité de von Mises vaut (voir §
[Link]) :
p(in 11) p(i)
n 1 p n
i i
1/2
2 f f
n (i)
n 1 : (i)
n 1
3
R
Notons aussi que : (i)
n 1 h (i)
n 1 coefficient d’écrouissage.
p
R
En injectant les expressions de , p i et
(i)
dans la condition de cohérence et en
p
linéarisant la fonction de charge, vient l’expression du multiplicateur de charge :
f n(i)1
i
f (i) f
n 1 : D :
(i)
n 1 h (i)
n 1n
Le multiplicateur plastique permet de calculer les variables qui caractérisent l’état plastique :
f
n 1 n 1 D :
i 1 i i (i)
n 1
R n 1 R n 1 hdp
i 1 i
Les itérations vont être poursuivies jusqu’à ce que la fonction de charge soit suffisamment
proche de zéro.
f n(i)1 tolérance 0
L’interprétation géométrique de cet algorithme est représentée sur la Figure 5.9 dans le
cas d’un comportement plastique avec écrouissage.
n 1 n 1 H n H n 1
e
f n 1 , A n 1 0
115
Fig. 5.9. Interprétation géométrique du schéma d’intégration «Cutting plane»
Le paramètre représenté sur la Figure 5.10 détermine la part du retour suivant la normale
initiale et finale.
0 1
Ce schéma correspond donc à une discrétisation temporelle de type semi-implicite. Le schéma
est purement implicite si = 1 et purement explicite si = 0 (donc pas de nécessité de
méthode itérative).
L’état de contrainte de la prédiction élastique essai
n 1
est ramené sur la surface de charge en fin
d’incrément par la projection suivant la direction de la normale déterminée par l’expression :
1 n n n n 1
C’est une combinaison linéaire des normales à la surface de charge aux temps tn et tn+1.
L’équation d’écoulement plastique correspondante s’écrit :
n 1 n 1 n n n n 1
p p
n 1 D 1 n n n n 1
n 1 essai
116
Fig. 5.10. Méthode des trapèzes généralisés
n 1 n Hn
f n 1 , A n 1 0
Les lois d’état au point milieu sont :
n (1 )n 1 n
An (1 )An An 1
L’état de contrainte final en fonction de l’état de contrainte prédit est défini par l’expression :
n 1 n 1 D : n n
essai
Il est possible de remplacer la condition de cohérence f n 1 , A n 1 0 par :
f n , A n 0
117
L’interprétation géométrique de l’algorithme du point milieu généralisé est représentée sur la
Figure 5.11.
La prédiction élastique peut conduire à l’un des états suivants : état élastique (Figure 5.12a),
état sur la frontière du domaine élastique (Figure 5.12b), état plastiquement inadmissible
(Figures 5.12c et d) ou une décharge élastique (Figure 5.12e).
Dans le cas où l’état prédit est plastiquement inadmissible, la projection sur le critère conduit
à une évolution de la limite élastique en fonction de la déformation plastique obtenue au
temps tn :
ny 0y h np
118
a) b) c) d) e)
Fig. 5.12. Illustration des cas de prédiction élastique
où 0y est la limite élastique initiale ;
h est le module d’écrouissage isotrope (voir § 2.2.2) :
ET
h
E
1 T
E
L’incrément de contrainte d’essai essai est décomposé en une partie élastique Ressai et
une partie inadmissible (1-R) essai (Figure 5.14).
La fonction de charge ou le critère s’écrit :
f essai essai ny n ny (1 R)E
119
essai ny (E h)p 0
d’où l’on tire :
essai ny f essai E
p
(1 R)
(E h) Eh Eh
A noter que : p 0
pn 1 np p
120
Ecrouissage cinématique
Considérons un état de contrainte en un point (M) vérifiant le critère de plasticité dont on
connait à l’instant tn : n, n, Xn, y0.
Xn étant la contrainte interne (back stress) à l’instant tn.
On applique un incrément de contrainte élastique essai E. Il en résulte :
=-X
=X
121
d’où l’on tire :
essai X essai 0y f essai
p
(E H) EH
Ce qui permet de mettre à jour les contraintes :
n 1 essai signe(essai X essai )E p
X n 1 X n signe(essai X essai )H p
122
Exemples d’applications
Exemple 2.1
On réalise un essai de traction sur une éprouvette d’un matériau de section circulaire (D = 6
mm, l0 = 100 mm). Les données expérimentales enregistrées au cours de cet essai de traction
dont on a besoin pour cet essai sont comme suit :
• sous une force F1 de 5 kN, la distance entre les points de référence est égale à 100,6 mm ;
• si la force est supprimée, la distance entre les points de référence est égale à 100,2 mm ;
• au maximum de la courbe brute de traction, la force Fmax est égale à 7 kN et la distance entre
les points de référence est égale à 112 mm ;
• la rupture se produit sous une force Ff de 6,5 kN et, après la rupture de l’éprouvette, la
distance entre les points de référence (obtenue en mettant bout-à-bout les deux morceaux
d’éprouvette) est égale à 114 mm.
À partir de ces données, on demande de caractériser le comportement en traction de ce
matériau et de calculer :
a) sa limite conventionnelle d’élasticité Re0,2 (en MPa).
b) son module de Young E (en GPa).
c) sa résistance à la traction Rm (en MPa).
d) son allongement après rupture A (en %).
e) la déformation AM (en %) qui est la borne supérieure du domaine de déformation plastique
homogène de ce matériau ;
f) l’énergie élastique Wél (en kJ/m3) qui était emmagasinée par unité de volume de matériau
juste avant sa rupture.
123
Force F (kN) Contrainte (MPa) Déformation (%)
5 176.8 Sous charge 0.6
Après décharge 0.2
7 247.6 Sous charge 1.2
6.5 229.9 Après décharge 1.4
a) Limite d’élasticité
Par définition, la limite conventionnelle d’élasticité à 0,2% (Re0,2) correspond à la contrainte
qui, après suppression, a causé une déformation ε permanente de 0,2% du matériau. Dans le
cas présent, cette limite d’élasticité Re0,2 correspond donc au point A de la courbe de traction.
Re0,2 = 176.8 MPa
b) Module d’Young du matériau
Le module de Young E du matériau est égal à la pente de la droite AC de la Figure ci-dessus.
Les dimensions du triangle ABC permet de déduire la valeur de E :
176.8
E 44200 MPa 44.2 GPa
(0.6 0.2) 10 2
c) Résistance à la traction
Par définition, la résistance à la traction Rm du matériau est la contrainte maximale σmax
atteinte durant l’essai de traction (point D de la courbe de traction).
Rm = 247.6 MPa
d) Allongement A (%) après rupture du matériau
Cet allongement après rupture est égal à la déformation permanente résiduelle du matériau
une fois qu’il est rompu. Sur la courbe de traction, cette déformation est caractérisée par le
point H, soit :
114 100
A(%) 14 %
l00
124
e) Déformation Am (%) délimitant la déformation plastique homogène du matériau
Cette déformation Am correspond à celle atteinte au maximum de la courbe de traction, alors
que le phénomène de striction (déformation inhomogène localisée) va apparaître au cours de
l’essai de traction. Cette déformation Am est caractérisée par le point E de la courbe de
traction :
Am(%) = 12%.
f) Énergie élastique Wel emmagasinée par unité de volume de matériau juste avant
rupture
Par définition, cette énergie est égale à la surface du triangle HGF de la courbe de traction,
soit :
1
Wel f elf
2
avec : f : contrainte à la rupture = 229,9 MPa (point F de la courbe) ;
εelf : déformation élastique du matériau juste avant rupture = longueur du segment HG.
Comme la valeur de la déformation totale au point G n’est pas connue, on applique la loi de
Hooke pour trouver la valeur de εelf, puisque la droite HF a pour pente le module de Young :
f
elf
E
d’où :
2
f2 229.6 106
Wel 598 103 N / m 2 598 103 J / m3
2E 2 44.2 10 9
598 kJ / m3
1 Joule = 1 N.m
Exemple 2.2
Pour un alliage d’aluminium, on effectue un essai de traction sur une éprouvette prélevée
d’une plaque de 3 mm d’épaisseur. L’installation de l’éprouvette dans les mors de la machine
de traction et les paramètres du diagramme contrainte-déformation sont illustrés sur la Figure
ci-dessous. Dans la zone utile, l'éprouvette a une largeur L1 = 20 mm et les rayons de
courbure de l’éprouvette sont suffisamment grands pour négliger l’effet de concentration de
contraintes. Les propriétés de l’alliage sont les suivantes :
E = 71 GPa Re = 310 MPa Re0,2 = 330 MPa Rm = 450 MPa A = 18%
a) Pour quelle force appliquée F (en kN) aura-t-on le début de la déformation plastique dans la
zone utile de l'éprouvette ?
b) Quelle est la valeur de la déformation élastique el (en %) à la limite conventionnelle
d’élasticité Re0,2 ?
c) Sachant que la longueur initiale AB, sur laquelle est mesuré l'allongement de l'éprouvette,
est égale à 10 cm et que la contrainte f à la rupture est égale à 400 MPa, quelle sera la
valeur de AB (en cm) juste avant la rupture de l'éprouvette ?
d) Quelle valeur minimale doit avoir la largeur L2 afin d'éviter qu'une déformation plastique
n'apparaisse dans les régions saisies par les mors au cours de l'essai complet de traction ?
125
a) Force F au début de la déformation plastique
Cette force correspond à celle appliquée quand on atteint la limite proportionnelle d’élasticité
Re de l’alliage.
Fel S0 R e 20 3 310 18600 N 18.6 kN
Notons que la limite de proportionnelle d’élasticité ne correspond pas toujours à la limite
élastique
b) Déformation élastique εel à Re0,2
Pour toute valeur de la contrainte, qu’elle soit située dans le domaine élastique ou dans le
domaine plastique, on aura une déformation élastique εel associée qui est donnée par la loi de
Hooke : εel = σ/E. Donc ici nous aurons :
R e,0.2 330
el 0.465%
E 71000
c) Longueur A’B’ juste avant la rupture
L’allongement final A après rupture est égal à 18%. Pour connaître la déformation totale εt
juste avant la rupture, il faut lui ajouter la déformation élastique εelr qui existe dans le
matériau avant sa rupture, donc :
t A elr A f / E 0.18 400 / 71000 18.56 %
La longueur A’B’ sera donc égale à (1 + t)AB = 1,1856 AB = 1,1856 x 10 cm = 11,856 cm.
d) Valeur minimale de la largeur L2
Au cours de l’essai de traction, la force maximale appliquée Fmax à l’éprouvette sera celle
correspondant à Rm. La valeur maximale atteinte par la contrainte vaut :
Fm R m S0 R meL1
où : S0 = e.L1 est la section initiale de la partie utile de l’éprouvette.
Pour qu’il n’y ait pas de déformation plastique dans les têtes de l’éprouvette, il faut que la
contrainte qui y est appliquée ne dépasse pas la limite proportionnelle d’élasticité Re du
126
matériau. La section des têtes S = e.L2 (où e = 3 mm est l’épaisseur de l’éprouvette) doit donc
être telle que :
Fm F
Re m
S e L2
d’où
Rm 450
L 2 L1 20 29.03mm
Re 310
Exemple 2.3
Un essai de traction a été réalisé sur une éprouvette en acier afin de déterminer les
caractéristiques mécaniques et la courbe contrainte-déformation pour en déterminer la loi de
comportement expérimentale. Le déplacement de la traverse se déroule à vitesse constante
(2mm/s). L’essai est réalisé sur une éprouvette « plate » en acier de section rectangulaire (1.02
mm x 8.10 mm).
S0 = 8.26 mm2
A partir de l’enregistrement dans la zone de plasticité, nous allons établir le diagramme
rationnel de la forme : σ = f(ε), soit la contrainte vraie en fonction de la déformation vraie
(voir § [Link]). Ce diagramme nous permettra ensuite d’établir la loi de comportement
plastique du matériau. La portion parabolique de la courbe déformation plastique vraie-
contrainte vraie sera représentée par la loi de Hollomon (voir § [Link]) :
n
vh k pv
F
v 1
S0
pv Ln 1 p
avec :
pv est la déformation plastique rationalisée : pv v e v
E
Le paramètre n représente le coefficient d’écrouissage (strain-hardening exponent), il rend
compte de la capacité du matériau à se déformer plastiquement.
Le paramètre k représente le coefficient de résistance du matériau (strength coefficient
exprimé en MPa).
Le Tableau ci-dessous donne les valeurs des déformations vraies et des contraintes vraies
obtenues à partir des valeurs bruts de l’essai.
Ln v Ln K n Ln pv
n
v k pv
127
donnera les paramètres n et K.
Ln(266) = b + aLn(0.0227) 5.58 = b -3.78a
Ln(364) = b + aLn(0.11) 5.9 = b -2.21a
5.9 5.58
a n 0.20
3.78 2.21
Ln(k) 6.336 k e6.336 565 MPa
0.20
v 565 pv
La Figure ci-dessous compare entre les courbes contrainte-déformation plastique
expérimentale, rationnelle et celle de la loi de Hollomon. On constate que l’écart entre les
courbes nominale et rationnelle augmente au fur et à mesure que la déformation augmente.
380
350
320
Contrainte (MPa)
200
0,0000 0,0300 0,0600 0,0900 0,1200
Déformation plastique
128
Exemple 2.4
Soit la structure à trois barreaux articulés à leurs deux extrémités ayant tous la même section
(S1 = S2 = S3 = 400 mm2), de longueur L1 = 2 m et L2 = L3 = 4 m et sont réalisés avec l’acier
de nuance S235. La charge F est supposée appliquée lentement afin que le problème soit
quasi-statique. Nous considérons deux types de comportement :
- comportement élasto-plastique parfait : E = 200000 MPa y = 235 MPa ;
- comportement élasto-plastique avec écrouissage linéaire : la courbe d’écrouissage est bi-
linéaire et caractérisée par les paramètres suivants : u = 360 MPa u = 26%.
Déterminer les efforts dans les barres pour les deux lois de comportement.
2m
2m S3
Les efforts normaux dans les barres sont les inconnues du problème.
Nombre d’inconnues = 3.
Nombre d’équations = 2.
Le système est une fois hyperstatique.
I) Comportement élasto-plastique parfait
a) Phase élastique (F Fy1)
N1 N 2 N 3 F
Équations d’équilibre :
N 2 N3
N1l1
u1 ES
Déplacements :
u u N 2 l 2 N 3 l3
2 3
ES ES
La symétrie du système impose la condition de compatibilité suivante : u1 = u2 = u3, d’où :
N1 2N 2 2N3
F F
N1 2N 2 4N 2 F N 2 N 3 et N1
4 2
C’est donc le barreau 1 qui se plastifie en premier lorsque la force appliquée a pour valeur :
Fy1
1 y 235 MPa
2S1
129
On en déduit :
Fy1 188000 N
N1 94000 N
N 2 N 3 47000 N
u 188000 2000 2.35 mm
2 2 105 400
N2 N3 2N2 yS F
282000 4000
u 4.7 mm
2 2 105 400
La Figure ci-dessous montre l’évolution de la force appliquée, des efforts normaux et des
déplacements durant les phases élastique et élasto-plastique.
300
282
250
200
188
Force F (kN)
150 F
N1
94
100 N2
50
2,35 4,7
0
0 1 2 3 4 5
Déplacement u (mm)
130
1
En fin de phase élastique nous avons : 2
2
b) Phase de comportement élasto-plastique
Pour des valeurs de la charge F supérieures au chargement limite Fy1, le barreau 1 subit des
déformations plastiques et une contrainte supérieure à y. La contrainte dans le barreau 1 peut
être exprimée sous la forme :
1 'y ET 1
avec (voir § 2.2.2) :
360 235
ET 483MPa
235
0.26
200000
y y'
E
E T 'y y 1 T y 235 MPa
y E
E
1
y ET
𝜎𝑦′
E
1
Après plastification, l’effort normal N1 s’écrit sous la forme :
N1 1S y E T 1 S yS E T 1S yS 2E T 2S
2E T N 2
yS
E
2N 2
1 2 2
ES
L’équation d’équilibre donne :
2E N
F yS T 2
F N1 E
N2
2 2
D’où l’on tire :
2E
F yS 1 T N2
E
les barreaux 2 et 3 se plastifient lorsque :
E
N 2 yS 94000 N Fy2 2 yS 1 T yS 285610 N
E
131
N1 285610 2 94000 97610 N
3. Phase de plasticité
L’écrouissage permet de poursuivre le chargement pour des valeurs de F supérieurs à Fy2.
L’évolution est alors complètement plastique.
L’état ultime est atteint lorsque : 1 = 0.26 et 1 = 360 MPa, 2 = 0.13, d’où :
N1 1S 144000 N N2 N3 yS ETS2 119116 N
La charge ultime est déterminée à partir de l’équation d’équilibre :
Fu N1 2N 2 382232N
La Figure ci-dessous montre la variation des efforts normaux dans les barres en fonction de la
force F.
160
140
120
100
N1 (kN)
80
60
ETET=0
=0
40
ETET=483
= 483 MPa
MPa
20
0
0 100 200 300 400
Force F (kN)
120
100
80
N2 (kN)
60
40
20 ET =483 MPa
ET =0
0
0 100 200 300 400
Force F (kN)
On peut remarquer que le comportement élasto-plastique avec écrouissage n’a d’influence que
lorsque toutes les barres sont plastifiées. Pour la structure étudiée, la considération d’une loi
de comportement élasto-plastique avec écrouissage et un module tangent ET = 483 MPa fait
augmenter la charge ultime de 35% par rapport à un calcul avec un comportement élasto-
plastique parfait.
Exemple 3.1
Considérons l’état de contrainte plan en un point M d’un élément en bronze dont la contrainte
limite est u = y = 345 MPa et dont les valeurs des contraintes normales et tangentielles sont
indiquées sur la Figure ci-dessous.
132
Vérifier la résistance du matériau avec les critères de Rankine, de Tresca et de Von mises?
80 125 0
125 190 0
0 0 0
Le signe (-) indique une compression
Les contraintes principales ont pour valeurs :
133
La condition de résistance n’est donc pas vérifiée par le critère de Tresca.
Critère de von Mises (voir § 3.1.3)
Le critère de von Mises s’écrit :
eq (239) 2
(129)2 (239)(129)
1
190 80 80 190 6 125 323MPa
2 2 2 2
2
Le coefficient de sécurité FS vaut :
345
FS 1.067 1
323
La condition de résistance est vérifiée par le critère von Mises. Le critère de Tresca est plus
sévère que le critère de von Mises.
Exemple 3.2
Considérons l’état de contrainte plan en un point M d’un élément en fonte dont les
contraintes limites en compression et en traction sont respectivement : uc = 970 MPa et ut = 295
MPa. Les valeurs des contraintes normales et tangentielles sont indiquées sur la Figure Ci-dessous.
Vérifier la résistance du matériau par le critère de Mohr-coulomb simplifié.
Le critère de Mohr-Coulomb simplifié peut être utilisé pour les matériaux ductiles ayant les
limites élastiques différentes en compression et en traction. Ceci n’est pas le cas des critères
de Tresca et von Mises qui sont des critères isotropes. Le critère de résistance s’écrit :
UT
FS 1 si 1 et 2 sont des contraintes de traction ;
max(1 , 2 )
UC
FS 1 si 1 et 2 sont des contraintes de compression ;
max( 1 , 2 )
UT UC
FS 1 si 1 et 2 sont de signes opposés.
1 2
134
1
1 2
150 100
2
150 100
50 159.6 MPa
2
1
2 2
150 100
2
150 100
50 109.6 MPa
2
2
2 2
La représentation graphique de Mohr de l’état de contrainte est montré sur la Figure ci-
dessous.
135
Les contraintes extrêmes dans toutes les sections du barreau ont pour valeurs :
P 100000 Mr 2300000 30
x 35.36 MPa x 108.46 MPa
A 60 2
I
30
4
4 4
Tr 5000000 30
xy 117.9 MPa x 35.36 108.46 143.8 MPa
J
30
4
2
Contraintes principales :
2
143.8 143.8
117.9 210 MPa
2
1
2 2
2
143.8 143.8
117.9 66.2 MPa
2
2
2 2
- Critère de Tresca :
1 2 210 66.2
max 138.1 MPa
2 2
300
FS 2 1.08 1
138.1
La condition de résistance est vérifiée par le critère de Tresca.
- Critère de von Mises :
eq (210) 2
(66.2 210) (66.2) 2 (143.8) 2
3(117.9 ) 2 250 MPa
300
FS 1.2 1
249.76
La condition de résistance est vérifiée par le critère de von Mises.
Les résultats obtenus sont représentés graphiquement sur la Figure ci-dessous.
136
2
210 MPa
1
66.2 MPa
Exemple 3.4
L’état de contrainte en un point M est défini dans la base e1 , e 2 , e3 par le tenseur suivant :
40 20 10
20 80 5 MPa
10 5 60
40 20 10
I3 ij Det(ij ) 20 80
5 11 22 33 212 23 31 11 23
2
22 31
2
33 122
10 5 60
40 80 60 +2 20 (5) 10 – 40 52 (80) 102 60 (20) 2
211000 (MPa )3
Les racines de l’équation du troisième degré représentent les contraintes principales au point
M:
3 20 2 1085 227000 0
La solution de l’équation du troisième degré fournie les contraintes principales :
1 64.21MPa 2 39.33MPa 3 83.54 MPa
137
- Contrainte équivalente de Tresca :
eq 1 3 64.21 83.54 147.75MPa
eq
2
40 80 80 60 40 60
2 2 2
6 (20) 2 (5) 2 (10) 2
2
137 MPa
Exemple 3.5
Dans le plan traction-cisaillement, l’état de contrainte en un point M est défini par le tenseur
suivant :
0
0 0
0 0 0
Démontrer que dans ce plan, le critère de Tresca et le critère de von Mises ont pour
expressions :
Tresca : 2 4 2 y
von Mises : 2 32 y
3p p2 2 p p p2 p 2 0
Les contraintes principales sont :
2 4 2 2 4 2
1 2 0 3 0
2 2
1 3 y 2 42 y
- Critère de von Mises
Dans la base principale le critère de von Mises s’écrit :
138
f 1 2 2 3 1 3 22y 0
2 2 2
2 2
2 4 2 2 42
2
2
2 4 2
2 2 y
2 32 y
Sous cette forme, le critère de von Mises est utilisé dans l’EC3 pour le calcul de la résistance
plastique des sections à l’effort tranchant ou au moment de flexion en présence de l’effort
tranchant (voir les clauses 6.2.6 et 6.2.8 de l’EC3).
Dans une base quelconque le critère de von Mises s’écrit :
f 11 22 22 33 11 33 6 12
2 2 2 2
23
2
31
2
2 2y 0
f 2 2 62 2 2y 0
2 32 y
Exemple 3.6
Soit une barre fabriquée à partir d’un matériau ayant les propriétés mécaniques suivantes :
limite élastique : y = 400 MPa ; coefficient d’écrouissage h = 200 MPa. La barre subit une
déformation plastique 11p
0.1 .
On demande de calculer la contrainte normale dans la barre dans les cas suivants :
- écrouissage isotrope linéaire ;
- écrouissage cinématique linéaire ;
- écrouissage mixte.
Le critère de von Mises dans le cas d’un comportement plastique sans écrouissage peut
s’écrire sous la forme suivante (voir § 3.1.3) :
2
sij y 0
3
Pour un état de contrainte unidimensionnel :
2
3 0 0
0 0
ij 0 0 0 sij 0 0
3
0 0 0
0
0
3
4 2 1 2 1 2 2
sij
9 9 9 3
Déformation plastique cumulée :
p 11
p
0.1
139
Le critère de plasticité de von Mises s’écrit :
s
2
3
y R s
2
3
y hp 2
3
2
3
(400 200 0.1) 0 420 MPa
2 2 2 2
hp y 1 hp 0
3 3 3 3
y hp 420MPa 0 420MPa
Exemple 3.7
Le comportement élasto-plastique du matériau d’une structure en acier est représenté par la
courbe rationnelle illustrée sur la Figure ci-dessous, le critère de plasticité de von Mises et la
loi d’écoulement plastique associée. Les propriétés élastiques sont caractérisées par : le
module de Young E = 200 GPa et le coefficient de Poisson = 0.3.
Déterminer les états de déformations dans les cas suivants :
- état de contrainte uniaxiale 11= ;
- état plan de déformation 11et 22 0 et 22 = 0.
550 550
400
400
p
0.002 0.25 0.25-550/200000 = 024725
140
1) Etat de contrainte unidimensionnel
- Tenseur des contraintes et des contraintes déviatoriques :
2
3 0 0
0 0
ij 0 0 0 sij 0 0
3
0 0 0
0
0
3
- Comportement élastique
eq = 11= y0 = 400 MPa
Tenseur des déformations élastiques :
E 0 0
ij 0
e
0
E
0 à
E
i, j: ijp 0
- Comportement élasto-plastique avec écrouissage
eq = 11= > y0 = 400 MPa
Dans le cas d’un état de contrainte uniaxiale, la déformation cumulée p est égale à la
déformation longitudinale, p 11
p
.
Le module d’écrouissage vaut :
d 150
h = 606.673MPa
dp 0.24725
L’évolution de la limite élastique en fonction de la déformation plastique est déterminée à
partir de la relation :
d
dp
h
d’où :
d 400
p
400 h 606.673
141
Sous chargement uniaxial, le tenseur vitesse de déformation plastique s’écrit : (voir §
[Link]) :
0 0
p 0 0 h
ij 0 0.5p 0 0 0.5 0
p
h
0 0.5p
0
0 0 0.5
h
Le tenseur de déformation plastique :
400
606.673 0 0
400
ij 0
p
0
1213.346
0 400
0
1213.346
142
En remplaçant 22 par 11, la résolution de l’équation du second degré donne :
400
11 450 MPa
1 2
soit :
22 0.311 135MPa
Pour un comportement plastique à écrouissage isotrope (eq > 400 MPa), la condition de
cohérence s’écrit :
f f
dij dR 0
ij R
que l’on peut écrire sous la forme suivante :
3 s11 3 s 22
11 22 hp 0
2 eq 2 eq
On en déduit :
p
211 22 11 222 11 22
2heq
Par ailleurs, la règle de partition s’écrit :
22 f
22 e22 22
p
0 11 p 0
E E 22
En arrangeant les deux équations précédentes, on obtient :
E 11 222 211 22 4heq
2
22 11
E 11 222 4heq
2 2
Une fois le tenseur des contraintes et de l’incrément de déformation cumulée déterminés, les
composantes du tenseur taux de déformation sont déterminées à l’aide de la relation :
3sij
ijp p
2eq
Les contraintes et les déformations peuvent être déterminées en effectuant un calcul
incrémental en plusieurs pas.
(i 1)
(i)
ij ij ij
ij (i 1) ij ij
(i)
Dans le Tableau ci-dessous on se limite à présenter les résultats des calculs relatifs aux
contraintes et aux déformations obtenus aux trois premiers pas et aux trois derniers pas.
L’incrément 11 doit être choisi suffisamment faible de manière à minimiser les erreurs de
calcul. Pour cet exemple nous avons adopté une valeur de 11 = 0.5 MPa. Les calculs sont
réalisés avec Excel.
Pas ′i′ 11
(i-1)
22
(i-1)
eq
(i-1)
22 11
(i)
22
(i)
eq
(i)
p 105 (i)
s11 (i)
s22 (i)
s33 11
p
105 p22 105 33
p
105
5,49882 254,66 - 58,82 - 195,84 5,25131 - 1,213 - 4,0383
5,70501 254,31 - 57,63 - 196,69 5,44034 - 1,2328 - 4,2076
5,92676 253,95 - 56,41 - 197,55 5,64335 - 1,2535 - 4,3899
(i)
11
p (i)
p22 (i)
33
p
400 0,15
ep11
300
0,1
200
100 0,05
0 0
0 100 200 300 400 500 600 0 100200300400500600700
11 (MPa) 11 (MPa)
-0,00035
ep22
-0,00045 -0,15
-0,00055 -0,2
-0,00065
-0,00075 -0,25
11 (MPa) 11 (MPa)
144
Exemple 3.8
Dans le plan des contraintes principales, l’état de contrainte à l’instant tn en un point M est
défini par les tenseurs de contrainte et de déformation suivants :
151.8 0 0 0.0004 0 0
ij 0 496.8 0 MPa
ij 0 0.00189 0
0 0 276.8 0 0 0.00041
Les propriétés élastiques sont caractérisées par : le module de Young E = 200 GPa et le
coefficient de Poisson = 0.3.
Le comportement élasto-plastique du matériau est représenté par la courbe contrainte-
déformation plastique illustrée sur la Figure ci-dessous, le critère de plasticité de von Mises et
la loi d’écoulement plastique associée.
500 300
Le coefficient d’écrouissage : h 1000 MPa 0.005E
0.2
eq
500 MPa
300 MPa
p p
0 0.2
Déterminer le tenseur des déformations résultant d’un accroissement des contraintes défini
par le tenseur suivant :
68.2 0 0
ij 0 223.2 0 MPa
0 0 124
ijp 0
145
Le tenseur des contraintes après accroissement :
226.67 0 0
sij 0 273.33 0 MPa
0 0 46.67
2 eq 2 eq
226.67 0 0 0.1054 0 0
3
ijp
0.136 0 273.33 0 0 0.1271 0
2 438.6
0 0 46.67 0 0 0.0217
146
1 1
e22 22 11 33 223.2 0.3 68.2 124 0.00083
E 200000
1 1
33
e
33 11 22 124 0.3 68.2 223.2 0.00018
E 200000
0.00018 0 0 0.1054 0 0
ij 0 0.00083 0 0 0.1271 0
0 0 0.00018 0 0 0.0217
0.1056 0 0
0 0.1280 0
0 0 0.0215
0.0004 0 0 0.1056 0 0
ij 0 0.00189
0 0 0.1280 0
0 0 0.00041 0 0 0.0215
0.1052 0 0
0 0.1289 0
0 0 0.0211
Remarquons que l’influence des déformations élastiques est négligeable sur l’état de
déformation final.
Exemple 3.9
Considérons un point M soumis à l’état de contrainte bidimensionnel suivant :
11 12 0
ij 21 0 0
0 0 0
Les propriétés élastiques sont caractérisées par : le module de Young E = 200 GPa, le
200000
coefficient de Poisson = 0.3 et le module de cisaillement G 76923MPa .
2(1 0.3)
Le comportement du matériau est représenté par un comportement élasto-plastique à
écrouissage cinématique non linéaire d’Armstrong Fredericks en petites transformations, le critère
de plasticité de von Mises et la loi d’écoulement plastique associée. Le Tableau ci-dessous
donne quelques valeurs de la courbe contrainte-déformation de l’essai de traction uniaxiale.
0 0.0009 0.018 0.02045 0.04905 0.09495 0.18235 0.23975
[MPa] 0 180 360 395 450 510 570 590
147
Déterminer les composantes des tenseurs des déformations et des contraintes internes Xij en
considérant les étapes de chargement suivants :
1) L’élément est soumis progressivement à un cisaillement pur jusqu’à ce que les contraintes
12 et 21 soient égales à 300 MPa.
2) On supprime les contraintes 12 et 21.
3) On augmente progressivement la contrainte 11 jusqu’à 500 MPa.
L’essai de traction simple montre que le matériau commence à se plastifier pour une
contrainte y = 360 MPa. Dans le cas d’un écrouissage cinématique, la limite élastique y =
11-X11 = cst = 360 MPa lorsque la déformation plastique évolue. La variation de la
déformation plastique cumulée p 11
p
en fonction de 11 en traction uniaxiale est déterminée
à partir de la règle de partition :
11
p 11
p
11
E
Le Tableau ci-dessous fourni des valeurs de la déformation plastique cumulée p et de la
contrainte interne X11 en fonction de 11 obtenues expérimentalement.
11 [MPa] 0 180 360 395 450 510 570 590
p 0 0 0 0.0185 0.0468 0.0924 0,1795 0,2368
X11 [MPa] 0 0 0 35 90 150 210 230
X11
C
1 e p
Les paramètres C et sont ajustées à partir de la courbe expérimentale de traction en utilisant un
logiciel approprié : C = 2352.3 MPa, = 8.94. On en déduit :
X11 263.09 1 e8.94p
La Figure ci-dessous montre les courbes expérimentales représentant l’évolution de la contrainte 11
et la contrainte interne X11, obtenues expérimentalement et la courbe obtenue à partir de la formule
d’Armstrong et Frédéricks en fonction de la déformation plastique cumulée.
- 1ère Phase
Tenseur des contraintes après plastification :
0 12 X12 0
ij Xij sij Xij 21 X 21 0 0
0 0 0
148
600
500
400
360MPa X Expérimentale
Série1
, X [MPa]
300
X Armstrong et
Série2 Fréderick
200 - Expérimentale
Série3
100
X
0 p
0 0,025 0,05 0,075 0,1 0,125 0,15 0,175 0,2 0,225 0,25
Avant plastification du matériau, les contraintes internes X12 et X21 sont nulles. La
plastification du matériau commence lorsque l’état de contrainte vérifie le critère de von
Mises :
3 3 2
eq sij : sij (12 21
2
) y 360 MPa
2 2
12 = 21, on en déduit:
3
eq s12 X12 : s12 X12 3 12 X12 y 360 MPa
2
X12 = X21.
La condition de cohérence dans le cas d’un écrouissage cinématique s’écrit (voir § [Link]) :
f f
df d ij dXij 0
ij Xij
149
Cette équation contient l’incrément de déformation plastique équivalente qu’on retrouve dans
le comportement des matériaux à écrouissage isotrope. Seules les contraintes 12 et 21 sont
non nulles, d’où :
f f f f
f 12 21 X12 X 21 0
12 21 X12 X 21
f f f f
f
12
12
21
21
X12
C12
p
X12 p
X 21
C 21
p
X 21p 0
.
Les gradients de f par rapport à ont pour valeurs :
f
12
12 12
p
f
C12p X12 2352.312p 8.94X12
X12
Dans cet exemple, les valeurs des incréments 12 et X12 sont prises égales à 20 MPa.
Les incréments des déformations plastiques cumulées p obtenues à la fin de chaque
incrément des contraintes 12 et X12 sont résumées dans le tableau ci-dessous.
Une fois les incréments de déformations plastiques cumulées calculés, on détermine les
incréments de déformation plastique (i) 12
p
en utilisant la règle de normalité associée (voir §
[Link]) puis les composantes des tenseurs des déformations plastiques et totales en utilisant
les expressions suivantes :
150
f 12
(i)
12
p
p(i) 0.866p (i) 12
e
e21 (i)
12
p
(i 1)
p21 (i 1) p21
(i) 2 76923
(i 1)
(i)
p p (i 1) p (i)
12 (i) 12
e
(i) p21
′i′ p
(i) (i)
12
p
(i)
p21 12 (i)
12
e
(i)
e21 (i)
12
p
(i)
p21
(i)
12
1 0,0000 0,0128 207.85 0,00135 0 0,0014
2 0,0147 0,0147 227.85 0,00148 0,0128 0,0142
3 0,0317 0,0173 247.85 0,00161 0,0274 0,0291
4 0,0517 0,0211 267.85 0,00174 0,0448 0,0465
5 0,0760 0,0164 287.85 0,00187 0,0658 0,0677
6 0,0949 300 0,00195 0,0822 0,0841
- 2ème Phase
La suppression des contraintes de cisaillement 12 et 21 entraine une décharge élastique et la
contrainte équivalente devient négative :
- Phase 3
Tenseurs des contraintes après plastification :
2X11 2
3 X12 0 3 11 X11 X12 0
Xij X 21 0 s(ij X ij )
X11 X
X 21 11 11 0
3 3 3
0 X X
0 11 0 0 11 11
3 3 3
3 211 2X11
2
11 X11
2
2 X12
2
eq 3
2 3 3 3 3
3 2
2
11 X11
2
2X11
3X 11
2
12 3
2 3 3 3 3
Lorsque l’état de contrainte se trouve sur la frontière du domaine élastique : X12 = X21 = 92.15
151
MPa, on calcule 11 à partir du critère de von Mises :
3 2
2 2
2 3 3
La plastification commence de façon prématurée pour une contrainte 11 = 322.664 MPa <
360 MPa à cause du chargement considéré dans les phases 1 et 2 (histoire du chargement).
Condition de cohérence :
f f
3 3
f ij Cijp Xijp 0
i, j ij i, j X ij
f f f f
f 11 X11 X12 X 21 0
11 X11 X12 X 21
f f (11 X11 ) f (12 X12 )
f 11 C p X11p C p X12 p
11 X11 y X12 y
f ( 21 X 21 )
C p X 21p 0
X 21 y
12 = 21 , X12 = X21,
2
3 s11 3 3 11
X11
n11
f
f
X11
11
11 X11 2 eq 2 eq 360
f 3s 3 X 21
n12 n 21 12
X12 2 eq 2 360
f
f (11 X11 ) f (X12 )
f 11 p C X11 2 C X12 0
11
X11
y X12 y
f
11
11
p
f (11 X11 ) f (X12 )
C X11 2 C X12
X11 y X12 y
les formules :
(X12 )
(1)
X11 C 11 p (1)
X12 C X12 p
y y
Pour un incrément (i), les formules ci-dessus deviennent :
(i)
X11 = (i-1) X11+ (i-1) X11 (i)
X12 = (i-1) X12 + (i-1) X12
152
avec :
Une fois les incréments de déformations plastiques cumulées calculés, on calcule à chaque
incrément de contrainte 11 les composantes des tenseurs des incréments de déformation
plastique 11p
et 12
p
, des déformations plastiques 11
p
, p22 33
p
et 12
p
, des déformations
élastiques 11
e
en utilisant les expressions suivantes :
f f
11
p
p 12
p
21
p
p (1)
p 0, 0949 (i)
p (i 1) p (i 1) p
11 12
(i)
11
p
(i)
11
p
(i 1) 11
p
(i 1) 11
p
0.084
(1) p
12 (i 1) 12
(i) p
12
p
(i 1) 12
p
(i) 33
(i) p
22
p
2
11 e
11
e
22 33
e
11
e
E
Les résultats des calculs sont résumés dans le Tableau ci-après.
′i′ 11
p
12
p
p 11
p
12
p
p22 11
e
e22 33
e
153
Exemple 3.10
Considérons l’état de contrainte en un point M d’un élément
400 200 0
ij 200 450 80 MPa
0 80 186
800
Contrainte MPa
600
400
200
0
0 0,01 0,02 0,03 0,04 0,05 0,06
Déformation plastique p
Contrainte moyenne :
400 450 186
p 45.33MPa
3
354.67 200 0
sij 200 495.33 80 MPa
0 80 140.67
154
3 3 3 3
eq
2
sij : sij
2 i 1 j1
sij : sij 851.82 MPa 0y 150 MPa
L’état de contrainte considéré engendre des déformations plastiques et admissible (eq < u =
1000 MPa)
Le critère de von Mises peut s’écrire sous la forme suivante :
2 2
sij : sij y
3
L’état de contrainte situé sur la surface de charge est défini comme suit :
2 2 2
y 1502
2 3 3
2sij : sij 2y 1
0.1760937
3
s ij : sij 2 483730.6667
70.44 35.22 0
ij 35.22 79.24 14.09 MPa
0 14.09 32.73
Les composantes du tenseur des déformations élastiques sont définies comme suit :
x y z z x y
1 1 1
ex ey y x z ez
E E E
xy yx xy xz yz
exy xy yx 2 xy exz eyz
2G 2G G G G
0.000422 0.000458 0
ij 0.000458 0.000551 0.000183
e
ij
0 0.000183 0.000177
Posons r = (1-)
Le tenseur des déformations induit par le tenseur des contraintes r ij , dirigé vers l’extérieur
de la surface de charge, est déterminé de manière incrémentale. Pour cet exemple, les
contraintes du champ r ij sont appliquées en 100 incréments.
155
400 200 0 70.44 35.22 0 329.56 164, 78 0
rij 200 450 80 35.22
79.24 14.09 164, 78 370.76 65.91 MPa
0 80 186 0 14.09 32.73 0 65.91 153.24
329.56 164, 78 0 3.3 1.65 0
1
rij 164, 78 370.76 65.91 1.65 3.71 0.66 MPa
100
0 65.91 153.24 0 0.66 1.53
Incrément i = 1 : l’état de contrainte initial est situé sur la surface de charge initiale :
3.3 1.65 0
ij 1.65 3.71 0.66 MPa
0 0.66 1.53
156
On en déduit :
sij : dij
dijp s
2 2 ij
y
3
11
p
12p
13
p
px 2 pxy 2 pxz
p
A noter que : ijp p21 p22 23 pxy 2 pyx py 2 pyz
31
p
32p p
33 2 pzx 2 pzy pz
Dans le cas particulier d’un chargement uniaxial, nous avons :
x 0 0 d x 0 0 1 0 0
2 x
ij 0 0 0 dij 0 0 0 s ij 0 0.5 0
3
0 0 0 0 0 0 0 0 0.5
0 9.28 37.96
mn : d mn
s(1) mn : d mn
3 1 s (1)
h 104435748,8 MPa f 6.72011010
2 2 2 h 2 2
eq eq
3 3
4.197 2.367 0
2.367 5.861 0.947 108 p
2 p
(1)
f s
p
ij
(1)
ij ij : ijp 6.72 108
3
0 0.947 1.665
157
Déformation plastique cumulée :
(1)
p (1)
p 6.7201108
1.974 1.071 0
(1)
1.071 2.578 0.428 105
e
ij
0 0.428 0.828
ij
ije pour i j
2G
L’incrément du tenseur des déformations élastiques ne varie pas au cours de l’évolution des
déformations plastiques.
1.978 1.073 0
(1)
ij 1.073 2.583 0.429 4 105
(1)
e
ij
(1) p
ij
0 0.429 0.83
1.978 2.147 0
(1)
xy ij ij 2.147 2.584 0.859 105
(1) e (1) p
0 0.859 0.83
Incrément i = 2
65.38 36.87 0
s (2)
ij 36.87 91.31 14.75 MPa
0 14.75 25.93
215.46 60.75 0
s (2)
ij : ij 60.75 338.52 9.72 MPa 2
0 9.72 39.73
158
5.438 3.067 0
(2)
f s
p
ij
(2)
ij 3.067 7.595 1.227 108
0 1.227 2.157
(2)
p 8.7078 108 ( 2)
p 6.7201108 8.7078 108 1.5428 107
1.9741 1.0711 0
(2)
1.0711 2.578 0.428 105
e
ij
0 0.428 0.828
1.979 1.074 0
(2)
ij (2)
e
ij
(2)
1.074
e
ij 2.586 0.429 105
0 0.429 0.831
1.979 2.148 0
(2)
xy (2)
e
xy
(2)
2.148 2.586 0.859 105
p
xy
0 0.859 0.83
118.07 66.6 0
66.6 164.67 26.64 10
5
(101) p
ij
0 26.64 46.8
(101)
p 189.0375 105 (101)
p 0.02426244
0.01757 0.01971 0
(101)
xy 0.01971 0.02432 0.00788
0 0.00788 0.00702
159
1000
Contrainte équivalente e
800
Contrainte [MPa]
600
Série2
eq= f(p)
400
= f( p)
Série1
200
0
0 0,005 0,01 0,015 0,02 0,025 0,03
Déformation plastique p
Déformation plastique cumulée p
3 sij
ijp p
2 eq
La déformation plastique cumulée est calculée à partir de l’équation d’Osgood- Ramberg :
1
eq 0.15
p pour eq 851.82 MPa p 0.023 0.024
1500
0.024 étant la valeur de p que nous avons déterminée avec la méthode incrémentale.
0.01436 0.00810 0
3 sij
p
p
0.00810 0.02006 0.00324
2 eq
ij
0 0.00324 0.00570
160
Exemple 4.1
Déterminer le moment plastique de la section représentée sur la Figure ci-dessous.
1
A/2 240 15 180 15 400 8 4750 mm 2
2
240 15 (y plastic 15) 8 4750 mm 2
4750 15 8 240 15
y plastic 158.75 mm
8
yplastic 158.75
ANP
15
G
143.75
y1 132
161
y2 205.22
256.25 G
15
180
Exemple 4.2
Déterminer la charge de ruine Pu de la poutre représentée sur la Figure ci-dessous. Le
moment plastique est constant le long de la poutre.
Mp = Const
- Degré d’hyperstaticité = 1.
- Nombre de rotules plastique nécessaire pour avoir un mécanisme = 2
- Sections critiques : A, C, D
a) Méthode cinématique
- Mécanisme de ruine 1 : rotules plastiques en A et D
d’où :
5M P
Pu
4L
On admet : tan (comportement en petites rotations).
- Mécanisme de ruine 2 : rotules plastiques en A et C
L’équation d’équilibre des travaux virtuels s’écrit :
L
Pu L 2Pu M P M P 2
2
d’où :
3M P
Pu
2L
162
Mécanisme de ruine 3 : rotules plastiques en C et D
L’équation d’équilibre des travaux virtuels s’écrit :
L
2Pu 3M P
2
d’où :
3M P
Pu
L
163
A
On peut aussi procéder autrement. Posons d1d2 = Mp, étant donnée qu’une rotule
plastique se forme au point D. Ainsi, on peut écrire la relation suivante :
Mp Pu L Mp
4
5M p
Pu
4l
La méthode cinématique et la méthode statique fournissent des charges ultimes égales.
En conséquence, la valeur calculée est la charge ultime exacte.
Si on admet qu’au point C se forme une rotule plastique (c1c2 = Mp), le moment au
point D sera supérieur au moment Mp et par suite la condition de résistance n’est pas
satisfaite. Le mécanisme 2 ne peut pas donc être le mécanisme exact.
Exemple 4.3
Déterminer la charge de ruine Pu de la poutre bi-encastrée en A et B représentée sur la Figure
ci-dessous.
P 2P
A B
Mp = const
L C L D L
- Degré d’hyperstaticité n = 2.
- Nombre de rotules plastique nécessaire pour avoir un mécanisme = 3
- Sections critiques : A, B, C, D.
164
- Nombre de mécanismes possibles = 4
a) Méthode Cinématique
- Mécanisme de ruine 1 : Rotules plastiques en A, C et B
P 2 L 2P L M 2 M P 3 M P
u u P
3M p
Pu
2L
Pu 2Pu
L L
L
A B
C D
2 L
2L
3
P L 2P 2 L M P M P 3 M P 2
u u
6M p
Pu
5L
Pu 2Pu
L L L
A B
C D
2
L
2L
2P
u
L M P M P 2 M P
2M p
Pu
L
L Pu L 2Pu L
A C D B
L
165
Pu 2Pu
L L L
A B
C D
L
2
A B
L C L D L
d2
c2
a2 b2 3PuL
c1 d1
Mp Mp
a1 b1
Si le mécanisme de ruine 2 est le mécanisme exact, on doit avoir :
d1d 2 3Pu L M p M p
2M p
Pu
L
On retrouve un résultat identique à celui obtenu avec la méthode cinématique. La condition de
résistance est vérifiée dans toutes les sections de la poutre. On conclut que le mécanisme de
ruine 2 est le mécanisme exact.
Exemple 4.4
Déterminer la charge ultime qu de la poutre encastrée en A et simplement appuyée B
représentée sur la Figure ci-dessous.
166
A Mp = const q B
a) Méthode cinématique
La poutre est une fois hyperstatique, deux rotules plastiques sont donc nécessaires pour
produire un mécanisme de ruine. La première rotule plastique se forme à l’encastrement. En
augmentant la charge q, une deuxième rotule plastique se forme en travée courante et
provoque la ruine de la poutre. Cependant, contrairement aux charges concentrées, on ne
connait pas à priori l’emplacement des rotules plastiques. On suppose donc en premier lieu
que la rotule plastique est située à une distance x de l’encastrement comme il est montré sur la
Figure ci-dessous.
A q Mp = const B
x L-x
A
x/(L-x)
dx B
x
L/(L-x)
C
Le travail virtuel d’une force élémentaire qdx vaut :
dWe q u dx
est le déplacement virtuel d’un point d’abscisse x.
Le travail virtuel extérieur total vaut :
L 1
We q u dx Aire du triangle ABC q u q u x
o 2
L’équation des travaux virtuels s’écrit :
qu L L 2L x
x M p M p M p
2 Lx Lx
2M p 2L x
qu
L Lx x 2
A chaque valeur de x correspond un mécanisme de ruine plausible. Le mécanisme exacte
correspond à la charge ultimé la plus faible. qu est minimum lorsque :
dq u 2M p Lx x 2 L 2x 2L x
0
dx L Lx x
2
Ce qui conduit à résoudre l’équation :
x 2 4xL 2L2 0
167
L’équation à deux racines :
4L 16L2 8L2
x
2
x1 2L L 2 0.558L
x 2 2L L 2 3.41L L
On retient donc la valeur x1= 0.5858L.
La charge q à l’état ultime vaut donc :
2M p 2L 0.5858L 11.66M p
qu
L 0.5858L2 0.5858L 2 L2
b) Vérification par la méthode statique
L’expression du moment fléchissant associée au mécanisme de ruine est déterminée par
superposition de l’expression du moment due à un moment concentré Mp appliqué en A et
celle due à la charge répartie dans la poutre isostatique associée.
qu x 2 Mp qu L
M(x) x Mp
2 L 2
L’effort tranchant a pour expression :
dM(x) Mp qu L
T(x) q u x
dx L 2
L’abscisse du moment fléchissant maximum s’obtient à partir de l’équation :
dM(x) Mp qu L
T(x) q u x 0
dx L 2
Ce qui donne pour le mécanisme de ruine considéré :
L L
x 0.5858L
2 11.66
abscisse de la rotule plastique déterminée à l’aide de la méthode cinématique. Si le
mécanisme de ruine est statiquement admissible, le moment à l’abscisse x = 0.5858L doit être
égal à Mp :
11.66 M p 0.5858L
2
M p 11.66M p L
M(x) (0.5858L) M p M p
2L2 L 2L2
Le moment maximum est obtenu au point C et égal à Mp, la condition de résistance est
vérifiée. Le diagramme des moments de flexion à l’état d’équilibre ultime est représenté sur la
Figure ci-dessous.
168
quLx/2 - qux/2
Mp
Mp Mpx/L
0.5858L
Exemple 4.5
Déterminer pour la poutre continue représentée sur la Figure ci-dessous le moment plastique
minimum nécéssaire pour résister aux charges auxquelles elle est soumise.
225 kN 225 kN
Pour une poutre continue non symmétrique, l’état de ruine se produit selon un mécanisme
partiel de l’une des travées. Deux rotules plastiques sont nécessaires pour produire un
mécanisme de ruine dans une travée de rive simplement appuyée à son extrémité. Les rotules
plastiques se forment sur les appuis intermédiaires et en travée courante. Le nombre de
mécanismes plausibles est donc au moins égal au nombre des travées de la poutre. Pour la
poutre considérée dans cet exemple, trois mécanismes de ruine sont plausibles.
Mécanisme de ruine de la travée AB
On admet que les rotules plastiques se forment au point F à une distance x de l’appui A et sur
l’appui intermédiaire. Au lieu d’utiliser l’équation des travaux virtuels comme dans les
exemples 4.2 et 4.3, on détermine la position des rotules plastiques à partir des équations
d’équilibres des consoles AF et FB (méthode de l’équilibre).
A 60 kN/m 2Mp B
x 12-x
A A B B
Ax
A + B
169
A F
60 kN/m 60 kN/m B Mp
2Mp
x 12-x
- Equilibre de la console AF :
x2
M / A 2M p 60 0
2
M p 15x 2
- Equilibre de la console FB :
60 12 x
2
M/ B 3M p 0
2
M p 10 12 x
2
En égalisant entre les deux équations d’équilibres, on aboutit à une équation du second degré
pour laquelle on retient la solution x = 5.394 m.
Le moment plastique ne doit donc pas être inférieur à la valeur suivante :
Mp 15 5.394 436.43kN.m
2
2 2 12
dM(x)
60x 323.63 0 x 5.394 m
dx
M(5.394m) 30 5.394 323.63 5.394 872.8kN.m 2M p
2
La Figure ci-dessous montre que les conditions de mécanisme et de résistance sont vérifiées .
On peut retrouver la valeur du moment plastique en utilisant la méthode cinématique.
Il est aussi possible de retrouver la valeur du moment plastique à l’aide de la condition
d’équilibre suivante :
M D 2M p M iso 0.45M P 360x 30x 2 0.45M P
360x 30x 2
Mp
2, 45
Pour x 5.394m, Mp 436.32KN.m
Mécanismes de ruine de la travée BC
Les deux mécanismes de ruine plausibles sont :
- Rotules plastiques en B et D ;
- Rotules plastiques en B et E.
170
2Mp
Mp
0.45Mp
A B
5.394 m
Pour déterminer le mécanisme exacte, on utilise la méthode statique ou la méthode d’équibre.
Le diagramme des moments associé à un mécanisme choisi, résulte de la superposition du
diagramme des moments provoqués par les deux charges concentrées et celui des moments
provoqués par un moment concentré appliqué au point B dans la poutre isostatique associée
Dans les deux cas, on obtient le diagramme représenté sur la Figure ci-dessous.
Mp
Mp 675 kN.m
Mp /3
B D E F
Pour que la condition de résistance soit vérifiée, le moment au point E doit être inférieur à Mp,
puisque ME > MD. Le mécanisme de avec des rotules plastique en B et E est donc le
mécanisme plausible.
La position des rotules plastiques étant connue, on peut écrire :
Mp 3 675
M E M p 225 3 Mp 506.25 kN.m
3 4
Le mécanisme exacte est le mécanisme qui conduit au moment plastique le plus grand, soit
Mp = 506. 25 kN.m (> 436.32 kN.m). L’état de ruine de la poutre se produit lorsqu’il se forme
deux rotules plastiques aux points B et E de la travée BF (mécanisme partiel).
Le diagrramme des moments de flexion associé au mécanisme exact est représenté sur la
Figure ci-dessous.
L’équation du moment de flexion dans la travée AB associé au mécanisme exact s’écrit :
60 2 60 12 M p
M(x) x x 30x 317.816x
2
2 2 12
dM(x)
60x 317.81 0 x 5.3m
dx
171
MC 30 5.3 317.816 5.3 841.69kN.m
2
2 506.25
M D 675 337.5 kN.m
3
841.69 < 2Mp = 872.8 kN.m
A 5.3 m C B D E F
Exemple 4.6
Déterminer l’expression de la charge ultime d’une poutre continue à section efficace variable
le long de l’axe longitudinal soumise à une charge uniformément répartie.
Dans la pratique, les poutres continues des batiments et des ponts sont à sections variables.
L’analyse globale s’effectue en considérant un moment plastique sous moments positifs pour
les sections courantes et des moments plastiques sur appuis intermédiares sous moments
négatifs.
a) Travée de rive
Soient :
Mpi : moment de plastification négatif sur appui intermédiaire i ;
Mpt : moment de pastification positif en travée i ;
: déplacement virtuel au droit de la rotule plastique en travée courante.
En se référant aux résultats obtenus dans les exemples précédents, on peut écrire :
0 1
x L1 x
L’équation issue du principe des travaux virtuels s’écrit alors :
q L1
u
M M p1
2
pt x L1 x L1 x
d’où :
q L1 1 1 1
u
M M
2
pt x L1 x p1 L x
1
A chaque valeur de x correspond un mécanisme de ruine plausible. Le mécanisme exacte
correspond à la charge ultime la plus faible pour laquelle :
dq u d 1 1 1 M M M p1
0
pt pt
pt
M M p1
2 2 2
dx dx x L1 x 1
L x x L1 x L1 x
172
Ce qui ramène à résoudre l’équation du second degré suivante :
M x 2 2L1M x L21M 0
p1 pt pt
M M 2pt M M
pt pt p1
x L1
M
p1
A noter que le moment de plastification Mp1 doit être introduit dans ces équations avec le
signe négatif.
dq u M Pt M pi M M pi 1
2 2 0
pt
Li x Li x
2 2
dx x x
M pi
M
pi 1 x 2L M
2
i
pt
M
pi x L M 2
i
pt
M
pi 0
La racine retenue est :
M M
2
M M M M M M
pt pi pt pi pt pi p1 pi 1
x Li
M M
pi pi 1
173
Exemple 4.7
Soit un portique bi-encastré soumis à une charge concentrée appliquée au point C de la poutre et une
charge latérale concentrée au nœud B. Déterminer le multiplicateur de ruine des charges
appliquées.
- Sections critiques : A, B, C, D et E
- Nombre de section critiques n = 5
- Degré d’hyperstaticité p= 3
- Nombre des mécanismes indépendants plausibles k = n-p = 5 - 3 = 2
- Nombre de mécanismes (indépendants + combinés) = 2k -1 = 22 -1= 3.
a) Méthode cinématique
Mécanisme 1 : Mécanisme de poutre
Equation des travaux virtuels :
40 4 60 6 2.25
174
Mécanisme 3 : Mécanisme combiné
Le mécanisme 3 résulte de la combinaison des mécanismes de poutre et de panneau.
Equation des travaux virtuels :
30 3 40 4 60 3 2 1.92
Le mécanisme combiné fournit la charge ultime la plus faible. Il est donc le mécanisme
probablement exact.
b) Vérification par la méthode statique
On construit le diagramme des moments en disposant à l’emplacement des rotules plastiques
une articulation mécanique et un moment concentré égal au moment plastique. La Figure ci-
dessous montre la structure isostatique associée au mécanisme combiné soumise aux forces
extérieures, aux réactions d’encastrement et aux efforts intérieurs au droit des rotules
plastiques.
- Equilibre des efforts dans le poteau ED donne :
Efforts tranchants :
60 60
VE VD 40 kN
3
On en déduit l’effort de compression dans le tronçon CD en considérant l’équilibre des efforts
horizontaux au nœud D et dans le tronçon CD.
L’équilibre des moments dans le tronçon CD donne :
60 60
Vc VD 60 kN
2
L’équilibre des efforts verticaux dans le tronçon ABC donne :
H A 76.8 60 16.8kN
L’équilibre des efforts horizontaux dans le tronçon ABC donne :
VA 57.6 40 17.6 kN
On en déduit le moment au nœud B :
M B 60 17.6 3 7.2 kN.m
Le diagramme des moments de flexion associé au mécanisme de ruine combiné est représenté
sur la Figure ci-dessous.
175
40x1.92= 76.8 kN
60 kN.m
60 kN.m
30 x1.92= 57.6 kN C 40 kN C 60 kN.m
40 kN.m
B 40 kN D
60 kN 60 kN
60 kN
60 kN 60 kN.m
D
40 kN
60 kN.m
A 17.6 kN 40 kN
60 kN.m
76.8 kN E
60 kN
Le diagramme des moments montre que le moment de flexion dans toutes les sections du
portique est inférieur au moment plastique. La condition de résistance est donc vérifiée. Le
mécanisme combiné est donc le mécanisme exact.
Exemple 4.8
Soit un portique bi-encastré en A et E soumis à une charge ponctuelle P = 70 au milieu de la traverse
et une force latérale 30 appliquée en B. Déterminer le multiplicateur de de ruine des charges
appliquées.
176
a) Méthode cinématique
Mécanisme 1 : Mécanisme de poutre
Equation des travaux virtuels :
3M P
70 4 2M 2 M p M P
P 140
Les moments plastiques aux nœuds B et D sont pris égaux à la petite valeur des moments
plastiques des éléments assemblés à ces nœuds, soit Mp dans le cas de cet exemple.
177
Le mécanisme combiné est donc le mécanisme probablement exact
On commence en premier lieu par étudier l’équilibre des moments dans le poteau DE :
Mp
H E 4 2M p H E
2
Une fois la réaction horizontale HE déterminée, on calcule les autres réactions à partir des
équations d’équilibre de la structure complète.
178
L’équilibre des efforts verticaux donne :
63M p 138M p 106.5M p
VA 0 VA
46 23 23
2M p 11
MB 6 Mp M p 0.48M p
23 23
Le diagramme des moments associé au mécanisme de ruine est représenté sur la Figure ci-
dessous. On en déduit que la condition de résistance est vérifiée. Le mécanisme combiné est
donc le mécanisme exact.
Exemple 4.9
Soit un portique bi-encastré en A et E soumis à une charge ponctuelle 2P au milieu de la
traverse et une force latérale P appliquée en B. Déterminer la charge ultime exacte.
a) Méthode cinématique
Mécanisme 1 : Mécanisme de poutre
20Pu 30 20 2 55 Pu 8kN
179
Mécanisme 2 : Mécanisme de panneau
40Pu 30(4) 20(2) Pu 4 kN
180
Le mécanisme de panneau fourni la charge ultime la plus faible, il est donc le mécanisme
probablement exact avec Pu = 4.5 kN.
b) Vérification par la méthode statique
On commence par déterminer les réactions en considérant la structure complète ou un élément
isolé (poteau ED), telle qu’il est représenté sur la Figure ci-dessous.
181
Exemple 4.10
Soit un portique bi-encastré en A et E soumis à une charge uniformément répartie sur la
traverse BC, q= 4 kN/m, et une force latérale en B, H=24 kN. Les sections des montants
et des traverses ont un moment plastique Mp = 80 kN.m.
Déterminer le multiplicateur de ruine des charges appliquées.
4 kN/m
24 kN
B D
5m
12 m
A E
Le travail virtuel externe produit par une charge uniformément répartie sur toute la longueur
d’un élément peut être pris égal au travail externe produit par trois charges concentrées,
comme il est représenté sur la Figure ci-dessous.
Les forces d’extrémités n’ont aucun effet si leur ligne d’action passe par la ligne moyenne
d’un poteau.
a) Méthode statique
Mécanisme de poutre
Equation des travaux virtuels :
4 6 6 80 4 2.22
4 kN/m
24 kN
B D
6
2
6m 6m 5m
A E
182
4 kN/m
24 kN 5 5
5m
x 12 - x
24 kN 5 5
x/12-x
x
12/12-x
5m
183
b) Vérification par la méthode statique
7.224 kN/m VD
43.344 kN
B D HD D 80 kN.m
5m
12 m
A E E HE
80 kN.m HA HE
80 kN.m
VA 80 kN.m VE
VE
- Poteau DE :
M / D HE 5 2 80 0 HE 32 kN
- Portique :
FH HA HE 43.344 0 H A 11.344 kN
122
M / A VE 12 2 80 7.224 2 43.344 5 0 VE 48.071kN
FV 0 VE VA 7.224 12 0 VA 38.617 kN
M B 80 11.344 5 23.28kN.m
- Moment fléchissant de la traverse BD :
23.28 kN.m
7.224 kN/m
x
38.617 kN
7.224x 2
M(x) 23, 28 38.617x
2
On obtient le moment fléchissant maximum lorsque :
dM(x)
38.617 7.224x 0 x 5.345 m
dx
7.224 (5.345) 2
M max M(5,345) 23, 28 38.617(5.345) 79,94 M p 80 kN.m
2
Ce résultat est similaire à celui obtenu avec la méthode cinématique en considérant un
mécanisme de ruine combiné.
Le diagramme des moments de flexion associé au mécanisme de ruine combiné est représenté
sur la Figure ci-dessous. On en déduit que la condition de résistance est vérifiée. Le
mécanisme combiné est donc le mécanisme exact.
184
Exemple 4.11
Déterminer à l’aide la méthode cinématique la charge ultime probablement exacte du portique
à deux travées représenté sur la Figure ci-dessous.
- Degré d’hyperstaticité : p = 6.
- Nombre de rotules plastiques maximum pour transformer le portique en un mécanisme
complet = 6 + 1 = 7.
- Nombre de sections critiques n = 10 (voir Figure ci-dessus).
- Nombre des mécanismes indépendants plausibles k = n-p = 4
- Nombre de mécanismes (indépendants + combinés) = 24 – 1= 15
A Mécanismes indépendants
- Mécanisme (1) : Mécanisme de poutre
2M p
2Pu L M p 4 Pu
L
185
- Mécanisme (2) : Mécanisme de poutre
4M p 2M p
Pu L M p 4 Pu
L L
Ce mécanisme ne peut pas donc être le mécanisme exacte.
- Mécanisme (3) : Mécanisme de panneau
3M p 2M p
Pu 2L M p 6 Pu
L L
Ce mécanisme ne peut pas donc être le mécanisme exacte.
Mécanisme (4) : Mécanisme de nœud
We 0 Wi 3M p
B Mécanismes combinés
Mécanisme (5) = Mécanisme (1) + Mécanisme (2)
186
Mécanisme We Wi Pu
(1) PL 4Mp
2.67M p 2M p
(2) 2PL 4Mp
L L
(5) = (1) + (2) 3PL 8 Mp
Mécanisme Sections
We Wi,ind critiques non -Wi,2 Wi,comb Pu
plastifiées
(1) PuL 4Mp 2 - Mp
(3) 2PuL 6Mp 2 - Mp 2.67M p 2M p
(1) + (3) 3PuL 10Mp 2 - 2Mp L L
(6) 3PuL 2 8Mp
Les rotules plastiques qui forment un mécanisme de ruine combiné ne comprend pas
nécéssairement toutes les rotules plastiques qui forment les mécanismes indépendants à partir
desquels il a été obtenu. Dans le cas du mécanisme (6), la section critique 2 ne comporte pas
de rotule plastique comme dans les cas des mécanismes (l) et (2). Le travail virtuel des forces
intérieures issu du mécanisme (6) est égal à la somme des travaux virtuels issus des
mécanismes (1) et (2) mais sans considération du travail des rotules plastiques qui se forment
dans la section critique 2. La colonne 4 du tableau ci-dessus indique le numéro de la section
critique qui n’est pas considérée dans le mécanisme combiné et dans la colonne suivante le
travail virtuel qui lui est associé.
- Mécanisme (7) = Mécanisme (1) + Mécanisme (4)
187
Mécanisme Sections -Wi,4 Wi,comb Pu
We Wi,ind critiques non
plastifiées
(1) PuL 4Mp 4 - Mp
5M p 2M p
(4) 0 3Mp 4 - Mp
L L
(1) + (4) PuL 7Mp 4 - 2Mp
(7) PuL 4 5Mp
188
- Mécanisme (10) = Mécanisme (3) + Mécanisme (4)
189
Mécanisme Sections -Wi,4 Wi,com Pu
We Wi,ind critiques non
plastifiées
(1) PuL 4Mp 4 - Mp
(2) 2PuL 4Mp 4 0 3M p 2M p
(4) 0 3Mp 4 - Mp L L
(1) + (2) + (4) 3PuL 11Mp 4 - 2Mp
(12) 3PuL 4 9Mp
190
Mécanisme Sections -Wi,5,6 Wi,comb Pu
We Wi,ind critiques non
plastifiées
(2) 2PuL 4Mp 5 - Mp
(3) 2PuL 6Mp 6 - Mp 2.25M p 2M p
L L
(4) 0 3Mp 5,6 - Mp
(2) + (3) + (4) 4PuL 13Mp 4 - 4Mp
(14) 4PuL 9Mp
- Mécanisme (15) = Mécanisme (1) + Mécanisme (2) +Mécanisme (3) +Mécanisme (4)
2M p
Les mécanismes (1) et (8) fournissent la charge ultime probablement exacte Pu .
L
Exemple 4.12
Soit le portique ABDE encastré à la base et soumis à une charge uniformément répartie sur la traverse
BD et une force horizontale appliquée en B. Les montants et la traverse sont des HE 100A. Déterminer
à l’aide de la méthode pas à pas le multiplicateur de ruine des charges appliquées.
191
Le Tableau ci-dessous donne les moments fléchissants dans les sections critiques obtenus
pour les valeurs de = 1 et 1 = 1.446.
Elément Nœud M M1
=1 1 = 1.446
AB A - 4.61 - 6.67
B - 6.98 - 10.1
BD B - 6.98 - 10.1
C 15.449 22.335
x = 2.6 m
D - 22.96 - 33.2
DE D - 22.96 - 33.2
E 19.411 - 28.07
La répartition des moments de flexion, à la fin de la phase élastique (1 =1.446), est montré
sur la Figure ci-dessous. A noter, que contrairement à la convention de signe généralement
adoptée pour construire les diagrammes des moments, le logiciel LUSAS affecte un signe
négatif au moment en travée courante et un signe positif aux encastrements et aux nœuds des
portiques.
192
hyperstatique. Le moment de flexion maximal, obtenu dans le nouveau système statique, se
produit au nœud E et vaut 23.09 kN.m. On calcule l’accroissement du multiplicateur de
charge d2 nécessaire pour entraîner, dans le nouveau système statique, la plastification de la
section E en tenant compte du moment obtenu durant la phase précédente. Le moment de
flexion au point E qu’on obtient avec le multiplicateur de charge d2, pour que cette section
soit plastifiée à la fin de cette phase, doit être égal à : 33.2 – 28.07 = 5.13 kN.m.
L’accroissement du multiplicateur de charge d2 vaut donc :
5.13
d 2 0.22 M E 28.07 23.09*0.22 33.2 kN.m
23.309
A la fin de la deuxième phase, le multiplicateur de charge vaut donc :
2 1.446 0.22 1.666
Le Tableau ci-dessous donne les moments fléchissants dans les sections critiques obtenus
pour les valeurs de = 1 et 2 = 1.666.
Elément Nœud M M2 = M1 + dM2
=1 2 = 1.666
AB A - 21.136 - 11.32
B - 4.454 - 11.08
BD B - 4.454 - 11.08
Section 27.773 28.445
courante de BD
x = 2.67 m
D 0 - 33.2
DE D 0 - 33.2
E 23.309 - 33.2
La répartition des moments de flexion obtenue avec un multiplicateur de charge d2 = 0.22 et
une rotule mécanique en D est représentée sur la Figure ci-dessous.
193
élastique linéaire en considérant un multiplicateur de charge = 1. La structure devient donc
une fois hyperstatique. Le moment fléchissant maximum se produit dans la section C
d’abscisse x = 2.7 m de la traverse BD et vaut 31.667 kN .m. On calcule l’accroissement du
multiplicateur de charge d3 nécessaire pour entraîner, dans le nouveau système statique, la
plastification de la section C en tenant compte du moment obtenu durant la phase précédente.
Le moment de flexion dans la section C qu’on obtient avec le multiplicateur de charge d3,
pour que cette section soit plastifiée à la fin de cette phase, doit être égal à : 33.2 – 28.445 =
4.75 kN.m. L’accroissement du multiplicateur de charge d3 vaut donc :
4.75
d3 0.15 M c 28.445 31.667 *0.15 33.2 kN.m
31.667
A la fin de la troisième phase, le multiplicateur de charge vaut :
3 1.666 0.15 1.806
Le Tableau ci-dessous donne les moments fléchissants dans les sections critiques obtenus
pour les valeurs de = 1 et 3 = 1.806.
La répartition des moments de flexion, obtenue avec un multiplicateur de charge d3 = 0.15
et deux rotules mécaniques en E et D, est représentée sur la Figure ci-dessous.
194
- Phase élasto-plastique : 1.806 < < 4
On dispose une rotule mécanique supplémentaire dans la section C de la traverse BD localisée
à une distance de 2.7 m du nœud B et on effectue une analyse élastique linéaire en considérant
un multiplicateur de charge = 1. La structure devient isostatique. Le moment fléchissant
maximum dans la section C vaut 94.022 kN.m. On calcule l’accroissement du multiplicateur
de charge d4 nécessaire pour entraîner, dans le nouveau système statique, la plastification de
la section A en tenant compte du moment obtenu durant la phase précédente. Le moment de
flexion dans la section A qu’on obtient avec le multiplicateur de charge d4, pour que cette
section soit plastifiée à la fin de cette phase, doit être égal à : 33.2 – 16.82 = 16.38 kN.m.
L’accroissement du multiplicateur de charge d3 vaut donc :
16.38
d 4 0.174 M A 16.82 94.022*0.174 33.2 kN.m
94.022
Après la formation de la rotule plastique au point A, le portique devient un mécanisme et ne
peut supporter aucune charge supplémentaire. Le multiplicateur de charge à l’état de rupture a
donc pour valeur :
4 1.806 0.174 1.98
Le Tableau ci-dessous donne les moments fléchissants dans les sections critiques obtenus
pour les valeurs de = 1 et 4 = 1.98.
La répartition des moments de flexion, obtenue avec un facteur de charge d4 = 0.174 et trois
rotules mécaniques en E, D et C, est représentée sur la Figure ci-dessous.
Elément Nœud M4 M4= M3 + dM4
4 = 1 4 = 1.98
AB A - 94.022 - 33.2
B - 54.09 - 21.262
BD B - 54.09 - 21.292
Section courante 0 33.2
x = 2.7 m
D 0 - 33.2
DE D 0 - 33.2
E 0 -33.2
Le mécanisme de ruine est similaire à celui du mécanisme combiné. Nous allons à présent
utiliser la méthode cinématique pour calculer la charge ultime.
Considérons, à priori, que la rotule plastique à une distance x du nœud B. L’équation
d’équilibre des travaux virtuels s’écrit :
6 x
10 4 20 x 33.2
6 x 6 x
12
33.2 2
6 x
6 12 x
2 x 3.32 1 3.32
6 x 6x
195
12 x 12 x
3.32 3.32 2
6 x 2 x x 4x 12
d
La valeur de est minimale lorsque : 0
dx
d’où :
x 2 24x 60
33.2 2
0
2 x 6 x
x1 = 21.16 m x2 = 2.83 m
On retient la deuxième solution, soit x = 2.83 m.
Le multiplicateur plastique de ruine a donc pour valeur :
12 2.83
3.32 1.988
2 2.83 6 2.83
6.67 kN/m
4 B 4
10 kN D
x
x C
A E
196
Ce résultat est quasiment identique à celui fourni par la méthode pas-à-pas mais avec une
légère différence concernant la position de la rotule plastique dans la traverse BD.
Exemple 5.1
Soit une barre en acier ayant un comportement élasto-plastique avec écrouissage isotrope
linéaire. Les propriétés mécaniques de l’acier sont :
E = 200 000MPa Module d’écrouissage : h =25000 MPa Limite élastique
initiale 0Y 250MPa . Le chargement est appliqué de manière incrémentale. Au pas ′n′, les
valeurs suivantes ont été obtenues : n = 150 MPa, déformation permanente pn = 0.0001.
Déterminer, à l’aide de la méthode du retour radial (voir § [Link]), la contrainte n+1 et la
déformation permanente pn 1 qui résultent d’un acroissement de déformation = 0.002.
Exemple 5.2
Soit une barre en acier ayant un comportement élasto-plastique avec écrouissage cinématique
linéaire. Les propriétés mécaniques de l’acier sont :
E = 200 000MPa Module d’écrouissage : h =25000 MPa
Limite élastique initiale : 0Y 200MPa .
Les paramètres évalués à l’instant tn sont :
La contrainte normale : n = 150 MPa le paramètre d’écrouissage Xn = 50 MPa.
Déterminer, à l’aide de la méthode du retour radial, la contrainte n+1 et la déformation
permanente pn 1 qui résultent d’un acroissement de déformation = - 0.002.
n X n 100 MPa 200 MPa le comportement de l’élément à l’instant tn est élastique.
Incrément de contrainte d’essai essai E 400 MPa
La contrainte d’essai vaut : essai n essai 250MPa
Xessai Xn essai Xn 300MPa
essai 0Y 100 MPa 0 l’état de contrainte correspond à un comportement plastique en
compression.
f essai 100
Incrément de déformation plastique : p 0.000444
E H 200000 25000
197
Contrainte actualisée : n 1 essai signe(essai X essai )E p 161MPa
Variable interne Xn 1 Xn signe(essai Xessai ) H p 48.9MPa
Xn+1 Xn
198
Bibliographie
[1] Coates, R.C. , Coutie, M.G., and Kong, F.K. Structural Analysis, Nelson,1980.
[2] Majid, K.I. Non-Linear Structures, Butterworths, 1972.
[3] Ghali, A. and Neville. A.M, Structural Analysis, Chapman and Hall, 1983.
[4] Bailon, J-P. et Dorlot, J-M. Des matériaux. Presses internationales Polytechnique.
[Link]
[5] Lusas 13. Lusas manual. Finite Element Analysis Ltd. Version 13.5, UK, 2003.
[6] Salençon, J. De l’Élasto-plasticité au calcul à la rupture. Éditions de l’École
polytechnique 2002 - 262 pages.
[Link]
[7] Salençon, J. Viscoélasticité pour le calcul des structures. Éditions de l’École
polytechnique et Presses de l’École nationale des ponts et chaussées 2009 - 158
pages.
[Link]
[8] Neal, B.G. The plastic methods of structural analysis. Chapman and Hall, 1977.
[Link]
[9] Courbon, J. Plasticité appliquée au calcul des structures. Techniques de l’ingénieur C
2050-30.
[Link] › wp-content › uploads.
[10] Walton, G., Hedayat, A. and Rahjoo, M. Relating Plastic potential function to
experimentally obtained dilatancy observations for geomaterials with a confinement-
dependent dilation angle. Int. J. Geomechanics ASCE, 2019, 19(8).
[11] Liingaard, M., Augustesen, A. and Lade, P. Characterization of models for time-
dependent behavior of soils. Journal of Geomechanics, Vol. 4, No. 3, ASCE
September 1, 2004.
[Link]
[12] Maranha, J.R., Maranha das Neves E. The experimental determination of the angle of
dilatancy in soils. 17th International Conference on Soil Mechanics and Geotechnical
Engineering (Alexandria).
[Link]
[13] Nutor R.K and Adomako N.K. Using the Hollomon Model to Predict Strain-
Hardening in Metals. American Journal of Materials Synthesis and Processing. Vol.
2, No. 1, 2017, pp. 1-4.
[Link] › 3216...
[14] Dehard, J. Cours de charpente métallique Analyse des structures (1) en acier. Institut
Hemes Gramme, 2008.
[Link]
[15] Dehard, J., Cours de charpente métallique Vérification des sections (2). Institut
Hemes Gramme, 2008.
[Link]
[16] Chapter 2 Failure/Fracture Criterion.
[Link] › files.
[17] Cours d'introduction à l'analyse numérique. Méthodes de Runge-Kutta.
[Link]
[18] Feng, Z. Mécanique non linéaire. UFR Science et Technologies Université d’Evry
Val d’Essonne.
[Link] › M2_M...
199
[19] Heuzé T. Plasticité des structures. École d’ingénieur. Ecole Centrale de Nantes,
France. 2015.
[Link] › cel-02059317 › document.
[20] François, M. Formalisme général, détail des modèles ECL, EI, ECNL et bases des
schémas d’intégration. Institut de Recherches en Génie Civil et Mécaniques
Université de Nantes.
[21] Kim, N. FEA for Elastoplastic Problems. University of florida.
[Link] › nkim › egm6352 › Chap4
[22] Mohr, D. Integration Algorithms for Rate-independent Plasticity (1D). ETH Zurich.
[Link] › Studium › Downloads › Lecture4.
[23] Kelly PA. Mechanics Lecture Notes : An introduction to Solid Mechanics. University
of Auckland.
[Link]
html.
[24] Gutierrez, ER. Caractérisation mécanique de la matrice de composites type BMC.
Arts et Métiers ParisTech.
[Link] › bitstream.
[25] Baes JM, Bourcart M., Lagoutte C. Essai de traction. Institut Technique d’Ingénieur
d’Industrie, Mai 2004.
[Link] › 19848661-Essai.
[26] Marc F. Plasticité 3D. Institut de Recherche en Génie Civil et Mécanique. Université
de Nantes.
[Link]
[27] Charmet, JC. Mécanique du solide et des matériaux Elasticité-Plasticité-Rupture.
ESPCI Paris - Laboratoire d’Hydrodynamique et Mécanique Physique.
[Link]
[28] Cailletaud, G., Tijani M. Bletry, M., Heripre,E. et Rouabhi, A. Mécanique des matériaux
solides. Support de cours École Nationale Supérieure des Mines de Paris, 2007.
[Link]
[29] Vrouwenvelder, A.C.W.M. The plastic behaviour and the calculation of beams and
frames subjected to bending. Technical University Delft Faculty of Civil Engineering
and Geosciences, 2003.
[Link] .
[30] Baes JM, Bourcart M., Lagoutte C. Essai de traction. Institut Technique d’Ingénieur
d’Industrie, Mai 2004.
[Link] › 19848661-Essai.
[31] Haddag, B. Contribution à la modélisation de la mise en forme des tôles métalliques :
Application au retour élastique et à la localisation. Thèse de doctorat. École Nationale
Supérieure d'Arts et Métiers, 2007.
[Link]
[32] Bouby, C. Adaptation élastoplastique de structures sous chargements variables avec
règle d'écrouissage cinématique non linéaire et non associée. Thèse de Doctorat,
Université des Sciences et Technologies de Lille, 2006.
Archive ouverte HAL tel-00109296 › document [Link] ›.
[33] EN 1992-1-1 : Eurocode 2, Calcul des structures en béton. Partie 1-1 : Règles
générales et règles pour les bâtiments. Octobre 2005.
[34] EN 1992-2 : Eurocode 2, Calcul des structures en béton. Partie 2 : Ponts en béton,
Calcul des dispositions constructives. Mai 2006.
[35] EN 1993-1-1: Eurocode 3. Calcul des structures en acier. Partie 1-1 : Règles
générales et règles pour les bâtiments. Octobre 2005.
200
[36] EN 1994-1-1 : Eurocode 4, Calcul des structures mixtes acier/béton. Partie 1- 1 :
Règles générales et règles pour les bâtiments. Juin 2005.
McKenzie, W.M.C. Examples in structural analysis. Taylor & Francis, 2006.
[Link]
[37] Plastic analysis. University of Marylands.
[Link]
201