Zbigniew Starczewski
Drgania mechaniczne
Warszawa 2010
Politechnika Warszawska
Wydział Samochodów i Maszyn Roboczych
Kierunek "Edukacja techniczno informatyczna"
02-524 Warszawa, ul. Narbutta 84, tel (22) 849 43 07, (22) 234 83 48
ipbmvr.simr.pw.edu.pl/spin/, e-mail: sto@simr.pw.edu.pl
Opiniodawca: prof. nzw. dr hab. Zbigniew SKUP
Projekt okładki: Norbert SKUMIAA, Stefan TOMASZEK
Projekt układu graficznego tekstu: Grzegorz LINKIEWICZ
Skład tekstu: Janusz BONAROWSKI, Piotr KORCZAK-KOMOROWSKI
Publikacja przeznaczona jest dla studentów kierunku
"Edukacja techniczno informatyczna"
Copyright © 2010 Politechnika Warszawska
Utwór w całości ani we fragmentach nie może być powielany
ani rozpowszechniany za pomocą urządzeń elektronicznych, mechanicznych,
kopiujÄ…cych, nagrywajÄ…cych i innych bez pisemnej zgody posiadacza praw
autorskich.
ISBN 83-89703-45-9
Druk i oprawa: Drukarnia Expol P. Rybiński, J. Dąbek Spółka Jawna,
87-800 Włocławek, ul. Brzeska 4
Spis treści
Wstęp..................................................................... 5
1. Wprowadzenie................................................... 7
2. Kinematyka drgań .......................................... 11
2.1 Pojęcia podstawowe ..................................................................... 12
3. Składanie ruchów harmonicznych .................. 15
3.1 Składanie drgań o takich samych częstościach .......................... 16
3.2 Składanie drgań o różnych częstościach..................................... 17
4. Elementy analizy harmonicznej ...................... 31
4.1 Przekształcenie Fourier a ........................................................... 32
5. Modelowanie układów drgających.................. 39
6. Układanie równań ruchu ................................. 43
7. Siły w ruchu drgającym................................... 55
8. Krótka klasyfikacja drgań ............................... 63
9. Drgania swobodne liniowego układu
drgajÄ…cego o jednym stopniu swobody
(bez tłumienia)................................................ 67
10. Drgania swobodne układu
o jednym stopniu swobody
tłumione tarciem wiskotycznym .................... 81
11. Drgania wymuszane układu
o jednym stopniu swobody bez tłumienia... 93
12. Drgania wymuszane liniowego układu
drgajÄ…cego o jednym stopniu swobody
z tłumieniem wiskotycznym ......................... 103
13. Drgania liniowe układu o jednym stopniu
swobody wymuszane bezwładnościowo
(z tłumieniem) .............................................. 113
14. Drgania układów o jednym stopniu
swobody przy wymuszeniu kinematycznym
(z tłumieniem) .............................................. 121
15. Amortyzacja drgań ....................................... 131
16. Rejestracja drgań......................................... 137
17. Drgania swobodne układu liniowego
o dwóch stopniach swobody
bez tłumienia ............................................. 143
18. Drgania wymuszane układów o dwóch
stopniach swobody, tłumienie dynamiczne . 151
19. Literatura...................................................... 163
Wstęp
Niniejsze materiały zostały opracowane w ramach realizacji Programu
Rozwojowego Politechniki Warszawskiej współfinansowanego ze środ-
ków PROGRAMU OPERACYJNEGO KAPITAA LUDZKI. Przezna-
czone są dla studentów studiów inżynierskich kierunku Edukacja tech-
niczno-informatyczna prowadzonych na Wydziale Samochodów i Ma-
szyn Roboczych Politechniki Warszawskiej
Niniejsze opracowanie przygotowano dla przedmiotu pt. DRGANIA
MECHANICZNE . Jego zawartość merytoryczna w pełni odpowiada
zakresowi opisanemu w sylabusie opracowanym dla tego przedmiotu.
Całość opracowanych materiałów dydaktycznych dla ww przedmiotu za-
warta została w 18 rozdziałach.
Rozdział 1 został poświęcony ogólnym pojęciom z zakresu drgań, opisa-
ny jest cel badania drgań, możliwości zastosowania ruchów drgających,
a także definicji ruchu drgającego.
Rozdział 2 dostarczy kinematyki ruchów drgających, opisane są podsta-
wowe pojęcia takie jak przemieszczenia w ruchu drgającym okres, drgań
częstość drgań, faza początkowa, definicja ruchu okresowego.
Rozdział 3 poświęcony jest składaniu ruchów drgających o tych samych
częstościach i różnych amplitudach, o różnych częstościach drgań skła-
dowych, omówione zostało pojęcie dudnienia.
Rozdział 4 zawiera elementy analizy harmonicznej co związane jest
z rozwinięciem funkcji okresowej w szereg Fourier`a.
Rozdział 5 poświęcony jest problemom modelowania rzeczywistych
układów drgających. Naturalną konsekwencją procesu modelowania jest
opis matematyczny modelu. Wiąże się to z układaniem równań ruchu.
Metodom układania równań ruchu poświęcony jest rozdział 6.
W równaniach ruchu występują określone siły związane z ruchami
drgającymi. Siły te są opisane w rozdziale 7.
Rozdział 8 zawiera krótką podstawową klasyfikacje drgań. Autor
posłużył się tu podziałem zaproponowanym przez prof. Zbigniewa
Osińskiego.
Rozdziały 9, 10 poświęcone są drganiom swobodnym układów o jednym
stopniu swobody bez tłumienia i z tłumieniem.
Rozdziały 11, 12, 13, 14 poświęcone są drganiom układów o jednym
stopniu swobody z różnymi rodzajami wymuszeń (siłowym, bezwład-
nościowym i kinematycznym). Rozpatrzono pojęcie współczynnika
uwielokrotnienia amplitudy drgań oraz krzywych rezonansowych.
Rozdział 15 przedstawia problematykę amortyzacji drgań, to znaczy
ochronę otoczenia przed skutkami drgań obiektu i ochronę obiektu przed
skutkami drgań otoczenia.
Ważnemu problemowi pomiaru parametrów układów drgających (czę-
stość, amplituda, miejsca występowania) poświęcony jest rozdział 16.
Rozdziały 17 i 18 poświęcone są drganiom układów o dwóch stopniach
swobody bez tłumienia (swobodnych) oraz z wymuszeniem harmonicz-
nym. Czytelnik jest wprowadzony w pojęcie tak zwanego tłumienia
dynamicznego.
Należy podkreślić iż do każdego rozdziału wprowadzone są przykłady
zadaniowe pokazujące zastosowanie przedstawionego materiału teore-
tycznego. Przewidziano także dwa ćwiczenia laboratoryjne (krzywe re-
zonansowe belki z wymuszeniem bezwładnościowym oraz badanie
dynamicznego eliminatora drgań) które lepiej utrwalą przedstawiony
materiał teoretyczny i zadaniowy. Myślę, że tak skonstruowane mater-
iały dydaktyczne pomogą słuchaczowi w nabyciu teoretycznych i prak-
tycznych umiejętności z zakresu przedstawionego materiału.
Zajęcia dydaktyczne zdecydowanej większości przedmiotów składają-
cych się na program studiów będą realizowane, oprócz wykładu, także
w formie ćwiczeń laboratoryjnych prac projektowych. Dlatego istotną
częścią tych materiałów, oprócz prezentacji materiału teoretycznego, są
opisy przebiegu ćwiczeń wykonywanych podczas zajęć dydaktycznych
oraz propozycje zadań do samodzielnego wykonania przez słuchaczy.
Tak skonstruowane materiały dydaktyczne pomogą słuchaczom w naby-
ciu praktycznych umiejętności z zakresu posługiwania się technikami
komputerowymi niezbędnych w realizacji współczesnych procesów
projektowo wytwórczych.
Wprowadzenie
1
ROZDZIAA 1
Zachowanie się układów mechanicznych w trakcie drgań jest nieustannie
przedmiotem zainteresowań wielu badaczy i instytucji naukowych.
Chodzi o zbadanie, jaki wpływ mają drgania na wytrzymałość i żywot-
ność, a co za tym idzie niezawodność elementów i maszyn, oraz jakie są
przyczyny, zródła drgań, i jak ochronić się przed nimi.
W wyniku drgań elementów maszyn pojawiają się negatywne zjawiska
z których do najważniejszych zaliczamy:
1. Zakłócenia prawidłowości działania maszyn.
Nadmierne drgania mogą spowodować wadliwą, nierównomierną pracę
maszyn i urządzeń. Np. w obrabiarkach mogą utrudnić uzyskanie odpo-
wiedniej dokładności obróbki. w elementach złącznych gwintowych,
zaciskowych mogą być przyczyną ich rozłączania się.
2. Zmniejszenie trwałości maszyn i urządzeń.
Zjawisko drgań powoduje powstawanie w elementach maszyn zmien-
nych naprężeń co prowadzi poprzez procesy zmęczeniowe do szybszego
ich zużycia. Szczególnie grozne jest to w przypadku wałów maszy-
nowych, łożysk ślizgowych i tocznych, łopatek wirników, wszelkiego
rodzaju elementów zawieszeń.
3. Niekorzystny wpływ drgań na organizm człowieka.
Generalnie wszelkie postacie drgań mają wpływ szkodliwy dla organiz-
mu ludzkiego. Drgania powstajÄ…ce w maszynach roboczych takich jak
młoty pneumatyczne, koparki, żurawie budowlane, walcarki i szereg
innych bywają bardzo często powodami t. zw. chorób zawodowych.
4. Hałas.
Drgania są przyczyną hałasu. yródła hałasu są różnorakie, są to zarówno
drgania ośrodka (gazy), jak i drgania elementów maszyn i urządzeń.
długotrwałe przebywanie w środowisku o podwyższonym hałasie
wywołuje uczucie zmęczenia, rozdrażnienia, występuje zjawisko stresu,
a często uszkodzenie organów człowieka (głuchota) lub w przypadku
infradzwięków (drgania o bardzo niskich częstotliwościach) wręcz
fizyczne nieodwracalne uszkodzenie tych organów.
Strona 8
8
8
8
WPROWADZENIE
Generalnie należy mówić o szkodliwości drgań, jednakże bywają one
wykorzystywane z pożytkiem dla człowieka. Występuje to w przypadku
wszelkiego rodzaju przenośników wibracyjnych, przesiewaczy, zagęsz-
czaczy.
Cały oddzielny rozdział to muzyka. Zarówno ta poważna, jak i rozryw-
kowa. Któż z nas nie podziwiał wspaniałych solówek wykonanych na
instrumentach dętych, strunowych czy perkusyjnych.
Zdefiniujmy zatem co to jest drganie zwane niekiedy ruchem drgajÄ…cym.
Według Osińskiego definicja ta ma następujące brzmienie: DRGANIEM
lub RUCHEM DRGAJCYM nazywamy taki ruch w którym badana
współrzędna na przemian zbliża się i oddala od pewnej wartości prze-
ciętnej. Wartość ta może być ustalona w czasie. Zwykle przyjmuje się ją
zerową w przyjętym układzie współrzędnych. Wartość przeciętna może
też być zmienna w czasie w dowolny sposób.
a) b)
Rysunek 1.1 Ruchy drgające. a) wartość przeciętna równa zero,
b) wartość przeciętna zmienna w czasie
Strona 9
9
9
9
ROZDZIAA 1
Strona 10
10
10
10
`
Kinematyka drgań
2
ROZDZIAA 2
Pojęcia podstawowe
Podstawowe pojęcia związane z drganiami oparte są na opisie ruchu
harmonicznego prostego. Ruch taki opisany jest równaniem.
x = a cos(Ét + Õ0 ) (2.1)
gdzie:
x współrzędna ruchu drgającego,
a amplituda drgań,
É czÄ™stość kÄ…towa drgaÅ„,
Õ0 faza poczÄ…tkowa drgaÅ„ (przesuniÄ™cie fazowe),
t czas,
Rysunek 2.1 Ilustracja przebiegu drgań w ruchu harmonicznym prostym
i podstawowe parametry tego ruchu
x0 amplituda poczÄ…tkowa,
T okres drgań.
Pomiędzy częstością f wyrażoną w hercach (wielkość ta zwana jest
przez elektrotechników częstotliwością) i okresem T zachodzą
zależności:
2Ä„ 2Ä„ 1 É
T = ; É = = 2Ä„f ; f = = (2.2)
É T T 2Ä„
Strona 12
12
12
12
KINEMATYKA DRGAC
Zapis ruchu harmonicznego może być przedstawiony w postaci:
x = AcosÉt + B sinÉt (2.3)
Jeżeli dokonamy podstawienia:
A = a cosÈ ; B = -a sinÈ (2.4)
x = a cosÈ cosÉt - asinÈ sinÉt (2.5)
cosÉt cosÈ - sinÉt sinÈ = cos(Ét +È )
Ostatecznie:
x = a cos(Ét +È ) (2.6)
gdzie w zależności (2.4)
B
B
öÅ‚
a = A2 + B2 ; tgÈ = - ; È = arctgëÅ‚- ÷Å‚
(2.7)
ìÅ‚
A
A
íÅ‚ Å‚Å‚
Ruch będziemy nazywać okresowym wtedy, gdy spełniona będzie
zależność:
x(t + T ) = x(t) (2.8)
Należy pamiętać iż każdy ruch harmoniczny jest ruchem okresowym,
natomiast nie każdy ruch okresowy jest ruchem harmonicznym.
Strona 13
13
13
13
ROZDZIAA 2
Strona 14
14
14
14
`
Składanie
ruchów harmonicznych
3
ROZDZIAA 3
3.1. Składanie drgań o takich samych częstościach
Drganie wypadkowe x(t) jest sumą dwóch drgań harmonicznych o tej
samej czÄ™stoÅ›ci koÅ‚owej É i różnych amplitudach a i różnych przesuniÄ™-
ciach fazowych Õ .
x(t) = a1 sin(Ét + Õ1) + a2 sin(Ét + Õ2 )
sin(Ét + Õ1) = sinÉt cosÕ1 + cosÉt sinÕ1
Å„Å‚
(3.1.1)
òÅ‚
ółsin(Ét + Õ2 ) = sinÉt cosÕ2 + cosÉt sinÕ2
x(t) = a1 sinÉt cosÕ1 + a1 cosÉt sinÕ1 +
(3.1.2)
a2 sinÉt cosÕ2 + a2 cosÉt sinÕ2
GrupujÄ…c wyrazy z cosÉt i sin Ét otrzymujemy:
x(t) = (a1 sinÕ1 + a2 sinÕ2 ) cosÉt
(3.1.3)
+ (a1 cosÕ1 + a2 cosÕ2 )sinÉt
Wykonujemy podstawienie:
a1 sinÕ1 + a2 sinÕ2 = a sinÈ
Å„Å‚
(3.1.4)
òÅ‚
cosÕ1 + a2 cosÕ2 = a cosÈ
óła1
x(t) = a cosÉt sinÈ + a sinÉt cosÈ (3.1.5)
Ostatecznie:
x(t) = asin(Ét +È ) (3.1.6)
gdzie z (2.1.2)
2 2
a = (a1 sinÕ1) + (a2 sin Õ2 ) (3.1.7)
ëÅ‚ a1 sinÕ1 + a2 sinÕ2 öÅ‚
ìÅ‚ ÷Å‚
È = arctgìÅ‚ (3.1.8)
a1 cosÕ1 + a2 cosÕ2 ÷Å‚
íÅ‚ Å‚Å‚
Strona 16
16
16
16
SKAADANIE RUCHÓW HARMONICZNYCH
Wynika stąd iż drganie wypadkowe będzie też drganiem harmonicznym
o czÄ™stoÅ›ci É . Przypadek ten można uogólnić na sumÄ™ n drgaÅ„ harmo-
nicznych o czÄ™stoÅ›ci É .
i=n
x(t) = sin(Ét + Õi ) = a sin(Ét +È ) (3.1.9)
"ai
i=1
gdzie:
i=n i=n
a = (3.1.10)
i i
"(a sin Õi )2 + "(a cosÕi )2
i=1 i=1
i=n
ëÅ‚ öÅ‚
ìÅ‚ sinÕi ÷Å‚
"ai
÷Å‚
È = arctgìÅ‚ i=1 (3.1.11)
i=n
ìÅ‚ ÷Å‚
cosÕi ÷Å‚
ìÅ‚
"ai
íÅ‚ i=1 Å‚Å‚
3.2. Składanie ruchów drgań o różnych częstościach
Rozpatrzmy przypadek, gdy wypadkowe drganie x(t) jest sumą dwóch
drgań harmonicznych o różnych częstościach.
x(t) = a1 sin(É1t + Õ1) + a2 sin(É2t + Õ2 ) (3.2.1)
Dla tej postaci równania możemy wyróżnić trzy przypadki:
a) Częstość jednego z drgań składowych jest dużo większa od częstości
drgań drugiego.
x1(t) = a1 sin(É1t + Õ1) ; x2 (t) = a2 sin(É2t + Õ2 )
załóżmy że a1 < a2 ; É1 << É2 ; Õ1 = Õ2 = 0
Strona 17
17
17
17
ROZDZIAA 3
Rysunek 3.2.1 Przypadek (a), a1 < a2 , É1 << É2 , Õ1 = Õ2 = 0
Strona 18
18
18
18
SKAADANIE RUCHÓW HARMONICZNYCH
gdy a1 > a2 , É1 << É2 otrzymujemy:
Rysunek 3.2.2 Przypadek (a), É1 << É2 , a1 > a2
Strona 19
19
19
19
ROZDZIAA 3
b) Częstości drgań składowych różnią się nieznacznie od siebie.
x1(t) = a1 sin(Ét + Õ) , x2 (t) = a2 sin(Ét + "Ét) (3.2.2)
Należy pamiętać iż:
"É << É , Õ - przesuniÄ™cie poczÄ…tkowe x1(t) wzglÄ™dem x2 (t) .
Drgania wypadkowe otrzymujemy w postaci:
x(t) = x1(t) + x2 (t) = a1 sin(Ét + Õ) + a2 sin(Ét + "Ét) =
a1 sinÉt cosÕ + a1 cosÉt sinÕ + a2 sinÉt cos "Ét +
(3.2.3)
+ a2 cosÉt sin "Ét = (a1 sinÕ + a2 sin "Ét) cosÉt +
+ (a1 cosÕ + a2 cos "Ét)sinÉt
WprowadzajÄ…c oznaczenia:
a1 sinÕ + a2 sin "Ét = A(t)sin ¨(t)
(3.2.4)
a1 cosÕ + a2 cos "Ét = A(t)cos ¨(t)
Otrzymujemy ostatecznie:
x(t) = A(t)sin(Ét + ¨(t)) (3.2.5)
gdzie:
A(t) = (a1 sinÕ + a2 sin "Ét)2 + (a1 cosÕ + a2 cos "Ét)2
a1 sinÕ + a2 sin "Ét
¨(t) = arctg
a1 cosÕ + a2 cos "Ét
Tak więc amplituda zmienia się okresowo od Amax do Amin , gdzie:
Amax = a1 + a2
(3.2.6)
Amin = a1 - a2
Drgania mają postać jak na rysunku 3.2.3 i tę postać drgań nazywamy
DUDNIENIEM.
Strona 20
20
20
20
SKAADANIE RUCHÓW HARMONICZNYCH
Rysunek 3.2.3. Przypadek (b), dudnienie
c) Stosunek częstości drgań składowych wyraża się przez niewielkie
liczby naturalne.
W tym przypadku przebiegi drgania wypadkowego x(t) w zależności od
stosunku amplitud, częstości oraz kątów przesunięcia fazowego może
przyjąć różne formy.
a1
Niech x1(t) = a1 sinÉt , x2 (t) = a2 sin(2Ét + Õ) oraz = 2
a2
Rysunek 3.2.4. Przebieg x1(t) + x2 (t) dla Õ = 0o
Strona 21
21
21
21
ROZDZIAA 3
Rysunek 3.2.5 Przebieg x1(t) + x2 (t) dla Õ = 90o
W przypadku gdy É1 i É2 nie sÄ… współmierne, to drganie wypadkowe
jest nieokresowe.
Przykład 1.
Ruch punktu opisany jest superpozycją dwóch ruchów opisanych
równaniami:
Ä„
öÅ‚
x1 = 5sinÉt , x2 = 3sinëÅ‚Ét - ÷Å‚
ìÅ‚
4
íÅ‚ Å‚Å‚
Znalezć amplitudę i przesunięcie fazowe ruchu wypadkowego.
Ä„
öÅ‚
x1 = 5sin Ét Ò! x1 = 5cosëÅ‚Ét - ÷Å‚
ìÅ‚
2
íÅ‚ Å‚Å‚
Amplituda A :
i=2 i=2
A = sinÕi )2 + cosÕi )2 =
"(ai "(ai
i=1 i=1
2 2
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
ìÅ‚5 Å" (-1) + 3Å" 2 ÷Å‚ ìÅ‚5 Å" 0 + 3Å" 2 ÷Å‚
+ E" 7,43
ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚
2 2
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
Strona 22
22
22
22
SKAADANIE RUCHÓW HARMONICZNYCH
i=2
sinÕi 5Å" (-1) - 3Å" 2
"ai
i=1 2
tgÈ = = = -3,35
i=2
cosÕi 5 Å" 0 + 3Å" 2
"ai
i=1 2
Przykład 2.
Ruch punktu opisany jest równaniami :
x1 = 5cosÉt ; x2 = 3 - 5sinÉt
Znalezć tor punktu.
Równanie parametryczne na x1 i x2 możemy przedstawić w postaci:
x1 x2 - 3
= cosÉt ; = -sinÉt
5 5
PodnoszÄ…c obustronnie do kwadratu i dodajÄ…c do siebie stronami
otrzymujemy:
x12 (x2 - 3)2
+ = cos2 Ét + sin2 Ét
25 25
Czyli ostatecznie:
x12 (x2 - 3)2
+ = 1
25 25
Rysunek 3.2.6 Trajektoria punktu którego ruch opisują równania x1i x2
Strona 23
23
23
23
ROZDZIAA 3
Jest to równanie okręgu o promieniu r = 5 i środku przesuniętym o 3
wzdłuż osi x2 .
Przykład 3.
Znalezć tor punktu opisanego równaniami:
Ä„ Ä„
öÅ‚ öÅ‚
x = 4cosëÅ‚Ét + ; y = 7 cosëÅ‚Ét - ÷Å‚
ìÅ‚ ÷Å‚ ìÅ‚
2 2
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
Ä„ Ä„ Ä„
öÅ‚
cosëÅ‚Ét + = cosÉt cos - sinÉt sin = -sinÉt
ìÅ‚ ÷Å‚
2 2 2
íÅ‚ Å‚Å‚
Ä„ Ä„ Ä„
öÅ‚
cosëÅ‚Ét - ÷Å‚
= cosÉt cos + sinÉt sin = sinÉt
ìÅ‚
2 2 2
íÅ‚ Å‚Å‚
Zatem:
x = -4sinÉt ; y = 7sinÉt
StÄ…d:
7
y = - x
4
Jest to równanie prostej przedstawionej na rysunku 3.2.7.
Rysunek 3.2.7 Prosta opisana równaniami x i y
Przykład 4.
Znalezć tor punktu opisany równaniami:
x = 3 + 5sinÉt ; y = 7 + 3cosÉt
Strona 24
24
24
24
SKAADANIE RUCHÓW HARMONICZNYCH
Równania możemy łatwo przekształcić w postać:
x - 3 y - 7
= sinÉt ; = cosÉt
5 3
PodnoszÄ…c obustronnie do kwadratu i dodajÄ…c stronami, otrzymujemy:
(x - 3)2 (y - 7)2
+ = 1
25 9
Jest to równanie elipsy o środku określonym współrzędnymi x = 3,
y = 7 i ramionach 3 i 5, przedstawionej na rysunku 3.2.8.
Rysunek 3.2.8 Elipsa opisana równaniami x i y
Przykład 5.
Znalezć tor punktu poruszającego się zgodnie z równaniami:
x = a cos(Ét + Ä… )
y = bsin(Ét + ² )
Równanie drugie możemy przedstawić w formie:
y = bsin(Ét + ² + Ä… - Ä… )
Zauważmy iż w wyrażeniu w nawiasach dodaliśmy i odjęliśmy to samo
wyrażenie ą , zatem:
y = bsin[(Ét + Ä…) + (² - Ä… )] = b[sin(Ét + Ä…) cos(² - Ä… ) +
+ cos(Ét + Ä…)sin(² -Ä…)]
Strona 25
25
25
25
ROZDZIAA 3
Z równania pierwszego otrzymujemy:
x
= cos(Ét + Ä… )
a
Zatem:
2
x
ëÅ‚ öÅ‚
sin(Ét + Ä…) = 1- ìÅ‚ ÷Å‚
a
íÅ‚ Å‚Å‚
Uwzględniając ostatecznie wyrażenie w równaniu na y mamy:
2
îÅ‚ Å‚Å‚
x x
ëÅ‚ öÅ‚
ïÅ‚ śł
y = b 1- ìÅ‚ ÷Å‚
cos(² -Ä…) + sin(² - Ä… )
a a
ïÅ‚ íÅ‚ Å‚Å‚ śł
ðÅ‚ ûÅ‚
StÄ…d:
2
y x x
ëÅ‚ öÅ‚
- sin(² -Ä…) = 1- ìÅ‚ ÷Å‚
cos(² - Ä…)
b a a
íÅ‚ Å‚Å‚
PodnoszÄ…c obustronnie do kwadratu otrzymujemy:
2 2
y 2xy x
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
2
ìÅ‚ ÷Å‚ - sin(² - Ä…) + sin (² - Ä… ) =
ìÅ‚ ÷Å‚
b ab a
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
2
ëÅ‚ öÅ‚
ìÅ‚1 - ëÅ‚ x öÅ‚ ÷Å‚
cos2 (² - Ä…)
ìÅ‚ ÷Å‚
ìÅ‚ ÷Å‚
a
íÅ‚ Å‚Å‚
íÅ‚ Å‚Å‚
Po uproszczeniu:
2 2
y xy x
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
ìÅ‚ ÷Å‚ - 2 sin(² -Ä…) + [sin2 (² -Ä…) + cos2 (² -Ä…)]=
ìÅ‚ ÷Å‚
b ab a
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
cos2 (² -Ä…)
Uwzględniając jedynkę trygonometryczną otrzymujemy ostatecznie:
2 2
y xy x
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
ìÅ‚ ÷Å‚ - 2 sin(² - Ä… ) + = cos2 (² -Ä…)
ìÅ‚ ÷Å‚
b ab a
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
Strona 26
26
26
26
SKAADANIE RUCHÓW HARMONICZNYCH
Jak widać charakter trajektorii bÄ™dzie zależaÅ‚ od wyrażenia (² -Ä… ) .
a) gdy: ² - Ä… = k Å"Ä„ ; k = 0,2,4,.....
sin(² - Ä… ) = 0 ; cos(² - Ä… ) = 1
Dla takich wartości funkcji równanie trajektorii ma postać:
2 2
y x
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
+ = 1
ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚
b a
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
Jest to równanie elipsy o środku umieszczonym w punkcie (0;0)
przyjętego kartezjańskiego układu współrzędnych,
Ä„
b) gdy: ² -Ä… = + k Å"Ä„ ; k = 0,2,4,.....
2
sin(² - Ä…) = 1; cos(² -Ä… ) = 0
Równanie trajektorii przyjmuje postać:
2 2 2
y xy x y x
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
= 0
ìÅ‚ ÷Å‚ - 2 + = 0 , czyli - ÷Å‚
ìÅ‚ ÷Å‚ ìÅ‚
b ab a b a
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
Ostatecznie:
b
y = x
a
b
Jest to równanie prostej o współczynniku kierunkowym i
a
przechodzącej przez środek układu współrzędnych,
Ä„
c) gdy: ² -Ä… = + k Å"Ä„ ; k = 1,3,5,.....
2
sin(² -Ä… ) = -1; cos(² -Ä… ) = 0
Strona 27
27
27
27
ROZDZIAA 3
a)
b)
c)
Rysunek 3.2.9 Trajektorie punktu opisane równaniem końcowym
z przykładu 5
Równanie trajektorii przyjmuje postać:
2 2 2
y xy x y x
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
+ 2 + = 0 , czyli - ÷Å‚
= 0
ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚ ìÅ‚
b ab a b a
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
ostatecznie:
Strona 28
28
28
28
SKAADANIE RUCHÓW HARMONICZNYCH
b
y = - x
a
b
Jest to równanie prostej o współczynniku kierunkowym - i
a
przechodzącej przez środek układu współrzędnych.
Przypadki a), b), c) obrazuje rysunek 3.2.9.
Przykład 6.
Ruch drgający punktu jest wypadkową następujących składowych:
x1 = 4sin(15t) ; x2 = 4sin(16t)
Ponieważ mamy do czynienia niewielką różnicą prędkości kątowych na
pewno wystąpi zjawisko dudnienia. Należy wyznaczyć maksymalne i
minimalne wartości amplitud, częstość oraz okres dudnień.
amax = a1 + a2 = 4 + 4 = 8
amin = a1 - a2 = 4 - 4 = 0
Częstość dudnień:
1
Éd = É2 - É1 =16 -15 =1îÅ‚ Å‚Å‚
ïÅ‚ śł
s
ðÅ‚ ûÅ‚
Okres dudnień:
2Ä„ 2Ä„
Td = = = 2Ä„[s]
Éd 1
Przebieg wypadkowy ilustruje rysunek 3.2.10.
Strona 29
29
29
29
ROZDZIAA 3
Rysunek 3.2.10 Wypadkowa trajektoria punku z przykładu 6
Strona 30
30
30
30
`
Elementy analizy
harmonicznej
4
ROZDZIAA 4
4.1. Przekształcenie Fourier a
Dowolny przebieg drgań okresowych można rozłożyć na sumę składo-
wych harmonicznych. Analiza harmoniczna polega na rozwinięciu
funkcji x(t) o okresie T w tak zwany szereg Fourier`a.
Szereg ten możemy przedstawić w następującej formie:
a0 i="
x(t) = + cos nÉt + bn sin nÉt) (4.1)
"(an
2
i=1
Przy czym:
2Ä„
É = (4.2)
T
Poszczególne współczynniki szeregu Fourier`a wyrażone są następujący-
mi zależnościami:
T
2
a0 = x(t)dt (4.3)
+"
T
0
T
2
an = x(t) cos nÉtdt , n = 0,1,2,3...... (4.4)
+"
T
0
T
2
bn = x(t)sin nÉtdt , n = 1,2,3...... (4.5)
+"
T
0
Kolejny wyraz szeregu Fourier`a jest nazywany n tÄ… harmonicznÄ…
drgań okresowych. Wyraz wolny a0 nazywamy składową stałą drgań.
Pierwsza harmoniczna nazywana jest harmonicznÄ… podstawowÄ….
Dla scharakteryzowania składowych harmonicznych drgań okresowych
stosuje się tak zwane widmo funkcji będące zbiorem par liczb, a miano-
wicie kolejnych czÄ™stoÅ›ci Én oraz sumy kwadratów odpowiadajÄ…cych
im amplitud.
an 2 + bn 2 = An 2
Strona 32
32
32
32
ELEMENTY ANALIZY HARMONICZNEJ
Rysunek 4.1.1 Widmo funkcji
Najczęściej zdarza się, że harmoniczne wyższego rzędu mają małe
amplitudy, wtedy można przybliżoną funkcję wyrazić przez tak zwany
wielomian Fouriera`a.
N
x(t) = cos nÉt + bn sin nÉt) (4.6)
"(an
1
Oczywiste jest, że im więcej harmonik uwzględniamy tym bardziej
dokładnie rozwijana funkcja w szereg Fouriera`a oddaje oryginał.
Przykład 1.
Znalezć widmo funkcji przedstawionej na rysunku 4.1.2.
Rysunek 4.1.2 Przebieg badanej funkcji
Aatwo zauważyć iż okres funkcji wynosi 2Ą , a częstość 1. Poszczególne
współczynniki szeregu Fouriera`a obliczamy z zależności (4.3), (4.4),
(4.5).
Strona 33
33
33
33
ROZDZIAA 4
T Ä„ 2Ä„
2 2 2 Ä„
a0 = x(t)dt =
+" +"Ä„dt + +"0dt = t = Ä„
0
T 2Ä„ 2Ä„
0 0 Ä„
T Ä„ 2Ä„
2 2 2
an = x(t) cos nÉtdt = cos nÉtdt + cos nÉtdt =
+" +"Ä„ +"0
Ä„ 2Ä„ 2Ä„
0 0 Ä„
Ä„ Ä„ Ä„ Ä„
2Ä„ 1
+"cos nÉtdt = +"cos n tdt = +"cos ntdt = - sin nt = 0
2Ä„ n
0
0 0 0
T Ä„ 2Ä„
2 2 2
bn = x(t) sin nÉtdt = sin nÉtdt +
+" +"Ä„ +"0sin nÉtdt =
Ä„ 2Ä„ 2Ä„
0 0 Ä„
Ä„ Ä„ Ä„ Ä„
2Ä„ 1 1 1
+"sin nÉtdt = +"sin n tdt = +"sin ntdt = cos nt = - cos nÄ„ + =
2Ä„ n n n
0
0 0 0
1
(1 - cosĄ )
n
2
Gdy n parzyste, wtedy bn = 0 , gdy n nieparzyste, wtedy bn = .
n
Zatem poszczególne współczynniki An 2 będą wynosić:
2
A0 2 = a0 2 = Ä„
A12 = b12 = 4
A3 2 = b3 2 E" 0,44
Rysunek 4.1.3 Widmo funkcji wyznaczone w przykładzie 1
Strona 34
34
34
34
ELEMENTY ANALIZY HARMONICZNEJ
Przykład 2.
Wyznaczyć współczynniki widma funkcji określonej zależnością:
Ä„
ëÅ‚Ét öÅ‚
4
x(t) = sin +
ìÅ‚ ÷Å‚
3
íÅ‚ Å‚Å‚
Powyższe współrzędne możemy podać jako iloczyn:
Ä„ Ä„
öÅ‚ öÅ‚
x(t) = sin2 ëÅ‚Ét + Å"sin2 ëÅ‚Ét +
ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚
3 3
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
Ale:
Ä„
öÅ‚
1- cos 2ëÅ‚Ét +
ìÅ‚ ÷Å‚
Ä„ 3
öÅ‚
íÅ‚ Å‚Å‚
sin2 ëÅ‚Ét + =
ìÅ‚ ÷Å‚
3 2
íÅ‚ Å‚Å‚
Zatem:
ëÅ‚ Ä„ öÅ‚ ëÅ‚ Ä„ öÅ‚
öÅ‚÷Å‚ öÅ‚÷Å‚
ìÅ‚1- cos 2ëÅ‚Ét + ìÅ‚1- cos 2ëÅ‚Ét +
ìÅ‚ ÷Å‚÷Å‚ ìÅ‚ ìÅ‚ ÷Å‚÷Å‚
ìÅ‚
3 3
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
x(t) = Å" =
2 2
îÅ‚ Å‚Å‚
1 îÅ‚ Ä„ Å‚Å‚ îÅ‚ Ä„ Å‚Å‚
öÅ‚ öÅ‚
= cos 2ëÅ‚Ét + Å" cos 2ëÅ‚Ét + =
ïÅ‚
ïÅ‚1- ìÅ‚ 3 ÷łśł ïÅ‚1- ìÅ‚ 3 ÷łśłśł
4
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚
ðÅ‚ ûÅ‚
1 ëÅ‚ Ä„ Ä„ öÅ‚
öÅ‚ öÅ‚÷Å‚
= ìÅ‚1- cos 2ëÅ‚Ét + + cos2 2ëÅ‚Ét + =
ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚÷Å‚
ìÅ‚
4 3 3
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
íÅ‚ Å‚Å‚
Strona 35
35
35
35
ROZDZIAA 4
Ä„
öÅ‚
1+ cos 4ëÅ‚Ét +
ìÅ‚ ÷Å‚
Ä„ 3
öÅ‚
íÅ‚ Å‚Å‚
= cos2 2ëÅ‚Ét + = =
ìÅ‚ ÷Å‚
3 2
íÅ‚ Å‚Å‚
1 1 Ä„ 1 îÅ‚ Ä„ Å‚Å‚
öÅ‚ öÅ‚
= - cos 2ëÅ‚Ét + + cos 4ëÅ‚Ét + =
ìÅ‚ ÷Å‚
ïÅ‚1+ ìÅ‚ 3 ÷łśł
4 2 3 8
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
ðÅ‚ ûÅ‚
1 1 2 1 1 4
öÅ‚ öÅ‚
= - cosëÅ‚2Ét + Ä„ + + cosëÅ‚4Ét + Ä„ =
ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚
4 2 3 8 8 3
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
3 1 2 1 4
öÅ‚ öÅ‚
= - cosëÅ‚2Ét + Ä„ + cosëÅ‚4Ét + Ä„
ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚
8 2 3 8 3
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
Ale:
2 2 2
öÅ‚
cosëÅ‚2Ét + Ä„ = cos 2Ét cos Ä„ - sin 2Ét sin Ä„
ìÅ‚ ÷Å‚
3 3 3
íÅ‚ Å‚Å‚
4 4 4
öÅ‚
cosëÅ‚4Ét + Ä„ = cos 4Ét cos Ä„ - sin 4Ét sin Ä„
ìÅ‚ ÷Å‚
3 3 3
íÅ‚ Å‚Å‚
Tak więc:
3 1 2 2
îÅ‚cos Å‚Å‚
x(t) = - 2Ét cos Ä„ - sin 2Ét sin Ä„ +
ïÅ‚ śł
8 2 3 3
ðÅ‚ ûÅ‚
1 4 4
îÅ‚ Å‚Å‚
+ 4Ét cos Ä„ - sin 4Ét sin Ä„ =
ïÅ‚cos śł
8 3 3
ðÅ‚ ûÅ‚
îÅ‚ Å‚Å‚
3 1 1 3
= - cos 2Ét - sin 2Étśł +
ïÅ‚-
8 2 2 2
ðÅ‚ ûÅ‚
îÅ‚ Å‚Å‚
1 1 3
+ cos 4Ét + sin 4Étśł =
ïÅ‚-
8 2 2
ðÅ‚ ûÅ‚
3 1 3 1 3
= + cos 2Ét + sin 2Ét - cos 4Ét + sin 4Ét
8 4 4 16 16
Z powyższej analizy:
a0 3 3
= Ò! a0 =
2 8 4
Strona 36
36
36
36
ELEMENTY ANALIZY HARMONICZNEJ
a1 = 0 ; b1 = 0
1 3
a2 = ; b2 =
4 4
a3 = 0 ; b3 = 0
1 3
a4 = - ; b4 =
16 16
Rysunek 4.1.4 Widmo funkcji analizowanej w przykładzie 2
Strona 37
37
37
37
ROZDZIAA 4
Strona 38
38
38
38
`
Modelowanie układów
drgajÄ…cych
5
ROZDZIAA 5
Gdy przystępujemy do analizy drgań konkretnego układu musimy przed-
stawić układ rzeczywisty w postaci modelu o mniejszym lub większym
stopniu komplikacji. W skład takiego modelu wchodzą punkty ma-
terialne, ciała sztywne, ciała odkształcalne o masach różnych od zera,
ciała odkształcalne o masach przyjmowanych jako zerowe. Proces
modelowania polega na wprowadzaniu pewnych uproszczeń w stosunku
do rzeczywistego układu drgającego. Gdy do analizy przyjęlibyśmy
układ rzeczywisty okazałoby się iż bardzo skomplikowana (niekiedy
wręcz niemożliwa) analiza dawała by niewiele lepsze rezultaty niż jej
uproszczony model.
Położenie modelu określa się współrzędnymi uogólnionymi. Jeżeli anali-
zowany układ składa się ze skończonej liczby punktów materialnych lub
ciał sztywnych, to liczba współrzędnych uogólnionych jest skończona.
Gdy mamy do czynienia z ciałami odkształcalnymi o masach rozłożo-
nych w sposób ciągły, wtedy liczba współrzędnych uogólnionych jest
nieskończenie wielka i mówimy że analizowany układ ma nieskończenie
wielką liczbę stopni swobody. W naszych rozważaniach będziemy zaj-
mować się układami o skończonej liczbie stopni swobody. Szczególnym
przypadkiem tych układów jest układ o jednym stopniu swobody.
Ewolucje procesu modelowania prześledzimy na prostym przykładzie
n.p. samochodu.
Rysunek 5.1 Model pojazdu jako układ o jednym stopniu swobody
Strona 40
40
40
40
MODELOWANIE UKAADÓW DRGAJCYCH
Pojazd przedstawiony został jako ciało o masie m z elementami zawie-
szenia (sztywność k , oraz tłumik o współczynniku tłumienia c ),
wykonujące drgania pionowe wywołane oddziaływaniem funkcji opisa-
nej drogą f (s) . Aatwo zauważyć iż ten najprostszy model ma jeden sto-
pień swobody, a ruch ciała o masie m , określany jest jedną współrzędną
x , będącą jego przemieszczeniem.
Nawet kompletna noga techniczna zauważy iż model przedstawiony
na rysunku 5.1 ma mało wspólnego z rzeczywistym pojazdem. Spróbuj-
my zatem nieco skomplikować badany układ samochodu.
Rysunek 5.2 Model pojazdu jako układ o dwóch stopniach swobody
Widać iż przedstawiony na rysunku 5.2 model pojazdu trochę zbliżył się
do rzeczywistości. Ma on teraz dwa stopnie swobody, resorowana masa
pojazdu może poruszać się niezależnie pionowo i jednocześnie wykony-
wać ruch obrotowy wokół środka masy z . Tak więc jego ruch opisany
jest dwoma współrzÄ™dnymi, przemieszczeniem x i kÄ…tem obrotu Õ . Ko-
lejny etap przybliżania modelu do rzeczywistości obrazuje rysunek 5.3.
W modelu tym uwzględniono sztywności kop , tłumienie cop opon pojaz-
du i masy kół mk .
Nastąpiło dalsze powiększenie liczby stopni swobody modelu. Zwięk-
szyła się liczba współrzędnych opisujących ruch masy resorowanej.
Zwiększyła się zatem liczba równań opisujących ten model.
Strona 41
41
41
41
ROZDZIAA 5
Rysunek 5.3 Model pojazdu w którym uwzględniono masę kół, oraz
współczynniki sztywności i tłumienia opon
Proces modelowania powinniśmy zakończyć na tym stopniu komplika-
cji, który pozwoli poznać najwłaściwsze parametry drganiowe analizo-
wanego układu rzeczywistego.
Strona 42
42
42
42
`
Układanie równań
ruchu
6
ROZDZIAA 6
Mając stworzony mniej lub bardziej przybliżony model rzeczywistego
układu drgającego możemy pokusić się o jego opis matematyczny co
pozwoli na dalszą analizę jego parametrów drganiowych. Równania
ruchu stworzonego układu materialnego można wyprowadzić za pomocą
dowolnej z metod poznanej z wykładu z mechaniki. w szczególnie
w prostych przypadkach gdy w grę wchodzą układy o jednym stopniu
swobody można zastosować bezpośrednio II zasadę dynamiki Newtona
lub metodą energetyczną. Dla układów bardziej złożonych wygodniej
posłużyć się równaniami Lagrange`a drugiego rodzaju. Przybliżmy te
trzy metody.
a) Metoda Newtona
Rozpatrzmy układ o jednym stopniu swobody składający się z ciała
o masie m mogącego poruszać się pionowo, pobudzanego do drgań siłą
P , podpartego elementami sztywnymi o sztywności k , i elementami
tłumiącymi o współczynniku tłumienia c . Układ przedstawiony jest na
rysunku 6.1.
Rysunek 6.1 Rozpatrywany model i układ sił działających na ciało
o masie m
Zgodnie z przyjętym układem współrzędnych równanie równowagi sił
będzie następujące:
P + G + B - S - R = 0 (6.1)
Strona 44
44
44
44
UKAADANIE RÓWNAC RUCHU
gdzie:
P - siła zewnętrzna,
B - siła bezwładności,
G - obciążenie układu,
S - siła indukowana w elemencie sprężystym,
R - siła oporu.
Podstawiając pod poszczególne oznaczenia konkretne zależności, rów-
nanie (6.1), stanie się równaniem ruchu analizowanego układu.
b) Metoda energetyczna
Dla układów zachowawczych, to znaczy takich, w których całkowita
energia układu pozostaje niezmienna w czasie ruchu, możemy układać
równania ruchu w oparciu o zasadę:
EK + EP = const (6.2)
gdzie:
EK - energia kinetyczna,
EP - energia potencjalna.
Różniczkując po czasie zależność (6.2) otrzymujemy:
d
(EK + EP ) = 0 (6.3)
dt
W zastosowaniu do układów drgających, równanie (6.3) staje się równa-
niem ruchu.
c) Metoda z zastosowaniem równania Lagrange`a II rodzaju
Równanie Lagrange`a dla układów holonomicznych ma postać:
ëÅ‚ öÅ‚
d "EK "EP "EK
ìÅ‚ ÷Å‚
+ - = Qj - R (6.4)
j
ìÅ‚ ÷Å‚
&
dt "q "q "q
j j j
íÅ‚ Å‚Å‚
Strona 45
45
45
45
ROZDZIAA 6
gdzie:
EK - energia kinetyczna układu drgającego,
EP - energia potencjalna układu drgającego,
Q = Qj (t) - zewnętrzna siła uogólniona odpowiadająca
j
współrzędnej uogólnionej q .
j
R - uogólniona siła oporu, odpowiadająca współrzędnej
j
uogólnionej q skierowana przeciwnie doQ .
j j
Po wykonaniu operacji różniczkowania zgodnie z zależnością (6.4) uzys-
kujemy bezpośrednio równania ruchu układu drgającego.
Przykład 1.
Wyznaczyć sztywność zastępczą sprężyn oraz wyznaczyć równania
ruchu układu jak na rysunku 6.2 stosując metodę Newtona, energetyczną
i równania Lagrange`a II rodzaju. Dane: k1 , k2 , k3 , k4 , k5 , k6 , G .
Rysunek 6.2 Model układu analizowany w przykładzie 1
W pierwszym etapie należy znalezć sztywność zastępczą elementów
sprężystych. Sprężyny k1 i k2 są połączone szeregowo. Ich sztywność
zastępcza k12 wynosi:
Strona 46
46
46
46
UKAADANIE RÓWNAC RUCHU
1 1 1 1 k1 + k2
= + stÄ…d =
k12 k1 k2 k12 k1 Å" k2
Czyli:
k1 Å" k2
k12 =
k1 + k2
Tak samo postępujemy w przypadku elementów sprężystych k3 i k4 .
k3 Å" k4
k34 =
k3 + k4
Elementy k12 , k34 , k5 i k6 stanowią sobą połączenie równoległe, zatem
sztywność wypadkowa kW jest równa sumie poszczególnych sztyw-
ności:
kW = k12 + k34 + k5 + k6
Zgodnie z zależnością (6.1) na układ działają siły:
G
&&
B = -mx = - &&
x - siła bezwładności,
g
G = mg - obciążenie układu,
S = kW (x + lst ) - reakcja elementów sprężystych,
gdzie:
G
lst = - ugięcie statyczne elementów sprężystych wywołane
kW
obciążeniem G ,
Zatem:
B + G - S = 0
Strona 47
47
47
47
ROZDZIAA 6
ëÅ‚ öÅ‚
G
&& ìÅ‚ ÷Å‚
- mx + mg - kW ìÅ‚ x + = 0
kW ÷Å‚
íÅ‚ Å‚Å‚
&&
- mx + mg - kW x - mg = 0
Ostatecznie:
G
&& &&
mx + kW x = 0 ; x + kW x = 0
g
Jest to równanie ruchu układu z rysunku 6.2 określone metodą Newtona.
Dla zastosowania metody energetycznej konieczne jest określenie
energii kinetycznej i potencjalnej.
Energia kinetyczna układu:
1 1 G
& &
EK = mx2 = x2
2 2 g
Energia potencjalna układu:
1
EP = kW x2
2
Korzystając z równania (6.3) otrzymujemy:
ëÅ‚
d 1 G 1 öÅ‚
&
ìÅ‚ x2 + kW x2 ÷Å‚ = 0
ìÅ‚ ÷Å‚
dt 2 g 2
íÅ‚ Å‚Å‚
1 G 1
&&& &
2 xx + 2kW xx = 0
2 g 2
G
&&
x + kW x = 0
g
Metoda równania Lagrange`a II rodzaju.
Korzystając z równania (6.4) otrzymujemy:
Strona 48
48
48
48
UKAADANIE RÓWNAC RUCHU
ëÅ‚ öÅ‚
d "EK d G G
ëÅ‚ öÅ‚
& &&
= ìÅ‚ x÷Å‚ = x
ìÅ‚ ÷Å‚
ìÅ‚ ÷Å‚
&
dt "x dt g g
íÅ‚ Å‚Å‚
íÅ‚ Å‚Å‚
"EP
= kW x
"x
"EK
= 0
"x
Q = 0
j
R = 0
j
Zatem:
G
&&
x + kW x = 0
g
We wszystkich trzech przypadkach otrzymaliśmy to samo równanie
opisujące ruch układu.
Jak wspomnieliśmy metoda Newtona sprawdza się przy stosunkowo
prostych układach drgających (zwykle o jednym stopniu swobody),
metoda energetyczna ma silne ograniczenia w postaci takiej iż układ
musi być autonomiczny. Najwygodniej jest stosować metodę równania
Lagrange`a II rodzaju, i w następnych przykładach będziemy właśnie ją
stosować.
Przykład 2.
Określić za pomocą równania Lagrange`a II rodzaju równanie ruchu
układu jak na rysunku 6.3.
Dane: sztywności elementów sprężystych k1 , k2 , k3 , k4 , k5 , k6 ,
Współczynniki tłumienia c1 , c2 , ciężar ciała G , P(t) siła wymusza-
jąca, x przemieszczenie ciała (współrzędna uogólniona).
Strona 49
49
49
49
ROZDZIAA 6
Rysunek 6.3 Model układu analizowany w przykładzie 2
Sztywność zastępcza kW :
1 1 1 1 k2k3 + k1k3 + k1k2
= + + =
k123 k1 k2 k3 k1k2k3
czyli:
k1k2k3
k123 =
k2k3 + k1k3 + k1k2
kW = k4 + k5 + k123
Zastępczy współczynnik tłumienia:
cW = c1 + c2
Energia kinetyczna rozpatrywanego układu:
1 1 G
& &
EK = mx2 = mx2
2 2 g
Energia potencjalna:
1
EP = kW x2
2
Strona 50
50
50
50
UKAADANIE RÓWNAC RUCHU
Siła oporu (rozpraszająca):
&
R = cW x
j
Siła wymuszająca:
Qj = P(t)
Wyznaczamy poszczególne człony równania Lagrange`a II rodzaju.
ëÅ‚ öÅ‚
d "EK d G G
ëÅ‚ öÅ‚
& &&
= ìÅ‚ x÷Å‚ = x
ìÅ‚ ÷Å‚
ìÅ‚ ÷Å‚
&
dt "x dt g g
íÅ‚ Å‚Å‚
íÅ‚ Å‚Å‚
"EK
= kW x
"x
Podstawiając do równania (6.4) otrzymujemy:
G
&& &
x + cW x + kW x = P(t)
g
Ostatecznie:
ëÅ‚ öÅ‚
G
ìÅ‚k4 ëÅ‚ k1k2k3 öÅ‚÷Å‚x = P(t)
&& & ìÅ‚ ÷Å‚÷Å‚
x + (c1 + c2 )x + + k5 +
ìÅ‚
ìÅ‚
g k2k3 + k1k3 + k1k2 ÷Å‚Å‚Å‚
íÅ‚ Å‚Å‚
íÅ‚
Przykład 3.
Określić równania ruchu układu jak na rysunku (6.4). Metodą równania
Lagrange`a II rodzaju.
Dane: G1 , G2 , k2 , k3 , P(t) , rysunek 6.4.
Układ ma dwa stopnie swobody, czyli może wykonywać dwa niezależne
od siebie ruchy. Przemieszczenia ciał o ciężarach G1 , G2 wynoszą x1,
x2 .
Strona 51
51
51
51
ROZDZIAA 6
Rysunek 6.4 Model układu analizowany w przykładzie 3
Energia kinetyczna układu:
1 1
& &
EK = m1x12 + m2 x2 2
2 2
G1 G2
m1 = ; m2 =
g g
Zatem energia kinetyczna wynosi:
1 G1 1 G2
& &
EK = x12 + x2 2
2 g 2 g
Energia potencjalna indukowana w sprężynach wynosi:
1 1
EP = k2 (x1 - x2 )2 + k3x2 2
2 2
Poszczególne człony równania (6.4):
Strona 52
52
52
52
UKAADANIE RÓWNAC RUCHU
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
d "EK G1 d "EK G2
ìÅ‚ ÷Å‚ &&1 ìÅ‚ ÷Å‚ &&2
= x ; = x
ìÅ‚ ìÅ‚
& &
dt "x1 ÷Å‚ g dt "x2 ÷Å‚ g
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
Energię potencjalną możemy przedstawić w wygodnej formie:
1 1 1
EP = k2 x12 - k2 x1x2 + k2x2 2 + k3x2 2
2 2 2
zatem:
"EP
= k2x1 - k2 x2 = k2 (x1 - x2 )
"x1
"EP
= -k2 x1 + k2 x2 + k3x2 = k2 (x2 - x1) + k3x2
"x2
Ostatecznie równania ruchu przyjmują postać:
G1
&&1
x + k2 (x1 - x2 ) = P(t)
g
G2
&&2
x + k2 (x2 - x1) + k3x2 = 0
g
Przykład 4.
Określić równania ruchu dla układu jak na rysunku 6.5, metodą równa-
nia Lagrange`a II rodzaju. Układ posiada trzy stopnie swobody.
Rysunek 6.5 Model układu analizowany w przykładzie 4
Dane: Momenty bezwładności krążków I , I , I , sztywności na skrę-
1 2 3
canie wałów łączących k1, k , k , moment wymuszający drgania M(t).
2 3
Strona 53
53
53
53
ROZDZIAA 6
Energia kinetyczna układu:
1 1 1
& & &
EK = I1Õ12 + I2Õ2 2 + I3Õ3 2
2 2 2
gdzie: Õ1, Õ , Õ - kÄ…ty skrÄ™cania poszczególnych wałów o momen-
2 3
tach bezwładności tarcz I , I , I .
1 2 3
Energia potencjalna układu:
1 1
EP = k2 (Õ2 - Õ1)2 + k3 (Õ3 -Õ2 )2
2 2
Energię potencjalną przedstawmy w wygodnej do różniczkowania
formie:
1 1 1 1
EP = k2Õ2 2 + k2Õ1Õ2 + k2Õ12 + k3Õ3 2 + k3Õ2Õ3 + k3Õ2 2
2 2 2 2
Poszczególne człony równania (6.1).
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
d "EK && d "EK d "EK
ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚ && ìÅ‚ ÷Å‚ &&
= I1Õ1 ; = I2Õ2 ; = I3Õ3
ìÅ‚ ìÅ‚ ìÅ‚
& & &
dt "Õ1 ÷Å‚ dt "Õ2 ÷Å‚ dt "Õ3 ÷Å‚
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
"EP "EP
= -k2Õ2 + k2Õ1 ; = k2Õ2 + k2Õ1 - k3Õ3 + k3Õ2
"Õ1 "Õ2
"EP
= k3Õ3 + k3Õ2
"Õ3
Podstawiając do równania (6.4) otrzymamy:
&&
I1Õ1 + k2 (Õ1 - Õ2 ) = 0
&&
I2Õ2 + k2 (Õ2 -Õ1) + k3 (Õ2 -Õ3 ) = M (t)
&&
I3Õ3 + k3 (Õ3 -Õ2 ) = 0
Powyższy układ trzech równań opisuje ruch rozpatrywanego układu.
Strona 54
54
54
54
`
Siły w ruchu drgającym
7
ROZDZIAA 7
Ogólne równanie różniczkowe drgań układu o jednym stopniu swobody
możemy zapisać w formie:
&& &
mx + F(x, x,t) = 0 (7.1)
&
Dla wielu przypadków siła F(x, x,t) może być przedstawiona jako
superpozycja składników z których każdy będzie zależał od jednej
z wymienionych wielkości:
& &
F(x, x,t) = S(x) + R(x) + G(t) (7.2)
Uwzględniając ostatecznie równanie ruchu drgającego może być
przedstawione w następującej formie:
&& &
mx + R(x) + S(x) = f (t) (7.3)
gdzie: f (t) = -G(t)
&
Gdy siła sprężysta S(x) i siła tłumienia R(x) są liniowymi funkcjami
&
przemieszczenia x i prędkości x , równanie (7.3) możemy przedstawić
w formie:
&& &
mx + cx + kx = f (t) (7.4)
Siła która jest zależna od przemieszczenia , jako funkcja S(x)
nazywana jest siłą restytucyjną lub wznawiającą. Możemy
wyróżnić dwa rodzaje sił restytucyjnych:
Strona 56
56
56
56
SIAY W RUCHU DRGAJCYM
a) grawitacyjna
Rysunek 7.1 Charakter siły restytucyjnej grawitacyjnej
SkÅ‚adowa mg cosÕ , napiÄ™cie nici o dÅ‚ugoÅ›ci l . SkÅ‚adowa
mg sinÕ jest siÅ‚Ä… wznawiajÄ…cÄ… grawitacyjnÄ….
b) sprężysta
Powstawanie sił restytucyjnych sprężystych jest związane z
właściwościami sprężystymi zastosowanych materiałów
konstrukcyjnych. Przykłady sił restytucyjnych sprężystych
przedstawiono na rysunku 7.2
W rozpatrywanych przez nas układach liniowych zależność pomiędzy
siłą sprężystą i przemieszczeniem ciała jest linią prostą. Zależność ta
występuje dla ciał które spełniają prawo Hook`a oraz przy małych
odkształceniach.
Siły zależne od prędkości są w drganiach siłami oporu (rozpraszają
energię). Skierowane są przeciwnie do zwrotu prędkości. Siły
rozpraszające powodują tłumienie drgań. Charakterystyka tłumienia to
zależność siły oporu od prędkości.
Strona 57
57
57
57
ROZDZIAA 7
l - długość belki
E - moduł Young`a
J - moment bezwładności
G - moduł Kirchoffa
Rysunek 7.2 Przykłady siły restytucyjnej sprężystej
W układach drgających liniowych siła oporu jest zależna liniowo od
prędkości. Mówimy wtedy o tak zwanym tarciu (tłumieniu)
wiskotycznym.
Strona 58
58
58
58
SIAY W RUCHU DRGAJCYM
& &
R(x) = cx (7.5)
gdzie: c - współczynnik tłumienia
Należy podkreślić, że siła oporu wiskotycznego występuje przy ruchu
ciała w płynie lepkim. Musi też być zachowany przepływ laminarny
(warstwowy) cieczy. Występuje to zwykle przy małych prędkościach
ciała.
Siły zależne tylko od czasu, a niezależne od przemieszczenia i prędkości
nazywamy siłami wymuszającymi. Siły te mogą mieć charakter
okresowy oraz krótkotrwały (impulsowy). Siły impulsowe
wyprowadzają układ z położenia równowagi, po czym drga on z
częstotliwością drgań własnych zależną od parametrów układu. Możemy
wyróżnić następujące typowe siły wymuszające:
a) Siła okresowa harmoniczna o stałej amplitudzie.
f (t) = Asin½t
gdzie:
A - stała amplituda,
½ - czÄ™stość siÅ‚y wymuszajÄ…cej,
t - czas.
b) Siła okresowa wynikająca z niewyważenia wirującego ciała względem
osi obrotu (wymuszenie bezwładnościowe).
2
f (t) = md r½ sin½t
gdzie:
md - niewyważone ciało o masie md ,
r - promień niewyważenia,
½ - czÄ™stość wymuszenia,
t - czas.
Siła ta jest szczególnie niebezpieczna, jej amplituda jak widać zależy od
kwadratu prędkości.
Strona 59
59
59
59
ROZDZIAA 7
c) Wymuszenie kinematyczne.
Rysunek 7.3 Idea wymuszenia kinematycznego
Punkt zamocowania sprężyny wykonuje ruch okresowy opisany
zależnością:
u = Asin½t
gdzie:
A - amplituda przemieszczenia punktu zamocowania sprężyny,
½ - czÄ™stość drgaÅ„ punktu zamocowania.
Zatem caÅ‚kowite odksztaÅ‚cenie sprężyny (¾ ) , bÄ™dzie różnicÄ… przemiesz-
czeń dolnego i górnego końca.
¾ = x - u = x - Asin½t
Siła sprężysta indukowana w sprężynie:
S(x) = -kx - kAsin½t
Widać że siła sprężysta może być rozdzielona na siłę zależną od
przemieszczenia ciała x i siłę zewnętrzną zależną od czasu.
d) Wymuszenie impulsowe.
W tym typie wymuszenia może to być jeden krótkotrwały impuls wytrą-
cający układ drgający z położenia równowagi lub seria impulsów
następujących po sobie.
Strona 60
60
60
60
SIAY W RUCHU DRGAJCYM
Rysunek 7.4 Przykłady wymuszeń impulsowych
Strona 61
61
61
61
ROZDZIAA 7
Strona 62
62
62
62
`
Krótka klasyfikacja
drgań
8
ROZDZIAA 8
Drgania klasyfikujemy w różny sposób. Przytoczmy klasyfikację
zaproponowaną przez Z. Osińskiego.
Według Z. Osińskiego możemy rozważać:
a) drgania o jednym stopniu swobody,
b) drgania o skończonej liczbie stopni swobody,
c) drgania układów o masach rozłożonych w sposób ciągły (nieskoń-
czenie wielka liczba stopni swobody).
Drgania mogą być:
a) swobodne, gdy nie ma siły wymuszającej, wymuszenie jest poprzez
warunki poczÄ…tkowe (poczÄ…tkowe przemieszczenie, poczÄ…tkowa
prędkość, zadane układowi drgającemu),
b) wymuszone, gdy układ drgający poddany jest działaniu jednej
z omawianych w poprzednim punkcie sił wymuszających,
c) samowzbudne, gdy układ nie jest poddany jawnemu działaniu siły
zewnętrznej, ale istnieje doprowadzenie energii sterowane przez sam
układ drgający.
Układy na które nie działają siły zewnętrzne nazywamy autonomiczny-
mi, a te na które działają siły zewnętrzne nieautonomicznymi.
Drganiami parametrycznymi nazywamy drgania układów w których
parametry takie jak masa lub sztywność zależą od czasu (najczęściej
w sposób okresowy). Układy te są opisane równaniami różniczkowymi
o zmiennych współczynnikach.
Jeżeli drgania opisane są przez równania różniczkowe liniowe, to
mówimy o drganiach liniowych. Ich charakterystyki sprężyste i tłumie-
nia sÄ… liniami prostymi.
Jeżeli charakterystyki sprężyste i tłumienia są nieliniowe, układ drgający
jest opisany równaniami różniczkowymi nieliniowymi i mamy wtedy do
czynienia z drganiami nieliniowymi.
Strona 64
64
64
64
KRÓTKA KLASYFIKACJA DRGAC
Drgania nazywamy tłumionymi, jeżeli w układzie drgającym występuje
rozproszenie energii, oraz nietłumionymi gdy nie ma rozproszenia
energii.
Strona 65
65
65
65
ROZDZIAA 8
Strona 66
66
66
66
`
Drgania swobodne
liniowego układu
drgajÄ…cego o jednym
9
stopniu swobody
(bez tłumienia)
ROZDZIAA 9
Rozpatrzmy układ przedstawiony na rysunku 9.1.
Rysunek 9.1 Rozpatrywany model układu
Układ wykonuje drgania pionowe. Dane: element sprężysty o sztywności
k , ciało o masie m .
Równanie ruchu układu:
&&
mx + kx = 0 (9.1)
Jeżeli podzielimy obie strony równania przez masę, otrzymamy:
&&
x + É0 2 x = 0 (9.2)
gdzie:
k
É0 = (9.3)
m
Zależność (9.3) nazywamy częstością drgań własnych.
Rozwiązanie równania (9.2) przewidujemy w postaci:
x = C1 cosÉ0t + C2 sinÉ0t (9.4)
Strona 68
68
68
68
DRGANIA SWOBODNE LINIOWEGO UKAADU DRGAJCEGO O JEDNYM STOPNIU SWOBODY
(BEZ TAUMIENIA)
Stałe C1 i C2 wyznaczamy z warunków początkowych albowiem
&
w chwili t = 0 , x(t = 0) = x0 , a prędkość x(t = 0) = V0 .
&
x = -C1É0simÉ0t + C2É0 cosÉ0t (9.5)
Stosując warunki początkowe na przemieszczenie z równania (9.4)
otrzymujemy:
x0 = C1 cosÉ0 0 + C2 sinÉ0 0 (9.6)
stÄ…d:
C1 = x0 (9.7)
Stosując warunek początkowy na prędkość otrzymujemy:
V0 = -C1É0 sin É0 0 + C2É0 cosÉ0 0 (9.8)
StÄ…d:
V0
C2 = (9.9)
É0
Zatem ostatecznie rozwiązanie z uwzględnieniem stałych C1 i C2 ma
postać:
V0
x = x0 cosÉ0t + sin É0t (9.10)
É0
Formułę (9.10) możemy zapisać w postaci:
x = asin(É0t +È ) (9.11)
V0
x0 = asinÈ ; = a cosÈ (9.12)
É0
Wyrażenia (9.12) podniesione do kwadratu i dodane stronami dają:
V0 2
a = x0 2 + (9.13)
É0 2
Strona 69
69
69
69
ROZDZIAA 9
x0É0
tgÈ = (9.14)
V0
Ciało będzie wykonywać ruch harmoniczny o stałej amplitudzie i fazie,
zależnej od warunków poczÄ…tkowych i czÄ™stoÅ›ci É0 zależnej od parame-
trów układu.
Przykład 1.
Wskazówka przyrządu pomiarowego ma masę m i zamocowana jest jak
na rysunku. Wskazówka wykonuje małe drgania wokół punktu 0 na
skali. Wyznaczyć częstość drgań własnych jeżeli sztywności sprężyn
podtrzymujących ją mają sztywność k , a sztywność sprężyny na
skręcanie w punkcie zamocowania wynosi kskr . Długość wskazówki
wynosi l .
Rysunek 9.2 Układ rozpatrywany w przykładzie 1
WspółrzÄ™dnÄ… okreÅ›lajÄ…cÄ… przemieszczenie koÅ„ca wskazówki jest kÄ…t Õ .
Energia kinetyczna wskazówki (ruch obrotowy wokół punktu A ):
1
2
&
EK = IÕ
2
gdzie: I - moment bezwładności wskazówki.
1
2
I = ml
3
Zatem ostatecznie energia kinetyczna wskazówki:
1
2 2
&
EK = ml Õ
6
Strona 70
70
70
70
DRGANIA SWOBODNE LINIOWEGO UKAADU DRGAJCEGO O JEDNYM STOPNIU SWOBODY
(BEZ TAUMIENIA)
Energia potencjalna związana z wychyleniem będzie magazynowana
w sprężynach k , k i kskr . Przemieszczenie sprężyn o współczynnikach
sztywności k , dla małych wychyleń wynosi:
x = lÕ
Zatem energia potencjalna:
1 1 1 1
2 2 2 2 2 2 2 2
EP = kl Õ + kl Õ + kskrÕ = kl Õ + kskrÕ
2 2 2 2
Wyliczając poszczególne człony równania Lagrange`a II rodzaju mamy:
d ëÅ‚ "EK öÅ‚ 1
2
&&
ìÅ‚ ÷Å‚ = ml Õ
ìÅ‚ ÷Å‚
&
dt "Õ 3
íÅ‚ Å‚Å‚
"EP 2
2
= 2kl Õ + kskrÕ = Õ(2kl + kskr )
"Õ
Ostateczne:
1
2 2
&&
ml Õ + Õ(2kl + kskr ) = 0
3
1
2
Po podzieleniu obu członów ostatniego równania przez wyrażenie ml
3
otrzymamy:
2
3(2kl + kskr )
&&
Õ + Õ = 0
2
ml
stÄ…d:
2
3(2kl + kskr )
É0 =
2
ml
Ponieważ mamy do czynienia z układem autonomicznym identyczny
wynik otrzymamy stosując metodę Newtona układania równań ruchu jak
i energetycznÄ….
Strona 71
71
71
71
ROZDZIAA 9
Przykład 2.
Dla układu jak na rysunku określić równania ruchu oraz częstość drgań
ciężaru G1 . Dane: G1 - ciężar drgający, G2 - ciężar krążka, y - współ-
rzędna określająca punkt zamocowania sprężyny, d - średnica krążka,
k - sztywność elementu sprężystego wznawiającego drgania.
Rysunek 9.3 Rozpatrywany układ drgający
Równania ruchu układamy korzystając z równania Lagrange`a II
rodzaju.
Ciężar G1 porusza się ruchem postępowym, zatem jego energia kine-
tyczna wyznaczona będzie zależnością:
1 1 G1 d
2
& & &
EK 1 = mV = x2 ; x = V = Õ
2 2 g 2
1 G1 2 & 2
EK 1 = d Õ
8 g
Wykonujący ruch obrotowy krążek o ciężarze G2 posiada energię
kinetycznÄ…:
1
2
&
EK 2 = IÉ ; É = Õ
2
gdzie:
Strona 72
72
72
72
DRGANIA SWOBODNE LINIOWEGO UKAADU DRGAJCEGO O JEDNYM STOPNIU SWOBODY
(BEZ TAUMIENIA)
2
m2r G2 2
d
I = =
2 g 8
Ostatecznie:
1 G2 2 & 2
EK 2 = d Õ
16 g
Zatem całkowita energia kinetyczna rozpatrywanego układu:
1 G1 2 2 G2 2 2 2
1 1 d
2
& & &
EK = EK 1 + EK 2 = d Õ + d Õ = (G2 + 2G1)Õ
8 g 16 g 16 g
Energia potencjalna zmagazynowana w sprężynie o sztywności k
wynosi:
x = yÕ - przemieszczenie zamocowanego koÅ„ca sprężyny do ciaÅ‚a
o ciężarze G2 .
1 1
2
EP = kx2 = ky2Õ
2 2
Poszczególne współczynniki równania Lagrange`a II rodzaju są:
2
d ëÅ‚ "EK öÅ‚ 1 d
&&
ìÅ‚ ÷Å‚ = (G2 + 2G1)Õ
ìÅ‚ ÷Å‚
&
dt "Õ 8 g
íÅ‚ Å‚Å‚
"EP
= ky2Õ
"Õ
Ostateczne:
2
1 d
&&
(G2 + 2G1)Õ + ky2Õ = 0
8 g
Po podzieleniu przez wyrażenie przy drugiej pochodnej otrzymujemy:
2
8ky g
&&
Õ + Õ = 0
2
d (G2 + 2G1 )
Strona 73
73
73
73
ROZDZIAA 9
stad:
8ky2 g
É0 =
2
d (G2 + 2G1)
Przykład 3.
Krążek którego walce o średnicy d wykonują drgania wokół najniższego
punktu toru będącego wycinkiem okręgu o średnicy D . Wyznaczyć
równanie ruchu krążka, przyjmując, że jego moment bezwładności
wynosi I , a ciężar G .
Rysunek 9.4 Układ analizowany w przykładzie 3
WspółrzÄ™dnÄ… okreÅ›lajÄ…cÄ… poÅ‚ożenie krążka bÄ™dzie kÄ…t Õ . Krążek bÄ™dzie
się poruszał ruchem obrotowym wokół osi w punkcie A i jednocześnie
będzie się przemieszczał po wycinku okręgu o średnicy D . Przemiesz-
czenie liniowe punktu A równe x wynosi:
D - d
x = Õ
2
a prędkość:
D - d
& &
x = Õ
2
Prędkość kątowa walca o średnicy d , wynosi:
Strona 74
74
74
74
DRGANIA SWOBODNE LINIOWEGO UKAADU DRGAJCEGO O JEDNYM STOPNIU SWOBODY
(BEZ TAUMIENIA)
&
x D - d 2 D - d
& &
É = = Õ = Õ
d
2 d d
2
Całkowita energia kinetyczna krążka wynosi:
1 1 1 (D - d)2 1 G (D - d)2
2 2 2 2
& &
EK = IÉ + mV = I Õ + Õ
2
2 2 2 2 g 4
d
Energia potencjalna:
G G
EK = mgh - mgh1 = g(D - d) - g(D - d) cosÕ =
g g
= G(D - d) - G(D - d) cosÕ = G(D - d)(1- cosÕ)
Współczynniki równania Lagrange`a II rodzaju:
d ëÅ‚ "EK öÅ‚ 1 (D - d)2 G (D - d)2
&&2 &&2
ìÅ‚ ÷Å‚ = I Õ + Õ =
ìÅ‚ ÷Å‚ 2
&
dt "Õ 2 d g 4
íÅ‚ Å‚Å‚
ëÅ‚ - d)2 G(D - d)2
öÅ‚
I (D
&&ìÅ‚
= ÕìÅ‚ 2 + ÷Å‚
÷Å‚
d 4
íÅ‚ Å‚Å‚
"EP
= G(D - d)sinÕ
"Õ
ëÅ‚ - d)2 G(D - d)2
öÅ‚
I(D
&&ìÅ‚ + G(D
ÕìÅ‚ 2 + ÷Å‚ - d)sinÕ = 0
÷Å‚
d 4
íÅ‚ Å‚Å‚
Po uproszczeniu oraz przyjÄ™ciu, że dla maÅ‚ych kÄ…tów Õ , sinÕ = Õ .
2
4Ggd
&&
Õ + Õ = 0
2
(D - d )(4gI + Gd )
CzÄ™stość drgaÅ„ wÅ‚asnych É0 wynosi:
2
4Ggd
É0 =
2
(D - d)(4gI + Gd )
Strona 75
75
75
75
ROZDZIAA 9
Przykład 4.
Obliczyć amplitudę drgań swobodnych podłużnych ciężaru Q zawieszo-
nego na końcu nieważkiego pręta pryzmatycznego o średnicy d . Prze-
mieszczenie początkowe x0 = 0,0003m , prędkość początkowa
m
V0 = 0,05
, ciężar Q = 1500N , długość pręta l = 1,25m , moduł
sek
N
Young`a E = 2 Å"1011 .
m2
Rysunek 9.5 Układ analizowany w przykładzie 4
Mając ciężar ciała możemy wyznaczyć masę m :
Q
m =
g
Sztywność pręta pryzmatycznego określona jest zależnością:
EF
k =
l
gdzie: F - pole przekroju poprzecznego.
Dla pręta o przekroju okrągłym
2
Ä„d
F =
4
Strona 76
76
76
76
DRGANIA SWOBODNE LINIOWEGO UKAADU DRGAJCEGO O JEDNYM STOPNIU SWOBODY
(BEZ TAUMIENIA)
Zatem ostatecznie:
2
EÄ„d
k =
4l
Częstość drgań własnych układu autonomicznego o jednym stopniu
swobody:
2
k EÄ„d g
É0 = =
m 4l Q
Zatem amplituda a :
V0 2 V0 2 4lQ
a = x0 2 + = x0 2 +
2
EÄ„d g
É0 2
Obliczenie wyniku pozostawiam czytelnikowi.
Przykład 5.
Określić moduł Kirchoffa G materiału metodą drgań skrętnych na pod-
stawie danych: długość pręta l = 1m , średnica d = 0,0125m , średnica
krążka D = 0,3m , ciężar krążka Q = 45N , zmierzona częstotli-
cykli
wość drgań swobodnych f = 10 = 10Hz .
sek
Rysunek 9.6 Badany układ
Strona 77
77
77
77
ROZDZIAA 9
Równanie ruch drgań skrętnych swobodnych:
&&
IÕ + kÕ = 0
DzielÄ…c przez I otrzymujemy:
&&
Õ + É0 2Õ = 0
k
É0 2 =
I
É0 = 2Ä„f
É0 1 k
f = =
2Ä„ 2Ä„ I
Sztywność k wynosi (skręcanie)
GJ0
k =
l
gdzie:
4
Ä„d
J0 = - moment bezwładności pręta.
32
czyli:
4
Ä„d G
k =
32l
Moment bezwładności krążka o średnicy D :
2
mr Q D2
I = =
2 g 8
czyli:
4
1 Ä„d G 8g
f =
2Ä„ 32l QD2
Strona 78
78
78
78
DRGANIA SWOBODNE LINIOWEGO UKAADU DRGAJCEGO O JEDNYM STOPNIU SWOBODY
(BEZ TAUMIENIA)
PodnoszÄ…c obustronnie do kwadratu:
2 4
4Ä„f 32lQD2 = d G8g
stÄ…d ostatecznie:
D2
2
G = 16Ä„f lQ
4
gd
Obliczenia pozostawiam czytelnikowi.
Strona 79
79
79
79
ROZDZIAA 9
Strona 80
80
80
80
`
Drgania swobodne
układu o jednym
stopniu swobody
10
tłumione tarciem
wiskotycznym
ROZDZIAA 5
Rozpatrzmy układ jak na rysunku.
Rysunek 10.1 Analizowany układ
k - sztywność elementu sprężystego,
c - współczynnik tłumienia wiskotycznego,
x - współrzędna określająca położenie ciała o masie m ,
g - przyspieszenie ziemskie.
Równanie ruchu ma postać:
&& &
mx + cx + kx = 0 (10.1)
Po podzieleniu obustronni przez m uzyskujemy:
&& & (10.2)
x + 2hx + É0 2 x = 0
gdzie:
c
= h - zredukowany współczynnik tłumienia,
2m
Strona 82
82
82
82
DRGANIA SWOBODNE UKAADU O JEDNYM STOPNIU SWOBODY
TAUMIONE TARCIEM WISKOTYCZNYM
k
É0 2 = - czÄ™stość drgaÅ„ wÅ‚asnych.
m
Rozwiązanie równania (10.2) Przewidujemy w postaci:
x = e-ht¾ (t) (10.3)
¾ (t) - nieznana funkcja której bÄ™dziemy poszukiwać.
Różniczkując dwustronnie zależność (10.3) otrzymujemy:
&
x = -he-ht¾ (t) + e-ht¾&(t) (10.4)
&
&&
x = h2e-ht¾ (t) - ¾&(t)e-hth - he-ht¾&(t) + e-ht¾&(t) (10.5)
Podstawiając wyrażenia (10.5), (10.4), (10.3) do równania (10.2)
otrzymujemy:
&
h2e-ht¾ (t) - ¾&(t)e-hth - he-ht¾&(t) + e-ht¾&(t) +
(10.6)
- 2h2e-ht¾ (t) + 2he-ht¾&(t) + É0 2e-ht¾ (t) = 0
Po podzieleniu równania (10.6) przez e-ht otrzymujemy:
&
h2¾ (t) - 2¾&(t)h + ¾&(t) - 2h2¾(t) + 2h¾&(t) + É0 2¾ (t) = 0 (10.7)
Ostatecznie po uproszczeniu mamy:
&
¾&(t) + (É0 2 - h2 )¾ (t) = 0 (10.8)
OznaczajÄ…c:
(É0 2 - h2 ) = p2 (10.9)
Otrzymujemy:
&
¾&(t) + p2¾ (t) = 0 (10.10)
Jest to klasyczne równanie jak dla drgań swobodnych nie tłumionych z
nowÄ… czÄ™stoÅ›ciÄ… p = É0 2 - h2 .
Rozwiązaniem równania (10.10) będzie wyrażenie:
Strona 83
83
83
83
ROZDZIAA 5
¾ (t) = C1 cos pt + C2 sin pt (10.11)
Zatem ogólne rozwiązanie równania (10.2) ma postać:
x(t) = ¾ (t)e-ht = e-ht (C1 cos pt + C2 sin pt) (10.12)
Stałe C1 i C2 wyznaczamy z warunków początkowych.
t = 0 ; x(t = 0) = x0 (10.13)
Uwzględniając ten warunek w równaniu (10.12) otrzymujemy:
x0 = 1(C1 cos p0 + C2 sin p0) (10.14)
StÄ…d :
C1 = x0
Różniczkując wyrażenie (10.12) otrzymujemy:
&
x(t) = -he-ht (C1 cos pt + C2 sin pt) +
(10.15)
+ e-ht (-C1 cos pt + C2 sin pt)
Uwzględniając drugi warunek początkowy w formie:
&
t = 0 ; x(t = 0) = V0 (10.16)
Otrzymujemy:
V0 = -h(C1 + C2 sin p0) + (-C1 cos p0 p + C2 sin p0 p) (10.17)
stÄ…d:
V0 = -hx0 + C2 p (10.18)
Ostatecznie:
V0 + hx0
C2 =
p
Pełne rozwiązanie równania (10.2) przyjmuje formę:
Strona 84
84
84
84
DRGANIA SWOBODNE UKAADU O JEDNYM STOPNIU SWOBODY
TAUMIONE TARCIEM WISKOTYCZNYM
V0 + hx0
x(t) = e-ht (x0 cos pt + sin pt) (10.19)
p
Z ostatniej zależności jednoznacznie wynika że postacie drgań swobod-
nych tłumionych zależą od charakteru tłumienia.
a) małe tłumienie, wtedy:
h < É0 oraz p2 = É0 2 - h2 < É0 (10.20)
PrzyjmujÄ…c:
V0 + hx0
x0 = asinÕ ; = a cosÕ (10.21)
p
Otrzymujemy:
x(t) = e-hta(cos pt sinÕ + sin pt cosÕ) = ae-ht sin( pt + Õ) (10.22)
gdzie:
(V0 + hx0 )2
a = x0 2 + (10.23)
p2
x0 p
tgÕ = (10.24)
hx0 +V0
Rysunek 10.2 Przebieg rozwiązania równania (10.2)
Strona 85
85
85
85
ROZDZIAA 5
Z ogólnego rozwiązania wynika że dla t " , x(t) 0 , to znaczy,
że drgania wygasają całkowicie po nieskończenie długim czasie.
Wielkość:
2Ä„ 2Ä„
Th = =
(10.25)
p
É0 2 - h2
nazywamy okresem drgań tłumionych.
Okres drgań nie tłumionych:
2Ä„
T = (10.26)
É0
É0 2 - h2 < É0 , zatem Th > T
Stosunek:
x(t)
h
= ehT (10.27)
x(t + Th )
Jest niezależny od czasu i jest równy stosunkowi kolejnych maksymal-
nych wychyleń w czasie jednego okresu drgań.
Wielkość:
x(t)
´ = ln = hTh (10.28)
x(t + Th )
nazywamy logarytmicznym dekrementem tłumienia i jest miarą
tłumienia w układzie.
b) tłumienie krytyczne:
h = hkr = É0 ; p = 0
wtedy współczynnik tłumienia ckr ma postać:
ckr = 2hkrm = 2mÉ0 = 2 mk (10.29)
W tym przypadku postać ruchu swobodnego:
Strona 86
86
86
86
DRGANIA SWOBODNE UKAADU O JEDNYM STOPNIU SWOBODY
TAUMIONE TARCIEM WISKOTYCZNYM
V0 + hx0
x(t) = e-ht (x0 cos pt + sin pt) =
pt
îÅ‚ sin pt Å‚Å‚
= limïÅ‚e-ht (x0 cos pt + (V0 + hx0 )t )śł = (10.30)
p0
pt
ðÅ‚ ûÅ‚
= e-ht (x0 + (V0 + hx0 )t)
Jest to ruch niedrgajÄ…cy zanikajÄ…cy z czasem. TÅ‚umienie krytyczne
wyznacza granicę pomiędzy drganiami harmonicznymi a ruchem
niedrgajÄ…cym.
c) duże tłumienie:
h > É0 ; p2 = É0 2 - h2 < 0
p , jest zatem wartością urojoną.
p = i É0 2 - h2 = ip (10.31)
RozwiÄ…zanie uzyskamy stosujÄ…c podstawienie:
cosix = cosh x ; sin ix = sinh x (10.32)
Tak więc:
V0 + hx0
x(t) = e-ht (x0 cosh pt + sinh pt) (10.33)
p
Jest to też ruch aperiodyczny co ilustruje rysunek 10.3.
Rysunek 10.3 Ilustracja ruchu, dla h > É0
Strona 87
87
87
87
ROZDZIAA 5
Przykład 1.
Wyznaczyć częstość i okres drgań układu mechanicznego przedstawio-
nego na rysunku 10.4. Dane są wielkości a i l wyznaczające zamoco-
wania tłumika i elementu sprężystego oraz masa pręta m . Masa skupio-
2
na ciała M = m . Znany jest współczynnik sztywności k sprężyny
3
i wiadomo że siła tarcia jest proporcjonalna do pierwszej potęgi prędkoś-
&
ci (tarcie wiskotyczne), R = Ä… Å" x , Ä… - współczynnik proporcjonal-
ności.
Rysunek 10.4 Ilustracja do przykładu 1
Równanie ruchu ma postać:
2
&& &
IÕ + Ä…a2Õ + kl Õ = 0
gdzie:
2 2
ml ml 2
2 2 2
I = + Ml = + ml = ml
3 3 3
Po podzieleniu równania ruchu przez I otrzymujemy:
&& &
Õ + 2hÕ + É0 2Õ = 0
gdzie:
Ä…a2 Ä…a2
2h = czyli h =
2 2
ml 2ml
Strona 88
88
88
88
DRGANIA SWOBODNE UKAADU O JEDNYM STOPNIU SWOBODY
TAUMIONE TARCIEM WISKOTYCZNYM
CzÄ™stość drgaÅ„ wÅ‚asnych É0 :
k
É0 =
m
Znając częstość drgań własnych oraz współczynnik tłumienia można
wyznaczyć częstość drgań tłumionych oraz okres drgań:
2Ä„ 2Ä„
p = É0 2 - h2 ; Th = =
p
É0 2 - h2
Aby wystąpił ruch aperiodyczny (nieokresowy) musi być spełniona
zależność:
h e" É0
czyli:
Ä…a2 k
e"
2
2ml m
Tak więc:
2
2l
Ä… e" km
a2
Przykład 2.
Ciężar Q zawieszony na sprężynie o sztywności k i zanurzony w ieczy
stawiającej opór wiskotyczny wykonuje drgania pionowe. Doświadczal-
nie zmierzono iż amplituda tych drgań po czterech wahnięciach zmalała
12-krotnie. Obliczyć okres drgań tłumionych i wyznaczyć logarytmiczny
N
dekrement tłumienia. Dane: Q = 50N , k = 2000 . Analizowany
m
układ przedstawia rysunek 10.5.
Strona 89
89
89
89
ROZDZIAA 5
Rysunek 10.5 Układ ilustrujący przykład 2
Amplitudowe wymuszenie w n tym okresie możemy przedstawić
zależnością:
an = ae-h(t +nTh )
A po czterech okresach:
h
an+4 = ae-h(t+(n+4)T )
gdzie: Th - okres drgań tłumionych.
Utwórzmy stosunek:
an+4
h
= e-4hT
an
LogarytmujÄ…c obustronnie otrzymujemy:
an+4 an
- 4hTh = ln = -ln
an an+4
czyli :
1 an
hTh = ´ = ln
4 an+4
PodstawiajÄ…c dane otrzymujemy:
Strona 90
90
90
90
DRGANIA SWOBODNE UKAADU O JEDNYM STOPNIU SWOBODY
TAUMIONE TARCIEM WISKOTYCZNYM
1
hTh = ´ = ln12 E" 0.622
4
Z ostatniej zależności możemy wyznaczyć współczynnik tłumienia h :
´ 0.622 1
îÅ‚ Å‚Å‚
h = =
śł
Th Th ïÅ‚ s
ðÅ‚ ûÅ‚
Podstawiając do wzoru na okres drgań tłumionych w formie:
2Ä„
Th =
É0 2 - h2
otrzymujemy:
1
2
Th = 4Ä„ + 0.6222
É0
gdzie:
k kg 2000 Å"9.81 1
É0 = = = = 19.8îÅ‚ Å‚Å‚
ïÅ‚s śł
m Q 50
ðÅ‚ ûÅ‚
Po podstawieniu do zależności na okres drgań tłumionych otrzymujemy:
1
2
Th = 4Ä„ + 0.6222 E" 0.319[s]
19.8
Strona 91
91
91
91
ROZDZIAA 5
Strona 92
92
92
92
`
Drgania wymuszane
układu o jednym
stopniu swobody
11
bez tłumienia
ROZDZIAA 11
Rozpatrywany układ drgający przedstawiono na rysunku 11.1. Układ
wykonuje drgania pionowe. Do ciała o masie m zawieszonego na ele-
mencie sprężystym o sztywności k przyłożona jest siła zależna od
czasu P(t) = Acos½t , gdzie A - amplituda siÅ‚y , ½ - czÄ™stość siÅ‚y
wymuszajÄ…cej.
Rysunek 11.1 Analizowany układ
Równanie ruchu układu:
&&
mx + kx = Acos½t (11.1)
Po podzieleniu przez m otrzymujemy:
&&
x + É0 2 x = q cos½t (11.2)
gdzie:
k
É0 =
- częstość drgań własnych,
m
Strona 94
94
94
94
DRGANIA WYMUSZANE UKAADU O JEDNYM STOPNIU SWOBODY BEZ TAUMIENIA
A
q =
m
Równanie (11.2) jest niejednorodnym równaniem różniczkowym. Jego
rozwiązanie jest superpozycją rozwiązania ogólnego równania jednorod-
nego xO oraz rozwiązania szczególnego równania niejednorodnego xS ,
czyli:
x = xO + xS (11.3)
Rozwiązanie ogólne ma postać:
xO = C1 cosÉ0t + C2 sinÉ0t (11.4)
Rozwiązania szczególnego poszukujemy w postaci:
xS = H cos½t (11.5)
gdzie: H - należy wyznaczyć.
Podstawiając wyrażenie (11.5) do równania (11.2) otrzymujemy:
&
xS = -H½ sin½t
(11.6)
2
&&S -H½ cos½t
x =
2
- H½ cos½t + É0 2H cos½t = q cos½t (11.7)
Po uproszczeniu:
2
[H (É0 2 +½ ) - q]cos½t = 0 (11.8)
Aby równanie (11.8) było spełnione dla wszystkich wartości t :
2
H (É0 2 +½ ) - q = 0 (11.9)
czyli:
q
H = (11.10)
2
É0 2 -½
Tak więc poszukiwane rozwiązanie szczególne ma postać:
Strona 95
95
95
95
ROZDZIAA 11
q
xS = cos½t (11.11)
2
É0 2 -½
Zatem całkowite rozwiązanie według (11.3) wygląda:
q
x = C1 cosÉ0t + C2 sinÉ0t + cos½t (11.12)
2
É0 2 -½
Widać że ruch ciała o masie m jest sumą dwóch ruchów harmonicz-
nych. Drgań swobodnych nietłumionych i drgań wymuszonych o czę-
stoÅ›ci ½ . Amplituda drgaÅ„ wymuszonych H opisana jest zależnoÅ›ciÄ…:
q 1 1
H = = lst 2 (11.13)
É0 2 2
½ ½
1- 1-
É0 2 É0 2
A
A q
m
lst = = =
(11.14)
k
k
É0 2
m
gdzie lst jest to przemieszczenie statyczne badanego ciała o masie m ,
pod wpływem siły o amplitudzie siły wymuszającej działającej w sposób
statyczny. Wprowadzmy pojęcie współczynnika uwielokrotnienia ampli-
tudy µ . Jest to stosunek amplitudy drgaÅ„ H do statycznego przemiesz-
czenia lst jakie wywołała by statycznie przyłożona do układu siła, równa
amplitudzie siły wymuszającej, czyli:
H 1
µ = = (11.15)
lst ½ 2
1-
É0 2
Przebiegi współczynnika uwielokrotnienia amplitudy µ w funkcji sto-
sunku częstości siły wymuszającej do częstotliwości drgań własnych ob-
razuje rysunek 11.2.
Strona 96
96
96
96
DRGANIA WYMUSZANE UKAADU O JEDNYM STOPNIU SWOBODY BEZ TAUMIENIA
Rysunek 11.2 Przebieg współczynnika uwielokrotnienia µ
½
w funkcji stosunku
É0
Widać, że dla czÄ™stoÅ›ci siÅ‚y wymuszajÄ…cej ½ = 0 , µ = 1, dla stosunku
½
= 1, µ dąży do nieskoÅ„czonoÅ›ci, wystÄ™puje zjawisko tzw. rezonan-
É0
½
su. Gdy dąży do nieskoÅ„czonoÅ›ci µ dąży do zera.
É0
Zajmijmy siÄ™ przypadkiem, gdy ½ = É0 , czyli przypadkiem rezonansu.
Rozwiązanie szczególne równania ruchu (11.2) ma wtedy formę:
xS = Ht sinÉ0t (11.16)
Różniczkując dwukrotnie otrzymujemy:
&&S -HÉ0 2t sin É0t + 2HÉ0 cosÉ0t
x =
(11.17)
Po podstawieniu do równania (11.2) otrzymujemy:
2HÉ0 cosÉ0t = q cosÉ0t
(11.18)
Po uproszczeniu:
Strona 97
97
97
97
ROZDZIAA 11
2HÉ0 cos É0t - q cos É0t = 0 (11.19)
(2HÉ0 - q) cosÉ0t = 0 (11.20)
2HÉ0 = q (11.21)
czyli:
q
H = (11.22)
2É0
Zatem ostatecznie:
q
xS = t sinÉ0t (11.23)
2É0
Drgania wymuszone dla omawianego przypadku nie sÄ… harmonicz-
ne, można je traktować jako drgania okresowe o narastającej am-
plitudzie proporcjonalnie do czasu.
Rysunek 11.3 Przebieg drgaÅ„ wymuszonych w przypadku ½ = É0
Stałe C1 i C2 występujące w równaniu (11.12) należy wyznaczyć
z warunków początkowych w formie; dla t = 0 ; x(t = 0) = 0 i dla
&
t = 0 ; x(t = 0) = 0 . Wyznaczenie stałych pozostawiam czytelnikowi.
Strona 98
98
98
98
DRGANIA WYMUSZANE UKAADU O JEDNYM STOPNIU SWOBODY BEZ TAUMIENIA
Przykład 1.
Ciało o ciężarze Q , zawieszone na nieważkim pręcie pryzmatycznym,
wykonuje drgania podłużne wymuszone siłą sinusoidalnie zmienną
w czasie. Dane Q = 50N , wydłużenie statyczne pręta pod wpływem tej
siÅ‚y lst = 0.025m , P(t) = P0 sin½t , P0 = 10N , liczba cykli siÅ‚y wy-
muszającej f = 5Hz . Obliczyć:
a) współczynnik uwielokrotnienia amplitudy µ ,
b) całkowite przemieszczenie ciężaru Q po upływie czasu t = 1s od
chwili poczÄ…tkowej ruchu. Warunki poczÄ…tkowe: dla t = 0 ;
&
x(t = 0) = 0 i dla t = 0 ; x(t = 0) = 0 .
Rysunek 11.4 Analizowany układ drgający
CzÄ™stość drgaÅ„ wÅ‚asnych É0 wynosi:
k
É0 =
m
Q
m =
g
Wydłużenie statyczne określone jest zależnością :
Q Q
lst = skÄ…d k =
k lst
Strona 99
99
99
99
ROZDZIAA 11
Podstawiając obliczoną sztywność i masę do wzoru na częstość drgań
własnych otrzymujemy:
g 9.81 1
É0 = = E" 19.8
lst 0.025 s
Częstość kątowa siły wymuszającej:
1
½ = 2Ä„f = 2Ä„ 5 = 10Ä„ = 31.4
s
Współczynnik uwielokrotnienia amplitudy µ :
1 1
µ = = E" 0.66
2
19.82
½
1-
1-
É0 2 31.42
Pod wpływem Q , ciężar przemieszcza się o lst , a pod wpływem P0
o ´0 , czyli:
Q P0
k = =
lst ´0
StÄ…d:
P0
´0 = lst
Q
Po podstawieniu danych liczbowych:
10
´0 = Å" 0.025 = 0.005m
50
Zatem amplituda drgań H wynosi:
H = µ´0 = 0.66 Å" 0.005 E" 0.0033m
Strona 100
100
100
100
DRGANIA WYMUSZANE UKAADU O JEDNYM STOPNIU SWOBODY BEZ TAUMIENIA
Chcąc obliczyć przemieszczenie x w danej chwili musimy w rozwiąza-
niu równania ruchu (11.12) wyznaczyć stałe C1 i C2 .
Po zastosowaniu przyjętych warunków początkowych:
H½
C1 = - ; C2 = 0
É0
Wobec tego równanie (11.12) przyjmuje postać:
H½
x = - cosÉ0t + H cos½t
É0
StÄ…d:
îÅ‚ Å‚Å‚
½
x = H - cosÉ0tśł
ïÅ‚cos½t
É0 ûÅ‚
ðÅ‚
Ponieważ:
1 1 ½
½ = 10Ä„ ; É0 = 19.8 H" 6.3Ä„ ; E" 1.59
s s É0
Więc dla t = 1sek , otrzymamy:
x(t = 1s) = H (cos10Ä„ Å"1-1.59cos 6.3Ä„ Å"1) =
= 0.0033(1-1.59 Å" cos0.3Ä„ ) H" -0.00422m
Strona 101
101
101
101
ROZDZIAA 11
Strona 102
102
102
102
`
Drgania wymuszone
liniowego układu
drgajÄ…cego o jednym
stopniu swobody
12
z tłumieniem
wiskotycznym
ROZDZIAA 12
Analizowany układ drgający przedstawia rysunek 12.1.
Rysunek 12.1 Analizowany układ drgający
k - sztywność elementu sprężystego,
c - współczynnik tłumienia wiskotycznego,
x - przemieszczenie,
m - masa ciała,
g - przyspieszenie ziemskie,
P(t) = Asin½t - siÅ‚a wymuszajÄ…ca o amplitudzie A i czÄ™stoÅ›ci ½ .
Równanie ruchu układu:
&& &
mx + cx + kx = Asin½t (12.1)
Po podzieleniu przez m otrzymujemy:
&& &
x + 2hx + É0 2 x = qsin½t (12.2)
gdzie:
c A
k
h = ; q = ; É0 =
m
2m m
Strona 104
104
104
104
DRGANIA WYMUSZANE LINIOWEGO UKAADU DRGAJCEGO O JEDNYM STOPNIU SWOBODY
Z TAUMIENIEM WISKOTYCZNYM
Podobnie jak w rozdziale 11 rozwiązanie równania (12.2) jest superpo-
zycją dwóch rozwiązań równań, ogólnego i szczególnego.
x = xO +xS (12.3)
Rozwiązanie ogólne:
xO = e-ht (C1 cos pt + C2 sin pt) (12.4)
Rozwiązanie szczególne przyjmujemy w formie:
xS = C3 sin½t + C4 cos½t (12.5)
Poszczególne pochodne mają postać:
&
xS = C3½ cos½t - C4½ sin½t (12.6)
2 2
&&S -C3½ sin½t - C4½ cos½t (12.7)
x =
Podstawiając zależności (12.5), (12.6), (12.7) do równania (12.2)
otrzymujemy:
2 2
- C3½ sin½t - C4½ cos½t + 2h½C3 cos½t - 2h½C4 sin½t +
(12.8)
+ É0 2C3 sin½t + É0 2C4 cos½t - q sin½t = 0
Zgrupujmy wyrażenia z cos½t i sin½t .
2
(C3 (É0 2 -½ ) - 2h½C4 - q)sin½t = 0 (12.9)
2
(C3 2h½ - (É0 2 -½ )C4 )cos½t = 0 (12.10)
Wyrażenia powyższe będą się zerować dla każdego t , gdy elementy
w nawiasach będą równe zeru, czyli:
2
C3 (É0 2 -½ ) - 2h½C4 - q = 0 (12.11)
2
C3 2h½ - (É0 2 -½ )C4 = 0 (12.12)
Z równania (12.12) wyznaczamy wielkość C4 .
Strona 105
105
105
105
ROZDZIAA 12
C3 2h½
C4 = (12.13)
2
(É0 2 -½ )
Po podstawieniu do równania (12.11) otrzymujemy:
2
C3 4h2½
2
C3 (É0 2 -½ ) - - q (12.14)
2
(É0 2 -½ )
Stad:
2
q(É0 2 -½ )
C3 = (12.15)
2 2
(É0 2 -½ ) + 4h2½
Mając wyrażenie C3 łatwo wyprowadzić z zależności (12.13) stałą C4 .
- 2h½q
C4 = (12.16)
2 2
(É0 2 -½ ) + 4h2½
Zatem rozwiązanie szczególne przyjmuje postać:
2
q(É0 2 -½ ) 2h½q
xS = sin½t - cos½t (12.17)
2 2 2 2
(É0 2 -½ ) + 4h2½ (É0 2 -½ ) + 4h2½
PrzyjmujÄ…c:
2
q(É0 2 -½ )
= H cosÕ (12.18)
2 2
(É0 2 -½ ) + 4h2½
- 2h½q
= H sinÕ (12.19)
2 2
(É0 2 -½ ) + 4h2½
i podstawiając do wyrażenia (12.17) otrzymujemy:
xS = H cosÕ sin½t + H sinÕ cos½t = H sin(½t + Õ) (12.20)
gdzie:
Strona 106
106
106
106
DRGANIA WYMUSZANE LINIOWEGO UKAADU DRGAJCEGO O JEDNYM STOPNIU SWOBODY
Z TAUMIENIEM WISKOTYCZNYM
2 2
q2 (É0 2 -½ )2 + 4h2½ q2
H = C3 2 + C4 2 = =
2 2
((É0 2 -½ ) + 4h2½ )2
(12.21)
2 2
(É0 2 -½ )2 + 4h2½ q
= q =
2 2
((É0 2 -½ ) + 4h2½ )2 (É0 2 2 ) + 4h2½ 2
-½
Wydłużenie statyczne:
A
A q
m
lst = = = (12.22)
k
k
É0 2
m
Wyrażenie (12.21) możemy przedstawić w postaci:
q
q É0 2
H = =
(12.23)
2 2
2 2 2 2
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
½ 4h2½ 4h2½
ìÅ‚1- ½ ÷Å‚
É0 2 ìÅ‚1- ÷Å‚ + +
ìÅ‚ ÷Å‚
É0 2 É0 4 ìÅ‚ É0 2 ÷Å‚ É0 4
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
Zatem:
lst
xS = sin(½t + Õ) (12.24)
2
2 2
ëÅ‚ öÅ‚
4h2½
ìÅ‚1- ½ ÷Å‚
+
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
Przesunięcie fazowe:
2 2
C4 - 2h½q (É0 2 -½ ) + 4h2½
tgÕ = = Å" =
2 2
C3 2 -½ ) + 4h2½ q(É0 2 -½ )
(É0 2
2h½
(12.25)
2h½ É0 2
= =
2
2
½
É0 2 -½
1-
É0 2
Strona 107
107
107
107
ROZDZIAA 12
Współczynnik uwielokrotnienia amplitudy µ :
H lst 1
µ = = Å" =
lst ëÅ‚ 2 öÅ‚2 lst
2
4h2½
ìÅ‚1- ½ ÷Å‚
+
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
(12.26)
1
=
2
2 2
ëÅ‚ öÅ‚
4h2½
ìÅ‚1- ½ ÷Å‚
+
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
Rysunek 12.2 Zmiana współczynnika uwielokrotnienia amplitudy
½
w funkcji (krzywe rezonansowe)
É0
Jak widać wprowadzone tłumienie h `" 0 powoduje ograniczenie
amplitud w rezonansie. Jednocześnie wierzchołki krzywych rezonanso-
½
wych przesuwają się w kierunku niższych stosunków .
É0
Zatem całkowite rozwiązanie równania (12.2) przyjmuje postać:
Strona 108
108
108
108
DRGANIA WYMUSZANE LINIOWEGO UKAADU DRGAJCEGO O JEDNYM STOPNIU SWOBODY
Z TAUMIENIEM WISKOTYCZNYM
x = xO + xS = e-ht (C1 cos pt + C2 sin pt) +
q
É0 2
(12.27)
+
2
2 2
ëÅ‚ öÅ‚
4h2½
ìÅ‚1- ½ ÷Å‚
+
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
Stałe C1 i C2 w równaniu trzeba wyznaczyć z przyjętych warunków
początkowych. Człon pierwszy równania (12.27) opisuje drgania
tłumione z częstością p , drugi zaś opisuje drgania wymuszone z czę-
stoÅ›ciÄ… siÅ‚y wymuszajÄ…cej ½ .
Przykład 1.
Ciężar Q zawieszony na końcu sprężyny o współczynniku sztywności
k wykonuje drgania podÅ‚użne wymuszone siÅ‚Ä… P(t) =P0sin½t i tÅ‚u-
mione oporem wiskotycznym . Obliczyć czÄ™stość siÅ‚y wymuszajÄ…cej ½
przy której zachodzi rezonans. W przypadku rezonansu obliczyć
amplitudę drgań wymuszonych oraz wartość współczynnika uwielokrot-
N
nienia amplitudy. Dane: Q = 4900N , k = 50000 , P0 = 2000N ,
m
&
siła oporu wiskotycznego tarcia R = 0.2 kmx .
Rysunek 12.3 Analizowany w przykładzie 1 układ drgający
Strona 109
109
109
109
ROZDZIAA 12
Z porównania z zapisem (7.5) siły tarcia wiskotycznego wynika,
c
że c = 0.2 km . Zarazem współczynnik tłumienia h = , stąd:
2m
0.2 km k
h = = 0.1
2m m
Ponieważ częstość drgań własnych:
k
É0 2 =
m
więc:
h = 0.1É0
Stosunek:
h É0
= 0.1, czyli = 10
É0 h
Obliczona częstość drgań własnych:
k kg 50000 Å"9.81 1
É0 = = = H" 10
m Q 4900 s
Przypadek rezonansu zachodzi gdy É0 =½ , czyli:
1
½ = 10
s
W przypadku rezonansu amplituda drgań wymuszonych wynosi:
lst lst lst
H = = =
2
2 2
4h2 k
2
ëÅ‚ öÅ‚
2
4h2½
(1-1) +
ìÅ‚1- ½ ÷Å‚
+
ìÅ‚ ÷Å‚ É0 2 É0
É0 2 É0 4
íÅ‚ Å‚Å‚
P0 2000
lst = = = 0.04m
k 50000
Strona 110
110
110
110
DRGANIA WYMUSZANE LINIOWEGO UKAADU DRGAJCEGO O JEDNYM STOPNIU SWOBODY
Z TAUMIENIEM WISKOTYCZNYM
Po podstawieniu danych otrzymujemy:
lst É0 0.04 Å"10
H½ =É0 = = = 0.2m
2 h 2
Przykład 2.
Ciężar Q umieszczony w środku belki o sztywności k wykonuje drga-
nia poprzeczne wymuszone siÅ‚Ä… zmiennÄ… w czasie P(t) =P0sin½t
i tłumione oporem wiskotycznym. Obliczyć ile zmaleje amplituda drgań
własnych w przypadku 3 krotnego wzrostu oporów tłumienia dla
N
danych: Q = 15000N , k = 1000000 , częstość siły wymuszającej
m
1
½ = 25.6 .
s
Rysunek 12.4 Rozpatrywany w przykładzie układ drgający
Pierwsza siła oporu wiskotycznego wynosi R1 . Zatem przy trzykrotnym
wzroście oporu tłumienia nowa siła R2 = 3R1 . Na podstawie zależności
(7.5) wynika ze w tym samym stosunku wzrośnie współczynnik
tłumienia h .
h2 = 3h1
Amplituda pierwotna określona jest zależnością:
lst
H1 =
2
2 2
ëÅ‚ öÅ‚
4h12½
ìÅ‚1- ½ ÷Å‚
+
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
Strona 111
111
111
111
ROZDZIAA 12
Po trzykrotnym wzroście oporu wiskotycznego amplituda wynosi:
lst
H =
2
2
2 2
ëÅ‚ öÅ‚
4h2 2½
ìÅ‚1- ½ ÷Å‚
+
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
Zatem stosunek amplitud H do H1 określony jest zależnością:
2
2
2 2
ëÅ‚ öÅ‚
4h2 2½
ìÅ‚1- ½ ÷Å‚
+
É0 2 É0 4
H2 ìÅ‚ ÷Å‚
íÅ‚ Å‚Å‚
=
2 2
H1 ëÅ‚ ½ öÅ‚2
4h12½
ìÅ‚1- ÷Å‚
+
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
Wyznaczona częstość drgań własnych wynosi:
k kg 1000000 Å"9.81 1
É0 = = = H" 25.6
m Q 15000 s
Widać że czÄ™stość drgaÅ„ wÅ‚asnych É0 jest równa czÄ™stoÅ›ci siÅ‚y wymu-
szajÄ…cej, czyli mamy do czynienia z przypadkiem rezonansu.
Zatem:
2
ëÅ‚ öÅ‚
h1
ìÅ‚ ÷Å‚
4ìÅ‚
É0 Å‚Å‚
H2 íÅ‚ ÷Å‚ h1 1
= = =
H1 ëÅ‚ h2 öÅ‚2 h2 3
ìÅ‚ ÷Å‚
4ìÅ‚ ÷Å‚
É0
íÅ‚ Å‚Å‚
czyli:
1
H2 = H1
3
Strona 112
112
112
112
`
Drgania liniowe
układu o jednym
stopniu swobody
wymuszane
13
bezwładnościowo
(z tłumieniem)
ROZDZIAA 13
Rozpatrzmy układ jak na rysunku 13.1.
Rysunek 13.1 Analizowany układ
Ciało o masie m pobudzane jest do drgań poprzez siłę Fy pochodzącą
od niewyważenia wzglÄ™dem osi 0 . droga kÄ…towa Õ okreÅ›lona jest
zależnością:
Õ =½ Å" t (13.1)
gdzie:
½ - prÄ™dkość kÄ…towa wirujÄ…cego krążka.
md - dodatkowe ciało wywołujące niewyrównoważenie względem
osi 0 .
Strona 114
114
114
114
DRGANIA LINIOWE UKAADU O JEDNYM STOPNIU SWOBODY
WYMUSZANE BEZWAADNOÅšCIOWO (Z TAUMIENIEM)
Ciało o masie m może wykonywać tylko drgania pionowe wywołane siłą
F .
y
Fy = F sinÕ = F sin½t (13,2)
2 2 2
mdV md½ r
2
F = = = md½ r (13.3)
r r
Zatem:
2
Fy = md r½ sinÕ (13.4)
Równanie ruchu analizowanego układu będzie miało postać:
2
&& &
mx + cx + kx = md r½ sin½t (13.5)
Po podzieleniu przez m , otrzymujemy:
2
&& & (13.6)
x + 2hx + É0 2x = q0½ sin½t
gdzie:
c k md r
h = , É0 = , q0 = (13.7)
2m m m
Rozwiązanie szczególne równania (13.6) ma postać:
2
q0½
xS = sin½t =
2
2 2
(É0 2 -½ ) + 4h2½
2
½
q0 2
(13.8)
É0
= sin½
2
2 2
ëÅ‚ öÅ‚
4h2½
ìÅ‚1- ½ ÷Å‚
+
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
Współczynnik uwielokrotnienia amplitudy µ :
Strona 115
115
115
115
ROZDZIAA 13
2 2
½ ½
q0 2
H É0 É0 2
µ = = == (13.9)
2
q0 ëÅ‚ 2 öÅ‚2
2 2 2
ëÅ‚ öÅ‚
4h2½ ½ 4h2½
ìÅ‚1- ½ ÷Å‚
+ q0 ìÅ‚1- ÷Å‚ +
ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚
É0 2 É0 4 É0 2 É0 4
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
Przebieg zmian współczynnika uwielokrotnienia amplitudy µ przedsta-
wia rysunek 13.2.
Rysunek 13.2 Zmiana współczynnika uwielokrotnienia amplitudy w
½
funkcji
É0
Wierzchołki rezonansu wraz ze wzrostem częstości siły wymuszającej
przesuwajÄ… siÄ™ w prawo .
Przykład 1.
Silnik elektryczny zamocowano za sprężystej belce przy czym jej
strzałka ugięcia jest y . Mimośród wirnika obracającego się z prędkością
kÄ…towÄ… ½ jest równy r , a masa wirnika m . Masa silnika elektrycznego
wraz z wirnikiem wynosi M . Znalezć amplitudy drgań pionowych wy-
muszonych silnika . Przy jakiej wartoÅ›ci ½ może wystÄ…pić rezonans.
Strona 116
116
116
116
DRGANIA LINIOWE UKAADU O JEDNYM STOPNIU SWOBODY
WYMUSZANE BEZWAADNOÅšCIOWO (Z TAUMIENIEM)
Rysunek 13.3 Badany układ niewyważonego wirnika
Masa niewyważonego wirnika wynosi m , masa wirnika plus masa
obudowy i belki wynosi M , g - przyspieszenie ziemskie.
Równanie ruchu układu:
&&
My = -S + F cosÕ
2
gdzie: F = mr½ , Õ =½t , S = ky
Pod wpÅ‚ywem caÅ‚kowitego obciążenia M Å" g belka wychyli siÄ™ o war-
tość y , zatem:
ky = Mg
Stąd możemy wyznaczyć sztywność belki k :
Mg
k =
y
Podstawiając do równania ruchu otrzymujemy:
Mg
2
&&
My = - x + m½ r cosÕ
y
i ostatecznie:
Mg
2
&&
My + y = m½ r cos½t
y
Po podzieleniu przez M mamy:
Strona 117
117
117
117
ROZDZIAA 13
g m
2
&&
x + x = ½ r cos½t
y M
Rozwiązanie szczególne przewidujemy w postaci:
xS = H cos½t
2
&&S -H½ cos½t
x =
Po podstawieniu do równania ruchu otrzymujemy:
g m
2 2
- H½ cos½t + H cos½t = ½ r cos½t
y M
Po uproszczeniu:
g m
2 2
- H½ + H = ½ r
y M
stÄ…d:
2
m½ ry
H =
2
M (g -½ y)
Ostatecznie rozwiązanie szczególne:
2
m½ ry
xS = cos½t
2
M (g -½ y)
g
2
Przypadek rezonansu wystÄ…pi gdy =½ .
f
RozwiÄ…zanie przyjmujemy w postaci:
xS = Ht sin½t
2
&&S
x = H½ cos½t + H½ cos½t - Ht½ sin½t
&&S
Ostatecznie x :
2
&&S
x = 2H½ cos½t - Ht½ sin½t
Strona 118
118
118
118
DRGANIA LINIOWE UKAADU O JEDNYM STOPNIU SWOBODY
WYMUSZANE BEZWAADNOÅšCIOWO (Z TAUMIENIEM)
Podstawiając do równania ruchu otrzymujemy:
g m
2 2
2H½ cos½t - Ht½ sin½t + Ht sin½t = r½ cos½t
f M
Po uproszczeniu:
m
2
2H½ cos½t = r½ cos½t
M
m
2
2H½ = r½
M
Ostatecznie amplituda w rezonansie wynosi:
2
mr½ mr½
H = =
2M½ 2M
a rozwiązanie szczególne dla przypadku rezonansu:
mr½
xS = t sin½t
2M
Doświadczenie 1.
Celem eksperymentu jest praktyczne zaznajomienie czytelnika z analizÄ…
drgań belki wymuszonych bezwładnościowo. Schemat stanowiska
przedstawiony jest na rysunku 13.4.
Rysunek 13.4 Schemat stanowiska
Strona 119
119
119
119
ROZDZIAA 13
Składa się on z następujących zespołów:
1 - Belki stalowej (1) utwierdzonej jednym końcem w podstawie.
2 - Układu wymuszającego składającego się z dwóch niewyważonych
względem osi obrotu kół zębatych (2), napędzanych silnikiem (3), patrz
rysunek 13.5.
3 - TÅ‚umika olejowego (5).
4 - Układu rejestrującego (6) służącego do rejestracji drgań belki z wy-
korzystaniem tensometru.
Rysunek 13.5 Schemat zespołu wymuszającego
W ćwiczeniu należy zarejestrować zmiany amplitudy drgań belki (1)
w funkcji zmiany częstości wymuszenia. Zmieniając wartości siły oporu
tłumienia można doświadczalnie wyznaczyć rodzinę krzywych rezonan-
sowych przedstawionych jako teoretyczne na rysunku 13.2.
Strona 120
120
120
120
`
Drgania układu
o jednym stopniu
swobody przy
wymuszeniu
14
kinematycznym
(z tłumieniem)
ROZDZIAA 14
Rozpatrzmy układ jak na rysunku 14.1.
Rysunek 14.1 Schemat analizowanego układu
Widać iż w tym przypadku wymuszenie jest przekazywane od profilu
drogi ¾ (t) poprzez element sprężysty i tÅ‚umik na ciaÅ‚o o masie m .
Zbadajmy przemieszczenie x ciała o masie m .
Równanie ruchu ciała o masie m ma postać:
&& &
mx + c(x - ¾&) + k(x - ¾ ) = 0 (14.1)
gdzie, zakładany profil drogi przyjmujemy w formie:
¾ = Asin½t (14.2)
¾& = A½ cos½t (14.3)
Podstawiając do równania (14.1) otrzymujemy:
&& &
mx + cx + kx = cA½ cos½t + kAsin½t (14.4)
Dzieląc równanie (14.4) przez m otrzymujemy:
Strona 122
122
122
122
DRGANIA UKAADU O JEDNYM STOPNIU SWOBODY PRZY WYMUSZENIU KINEMATYCZNYM
(Z TAUMIENIEM)
cA½ kA
&& &
x + 2hx + É0 2 x = cos½t + sin½t (14.5)
m m
gdzie:
* *
kA = H sinÕ ; c½A = H cosÕ (14.6)
PodstawiajÄ…c (14.6) do (14.4) otrzymujemy:
* *
H H
&& &
x + 2hx + É0 2 x = cos½t cosÕ + sin½t sinÕ (14.7)
m m
Zatem:
*
H
&& &
x + 2hx + É0 2 x = cos(½t - Õ) (14.8)
m
Korzystając z zależności (14.6) mamy:
* 2 2 2 2
H = k A2 + c2½ A2 = A k + c2½ (14.9)
kA k
tgÕ = = (14.10)
c½A c½
2
k c2
UwzglÄ™dniajÄ…c, że = É0 4 oraz = 4h2 otrzymujemy:
m2 m2
2
&& &
x + 2hx + É0 2 x = A É0 4 + 4h2½ cos(½t -Õ) (14.11)
Zatem amplituda H , i rozwiązanie szczególne wynosi:
2
A É0 4 + 4h2½
xS = cos½t (14.12)
2
2 2
(É0 2 -½ ) + 4h2½
gdzie:
Strona 123
123
123
123
ROZDZIAA 14
2
4h2½
AÉ0 2 1 +
2
A É0 4 + 4h2½ É0 4
*
H = = =
2 2
2 2
2 2
ëÅ‚ öÅ‚
(É0 2 -½ ) + 4h2½
½ h2½
÷Å‚
É0 2 ìÅ‚1 - + 4
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
(14.13)
2
4h2½
A 1 +
É0 4
=
2
2 2
ëÅ‚ öÅ‚
h2½
ìÅ‚1 - ½ ÷Å‚
+ 4
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
A wiÄ™c współczynnik uwielokrotnienia amplitudy µ wynosi:
2
4h2½
1 +
*
É0 4
H
*
µ = = (14.14)
2
A
2 2
ëÅ‚ öÅ‚
h2½
ìÅ‚1 - ½ ÷Å‚
+ 4
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
½
Przebieg µ w funkcji obrazuje rysunek 14.2.
É0
Strona 124
124
124
124
DRGANIA UKAADU O JEDNYM STOPNIU SWOBODY PRZY WYMUSZENIU KINEMATYCZNYM
(Z TAUMIENIEM)
½
Rysunek 14.2 Przebiegi zmian µ funkcji
É0
W przypadku drgań wymuszanych kinematycznie, bardziej interesujący
jest przypadek gdy badamy przemieszczenia względne ciała o masie m ,
to znaczy ugiÄ™cie elementu sprężystego y = (x - ¾ ) .
Przy wykorzystaniu takiego oznaczenia mamy:
y = x - ¾ (14.15)
& &
y = x - ¾& (14.16)
&
&& && - ¾&
y = x (14.17)
Równanie ruchu przyjmuje postać:
&& &
mx + cy + ky = 0 (14.18)
Ale:
&
&& &&
x = y + ¾& (14.19)
Więc:
&
&& &
my + cy + ky = -m¾& (14.20)
Strona 125
125
125
125
ROZDZIAA 14
2
&
Jeżeli ¾ = Asin½t , to ¾& = -A½ sin½t .
Po podstawieniu do zależności (14.20) otrzymujemy:
2
&& &
my + cy + ky = mA½ sin½t (14.21)
Dzieląc obydwie strony równania (14.21) przez m , mamy:
2
&& &
y + 2hy + É0 2 y = A½ sin½t
Stosując analizę z rozdziałów 12 i 13 otrzymujemy:
2
½
A
2
A½ É0 2
*
H = = (14.22)
2 2
2 2
2 2
ëÅ‚ öÅ‚
(É0 2 -½ ) + 4h2½
h2½
ìÅ‚1 - ½ ÷Å‚
+ 4
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
Zatem współczynnik uwielokrotnienia amplitudy µ :
2
½
*
H É0 2
*
µ = = (14.23)
2
A
2 2
ëÅ‚ öÅ‚
h2½
ìÅ‚1 - ½ ÷Å‚
+ 4
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
Przebieg zmian współczynnika uwielokrotnienia amplitudy µ przedsta-
wia rysunek 14.3.
Strona 126
126
126
126
DRGANIA UKAADU O JEDNYM STOPNIU SWOBODY PRZY WYMUSZENIU KINEMATYCZNYM
(Z TAUMIENIEM)
½
Rysunek 14.3 Przebiegi zmian µ funkcji dla różnych
É0
współczynników tłumienia
Wierzchołki krzywych rezonansowych przesuwają się w prawo
ze wzrostem stosunku czÄ™stoÅ›ci siÅ‚y wymuszajÄ…cej ½ i współczynnika
tłumienia h .
Przykład 1.
Ciężar jak na rysunku 14.4 porusza się ze stałą prędkością V wzdłuż
nierównej drogi. Wzdłużny profil drogi przedstawiony jest równaniem
S(t) . Znając masę ciężaru równą m , sztywność sprężyny k
o współczynnik tłumienia c , znalezć równanie drgań ciężaru;
Ä„¾
S(t) = s0 sin2 , długość fali jak na rysunku.
l
Strona 127
127
127
127
ROZDZIAA 14
Rysunek 14.4 Profil drogi i analizowany w przykładzie układ
Równanie drgań pionowych ciała o masie m :
&& & &
mz + c(z - s) + k(z - s) = 0
z - s = y
& & &
z - s = y
&& &
mz + cy + ky = 0
&&- && && && && &&
z s = y Ò! z = y + s
czyli:
&& && &
m( y + s) + cy + ky = 0
Zetem:
&& & &&
my + cy + ky = -ms
Po podzieleniu przez m otrzymujemy:
&& &
y + 2hy + É0 2 y = -&&
s
c k
gdzie: h = , É0 2 = , ¾ = Vt
2m m
Strona 128
128
128
128
DRGANIA UKAADU O JEDNYM STOPNIU SWOBODY PRZY WYMUSZENIU KINEMATYCZNYM
(Z TAUMIENIEM)
Zatem:
Ä„Vt
S(t) = s0 sin2
l
Ä„Vt Ä„Vt Ä„V Ä„V 2Ä„Vt
&
s = 2s0 sin cos = s0 sin
l l l l l
stÄ…d:
2 2
Ä„V 2Ä„Vt 2Ä„V Ä„ V 2Ä„Vt
&&
s = s0 cos = 2s0 2 cos
l l l l l
Zatem równanie ruchu przyjmuje postać:
2 2
Ä„ V 2Ä„Vt
&& &
y + 2hy + É0 2 y = -2s0 2 cos
l l
lub
&& &
y + 2hy + É0 2 y = - p0 cos½t
2 2
Ä„ V 2Ä„V
gdzie: p0 = 2s0 2 , ½ = .
l l
Rozwiązanie ostatniego równania ma postać:
y(t) = Ae-ht sin pt + Be-ht cos pt - C cos(½t -Õ)
Stała C dla rozwiązania szczególnego wynosi:
2 2
p0
2s0Ä„ V m
2
É0 2
l k
C = =
2 2
2 2 2 2
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
h2½ h2½
ìÅ‚1- ½ ÷Å‚ ìÅ‚1- ½ ÷Å‚
+ 4 + 4
ìÅ‚ ÷Å‚
É0 2 É0 4 ìÅ‚ É0 2 ÷Å‚ É0 4
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
2k½
tgÕ =
É0 2
Strona 129
129
129
129
ROZDZIAA 14
Stałe A i B wyznaczamy z warunków początkowych. Dla t = 0 ;
y(t = 0) = 0 .
0 = B - C cosÕ
stÄ…d:
B = C cosÕ
&
Dla t = 0 ; y(t = 0) = 0 :
0 = Ap - hB - C½ cosÕ
stÄ…d:
hB - C½ sinÕ C
A = = (h cosÕ +½ sinÕ)
p p
Zatem ostatecznie równanie drgań przyjmuje postać:
1
y(t) = Ce-ht (cosÕ cos pt + (hcosÕ +½ sinÕ))sin pt - C cos(½t -Õ)
p
Strona 130
130
130
130
`
Amortyzacja drgań
15
ROZDZIAA 15
Celem amortyzacji jest złagodzenie skutków drgań. Możemy tu wyróż-
nić dwa problemy (przypadki) tego zagadnienia:
a) Ochronę otoczenia przed skutkami drgań maszyny.
b) Ochronę maszyny od skutków drgań otoczenia.
Rozpatrzy problem a .
Mając układ jak na rysunku 15.1.
Rysunek 15.1 Analizowany układ
Wyznaczmy siłę przenoszoną na podłoże:
Równanie ruchu:
&& &
mx + cx + kx = H sin½t (15.1)
Rozwiązanie szczególne ma postać:
xS = Asin(½t + Å‚ ) (15.2)
gdzie:
Strona 132
132
132
132
AMORTYZACJA DRGAC
lst
A = (15.3)
2
2 2
ëÅ‚ öÅ‚
h2½
ìÅ‚1- ½ ÷Å‚
+ 4
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
Siła działająca na podłoże składa się z sił przenoszonych przez sprężynę
i przez tłumik. Siła sprężysta wynosi:
R = kx = kAsin(½t + Å‚ ) (15.4)
Siła przenoszona przez tłumik:
&
S = cx = c½Acos(½t + Å‚ ) (15.5)
Maksymalna wartość siły jest równa:
2
c½
ëÅ‚ öÅ‚
Pmax = (kA)2 + (c½A)2 = kA 1+ (15.6)
ìÅ‚ ÷Å‚
k
íÅ‚ Å‚Å‚
Podstawiając zależność na A otrzymujemy:
2
kst c½
ëÅ‚ öÅ‚
Pmax = 1 + =
ìÅ‚ ÷Å‚
2
k
íÅ‚ Å‚Å‚
2 2
ëÅ‚ öÅ‚
h2½
ìÅ‚1 - ½ ÷Å‚
+ 4
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
2 (15.7)
c½
ëÅ‚ öÅ‚
kst 1 +
ìÅ‚ ÷Å‚
k
íÅ‚ Å‚Å‚
= = kstµ
2
2 2
ëÅ‚ öÅ‚
h2½
ìÅ‚1 - ½ ÷Å‚
+ 4
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
kst = P (15.8)
st
WiÄ™c µ oznacza stosunek maksymalnej siÅ‚y przenoszonej na podÅ‚oże
w czasie drgań do siły statycznej i nazywa się go współczynnikiem
½
przenoszenia. Dla stosunku > 2 występuje zmniejszenie sił prze-
É0
noszonych na podłoże.
Strona 133
133
133
133
ROZDZIAA 15
½
Rysunek 15.2 Zmiana współczynnika przenoszenia µ w funkcji
É0
Aby uzyskać właściwą amplitudę należy sztywność układy k dobrać tak
by były spełnione warunki:
½ ½ m
= > 2 (15.9)
É0
k
Jak widać właściwie nie ma potrzeby stosowania tłumika, gdyż nie-
znacznie powiększa on siły. Tłumik niezbędny jest jednak w strefie
rezonansu, gdzie poważnie wpływa na zmniejszenie amplitudy sił. Przy
rozruchu maszyny często przechodzimy przez strefę rezonansu.
Rozpatrzmy przypadek b
Przeanalizujmy układ jak na rysunku 15.3.
Rysunek 15.3 Analizowany układ
Strona 134
134
134
134
AMORTYZACJA DRGAC
Załóżmy że podłoże porusza się ruchem harmonicznym:
s = Acos½t (15.10)
Równanie ruchu:
&& & &
mx + c(x - s) + k(x - s) = 0 (15.11)
&
s = -A½ sin½t (15.12)
PodstawiajÄ…c otrzymujemy:
&& &
mx + cx + kx = -Ac½ sin½t + kAcos½t (15.13)
&& &
mx + cx + kx = A(k cos½t - c½ sin½t) (15.14)
PodstawiajÄ…c:
a cosÕ = k ; a sinÕ = -c½ (15.15)
otrzymujemy:
2 2
&& &
mx + cx + kx = a k +½ c2 cos(½t + Õ) (15.16)
c½
2 2
a = k +½ c2 ; tgÕ = - (15.17)
k
Ostatecznie rozwiązanie równania (15.16) ma postać:
2
c½
ëÅ‚ öÅ‚
1+
ìÅ‚ ÷Å‚
k
íÅ‚ Å‚Å‚
D =
(15.18)
2
2 2
ëÅ‚ öÅ‚
h2½
ìÅ‚1- ½ ÷Å‚
+ 4
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
Stosunek maksymalnej amplitudy drgań własnych do przemieszczenia
statycznego jest taki sam jak w poprzednim przypadku. Tak więc zasady
amortyzacji sÄ… takie same.
Strona 135
135
135
135
ROZDZIAA 15
Strona 136
136
136
136
`
Rejestracja drgań
16
ROZDZIAA 16
Badanie i rejestracja drgań pozwala określić zródła ich powstania,
określić ich szkodliwość, a poprzez pomiar parametrów takich jak czę-
stość, amplituda wartości przyspieszeń oraz sił przewidzieć sposoby ich
zmniejszenia lub wręcz usunięcia. W przypadku gdy drgania są wyko-
rzystywane, wtedy ich badania pozwalają określić optymalne warunki
pracy.
Pomiary i rejestrację drgań dokonujemy w oparciu o dwie zasady:
a) Pomiar drgań badanego obiektu względem nieruchomego układu
odniesienia.
Rysunek 16.1 Ilustracja zasady pomiarowej a
b) Umieszczenie dodatkowego układu drgającego na obiekcie badanym i
pomiar drgań dodatkowego ciała względem jego obudowy związanej z
badanym obiektem.
Strona 138
138
138
138
REJESTRACJA DRGAC
Rysunek 16.2 Ilustracja zasady pomiarowej b
W przypadku a pomiar i rejestracja jest prosta , ale występują duże
trudności z utrzymaniem stałego położenia układu nieruchomego
(rejestrującego). W przypadku b nie trzeba stałego układu odniesienia
i stąd popularność tej metody rejestracji i pomiaru. Przyrządy pracujące
na zasadzie b nazywane sÄ… przyrzÄ…dami inercyjnymi lub sejsmicz-
nymi.
Bardzo ważnym problemem jest właściwe dobranie przyrządu sejsmicz-
nego tak by wskazywał rzeczywiste przemieszczenie lub przyspieszenie
badanego obiektu.
Rysunek 16.3 Schemat przyrzÄ…du sejsmicznego
Dodatkowe ciało o masie m w czujniku sejsmicznym w skutek ruchu
obudowy poddane jest wymuszeniu kinematycznemu. Ruch bezwzględ-
ny ciała opisany jest równaniem:
&& & &
mx + c(x - s) + k(x - s) = 0 (16.1)
gdzie: s - przemieszczenie podstawy.
Względne przemieszczenie ciała i podstawy określone jest zależnością:
Strona 139
139
139
139
ROZDZIAA 16
y = x - s (16.2)
& & &
y = x - s (16.3)
&& && - && && && && (16.4)
y = x s Ò! x = y + s
Zatem równanie ruchu względnego:
&& && &
my + ms + cy + ky = 0 (16.5)
Po podzieleniu przez m :
&& && &
y + s + 2hy + É0 2 y = 0 (16.6)
Gdy s jest znaną funkcją czasu, równanie sprowadzamy do postaci:
&& & s(t)
y + 2hy + É0 2 y = -&& (16.7)
Załóżmy że obiekt badany, a wraz z nim podstawa czujnika porusza się
ruchem harmonicznym o postaci:
s = H sin½t (16.8)
Równanie ruchu przyjmuje postać:
2
&& &
y + 2hy + É0 2 y = H½ sin½t (16.9)
Drgania wymuszone , ustalone badanego ciała mają postać:
y(t) = BAH sin(½t + ´ ) (16.10)
Wychylenie względne ciała o masie m , czyli wskazania czujnika są
więc proporcjonalne do przemieszczeń obiektu badanego i przesunięte
´
w czasie o czas Ä = .
½
y(t) = BAH sin(t +Ä ) (16.11)
Współczynnik BA nazywamy współczynnikiem czułości wibrometru
i wynosi on:
Strona 140
140
140
140
REJESTRACJA DRGAC
2
½
É0 2
BA = (16.12)
2
2 2
ëÅ‚ öÅ‚
h2½
ìÅ‚1- ½ ÷Å‚
+ 4
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
Rysunek 16.4 Przebiegi współczynnika czułości wibrometru
Jak widać przy częstościach drgań mierzonych, wyższych kilkukrotnie
od częstości drgań własnych czujnika, współczynnik czułości BA jest
bliski jedności, a więc czujnik mierzy przemieszczenie badanego
obiektu. Wniosek, iż czujniki przemieszczeń powinny być tak skonstruo-
wane, aby jego częstości drgań własnych była niska. Jednocześnie pro-
wadzi to do utrudnień gdyż czujniki o niskich częstościach muszą być
duże i ciężkie.
Te same przyrządy używane są do pomiarów przyspieszeń. Rozwiązanie
równania możemy przedstawić w następującej formie:
BP 2 BP
y(t) = H½ sin(½t + ´ ) = (-&&
s(t +Ä )) (16.13)
É0 2 É0 2
Tak więc wychylenie względne ciała i wskazania czujnika są proporcjo-
nalne do przyspieszenia obiektu drgającego. Współczynnik proporcjo-
nalności BP nazywamy czułością przyspieszeniomierza.
Strona 141
141
141
141
ROZDZIAA 16
1
BP = (16.14)
2
2 2
ëÅ‚ öÅ‚
h2½
ìÅ‚1- ½ ÷Å‚
+ 4
ìÅ‚ ÷Å‚
É0 2 É0 4
íÅ‚ Å‚Å‚
Rysunek 16.5 Charakterystyka przyspieszeniomierza
W rozpatrywanym przypadku czujnik może być nastrojony inaczej.
Częstość drgań własnych czujnika powinna być dużo większa od czę-
stości drgań mierzonych., gdyż tylko w tedy wartość współczynnika
czułości przyspieszeniomierza BP jest równa jedności, a więc czujnik
mierzy przyspieszenie.
Strona 142
142
142
142
`
Drgania swobodne
układu liniowego
o dwóch stopniach
17
swobody
bez tłumienia
ROZDZIAA 17
Rozpatrzmy układ jak na rysunku 17.1.
Rysunek 17.1 Analizowany układ
Składa się on z dwóch ciał o masach m1 i m2 połączonych elementami
sprężystymi o współczynnikach sztywności k1 i k2 . Przemieszczenie
poszczególnych ciał opisane jest dwoma współrzędnymi x1 i x2 . Ko-
rzystając z równania Lagrange`a II rodzaju ułóżmy równania ruchu.
Energia kinetyczna układu EK :
1 1
& &
EK = m1x12 + m2 x2 2 (17.1)
2 2
Energia potencjalna układu EP :
Strona 144
144
144
144
DRGANIA SWOBODNE UKAADU LINIOWEGO O DWÓCH STOPNIACH SWOBODY
BEZ TAUMIENIA
1 1
2
EP = k1x12 + k (x2 - x1) =
2
2 2
(17.2)
1 1 1
k1x12 + k x2 2 - k2 x1x2 + k x12
2 2
2 2 2
Ponieważ analizowany układ jest autonomiczny, równanie Lagrange`a II
rodzaju przyjmujemy w postaci:
ëÅ‚ öÅ‚
d "EK "EP
ìÅ‚ ÷Å‚
+ = 0 (17.3)
ìÅ‚
&
dt "xi ÷Å‚ "xi
íÅ‚ Å‚Å‚
Poszczególne człony równania (17.3) mają postać:
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
d "EK d "EK
ìÅ‚ ÷Å‚ && ìÅ‚ ÷Å‚ &&2
= m1x1; = m2 x (17.4)
ìÅ‚ ìÅ‚
& &
dt "x1 ÷Å‚ dt "x2 ÷Å‚
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
"EP
= k1x1 - k2 x2 + k2 x1 = k1x1 + k2(x1 - x2 ) (17.5)
"x1
"EP
= k2 x2 - k2 x1 = k2(x2 - x1) (17.6)
"x2
Zatem po podstawieniu do równania Lagrange`a w przyjętej postaci,
uwzględniając indeksy, otrzymujemy:
&&
m1x1 + k1x1 + k2(x1 - x2 )= 0
(17.7)
&&2
m2 x + k2(x2 - x1)= 0
Rozwiązania równań (17.7) przewidujemy w postaci:
x1 = a1 sin Ét ; x2 = a2 sin Ét
& &
x1 = a1É cosÉt ; x2 = a2É cosÉt (17.8)
2
&&1 -a1É sin Ét ; x = -a2É sin Ét
&&2 2
x =
Uwzględniając (17.8) w równaniach ruchu (17.7) otrzymujemy:
2
- m1a1É sin Ét + k1a1 sin Ét + k2a1 sin Ét - k2a2 sin Ét = 0
(17.9)
2
- m2a2É sin Ét + k2a2 sin Ét - k2a1 sin Ét = 0
Strona 145
145
145
145
ROZDZIAA 17
Po uporzÄ…dkowaniu otrzymujemy:
2
(a1(- m1É + k1 + k2)- k2a2)sin Ét = 0
(17.10)
2
(- k2a1 + a2(- m2É + k2))sin Ét = 0
Aby układ równań (17.10) zerował się dla każdego czasu t , muszą być
spełnione warunki:
2
a1(- m1É + k1 + k2)- k2a2 = 0
(17.11)
2
- k2a1 + a2(- m2É + k2)= 0
Aby ostatnie równania miały niezerowe rozwiązania musi być spełniony
warunek:
2
- m1É + k1 + k2 - k2
= 0 (17.12)
2
- k2 - m2É + k2
RozwiÄ…zujÄ…c wyznacznik (17.12) otrzymujemy:
4 2 2 2
m1m2É - k2m1É - k1m2É + k1k2 - k2m2É + k2 2 - k2 2 = 0 (17.13)
Po uproszczeniu i uporządkowaniu otrzymujemy równanie częstości:
4 2
m1m2É - (k2m1 + k1m2 + k2m2 )É + k1k2 = 0 (17.14)
Rozwiązując równanie (17.14) otrzymujemy dwie częstości drgań
własnych analizowanego układu.
2
" = (k2m1 + k1m2 + k2m2 ) - 4m1m2k1k2 (17.15)
2
É1,2 = (k2m1 + k1m2 + k2m2 )m (k2m1 + k1m2 + k2m2 ) - 4m1m2k1k2
Zatem układ będzie drgał w ogólnym przypadku z dwoma częstościami
É1 i É2 zależnymi od jego parametrów.
Przykład 1.
Wyznaczyć częstości drgań własnych dla układu jak na rysunku 17.2.
Strona 146
146
146
146
DRGANIA SWOBODNE UKAADU LINIOWEGO O DWÓCH STOPNIACH SWOBODY
BEZ TAUMIENIA
Rysunek 17.2 Analizowany układ drgający
WspółrzÄ™dne skrÄ™cenia waÅ‚u poprzez krążki odpowiednio Õ1 i Õ2 .
Badany układ składa się z dwóch krążków o momentach bezwładności
I1 i I2 połączonych ze sobą wałem o sztywności k na skręcanie.
Ułóżmy równania ruchu układu:
1 1
& &
EK = I Õ12 + I Õ2 2
1 2
2 2
1 1 1
2
EP = k(Õ2 - Õ1) = kÕ2 2 - kÕ1Õ2 + kÕ12
2 2 2
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
d "EK d "EK
ìÅ‚ ÷Å‚ && ìÅ‚ ÷Å‚ &&
= I1Õ1 ; = I Õ2
ìÅ‚ ìÅ‚
& &
dt "Õ1 ÷Å‚ dt "Õ2 ÷Å‚ 2
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
"EP
= -kÕ2 + kÕ1 = k(Õ1 - Õ2 )
"Õ1
"EP
= kÕ2 - kÕ1 = k(Õ2 - Õ1)
"Õ2
&&
I1Õ1 + k(Õ1 - Õ2 )= 0
&&
I2Õ2 + k(Õ2 - Õ1) = 0
Õ1 = Ä…1 sin Ét ; Õ2 = Ä…2 sin Ét
Strona 147
147
147
147
ROZDZIAA 17
& &
Õ1 = Ä…1É cosÉt ; Õ2 = Ä… É cosÉt
2
2 2
&& &&
Õ1 = -Ä…1É sin Ét ; Õ2 = -Ä…2É sin Ét
Podstawiając do równania ruchu:
2
- I1Ä…1É sin Ét + kÄ…1 sin Ét - kÄ… sin Ét = 0
2
2
- I Ä…2É sin Ét + kÄ… sin Ét - kÄ…1 sin Ét = 0
2 2
Po uproszczeniu.
2
[Ä…1(- I1É + k)- kÄ… ]sin Ét = 0
2
2
[- kÄ…1 + (- I2É + k)Ä… ]sin Ét = 0
2
Aby równania się zerowały dla każdego czasu t wyrażenia w nawiasach
kwadratowych muszą być równe zeru.
2
Ä…1(- I1É + k)- kÄ… = 0
2
2
- kÄ…1 + (- I É + k)Ä… = 0
2 2
Aby istniały niezerowe rozwiązania na ą1 i ą musi być spełniony
2
warunek:
2
- I1É + k - k
= 0
2
- k - I É + k
2
RozwiÄ…zujÄ…c wyznacznik otrzymujemy:
4 2 2 2 2
I1I2É - kI1É - kI2É + k - k = 0
Ostatecznie:
4 2
I1I2É - (kI1 + kI2 )É = 0
Równanie to można przedstawić w postaci:
2 2
É [I1I2É - k(I1 + I2 )]= 0
É1 = 0
Strona 148
148
148
148
DRGANIA SWOBODNE UKAADU LINIOWEGO O DWÓCH STOPNIACH SWOBODY
BEZ TAUMIENIA
k(I1 + I2 )
É2 =
I1I
2
CzÄ™stość É1 odpowiada toczeniu siÄ™ rozpatrywanego ukÅ‚adu, natomiast
czÄ™stość É2 wynika z wzajemnych drgaÅ„ wzglÄ™dem siebie krążków
o momentach bezwładności I1 i I .
2
Przykład 2.
Dla układu z przykładu 1 znalezć współrzędne przekroju pręta o sztyw-
ności k , który nie będzie podlegał skręcaniu. Pręt okrągły.
Rysunek 17.2 Przypadek analizowany w przykładzie 2
Częstość drgań krążka o momencie bezwładności I1 wynosi:
ka
Éa =
I1
gdzie: ka - sztywność na skręcanie odcinka a pręta.
Częstość drgań krążka o momencie bezwładności I wynosi:
2
kb
Éb =
I
2
gdzie: kb - sztywność na skręcanie odcinka b pręta.
Sztywności odcinków pręta a i b wynoszą odpowiednio:
GJ GJ0
0
ka = ; kb =
a b
Strona 149
149
149
149
ROZDZIAA 17
gdzie:
G - moduł Kirchoffa materiału pręta o sztywności k ,
J0 - osiowy moment bezwładności.
Zatem:
GJ0 GJ0
Éa = ; Éb =
aI1 bI2
CzÄ™stoÅ›ci Éa i Éb muszÄ… być sobie równe.
Éa = Éb
czyli:
GJ0 GJ
0
=
aI1 bI
2
Zatem:
1 1
=
aI1 bI2
Tak więc:
bI2 = aI1
Po uporzÄ…dkowaniu:
a I
2
=
b I1
Szukany przekrój nie podlegający skręcaniu zależny jest od stosunku
momentów bezwładności I1 i I krążków.
2
Strona 150
150
150
150
`
Drgania wymuszane
układów o dwóch
stopniach swobody,
18
tłumienie dynamiczne
ROZDZIAA 18
Rozpatrzmy liniowy układ o dwóch stopniach swobody, pozbawiony
tłumienia, jak na rysunku 18.1.
Rysunek 18.1 Rozpatrywany układ o dwóch stopniach swobody
Ciało o masie m1 zawieszone jest na sprężynie o sztywności k1 . Z tym
ciałem za pośrednictwem sprężyny o sztywności k2 połączone jest
drugie ciało o masie m2 . Na ciało o masie m1 działa harmoniczna siła
wymuszajÄ…ca w postaci f (t) = P sin½t . UkÅ‚ad wykonuje drgania
podłużne.
Energia kinetyczna układu EK :
1 1
& &
EK = m1x12 + m2 x2 2 (18.1)
2 2
Energia potencjalna układu EP :
Strona 152
152
152
152
DRGANIA WYMUSZONE UKAADÓW O DWÓCH STOPNIACH SWOBODY, TAUMIENIE DYNAMICZNE
1 1
2
EP = k1x12 + k (x1 - x2 ) =
2
2 2
(18.2)
1 1 1
k1x12 + k x12 - k2 x1x2 + k x2 2
2 2
2 2 2
Korzystając z równania Lagrange`a II rodzaju możemy napisać:
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
d "EK d "EK
ìÅ‚ ÷Å‚ && ìÅ‚ ÷Å‚ &&2
= m1x1; = m2 x (18.3)
ìÅ‚ ìÅ‚
& &
dt "x1 ÷Å‚ dt "x2 ÷Å‚
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
"EP
= k1x1 + k2 x1 - k2 x2 = k1x1 + k2(x1 - x2 ) (18.4)
"x1
"EP
= k2 x2 - k2 x1 = k2(x2 - x1) (18.5)
"x2
Podstawiając do równania Lagrange`a i uwzględniając siłę wymuszającą,
otrzymujemy układ równań, opisujących drgania badanego układu.
&&
m1x1 + k1x1 + k2(x1 - x2 )= P sin½t
(18.6)
&&2
m2 x + k2(x2 - x1)= 0
Oczywiście x1 i x2 współrzędne liczone od położenia równowagi
statycznej.
Rozwiązania szczególne drgań wymuszonych przewidujemy w postaci:
x1 = A1 sin½t ; x2 = A2 sin½t
& &
x1 = A1½ cos½t ; x2 = A2½ cos½t (18.7)
2
&&1 -A1½ sin½t ; x = -A2½ sin½t
&&2 2
x =
Podstawiając do układu równań (18.6) otrzymujemy:
2
- m1 A1½ sin½t + k1A1 sin½t + k2 A1 sin½t - k2 A2 sin½t = H sin½t
(18.8)
2
- m2 A2½ sin½t + k2 A2 sin½t - k2 A1 sin½t = 0
Na podstawie (17.11) możemy zapisać:
Strona 153
153
153
153
ROZDZIAA 18
2
(- m1½ + k1 + k2)A1 - k2 A2 = H
(18.9)
2
- k2 A1 + (- m2½ + k2)A2 = 0
Z równania drugiego wyznaczmy A2 :
k2 A1
A2 = (18.10)
2
(k2 - m2½ )
PodstawiajÄ…c do pierwszego otrzymujemy:
k2 2 A1
2
(k1 + k2 - m1½ )A1 - (18.11)
2
(k2 - m2½ )= P
stÄ…d A1 :
2
P(k2 - m2½ )
A1 = (18.12)
2 2
(k1 + k2 - m1½ )(k2 - m2½ )- k2 2
Podstawiając do wyrażenia na A2 otrzymujemy:
2
k2 P(k2 - m2½ )
A1 =
2 2 2
(k2 - m2½ )[(k1 + k2 - m1½ )(k2 - m2½ )- k2 2]=
(18.13)
k2 H
=
2 2
(k1 + k2 - m1½ )(k2 - m2½ )- k2 2
Wyznaczmy współczynniki uwielokrotnienia amplitudy µ1 i µ2 :
A1 A2 P
µ1 = ; µ2 = ; st = (18.14)
st st k1
Podstawiając wyrażenie na A1, otrzymujemy:
Strona 154
154
154
154
DRGANIA WYMUSZONE UKAADÓW O DWÓCH STOPNIACH SWOBODY, TAUMIENIE DYNAMICZNE
2
A1 P(k2 - m2½ ) k1
µ1 = = =
2
st + k2 - m1½ 2
(k1 )(k2 - m2½ )- k2 2 H
2
ëÅ‚ öÅ‚
½
÷Å‚
k1k2 ìÅ‚1 -
ìÅ‚ ÷Å‚
É2 2
íÅ‚ Å‚Å‚
= = (18.15)
ëÅ‚
k2 2 öÅ‚ ëÅ‚ 2 öÅ‚ 2
½
÷Å‚k2 ìÅ‚1 - ½ ÷Å‚
k1ìÅ‚1 + - - k2
ìÅ‚
k1 É12 ÷Å‚ ìÅ‚ É2 2 ÷Å‚
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
2
ëÅ‚ öÅ‚
ìÅ‚1 - ½ ÷Å‚
ìÅ‚ ÷Å‚
É2 2
íÅ‚ Å‚Å‚
=
2 2
ëÅ‚ öÅ‚ëÅ‚ öÅ‚
k2
ìÅ‚1 + k2 ½ ÷Å‚ìÅ‚1 - ½ ÷Å‚
- -
ìÅ‚
k1 É12 ÷Å‚ìÅ‚ É2 2 ÷Å‚ k1
íÅ‚ Å‚Å‚íÅ‚ Å‚Å‚
Podstawiając wyrażenie na A2 , otrzymujemy:
A2 k2 P k1
µ2 = = =
2
st + k2 - m1½ 2
(k1 )(k2 - m2½ )- k2 2 P
k1k2
= = (18.16)
ëÅ‚
k2 2 öÅ‚ ëÅ‚ 2 öÅ‚ 2
½
÷Å‚k2 ìÅ‚1 - ½ ÷Å‚
k1ìÅ‚1 + - - k2
ìÅ‚
k1 É12 ÷Å‚ ìÅ‚ É2 2 ÷Å‚
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
1
=
2 2
ëÅ‚ öÅ‚ëÅ‚ öÅ‚
k2
ìÅ‚1 + k2 ½ ÷Å‚ìÅ‚1 - ½ ÷Å‚
- -
ìÅ‚
k1 É12 ÷Å‚ìÅ‚ É2 2 ÷Å‚ k1
íÅ‚ Å‚Å‚íÅ‚ Å‚Å‚
Przebieg współczynnika µ1 przedstawia rysunek 18.2.
Strona 155
155
155
155
ROZDZIAA 18
½
Rysunek 18.2 Przebiegi współczynnika µ1 w funkcji parametru
É1
W przypadku gdy:
½ = É1 = É2 (18.17)
Wtedy µ1 = 0 .
Następuje całkowity zanik drgań ciała o masie m1 . Takie zjawisko nazy-
wamy tłumieniem dynamicznym i jest ono bardzo często wykorzystywa-
ne w technice. Gdy chcemy aby ciało o masie m1 pozostawało w spo-
czynku, (mimo iż przyłożona jest doń siła wymuszająca f (t) ), to
dodajemy do niego dodatkowe ciało o masie m2 na sprężynie o sztyw-
noÅ›ci k2 , tak dobranej aby byÅ‚ speÅ‚niony warunek ½ = É1 = É2 . Tak
dobrane ciało i sprężyna nazywa się tłumikiem dynamicznym.
Tłumik wykonuje ruch określony zależnością:
ëÅ‚ öÅ‚
P k1 P
ìÅ‚- ÷Å‚ - sin½t (18.18)
x2 = A2 sin½t = st µ2 sin½t = sin½t =
k1 ìÅ‚ k2 ÷Å‚ k2
íÅ‚ Å‚Å‚
Siła przenoszona przez sprężynę:
P
R = k2 x2 = -k2 sin½t = -P sin½t (18.19)
k2
Strona 156
156
156
156
DRGANIA WYMUSZONE UKAADÓW O DWÓCH STOPNIACH SWOBODY, TAUMIENIE DYNAMICZNE
Tak więc reakcja sprężyny jest równa co do wartości sile wymuszającej
i ma zwrot przeciwny. Równoważenie się tych sił powoduje iż ciało
o masie m1 do którego przyłożona jest siła wymuszająca pozostaje
w spoczynku.
Taki tłumik dynamiczny posiada podstawową wadę. Jak widać z prze-
biegu współczynnika µ1 drgania sÄ… wytÅ‚umione do zera praktycznie tyl-
ko dla jednej częstości siły wymuszającej. Przy reaktywnych krzywych
µ1 najmniejsza zmiana czÄ™stoÅ›ci siÅ‚y wymuszajÄ…cej powoduje już
bardzo duży wzrost współczynnika uwielokrotnienia, a co za tym idzie
znaczny wzrost amplitudy ciała o masie m1 . Taki tłumik jest więc mało
efektywny.
Jeżeli wprowadzimy tłumik pomiędzy ciałem głównym o masie m1
a ciałem m2 , czyli w układzie eliminatora pojawi się tłumik jak na
rysunku 18.3.
Rysunek 18.3 Dynamiczny eliminator drgań z dodatkowym tłumikiem
W równaniach ruchu pojawi się dodatkowy człon z tłumieniem:
&& & &
m1x1 + k1x1 + k2(x1 - x2 )+ c2(x1 - x2 )= P sin½t
(18.20)
&&2 & &
m2 x + k2(x2 - x1)+ c2(x2 - x1)= 0
Współczynnik uwielokrotnienia amplitudy ma w tym przypadku
przebieg.
Strona 157
157
157
157
ROZDZIAA 18
Rysunek 18.4 Przebiegi współczynnika µ1 ze zÅ‚agodzeniem
wpływu rezonansów
Widać wyraznie że następuje złagodzenie wpływu rezonansów. Działa-
nie tÅ‚umika jest jednak osÅ‚abione gdyż w punkcie ½ = É1 = É2 współ-
czynnik uwielokrotnienia jest większy od zera.
Można tak dobrać parametry układu aby styczne do krzywych rezonan-
sowych w punktach przejścia P i Q były poziome.
Rysunek 18.5 Przebiegi współczynnika µ1 z poziomymi stycznymi
w punktach P i Q
Strona 158
158
158
158
DRGANIA WYMUSZONE UKAADÓW O DWÓCH STOPNIACH SWOBODY, TAUMIENIE DYNAMICZNE
Przykład 1.
Podać warunki, stanowiące podstawę doboru amortyzatora drgań skręt-
nych krążka o masie m1 , wymuszanych momentem sinusoidalnie zmien-
nym M (t) = M sin½t . Schemat ukÅ‚adu przedstawiono na rysunku 18.6.
O
Rysunek 18.6 Analizowany układ w przykładzie 1
Krążek o masie m1 , obciążony momentem sił zewnętrznych M (t) , oraz
krążek o masie m2 amortyzatora tworzą układ o dwóch stopniach
swobody. Przemieszczenia krążków opisywać będą odpowiednio kąty
skrÄ™cenia Õ1 i Õ2 .Równania ruchu ukÅ‚adu majÄ… postać:
&&
I1Õ1 - k2Õ2 + (k1 + k2 )Õ1 = M sin½t
0
&&
I Õ2 + k2Õ2 - k2Õ1 = 0
2
Lub:
&&
I1Õ1 - k2Õ2 + kÕ1 = M sin½t
0
&&
I Õ2 + k2Õ2 - k2Õ1 = 0
2
gdzie: k = k1 + k2
Przewidujemy rozwiązania równań w postaci:
Õ1 = Ä…1 sin½t i Õ2 = Ä… sin½t
2
gdzie:
ą1 ; ą - odpowiednio amplitudalne wychylenia krążków o ma-
2
sach m1 i m2 .
Po wprowadzeniu tych funkcji i ich drugich pochodnych do równań
ruchu otrzymujemy po uproszczeniu:
Strona 159
159
159
159
ROZDZIAA 18
2 2
(k - I1½ )Ä…1 - k2Ä… = M ; (k2 - I ½ )Ä… - k2Ä…1 = 0
2 O 2 2
Z powyższego równania wyznaczamy amplitudalne wychylenia krążków
Ä…1 ; Ä…2 :
2
ëÅ‚ öÅ‚
ìÅ‚1 - ëÅ‚ ½ öÅ‚ ÷Å‚Õ
ìÅ‚ ÷Å‚
ìÅ‚ ÷Å‚
ìÅ‚
É2 ÷Å‚ st
íÅ‚ Å‚Å‚
íÅ‚ Å‚Å‚
Ä…1 =
2 2
îÅ‚ Å‚Å‚îÅ‚ Å‚Å‚
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚ ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
k2 ½ ½ k2
ïÅ‚ ìÅ‚ ÷Å‚ - ìÅ‚ ÷Å‚ śłïÅ‚ - ìÅ‚ ÷Å‚ śł - ìÅ‚ ÷Å‚
1 + 1
ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚
k1 É1 ûÅ‚ðÅ‚ É2 ûÅ‚ k1
ïÅ‚ śłïÅ‚ śł
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
ðÅ‚
Õst
Ä… =
2
2 2
îÅ‚ Å‚Å‚îÅ‚ Å‚Å‚
ëÅ‚ öÅ‚ ëÅ‚ öÅ‚ ëÅ‚ öÅ‚ ëÅ‚ öÅ‚
k2 ½ ½ k2
ïÅ‚ ìÅ‚ ÷Å‚ - ìÅ‚ ÷Å‚ śłïÅ‚ - ìÅ‚ ÷Å‚ śł - ìÅ‚ ÷Å‚
1 + 1
ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚ ìÅ‚ ÷Å‚
k1 É1 ûÅ‚ðÅ‚ É2 ûÅ‚ k1
ïÅ‚ śłïÅ‚ śł
íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚
ðÅ‚
gdzie:
k1
É1 = - czÄ™stość drgaÅ„ wÅ‚asnych skrÄ™tnych krążka m1 ,
I1
k2
É2 = - czÄ™stość drgaÅ„ wÅ‚asnych skrÄ™tnych krążka m2 ,
I2
M
O
Õst = - statyczny kÄ…t skrÄ™cenia waÅ‚u pierwszego pod
k1
wpływem amplitudalnej wartości momentu
wymuszajÄ…cego.
Dobór amortyzatora m2 opierać się będzie na następujących dwóch
warunkach:
k2
½ = É2 =
I2
Oraz:
k2
½ = É2 =
I2
Strona 160
160
160
160
DRGANIA WYMUSZONE UKAADÓW O DWÓCH STOPNIACH SWOBODY, TAUMIENIE DYNAMICZNE
Stosując dynamiczny eliminator drgań można także tłumić drgania
układów o większej liczbie stopni swobody niż dwa. Przypadek taki
obrazuje przykład 2.
Przykład 2.
Podać warunki, stanowiące podstawę doboru amortyzatora drgań skręt-
nych krążka o masie m3 , drgań skrętnych układu z dwoma krążkami m1
i m2 wymuszanych momentem przyłożonym do krążka m2 sinusoidal-
nie zmiennym M (t) = M sin½t . Schemat ukÅ‚adu przedstawiono na
O
rysunku 18.7.
Rysunek 18.7 Analizowany układ w przykładzie 2
Oba krążki m1 i m2 oraz amortyzator m3 tworzą układ o trzech stop-
niach swobody tego samego typu jak układ analizowany w przykła-
dzie 2. Wobec tego przystosujemy tutaj wyjściowe równania różniczko-
we z poprzedniego przykładu, uzupełniając drugie z nich czynnikiem,
wyrażającym moment M (t) sił zewnętrznych.
W ten sposób będziemy mieli:
&&
I1Õ1 - k2(Õ2 - Õ1)= 0
&&
I Õ2 + k2(Õ2 - Õ1)- k3(Õ3 - Õ2 )= M sin½t
2 0
&&
I3Õ3 + k3(Õ3 - Õ2 )= 0
Strona 161
161
161
161
ROZDZIAA 18
Doświadczenie 1.
Celem ćwiczenia jest zapoznanie się z zasadą działania dynamicznych
eliminatorów drgań oraz porównanie krzywych rezonansowych dla bez-
władnościowo wymuszanych drgań układu głównego: bez eliminatora
oraz z tłumionym dynamicznym eliminatorem drgań.
Schemat stanowiska wykorzystywanego w doświadczeniu przedstawia
rysunek 18.8.
Rysunek 18.8 Schemat stanowiska
Stanowisko zawiera następujące elementy:
1 - Układ główny o masie M podwieszony na płaskich sprężynach (1).
2 - Eliminator o masie m podwieszony na sprężynach (3) i zawierający
tłumik (5).
6 - Wzbudnik drgań.
W ćwiczeniu należy zarejestrować przebiegi drgań układu głównego (1)
ze zblokowanym eliminatorem przy płynnej zmianie prędkości obroto-
wej wibratora oraz przy różnych ustalonych prędkościach obrotowych
wibratora. Wykonać analogiczne pomiary dla układu z odblokowanym
eliminatorem dynamicznym. Wyznaczyć rodziny krzywych rezonanso-
wych dla obydwu przypadków przedstawione jako teoretyczne na rysun-
kach 18.2 i 18.4.
Strona 162
162
162
162
`
Literatura
19
ROZDZIAA 6
1. Osiński. Z., Teoria drgań., PWN, Warszawa, 1980.
2. Piszczek K., Walczak J., Drgania w budowie maszyn, PWN,
Warszawa, 1972.
3. Kaliski S. i zespół, Drgania i fale, PWN, Warszawa, 1966.
4. Nizioł J., Podstawy drgań w maszynach., Politechnika Krakow-
ska, Kraków, 1996.
5. Giergiel J., Tłumienie drgań mechanicznych, PWN, Warszawa,
1990.
6. Osiński Z., Tłumienie drgań mechanicznych, PWN, Warszawa,
1976.
7. Kamiński E., Podstawy dynamiki maszyn, Wydawnictwo Poli-
techniki Warszawskiej , Warszawa, 1980.
8. Giergiel J., Drgania mechaniczne układów dyskretnych, Wydaw-
nictwo Politechniki Rzeszowskiej, Rzeszów, 2004.
9. Praca zbiorowa, Teoria drgań zbiór zadań, Wydawnictwo Po-
litechniki Warszawskiej, Warszawa, 1985.
10. Nizioł J., Metodyka rozwiązywania zadań z mechaniki, WT,
Warszawa, 2002.
11. Praca zbiorowa, Drgania mechaniczne, (laboratorium), Oficyna
wydawnicza PW, Warszawa, 1995.
Strona 164
164
164
164
Wyszukiwarka
Podobne podstrony:
dziurdz drgania mechaniczne i halasDrgania mechaniczne05 Drgania mechaniczneHałas i drgania mechaniczneOchrona przed hałasem i drganiami mechanicznymi w środowisku pracyjach,fizyka środowiska pracy, drgania mechanicznelekcja 4 drgania mechan wibracje hałas w srod pracyDrgania mechaniczne (Wibracje)21 mechanika budowli wykład 21 drgania wymuszone nietlumionewykl mechanika budowli opis ruchu drgania wlasne tlumione20 mechanika budowli wykład 20 drgania pretow pryzmatycznych?więcej podobnych podstron