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.
En la Figura 1 la <strong>de</strong>finición <strong>de</strong> la celda o volumen finito se ha hecho en función <strong>de</strong> las “medianas”<br />
<strong>de</strong> cada triángulo que parten <strong>de</strong> la mitad <strong>de</strong> cada lado y se cortan en el centro <strong>de</strong>l elemento. Esto<br />
hace que a lo largo <strong>de</strong> cada mediana dos <strong>de</strong> las coor<strong>de</strong>nadas triangulares (asociadas a los nudos<br />
extremos <strong>de</strong>l lado) tengan el mismo valor ( 1 sobre el lado y 1 en el centro <strong>de</strong>l elemento), en tanto<br />
2 3<br />
que la coor<strong>de</strong>nada triangular restante varíe entre 0 (sobre el lado) y 1 en el centro <strong>de</strong>l elemento.<br />
3<br />
Las coor<strong>de</strong>nadas <strong>de</strong>l centro <strong>de</strong>l triángulo son<br />
C x = 1 3<br />
(<br />
x 1 + x 2 + x 3)<br />
en tanto que las <strong>de</strong>l centro <strong>de</strong> cada lado son sencillamente (notar que con un supraíndice a la<br />
<strong>de</strong>recha se indican valores nodales, como se ha hecho habitualmente, y con supraíndice izquierdo<br />
se indican las coor<strong>de</strong>nadas <strong>de</strong>l centro <strong>de</strong>l lado, opuesto al nudo correspondiente)<br />
1 x = 1 (<br />
x 2 + x 3)<br />
2<br />
2 x = 1 (<br />
x 3 + x 1)<br />
2<br />
3 x = 1 (<br />
x 1 + x 2)<br />
2<br />
A partir <strong>de</strong> estas coor<strong>de</strong>nadas, cada interfaz queda <strong>de</strong>finida por el vector orientado <strong>de</strong>l centro<br />
<strong>de</strong>l triángulo al centro <strong>de</strong>l lado<br />
s 1 = 1 x− C x = − 1 3 x1 + 1 6 x2 + 1 6 x3<br />
s 2 = 2 x− C x = − 1 3 x2 + 1 6 x3 + 1 6 x1<br />
s 3 = 3 x− C x = − 1 3 x3 + 1 6 x1 + 1 6 x2<br />
A la longitud <strong>de</strong> cada interfaz la <strong>de</strong>nominarse como |s I | en tanto que el vector normal a la<br />
misma se expresa como<br />
n I = 1 ( ) (<br />
−s<br />
I<br />
|s I | 2 , s I 1 = n<br />
I<br />
1 , n2)<br />
I (9.1)<br />
<strong>de</strong> tal forma que el producto s I × n I sea saliente al plano.<br />
Figura 3<br />
η como función <strong>de</strong> ξ en cada interfaz interna<br />
Observando el elemento maestro y el elemento en el plano real resulta<br />
151