1 - Determine as equaes do movimento, as freqncias naturais e os vetores

modais do sistema abaixo. Escreva a expresso geral das respostas livres e plote as
resposta utilizando Matlab. (k=50.000 N/m; m=5 kg).x(0)=0,6 mm; v(0)=20 m/s.

Em srie
k 1=3/4 K

Em paralelo
k 1=3 K

Figura 1. Problema a ser determinado.

Dados:
k =50000 N /m
m=5 kg
x 0=0,0006 m
x 0=20 m/s

As equaes dinmicas do sistema so:

M x1 +3/4 K x 1 +3 K ( x 1x 2 ) =0
M x23 K ( x 1x 2 )=0
Na forma matricial, obtm-se a equao de estado:

]{ } [

M
0

]{ } { }

0 x 1 + 3,75 K 3 K x1 = 0
M x 2
3 K
3 K x2
0

Sabendo que a resposta do sistema pode ser descrita pela equao abaixo, a acelerao
obtida com a sua segunda derivada:
x i=X i sen ( n t )
x i=2n X i sen ( n t )
Substituindo os termos acima na equao de estado e eliminando o seno, teremos:

2n M
0

]{ } [

]{ } { }

X 1 + 3,75 K 3 K X 1 = 0
3 K
3 K X2
0
M X 2
2
n

Organizando a equao acima, teremos:

]{ } { }

2n M +3,75 K
3 K
X1 = 0
0
3 K
2n M +3 K X 2

Resolvendo a equao, igualando o determinante esquerda do vetor de amplitude

zero, teremos:
det

2n M +3,75 K

3 K
2
n

3 K

M +3 K

= M 2 4n6,75 KM 2n +2,25 K 2=0

Fazendo que s=2n :

2 2

M s 6,75 KMs+ 2,25 K =0

Resolvendo a equao de segundo grau, teremos:
s=

6,75 KM 6,75 2 K 2 M 2 42,25 K 2 M 2

2 M2
2 M2

Simplificando, teremos:
s=

K
( 3,375 3,023 )
M

Substituindo os valores das constantes, temos:

s=

50000
( 3,375 3,023 )
5

Obtm-se que s=63983,4666 e s=3516,5334, logo n1 =59,3 rad/s e n2 =252,95

rad/s.Substituindo a primeira frequncia natural na equao dinmica, obtm-se que:

(59,325+3,7550000 ) X 1( 350000 ) X 2=0

X1
0,883
=
. Substituindo a segunda frequncia natural na
X 2 =
1
equao dinmica, obtm-se a equao abaixo:
Obtm-se que

( )

n1

(252,9525+3,7550000 ) X 1( 350000 ) X 2=0

Obtm-se que

1= 0,8828
1

X1
X2

( )

n= n

=
2

1,133
. Portanto, os autovetores do sistema so obtidos:
1

2= 1,1328
1

A resposta do sistema ento dada por:

{ }

x 1 ( t ) = A sen t + + B sen t+
( n1
( n2
1
1)
2
2)
x2 ( t )

Derivando a equao acima, obtm-se que:

{ }

x1 ( t )
= A n1 1 cos ( n1 t + 1 ) + B n2 2 cos ( n2 t + 2)
x2 ( t )

Fazendo que t=0 e substituindo os valores iniciais do sistema, as equaes ficam como:

{ }{ }{ }
{ }{ } { }

x 1 ( 0 ) = 0,0006 = 0,8828 A sen + 1,1328 B sen

( 1)
( 2)
0,0006
1
1
x2 ( 0)

x1 ( 0 )
= 20 =59,3 0,8828 A cos ( 1 ) +252,95 1,1328 B cos ( 2 )
20
1
1
x2 ( 0 )

Fazendo que X =A sen ( 1 ) , Y =B sen ( 2) , Z =A cos ( 1 ) , W = A cos ( 1 ) ,

obtemos o sistema:

][ ] [ ]
[ ][ ]

0,8828 1,1328
0
0
X
0,0006
1
1
0
0
Y = 0,0006
0
0
52,35 286,54 Z
20
0
0
59,3
252,95 W
20

Obtemos ento que:

Elevando

X
6,3489E-4
Y = 3,4894E-5
Z
0,3569
W
0,0046

A sen ( 1 ) e A cos ( 1 ) ao quadrado e somando os dois, obtm-se que:

