B.9.2 Loi de comportement en 3D

Utilisé par les éléments 3D (NPAR(1) = 3).
On a pour ce modèle MODEL=NPAR(15)=30

Le modèle 30 est un modèle nonlinéaire. Le nombre des propriétés (NCON) et la taille du tableau histoire (JDETAT) doivent être donnés par l'utilisateur au niveau de la carte 7.2.2.1.
6 Contraintes (SIGxx, SIGyy, SIGzz,SIGyz,SIGxz, SIGxy) et 6 déformations (EPSyy, EPSxx, EPSzz, EPSyz, EPSxz, EPSxy) et la pression interstitielle sont déjà prévues dans le tableau d?histoire. Donc, comme valeur de JDETAT, il faut entrer la taille du tableau histoire que l?utilisateur a introduit dans sa loi à part les contraintes, les déformations et la pression interstitielle.

NCON = NPAR(17)

JDETAT = NPAR(18)

La valeur de JDETAT ne doit pas dépasser 30.
Le nombre de lignes de données dépend du nombre des propriétés.

lecture dans EFNL3
	READ (LECG,1010) (PROP(j,n), j=1,NCON)
1010	FORMAT(8F10.0)

Les six premières valeurs du tableau des propriétés sont déjà affectées aux paramètres suivants:

note	colonnes	variable	définition
ligne 1:
1	1-10	PROP(1,n)	(AK0) coef.de poussée des terres initial ($K_0$)
	11-20	PROP(2,n)	(EAUXA) module d'élasticité utilisé pour calculer la matrice de rigidité
	21-30	PROP(3,n)	(XNUAUX) coeff. de Poisson utilisé pour calculer la matrice de rigidité
2	31-40	PROP(4,n)	(AKAUX2(1)) perméabilité saturée à la porosité POROS0 dans la direction $o_x$
2	41-50	PROP(5,n)	(AKAUX2(1)) perméabilité saturée à la porosité POROS0 dans la direction $o_y$
2	51-60	PROP(6,n)	(AKAUX2(1)) perméabilité saturée à la porosité POROS0 dans la direction $o_z$
	61-70	PROP(7,n)	propriété définie par l'utilisateur
	71-80	PROP(8,n)	propriété définie par l'utilisateur
ligne 2:
	$\ldots$	$\ldots$	$\ldots$
	$\ldots$	$\ldots$	$\ldots$

notes

1. $K_0$ ne sert qu'à déterminer les contraintes initiales statiques à l'aide de la topographie définie par les lignes 7.5.2.3 ("Initialisation des contraintes effectives").
   Les contraintes initiales calculées sont définies, à la cote $z$, par le poids des terres supportées pour $\sigma_z$ et par $\sigma_x = \sigma_y = K_0 .\sigma_z$
2. Les perméabilités doivent être données dans des unités compatibles avec les unités des autres données (pas de temps par exemple). Elles ont toujours la dimension d'une unité de longueur divisée par une unité de temps.

Les sous programmes concernés par ce modèle de comportement sont:

• INITMP, appelé par IETA3D qui sert à initialiser les valeurs de AKO et AKAUX.
  Normalement l'utilisateur n'a pas à modifier ce sous-programme sauf si il a des propriétés qui sont ensuite stockées dans le tableau histoire.
• HINT3D, appelé par INIT3D qui sert à initialiser les variables d'histoire. Dans le cas d'existence des variables d'histoire initiales non nulles, il faudra compléter ce sous-programme, sinon ce sous programme se charge d?initialiser les contraintes et les déformations ainsi que la pression interstitielle.
• LOITP3, appelé par SIGEP3 fait appel aux sous-programmes CALDEP, CALAUX, LAWMP et REMISE. Normalement l'utilisateur n'a pas à modifier ce sous-programme.
• CALAUX, appelé par LOITP2 calcule la matrice d'élasticité auxiliaire qui sert à calculer la matrice de rigidité. Les valeurs de la matrice auxiliaire sont calculées avec les propriétés numéro 2 et 3.
• CALDEP, appelé par LOITP2 calcule l?incrément de déformations à une étape donnée.
• LAWTMP, appelé par LOITP3 est le programme qui doit contenir la loi de comportement proprement dit. Ce sous programme a comme argument d?entrée l?état initial (contraintes, déformations, variables d?etat) et l?incrément de déformations. A la sortie, il doit fournir l?état final, c?est à dire les contraintes et les variables d?état finaux. Ce sous programme doit être complété par l?utilisateur.
  Le tableau de contraintes est rangé de manière suivante:
  SIGxx SIGyy SIGzz SIGxy SIGxz SIGxy
  Le tableau de déformations est rangé de manière suivante:
  EPSxx EPSyy EPSzz EPSxy EPSxz EPSyz
• REMISE, appelé par LOITMP, s?occupe de la mise à jour des variables d?histoire. Normalement l?utilisateur n?a pas à le modifier.