1.2 Líneas de control (definición del problema)

lectura en IGNRAL

READ (LECG,2010)

NUMNP,(IDDL(I),I=1,7),NEGNL,MODEX,IDGEOM,NODLGN,(IPREV(J),J=1,3))

FORMAT(I5,1X,7I1,I2,I5,5X,2I5,2X,3I1)

nota	columnas	variables	definición
1	1-5	NUMNP	cantidad de nodos
2	7-13	IDDL	códigos del tipo de problema
3	14-15	NEGNL	cantidad de grupos de elementos $1 \leq \mbox{\texttt{NEGNL}} \leq 50$
4,7	16-20	MODEX	código de control de ejecución = 0, control de datos = 1, ejecución normal = 2, retoma sin reformar los operadores auxiliares, no disponible a partir de la revisión 7. = 3, retoma reformando todo = 4, retoma de un cálculo dinámico a partir de un estudio estático, reformando todo. = 5, retoma en MODEX = 3 con puesta a cero de los desplazamientos y deformaciones = 6, retoma en MODEX = 4 con puesta a cero de los desplazamientos y deformaciones
5	26-30	IDGEOM	código de la geometría del problema tratado. = 2, bidimensional (plano $y-z$) = 3, tridimensional = 0, por defecto igual a 3
6	31-35	NODLGN	cantidad de nodos (borde superior de elementos) que definen la línea para el cálculo de las resultantes (sólo 2D)
7	38-40	IPREV(1-3)	número de la revisión que ha sido empleada para el cálculo precedente. Se debe entregar 3 cifras; e.g. para la revisión 5.1.0 se entrega IPREV = 510

notes

1. La cantidad total de nodos (NUMNP) controla el número de datos a leer en §3 (``información nodal''). Si NUMNP=0 el programa se detiene.

2. Se debe escoger el código en función del problema a tratar, de acuerdo a la tabla siguiente:

   tipo de problema	IDDL
   2D mecánico sin giros	1001111
   2D mecánico con giros	1000111
   2D Biot simple sin giros	1001110
   2D Biot simple con giros	1000110
   2D Biot completo sin giros	1000011
   2D Biot completo con giros	1000001
   3D mecánico sin giros	0001111
   3D mecánico con giros	0000001
   3D Biot simple sin giros	0001110
   3D Biot simple con giros	0000000
   hidráulico puro 2D o 3D	1111110

   La ejecución directa de un problema de hidráulica pura (IDDL=1111110) no es posible. Sin embargo, se puede tratar un problema de hidráulica pura seleccionado el código IDDL=1001110 bloqueando los grados de libertad mecánicos en la definición nodal. El problema tipo Biot simple corresponde a un modelo acoplado mecánico/hidráulico en estática o en dinámica. Para los problemas dinámicos bifásico, existen tres tipos de formulaciones posibles:

   1. Biot formulación simplificada 1.
      Se trata de una formulación denominada $u-p$, es decir, que las incógnitas son los desplazamientos y rotaciones del esqueleto sólido así como la presión de poros. además, la ley de Darcy es no generalizada.
   2. Biot formulación simplificada 2.
      Es idéntica a la formulación simplificada 1, pero emplea la ley de Darcy generalizada, es decir, toma en cuenta los términos inerciales.
   3. Biot formulación completa.
      Las incógnitas son los desplazamientos nodales absolutos del esqueleto sólidos así como los desplazamientos relativos del fluido intersticial.
      La formulación toma en cuenta tanto la compresibilidad del fluido como la del sólido.

3. Un grupo de elementos queda definido por un conjunto de características comunes (§7.1, vector NPAR). Las características principales son las siguientes:
   • mismo tipo de elementos (volúmenes, vigas, interfaces, etc.)
   • mismo tipo de análisis no lineal
   • mismo modelo reológico
   • mismo orden de integración
   • mismos grados de libertad locales
   • mismas etapas de construcción/excavación
   • mismo tipo de integración en el tiempo
   El número máximo de grupos de elementos es 50. Todos los grupos pueden se no lineales en esta versión. Los grupos no lineales toman en cuenta sólo la no linealidad material.

