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 ...

efn.uncor.edu
from efn.uncor.edu More from this publisher
15.07.2014 Views

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

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!