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.

Notar que K es no simétrica aun cuando se utilizó el método <strong>de</strong> Galerkin. La solución es<br />

a 1 = 0,068 a 2 = 0,632 a 3 = 0,226<br />

La convergencia en la aproximación a las condiciones <strong>de</strong> bor<strong>de</strong> en x = 0 y x = 1 se observa en la<br />

siguiente tabla con aproximaciones usando 1,2 o 3 términos<br />

N ◦ términos 1 2 3 exacto<br />

x=0 1/3 −0,095 0,068 0<br />

x=1 1/3 0,762 0,925 1<br />

Ejemplo N ◦ 4: resolveremos nuevamente el ejemplo 2 utilizando nuevamente funciones pares<br />

pero relajando la condición <strong>de</strong> que las funciones <strong>de</strong> <strong>prueba</strong> satisfagan las condiciones <strong>de</strong> bor<strong>de</strong>.<br />

Elegimos el conjunto <strong>de</strong> funciones<br />

φ 1 = 4 − y 2 ; φ 2 = x 2 φ 1 ; φ 3 = y 2 φ 1<br />

φ 4 = x 2 y 2 φ 1 ; φ 5 = x 4 φ 1 . . .<br />

que para una aproximación en 5 términos conduce a<br />

ˆϕ = ( 4 − y 2) ( a 1 + a 2 x + a 3 x 2 + a 4 x 3 + a 5 x 4)<br />

que satisface las condiciones <strong>de</strong> bor<strong>de</strong> en y = ±2 pero no en x = ±3. En consecuencia, esta<br />

condición <strong>de</strong>be incorporarse al problema a través <strong>de</strong>l enunciado <strong>de</strong> <strong>residuos</strong> pon<strong>de</strong>rados<br />

∫ 3 ∫ 2<br />

( )<br />

∂ 2 ϕ<br />

∂x + ∂2 ϕ<br />

2 ∂y + 2 W 2 l dydx+<br />

∫ 2<br />

−2<br />

−3<br />

−2<br />

[ ] ∫ 2<br />

¯Wlˆϕ dy − [ ¯Wl ˆϕ ] dy = 0<br />

x=3 x=−3<br />

Adoptando W l = φ l ; ¯W l = φ l | Γ<br />

, resulta<br />

⎡<br />

⎤ ⎡<br />

6,7 −44 43,7 80 813,6<br />

−12 −333,6 105,6 355,2 −5748,7<br />

K =<br />

⎢ 11,7 −16 159,1 389,5 −547,2<br />

⎥<br />

⎣ 9,6 −163,2 433,4 1971,6 −3932,4 ⎦ ; f = ⎢<br />

⎣<br />

−237,6 −3156,7 432 1473,7 −46239,4<br />

−2<br />

12<br />

36<br />

16<br />

48<br />

194,4<br />

A continuación se presentan la máxima tensión <strong>de</strong> corte τ y momento actuante T obtenidos con<br />

2 y 5 términos<br />

N ◦ términos 2 5 exacto<br />

T 58,94 73,6 76,40<br />

|τ| 3,52 3,33 2,96<br />

1.4.2.2. Ejercicios.<br />

Ejercicio N ◦ 11: graficar la solución obtenida en el ejemplo 3 conjuntamente con la exacta y las<br />

obtenidas en el ejemplo 1.<br />

Ejercicio N ◦ 12: con los datos <strong>de</strong>l ejemplo 4 verifique el valor <strong>de</strong>l momento T que figura en la<br />

tabla y obtenido con 5 términos.<br />

Ejercicio N ◦ 13: Utilizar el método <strong>de</strong> Galerkin para resolver el siguiente ejercicio y comparar<br />

con la solución exacta<br />

d 2 u<br />

dx 2 + u = 0 con { u = 1 en x = 0<br />

u = 0 en x = 1<br />

⎤<br />

⎥<br />

⎦<br />

13

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

Saved successfully!

Ooh no, something went wrong!