17.04.2013 Views

Modelagem Física e Computacional de um Escoamento Fluvial

Modelagem Física e Computacional de um Escoamento Fluvial

Modelagem Física e Computacional de um Escoamento Fluvial

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

FIGURA 6.10 – QUADRATURA DE GAUSS COM NOVE PONTOS PARA ELEMENTO QUADRAN-<br />

GULAR<br />

FONTE: FROEHLICH, (2002)<br />

Ponto Coor<strong>de</strong>nadas Fator <strong>de</strong><br />

ξ η pon<strong>de</strong>ração wi<br />

1 -0,77459667 0,77459667 0,07716049<br />

2 0 0,77459667 0,12345679<br />

3 0,77459667 0,77459667 0,07716049<br />

4 -0,77459667 0 0,12345679<br />

5 0 0 0,19753086<br />

6 0,77459667 0 0,12345679<br />

7 -0,77459667 -0,77459667 0,07716049<br />

8 0 -0,77459667 0,12345679<br />

9 0,77459667 -0,77459667 0,07716049<br />

O número e a posição dos pontos <strong>de</strong> integração <strong>de</strong>termina a or<strong>de</strong>m do polinômio<br />

que po<strong>de</strong> ser integrado com exatidão. Por exemplo, a integração <strong>de</strong> nove pontos do<br />

elemento quadrangular permite a integração exata <strong>de</strong> <strong>um</strong> polinômio <strong>de</strong> quinta or<strong>de</strong>m<br />

em x e y. O RMA2 utiliza como padrão para elementos <strong>de</strong> lados retos a precisão <strong>de</strong><br />

quinta or<strong>de</strong>m. Para elementos com lados curvos a grau <strong>de</strong> precisão é a<strong>um</strong>entado para<br />

sétima or<strong>de</strong>m para acomodar os termos extras que aparecem durante a transformação<br />

para elementos com lados retos (KING, 1993).<br />

6.6 ESQUEMA DE SOLUÇÃO DO SISTEMA DE EQUAÇÕES<br />

Na aplicação do método dos elementos finitos, a matriz quadrada <strong>de</strong> coeficientes<br />

K do sistema <strong>de</strong> equações K U = f é geralmente esparsa 6 , dado que os termos fora da<br />

diagonal principal acoplando duas incógnitas são nulos a menos que essas incógnitas sejam<br />

comuns a <strong>um</strong> <strong>de</strong>terminado elemento.<br />

78<br />

É essencial que essa esparsida<strong>de</strong> seja totalmente<br />

explorada com o intuito <strong>de</strong> reduzir a necessida<strong>de</strong> <strong>de</strong> armazenamento computacional e o<br />

número total <strong>de</strong> operações matriciais necessárias à solução do sistema <strong>de</strong> equações (LEE;<br />

FROEHLICH, 1986).<br />

6 Uma matriz é dita esparsa quando o número <strong>de</strong> coeficientes nulos é muito superior ao número <strong>de</strong><br />

coeficientes não-nulos.

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

Saved successfully!

Ooh no, something went wrong!