3 Construcción del modelo

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 $(0,0)$ y $(1.000, 1 \times 10^{-2} )$, 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 $G$, pero el modelo requiere el módulo de Young E. Para obtener el valor de $E$ basta con evaluar:
  
  $$E = 2 G(1 + \nu)$$
  
  

• 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).