15.07.2014 Views

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

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Un nudo <strong>de</strong>finido sobre el elemento maestro con coor<strong>de</strong>nadas ( ξ I , η I) se correspon<strong>de</strong>rá en el<br />

plano con el par coor<strong>de</strong>nado x I .<br />

Las coor<strong>de</strong>nadas <strong>de</strong>l contorno quedan <strong>de</strong>finidas exclusivamente por las coor<strong>de</strong>nadas <strong>de</strong> los<br />

nudos que forman el lado, lo que asegura continuidad <strong>de</strong> la geometría entre elementos sin<br />

solapamientos ni brechas<br />

Figura 4 Elemento finito Ω e en el plano (x, y) obtenido como la imagen <strong>de</strong>l mapeamiento T e <strong>de</strong>l<br />

correspondiente elemento maestro ¯Ω en el plano (ξ, η). También se indica el mapeamiento inverso Te<br />

−1<br />

<strong>de</strong> Ω e a ¯Ω.<br />

En el caso <strong>de</strong>l triángulo (lineal) <strong>de</strong> 3 nodos esto no representa ningún cambio práctico pero<br />

si <strong>de</strong>s<strong>de</strong> el punto <strong>de</strong> vista formal, en tanto que para el cuadrilátero bilineal si hay un cambio<br />

substancial que permite utilizar ahora elementos cuadriláteros <strong>de</strong> forma arbitraria (aunque veremos<br />

más a<strong>de</strong>lante que los ángulos interiores no <strong>de</strong>ben superar los 180 o en ningún caso) y ya no<br />

exclusivamente rectangulares.<br />

Resulta entonces posible ahora utilizar elementos <strong>de</strong> lados curvos (en elementos con más <strong>de</strong><br />

dos nudos en cada lado), en particular esto resulta útil en los contornos exteriores <strong>de</strong> la geometría<br />

<strong>de</strong>l problema, ya que para la interface entre elementos es conveniente utilizar contornos rectos.<br />

Veamos como interviene esta parametrización <strong>de</strong> la geometría (el dominio) en la generación <strong>de</strong><br />

las ecuaciones <strong>de</strong>l problema. Para la obtención <strong>de</strong>l gradiente <strong>de</strong> la variable tendremos ahora que:<br />

( ∂u<br />

∇u = , ∂u )<br />

=<br />

∂x 1 ∂x 2<br />

( NN∑<br />

I=1<br />

don<strong>de</strong> la aplicación <strong>de</strong> la regla <strong>de</strong> la ca<strong>de</strong>na supone<br />

)<br />

∂N I (ξ, η)<br />

NN∑<br />

u I ∂N I (ξ, η)<br />

,<br />

u I<br />

∂x 1 ∂x<br />

I=1 2<br />

112<br />

∂N I (ξ, η)<br />

= ∂N I (ξ, η) ∂ξ<br />

+ ∂N I (ξ, η) ∂η<br />

∂x 1 ∂ξ ∂x 1 ∂η ∂x 1

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

Saved successfully!

Ooh no, something went wrong!