Loi de comportement développée Chapitre VII Loi de comportement développée : loi MMC 201 Loi de comportement développée Table des matières I. INTRODUCTION ................................................................................................................................... 203 II. FONDEMENTS DE LA LOI MMC ...................................................................................................... 203 III. DESCRIPTION DES COMPOSANTES DU MODÈLE DÉVELOPPÉ ............................................ 203 III.1. COMPORTEMENT ÉLASTIQUE INCRÉMENTAL..................................................................................... 204 III.2. CRITÈRE DE RUPTURE ....................................................................................................................... 205 III.3. RÈGLE D’ÉCOULEMENT..................................................................................................................... 205 III.4. DÉFINITION DE LA SURFACE DE CHARGE........................................................................................... 206 III.5. LA SURFACE DE CHARGE DE LA LOI MMC ........................................................................................ 207 III.6. CALCUL DE LA DÉFORMATION EN CHARGEMENT PRIMAIRE .............................................................. 211 III.7. CALCUL DU MULTIPLICATEUR ...
I. INTRODUCTION ................................................................................................................................... 203 II. FONDEMENTSDE LA LOI MMC ...................................................................................................... 203 III. DESCRIPTION DES COMPOSANTES DU MODÈLE DÉVELOPPÉ ............................................ 203 III.1. COMPORTEMENT ÉLASTIQUE INCRÉMENTAL..................................................................................... 204 III.2. CRITÈRE DE RUPTURE....................................................................................................................... 205 III.3. RÈGLE D’ÉCOULEMENT..................................................................................................................... 205 III.4. DÉFINITION DE LA SURFACE DE CHARGE........................................................................................... 206 III.5. LA SURFACE DE CHARGE DE LA LOI 207MMC ........................................................................................ III.6. CALCUL DE LA DÉFORMATION EN CHARGEMENT PRIMAIRE.............................................................. 211 III.7. CALCUL DU MULTIPLICATEUR PLASTIQUE........................................................................................ 211 III.8. DÉTERMINATION DE L’ÉTAT DE DÉCHARGEMENT-RECHARGEMENT................................................. 213 III.8.1. en déchargement ............................................................................................................ 214 Module III.8.2. en rechargement............................................................................................................. 215 Module IV. IDENTIFICATION DES PARAMÈTRES DE LA LOI MMC........................................................... 216 IV.1. DÉTERMINATION DES PARAMÈTRES D’ÉLASTICITÉ............................................................................ 216 IV.1.1. K Paramètrec 216 .......................................................................................................................et n. IV.1.2. Paramètre Kd............................................................................................................................... 218 IV.1.3. Paramètreν................................................................................................................................. 219 IV.2. DÉTERMINATION DES PARAMÈTRES DE PLASTICITÉ.......................................................................... 220 IV.2.1. Paramètreϕ. ............................................................................................................................... 220 IV.2.2. Paramètresψ1,ψr....................................................................................................................... 221 IV.2.3. Paramètre SLψ............................................................................................................................ 222 IV.2.4. Paramètre C0............................................................................................................................... 222 V. CONCLUSIONS SUR LA LOI MMC................................................................................................... 222
202
Loi de comportement développée
I.Introduction Afin de représenter le comportement du sol analogique bidimensionnel de Schneebeli utilisé dans notre expérimentation une loi de comportement bidimensionnelle spécifique a été développée. Ce modèle de comportement a été implanté dans un code de calcul en différences finies FLAC. (Cundall [1980]). Vu la nature et le comportement du matériau utilisé, ainsi que la nature des sollicitations dans le creusement des tunnels la loi développée doit être strictement bidimensionnelle. Le comportement en cisaillement écrouissable et fortement dilatant doit correctement simuler les cycles de déchargement et de rechargement. A la rupture il n’y a pas de comportement radoucissant. La loi développée sera appelée loi MMC. Les essais mécaniques (biaxial et oedométrique) réalisés avec des cycles de chargement et de déchargement sur les éprouvettes de matériaux analogiques nous ont permis de déterminer de manière expérimentale les relations caractérisant la rhéologie de ce matériau. II.Fondements de la loi MMC Le modèle que nous développons provient de deux origines différentes : l’élastoplasticité basée sur un seul mécanisme écrouissable avec un critère de rupture défini de type Mohr-Coulomb pour le chargement primaire et l’élasticité incrémentale pour le déchargement-rechargement. Le modèle de comportement développé est caractérisé par: Chargement primaire : •Pas de comportement élastique réversible •La dilatance doit pouvoir être prise en compte dès le début de la sollicitation •Un écrouissage isotrope •Une surface de charge ouverte avec un mécanisme plastique déviatorique •Une règle d’écoulement non-associée Déchargement-réchargement •charge le comportement est élastique non linéaireA l’intérieur de la surface de incrémental de type hypoélastique •Le comportement est distingué entre le déchargement et le rechargement, l’effet de rochet observé dans l’expérience est décrit. La rupture •Le critère de rupture de type Mohr-Coulomb III.Description des composantes du modèle développé En chargement primaire nous avons décomposé la déformation totale en la somme d’une déformation élastique et d’une déformation plastique (Figure 1):
203
Loi de comportement développée
dεij=dεeij+dεijp
Figure 1 : Illustration schématique des déformations élastiques et plastiques pour un essai triaxial Ces déformations sont calculées séparément. La déformation élastique est obtenue en utilisant la loi de Hooke. La déformation plastique, par contre, nécessite l’existence d’une surface de charge, d’une règle d’écoulement et d’une loi d’écrouissage durcissant. La formulation que nous proposons correspond à la logique de différences finies, où d’élément finis, dont la variable de base est le déplacement. Les déformations plastiques sont obtenues par dérivation et les contraintes sont calculées en utilisant la loi rhéologique. Afin de présenter les séquences logiques de développement nous décrivons les composantes dans l’ordre suivant : Comportement élastique, critère de rupture, règle d’écoulement, surface de charge, comportement en déchargement-rechargement. III.1.Comportement élastique incrémental L’élasticité est fondée sur l’élasticité de Hooke, avec un module de Young dépendant de la contrainte principale mineureσ2. Les nouvelles contraintes élastiques sont calculées à partir des déformations dans les axes globaux OX et OY : σxx= α1εxx+ α2εyy σyy= α1εyy+ α2εxxτxy=2Gεxy Avec : α1= E0le module oedométrique α2= K0E0, avec K0le coefficient des terres au repos Dans des conditions bidimensionnelles, on obtient les valeurs suivantes en élasticité isotrope : = E01E2K0=−ν La valeur deνest constante ( 0<ν<1,ν=1 correspond à un matériau incompressible).
204
Loi de comportement développée
Le module tangent initial E est variable. Il est calculé à partir de la formule proposée par Janbu [1963] : ⎛ σ ⎞n K2 Ei=⎜⎝Patm⎟⎠avec : K coefficient adimensionnel Patmla pression atmosphérique de référence Les 4 paramètres définissent l’élasticité isotrope : Kc, Kd,ν, n. Le paramètre Kdsera défini plus loin. III.2.Critère de rupture Dans la loi MMC le critère de rupture est celui de Mohr-Coulomb, il est parfaitement adapté à notre modèle 2D car il ne dépend que de deux contraintes principales. Il est exprimé sous la forme suivante: F(σij)= σ1(1−sinϕ)− σ2(1+sinϕ)≤0Avec : σ1etσ2respectivement les contraintes principales (σ1>σ2>0). ϕl’angle de frottement du matériau (paramètre) Dans le plan des axes principaux, la surface de rupture est présentée par deux droites qui sont asymétriques par rapport à l’axe hydrostatique. Le critère de rupture passe par l’origine des axes principaux ce qui représente l’absence de cohésion dans notre matériau (Figure 2). 100 90 Critère de ru 80ermocranoisspture70σ1/σ2=tg2(π/4+ϕ/2) 60 50 40 30 20σ1/σ2=tg2(π/4-ϕ/2)10Critère de rupture ar extension 0 0 10 20 30 40 50 σ2 (kPa) Figure 2 : Critère de rupture dans l’espace des axes principauxσ1etσ2.
III.3.Règle d’écoulement En utilisant la règle d’écoulement plastique, l’incrément du tenseur des déformations plastiques (dεip) est colinéaire aux dérivées partielles du potentiel plastique, il va donc être déterminé de la manière suivante :
205
Loi de comportement développée
= dεipdλp∂σGii=1,2∂ où G est le potentiel plastique dλpla longueur de la déformation plastiqueest un scalaire représentant à Le potentiel plastique G pour un matériau granulaire est en générale différent de la surface de charge et les règles d’écoulement non-associées sont appliquées (Figure 3). Le potentiel plastique utilisé s’écrit sous la forme suivante: G(σij)= σ1(1−sinψ)− σ2(1+sinψ)+Cteoùψl’angle de dilatance variable (0<ψ<ϕ) d’où les incréments de déformations plastiques: dεp1=dλs(1−sinψ)dεp2=dλs(1+sinψ)dε1petdε2pétant les incréments de déformation plastique principale. ψ=ψr=Cteà la rupture
Figure 3 : Direction de l’incrément de déformation plastique Les 2 paramètres de la loi définissant la rupture sont :ϕ,ψr. III.4.Définition de la surface de charge La construction de la surface de charge se fait généralement à partir de résultats expérimentaux en considérant certaines hypothèses simplificatrices. La plupart des modèles élastoplastiques développés pour les géomatériaux sont généralement basés sur deux notions théoriques : le travail plastique et la surface de charge. Tous les modèles de la famille Cam Clay sont basés sur la première notion. Le modèle MMC se base sur la deuxième notion. L’expression mathématique de la surface de charge peut être obtenue à partir de la surface de rupture. Au cours de la sollicitation, la surface de charge évolue de façon isotrope vers la surface de rupture jusqu’à ce que les deux surfaces soient confondues. Pour déterminer la surface de charge nous avons utilisé la forme de la surface de rupture définie par le critère de Mohr-Coulomb. L’idée est basée sur le fait qu’en utilisant la surface
206
Loi de comportement développée
de rupture, par exemple Fr(σij)=r, où r est une constante. Une fonction de charge Fc, peut être déterminée par : Fc(σi,α)=Fr(σi)−H(α)=0 Oùαreprésente la variable d’écrouissage et H une fonction qui vaut à la rupture H(α)=r. On augmente la taille de la surface de charge vers la surface de rupture selon une loi d’écrouissage isotrope. III.5.La surface de charge de la loi MMC Afin de déterminer notre surface de charge, nous avons adopté la démarche suivante : La variable SL (Niveau de contrainte) est le rapport du rayon du cercle de Mohr actuel et du rayon du cercle de Mohr à la rupture calculé avec la même contrainte principale mineureσ2. (Figure 4). 30 Cercle de Mohr actuel 25Cercle de Mohr à la rupture 20Critère MC SL=Ract/Rrupφ 15 10 5Rrup=(σ1−σ2)/2 Ract 0 0 10σ 30 40 50 602 = 20 σ(kPa) Figure 4 : Représentation de la variable SL dans le planτ-σSL(σ1− σ2)(1−sin )0≤SL≤1 = 2σ2sinϕ Avec SL niveau de contrainte actuel. SLmaxgrande valeur de SL obtenue au cours du chargement.est la plus A chaque incrément de calcul n la valeur de SLmaxde l’étape précédente est réajustée en fonction de la nouvelle valeur deσ2: −σ SLmaxnSLn1amx⎜⎛2nn−1⎞⎟m =⎜⎝σ2⎟ ⎠ m est une constante dépendant de la nature de matériau (dans la loi m=0.25 pour les petits rouleaux) La variableSL, définie précédemment, avec 0<SL<1, ne permet pas de faire la distinction entre l’extension et la compression et une forte valeur deSL obtenue en extension empêcherait le retour au chargement primaire en compression.
207
Loi de comportement développée
Critère de rupture pa compression
1/2=tg2(π/4+ϕ/2)
1/2=tg2(π/4-ϕ/2) Critère de rupture pa extensionFigure 5 : Chemin de contrainte et la courbe contrainte-déformation pour un essai triaxial Afin de distinguer le comportement en compression du comportement en extension (Figure 5) nous avons donc créé la variableSLV: (i j1()−sinϕ) SL= v2σ ϕmin(σi,j) sin oùσietσjsont les contraintes principales non classées SLvest comprise entre 1 et +1 : •SLV= -1 correspond à la rupture par extension avec une forte baisse du premier invariant de contrainte I1=σ1+σ2. •SLVcorrespond à la rupture par compression avec augmentation de la variable I= +1 1Le calcul du signe deSLVsera basé sur l’inversion brutale de l’inclinaison des axes principaux et contrôler par l’évolution de premier invariant. La Figure 6 montre la variation de la variable SL dans un essai triaxial. La variable SL ne distingue pas le comportement en extension. Déchargement Rechargement 1 0,8 0,6 0,4 0,2 0 0 10 20 30 40 50 60 70 80 90 100 -0,2 -0,4 -0,6 -0,8 Pas de calcul (*1E+03) -1 Figure 6 : Variation de SL correspondant à l’essai biaxial
208
Loi de comportement développée La variation de la variable SLvun essai biaxial est présentée sur la Figure 7.dans 1 0,8 0,6 0,4Déchargement Rechargement 0,2 0 0 2 0 10 20 30 40 50 60 70 80 90 100 , -0,4 --0,6 -0,8 Pas de calcul (*1E+03) -1 Figure 7 : Variation de SLvcorrespondant à l’essai biaxial 50
40 30
Rechargement
20 10Déchargement Pas de calcul (*1E+03) 0 0 10 20 30 40 50 60 70 80 90 100 Figure 8 : Variation de I1correspondant à l’essai biaxial Le sens de variation de SLVsera contrôlé par la variation de I1(Figure 8). Si I1diminue,SLVpeut devenir négatif et diminue même siSL (car augmente|SLv|=SL). De même si I1augmente,SLVaugmente. La surface de charge s’écrit sous la forme : n−1 m F=SLvmn−xa1σ⎜⎛σ⎜⎝n2⎟⎞⎟⎠2 Avec F=±1 a la rupture L’évolution de la surface de charge tenant compte de la compression et de l’extension présentée sur les Figure 9 et Figure 10. Dans les deux représentations les surfaces de charge sont des courbes et ne sont pas des droites.
209
Loi de comportement développée
80Critère MC en compression 60 Chemin de contrainte pour SLv=1 40σ2=30kPA SLv=0,5 20 SLv=0 0 0 20 40 60 80 100 -20SLv=-0,5 -40SLv=-1 -60Critère MC en extension -80 p=(σi+σj)/2 Figure 9 : L’évolution de la surface de charge en plan p-q avecσ2=Cte 100 Critère de rupture SLv=1exA par compression SLv=0,5 hydrostatique 80
60
SLv=-0,5
SLv=-1
40 20 Critère de rupture par extension 0 0 20 40 60 80 100 σj(kPa) Figure 10 : L’évolution de la surface de charge en planσi-σj. L’évolution de la surface de charge utilisée dans notre loi de comportement est en accord avec les résultats expérimentaux effectués par Tatsuoka et Ishihara [1974] pour les essais triaxiaux avec trois densités différentes. En analysant les résultats présentés sur la Figure 11 on constate que la surface de charge est convexe et contient l’origine. Dans le cas du sable lâche la surface de charge à tendance à se fermer sur l’axe hydrostatique, par contre ce phénomène est moins significatif pour le sable dense.
210
Loi de comportement développée
Figure 11 : Observation de la surface de charge pour le sable lâche (a), sable moyennement dense (b) et très dense (c). (D’après Tatsuoka et Ishihara [1974]) III.6.Calcul de la déformation en chargement primaire Si la condition d’irréversibilité est satisfaite:σijdεpji≥0 matériau est en chargement le primaire. La direction de déformation plastique est normale au potentiel plastique et la variation de déformation plastique est donnée par: C=C0+SL 1−C0)où C est un paramètre d’écrouissage. C0est le paramètre de la loi (0<C0<1). Il exprime du pourcentage de déformations plastiques dès le départ de la sollicitation, il est fixé pour l’état SL=0 La proportion de la déformation plastique C évolue avec SL. A la rupture C=1 quand SL=1. ce qui permet de trouver la longueur de la déformation plastique et de calculer dεp: dεjip=dλp(C)mG∂σij Où : m est fixé en fonction de la nature de sol. Ce qui définit complètement l’état de déformation à partir de l’état de contrainte. L’angle de dilatance diminue au cours du chargement deψ1au début du chargement àψrà la rupture. Le passage entre les deux angles est contrôlé par la valeur de SLψ de la (paramètre loi) compris entre 0 et 1, donnant la valeur de SL à partir de laquelle l’angle de dilatance passera progressivement deψ1deψr. Les trois paramètres du modèle MMC définissant la plasticité avant rupture : C0, SLψ,ψ1. III.7.Calcul du multiplicateur plastique Le multiplicateur plastique s’obtient à partir de la condition de consistance: Fσ +dσ,α +dα)=0