1. Problème de référence#

1.1. Géométrie#

_images/Object_1.svg

Coordonnées des points (en mètres):

\(A\)

\(B\)

\(C\)

\(D\)

\(x\)

0.5

\(y\)

0.5

\(z\)

0.5

1.2. Propriétés matériaux pour le modèle LETK#

PA = 0.1

NELAS = 0.

SIGMA_C = 12.

H0_EXT = 1.10292

GAMMA_CJS = 0.8

XAMS = 0.1

ETA = 0.04

A_0 = 0.25

A_E = 0.60

A_PIC = 0.4

S_0 = 0.0005

M_0 = 0.01

M_E = 2.

M_PIC = 6.

M_ULT = 0.61

XI_ULT = 0.365

XI_E = 0.028

XI_PIC = 0.015

MV_MAX = 3.

XIV_MAX = 0.0039

A = 1.5e-12

N = 4.5

SIGMA_P1 = 57.8

MU0_V = 0.1

XI0_V = 0.3

MU1 = 0.1

XI1 = 0.3

1.3. Propriétés matériaux pour le modèle LKR#

PA = .1E6

GAMMA = .85

M_0 = .5

F_P = 0.136047510046

M_1 = 9.69880017363

SIGMA_C = 10.9985715832E6

A_2 = 0.580184800258

Q_I = 100.000648048E6

V_1 = 1.5

V_2 = 1.5

XI_1 = 1.e-2

XI_2 = 1.8e-2

XI_5 = 1.6e-2

A = 1.e-14

N = 3.5

RHO_1 = 1.

RHO_2 = 0.1

RHO_4 = 1.10668567265

R_Q = 1.e-6

R_M = 1.e-6

R_S = 1.e-6

R_X1 = 1.e-6

R_X2 = 1.e-6

R_X5 = 1.e-6

Z = 1.e-6

COUPLAGE_P_VP = 1

A_SIGC = 0.155495602806

B_SIGC = 4.69721443803

1.4. Propriétés matériaux pour le modèle NLH_CSRM#

Les paramètres matériaux sont donnés dans le Système d’unités International (SI).

YoungModulus=7.0e9

PoissonRatio=0.3

IsoCompElasLim=50.0e6

IsoTensElasLim=0.1e6

MCCSlopeCSL=0.5

NLHIndex=1.0

mBigoCritCoef=10.0

aBigoCritCoef=0.75

IncompIndex=15.0

Tau=2.0e2

PerzynaExpo=2.0

NLHModulusP=7.0e9/2.5

NLHModulusV=0.01*7.0e9

1.5. Conditions initiales, conditions aux limites, et chargement#

Phase 1:

On amène l’échantillon à un état homogène: \({\sigma}_{xx}^{0}={\sigma}_{yy}^{0}={\sigma}_{zz}^{0}\) , en imposant la pression de confinement correspondante sur les faces avant, latérale droite et supérieure. Les déplacements sont bloqués sur les faces arrière (\({u}_{x}=0\) ), latérale gauche (\({u}_{y}=0\) ) et inférieure (\({u}_{z}=0\) ).

Phase 2:

On maintient les déplacements bloqués sur les faces arrière (\({u}_{x}=0\) ), latérale gauche (\({u}_{y}=0\) ) et inférieure (\({u}_{z}=0\) ). Sur toutes les faces, la pression de l’eau est nulle.

On applique un déplacement imposé sur la face supérieure de façon à obtenir une déformation \({\epsilon}_{zz}=-6\text{\%}\) (comptée à partir du début de la phase 2). Sur les faces avant et latérale droite, on impose une contrainte de \(5\mathit{MPa}\) .

Pour les modélisations E, F et G, une température croissante est imposée via la commande AFFE_CHAR_MECA entre 0°C et 50°C dans la phase 1 et 50°C et 100°C dans la phase 2.

Pour la modélisationH, la déformation en fin d‘essai est cent fois supérieure à celle obtenue à la fin de la phase de confinement: –0.00028571428571428416 * 100.