Mathcad - Exemplo 6-6
Mathcad - Exemplo 6-6
Mathcad - Exemplo 6-6
mcd
m 1000 k 400 k
M 2M
m 0 k k
M K
0 2 .m k k
2) O acoplamento presente no sistema deve ser observado na EM padrão. Portanto, vemos que a matriz
de massa é diagonal. Desta forma, NÃO temos acoplamento de inércia ou dinâmico (visto que não
temos amortecimento no sistema também. Como a matriz de rigidez é NÃO diagonal, TEMOS
acoplamento de rigidez ou estático.
1) Pelo método da rigidez, tem de se multiplicar a EM padrão [M]{a}+[K]{x}={0} pela matriz de massa
inversa. Portanto:
1
0 Se tentarmos obter a matriz de flexibilidade pela inversa de [K] = [A]
1 m flexibilidade, vemos que a matriz [K] é singular. O determinante |K| = 0. Por
M isto, não existe flexibilidade. Basta ver o conceito de flexibilidade acima.
1
0
( 2 .m )
1
0
m k k
B . 0.4 0.4
1 k k B= ORIGIN 1
0 0.2 0.2
( 2 .m )
3) A equação de movimento padrão pelo método da rigidez será dada pela eq. (a) abaixo:
Xg1 , 1 = 0.894
(k λ .m ) 1.
k
m m Xg2 , 1 = 0.447
Z
1 . 1 .(k 2 .λ .m )
k Xg2 , 1
.
( 2 m) 2 m = 0.5
Xg1 , 1
k λ .m . k.
0 X1 X2 0 ou seja:
m m
λ 3 . ω λ
k
( 2 .m )
k
autovalores freq. naturais (rad/s) m k
X1 X2
0 0 k λ .m m
X1 X2
λ= ω= k 0 .m
0.6 0.775 m
para λ1=0 m
autovetores => matriz modal 1
X1 X2 X1
1 2 1
X Usando a solução ao lado, chegamos a matriz
1 1 modal [X] e podemos calcular a massa e rigidez generalizadas.
T k
mm X .M .X matriz de massa generalizada
m
para λ2=3k/2m X1 X2
T 3 k.
kk X .K.X matriz de rigidez generalizada k m
2m
3 m
3 . 10 0 0 0
mm = kk = 3
0 6 . 10
3 0 3.6 . 10 2
X1 2 .X 2 X2
1
kk 0 0
= corresponde a matriz de autovetores, onde a diagonal é igual aos autovetores e os
mm 0 0.6 demais valores são zero.
[m]-normalizando os modos de vibração. Iremos ver que a massa e rigidez generalizados dependem do
vetor modal utilizado, porém, o vetor modal [m]-normalizado será único. Antes iremos mostrar as matrizes de
massa e rigidez generalizadas utilizando o vetor modal {Xg}
T T KK 0.6 0
MM Xg . M . Xg KK Xg . K . Xg = aparece na ordem que os vetores
MM 0 0 modais foram encontrados.
3
1.2 . 10 0 720 0
MM = KK =
3 0 0
0 1.5 . 10 0.6
λg eigenvals( B ) λg =
0
720 0
Λ1 Λ1 = 0.6
1200 λg sort( λg ) λg =
0.6
3
3 mm2 , 2 = 6 .10
mm1 , 1 = 3 .10
3 3
MM1 , 1 = 1.2 .10 MM2 , 2 = 1.5 .10