3.1.1 Archivo de comandos estático

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

initialization file: homogeneous soil, ECP model, dry case
   64 1001111 2    1         2
    0     
    5       0.2       0.0       1.0 
    1    1    0
99991    0    0 -300   11    .0010    .0100
    1                    
    0
    0
         1         1
        0.        0.     -9.81     1000.        0.
    0    0    0    0    0    0    
    0    0                                                                            * ICON/IDESEQ
         0
   2  30   1   0   3   2   4   1   0   2   0   0   1   2   2   1   0   0   0   0      * Homegenous layer (NPAR(8)=1)
    1     1755.      0.00        0.        1.        1.                               * Toyoura Dr=38%
   444.0E6  222.36E6       0.4    -31.00      -31.    -43.00  -1.800E6    4.0E-3
      0.20    1.0E-4       1.0    5.0E-3     3.E-2       0.8    6.0E-2       3.5
    3.0E-2     1.E-3       0.0       1.0       0.2   1000.E6    500.E6        -2
       100     1.E-4     1.E-4     1.E-4      0.50        0.        0.        0.
                  1.        0.     -1.E6 
    1    1    1    1 
         1         0                                                                  * NCOUCH
         1        0.     1500.
         1         0         0
   2   1   1   0   3   2   4   1   0   2   0   0   1   2   1   1   0   0   0   0      * Bedrock (NPAR(8)=1)
    1     2000.       0.00                
  5305.0E6       0.3       30.     1.E50       30.       0.7     1.E-2     1.E-2      * Elastic
     0.E-0
         2         0                                                                  * NCOUCH
         1      -30.     1800.
         2        0.     1500.
         1         0         0
123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890

Sólo describiremos en detalle algunos aspectos particulares del cálculo estático. El significado de los otros campos puede ser interpretado de acuerdo a la información del manual.

El primer aspecto relevante corresponde al código del tipo de condiciones iniciales. Como se trata del primer cálculo, consideraremos condiciones iniciales nulas: ICON=0 y supondremos que el estado suministrado mediante ICOUCH sólo sirve para inicializar el modelo (IDESEQ=0). Con esta opción, el programa NO guarda las fuerzas internas nodales asociadas al estado de esfuerzos suministrados cuando ensambla el vector de fuerzas nodales. Luego, buscará el equilibrio en las iteraciones según la carga suministrada (sólo la gravedad en este caso).

El parámetro NPAR(8)=1 en la definición de los elementos sólidos correspondientes al suelo indica que se activa la fuerza de gravedad. Como modelo de comportamiento emplearemos el modelo elasto-plástico del ECP (Hujeux) (NPAR(15)=2). Los parámetros suministrados corresponden a una calibración efectuada para una arena de Toyoura, con una densidad relativa del 38%. Para inicializar el modelo emplearemos la opción NCOUCH=1, es decir, definiremos un capa ficticia de densidad 1500 con ``techo'' en $z=0.0$. Se ha seleccionado a propósito un número inferior a la densidad del material de forma de ilustrar el efecto de las iteraciones durante la fase estática.

La definición del estrato rocoso es similar, salvo que se optó por el modelo de Mohr-Coulomb, imponiendo un valor de cohesión muy grande (1.E50) para garantizar el comportamiento elástico. La inicialización mediante NCOUCH es análoga, sólo que empleando dos estratos ficticios.