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
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