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.

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!