pratica5_NDOF aula com sistema.
mcd 1 1oSem2015
Arquivo que mostra como se calcula os parâmetros modais para um sistema com vários graus-de-liberdade:
ORIGIN 1 Mudança da origem do programa de 0 para 1.
Especificar matrizes de massa e rigidez do sistema, assim como o intervalo de linhas e colunas destas:
.5 0 0 0 0 0 0 0 0 .3 .1 .1 0 0 0 0 0 0
0 1.5 0 0 0 0 0 0 0 .1 .3 .1 .1 0 0 0 0 0
0 0 1 0 0 0 0 0 0 .1 .1 .3 .1 0 0 0 0 0
0 0 0 2 0 0 0 0 0 0 .1 .1 .3 .1 0 0 0 0
M 0 0 0 0 1.5 0 0 0 0 K 0 0 0 .1 .3 .1 .1 0 0 . 107 N rows ( M )
0 0 0 0 0 2 0 0 0 0 0 0 0 .1 .3 .1 .1 0
mt cols ( M )
0 0 0 0 0 0 .5 0 0 0 0 0 0 .1 .1 .3 .1 0
i 1 .. N
0 0 0 0 0 0 0 1 0 0 0 0 0 0 .1 .1 .4 .1
0 0 0 0 0 0 0 0 1.5 0 0 0 0 0 0 0 .1 .2 j 1 .. mt
val genvals ( K , M ) vec genvecs ( K , M ) Cálculos dos autovalores e autovetores:
T
sort csort augment vec , val , mt 1 Agrupar as matrizes acima para se ordenar do menor para o maior
autovalor
T
λ submatrix ( sort , 1 , N , mt 1 , mt 1) ψ submatrix ( sort , 1 , N , 1 , mt ) Extrair os autovalores e autovetores já
ordenados.
T
mm ψ . M. ψ Cálculo da matriz de massa generalizada para normalização dos autovetores
cho cholesky( Re ( mm) ) φ ψ . cho 1 Cálculo dos autovetores normalizados.
T
cond2 φ . M . φ = 1 Verificação da normalização em relação a matriz de massa
T
λ = 1.536 . 105 4.888 . 105 1.17 . 106 1.833 . 106 2.608 . 106 3.382 . 106 3.968 . 106 6.719 . 106 7.012 . 106
ωr λ T
ω r = 391.936 699.109 1.082 . 103 1.354 . 103 1.615 . 103 1.839 . 103 1.992 . 103 2.592 . 103 2.648 . 103
λ
freq T
.
2π freq = ( 62.378 111.267 172.174 215.477 257.02 292.672 317.015 412.541 421.44 )
Para se ver que a variação de uma massa apenas não afeta da mesma forma todas as frequências naturais, f.orig é o valor
das frequências se m6 = 1kg.
Se m6 = 1 kg: f orig ( 64.937 121.766 174.507 226.008 273.93 292.784 345.203 412.543 423.868 )
Vibrações Mecânicas (Grad.) DEMEC/UFMG Profa. Maria Lúcia M. Duarte
pratica5_NDOF aula com sistema.mcd 2 1oSem2015
Logo, a variação que existirá quando se comparam as frequências originais com aquela onde o m6 = 2kg, será:
T
f orig freqi
i T 4
erroi erro = 3.94 8.623 1.337 4.66 6.173 0.038 8.166 5.189 . 10 0.573 %
T
f orig
i
Especificando os dados para obtenção das curvas de FRF:
fst 10 fend 300 nfp 800 p 1 .. ( nfp 1)
fend fst
finc fp fst finc . ( p 1) ωp fp . 2 . π finc = 0.363
nfp
Especificando os pontos para obtenção da curva: ip = resposta, jp = força
Como a análise foi feita sem considerar amortecimento, será considerado um amortecimento modal de 1% em
cada modo:
Colocar a variação de ηi = 1%, 2% e 5% para ver variação do Nyquist.
ip 4 jp 2 ηi 1 %
Fazer jp = 5 e 8 para observar o comportamento da curva de FRF (Hp e Ht).
N r 1 .. 9 ip = 4 jp = 2
φ ip , r . φ ip , r 2
Hpp . ω Hp dB 20 . log Hpp
2. . 2 p p φ ip , r . φ jp , r = φ ip , r . φ ip , r =
r=1 ω rr 1 i η r
ω p 1 1
1 0.124 1 0.138
N
φ ip , r . φ jp , r 2 0.076 2 0.049
Htp . ω 2 Ht dB 20 . log Htp 3 -3.194·10 -3 3 9.303·10 -4
2 2 p p
r=1 ωr . 1 i . ηr ωp 4 -0.154 4 0.167
r
5 -0.083 5 0.133
6 0.047 6 0.01
50
7 3.324·10 -6 7 3.831·10 -4
8 -6.592·10 -3 8 2.053·10 -3
9 -7.335·10 -5 9 2.514·10 -4
Hp dB
p
0
Ht dB
p
50
0 50 100 150 200 250 300
f
p
FRF de ponto
FRF de transferência
0 20
5 10
Im( Hp ) 10
Im( Ht ) 0
15
10
20
10 5 0 5 10
20
Re( Hp ) 10 0 10
Re( Ht )
Vibrações Mecânicas (Grad.) DEMEC/UFMG Profa. Maria Lúcia M. Duarte
