B.6 Propriétés du matériau. Modèle viscoélastique KELVIN-VOIGT

Utilisé par les éléments 2D (NPAR(1) = 2).
On a pour ce modèle MODEL=NPAR(15)=5 $\longrightarrow$ loi viscoélastique Kelvin-Voigt
Le programme prend NCON=9 et JDETAT=8.

2 lignes de données

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

note	colonnes	variable	définition
ligne 1:
	1-10	PROP(1,n)	(YOUNG) module de Young ($E$)
	11-20	PROP(2,n)	(XNU) coefficient de Poisson ($\nu$)
1	21-30	PROP(3,n)	(ETAV) coefficient de l'amortissement volumique ($\zeta_v$)
	31-40	PROP(4,n)	(ETAD) coefficient de l'amortissement déviatoire ($\zeta_d$)
2	41-50	PROP(5,n)	(AK0) coeff. de poussée des terres initial ($K_0$)
1	51-60	PROP(6,n)	(AKS(2)) perméabilité saturée à la porosité POROS0 dans la direction $o_y$
1	61-70	PROP(7,n)	(AKS(3)) perméabilité saturée à la porosité POROS0 dans la direction $o_z$
ligne 2:
3	1-10	PROP(8,n)	(FACNAP) coefficient d'une éventuelle nappe phréatique initiale
3	1-10	PROP(9,n)	(ZNAPPE) altitude de la nappe phréatique

notes

1. Les perméabilités et les viscosités doivent être données dans des unités compatibles avec les unités des autres données (pas de temps par exemple). Les perméabilités ont toujours la dimension d'une unité de longueur divisée par une unité de temps. Les viscosités ont la dimension d'inverse de l'unité de temps.
2. $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$
3. La prise en compte de la présence d'une nappe phréatique n'est utile que pour initialiser les pressions interstitielles PINT aux points d'intégration des éléments du groupe, par:
   PINT = FACNAP * DENW * ABS(GZ) * (ZIPT-ZNAPPE)
   où

   ZIPT est la cote du point d'intégration

   DENW est la masse volumique du fluide

   GZ est la gravité dans la direction $z$.

   Ceci suppose que les inconnues nodales aient été préalablement initialisées par ailleurs (avec ICON = 10 ou 12)