Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...
Capítulo 1 Métodos de residuos ponderados Funciones de prueba ... Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...
Figura 3 Cuadrado Maestro, y rectángulo en el espacio coordenado donde los nudos están en correspondencia con las siguientes coordenadas locales Nudo ξ η 1 -1 -1 2 1 -1 3 1 1 4 -1 1 y x 0 son las coordenadas del centro del rectángulo y a y b son la las longitudes de sus lados en las direcciones x 1 y x 2 respectivamente. El elemento rectangular más sencillo resulta entonces el bilineal obtenido de multiplicar los polinomios lineales en ambas direcciones o englobando las cuatro en una única expresión N 1 (ξ, η) = 1 (1 − ξ) (1 − η) 4 N 2 (ξ, η) = 1 (1 + ξ) (1 − η) 4 N 3 (ξ, η) = 1 (1 + ξ) (1 + η) 4 N 4 (ξ, η) = 1 (1 − ξ) (1 + η) 4 N I (ξ, η) = 1 4 ( 1 + ξ I ξ ) ( 1 + η I η ) De la misma forma puede encontrarse el elemento cuadrático Lagrangeano de 9 nodos y el cúbico de 16 nodos. Que las funciones propuestas cumplen con las condiciones expresadas inicialmente es muy fácil de demostrar y se deja como ejercicio. La utilización de las funciones de forma para rectángulo en términos de las coordenadas locales, permite escribir las derivadas en términos de la regla de la cadena, es decir dado: u (x) = NN∑ I=1 N I (ξ, η) u I 110 entonces ∂u ∂x 1 = NN∑ I=1 [ ∂N I ∂ξ (ξ, η) 2 a ] u I
∂u ∂x 2 = NN∑ I=1 [ ∂N I ∂η (ξ, η) 2 b Notemos que a diferencia del triángulo, en donde cuando se genera un elemento (cuadrático por ejemplo) aparece la cantidad exacta de coeficientes necesarios para dicha aproximación (6 en el caso cuadrático), para los elementos rectangulares aparece una cantidad de coeficientes (9 en el caso cuadrático) mayor que el número indispensable, asociados con términos de orden superior (x 2 1x 2 , x 2 1x 2 2, x 1 x 2 2). Además dado que en la matriz global de coeficientes, los parámetros asociados a los nudos internos del elemento sólo tienen contribución del mismo elemento, muchas veces suelen eliminarse estos grados de libertad por “condensación”. Estos considerandos han llevado a desarrollar elementos cuadráticos de mayor orden con sólo nudos en el contorno, estos elementos se conocen como “serendípitos” y se obtienen por inspección de las funciones de forma. Por ejemplo el elemento rectangular de 8 nodos en el que se ha eliminado el nudo central del elemento de 9 nodos y con él el término x 2 1x 2 2, de forma que del triángulo de Pascal sobreviven los siguientes 1 x 1 x 2 x 2 1 x 1 x 2 x 2 2 x 2 1x 2 x 1 x 2 2 ] u I La forma estándar de encontrar estas funciones de forma es escribirlas de la forma N I (x 1 , x 2 ) = a 1 + a 2 x 1 + a 3 x 2 + a 4 x 2 1 + a 5x 1 x 2 + a 6 x 2 2 + a 7x 2 1 x 2 + a 8 x 1 x 2 2 e imponer las condiciones correspondientes de que N ( I x J) = δ IJ que conduce a invertir un sistema de 8 × 8, que nos da simultáneamente los 8 coeficientes de cada una de las 8 funciones de forma. La otra forma es armarlas directamente inspeccionando la forma que deberían tener (de allí su nombre ‘serendipity’ en inglés). Estas funciones resultan de esta forma ( N I (ξ, η) = 1 4 1 + ξ I ξ ) ( 1 + η I η ) ( ξ I ξ + η I η − 1 ) nudos esquina ( N I (ξ, η) = 1 2 1 − ξ 2 1 + η I η ) nudos medios η = ±1 ( N I (ξ, η) = 1 2 1 + ξ I ξ ) (1 − η 2 ) nudos medios ξ = ±1 6.5. Mapeamiento de la geometría Hasta ahora hemos considerado elementos con geometrías sencillas. En principio hemos definido las funciones de forma a partir de elementos “maestros” definidos sobre un dominio normalizado. En el caso del triángulo ha sido posible pasar fácilmente a un elemento triangular general de lados rectos, en tanto que para el caso de elementos cuadriláteros nos hemos restringido a elementos rectangulares. Desde el punto de vista práctico el elemento rectangular resulta muy limitado para el modelado de geometrías reales, para lo cual si resulta conveniente el elemento triangular que es mucho más versátil en ese aspecto, sin embargo en ambos casos debe aproximarse el contorno mediante segmentos de recta. Los inconvenientes anteriores pueden resolverse si se recurre a mapear al elemento “maestro” sobre el plano x 1 − x 2 en una forma más general que la usada hasta ahora. La idea es interpolar la geometría del elemento usando aproximaciones similares a las usadas para interpolar las variables nodales, es decir si describimos la geometría del elemento mediante x (ξ, η) = NN∑ I=1 N I (ξ, η) x I Haciendo uso del concepto ya conocido de las funciones de forma resulta que: 111
- 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: Para el caso cuadrático (ξ, η) =
- 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
Figura 3<br />
Cuadrado Maestro, y rectángulo en el espacio coor<strong>de</strong>nado<br />
don<strong>de</strong> los nudos están en correspon<strong>de</strong>ncia con las siguientes coor<strong>de</strong>nadas locales<br />
Nudo ξ η<br />
1 -1 -1<br />
2 1 -1<br />
3 1 1<br />
4 -1 1<br />
y x 0 son las coor<strong>de</strong>nadas <strong>de</strong>l centro <strong>de</strong>l rectángulo y a y b son la las longitu<strong>de</strong>s <strong>de</strong> sus lados en las<br />
direcciones x 1 y x 2 respectivamente.<br />
El elemento rectangular más sencillo resulta entonces el bilineal obtenido <strong>de</strong> multiplicar los<br />
polinomios lineales en ambas direcciones<br />
o englobando las cuatro en una única expresión<br />
N 1 (ξ, η) = 1 (1 − ξ) (1 − η)<br />
4<br />
N 2 (ξ, η) = 1 (1 + ξ) (1 − η)<br />
4<br />
N 3 (ξ, η) = 1 (1 + ξ) (1 + η)<br />
4<br />
N 4 (ξ, η) = 1 (1 − ξ) (1 + η)<br />
4<br />
N I (ξ, η) = 1 4<br />
(<br />
1 + ξ I ξ ) ( 1 + η I η )<br />
De la misma forma pue<strong>de</strong> encontrarse el elemento cuadrático Lagrangeano <strong>de</strong> 9 nodos y el cúbico<br />
<strong>de</strong> 16 nodos. Que las funciones propuestas cumplen con las condiciones expresadas inicialmente es<br />
muy fácil <strong>de</strong> <strong>de</strong>mostrar y se <strong>de</strong>ja como ejercicio.<br />
La utilización <strong>de</strong> las funciones <strong>de</strong> forma para rectángulo en términos <strong>de</strong> las coor<strong>de</strong>nadas locales,<br />
permite escribir las <strong>de</strong>rivadas en términos <strong>de</strong> la regla <strong>de</strong> la ca<strong>de</strong>na, es <strong>de</strong>cir dado:<br />
u (x) =<br />
NN∑<br />
I=1<br />
N I (ξ, η) u I<br />
110<br />
entonces<br />
∂u<br />
∂x 1<br />
=<br />
NN∑<br />
I=1<br />
[ ∂N<br />
I<br />
∂ξ (ξ, η) 2 a<br />
]<br />
u I