7.2.2.1 Características de un grupo

lectura en IELMNT
READ (LECG,1000) NPAR
FORMAT(20I4)

nota columnas variables definición
  1-4 NPAR(1) dar el número 3
  5-8 NPAR(2) (NUME) cantidad de elementos en este grupo
NUME > 0
1 9-12 NPAR(3) (INDNL) código que indica el tipo de análisis no lineal
= 0, por defecto igual a "1"
= 1, solamente no linealidad material
7 13-16 NPAR(4) (IBIOT) código que indica el tipo de formulación de Biot empleada :
= 0, perfectamente drenado
= 1, con ley de Darcy no generalizada
= 2, con ley de Darcy generalizada
=10, modelación hidromecánica no drenada
  17-20 NPAR(5) (NTYPEL) tipos de elementos pueden existir en este grupo
= 0, los 3 tipos de elementos
= 1, únicamente hexaedros
= 2, únicamente prismas
= 3, únicamente tetraedros
  21-24 NPAR(6) (NDIM) cantidad de "gl locales" en los nodos de los elementos de este grupo
= 3, mecánico puro $u_x$, $u_y$ y $u_z$
= 4, acoplado $u_x$, $u_y$, $u_z$ y $p$
= 0, por defecto igual 3
1 25-28 NPAR(7) (MXNODS) cantidad máxima de nodos empleados para describir cada elemento por grupo:
8 a 21, para los hexaedros
6 a 15, para los prismas
4 a 10, para los tetraedros
por defecto igual a 21
3 29-32 NPAR(8) (IFVOL) código que indica si se consideran las fuerzas de gravedad.
= 0, no hay fuerzas de gravedad
= 1, fuerzas de gravedad en la fase sólida y líquida
= 2, sólo fuerzas de gravedad en la fase líquida
= 3, sólo fuerzas de gravedad en la fase sólida
  33-36 NPAR(9) (KLAYER) determinada automáticamente por el programa
código que indica si hay una construcción por capas o una excavación que afecte a los elementos de este grupo.
4 37-40 NPAR(10) (NINTRS) orden de integración numérica en las direcciones locales $r$ y $s$ empleada para las cuadraturas gaussianas en los hexaedros
$1 \leq \mbox{\texttt{NINTRS}} \leq 4$
por defecto igual a 2
4 41-44 NPAR(11) (NINTT) orden de integración numérica en la dirección local $t$ empleada en las cuadraturas gaussianas para los hexahedros
$1 \leq \mbox{\texttt{NINTT}} \leq 4$
por defecto igual a 2
  45-48 NPAR(12) (NINTMX) determinada por el programa
Sirve para almacenar la cantidad máxima de puntos de integración
  49-52 NPAR(13) (IPOUTR) determinado por el programa, determina si existen elementos con grado de libertad de giro en el problemas
5 53-56 NPAR(14) (INITPR) código para la impresión de esfuerzos
= 0, no hay impresión para este grupo
= 1, impresión de esfuerzos iniciales
= 2, impresión de esfuerzos en fin de etapa
= 3, impresión de esfuerzos iniciales y finales de la etapa
6 57-60 NPAR(15) (MODEL) número del modelo seleccionado para el grupo
= 1, elasto-plástico Von-Mises
= 2, elasto-plástico Hujeux
= 3, visco-plástico VPMCC
= 4, VEEP-BB Hujeux-Modaressi
= 6, elasto-plástico de criterio orientado GEOMAT
= 7, elasto-plástico en esfuerzos principales Hujeux modificado
= 8, elasto-plástico Mohr-Coulomb
= 10, elasto-termo-plástico de Hujeux
= 11, elast-termo-visco-plástico de Hujeux
= 12, visco-plástico de Mohr-Coulomb
= 14, Hoeck y Brown
= 15, Hormigón armado
= 16, mecánica de rocas (HBO)
= 17, elasto-plástico de Hujeux en no saturado
= 30, modelo definido por el usuario
6 61-64 NPAR(16) (NUMMAT) cantidad de "juegos de parámetros" diferentes para el modelo seleccionado.
$\texttt{NUMMAT} \geq 1$
no hay valor por defecto
6 65-68 NPAR(17) (NCON) cantidad de constantes requeridas por el modelo, determinado por el programa salvo para MODEL=30
6 69-72 NPAR(18) (JDETAT) dimensión de la tabla necesaria para almacenar las variables internas del material.
Calculada por el programa salvo para MODEL=30
8 73-76 NPAR(19) (INTRED) indicador para la integración reducida en mono o bifásico
= 0, integración estándar (NINT definido por el usuario)
= 1, tratamiento de la incompresibilidad (no disponible)
= 2, integración reducida (1 punto por elemento)
  77-80 NPAR(20) (IEXPL) código que permite definir si el grupo se tratará implícitamente o explícitamente en un cálculo dinámico con código de integración IOPE= 10 (§1.7).
