3a Lista de Exerccios - PMR-2420
(Soluo de Equaes Diferenciais Parciais pelo Mtodo de
Elementos Finitos)
1) Aplicar o mtodo de elementos finitos (Galerkin) para calcular a distribuio de momento
fletor M(x) ao longo da viga mostrada na figura abaixo. Esta viga tem comprimento
unitrio e est sujeita carga distribuda F(x)=sen(x). Lembrar que a equao que rege o
problema :
d 2M
F ( x)
dx 2
e as condies de contorno so: M(0)=0 e M(1)=0. Considerar 4 elementos iguais e
montar as equaes em forma matricial.
F(x)
x=0 x=1
2) A equao diferencial que descreve o comportamento dinmico de uma barra sujeita a
trao ou compresso (barra de trelia) dada por:
2u u
A EA 0
t 2
x x
Supor E e A constantes nesse caso. Considere o elemento mostrado na figura onde os graus de
liberdade so representados por u1 e u2.
u1 u2
n 1 n 2
L
x
Utilizando a funo aproximadora linear:
u ( x) c1 c 2 x
a) determine as funes de forma N1(x) e N2(x) desse elemento
b) mostre que a matriz de massa e rigidez do elemento so iguais respectivamente :
AL 2 1 AE 1 1
[ M ]e [ K ]e
6 1 2 L 1
1
3) A equao diferencial que descreve o comportamento de uma viga dada por:
2v 2 2v
EI 2 q ( x, t )
t 2 x 2 x
Supor E e I constantes nesse caso. Considere o elemento mostrado na figura onde os graus de
v
liberdade so representados por u1, u2, 1 e 2, onde .
x
v1 v2
1 2
L
x
n 1 n 2
Utilizando a funo aproximadora representada pelo polinmio:
v ( x) c 0 c1 x c 2 x 2 c3 x 3
c) determine as funes de forma N1(x),N2(x),N3(x) e N4(x) desse elemento
d) mostre que a matriz de rigidez do elemento igual :
156 22 L 54 13L
2
3L2
AL 22 L 4L 13L
[ M ]e
420 54 13L 156 22 L
13L 3L 22 L
2
4 L2
12 6L 12 6L
2
6L 2 L2
EI 6 L 4L
[ K ]e 3
L 12 6L 12 6L
6L 2L 2
6L 4 L2
4) Conduo de calor num cilindro oco. A conduo de calor num meio contnuo
bidimensional regida pela seguinte equao elptica:
2T 2T
k2T Q k k Q 0 , (Equao de Poisson)
x 2 y 2
sendo T = T(x,y) a temperatura, k, a condutividade trmica do material (assumida constante) e
Q, o calor gerado internamente por unidade de rea do material (um valor positivo de Q indica
gerao de calor, um valor negativo implica extrao de calor).
O problema apresenta simetria axial e no ocorre variao de temperatura na direo
circunferencial, o que corresponde a adotar T/n=0, como indicado na figura acima. Alm
T
disso: q k , onde q o fluxo de calor.
n
Montar o sistema de equaes algbricas para a determinao das temperaturas nos ns 1, 2, 3
e 4, utilizando o Mtodo de Elementos Finitos.
5) Repita o problema acima para o caso onde a temperatura nas superfcies interna e externa
so fixadas em 200oC e 10oC respectivamente. Resolva o sistema e compare os resultados
com a soluo analtica:
190
T 200 ln r
ln 2
6) Problemas de toro elstica de barras prismticas so regidos pela equao:
2 2
2G
x 2 y 2
sendo G o mdulo elstico de cisalhamento e o ngulo de toro de cada seo por unidade
de comprimento da barra prismtica. x,y) a funo que representa a tenso mecnica, tal
que =0 nas fronteiras.
O momento de toro T dado por:
T 2 dx dy ,
e a tenso de cisalhamento numa direo n na seo pode ser obtida de
.
n
Pede-se determinar o valor de T e a tenso mxima de cisalhamento para um dado valor de
G e para a seo retangular mostrada na figura abaixo. Por consideraes de simetria,
suficiente obter a soluo apenas para um quarto da seo, como mostrado na figura.
Toro elstica de uma barra prismtica. (a) barra de seo quadrada; (b) malha de elementos
finitos usada.