A 2= X 2+ Z 2 . : A= X 2 +Z 2= ( 6,3489E-4 ) + 0,35692=0,3569
2

Realizando o mesmo processo, para B sen ( 2 ) e B cos ( 2 ) , obtm-se que

B=0,0046 . Para se determinar as fases, temos que:
A sen ( 1 ) =X .: 1 =sen1

=0,0018 rad
( XA )=sen ( 6,3489E-4
0,3569 )
1

Realizando o mesmo processo, para B sen ( 2 ) , obtm-se que 2=0,0076 rad .

A resposta pode ento ser obtida graficamente na figura abaixo.

0.4
Massa 1
Massa 2

0.3

Deslocamento (m)

0.2
0.1
0
-0.1
-0.2
-0.3
-0.4

0.05

0.1

0.15

0.2

0.25
0.3
Tempo (s)

0.35

0.4

0.45

0.5

Figura 2. Resposta do sistema no domnio do tempo.

2- Determine as equaes do movimento, as freqncias naturais e os vetores
modais do sistema abaixo. Escreva a expresso geral das respostas livres e plote as
resposta utilizando Matlab. (k=10.000 N/m; m=2 kg).x(0)=0,6 mm; v(0)=20 m/s.
Neste sistema, despreza-se o atrito.

Figura 3. Problema a ser determinado.

Dados:
k =10000 N /m
m=2 kg
x 0=0,0006 m
x 0=20 m/s

As equaes dinmicas do sistema so:

2 M x1 + K x 1 + K ( x 1x 2 )=0
M x2K ( x 1x 2 ) +3 K x 2=0
Na forma matricial, obtm-se a equao de estado:

]{ } [

2M
0

]{ } { }

0 x1 + 2 K K x 1 = 0
M x2 K 4 K x 2
0
Sabendo que a resposta do sistema pode ser descrita pela equao abaixo, a acelerao
obtida com a sua segunda derivada:
x i=X i sen ( n t )
x i=2n X i sen ( n t )
Substituindo os termos acima na equao de estado e eliminando o seno, teremos:

2n 2 M
0

]{ } [

]{ } { }

X 1 + 2 K K X 1 = 0
0
M X 2 K 4 K X 2
2
n

Organizando a equao, teremos:

]{ } { }

2n 2 M +2 K
K
X1 = 0
2
0
K
n M + 4 K X 2

Resolvendo a equao igualando o determinante zero, teremos:

det

2n 2 M +2 K
K

K
2
n

M +4 K

=2 M 2 4n10 KM 2n +7 K 2=0

Fazendo que s=2n :

2 M 2 s 210 KMs+7 K 2=0
Resolvendo a equao acima, teremos:
s=

10 KM 102 K 2 M 2427 K 2 M 2

4 M2
4 M2

Simplificando, teremos:
s=

K
( 2,5 1,6583 )
M

Substituindo os valores das constantes, temos:

s=

10000
( 2,5 1,6583 )
2

Obtm-se que s=20791,5 e s=4208,5, logo n1 =64,87rad/s e n2 =144,19rad/s.

Substituindo a primeira frequncia natural na equao de estado, obtm-se a Eq. :

(64,87222+210000 ) X 1( 10000 ) X 2=0

X1
3,157
=
. Substituindo a segunda frequncia natural na
X 2 =
1
equao de estado, obtm-se a equao abaixo:
Obtm-se que

( )

n1

(144,19222+ 210000 ) X 1 ( 10000 ) X 2=0

Obtm-se que

X1
X2

( )

n= n

0,1583
. Portanto, os autovetores do sistema so obtidos:
1

{ }

1= 3,157
1

2= 0,1583
1

A resposta do sistema ento dada por:

{ }

x1 ( t )
= A 1 sen ( n1 t + 1 ) + B 2 sen ( n2 t+ 2 )
x2 ( t )

Derivando a equao acima, obtm-se que:

{ }

x1 ( t ) = A cos t + + B cos t+
( n1
( n2
n1 1
1)
n2 2
2)
x2 ( t )

Fazendo que t=0 e substituindo os valores iniciais do sistema, as equaesficam como:

{ }{ }{ }
{ }{ } { }

x1 ( 0)
= 0,0006 = 3,157 A sen ( 1 ) + 0,1583 B sen ( 2 )
0,0006
1
1
x2 ( 0)

x1 ( 0 ) = 20 =64,87 3,157 A cos +144,19 0,1583 B cos

