2.2 Visualización del modelo desde Matlab y ejecución

Antes de efectuar el cálculo con GEFDyn, es una buena idea verificar la lectura del archivo desde Matlab. Este paso previo permite detectar algunos errores en el archivo de comandos.

En Matlab, está disponible una toolbox gefdyn que permite tanto visualizar como extraer los resultados generados por GEFDyn. Esta toolbox ha sido escrita sobre la base de SDTools por lo que se recomienda leer el manual en línea para explorar otras opciones. Por ahora nos contentaremos con leer el modelo y graficarlo:

>> MODEL=gefread('BXDMC1.in')
>> feplot(MODEL)

El comando gefread permite leer el modelo. Se debe incluir la extensión *.in con el nombre correspondiente. La estructura MODEL contiene una serie de campos:

MODEL = 
       name: 'ENSAYO BIAXIAL DRENADO, MOHR-COULOMB ASOCIADO'
       file: 'BXDMC1.mat'
    Project: [1x1 struct]
       Node: [4x7 double]
        DOF: [5x1 double]
      Stack: {4x3 cell}
        Elt: [2x7 double]

para efectos de visualización las más importantes son MODEL.Node y MODEL.Elt, que contiene las características de los nodos y la conectividad de los elementos (ver manual SDT). El comando feplot(MODEL) abre la interfaz gráfica y traza el modelo:

\epsfig{file=figures/model.eps,clip=,width=.8\textwidth}
Presionando el botón derecho del mouse sobre la figura aparecen opciones para cambiar las vistas, ver el eje de referencia, cursor de elementos o de nodos, etc.

A continuación se procede a ejecutar el modelo. Como el archivo se ha denominado BXDMC1.in, basta emplear el siguiente código de ejecución:

for fichier in BXDMC1 
do 
echo $fichier > prefix.in 
dyn7406c.exe  
rm fort.* 
rm FT00
rm *.l3d1
rm *.prxl
rm *.tbsv
rm *.resu
rm *.conv
done

notar que el nombre del archivo se escribe sin la extensión *.in en el archivo de comandos. La versión actualmente disponible es la dyn7406c. fichier es el nombre de la variable interna del procedimiento de ejecución, por lo que puede ser reemplazado por cualquier otro.

Si el archivo de ejecución se ha denominado gefdyn.bat, bastará con ejecutar:

./gefdyn.bat

desde una consola. El archivo BXDMC1.in debe existir en directorio actual y los archivos BXDMC1_SAVE y/o BXDMC1_TMP de una ejecución previa deben haber sido eliminados.

La ejecución en pantalla tendrá el aspecto siguiente:

         *******************************************************
         *                                                  
         *        Licensed to:                              
         *                                                  
         *       Departamento de Ingenieria Estructural     
         *       y Geotecnica, Pontificia Universidad       
         *       Catolica de Chile                          
         *       Santiago, Chile                            
         *                                                  
         *                 tel: +56(2)354 1273              
         *                 e-mail: esaez@ing.puc.cl         
         *           GEFDYN        RELEASE : 7.4.6          
         *                         VERSION : LINUX          
         *                                                  
         *                                                  
         *                                  Oct 2009        
         *                                                  
         *******************************************************


 
  Enter the Language you prefer ? F (French) or E (English)
  donner le prefixe du fichier de donnees de Gefdyn (*.in)
  
  Dans le cas d'une reprise suivi du prefixe du fichier 
  de donnees du calcul a partir duquel on veut reprendre 
  
  
  
  
  ENTRER LE NOM GENERIQUE (PREFIX) :
 Ketap =         1 Ite =         2 Cpu time = .20000E-01
 Ketap =         2 Ite =         2 Cpu time = .20000E-01
 Ketap =         3 Ite =         2 Cpu time = .20000E-01
 Ketap =         4 Ite =         2 Cpu time = .20000E-01
 Ketap =         5 Ite =         2 Cpu time = .30000E-01
.
.
.
 Ketap =       238 Ite =        77 Cpu time = .41800E+01
 Ketap =       239 Ite =        77 Cpu time = .42100E+01
 Ketap =       240 Ite =        77 Cpu time = .42400E+01

hasta completar los 240 pasos de ejecución. La variable Ketap indica la etapa de cálculo, Ite el número de iteraciones efetcuadas y Cpu time la duración del cálculo.

Esteban Saez 2010-10-28