El archivo de comandos para resolver el problema es el siguiente:
RIGID FOOTING * TITRE
189 1001111 1 1 2 * NUMNP/IDDL/NEGNL/MODEX/IDGEOM/NODLGN
0 0 0 * IDYNAM/IMASS/ITHERM
100 1e-2 0. 1.000 * NPASDT/DTINIT/DATDEB/DATFIN
2 10 0 * IPRI/ISAVE/NIVIMP
1111 0 1 550 11 .001 .01 * SREF/NUMREF/IEQUIT/ITEMAX/ITETYP/UTOL/RTOL
1 * IOPE
0 * NCHGDT
0 0 * NPB/IDC
1 1 * IPPGEF/IFORMA
0. 0. -9.81 1000. * GX0/GY0/GZ0/DENW/COMPRW
0 1 2 3 * NCHARG/NCOURB/NPTM/NDIMP
1 1 2 0.0 0 1 * LL/LTYP/NPTC/DATACT/LECHRG/ICHPRI
0. 0. 1.000 -1.e-2 * TEMPC/VALC
1 -2 1 1. * NOD(i)/IDIRN(i)/NCUR(i)/FAC(i) i=1..6
10 -2 1 1.
19 -2 1 1.
0 0 * ICON/IDESEQ
0 * JCOUCH
2 160 1 0 1 2 4 0 0 2 0 0 1 0 1 1 0 * NPAR: NPAR(1)=2 (Solid 2D), NPAR(8)=0 (no hay gravedad), NPAR(15)=1 (Mohr-Coulomb)
1 1763. 0.3 * N/DEN/POROS0
1.333E06 .333 20.0 20.e10 20.0 1. * E/XNU/PHI/C/PSI/AK0
0. 0. 2 * FACNAP/ZNAPPE/IMODEL
0 * NCOUCH
1 * IPPGEF
123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890
Se deja la interpretación detallada del archivo de comandos
al lector. A continuación se indican sólo algunos aspectos importantes:
- En la segunda línea, el parámetro correspondiente al número de nodos NUMNP se puede completar de a acuerdo a la información obtenida con el comando femesh info durante la generación de la malla.
- Se efectuarán NPASDT=100 pasos, con un incremento de DTINIT=1e-2.
- Se empleará la estrategia de Newton modificado reformando cada SREF=1111 etapas, empleando un máximo de ITEMAX=550 iteraciones.
- Se define NCOURB=1 curvas de carga, definidas con un máximo de NPTM=2 puntos y se imponen NDIMP=3 grados de libertad.
La curva de carga LL=1 se define como LTYP=2 y comienza en DATACT=0.0. La curva está definida por dos pares y
, los valores intermedios serán interpolados linealmente.
Se definen los 3 nodos (NOD) afectados, correspondientes a los nodos 1, 10 y 19, código IDIRN=-2: vertical y negativo (desplazamiento impuesto), empleando la curva número NCUR=1 afectada por un factor FAC=1.0.
- Consideraremos condiciones iniciales nulas (ICON=0) y no se itera para equilibrar el estado inicial (IDESEQ=0).
- Se define un único grupo de elementos de tipo volumen 2D (NPAR(1)=2), de PAR(1)=160 elementos. Consideraremos que no hay gravedad (NPAR(8)=0) y emplearemos el modelo de Druker Prager (NPAR(15)=1).
- Ya que no se considera la gravedad, los valores de densidad de masa DEN y porosidad inicial POROS0 son completamente arbitrarios pues no intervienen en el cálculo.
- GEFDyn no posee en su librería de materiales un modelo de comportamiento lineal elástico puro. Sin embargo, podemos ``fabricar'' un comportamiento elástico imponiendo una cohesión artificialmente grande (C=20.e10) para garantizar que nunca se alcance la superficie de fluencia.
En los datos del problema se define el módulo de corte , pero el modelo requiere el módulo de Young E. Para obtener el valor de basta con evaluar:
- Como no hay napa: FACNAP=0. y ZNAPPE=0.. Se empleará además la opción de superficie de fluencia tipo Mohr-Coulomb (IMODEL=2).
- Ya que no hay gravedad y el comportamiento es lineal elástico, no vale la pena inicializar los esfuerzos efectivos (NCOUCH=0).
- Finalmente, los elementos serán leídos en el archivo de geometría *.geom (IPPGEF=1).
Esteban Saez
2010-10-16