B.1 Propiedades del material. Modelo de comportamiento de Drucker Prager

Empleado por los elementos 2D (NPAR(1) = 2).
Para este modelo, se tiene MODEL=NPAR(15)=1.

• modelo elasto-plástico de Drucker-Prager (formulación 2D)
• o modelo elasto-plástico de Mohr-Coulomb(formulation 2D)

El modelo 1 es un modelo no lineal.
El programa considera NCON=12 y JDETAT=11.

Dos líneas de datos

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

nota	columna	variable	definición
1era línea :
	1-10	PROP(1,n)	(E) módulo de Young
	11-20	PROP(2,n)	(XNU) coeficiente de Poisson
1	21-30	PROP(3,n)	(PHI) ángulo de fricción
	31-40	PROP(4,n)	(C) cohesión
1	41-50	PROP(5,n)	(PSI) ángulo de dilatancia
2	51-60	PROP(6,n)	(AK0) coeficiente de empuje de tierras en reposo inicial ($K_0$)
3	61-70	PROP(7,n)	permeabilidad saturada para la porosidad POROS0 en la dirección $y$
3	71-80	PROP(8,n)	permeabilidad saturada para la porosidad POROS0 en la dirección $z$
2da línea :
4	1-10	PROP(9,n)	(FACNAP) coeficiente de una eventual napa freática inicial
4	11-20	PROP(10,n)	(ZNAPPE)altura de la napa freática
5	21-30	PROP(11,n)	(IMODEL) indicador del tipo de modelo de comportamiento = 1 modelo de comportamiento de Drucker-Prager = 2 modelo de comportamiento de Mohr-Coulomb por defecto igual a 1
6	31-40	PROP(12,n)	(MLU) cantidad de sub-incrementos por defecto igual a 30

notas

1. El modelo de comportamiento puede ser asociado (PHI=PSI) o no asociado (PHI<>PSI).
   Los ángulos deben ser entregados en grados.
2. $K_0$ sólo se emplea para determinar los esfuerzos iniciales estáticos de acuerdo a la topografía definida por las líneas 7.5.2.3 ("Inicialización de esfuerzos efectivos").
   Los esfuerzos iniciales calculados están definidos a una cota $z$, por el peso de tierras vertical $\sigma_z$ y por $\sigma_x =\sigma_y=K_0 \sigma_z$.
3. Las permeabilidades deben ser definidas en un sistema de unidades compatibles con los otros datos del problema (paso de tiempo por ejemplo). Esta debe tener siempre la dimensión de un unidad de longitud divida por una unidad de tiempo.
4. La activación de una napa freática sólo sirve par inicializar las presiones intersticiales PINT en los puntos de integración del grupo de elementos por medio de:
   PINT = FACNAP * DENW * ABS(GZ) * (ZIPT-ZNAPPE)
   donde
   ZIPT es la cota del punto de integración
   DENW es la densidad del fluido
   GZ es la aceleración de gravedad en la dirección $z$.
   Esta inicialización supone que las incógnitas nodales hayan sido inicializadas previamente por medio de ICON=10 ó 12
5. El modelo de Mohr Coulomb desarrollado aquí (IMODEL=2) es una formulación bidimensional regularizada en la punta del cono, la cual se escribe analíticamente. El modelo de Mohr Coulomb completo 2D y 3D con un límite de tracción y aristas no-regularizadas está disponible como modelo #8 (MODEL=NPAR(15)=8).
6. MLU se emplea por ambos modelos de comportamiento.