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