7.2.1.1 Características de un grupo

lectura en IELMNT

READ (LECG,1000) NPAR

FORMAT(20I4)

nota	columnas	variables	definición
	1-4	NPAR(1)	ingresar el número 2
	5-8	NPAR(2)	(NUME) cantidad de elementos 2D en este grupo NUME > 0
1	9-12	NPAR(3)	(INDNL) código correspondiente al tipo de análisis no lineal = 0, por defecto igual a "1" = 1, material no lineal únicamente = 2, formulación Lagrangiana total
2	13-16	NPAR(4)	(IBIOT) = 0 perfectamente drenado = 1 ley de Darcy no generalizada = 2 ley de Darcy generalizada =10 modelación hidromecánica no-drenada = 3 si NDIM =4 et IDYNAM = 2
3	17-20	NPAR(5)	(ITYP2D) código correspondiente al tipo de modelación adoptada para este grupo = 0, simetría de revolución = 1, deformación plana = 2, esfuerzos planos
4	21-24	NPAR(6)	(NDIM) cantidad de "gdl locales" de los nodos de los elementos de este grupo = 2 gdl $u_y$ y $u_z$ (mecánica pura) = 3 gdl $u_y$, $u_z$ y $p$ (Biot simplificado) = 4 gdl $u_y$, $u_z$, $w_y$ y $w_z$ (Biot completo) = 0 por defecto igual a 2
5	25-28	NPAR(7)	(MXNODS) cantidad máxima de nodos empleados para describir un elemento del grupo: $3 \leq \texttt{MXNODS} \leq 9$ por defecto igual a 9
	29-32	NPAR(8)	(IFVOL) código que indica si las fuerzas de gravedad son consideradas. = 0, no hay fuerzas de gravedad = 1, fuerzas de gravedad del sólido y del fluido = 2, fuerzas de gravedad del fluido = 3, fuerzas de gravedad del sólido
	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.
6	37-40	NPAR(10)	(NINT) orden de la integración numérica por defecto igual a 2 $1 \leq \mbox{\texttt{NINT}} \leq 4$
	41-44	NPAR(11)	(IPOUTR) determinada por el programa código que indica la existencia de elementos viga en el modelo
6	45-48	NPAR(12)	(NINTMX) determinada por el programa Empleado para almacenar el número máximo de puntos de integración
	49-52	NPAR(13)	(NTYPEL) tipos de elementos 2D = 1, cuadrángulo de 4 a 8 nodos únicamente = 2, triángulo de 3 à 6 nodos únicamente = 3, cuadrángulo de 9 nodos únicamente = 0, tipo definido para cada elementos, i.e. el grupo puede contener los tres tipos de elementos
7	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 al final de la etapa = 3 impresión de esfuerzos iniciales y al final de la etapa
8	57-60	NPAR(15)	(MODEL) número del modelo seleccionado para el material $1 \leq \texttt{MODEL} \leq 7$ ou $\texttt{MODEL=30}$ = 1, elasto-plástico Drucker-Prager o Mohr-Coulomb = 2, elasto-plástico Hujeux = 3, visco-plástico VPMCC = 4, VEEP-BB Hujeux-Modaressi = 5, visco-elástico Kelvin-Voigt = 6, elasto-plástico con criterio orientado Hujeux modificado = 7, elasto-plástico en esfuerzos principales Hujeux modificado = 8, elasto-plástico Mohr-Coulomb = 10, elasto-termo-plástico de Hujeux = 11, elasto-termo-visco-plástico de Hujeux = 12, visco-plástico de Mohr-Coulomb = 14, Hoeck y Brown = 15, Hormigón armado = 16, modelo para mecánica de rocas (HBO) = 17, elasto-plástico Hujeux no-saturado = 30, modelo de comportamiento definido por el usuario
8	61-64	NPAR(16)	(NUMMAT) cantidad de "juego de propiedades", definidos para el tipo de material seleccionado. $\texttt{NUMMAT} \geq 1$ no hay valor por defecto
8	65-68	NPAR(17)	(NCON) cantidad de constante por juego de propiedades determinado por el programa para todos los modelos salvo MODEL=30
8	69-72	NPAR(18)	(JDETAT) dimensión de la tabla para guardar la información de la historia del material. determinada por el programa para todos los modelos salvo MODEL=30
9	73-76	NPAR(19)	(INCOMP) indicador del tipo de integración en mecánica pura (NDIM=2) o de integración reducida en mono y bi-fásica. = 0, integración estándar (NINT definida por el usuario) = 1, tratamiento de la incompresibilidad en mecánica pura = 2, integración reducida (1 punto por elemento)
	77-80	NPAR(20)	(IEXPL) código define si el grupo será tratado 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

notas

1. Si NPAR(3)=1; los desplazamientos y deformaciones se asumen infinitamnte peque?os.
   Solamente se toman en cuenta las no-linealidades materiales.
   Si NPAR(3)=2; se emplea un formulación Lagrangiana total.
   En este caso se toman en cuenta los grandes desplazamientos y las grandes deformaciones.
   Esta opción sólo está disponible para los modelos MODEL=1 y MODEL=2.
   El estudio térmico (ITHERM>0 en los datos generales) no se puede efectuar en el caso de la formulación Lagrangiana total (INDNL = NPAR(3) = 2).
   Si se ingresa un valor cero, se toma por defecto igual a 1.
