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.

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

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

Saved successfully!

Ooh no, something went wrong!