21.08.2013 Views

o método dos elementos de contorno aplicado à ... - Sistemas SET

o método dos elementos de contorno aplicado à ... - Sistemas SET

o método dos elementos de contorno aplicado à ... - Sistemas SET

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.

7.4 Processo Iterativo<br />

No caso <strong>de</strong> não-linearida<strong>de</strong> física, procura-se o equilíbrio em um dado incremento <strong>de</strong><br />

carga, através do Método <strong>de</strong> Newton-Raphson, consi<strong>de</strong>rando-se o nível <strong>de</strong> carga (Mext)<br />

constante, como é mostrado na figura (7.6). Nesse trabalho será consi<strong>de</strong>rado apenas o<br />

Método <strong>de</strong> Newton-Raphson Modificado, <strong>de</strong>scrito a seguir na figura (7.7).<br />

Momento (M)<br />

M ext(i+1)<br />

∆M<br />

e( ) 1<br />

M ext(i) = M int(i)<br />

Ψ 1<br />

1<br />

M int( i+<br />

1 )<br />

(1/r)i<br />

1<br />

0<br />

[ST]<br />

[ST]<br />

∆(1/r) 1<br />

(1/r) 1<br />

on<strong>de</strong>: [ST] n é a matriz tangente na iteração n,<br />

∆(1/r) 2<br />

(1/r) 2<br />

[ST] 2<br />

Ψ 2<br />

Curvatura (1/r)<br />

(1/r)i+1<br />

FIGURA 7.6 - Método <strong>de</strong> Newton Raphson Padrão<br />

(1/r) n é o vetor <strong>de</strong> curvaturas da iteração n,<br />

Ψ n é o vetor <strong>de</strong> resíduo da iteração n,<br />

Mext(i+1) é o vetor <strong>de</strong> momentos externos do incremento i+1,<br />

∆M en ( ) é o incremento <strong>de</strong> momentos elásticos da iteração n,<br />

n<br />

M int( i+1<br />

)<br />

i+1,<br />

156<br />

é o vetor <strong>de</strong> momentos internos verda<strong>de</strong>iros da iteração n, do incremento<br />

(1/r)i é o vetor <strong>de</strong> curvaturas do incremento i.<br />

Na primeira iteração do incremento (figura 7.6), tem-se um incremento <strong>de</strong> momentos<br />

e( )<br />

∆M 1 , que foi calculado elasticamente, consi<strong>de</strong>rando-se as relações <strong>de</strong> equilíbrio da<br />

estrutura. Se o mesmo não for totalmente elástico, <strong>de</strong>ve-se fazer novas iterações, até que se<br />

encontre um estado <strong>de</strong> momentos internos aproximadamente igual ao estado <strong>de</strong> momentos

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

Saved successfully!

Ooh no, something went wrong!