3.1 Preparación del modelo con GEFDyn

El archivo de comandos para efectuar el cálculo con GEFDyn es el siguiente:

ENSAYO TRIAXIAL MONOTONO,MOHR COULOMB NO ASOCIADO NO DRENADO
    4 1001110 1    1         2                                                   * NUMNP/IDDL/NEGNL/MODEX/IDGEOM/NODLGN
    0    0
  240      .005       0.0      1.20
    1    1    0
99991    0    1 -300    1    .0001    .0001                                      * SREF/NUMREF/IEQUIT/ITEMAX/ITETYP/UTOL/RTOL
    1        0.        0.
    0
    0
         0         1
    1 1 1 1 1 1 1 0              0.        0.        0.    0                     * N/IDT(1:7)/X1/Y1/Z1/KN/IT/NEQUIV/IPL(N)
    2 1 1-1 1 1 1 0              0.        0.        1.    0
    3 1 0 1 1 1 1 0              0.        1.        0.    0
    4 1 0-1 1 1 1 0              0.        1.        1.    0
        0.        0.     -9.81     1000.         0        0.                     * GX0/GY0/GZ0/DENW/COMPRW
    0    1   20    0    2    0 
         1         2         2        00         0
        0.       1.2
       0.0      -1.0
         2        -2         1    10.E-2
         4        -2         1    10.E-2
    0    1
         0
   2   1   1   0   0   3   4   0   0   2   0   0   1   0   8   1   0   0   0   0 * NPAR: NPAR(6)=3 (Hidro-Mecanico)
    1     2700.      0.35        0.        1.        1. 
    10.0E6     0.300       0.0      -30.       -5.    10.0E6     0.300     1.E01
     1.E50     -1.E6       100    1.E-05    1.E-05    1.E-05       1.0        0.
        0.        0.        0.        0.           0
         2         0
         1        4.     0000.
         2       12.     1911.
         0         0         0
    11 4 1                   0    1    1    2    1    3    4
123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890

Similarmente al caso biaxial, se modifica el código del tipo de problema IDDL=1001110, para permitir la variación de la presión de poros.

Luego, se modifica el campo ITETYP a 1 de forma de emplear el esquema de integración adecuado para un problema hidro-mecánico.

En la definición de los nodos, se modifica el campo IDT(7)=0 para que los grados de libertad locales sean compatibles con el tipo de problema.

En la definición de las propiedades del fluido, se ha optado por emplear una compresibilidad COMPRW nula.

Finalmente, en la descripción general del grupo de elementos, se ha asignado al parámetro NPAR(6)=3, de forma de indicar a GEFDyn que existen 2 grados de libertad mecánicos ($u_y$ y $u_z$) y uno hidráulico ($u_w$) por nodo. Más abajo, se ha modificado el ángulo de dilatancia a $\psi=5°$ de acuerdo a los parámetros materiales definidos previamente.

No se requieren modificaciones adicionales para el caso no drenado.