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

Create successful ePaper yourself

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

on<strong>de</strong>: f kl<br />

β<br />

104<br />

∂w<br />

qP<br />

β θ<br />

∂x q<br />

πDr<br />

f<br />

*<br />

, mmkl<br />

( , ) kl ( )<br />

( ) =−<br />

1<br />

3 (5.75)<br />

4<br />

β<br />

( θ)<br />

são funções das <strong>de</strong>rivadas <strong>de</strong> r em relação <strong>à</strong>s direções x1 e x2 e, portanto,<br />

<strong>de</strong>pen<strong>de</strong>m apenas do ângulo θ.<br />

se:<br />

Substituindo-se (5.75) e (5.59) em (5.74) e fazendo-se a integração sobre r, obtém-<br />

θ 3<br />

q<br />

m 1<br />

⎡ ξ α α α ln R 2<br />

eijkl ( q) = ∫ fβkl<br />

( θ)<br />

− ⎢ + ( b cosθ + a senθ)<br />

4π⎣R2( θ)<br />

2A<br />

θ<br />

1<br />

θ 3<br />

q<br />

1 ⎡ ξ α α α ln R 1<br />

− ∫ f βkl<br />

( θ)<br />

⎢−<br />

+ ( b cos θ + a sen θ)<br />

4π⎣R1( θ)<br />

2A<br />

θ<br />

1<br />

( θ)<br />

( θ)<br />

⎤<br />

⎥ dθ<br />

+<br />

⎦<br />

⎤<br />

⎥dθ<br />

⎦<br />

(5.76)<br />

Estas integrais po<strong>de</strong>m ser feitas numericamente em relação a θ, a não ser quando<br />

R1(θ) for nulo. Assim, nesse caso, <strong>de</strong>ve-se verificar as duas parcelas da segunda integral <strong>de</strong><br />

(5.76), envolvendo R1(θ). A primeira parcela envolve a integral:<br />

θ 3<br />

∫<br />

θ1<br />

f<br />

βkl<br />

q<br />

ξ α<br />

( θ)<br />

− dθ<br />

R ( θ)<br />

1<br />

⎡ ⎤<br />

⎢ ⎥ , (5.77)<br />

⎣ ⎦<br />

que apresenta singularida<strong>de</strong>, quando o ponto q coincidir com o vértice α do triângulo e,<br />

q<br />

portanto, ξ α = 1. Consi<strong>de</strong>rando-se, porém, a proprieda<strong>de</strong> da função f βkl<br />

( θ),<br />

dada por<br />

(5.73), a parcela (5.77) é nula, pois:<br />

θ 3<br />

∫ fβkl ( θ) dθ<br />

= 0<br />

(5.78)<br />

θ1<br />

A segunda parcela <strong>de</strong> (5.76), envolvendo R1(θ), que é dada por:<br />

θ 3<br />

∫<br />

θ1<br />

⎡<br />

fβkl ( θ) ⎢ b cosθ a senθ<br />

⎣<br />

α α<br />

( + )<br />

ln R 1<br />

2A<br />

( θ)<br />

⎤<br />

⎥ dθ<br />

, (5.79)<br />

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

Saved successfully!

Ooh no, something went wrong!