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:
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