21.08.2013 Views

Modelos Não Lineares do Método dos Elementos de Contorno para ...

Modelos Não Lineares do Método dos Elementos de Contorno para ...

Modelos Não Lineares do Método dos Elementos de Contorno para ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

das gran<strong>de</strong>zas <strong>do</strong> problema físico. Dessa forma os parâmetros <strong>do</strong> problema analisa<strong>do</strong>,<br />

tais como <strong>de</strong>slocamentos e forças <strong>de</strong> superfície, não são utiliza<strong>do</strong>s diretamente e sua<br />

<strong>de</strong>terminação é obtida por integrações das funções fictícias <strong>de</strong>finidas na formulação.<br />

Essa formulação teve gran<strong>de</strong> aplicação em campos como eletromagnetismo e mecânica<br />

<strong>do</strong>s flui<strong>do</strong>s.<br />

5.1 – Soluções Fundamentais <strong>para</strong> os Problemas Bidimensionais e<br />

Tridimensionais<br />

Para o <strong>de</strong>senvolvimento da formulação <strong>do</strong> MEC torna-se necessária a obtenção<br />

<strong>de</strong> soluções fundamentais <strong>para</strong> as gran<strong>de</strong>zas envolvidas no problema. Para problemas<br />

elásticos estas soluções representam as respostas <strong>de</strong> um <strong>do</strong>mínio infinito sujeito a um<br />

carregamento unitário pontual localiza<strong>do</strong> em um ponto <strong>de</strong>nomina<strong>do</strong> ponto fonte,<br />

BREBBIA & DOMINGUEZ (1992).<br />

As soluções fundamentais <strong>de</strong> problemas elástico lineares são atribuídas a Kelvin<br />

e po<strong>de</strong>m ser obtidas por meio da equação <strong>de</strong> equilíbrio <strong>do</strong> corpo representada pela Eq.<br />

(I.2), <strong>do</strong> Anexo I. Inicialmente tem-se que a carga unitária pontual aplicada no ponto<br />

fonte é representada pela função Delta <strong>de</strong> Dirac, Anexo E. Levan<strong>do</strong> essa condição à<br />

equação <strong>de</strong> equilíbrio, Eq. (I.2), verifica-se que a função Delta <strong>de</strong> Dirac <strong>de</strong>ve ser igual<br />

ao carregamento <strong>de</strong> corpo b i .<br />

f<br />

bi ∆ ⋅<br />

Capítulo 5 – Méto<strong>do</strong> <strong>do</strong>s <strong>Elementos</strong> <strong>de</strong> <strong>Contorno</strong> __________________________________<br />

i<br />

92<br />

= e<br />

(5.1)<br />

f<br />

sen<strong>do</strong>: ∆ função <strong>de</strong>lta <strong>de</strong> Dirac aplicada no ponto fonte e e i o versor necessário <strong>para</strong><br />

<strong>de</strong>finição da direção <strong>de</strong> aplicação da carga.<br />

De posse <strong>de</strong>sse resulta<strong>do</strong> a equação <strong>de</strong> equilíbrio <strong>do</strong> corpo, Eq. (I.2), po<strong>de</strong> ser<br />

reescrita em termos <strong>de</strong> <strong>de</strong>slocamentos. Para realizar esse procedimento <strong>de</strong>ve-se<br />

empregar a Lei <strong>de</strong> Hooke generalizada, Eq. (I.5), a qual relaciona tensões a<br />

<strong>de</strong>formações, e em seguida aplicar as relações <strong>de</strong>formação-<strong>de</strong>slocamento, Eq. (I.7),<br />

ambas discutidas no Anexo I. Proce<strong>de</strong>n<strong>do</strong> <strong>de</strong>ssa forma a equação <strong>de</strong> equilíbrio em<br />

termos <strong>de</strong> <strong>de</strong>slocamento, <strong>para</strong> material isótropo sob esta<strong>do</strong> plano <strong>de</strong> <strong>de</strong>formação, po<strong>de</strong><br />

ser <strong>de</strong>finida como:<br />

f<br />

⎛ 1 ⎞<br />

∆ ⋅ei<br />

⎜ ⎟⋅<br />

i,<br />

ij + u j,<br />

ii + = 0<br />

⎝1−<br />

2⋅υ<br />

⎠<br />

µ<br />

u (5.2)

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

Saved successfully!

Ooh no, something went wrong!