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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

p = F (1.29)<br />

en la que k y F no <strong>de</strong>pen<strong>de</strong>n <strong>de</strong> u. A su vez, las condiciones <strong>de</strong> bor<strong>de</strong> se pue<strong>de</strong>n expresar como<br />

B (u) = M (u) + r = 0 en Γ (1.30)<br />

en la M es un operador lineal a<strong>de</strong>cuado y r es in<strong>de</strong>pendiente <strong>de</strong> u. Por ejemplo, las condiciones<br />

<strong>de</strong> bor<strong>de</strong> <strong>de</strong> Dirichlet y Neumann asociadas con la ecuación (1.9) se escriben como<br />

M (u) = u r = − _ u en Γ u<br />

M (u) = −k ∂u<br />

(1.31)<br />

r = − σ _ en Γ σ<br />

∂n<br />

Siguiendo las i<strong>de</strong>as anteriores, se propone una solución <strong>de</strong> la forma (1.12)<br />

u ≃ û = ψ +<br />

M∑<br />

a m φ m (1.32)<br />

m=1<br />

en la que la funciones ψ y φ m son tales que<br />

M (ψ) = −r<br />

M (φ m ) = 0<br />

⎫<br />

⎬<br />

⎭<br />

en Γ (1.33)<br />

y en consecuencia û satisface las condiciones <strong>de</strong> bor<strong>de</strong> <strong>de</strong> la ecuación (1.28) para todos los valores<br />

<strong>de</strong> a m . Entonces, asumiendo que las funciones son suficientemente continuas y diferenciables en<br />

todo Ω, po<strong>de</strong>mos escribir<br />

∂u<br />

∂x ≃ ∂û<br />

∂x = ∂ψ<br />

M<br />

∂x + ∑ ∂φ<br />

a m<br />

m<br />

∂x<br />

∂ 2 u<br />

∂x ≃ ∂2 û<br />

2 ∂x = ∂2 ψ<br />

2<br />

m=1<br />

∂x 2 + M<br />

∑<br />

m=1<br />

a m<br />

∂ 2 m φ<br />

∂x 2 (1.34)<br />

y asi sucesivamente.<br />

Las condiciones <strong>de</strong> que las funciones <strong>de</strong> <strong>prueba</strong> sean continuamente diferenciables será relajada<br />

más a<strong>de</strong>lante.<br />

En lo que sigue, <strong>de</strong>bemos asegurar que û satisfaga solamente a la ecuación diferencial (1.28) ya<br />

que la expansión û se construyó satisfaciendo las condiciones <strong>de</strong> bor<strong>de</strong> (1.30). Sustituyendo û en<br />

la ec.(1.28) resulta el siguiente residuo<br />

R Ω = A(û) = L (û) + p = L (ψ) +<br />

M∑<br />

a m L (φ m ) + p (1.35)<br />

Este residuo pue<strong>de</strong> minimizarse, utilizando el método <strong>de</strong> <strong>residuos</strong> pon<strong>de</strong>rados, a fin <strong>de</strong> lograr que<br />

R Ω ≃ 0 en todo punto <strong>de</strong> Ω<br />

∫<br />

∫<br />

]<br />

M∑<br />

W l R Ω dΩ ≡ W l<br />

[L (ψ) + a m L (φ m ) + p dΩ = 0 (1.36)<br />

Ω<br />

Ω<br />

Evaluando la integral (1.36) para l = 1, 2...M, se obtiene un sistema <strong>de</strong> M ecuaciones algebraicas<br />

lineales<br />

Ka = f (1.37)<br />

m=1<br />

m=1<br />

7

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

Saved successfully!

Ooh no, something went wrong!