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