( 1)
( 2)
20
1
1
x2 ( 0 )

Fazendo que X =A sen ( 1 ) , Y =B sen ( 2) , Z =A cos ( 1 ) , W = A cos ( 1 ) ,

obtemos o sistema:

][ ] [ ]

3,157 0,1583
0
0
X
0,0006
1
1
0
0
Y = 0,0006
0
0
204,79 22,825 Z
20
0
0
64,87
144,19 W
20

[ ][ ]
X
2,0955E-04
Y = 3,9045E-04
Z
0,1077
W
0,0903

Obtemos ento que:

Elevando

A sen ( 1 ) e A cos ( 1 ) ao quadrado e somando os dois, obtm-se que:

A 2= X 2+ Z 2 . : A= X 2 +Z 2= ( 2,0955E-04 ) + 0,10772=0,1077
2

Realizando o mesmo processo, para B sen ( 2 ) e B cos ( 2 ) , obtm-se que

B=0,0903 . Para se determinar as fases, temos que:
1

A sen ( 1 ) =X .: 1 =sen

( 2,0955E-04 )
X
=sen1
=0,0019 rad
A
0,1077

( )

Realizando o mesmo processo, para B sen ( 2 ) , obtm-se que 2=0,0043 rad . A

resposta pode ento ser obtida graficamente na figura abaixo.
0.4
Massa 1
Massa 2

0.3

Deslocamento (m)

0.2
0.1
0
-0.1
-0.2
-0.3
-0.4

0.05

0.1

0.15

0.2

0.25
0.3
Tempo (s)

0.35

0.4

0.45

0.5

Figura 4. Resposta do sistema no domnio do tempo.

3 - No sistema mostrado abaixo, a massa m1 excitada por uma fora harmnica
tendo um valor mximo de 50 N e uma freqncia de 2 Hz. Determine a expresso
da vibrao forada de cada massa, sendo m1=10 kg, m2=5kg, k1=8000N/m, e
k2=2000 N/m. Plote os grficos utilizando o Matlab.

Figura 5. Problema a ser determinado.

