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

More documents

Recommendations

Info

Descripción global del programa A continuación se desarrolla una descripción del significado de cada rutina del programa según el orden en que van apareciendo. El flujo global del programa se controla desde el programa principal donde se definen la variables indicadas en la sección anterior. Rutina OPENFI Esta rutina aglutina los comandos necesarios para abrir los archivos a ser usados. Los archivos son de tipo ASCII, es decir que pueden verse y editarse con cualquier editor. Lee interactivamente los nombres de los mismos, verifica su sintaxis y les asigna las características adecuadas. Los archivos que abre son los siguientes unit 3 Archivo de salida, el nombre de este archivo es ingresado por el usuario. Allí van a parar el eco de la entrada de datos, valores generados y resultados de las variables nodales. Posteriormente también se escriben allí las reacciones nodales, el valor del flujo en los puntos de integración y los valores suavizados en los nodos. unit 4 Archivo de salida donde se escriben algunos mensajes de advertencia o para escribir valores auxiliares en la fase de depuración del programa. unit 5 Archivo de datos (ingresado por el usuario) primero verifica su existencia y concatena (agrupa) los archivos en que pueden estar separados los datos en un único archivo GAMMA.DAT, a estos fines llama a las rutinas GENFIL, RANDWR unit 7-9 Archivos ASCII de salida orientados a ser usados como interfaces con programas de visualización (Tecplot, GiD). unit NN Archivo auxiliar Rutina MATPRO Lee de la unidad 5 las características de los materiales y las almacena en la variable PROPS(nprop,nmaty Esta rutina (como otras) se controla según el caracter que aparece en la primera columna: si ese carácter es una “e” o “E” entiende que se ha terminado con los datos de materiales. si ese carácter es una “m” o “M” entiende que se empezarán a leer datos de un nuevo material, y lee el número del material correspondiente. cualquier otro caracter hace que lea el valor de una característica (en forma consecutiva) del material a partir de la columna 31, sirviendo las 30 primeras columnas como un espacio para comentario. Una vez terminada la lectura la rutina verifica que se hayan leído valores para todos los materiales (arreglo EXIST(nmaty)) y si no imprime una advertencia. Rutina COORDG Lee de la unidad 5 las coordenadas de los nudos y las almacena en la variable COORD(dimen,nodes). 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 coordenadas nodales. si ese caracter es un espacio en blanco entiende que se leerán datos de un nodo, y lee el número del nodo y su DIMEN coordenadas correspondientes. si el caracter no es alguno de los anteriores asume que la línea es un comentario y pasa a la siguiente. 167
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
Page 1 and 2:
Capítulo 1 Métodos de residuos po
Page 3 and 4:
( ∂ k ∂u ) + ∂ ( k ∂u ) + F
Page 5 and 6:
δ (x − x l ) = 0 para x ≠ x l
Page 7 and 8:
p = F (1.29) en la que k y F no dep
Page 9 and 10:
La solución de los sistemas conduc
Page 11 and 12:
Una viga de longitud 1 está simple
Page 13 and 14:
Notar que K es no simétrica aun cu
Page 15 and 16:
integrando por partes el primer té
Page 17 and 18:
La sustitución de la aproximación
Page 19 and 20:
Figura 2 Ejercicio 21. a) Problema
Page 21 and 22:
Capítulo 2 El método de elementos
Page 23 and 24:
Figura 1 Aproximación de una funci
Page 25 and 26:
Figura 2 Comportamiento de tres fun
Page 27 and 28:
2.4.1. Propiedades de la matriz K y
Page 29 and 30:
Figura 5 Descripción local y globa
Page 31 and 32:
2.4.3.2. Cálculo del vector de car
Page 33 and 34:
de funciones. Una medida universalm
Page 35:
Figura 6 Curvas del error utilizand
Page 38 and 39:
eferiremos como ” ley de conserva
Page 40 and 41:
ambos miembros haciendo tender el p
Page 42 and 43:
1. Condiciones de borde esenciales,
Page 44 and 45:
2. Las condiciones de borde no pued
Page 46 and 47:
de tal modo que todo punto x en el
Page 48 and 49:
N 1 (ξ) = (ξ − ξ 2) (ξ 1 −
Page 50 and 51:
modo de ubicar un nodo en todo punt
Page 52 and 53:
Figura 5 Malla de N nudos y N-1 ele
Page 54 and 55:
3. Condiciones naturales de Neumann
Page 56 and 57:
(c) Considere las condiciones (ii)
Page 58 and 59:
los polinomios de interpolación in
Page 60 and 61:
La derivada segunda de u respecto a
Page 62 and 63:
4.3.3. Formulación a partir del Pr
Page 64 and 65:
¯K = ∫ L 0 BT D B d¯x 1 = ¯K =
Page 66 and 67:
Las cargas normales al eje de la vi
Page 68 and 69:
2. θ 1 (ξ) = ∑ 2 I=1 N I (ξ)
Page 70 and 71:
4.5.7. Cambio de base La expresión
Page 72 and 73:
Figura 5 Viga con deformación de c
Page 74 and 75:
4.6.1. Matriz de rigidez de una vig
Page 76 and 77:
Donde las funciones de forma son: N
Page 78 and 79:
(una para cada W J ) en función de
Page 80 and 81:
Ejercicio 1-Sea el problema de conv
Page 82 and 83:
o directamente las coordenadas noda
Page 84 and 85:
donde r es el residuo que se quiere
Page 86:
⎡ K 3−4 = ⎢ ⎣ K 4−5 = ⎡
Page 89 and 90:
Figura 1 Conducción del calor en 2
Page 91 and 92:
5.3. Forma variacional del problema
Page 93 and 94:
5.5. Flujo en un medio poroso El fl
Page 95 and 96:
ortotropía coinciden con las direc
Page 97 and 98:
5.6.5. Forma débil de la ecuación
Page 99 and 100:
Como φ es conocida sobre el contor
Page 101 and 102:
La formulación débil que resulta
Page 103 and 104:
5.9.2. Formulación débil usando r
Page 105 and 106:
⎡ σ = ⎢ ⎣ ⎤ σ 11 σ 22 σ
Page 107 and 108:
Figura 7 Teoría de placas clásica
Page 109 and 110:
Los esfuerzos de corte transversal
Page 111 and 112:
Capítulo 6 Elementos finitos en do
Page 113 and 114:
el punto p (ξ, η), y calculamos s
Page 115 and 116:
Para el caso cuadrático (ξ, η) =
Page 117 and 118:
∂u ∂x 2 = NN∑ I=1 [ ∂N I
Page 119 and 120:
Figura 5 Mapeamientos para funcione
Page 121 and 122: Elemento (L ( 1 , L 2 , L 3 ) w l L
Page 123 and 124: ahora dependiente sólo del valor d
Page 125 and 126: ∫ v ∫ σ ij δε ij dv = v δε
Page 127 and 128: 6.8.1. Funciones de interpolación,
Page 129 and 130: Siendo todas las variables constant
Page 131 and 132: (el I por ejemplo) con dicha fuente
Page 133 and 134: Capítulo 7 Aspectos generales asoc
Page 135 and 136: 7.2. Imposición de restricciones n
Page 137 and 138: entonces ¯k ij = k ij + 1 ǫ a ia
Page 139 and 140: Figura 2 viga flexible entre column
Page 141 and 142: D Matriz diagonal (definida positiv
Page 143 and 144: 7.5. Suavizado de Variables En el m
Page 145 and 146: 138
Page 147 and 148: Difícilmente se pueda encontrar un
Page 149 and 150: que para el caso de placas esbeltas
Page 151 and 152: cualquier variación de deformacion
Page 153 and 154: donde se ha introducido la matriz B
Page 155 and 156: 148
Page 157 and 158: en tanto que las velocidades (conoc
Page 159 and 160: Interfaz 1 ξ en [ 1, 0] η en [ 1,
Page 161 and 162: Parte Difusiva ⎡ 1 ⎣ 2A s 2 2
Page 163 and 164: Figura 4 Borde donde el flujo σ ν
Page 165 and 166: φ=0 grad φ=0 φ=1 simetría Figur
Page 167 and 168: 162
Page 169 and 170: de ecuaciones), de tal forma de hac
Page 171: 3 4 7 3 6 5 8 9 6 1 4 2 1 5 2 Figur
Page 175 and 176: LM(ndofe) arreglo que relaciona cad
show all

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

Create successful ePaper yourself

Delete template?

Save as template?