Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...

Esto esta así pensado para que luego puedan construirse alguna forma sencilla de generación de
coordenada nodales estructuradas. Lo que se controlaría con alguna letra adecuada. Por ejemplo se
podría leer en la primera columna una “L” y a continuación hacer una interpolación lineal entre dos
nudos previamente leídos, o una “C” e interpolar nudos sobre un arco de círculo, etc. Al respecto
en la versión actual del programa hay algunas opciones de generación sencillas.
Finalmente la rutina va verificando de que nudos se ha leído o generado coordenadas y en caso
de haber coordenadas faltantes (arreglo EXIST(nodes)) se imprime un mensaje de advertencia.
Rutina ELMDAT
Lee de la unidad 5 las conectividades de los elementos y las almacena en la variable CNODS-
(nnode,nelem). Donde NNODE es el número de nudos del elemento y NELEM es la cantidad de elementos
de este tipo. Lee también el tipo de material del que está compuesto en NUMAT(nelem).Esta
rutina se controla según el caracter que aparece en la primera columna:
si ese caracter es una “e” o “E” entiende que se ha terminado con los datos de conectividades
elementales.
si ese caracter es un espacio en blanco entiende que se leerán datos de un elemento, y lee
el número del elemento, el material del que está constituido y sus NNODE nodos en el orden
correspondiente.
si el caracter no es alguno de los anteriores asume que la línea es un comentario y pasa a la
siguiente.
Esto esta así pensado para que luego puedan construirse alguna forma sencilla de generación
de conectividades nodales estructuradas. Lo que se controlaría con alguna letra adecuada.
La rutina verifica que se hayan leído o generado las datos de todos los elementos (arreglo
EXIST(nelem)) y si no fuese así detiene la ejecución. En caso de que las coordenadas de sus nudos
medios falten se generan equidistantes de los vértices. Finalmente para cada nudo asociado a un
elemento se dan de alta los NDOFN parámetros nodales como potenciales grados de libertad del
problema en el arreglo IDNOD(ndofn,nodes)(condición que luego podrá cambiar al introducir las
condiciones esenciales de borde).
Rutina KNOWNV
Esta rutina lee los valores de la variable en aquellos puntos donde es conocida (condiciones
esenciales de borde). Estos valores se van ordenando en forma secuencial y para cada valor conocido,
el grado de libertad asociado se da de baja. El valor conocido (nulo o no) se guarda en el arreglo
KNOWN y en la posición correspondiente del arreglo IDNOD se guarda la posición en el arreglo anterior
pero con signo cambiado. De esta forma en IDNOD tendremos hasta ahora:
1 grados de libertad que no han sido dados de alta
0 grados de libertad dados de alta (efectivos)
-? grados de libertad con valores conocidos
Notar que las condiciones de borde aquí incluidas no son generales. Por ejemplo para problemas
de elasticidad los grados de libertad que se pueden imponer deben coincidir con las direcciones del
sistema global de coordenadas, es decir no es factible en la presente implementación incluir apoyos
inclinados respecto a dicho sistema. Tampoco es posible incluir restricciones conocidas como nodos
“maestros” o restricciones donde un grado de libertad depende de otro(s) en forma lineal.
168