Folha de Rosto - Sistemas SET - USP
Folha de Rosto - Sistemas SET - USP
Folha de Rosto - Sistemas 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.
108 Capítulo 4-Estratégias <strong>de</strong> Enriquecimento<br />
A matriz <strong>de</strong> rigi<strong>de</strong>z K acaba por apresentar as seguintes parcelas:<br />
K = K + K<br />
b h<br />
T<br />
Kb= VBˆ CBˆ<br />
on<strong>de</strong> (4.116)<br />
K = ∫ B CB rdA<br />
h<br />
A<br />
T<br />
h h<br />
A matriz constitutiva C relaciona tensões e <strong>de</strong>formações e apresenta-se no<br />
caso axissimétrico conforme (3.24).<br />
A K b representa, <strong>de</strong> modo análogo ao que se apresentou no item 4.1.2.1,<br />
uma matriz base. Tendo-se em vista que a matriz base possui modos espúrios <strong>de</strong><br />
<strong>de</strong>formação, a matriz K h é <strong>de</strong> dita <strong>de</strong> alta or<strong>de</strong>m e tem a finalida<strong>de</strong> <strong>de</strong> garantir um<br />
correto ‘rank’ a matriz <strong>de</strong> rigi<strong>de</strong>z.<br />
Des<strong>de</strong> que, <strong>de</strong>vido a axissimetria, há somente um movimento <strong>de</strong> corpo<br />
rígido livre, sendo este relativo à translação na direção axial (direção z ), e<br />
consi<strong>de</strong>rando-se que o elemento possui oito gruas <strong>de</strong> liberda<strong>de</strong>, a matriz <strong>de</strong> rigi<strong>de</strong>z K ,<br />
possui apenas sete auto-vetores não-nulos, isto é, a sua or<strong>de</strong>m (‘rank’) é sete. Como a<br />
matriz dos coeficientes elásticos é positiva <strong>de</strong>finida, e po<strong>de</strong>-se mostrar que ˆ B possui<br />
or<strong>de</strong>m quatro, <strong>de</strong>vendo, portanto, para a matriz K h apresentar or<strong>de</strong>m três.<br />
O vetor <strong>de</strong> forcas nodais equivalente, mostrado a seguir, é composto por<br />
duas parcelas, sendo a primeira oriunda das contribuições <strong>de</strong> forças eventualmente<br />
aplicadas no contorno do elemento e a segunda formada pela ação das forças<br />
volúmicas.<br />
∫ ∫ (4.117)<br />
f = NtrdΓ+ NρbrdA Γ<br />
A<br />
Resta que a formulação do elemento finito fica <strong>de</strong>pen<strong>de</strong>nte <strong>de</strong> uma<br />
a<strong>de</strong>quada escolha da matriz B h , <strong>de</strong> modo que garanta boas características <strong>de</strong><br />
aproximação e atenda a condição pétrea (equação (4.113)).<br />
4.2.5 A matriz B h<br />
A escolha a<strong>de</strong>quada da matriz B h rege a <strong>de</strong>terminação do elemento finito