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