Dados:
k 1=8000 N /m
k 2=2000 N /m
m1=10 kg
m2=5 kg
Fo =50 N
=4 rad / s
As equaes dinmicas do sistema so:
m1 x1+ k 1 x 1 +k 2 ( x 1x 2 )=F o sen ( t )
m2 x2k 2 ( x 1x 2 )=0
Na forma matricial, obtm-se a equao de estado:

]{ } [

]{ } {

m1 0 x1 k 1+ k 2 k 2 x 1
+
= Fo sen ( t )
0 m2 x2
k 2
k2 x2
0

Sabendo que a resposta do sistema pode ser descrita pela equao abaixo, a acelerao
obtida com a sua segunda derivada:
x i=X i cos ( t )
x i=2 X i cos ( t )
Substituindo as equaes acima na equao de estado, teremos:

]{ } {

2 m1 +k 1 +k 2
k 2
X 1 = F o sen ( t )
2
k 2
m2 + k 2 X 2
0

A inversa da matriz esquerda do vetor de amplitude igual a funo de transferncia

que :

2 m1 +k 1 +k 2
k 2
k 2
2 m2 + k 2

2
k2
1 m2+ k 2
2
D ()
k2
m1 +k 1 +k 2

onde
4

D ( )=m1 m2 [ k 2 m1 + ( k 1+ k 2 ) m2 ] + k 1 k 2
H

Substituindo os valores, tm-se que a funo de transferncia

H= 0,0002 0,0003
0,0003 0,0014

igual a:

Logo, pode-se calcular a resposta do sistema, multiplicando a matriz de transferncia

pela fora:

{ } {

X 1 =H 50 sen ( 4 t )
X2
0

Obtm-se finalmente que a resposta do sistema dada por:

{ }{

X 1 = 0,0098 sen ( 4 t )
X2
0,0161 sen ( 4 t )

A resposta pode ento ser obtida graficamente na figura abaixo.

0.02
Massa 1
Massa 2

0.015

Deslocamento (m)

0.01
0.005
0
-0.005
-0.01
-0.015
-0.02

0.2

0.4

0.6

0.8

1
1.2
Tempo (s)

1.4

1.6

1.8

Figura 6. Resposta do sistema no domnio do tempo.

4 - Uma mquina em repouso pesando 2,1352104 N observada defletir 3,0510-2
m quando colocada sobre os seus isoladores tipo mola. Em operao, uma fora
harmnica de amplitude mxima de 444,8 N induz ressonncia no sistema. (i)
Projete um neutralizador de vibraes sem amortecimento de forma que o mesmo

tenha uma amplitude mxima de 410-3 m; (ii) Qual o valor da relao das
massas?
Dados:
X m =4103 m
mm=2176,554 kg
F0 =444,8 N
A frequncia natural calculada:
n=

g
9,81
rad
=
=17,934
2
est
s
3,0510

A rigidez mxima do neutralizador obtida:

X 2=

F0
F
444,8
N
.:k 2= 0 =
=111200
3
k2
X 2 410
m

Considerando que o neutralizador deve possuir a mesma frequncia natural do sistema,

obtm-se a massa do mesmo:
n2 =

k2
k
. : m 2= 2n =345,7 kg
m2
n2

m2
345,7
=
=0,15
mm 2176,554
projeto satisfatrio.

com a razo de massa est entre 0,05 0,15 0,25 o

5 - Deduzir o modelo matemtico para o sistema da figura abaixo, usando as

coordenadas generalizadas x(t) e (t). Represente as equaes na forma matricial.
Adote (m=1,0 kg e k=1000 N/m)

Figura 7. Problema a ser resolvido.

Utilizando o mtodo da energia, tem-se que a energia cintica total do sistema consiste
da soma de parcelas relativas translao e rotao como segue.
2

1
1 d ( x +lcos ( ) ) 1 d ( lsen ( ) )
T cinetica= m x2 + m
+ m
2
2
dt
2
dt

(1.3)

se n + cos2

2 x lsen()
1
1
1
T cinetica= m x2 + m x2 + ml 2
2
2
2

(1.3)

Para pequenos deslocamentos de

reescrita como:

, sen ( )

ento a equao acima pode ser

2 1
2
T cinetica=m x + m l
2

(1.3)

A energia potencial consiste na energia potencial do mola e do pendulo, de modo que a

parcela referente ao bloco no necessria uma vez que ele se encontra em uma posio
de equilbrio esttico:
1
U= k x 2 +mgl (1cos ( ))
2

(1.3)

1 2
Considerando a relao geomtrica cos ( ) =1 :
2
1 2
2
U= k x +mgl
2

(1.3)

Considerando a equao de Lagrange para a energia de um sistema qualquer:

d dT
dT dU
+
=0
dt dx
dx dx

( )

(1.3)

Inserindo ento U e T na equao, temos as equaes dinmicas

2 m x + kx=0

ml 2 +mgl=0

(1.3)

Que podem ser reescritas na seguinte forma matricial:

]{ } [ ] { } { }

2m 0 x
+ k 0 x =0
2

0 ml 0 ml
0

(3.2)

6 - Uma ponte rolante consiste de uma viga (rigidez flexo EI= 17,22109 Nm2 e
comprimento L = 9,15 m), um carro (massa m1 = 3627,5 kg) e um cabo de ao (E

= 2,071011 N/m2, rea da seo reta A = 7,110-4 m2e comprimento l = 6,1 m),
conforme ilustra a figura (a). A figura (b) mostra o modelo adotado. Se a carga
suportada de m2 = 907 kg, calcular as frequncias naturais do sistema.

Figura 8. Problema a ser resolvido.

Dados:
m1=3627,5 kg
m2=907 kg
EA 2,07 x 1011 x 7,1 x 104
N
K cabo =k c =
=
=24,09 x 106
l
6,1
m
9
48 EI 48 x 17,22 x 10
9 N
K viga=k v = 3 =
=1,0789 x 10
3
m
l
9,15
As equaes dinmicas para o sistema so obtidas por:
m1 x1+ k v x 1+ k c ( x 1x 2 )=0
m2 x2 +k c ( x 2x1 ) =0
A equao de estado do sistema dada por:

]{ } [

]{ } { }

m1 0 x1 k c + k v k c x1
+
=0
0 m2 x2
k c
k c x2
0

Resolvendo a equao, a determinante tirada e igualada a zero:

2 m1 +k c + k v
k c
det
=0
k c
2 m2+ k c

m1 m2 ( m2 k v +m2 k c + m1 k c ) + k v k c =0
Resolvendo a equao acima, obtm-se que

1=25920,29. : n1 =160,99rad /s
2=302249,7157.: n2=549,77 rad /s
7 - Determinar as frequncias associadas ao movimento angular (arfagem) e linear
vertical (oscilao) e a localizao dos centros de oscilao de um automvel da
Figura abaixo com os seguintes dados:
massa = m = 1000 kg
raio de girao = r = 0,9 m
distncia entre eixo dianteiro e C.G. = l1 = 1,0 m
distncia entre eixo traseiro e C.G. = l2 = 1,5 m
rigidez das molas dianteiras = kf = 18 kN/m
rigidez das molas traseiras = kr = 22 kN/m

Figura 9. Problema a ser resolvido.

Dados:
m=1000 kg
r=0,9m
l 1=1 m
l 2=1,5 m
k 1=18000 N /m
k 2=22000 N /m
m
1000
2
2
J o= ( l 1 +l 2 ) +mr =
(1+1,5 )2+10000,92=1330,833 kg m2
12
12
As equaes dinmicas para o sistema so obtidas por:
m x +k 1 ( xl 1 ) +k 2 ( x +l 2 ) =0
1 ( xl 1 ) l 1+ k 2 ( x+l 2 ) l 2=0
J o k
A equao de estado do sistema dada por:

[ ]{ } [

]{ } { }

( k 1 l 1k 2 l 2 ) x
( k 1 +k 2 )
m 0 x +
=0
2
2

0 J o ( k 1 l 1k 2 l 2) ( k 1 l 1 +k 2 l 2)
0

Resolvendo a equao, a determinante tirada e igualada a zero:

det

2 m+ k 1+ k 2

( k 1 l1 k 2 l2 )

( k 1 l 1k 2 l 2 )

J o + ( k 1 l 1 +k 2 l 2 )

=0

m J o m ( k 1 l 1+ k 2 l 2 ) + J o ( k 1+ k 2 ) + k 1 k 2 ( l 1 l 2) =0

Resolvendo a equao acima, obtm-se que

1=89,893. : n1=9,48 rad /s
2=0,8275. : n2=0,9097 rad / s
8 - Determinar a resposta de regime permanente do sistema mostrado na Fig. 5.11
em que a massa m1 sofre a ao da fora F1(t) = F10 cos t. Traar a curva de
resposta em frequncia.

Figura 10. Problema a ser resolvido.

Como a questo no possui dados de massa e rigidez das molas, utiliza-se:
Dados:
k =10000 N /m
m=5 kg
Fo =50 N

As equaes dinmicas do sistema so:

m x1 +k x 1+ k ( x1 x2 ) =F o sen ( t )
m x2k ( x 1x 2 ) + k x 2=0
Na forma matricial, obtm-se a equao de estado:

[ ]{ } [

]{ } {

m 0 x1 + 2k k x1 = Fo sen ( t )
0 m x2 k 2 k x2
0

Sabendo que a resposta do sistema pode ser descrita pela equao abaixo, a acelerao
obtida com a sua segunda derivada:
x i=X i cos ( t )
2
x i= X i cos ( t )

Substituindo as equaes acima na equao de estado, teremos:

]{ } {

2 m+2 k
k
X 1 = F o sen ( t )
2
k
m+2 k X 2
0

A inversa da matriz esquerda do vetor de amplitude igual a funo de transferncia

que :

2 m+2 k
k
2
k
m+2 k

1 2 m+ 2 k
k
2
D ()
k
m+ 2 k

onde
D ( )=m2 44 km 2 +3 k 2
Logo, pode-se calcular a resposta do sistema (amplitude), multiplicando a matriz de
transferncia pela matriz de fora, portanto:

{ }

F o ( m+ 2 k )

2 4
2
2
X1
= m 4 km +3 k
X2
Fo k
2

m 4 km +3 k

A resposta pode ento ser obtida graficamente na figura abaixo.

10

Massa 1
Massa 2
0

Amplitude (m)

10

-1

10

-2

10

-3

10

-4

10

8
10
Frequncia (Hz)

12

14

16

Figura 11. Resposta do sistema no domnio da frequncia.

9 - Um motor diesel, pesando 3000 N, est montado sobre um pedestal. Observouse que o motor induz vibraes na rea atravs do pedestal quando opera a 6000
rpm. Determinar os parmetros do absorvedor que reduziro as vibraes do
motor quando montado no pedestal. A magnitude da fora que produz a vibrao
no motor 250 N, e a amplitude do movimento da massa auxiliar est limitada em
2mm.
Dados:
mm=3000 /9,81=305,81 kg
60002
nm=
=628,318 rad /s
60
F0 =250 N
max
X 0,002 m

A rigidez do sistema calculada:

nm=

km
.: k m=nm2 mm=628,3182305,81=1,2073E8 N /m
mm

As equaes dinmicas do sistema so:

( x1 x2 ) =F o sen ( t )

||

mm x1+ k m x1 + k

( x1 x2 ) =0

||

|| x

k
m
2

Na forma matricial, obtm-se a equao de estado:

||
||
||

k m + k k k k|| x 1 = F o sen ( t )
x2
0
mm 0 x1
+
0 m|| x2

{}{
{}

Sabendo que a resposta do sistema pode ser descrita pela equao abaixo, a acelerao
obtida com a sua segunda derivada:
x i=X i cos ( t )
x i=2 X i cos ( t )
Substituindo as equaes acima na equao de estado, teremos:
||
||
||
||+k||

X1
X2

2 mm + k m +k k k 2 m { }= F o sen ( t )
0

A inversa da matriz esquerda do vetor de amplitude igual a funo de transferncia

que :
||
||
||
||+k||

2 mm +k m +k k k 2 m
1
H=

D()
onde

||
||m
||

2 +k m k
|| m + k m
m
m
k
4
||
D ( ) =m m m

||+k

Logo, pode-se calcular a resposta do sistema (amplitude), multiplicando a matriz de

transferncia pela matriz de fora, portanto:

||
||m

||

2 mm + k m +k

F o ( k||F }
o

+k m k||
m + k m m
k
4
||
mm m
X m= 1
X||

||+k

|| m

||
||m
2

||+k

+k m k||
m +k m m
k
|| 4
mm m
k||F

||=
X

||m

Considerando que a frequncia natural do absorvedor igual a frequncia natural do

sistema, temos que:
||

m
k||
. :k
m||

||=

Substituindo a rigidez acima no clculo da amplitude do absorvedor, temos:

m||
Fo
||=

X
2

Substituindo os valores e fazendo a amplitude do absorvedor igual amplitude mxima,

obtm-se que:
m||628,318 =0,002
250

A massa do absorvedor ento igual a 0,3166 kg e a rigidez do mesmo 124988,46

N/m.
10 - Um conjunto motor-gerador, mostrado na Fig. 5.15, projetado para operar
nafaixa de 2000 a 4000 rpm. Entretanto, o conjunto vibra violentamente em 3000
rpm devido a um pequeno desbalanceamento no rotor. Prope-se colocar uma viga
engastada com uma massa na extremidade como absorvedor para eliminar o
problema. Quando a viga com uma massa de 2 kg, sintonizada a 3000 rpm, ligada
ao conjunto, as frequncias naturais do sistema obtidas foram de 2500 rpm e 3500
rpm. Projetar o absorvedor a ser ligado (especificando sua massa e sua rigidez) de
forma que as frequncias naturais do sistema fiquem fora da faixa de operao do
conjunto motor-gerador.

Figura 12. Problema a ser resolvido.

Para o limite superior:
1 =2=314,16 rad /s
Sendo:

m2

; r 1= r1 ; r 2= r 2
m1
2
2

r 1=0,833 e r 2=1,16
Logo,

( ) ( )

r 21= 1+

m2

2
=14,87 kg
1+ 1 , se r 1=0,833, =0,1345 e m1=
0,1345
2
2

Para o limite inferior:

r 1=209,44 rad / s
r 1=0,667 , para esse valor de r 1 ,m2=10,327 kg .
A segunda frequncia de ressonncia pode ser encontrada por:

( ) ( )

2
+ 1+ 1 =2,249
2
2

=4499,4
rpm. A rigidez encontrada por
Logo,
r2
2
2
k 2=2 m2=314,16 10,327=1,0188106 N / m
r 22= 1+

11 - Uma bomba centrfuga para gua de velocidade varivel, usada em

laboratrio de ensaio, vibra de forma violenta na velocidade de operao de 764
rpm que se torna impossvel realizar a leitura da presso no manmetro do tubo de
descarga. Para evitar o problema, sugere-se que se instale um absorvedor de
vibrao na carcaa da bomba, fazendo-se um furo com rosca na carcaa de modo
a fixar uma haste de ao de 7 mm de dimetro. Qual devero ser o comprimento da
haste e o peso a ser colocado em sua extremidade para que a leitura do manmetro
possa ser realizada? Como o sistema passar a ter duas freqncias naturais, ser
necessrio desmontar o absorvedor quando a bomba necessitar operar prximo a
estas freqncias. Assim, qual o valor das freqncias naturais do conjunto se a
massa da bomba de 103 kg? Dados: E = 206 GPa e = 7.800 kg/m3.
Dados:
7642
=80 rad / s
60
mb=103 kg
||=0,007 m
d
nb=

Para um bom projeto de absorvedor, a massa que deve ser colocada na extremidade do
mesmo deve ser entre 0,05 e 0,25 vezes a massa do sistema, que se deseja reduzir a
vibrao. Para este projeto utiliza-se uma massa 0,15 vezes a massa da bomba, ou seja,
||=15,45 kg
.Ento a massa equivalente do absorvedor calculada:
M

||+0,23m

0,237800 0,007
lv
v =15,45+
4
e q||=M
m

A rigidez do absorvedor tambm calculada:

3206E9 0,007
3 EI
64
e q||= 3 =
l
l 3v
k

Sabendo que a frequncia natural do absorvedor deve ser igual a frequncia natural da
bomba, tm-se que:

k eq
=80. :
me q
||

||

3206E9 0,0074
64
l 3v
2

0,237800 0,007
15,45+
lv
4

=80. :l v =0,0903m

12 - Para o sistema massa-mola de trs graus de liberdade representados pela

figura abaixo, determine:
a) Desenhe os diagramas de corpo livre;
b) Formule as equaes de movimento;
c) Defina as solues;
d) Faa m1 m2 m , m3 2m , k1 k2 k3 k e k4 k5 k6 2k e formule as
equaes das amplitudes;
e) Formule a equao caracterstica;
f) Calcule analiticamente, as freqncias naturais;
g) Substitua as freqncias naturais, (item f) nas equaes das amplitudes do
(item d), calcule os modos naturais de vibraes e desenhe as formas modais
para cada freqncia natural;
h) Formule as respostasparacadafreqncia natural, considerando as
x 0 1, 2 m x2 0 0 x&1 0 0 x&2 0 0
seguintes condies iniciais: 1
,
,
e
;
i) Trace as curvas respostas (item h) em funo dotempo usando o MATLAB.
j) Identifique a matriz massa e a matriz rigidez.

