download completo - SET - USP
download completo - SET - USP
download completo - SET - USP
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