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

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Figura 1<br />

Triángulo Maestro, y triángulo en el espacio coor<strong>de</strong>nado físico<br />

elementos tendrán como parámetros comunes las variables asociadas a sus nodos comunes<br />

(los que <strong>de</strong>finen geométricamente su contorno común) es necesario que el valor <strong>de</strong> la variable<br />

u a lo largo <strong>de</strong> dicho contorno común sólo <strong>de</strong>penda <strong>de</strong> las variables asociadas a los nudos<br />

que lo <strong>de</strong>finen. En consecuencia la función <strong>de</strong> interpolación <strong>de</strong>be valer 0 no sólo en los otros<br />

nudos (condición (a)) sino también a lo largo <strong>de</strong> el(los) lado(s) que no lo incluyan.<br />

La condición (b) asegura que si el valor <strong>de</strong> la variable es constante en todo los nodos, entonces<br />

es constante en todo el elemento. La condición (c) consecuencia <strong>de</strong> la anterior dice que, en<br />

tal caso, el gradiente en todo el elemento será cero.<br />

A<strong>de</strong>más resulta conveniente que las funciones <strong>de</strong> aproximación sean capaces <strong>de</strong> representar<br />

un estado <strong>de</strong> gradiente constante, que es el límite que <strong>de</strong>be alcanzarse cuando <strong>de</strong> refina la<br />

discretización.<br />

6.3. Elementos triangulares<br />

Empecemos viendo el elemento más sencillo para problemas planos que es el triángulo lineal.<br />

Definamos inicialmente un elemento maestro, en forma similar a como hicimos en el problema<br />

unidimensional. En este caso el elemento maestro se <strong>de</strong>fine como un triángulo rectángulo con<br />

el ángulo recto en el origen <strong>de</strong> coor<strong>de</strong>nadas, lados paralelos a los ejes y <strong>de</strong> longitud unitaria.<br />

Numeremos a<strong>de</strong>más sus vértices en la forma indicada, 1 <strong>de</strong> coor<strong>de</strong>nadas (ξ = 1, η = 0), 2 <strong>de</strong><br />

coor<strong>de</strong>nadas (ξ = 0, η = 1), y 3 <strong>de</strong> coor<strong>de</strong>nadas (ξ = 0, η = 0)<br />

En el elemento así <strong>de</strong>finido llamemos ξ al eje horizontal y η al eje vertical. Observemos las<br />

siguientes funciones lineales <strong>de</strong>finidas sencillamente como:<br />

L 1 (ξ, η) = ξ<br />

L 2 (ξ, η) = η<br />

L 3 (ξ, η) = 1 − ξ − η<br />

Notemos que estas funciones satisfacen todas las condiciones pedidas anteriormente. La (a) resulta<br />

inmediata <strong>de</strong> evaluar en los nudos, en tanto que para la (b) basta ver la <strong>de</strong>finición <strong>de</strong> la 3ra.<br />

función. Observemos a<strong>de</strong>más que el gradiente <strong>de</strong> la variable es constante para cualquier valor que<br />

tomen las parámetros nodales.<br />

Trataremos ahora <strong>de</strong> darle un significado geométrico a las funciones <strong>de</strong> forma L I , tomemos<br />

un punto cualquiera “p” <strong>de</strong> coor<strong>de</strong>nadas (ξ, η) <strong>de</strong>ntro <strong>de</strong> elemento y unamos este punto con los 3<br />

vértices lo que nos <strong>de</strong>fine 3 triángulos. Si observamos el triángulo inferior <strong>de</strong>finido por el eje ξ y<br />

106

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

Saved successfully!

Ooh no, something went wrong!