13 - Para o sistema massa-mola de trs graus de liberdade representados pela

figura abaixo, determine:
a) Desenhe os diagramas de corpo livre;
b) Formule as equaes de movimento;
c) Defina as solues;
d) Formule as equaes das amplitudes;
e) Formule a equao caracterstica;
f) Calcule analiticamente, as freqncias naturais;
g) Substitua as freqncias naturais, (item f) nas equaes das amplitudes do
(item d), calcule os modos naturais de vibraes e desenhe as formas modais
para cada freqncia natural;
h) Formule as respostasparacadafreqncia natural, considerando as
x 0 0 x2 0 0 x3 0 1 x&1 0 0
seguintes condies iniciais: 1
,
,
,
,
x&2 0 0 x&3 0 0
e
;
i) Trace as curvas respostas (item h) em funo do tempo usando o MATLAB.
j) Identifique a matriz massa e a matriz rigidez.

clearall; clc
k=50000;
m=5;
x0=0.0006;
v0=20;
K=[3.75*k -3*k;-3*k 3*k];
M=[m 0;0 m];
[AVEC,AVAL]=eig(K,M);
w1=sqrt(AVAL(1));
w2=sqrt(AVAL(4));
av1=[3*k/(-w1^2*m+3.75*k);1];
av2=[3*k/(-w2^2*m+3.75*k);1];
A2=linsolve([(av1.^2)+(w1^2).*(av1.^2),(av2.^2)+(w2^2).*(av2.^2)],
[v0^2+x0^2;v0^2+x0^2]);
A=sqrt(A2(1));
B=sqrt(A2(2));
alphasin=linsolve([A*av1,B*av2],[x0;x0]);
alpha=asin(alphasin);
t=linspace(0,0.5,2000);
x1=A.*av1(1).*sin(w1.*t+alpha(1))+B.*av2(1).*sin(w2.*t+alpha(2));
x2=A.*av1(2).*sin(w1.*t+alpha(1))+B.*av2(2).*sin(w2.*t+alpha(2));
plot(t,x1,t,x2)