= 0, grupo implícito
= 1, grupo explícito

notes

  1. En la versión actual, dar NPAR(3)=0 ó 1.
    Solamente las nolinealidades materiales pueden ser consideradas: los desplazamientos y las deformaciones se asumen pequeñas .
  2. NPAR(7) define el número máximo de nodos empleado para definir cualquier elemento de un grupo. Un mínimo de 4 nodos y un máximo de 21 son empleados para describir elementos 3D.
  3. El cálculo de las fuerzas de gravedad se efectúa con el valor definido en §3.
  4. En el caso de elementos lineales, un orden de integración "2" es suficiente. Si el elemento está fuertemente deformado, un orden de integración más elevado es recomendado. La cantidad de puntos de integración de cada elemento depende del tipo de elemento y del orden de integración.
    tipo de elemento cantidad de puntos de Gauss
    hexaedro NINTRS*NINTT
    prisma 21
    tetraedro 4
  5. Para un grupo dado, se puede elegir entre escribir o no en el archivo PREFIJO.lst los esfuerzos del grupo. Los esfuerzos se imprimen en los puntos de integración de los elementos.
  6. Un solo modelo de comportamiento (definido por NPAR(15)) puede emplearse en un mismo grupo de elementos.
    Los modelos de comportamiento siguientes: son modelos destinados a la investigación. El usuario deberá verificar que se adapten a sus requerimientos.
    Los valores de NPAR(18) y NPAR(17) característicos de un determinado modelo son definidos por el programa a excepción del MODEL=30 que es definido por el usuario.
    No definir ningún valor para NPAR(17) y NPAR(18) si NPAR(15)=30.
    Si se selecciona MODEL=30, existen 13 variables preestablecidas para almacenar los esfuerzos, deformaciones y la presión de poros.
    Se debe definir el número de variables de estado específicas al modelo de comportamiento definido por el usuario. En este caso JDETAT está limitada a 30: $\texttt{JDETAT}<30$.
  7. ver la nota (2) de §7.2.1.1 relativa a IBIOT para los elementos de volumen 2D. El valor de IBIOT se considera incluso en estática. La tabla siguiente define los valores posibles de IBIOT dependiendo de IDYNAM y el valor de NDIM (NPAR(6)):
      IDYNAM
    NDIM 0 (estático) 1 (modelo $u$ ó $u-p$) 2 (modelo $u-w$)
    3 (mecánica pura) 0 ó 10 0 ó 10 -
    4 (acoplado $u-p$) 0 ó 1 1 ó 2 -
    un estudio térmico (ITHERM>0) no es posible en el caso de una modelación no drenada (IBIOT=10).
    IBIOT=10 sólo es posible en mecánica pura (NDIM=2). Esta opción permite calcular las presiones de poros a apritr de las deformaciones de volumen en el caso perfectamente no drenado (sin flujo interno o externo), sin requerir el grado de liberad relativo a la presión de poros.
  8. Ver nota (3) §7.2.1.1

Esteban Saez 2010-12-30