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.
1.4.2. Aproximación simultánea <strong>de</strong> la ecuación diferencial y <strong>de</strong> las condiciones <strong>de</strong><br />
bor<strong>de</strong><br />
En esta sección se admitirá como función <strong>de</strong> <strong>prueba</strong> una aproximación que no satisfaga i<strong>de</strong>nticamente<br />
las condiciones <strong>de</strong> bor<strong>de</strong>, lo que permitirá ampliar el rango <strong>de</strong> funciones admisibles<br />
M∑<br />
u ≃ û = a m φ m (1.43)<br />
m=1<br />
La expansión (1.43) no satisface alguna o ninguna <strong>de</strong> las condiciones <strong>de</strong> bor<strong>de</strong>, por lo que el residuo<br />
en el dominio<br />
R Ω = A(û) = L (û) + p en Ω (1.44)<br />
<strong>de</strong>be complementarse con un residuo en el bor<strong>de</strong><br />
R Γ = B(û) = M (û) + r en Γ (1.45)<br />
A fin <strong>de</strong> reducir los <strong>residuos</strong> en el dominio y en el bor<strong>de</strong>, se propone el siguiente planteo en <strong>residuos</strong><br />
pon<strong>de</strong>rados<br />
∫<br />
∫<br />
__<br />
W l R Ω dΩ + W l R Γ dΓ = 0 (1.46)<br />
Ω<br />
en don<strong>de</strong> las funciones <strong>de</strong> peso, W l y __ W l , pue<strong>de</strong>n ser elegidas en forma in<strong>de</strong>pendiente. Resulta claro<br />
que a medida que aumenta el número <strong>de</strong> funciones para las cuales se satisface la expresión (1.46),<br />
mejor será la aproximación <strong>de</strong> û a u, asumiendo que la ec. (1.43) se ha elegido en forma a<strong>de</strong>cuada.<br />
Sustituyendo la aproximación (1.43) en la ec. (1.46) resulta un sistema idéntico a (1.37) en el<br />
que la matriz <strong>de</strong> coeficientes y el término in<strong>de</strong>pendiente resultan<br />
∫<br />
∫<br />
__<br />
K lm = W l L (φ m ) dΩ + W l M (φ m ) dΓ<br />
1.4.2.1. Ejemplos<br />
Ω<br />
∫<br />
f l = −<br />
Ω<br />
Γ<br />
Γ<br />
∫<br />
W l p dΩ −<br />
Γ<br />
__<br />
W l r dΓ (1.47)<br />
Ejemplo N ◦ 3: Resolveremos nuevamente el Ejemplo 1 pero en este caso no satisfaceremos las<br />
condiciones <strong>de</strong> bor<strong>de</strong> con las funciones <strong>de</strong> <strong>prueba</strong>:<br />
M∑<br />
û = a m φ m ; φ m = x m−1 m = 1, 2...<br />
0<br />
m=1<br />
El bor<strong>de</strong> Γ son dos puntos (x = 0, x = 1) y la ec.(1.46) queda<br />
∫ 1<br />
[ __ ] [ __ ]<br />
W l R Ω dx + W l R Γ + W l R Γ<br />
Si tomamos W l<br />
x=0<br />
x=1<br />
= 0<br />
= φ l ; W __<br />
l = −φ l | Γ<br />
, la anterior queda<br />
∫ 1<br />
( ) d2û φ l<br />
dx − û dx − [φ 2 l û] x=0<br />
− [φ l (û − 1)] x=1<br />
= 0<br />
0<br />
Con tres términos, û = a 1 + a 2 x + a 3 x 2 , que conduce a un sistema Ka = f don<strong>de</strong><br />
⎡<br />
3<br />
3 − 2 ⎤<br />
⎡ ⎤<br />
2 3<br />
1<br />
K =<br />
3 4 1<br />
2 3 4<br />
; f =<br />
⎢ 1<br />
⎥<br />
⎣ ⎦<br />
⎢<br />
⎥<br />
⎣<br />
4 5 8<br />
⎦<br />
1<br />
3 4 15<br />
12