2. IBIOT se emplea sólo en dinámica con la formulación Biot simplificada.
   En la formulación de Dary generalizada, se consideran los términos de inercia de la parte sólida para el cálculo de la velocidad de infiltración del fluido:
   $\partial_t w = -k \left( \nabla p + \rho_f g + \rho_f \partial_t^2 u \right)$
   con el valor de gravedad definido en §3. El valor de IBIOT entregado por el usuario es considerado por el programa sólo si IDYNAM= 1 (§1.3).
   El valor de IBIOT se considera también en estática. La tabla siguiente indica los valores posibles de IBIOT de acuerdo al valor de IDYNAM (ver §1.3) y el valor de NDIM (NPAR(6)):

   	IDYNAM
   NDIM	0 (estático)	1 (modelo $u$ o $u-p$)	2 (modelo $u-w$)
   2 (mecánica pura)	0 ó 10	0 ó 10	-
   3 (acoplado $u-p$)	0 ó 1	1 ó 2	-
   4 (acoplado $u-w$)	0 ó 1	-	impuesto en 3

   Un estudio térmico (ITHERM>0) no es posible en el caso no-drenado (IBIOT = 10).
   IBIOT = 10 es posible sólo en mecánica pura (NDIM = 2). Esta opción permite calcular las presiones de poros a partir de las deformaciones de volumen en el caso de un medio perfectamente no drenado (sin flujo interno o externo), si requerir el grado de libertad asociado a la presión de poros.
3. La formulación en esfuerzos planos sólo está disponible para el modelo de Drucker-Prager, es decir, sólo si MODEL=1.
   Algunas observaciones en relación al tratamiento de la incompresibilidad e integración.
4. En estática, los valores posibles son:
   • NDIM=2, en el caso de grupos tratados en mecánica pura (drenado)
   • NDIM=3, en el caso de grupos de elementos tratados con Biot simplificado (no drenado)
   En dinámica, los valores de NDIM deben ser compatible con el valor de IDYNAM definidos en §1.3 y de IBIOT:
   • NDIM=2, en el caso de grupos de elementos tratados en mecánica pura (drenado) ; se requiere que IDYNAM=1.
   • NDIM=3, en el caso de grupos de elementos tratados con Biot simplificado (no drenado); se requiere que IDYNAM = 1 y IBIOT = 1 ó 2
   • NDIM=4, en el caso de grupos de elemento tratado con Biot completo (y no drenado); se requiere que IDYNAM = 2.
5. NPAR(7) fija la cantidad máxima de nodos empleada para describir un elemento del grupo.
   Se requiere un mínimo de 3 nodos y un máximo de 9 nodos.
6. En el caso de elementos rectangulares, un orden de integración de 2 es en general suficiente. Si el elemento está fuertemente deformado, se debe emplear un orden de integración mayor.
   El número de puntos de integración de cada elemento depende del tipo de elementos y del orden de integración:

   Orden de integración	Cuadrángulo	Triángulo
   1	1	1
   2	4	3
   3	9	4
   4	16	7

   Para elementos triangulares, un orden de integración de 3 puede ser requerido en ocasiones. En efecto, si el orden de integración es de 2, los puntos de integración son los centros de los costados; si el orden de integración es de 3, los puntos de integración del triángulo son 4 puntos interiores al triángulo.
   El número máximo de puntos de integración por elemento de cada grupo es determinado automáticamente por el programa en función del orden de integración NINT y del tipo de elemento 2D, NTYPEL, empleado en cada grupo.
   La matriz de masas consistente es siempre calculada empleando un orden de integración de 3.
7. Para un grupo dado, es posible elegir si se escribe o no en el archivo PREFIJO.lst y los esfuerzos del elemento del grupo:
   • si INITPR = 0 no hay impresión
   • si INITPR = 1 impresión de esfuerzo iniciales
   • si INITPR = 2 impresión de esfuerzos obtenidos al final de la etapa en los puntos de integración de los elementos par los cuales IPS vale 1 (§7.2.1.6) y si se trata de una etapa para la cual se guardan los resultados (§1.5, IPRI).
   • si INITPR = 3 impresión idéntica a INITPR = 1 y 2.
   Los esfuerzos efectivos son impresos en los puntos de integración de los elementos.
8. Un único modelo de comportamiento (definido por medio del valor de NPAR(15)) es aceptado para cada grupo de elementos.
   Los modelos de comportamiento siguiente:
   • VPMCC (MODEL=3)
   • VEEPBB (MODEL=4)
   • GEOMAT (MODEL=6)
   • ISOGEO (MODEL=7)
   • TCYCMO (MODEL=10)
   • TVCYCMO (MODEL=11)
   • CYCMOH_NS (MODEL=17)
   son modelos destinados a la investigación. El usuario deberá verificar que corresponden a sus necesidades antes de emplearlos.
   Los valores de NPAR(18) y NPAR(17) característicos del modelos son determinadas por el programa, excepto para el modelo MODEL=30 que debe ser definido por el usuario.
   No entregar ningún valor para NPAR(17) y NPAR(18) si NPAR(15)=30.
   En el caso que se emplee MODEL=30, existen 9 variables pre-establecidas para que el programa pueda guardar los esfuerzos, las deformaciones y las presiones de poros.
   Por lo tanto, sólo se debe ingresar el número de variables de estado específicas al modelo de comportamiento definido por el usuario JDETAT, en este caso se limita a 30: $\texttt{JDETAT}<30$.
9. Si se selecciona esta opción, se calcula por separado la parte isótropa de la parte desviadora del los términos del lado derecho y de la matriz de rigidez. Esta estrategia está justificada para el cálculo de materiales incompresibles (hormigón bituminoso por ejemplo) y evita el bloqueo numérico del modelo.