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