2.2.1 Lectura de nodos con 4 grados de libertad

SECCIÓN ''LECTURA DE NODOS CON 4 GRADOS DE LIBERTAD''

saltar esta sección si IFORMA <> 0

Repetir esta línea cuantas veces sea necesaria (1)

lectura en INODAL

READ (LECNOD,1000) N,(IDT(i),i=1,3),IDT(7),X1,Y1,Z1,KN,IT,NEQUIV,IPL(N)

FORMAT(5I5,3F10.0,I5,A1,2I5)

nota	columnas	variables	definición
1	1-5	N	número del nodo $1 \leq \mbox{\texttt{N}} \leq \mbox{\texttt{NUMNP}}$
2	6-10	IDT(1)	código del grado de libertad del nodo N
2	11-15	IDT(2)	código del grado de libertad del nodo N
2	16-20	IDT(3)	código del grado de libertad del nodo N
2	21-25	IDT(4 ó 7)	código del grado de libertad del nodo N
3	26-35	X1	coordenada $x$ del nodo N (no se emplea en 2D)
3	36-45	Y1	coordenada $y$ del nodo N
3	46-55	Z1	coordenada $z$ del nodo N
4	56-60	KN	incremento del número de los nodos para la generación automática = 0 no hay generación
5	61	IT	indicador de coordenadas cilíndricas = C coordenadas cilíndrica = blanco coordenadas cartesianas
6	62-66	NEQUIV	número del nodo equivalente
7	67-71	IPL(N)	número de la referencia para los desplazamientos propia al nodo N

notas

1. Los nodos 1 al NUMNP deben ser definidos o generados. En este último caso, se puede emplear una generación automática; (ver nota 4).
   Si el número de los nodos se repite, sólo la última definición en considerada.
   Los datos de los nodos deben obligatoriamente terminar con el nodo NUMNP (la lectura de los nodos se detiene si N=NUMNP)
   Atención: Si se emplea la generación automática,verificar que todos los nodos han sido generados correctamente, con los grados de libertad deseados.
2. Las condiciones de borde de cada nodo N son definidas por IDT(m), m=1,2,3,4 con :
   • IDT(m)=0 un g.l. libre
   • IDT(m)=1 un g.l. bloqueado nulo
   • IDT(m)=-1 un g.l. impuesto no nulo
   • IDT(m)=1 un g.l. inutilizado
   Los valores impuestos a los g.l. no nulos están definidos en §"Solicitaciones nodales externas".
   El sentido de cada g.l. se define en la tabla siguiente, con las condiciones definidas al comienzo de este manual:

   tipo de problema	IDT(1)	IDT(2)	IDT(3)	IDT(4 ou7)
   2D mecánico sin elementos estructurales	-	$u_y$	$u_z$	-
   2D Biot simple sin elementos estructurales	-	$u_y$	$u_z$	$p$
   2D Biot completo sin elementos estructurales	$u_y$	$u_z$	$w_y$	$w_z$
   3D mecánico sin elementos estructurales	$u_x$	$u_y$	$u_z$	-
   3D Biot simple sin elementos estructurales	$u_x$	$u_y$	$u_z$	$p$
   Hidráulico puro en 2D ó 3D	-	-	-	$p$

3. El sistema de coordenadas es una referencia cartesiana.
   Si el problema en estudio es bidimensional, se trabaja en plano $y-z$.
4. Las líneas de nodos no tienen porqué ser ingresadas en secuencia. Sin embargo, todos los nodos de 1 a NUMNP deben ser definidos. Está disponible un método de generación automática:
   Los nodos n1, n1+1*kn1, n1+2*kn1, ...,n1+i*kn1, ...,n2 pueden ser generados mediante 2 líneas de datos:
   n1, ID(1,n1), ..., ID(7,n1), X(n1), ..., kn1
   n2, ID(1,n2), ..., ID(7,n2), X(n2), ..., kn2
   El primer nodo generado es n1+1*kn1,
   El i-ésimo nodo generado es n1+ i*kn1,
   etc ...
   La generación se detiene para i definido por:
   n1+i*kn1 < n2 < n1+(i+1)*kn1
   Notar que n2-n1 puede ser no divisible por kn1.
   Las condiciones de borde de los nodos generados son aquellas de n1.
5. Si bien el programa emplea una referencia cartesiana ($x-y-z$ en 3D, ó $y-z$ en 2D), se puede, a nivel de los datos del mallado, ingreso los nodos en coordenadas cilíndricas. De esta forma se obtiene:
   • en 2D: Y1 representa $r$ y Z1 representa $\theta$.
     En consecuencia $y = r cos \theta$ y $z = r sin \theta$
   • en 3D: X1 representa $r$, Y1 representa $\theta$ y Z1 representa $z$.
     En consecuencia $x = r cos \theta$, $y = r sin \theta$ y $z = z$
6. Es posible imponer valores nodales idénticos g.l. por g.l. entre dos nodos. Estos nodos se dicen equivalentes.
   Si el nodo N1 tiene por nodo equivalente al nodo N2, se debe precisar en la misma línea del nodo N2, que el nodo N2 es equivalente al nodo N1.
   Atención: la generación automática no funciona con los nodos equivalentes.
7. Si el nodo emplea la referencia global, dar 0.
   En caso contrario, entregar el número de la referencia de desplazamientos definida en §2.1.2.
   La referencia de desplazamientos es la referencia en la cual los g.l. (desplazamientos y rotaciones) (archivo PREFIJO.ndrs) del nodo serán calculados.
   Esta opción es aplicable solamente a aquellos nodos conectados a vigas anchas (NPAR(1) =25) o placas anchas (NPAR(1)=26).
   Se debe entregar también las fuerzas y desplazamientos impuestos al nodo en dicho referencial.

FIN DE LA SECCIÓN ''LECTURA DE NODOS CON 4 GRADOS DE LIBERTAD''