B.10 Propiedades del material.Modelo elasto-plástico de MOHR-COULOMB

Empleado para elementos volumen 2D y 3D (NPAR(1) = 2 ou 3).
Para este modelo se tiene MODEL=NPAR(15)=8 $\longrightarrow$ Modelo elaso-plástico de Mohr-Coulomb
El modelo 8 es un modelo no lineal
El programa considera NCON=20 y JDETAT=23.

3 líneas de datos

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

nota	columna	variable	definición
línea 1:
	1-10	PROP(1,n)	(XEI) módulo de Young a la presión de referencia: PREF=PROP(10,n)
	11-20	PROP(2,n)	(XNU) coeficiente de Poisson ($\nu$)
	21-30	PROP(3,n)	(XN) exponente de la ley elástica ($n$)
1	31-40	PROP(4,n)	(PHI) ángulo de fricción en plasticidad perfecta ($\phi$) (en grados)
1	41-50	PROP(5,n)	(PSI) ángulo de dilatancia ($\psi$) (en grados)
2	51-60	PROP(6,n)	(XEAUX) módulo de Young para el operador auxiliar $\texttt{XEAUX} \geq \texttt{XEI}$
2	61-70	PROP(7,n)	(XNUUX) coeficiente de Poisson para el operador auxiliar
	71-80	PROP(8,n)	(COHSN) cohesión ($c$)
línea 2:
	1-10	PROP(9,n)	(TRCLMT) tensión límite para el mecanismo de tracción
	11-20	PROP(10,n)	(PREF) presión de referencia por defecto igual a $-1.\times^6$
	21-30	PROP(11,n)	(INCMAX) número máximo de sub-incrementos
	31-40	PROP(12,n)	(AKS(1)) permeabilidad saturada a la porosidad POROS0 en la dirección $x$ (no empleada en 2D)
	41-50	PROP(13,n)	(AKS(2)) permeabilidad saturada a la porosidad POROS0 en la dirección $o_y$
	51-60	PROP(14,n)	(AKS(3)) permeabilidad saturada a la porosidad POROS0 en la dirección $o_z$
	61-70	PROP(15,n)	(AK0) coeficiente de empuje de tierras en reposo ($K_0$)
	71-80	PROP(16,n)	(SATR) saturación residual ($S_r$) en estática y para un problema acoplado en no saturado
línea 3:
	1-10	PROP(17,n)	(ALFAS) coeficiente de la ley $S(p)$ en estática y para un problema acoplado en no saturado
	11-20	PROP(18,n)	(DSDPMX) determinada por el programa, máximo $\frac{d(sat)}{dp}$
	21-30	PROP(19,n)	(FACNAP) coeficiente de una eventual napa freática inicial
	31-40	PROP(20,n)	(ZNAPPE) cota de la napa freática

notas

1. El modelo puede ser asociado (PHI=PSI) o bien no asociado (PHI<>PSI). Los ángulos deben ser ingresados en grados.
2. XEAUX y XNUUX son empleados para calcular la matriz tangente auxiliar, la que es empleada para resolver el sistema no lineal. De forma general, asignar a XEAUX un valor igual o más grande que el valor de XEI. El incremento de dichos valores puede mejorar la convergencia, pero también puede aumentar el tiempo de ejecución.