15.05.2013 Views

Cálculo Numérico - Engenharia Civil UEM

Cálculo Numérico - Engenharia Civil UEM

Cálculo Numérico - Engenharia Civil UEM

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.

CAPÍTULO 5. SOLUÇÃO DE SISTEMAS LINEARES: MÉTODOS ITERATIVOS 186<br />

Aplicando a lei de Kirchoff, que diz que a soma algébrica das diferenças de potencial em qualquer<br />

circuito fechado é zero, achamos para as correntes i1, i2, i3:<br />

⎧<br />

⎨<br />

⎩<br />

6i1 + 10(i1 − i2) + 4(i1 − i3) − 26 = 0<br />

5i2 + 5i2 + 5(i2 − i3) + 10(i2 − i1) = 0<br />

11i3 + 4(i3 − i1) + 5(i3 − i2) − 7 = 0<br />

a) É possível aplicar ao sistema acima o método de Gauss - Seidel com convergência assegurada?<br />

Justifique.<br />

b) Se possível, obtenha a solução com erro relativo < 10 −2 .<br />

5.3 - Suponha uma barra de metal homogêneo, como na figura a seguir, onde AB = CD = 4 : AC =<br />

BD = 3.<br />

y<br />

C<br />

0 o<br />

A<br />

✻<br />

R<br />

0 o<br />

1 o<br />

A temperatura ao longo de AB, AC, BD é mantida constante e igual a 0 o C, enquanto que ao longo<br />

de CD ela é igual a 1 o C. A distribuição do calor na barra R obedece à seguinte equação:<br />

com as condições de contorno:<br />

D<br />

B<br />

0 o<br />

✲<br />

x<br />

∂2u ∂x2 + ∂2u = 0 , (5.25)<br />

∂y2 u(x, y) = 1 para 0 < x < 4 ,<br />

u(x, 0) = 0 para 0 < x < 4 ,<br />

u(0, y) = 1 para 0 < y < 3 ,<br />

u(x, y) = 0 para 0 < y < 3 .<br />

A solução numérica desse problema pode ser obtida considerando-se uma divisão do retângulo ABCD<br />

em retângulos menores a partir de uma divisão de AB em intervalos iguais de amplitude h e de uma<br />

divisão de CD em intervalos iguais de amplitude k,como é mostrado na figura a seguir:

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

Saved successfully!

Ooh no, something went wrong!