04.07.2014 Views

download completo - SET - USP

download completo - SET - USP

download completo - SET - USP

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.

126<br />

Eduardo Toledo de Lima Junior & Wilson Sergio Venturini<br />

d<br />

Nas quais jk representa a parcela do tensor de tensões associada ao dano. Os núcleos<br />

*<br />

das integrais são as soluções fundamentais X do problema, ou as derivadas destas soluções. A<br />

variável T representa as forças de superfície. Os termos livres TL resultam da diferenciação dos<br />

*<br />

núcleos ijk (S,q) na primeira das equações (7). As soluções fundamentais utilizadas, bem como o<br />

desenvolvimento matemático dessa formulação são aqui omitidos, em virtude do formato compacto<br />

desta publicação.<br />

Deve-se representar o sistema de equações de forma discreta, sobre os elementos de<br />

contorno e sobre as células triangulares de domínio, a fim de se obter os valores aproximados das<br />

variáveis em questão. Define-se N<br />

n<br />

como o número de nós do contorno e N<br />

i<br />

o número de nós<br />

internos. A discretização do sistema (7), seguida de algumas manipulações algébricas características<br />

do MEC, conduz à seguinte equação condensada<br />

d<br />

<br />

<br />

E Ns Np QS<br />

<br />

3N i,3N 3N<br />

i<br />

i,1 3N i,1 3N 3N 3N i,3Ni<br />

i ,1<br />

i ,1<br />

1<br />

1 <br />

Npb<br />

QSIK I QP <br />

<br />

<br />

(i) <br />

Np<br />

3N i,3N 3N i,Ni<br />

3N i ,1<br />

i <br />

M t<br />

N i,N <br />

i N,1 i<br />

sendo<br />

1<br />

<br />

2 1 <br />

<br />

E b<br />

E QSIK I QP <br />

(i)<br />

QP <br />

(i) Tr<br />

<br />

<br />

<br />

<br />

3N 3N<br />

i,3Ni i,3N t<br />

<br />

i 3N 3N<br />

i,3Ni i,N t <br />

i<br />

<br />

N N<br />

i,N <br />

i N i,Ni<br />

i,3N<br />

<br />

<br />

M<br />

<br />

i<br />

<br />

<br />

Devido à presença dos termos de correção associados ao dano, a equação (8) se mostra nãolinear<br />

a cada incremento t . O procedimento de resolução segue uma metodologia do tipo Newton-<br />

Raphson. Define-se, para uma iteração i+1, o incremento de deformação da forma<br />

i1 i i<br />

feita através de uma expansão<br />

n n n<br />

, sendo a determinação do corretor i n<br />

i<br />

em série de Taylor, truncada na primeira variação, em torno de n<br />

equilíbrio na iteração i , como segue<br />

<br />

<br />

<br />

<br />

<br />

i<br />

n<br />

<br />

Y <br />

i<br />

i<br />

Yn n<br />

0<br />

<br />

i<br />

n<br />

<br />

<br />

<br />

<br />

ij<br />

<br />

<br />

(8)<br />

(9)<br />

Y , resíduo da equação de<br />

(10)<br />

4 RESULTADOS OBTIDOS OU ESPERADOS<br />

A formulação apresentada encontra-se em fase de implementação computacional. Os<br />

modelos isolados para análise de danificação em sólidos e de escoamento de fluidos em meios<br />

poroelásticos já foram implementados, apresentando bom desempenho e precisão nos resultados<br />

obtidos em alguns exemplos simulados. Espera-se em breve viabilizar a modelagem de problemas<br />

acoplados, com ênfase na aplicação à mecânica dos solos e difusão de contaminantes.<br />

Cadernos de Engenharia de Estruturas, São Carlos, v. 11, n. 53, p. 123-127, 2009

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

Saved successfully!

Ooh no, something went wrong!