Trabalho de Vibrações
Trabalho de Vibrações
Trabalho de Vibrações
1
Ministério da Educação
Universidade Federal do Paraná
Setor de Tecnologia
Departamento de Engenharia Mecânica
Curso de Graduação em Engenharia Mecânica
1 OBJETIVO
2 RESOLUÇÃO DA QUESTÃO 1
𝑘 5000
𝜔𝑛 = √ =√ = 18,26 𝑟𝑎𝑑/𝑠
𝑚 15
𝑥𝑝 (𝑡) = 𝑋 𝑠𝑖𝑛(𝜔𝑡 − ɸ)
1/2
𝑋 𝑘 2 + (𝑐𝜔)2
=[ ] ; 𝑿 = 𝟎, 𝟎𝟓𝟒𝟎𝟓 𝒎
𝑌 (𝑘 − 𝑚𝜔 2 )2 + (𝑐𝜔)2
𝑚𝑐𝜔3
ɸ = 𝑡𝑎𝑛−1 [ ] ; ɸ = 𝟎, 𝟏𝟑𝟐𝟖 𝒓𝒂𝒅
𝑘(𝑘 − 𝑚𝜔 2 ) + (𝜔𝑐)2
Tendo estes dados e a partir do sistema abaixo, foi possível obter X0 e ɸ0.
2
Ministério da Educação
Universidade Federal do Paraná
Setor de Tecnologia
Departamento de Engenharia Mecânica
Curso de Graduação em Engenharia Mecânica
Figura 1. Respostas total, particular e homogênea do sistema, assim como a comparação entre a função excitadora e a resposta.
3 RESOLUÇÃO DA QUESTÃO 2
De acordo com RAO et al [1], a força externa F(t) pode ser expandida em uma série de Fourier
quando esta é considerada periódica com um período igual a 2π/ω:
∞ ∞
𝑎0
𝐹(𝑡) = + ∑ 𝑎𝑗 cos 𝑗𝜔𝑡 + ∑ 𝑏𝑗 sin 𝑗𝜔𝑡
2
𝑗=1 𝑗=1
Onde:
2 𝜏
𝑎𝑗 = ∫ 𝐹(𝑡) cos 𝑗𝜔𝑡 𝑑𝑡 , 𝑗 = 0,1,2, …
𝜏 0
3
Ministério da Educação
Universidade Federal do Paraná
Setor de Tecnologia
Departamento de Engenharia Mecânica
Curso de Graduação em Engenharia Mecânica
2 𝜏
𝑏𝑗 = ∫ 𝐹(𝑡) sin 𝑗𝜔𝑡 𝑑𝑡 , 𝑗 = 0,1,2, …
𝜏 0
𝑏𝑗
( )
𝑘
𝑥𝑝 (𝑡) = sin(𝑗𝑤𝑡 − 𝜑𝑗 )
√(1 − 𝑗 2 𝑟 2 )2 + (2𝜁𝑗𝑟)²
Em que:
2𝜁𝑗𝑟
𝜑𝑗 = 𝑡𝑎𝑛−1 ( )
1 − 𝑗 2𝑟 2
e
𝜔
𝑟=
𝜔𝑛
∞ 𝑎𝑗 ∞
𝑏𝑗
𝑎0 ( ) ( )
𝑘 𝑘
𝑥𝑝 (𝑡) = +∑ cos(𝑗𝑤𝑡 − 𝜑𝑗 ) + ∑ sin(𝑗𝑤𝑡 − 𝜑𝑗 )
2𝑘 √(1 − 𝑗 2 𝑟 2 )2 + (2𝜁𝑗𝑟)² √(1 − 𝑗 2 𝑟 2 )2 + (2𝜁𝑗𝑟)²
𝑗=1 𝑗=1
4
Ministério da Educação
Universidade Federal do Paraná
Setor de Tecnologia
Departamento de Engenharia Mecânica
Curso de Graduação em Engenharia Mecânica
A solução da equação acima denota um sistema em regime permanente. Para obter uma
solução completa, desde as partes iniciais até as transitórias, é necessário obter as condições
de contorno x(0) e x’(0), o que nem sempre é trivial. Por meio das equações acima, é possível
estabelecer as relações de proporcionalidade da amplitude com as variáveis j e ω. São elas:
1
𝑥𝑝 ∝
𝑗
Se considerarmos que 𝑗𝜔𝑛 = 𝜔, ter-se-á uma amplitude grande, desde que os valores de j e
de 𝜁 sejam pequenos.
Solução:
A válvula é considerada como uma massa ligada a uma mola e a um amortecedor de um lado
sujeita a uma função forçante F(t) do outro lado, a qual é expressa como:
𝐹(𝑡) = 𝐴𝑝(𝑡)
𝜋(50)²
𝐴= = 625𝜋 𝑚𝑚2 = 0,000625𝜋 𝑚²
4
5
Ministério da Educação
Universidade Federal do Paraná
Setor de Tecnologia
Departamento de Engenharia Mecânica
Curso de Graduação em Engenharia Mecânica
e p(t) é a pressão que age sobre a válvula a qualquer instante t. Visto que p(t) é periódica com
período igual à 2 segundos e A é constante, F(t) também é uma função periódica com período
igual à 2 segundos. A frequência da função forçante é:
2𝜋
𝜔= = 𝜋 𝑟𝑎𝑑/𝑠
𝜏
F(t) pode ser expressa em uma série de Fourier como:
𝑎0
𝐹(𝑡) = + 𝑎1 cos 𝜔𝑡 + 𝑎2 cos 𝜔𝑡 + ⋯ + 𝑏1 sin 𝜔𝑡 + 𝑏2 sin 𝜔𝑡 + ⋯
2
onde 𝑎𝑗 e 𝑏𝑗 são fornecidos pelas equações já demonstradas. Já que a função F(t) é dada
por:
2 1 2
𝑎0 = ∫ 50000𝐴𝑡 𝑑𝑡 + ∫ 50000𝐴(2 − 𝑡)𝑑𝑡 = 50000𝐴
2 0 1
2 1 2
2 . 105 𝐴
𝑎1 = ∫ 50000𝐴𝑡 cos 𝜋𝑡 𝑑𝑡 + ∫ 50000𝐴(2 − 𝑡) cos 𝜋𝑡 𝑑𝑡 = −
2 0 1 π²
2 1 2
𝑏1 = [∫ 50000𝐴𝑡 sin 𝜋𝑡 𝑑𝑡 + ∫ 50000𝐴(2 − 𝑡) sin 𝜋𝑡 𝑑𝑡] = 0
2 0 1
2 1 2
𝑎2 = [∫ 50000𝐴𝑡 cos 2𝜋𝑡 𝑑𝑡 + ∫ 50000𝐴(2 − 𝑡) cos 2𝜋𝑡 𝑑𝑡] = 0
2 0 1
2 1 2
𝑏2 = [∫ 50000𝐴𝑡 sin 2𝜋𝑡 𝑑𝑡 + ∫ 50000𝐴(2 − 𝑡) sin 2𝜋𝑡 𝑑𝑡] = 0
2 0 1
2 1 2
2 . 105 𝐴
𝑎3 = [∫ 50000𝐴𝑡 cos 3𝜋𝑡 𝑑𝑡 + ∫ 50000𝐴(2 − 𝑡) cos 3𝜋𝑡 𝑑𝑡] = −
2 0 1 9π²
2 1 2
𝑏3 = [∫ 50000𝐴𝑡 sin 3𝜋𝑡 𝑑𝑡 + ∫ 50000𝐴(2 − 𝑡) sin 3𝜋𝑡 𝑑𝑡] = 0
2 0 1
6
Ministério da Educação
Universidade Federal do Paraná
Setor de Tecnologia
Departamento de Engenharia Mecânica
Curso de Graduação em Engenharia Mecânica
2 . 105 𝐴 2 . 105 𝐴
𝐹(𝑡) ≈ 25000𝐴 − cos 𝜔𝑡 − cos 3𝜔𝑡
π2 9π²
A resposta em regime permanente da válvula à função forçante da equação acima pode ser
expressa como:
105 𝐴 105 𝐴
25000𝐴 (2 . ) (2 . )
𝑘π2 9𝑘π2
𝑥𝑝 (𝑡) = − cos(𝜔𝑡 − 𝜑1 ) − cos(3𝜔𝑡 − 𝜑3 )
𝑘 √(1 − 𝑟 2 )2 + (2𝜁𝑟)2 √(1 − 9𝑟 2 )2 + (6𝜁𝑟)2
𝑘 3000
𝜔𝑛 = √ =√ = 100 𝑟𝑎𝑑/𝑠
𝑚 0,3
2π 2π
𝜔= = = π rad/s
𝜏 2
Desta forma, a razão de frequências pode ser calculada:
𝜔 π
𝑟= = = 0,031416
𝜔𝑛 100
E o fator de amortecimento:
𝑐 𝑐 15
𝜁= = = = 0,25
𝑐𝑐 2𝑚𝜔𝑛 2(0,30)(100)
2𝜁𝑟 2(0,25)(0,031416)
𝜑1 = 𝑡𝑎𝑛−1 ( 2
) = 𝑡𝑎𝑛−1 ( ) = 0,0157222 𝑟𝑎𝑑
1−𝑟 1 − 0,0314162
6𝜁𝑟 6(0,25)(0,031416)
𝜑3 = 𝑡𝑎𝑛−1 ( 2
) = 𝑡𝑎𝑛−1 ( ) = 0,04751056 𝑟𝑎𝑑
1 − 9𝑟 1 − 9(0,031416)2
105 𝐴 105 𝐴
25000𝐴 (2 . ) (2 . )
𝑘π2 9𝑘π2
𝑥𝑝 (𝑡) = − cos(𝜔𝑡 − 𝜑1 ) − cos(3𝜔𝑡 − 𝜑3 )
𝑘 √(1 − 𝑟 2 )2 + (2𝜁𝑟)2 √(1 − 9𝑟 2 )2 + (6𝜁𝑟)2
7
Ministério da Educação
Universidade Federal do Paraná
Setor de Tecnologia
Departamento de Engenharia Mecânica
Curso de Graduação em Engenharia Mecânica
105 (0,000625)
25000(0,000625π) (2 . )
3000π2
𝑥𝑝 (𝑡) ≈ − cos(π𝑡 − 0,0157222 )
3000 √(1 − 0,0314162 )2 + (2(0,25)(0,031416))2
105 (0,000625)
(2 . )
9(3000)π2
− cos(3π𝑡 − 0,04751056)
√(1 − 9(0,0314162 )2 + (6(0,25)(0,031416))2
4 REFERÊNCIAS BIBLIOGRÁFICAS