4. MODEX indica si se efectúa un control, una ejecución o la retoma de un cálculo anterior:
   • MODEX = 0 : Control de datos
     Este control de datos efectúa la generación automática y la verificación de datos, así como la selección de los valores por defecto.
   • MODEX = 1 : Ejecución estándar.
   • MODEX = 2 : Retoma de resultados y de matrices a partir de la última etapa del cálculo precedente. Esta etapa ha sido eliminada a partir de la revisión 7.
   • MODEX = 3 o 4: Retoma a partir de los resultados de una etapa precedente de un cálculo anterior.
     MODEX=3 se emplea para retomas estáticas-estáticas o dinámicas-estáticas
     MODEX=4 se emplea para una retoma estática-dinámica
     Esta retoma se efectúa a partir de los resultados de una etapa cualquiera guardada de un cálculo precedente.
     Las matrices se recalculan cuando se retoma un cálculo precedente.
     Los archivos necesarios para retomar un cálculo son: PREFIJO_SAVE o PREFIJO con el prefijo correspondiente al cálculo anterior.
     Las modificaciones a efectuar el archivo de comandos son las siguientes:
     1. corregir el tiempo inicial (DATDEB) ;
     2. poner ICON = 0 (§condiciones iniciales)
     Observación: una retoma dinámica-estática es imposible.
   • MODEX = 5 o 6 : retoma de un cálculo previo a partir de una etapa calculado poniendo en cero los desplazamientos y las deformaciones.
     Este tipo de retoma de un cálculo se emplea luego de una primera etapa de inicialización de esfuerzos. Cuando se efectúa este tipo de cálculo, se debe poner en cero los desplazamientos y las deformaciones obtenidas durante la etapa de inicialización ficticia.
   • MODEX = 5 efectúa una retoma de tipo MODEX=3 y luego pone en cero los desplazamientos y deformaciones.
   • MODEX = 6 efectúa una retoma de tipo MODEX=4 y luego pone en cero los desplazamientos y deformaciones.
     Por lo tanto, se procede de la misma forma que para MODEX = 3 o 4.
     Observación: Las operaciones a efectuar sobre los archivos antes de retomar un cálculo se describen en el anexo ``retoma de un cálculo'' y pueden ser tratadas mediante el procedimiento de ejecución del programa (§ejecución de GEFDyn)

5. Atención: en el caso de problemas 2D, la ejecución se efectúa en el plano $yz$. IDGEOM debe ser inicializado obligatoriamente como 2 ó 3, para evitar todo error.
   A partir de la revisión 7, si IDGEOM se define como cero, se toma el valor IDGEOM=2.

6. Las resultantes son calculadas sobre una superficie que se define por los puntos que pasan por el extremo superior de los elementos. Estos puntos son definidos por la línea de datos (§2.3).
   En 2D, se trata de las componentes de las fuerzas según los ejes $y$ y $z$ $\left(F_y,F_z\right)$, el momento con respecto al eje $x$ $\left(M_x\right)$ et el flujo a través dicha superficie $\left(\underline{V}.\underline{n}\right)$.
   Las fuerzas son calculadas empleando los esfuerzos efectivos en el caso acoplado.
   En 3D (no disponible), se trata de las componente de las fuerzas según $x$, $y$ y $z$ $\left(F_x,F_y,F_z\right)$, los momentos con respecto a esos ejes $\left(M_x,M_y,M_z\right)$ et el flujo a través dicha superficie $\left(\underline{V}.\underline{n}\right)$.

7. En el caso de retomar un cálculo con MODEX>2, si el cálculo precedente se ha efectuado con una revisión anterior a la 5.3.0., se deberá estipular el número de la revisión empleando la variable IPREV.