Wahadło
Wahadło – ciało zawieszone w jednorodnym polu grawitacyjnym w taki sposób, że może wykonywać drgania wokół poziomej osi nieprzechodzącej przez środek ciężkości zawieszonego ciała.
W mechanice wyróżnia się dwa podstawowe modele fizyczne wahadeł[1]:
- matematyczne (proste) – opisujące wahadło jako punkt materialny, zawieszony na nieważkiej nici,
- fizyczne – opisujące wahadło jako bryłę sztywną.
Ważną cechą wahadeł fizycznego i matematycznego jest niemal pełna niezależność ich okresu drgań od amplitudy, przy założeniu że amplituda drgań jest mała[a]. Własność ta, zwana izochronizmem drgań, została odkryta około 1602 roku przez Galileusza, który używał wahadła do pomiaru czasu. Zainspirowany tą zasadą Christiaan Huygens zbudował w 1656 roku pierwszy zegar wahadłowy[2]. Zegary wahadłowe były najdokładniejszymi urządzeniami do pomiaru czasu aż do skonstruowania w latach 30. XX wieku zegarów kwarcowych.
Ogólnie wahadło jest oscylatorem anharmonicznym, jego okres drgań i inne parametry zależą od amplitudy drgań. Rozwiązanie równania ruchu wahadła upraszcza się, jeśli założy się drgania o małych amplitudach. Jednak znalezienie rozwiązania dla dowolnie dużych amplitud wymaga użycia funkcji nieelementarnej, tzw. funkcji eliptycznej sn Jacobiego. Opis ruchu wahadła komplikuje się jeszcze bardziej, gdy trzeba uwzględnić efekt tłumienia ruchu wahadła. Wtedy efektywne rozwiązanie równania drgań dla dowolnych amplitud wymaga zastosowania metod numerycznych.
Historia
W XVII w. Galileusz w czasach swej młodości odkrył izochronizm wahadła, oraz że okres drgań zależy jako pierwiastek kwadratowy długości wahadła. Wykorzystywał wahadło do odmierzania czasu. W 1644 Marin Mersenne wyznaczył długość wahadła sekundowego (o okresie 2 sekund). Współczesny mu Stanisław Pudłowski proponował oprzeć miarę długości na zjawisku wahadła. W 1657 Huygens przedstawił i opatentował zegar wahadłowy; wynalazek szybko rozprzestrzenia się. W 1673 Huygens przedstawił teorię wahadła, w tym zależność okresu drgań wahadła od miejsca zawieszenia wahadła fizycznego (p. Jean Richer). W 1687 Isaac Newton w pracy Principia zauważył, że przyspieszenie ziemskie można wyrazić jako długość wahadła sekundowego. W 1737 Pierre Bouguer wykonał pomiary przyspieszenia ziemskiego między innymi w Andach. Zauważył, że do pomiaru przyspieszenia ziemskiego wygodniej jest określanie okresu wahania, a nie długości wahadła sekundowego (ta wcześniejszą metoda wymagała zmieniania długości wahadła tak, by otrzymać okres drgań równy 1 s). Używając wahadeł, porównał gęstość Ziemi z gęstością Kordylierów. Od 1735 Charles Marie de La Condamine prowadził eksperymenty z wahadłami, dopracowując i wykonując pomiary przyspieszenia ziemskiego w różnych miejscach. W 1792 we Francji długość wahadła sekundowego była proponowana jako jednostka długości. Pomysł nie został przyjęty w metrycznym systemie miar[3]. W latach 1825/27 Bessel udoskonalił układ pomiarowy oraz wprowadził układ optyczny do obserwacji ruchu wahadła. W 1817 Henry Kater skonstruował wahadło rewersyjne, dając impuls do dokładnych i bezwzględnych pomiarów przyspieszenia ziemskiego. W latach 1827–1840 Francis Baily skonstruował różne wahadła, w tym wahadło poruszające się w próżni[3].
Wahadło matematyczne
Definicja
Wahadłem matematycznym nazywa się punkt materialny poruszający się po okręgu w płaszczyźnie pionowej w jednorodnym polu grawitacyjnym[1]. Wahadło rzeczywiste, złożone z ciała zawieszonego na nici, może być traktowane jako wahadło matematyczne, jeżeli spełnione są następujące założenia[4]:
- rozmiary ciała są niewielkie w porównaniu z długością nici,
- nić jest nieważka,
- nić jest nierozciągliwa,
- wahadłu nadano prędkość początkową tak, że drga w płaszczyźnie pionowej,
- na ciało działają jedynie siła ciężkości oraz siła reakcji nici (pomijalne są inne siły, np. siła oporów ruchu ze strony powietrza).
Wahadło matematyczne stanowi szczególny przypadek wahadła fizycznego (patrz niżej)[4].
Równanie ruchu wahadła matematycznego
Równanie ruchu wahadła określa wzór[1][4]:
gdzie:
- – kąt odchylenia wahadła od pionu w chwili przy czym kąt ten przyjmują wartości dodatnie np. dla odchyleń w prawo, a ujemne dla odchyleń w lewo,
- – przyspieszenie ziemskie,
- – długość nici.
Wahadło jest odchylone od pionu o kąt θ. Na ciało działa siła ciężkości oraz siła naprężenia nici. Zawieszenie wahadła wymusza ruch po łuku w płaszczyźnie pionowej. Siłę ciężkości wahadła rozkłada się na styczną do kierunku ruchu i prostopadłą do niego. Składowa prostopadła do kierunku ruchu nie wpływa na wartość prędkości, zmienia jedynie jej kierunek. Składowa równoległa nadaje ciału przyspieszenie:
Przyspieszenie to powoduje przebycie drogi, która jest długością łuku:
Przybliżenie I równania ruchu wahadła - założenie małej amplitudy drgań
Równanie ruchu wahadła jest równaniem różniczkowym drugiego stopnia, nieliniowym. Dokładne rozwiązanie takiego równania jest niełatwe. Tu podano pierwszy sposób znalezienia rozwiązania przybliżonego, przy założeniu że wahadło wykonuje drgania o niewielkiej amplitudzie kątowej. Przy takim założeniu funkcję sinus można przybliżyć jej argumentem (wzór Maclaurina)[b][1]:
wówczas równanie ruchu wahadła upraszcza się do postaci
Powyższe równanie jest równaniem drgań harmonicznych. Rozwiązanie określa zależność kąta wahań od czasu i może być określone wzorem[5]:
gdzie:
- – amplituda drgań,
- – częstość kołowa małych drgań,
- – faza początkowa drgań.
Okres małych drgań jest związany z częstością wzorem
okres drgań wynosi[4]
Wynika stąd, że gdy wahadło wykonuje małe drgania, to:
- zależność kąta wychylenia wahadła od czasu jest funkcją harmoniczną,
- okres drgań nie zależy od amplitudy, a jedynie od długości wahadła i przyspieszenia grawitacyjnego.
Izochronizm wahadła odkrył już Galileusz w XVI w., własność ta jest słuszna dla małych drgań.
Przybliżenie II równania ruchu wahadła - dwa wyrazy szeregu Maclaurina
Zagadnienie ruchu wahadła można rozwiązać dla dużych amplitud drgań w sposób przybliżony, przybliżając funkcję sinus do dwóch początkowych wyrazów szeregu Maclaurina: . Wówczas równanie ruchu wahadła przyjmuje postać[6]:
Rozwiązanie tego równania ma postać (z dokładnością do wyrazów 3-go rzędu)
gdzie:
- – amplituda drgań o częstotliwości
- – amplituda drgań o częstotliwości
- - częstotliwość małych drgań
- - częstotliwość drgań wahadła o dużej amplitudzie;
Rozwiązanie powyższe wskazuje, iż wahadło matematyczne nie jest oscylatorem harmonicznym, lecz:
- ruch wahadła jest złożeniem dwóch drgań harmonicznych mających częstotliwości oraz i amplitudy odpowiednio równe oraz
- częstotliwość drgań zależy od amplitudy (nie występuje izochronizm drgań charakterystyczny dla małych amplitud)
- częstotliwość drgań jest mniejsza niż dla drgań o małej amplitudzie; częstotliwość ta maleje wraz ze wzrostem amplitudy
- amplituda wyższej harmonicznej zależy w trzeciej potędze od amplitudy
Z drugiej strony, dla dostatecznie małych wartości częstotliwość drgań zbliża się do wartości a amplituda wyższej harmonicznej staje się pomijalnie mała – otrzymuje się drganie harmoniczne o amplitudzie i częstotliwości która nie zależy od amplitudy.
Okres drgań wahadła dla dowolnie dużych amplitud
Zależność czasu drgań od kąta wychylenia wahadła
Rozwiązanie ogólnego równania ruchu wahadła dla dowolnej amplitudy można podać w postaci uwikłanej[7]:
Wykonując całkowanie w zakresie od 0 do otrzymuje się zależność czasu t od kąta wychylenia wahadła:
Całka w powyższym wzorze jest całką eliptyczną niezupełną pierwszego rodzaju.
Okres drgań o dowolnej amplitudzie
Przyjmując w powyższym wzorze wartość kąta w granicy górnej całki i mnożąc go przez 4 otrzyma się wzór na okres drgań wahadła o dowolnej amplitudzie drgań[8]:
gdzie K jest całką eliptyczną zupełną pierwszego rodzaju. Całki te są stablicowane.
Dla dużych amplitud wahań okres drgań zależy od amplitudy i rośnie wraz z jej wzrostem, co pokazuje wykres.
Rozwinięcie wzoru na okres drgań o dowolnej amplitudzie w szeregi
(1) Całkę eliptyczną w powyższym wzorze można rozwinąć w szereg Taylora względem , co prowadzi do wzoru[8][9]
(2) Dokładniejsze przybliżenie uzyskuje się rozwijając całkę eliptyczną w szereg Maclaurina względem (przy czym całka winna być najpierw wyrażona w bazie wielomianów Legendre'a)[10]
Gdy w powyższych wzorach pominie się wyrazy sumy poza pierwszym wyrazem, równym 1, to otrzyma się wzór na okres małych drgań wahadła (patrz wyżej). Kąt graniczny izochronizmu zależy od przyjętej dokładności; dla kąta 6° okres drgań jest o około 0,07% większy od minimalnego[11].
Wzór na okres drgań wahadła dla dowolnych kątów można wyprowadzić z zasady zachowania energii. Jeżeli ciało porusza się w dół od stanu spoczynku, to nabywa energię kinetyczną kosztem utraty energii potencjalnej grawitacji. Jeżeli nie ma strat energii, to powyższe dwie wielkości są sobie równe, czyli[12]
Stąd prędkość wahadła wynosi:
Ponieważ to z powyższych wzorów otrzymuje się prędkość kątową wahadła
Wysokość na jakiej znajduje się wahadło:
Zmiana wysokości jest różnicą wysokości dwóch położeń, to
Ostatecznie otrzymuje się[12]
albo
Okres wahań T otrzymuje się całkując powyższe równanie w granicach od 0 do i mnożąc całkę przez 4 (wahadło wraca do początku ruchu po 4 takich ruchach)[13]:
Całka występująca w powyższym wzorze jest całką eliptyczną. Aby przepisać ją do postaci Lagendre’a, której wartości są stablicowane, wyraża się θ w zależności od u, dokonując przekształceń, oraz podstawienia prowadzi do wzoru na okres drgań wahadła[13][14]
gdzie K jest zupełną całką eliptyczną pierwszego rodzaju zdefiniowaną jako
gdzie:
Całkę tę można rozwinąć w szereg[13]
Zostanie tutaj wyprowadzone równanie ruchu wahadła w oparciu o zasadę zachowania energii mechanicznej. Energia mechaniczna wahadła jest zachowana, gdyż zakłada się tutaj, że na układ działają jedynie siły zachowawcze, a pomija się opory ruchu. W ogólności zachowanie energii mechanicznej pozwala na wyprowadzenie równania ruchu dowolnego układu bez potrzeby odwoływania się do konkretnej postaci działających sił, co zazwyczaj prowadzi do prostszych obliczeń.
Różniczkując względem czasu wyprowadzone wyżej równanie (z użyciem reguły łańcuchowej), otrzyma się przyspieszenie kątowe
czyli
Stąd:
co jest tym samym równaniem które zostało wyprowadzone wcześniej z analizy sił.
Dokładne rozwiązanie równania ruchu wahadła dla dowolnych amplitud
Jawną zależności kąta wychylenia wahadła od czasu dla dowolnie dużych amplitud drgań uzyskuje się obliczając funkcję odwrotną do funkcji podaną w poprzednim rozdziale, przy czym wykorzystuje się funkcję sinus amplitudy Jacobiego[9].
(a) Jeżeli zada się warunki początkowe , (tj. w chwili wahadło w pionie oraz posiada zadaną prędkość kątową), to rozwiązanie ogólne wyraża się przez funkcję sinus amplitudy Jacobiego [15][9]:
- ,
gdzie .
(b) Jeżeli zada się warunki początkowe (tj. w chwili wahadło odchylone od pionu o maksymalny kąt ), to rozwiązanie ogólne wyraża się przez funkcję sinus amplitudy Jacobiego, ale trzeba dokonać przesunięcia na osi czasu [15]
- ,
gdzie:
- - okres drgań wahadła
Uwaga:
Zachodzi równość , jeżeli dla amplitudy drgań prędkość kątowa wahadła wynosi , gdy przechodzi przez położenie równowagi.
Generowanie wykresów drgań z użyciem funkcji sinus amplitudy Jacobiego
Zamieszczony tu program w języku Python generuje wykresy dla dowolnych amplitud ; program korzysta m.in. z biblioteki scipy.special do liczenia całek eliptycznych oraz funkcji sn amplitudy Jacobiego. Program można testować, korzystając np. z darmowego notatnika colab google online.
import numpy as np import matplotlib.pyplot as plt from scipy.special import ellipk, ellipj # Constants (can be changed) --------------------- theta_0_deg_list = [10, 133, 160, 179.99999] # List of amplitudes in degrees time_max = 3.2 # Time range in seconds time_steps = 1000 plt.figure(figsize=(10, 10)) # Size of a plot g = 9.81 # Acceleration due to gravity (m/s^2) L = 0.25 # Length of the pendulum (m) T_0 = 4 * np.sqrt(L / g) * ellipk(0) # Period for small amplitudes # Calculations --------------------- t = np.linspace(0, time_max, time_steps) # Time array for theta_0_deg in theta_0_deg_list: # Loop through each amplitude theta_0 theta_0 = np.radians(theta_0_deg) # Convert to radians k = np.sin(theta_0 / 2) # Elliptic modulus T = 4*np.sqrt(L/g) * ellipk(k**2)/T_0 # Period of the pendulum t_0 = 0 #T/4 # Time shift (can be changed) u = np.sqrt(g/L) * (t+t_0) # Time parameter for theta_0 sn, cn, dn, ph = ellipj(u, k**2) # The Jacobi elliptic functions theta = 2 * np.arcsin(k * sn) # Theta as a function of time # Plotting the result for theta_0 plt.plot(t/T_0, np.degrees(theta), label=r'$\theta_0 = {}^\circ$'.format(theta_0_deg)) # Plot customization plt.xlabel('time (s)', fontsize=20) plt.ylabel('θ(t) [${}^\circ$]', fontsize=20) plt.title('Nonlinear pendulum solutions for different amplitudes', fontsize=20) plt.grid(True) plt.legend(fontsize=18) plt.tick_params(axis='both', which='major', labelsize=18) plt.show()
Rozwinięcie w szereg Fouriera
Ogólne zależności
Znając zależność kąta wychylenia od czasu można rozwinąć ją w szereg Fouriera:
gdzie:
- - okres drgań wahadła (dla dowolnej amplitudy drgań),
- - podstawowa częstotliwość kołowa w rozwinięciu funkcji w szereg Fouriera
Szeregu Fouriera wyrażony za pomocą nomu eliptycznego
Szereg Fouriera ma postać[16][17]
gdzie:
- - amplituda drgań wahadła
- - moduł
- - całki eliptyczne zupełne 1. rzędu
- - tzw. nom eliptyczny
Definiując
można przybliżyć wyrażeniem
Ponieważ for , dlatego przybliżenie to jest dobre nawet dla dużych amplitud.
Szereg Fouriera dla ruchu od maksymalnego wychylenia
Dla wahadła rozpoczynającego ruch w chwili od maksymalnego wychylenia kątowego, zależność kąta odchylenia wahadła od pionu jest symetryczną funkcją czasu, tj. Wtedy współczynniki przed funkcją sinus w szeregu Fouriera są równe zero. Stąd otrzymuje się:
Ponadto, parzyste harmoniczne składowej cosinus są równe zero, w tym składowa . Wynika stąd, że w widmie Fouriera takiego ruchu wahadła występują składowe harmoniczne o częstotliwościach razy większych od częstotliwości podstawowej .
A także: Dla dużych amplitud drgań częstotliwość podstawowa jest mniejsza niż dla małych drgań, gdyż rośnie okres drgań ze wzrostem amplitudy, więc częstotliwość maleje.
Wnioski:
Z analizy tej wynika, że dokładne rozwiązanie równania ruchu wahadła zawiera nie jedną harmoniczną (jaką otrzymuje się w przybliżeniu funkcji sinus jednym wyrazem szeregu Maclaurina, tj. przyjmując - przy założeniu małych drgań), nie dwie harmoniczne (co otrzymuje się w rozwinięciu Maclaurina dwoma wyrazami, tj. przyjmując ), ale w ogólności składowych harmonicznych jest nieskończenie wiele. Amplitudy wyższych harmonicznych są określone przez współczynniki szeregu Fouriera. Wyższe harmoniczne mają coraz mniejsze amplitudy.
Całkowita energia mechaniczna wahadła. Punkty zwrotne. Krzywe fazowe
W opisie ruchu wahadła zamiast kąta maksymalnego wychylenia jako stały parametr ruchu można przyjąć całkowitą energię mechaniczną wahadła [7]. Energia ta determinuje punkty zwrotne ruchu wahadła, które charakteryzują kąty maksymalnego odchylenia.
Krzywe fazowe wahadła są to krzywe prezentujące zależność współrzędnej prędkości kątowej wahadła od kąta odchylenia wahadła od pionu. Niech oznacza energię potrzebną do odchylenia wahadła z położenia równowagi do pionu, tj. odchylenia o kąt Na podstawie wykresów fazowych można odróżnić poszczególne przypadki ruchu[13] (por. wykres obok): 1) gdy to krzywe fazowe są krzywymi zamkniętymi; 2) gdy to krzywe fazowe tworzą przecinające się linie; 3) gdy to krzywe fazowe są liniami otwartymi.
Poniżej zestawiono animacje pokazujące sposoby (mody) oscylacji wahadła matematycznego w zależności od jego energii całkowitej. Animacje pokazują, że okres drgań zależy od amplitudy. Małe wykresy powyżej wahadeł są wykresami fazowymi ruchu wahadeł.
- Kąt początkowy 0°, równowaga trwała.
- Kąt początkowy 45°
- Kąt początkowy 90°
- Kąt początkowy 135°
- Kąt początkowy 170°
- Kąt początkowy 180°, równowaga nietrwała.
- Wahadło o energii nieco większej niż potrzebna na pełny obrót.
- Wahadło o energii znacznie większej niż minimalna potrzebnej na pełny obrót.
Uniwersalność wzoru opisującego ruch wahadła. Co oznacza g?
Równanie ruchu wahadła jest uniwersalne - jest słuszne nie tylko dla drgań na Ziemi, ale dla drgań na dowolnym ciele niebieskim (jeśli takie wahadło można by tam zainstalować); np. na Księżycu, gdzie przyspieszenie grawitacyjne jest około 6 razy mniejsze niż na Ziemi, identyczne wahadło miałoby razy dłuższy okres drgań.
Jest też słuszne np. w windzie, spadającej z przyspieszeniem w kierunku Ziemi. Gdyby winda poruszała się z przyspieszeniem ziemskim, to przedmioty znajdujące się w niej byłyby w stanie nieważkości. Wyjaśnienie tego jest następujące: do opisu ruchu wahadła w układzie odniesienia związanym z windą, który z racji ruchu przyspieszonego względem Ziemi jest układem nieinercjalnym, trzeba wprowadzić oprócz siły grawitacji siłę bezwładności - jeżeli winda poruszałaby się z przyspieszeniem ziemskim, to siła ciężkości ciała byłaby równoważona przez siłę bezwładności, działającej na to ciało - ciało byłoby się w stanie nieważkości. W zależności od warunków początkowych ciało wahadła spoczywałoby (było w równowadze trwałej) albo poruszało się ruchem jednostajnym po okręgu[18].
W równaniu ruchu wahadła współczynnik należy więc rozumieć jako przyspieszenie ciał swobodnie spadających, mierzone w układzie nieruchomym względem punktu zamocowania nici wahadła. Łatwo pokazać zgodność powyższych efektów (obserwowanych np. w stacjach krążących wokół Ziemi) z opisem matematycznym: Z równania ogólnego ruchu wahadła dla otrzymamy proste równanie różniczkowe
z którego natychmiast wynika rozwiązanie:
- ,
gdzie: - położenie kątowe początkowe wahadła, - początkowa prędkość kątowa, nadana wahadłu.
Gdy , to wahadło będzie w pozycji nieruchomej; w przeciwnym razie będzie wykonywać ruch jednostajny po okręgu.
Opis ruchu wahadła w ramach mechaniki Lagrange'a
Z definicji wahadła prostego, jego ruch jest ograniczony przez więzy do ruchu po okręgu. Suma składowych sił działających na ciało prostopadłe do toru ruchu jest siłą dośrodkową, jej wartość określa wzór[19]
przy czym znak „minus” jest dlatego, że siła działa w stronę środka okręgu, przeciwnie do zwrotu współrzędnej układu współrzędnych biegunowych. Zależność tej siły od kąta można określić wyznaczając prędkość wahadła z zasady zachowania energii, co daje[19].
gdzie – kąt maksymalnego odchylenia wahadła.
Siłę napięcia nici określa wzór[19]
W przyjętym tu układzie współrzędnych biegunowych, który jest zgodny z więzami, wyznaczenie siły reakcji więzów jest niepotrzebne do opisu ruchu wahadła. Wyznaczenie tej siły byłoby konieczne, gdyby siły opisywać w układzie współrzędnych kartezjańskich. Dobór układu współrzędnych zgodnych z więzami stanowi podstawę sformułowania mechaniki klasycznej w ujęciu mechaniki Lagrange’a.
W ramach mechaniki Lagrange'a wyprowadza się równanie ruchu wahadła matematycznego (pomijamy to tutaj). Jednak metoda Lagrange'a, która zrodziła się w rozwiazywaniu problemów klasycznej fizyki, znajduje znacznie poważniejsze zastosowania we współczesnej teorii pola (w fizyce cząstek elementarnych), gdzie równania ruchu cząstek wyprowadza się poprzez konstrukcję odpowiedniego Lagranżjanu.
Wahadło matematyczne – tabela okresów drgań dla dowolnie dużych amplitud
2° | 1.0001 | 22° | 1.0093 | 42° | 1.0347 | 62° | 1.0785 | 82° | 1.1454 | 102° | 1.2439 | 122° | 1.3905 | 142° | 1.6238 | 162° | 2.0724 |
4° | 1.0003 | 24° | 1.0111 | 44° | 1.0382 | 64° | 1.0841 | 84° | 1.1537 | 104° | 1.2560 | 124° | 1.4090 | 144° | 1.6551 | 164° | 2.1453 |
6° | 1.0007 | 26° | 1.0130 | 46° | 1.0418 | 66° | 1.0898 | 86° | 1.1622 | 106° | 1.2686 | 126° | 1.4283 | 146° | 1.6884 | 166° | 2.2284 |
8° | 1.0012 | 28° | 1.0151 | 48° | 1.0457 | 68° | 1.0959 | 88° | 1.1711 | 108° | 1.2817 | 128° | 1.4485 | 148° | 1.7240 | 168° | 2.3248 |
10° | 1.0019 | 30° | 1.0174 | 50° | 1.0498 | 70° | 1.1021 | 90° | 1.1803 | 110° | 1.2953 | 130° | 1.4698 | 150° | 1.7622 | 170° | 2.4394 |
12° | 1.0027 | 32° | 1.0199 | 52° | 1.0540 | 72° | 1.1087 | 92° | 1.1899 | 112° | 1.3096 | 132° | 1.4922 | 152° | 1.8033 | 172° | 2.5801 |
14° | 1.0037 | 34° | 1.0225 | 54° | 1.0585 | 74° | 1.1155 | 94° | 1.1999 | 114° | 1.3244 | 134° | 1.5157 | 154° | 1.8478 | 174° | 2.7621 |
16° | 1.0049 | 36° | 1.0252 | 56° | 1.0632 | 76° | 1.1225 | 96° | 1.2103 | 116° | 1.3399 | 136° | 1.5405 | 156° | 1.8962 | 176° | 3.0193 |
18° | 1.0062 | 38° | 1.0282 | 58° | 1.0681 | 78° | 1.1299 | 98° | 1.2210 | 118° | 1.3560 | 138° | 1.5667 | 158° | 1.9492 | 178° | 3.4600 |
20° | 1.0077 | 40° | 1.0313 | 60° | 1.0732 | 80° | 1.1375 | 100° | 1.2322 | 120° | 1.3729 | 140° | 1.5944 | 160° | 2.0075 | 180° |
Okresy drgań wahadła matematycznego nietłumionego w zależność od amplitudy, zamieszczone w tabeli, obliczono korzystając ze wzoru[8]:
Dla mamy Dzieląc stronami oba wzory otrzyma się formułę, z której łatwo wyznaczyć szukane wielkości:
Wartości całek eliptycznych zupełnych pierwszego rodzaju są stabelaryzowane (patrz tu – należy odszukać kąt i odczytać wartość całki ).
Na podstawie tabeli można wyznaczyć okres drgań wahadła o zadanej długości mnożąc wartość odczytaną z tabeli przez
Tabela może służyć też do testowania algorytmów obliczania złożonych ruchów, co omówiono poniżej.
Drgania tłumione i wymuszone wahadła matematycznego
W omówionych tu zagadnieniach opisane były drgania swobodne wahadła matematycznego i fizycznego, tj. drgania odbywające się jedynie pod wpływem siły ciężkości oraz siły reakcji nici czy podpory. Siły grawitacji są siłami zachowawczymi. Podobnie zachowawcze są rozważane tu siły reakcji, gdy punkt zaczepienia nici jest nieruchomy (siła reakcji nie zależy wtedy jawnie od czasu i działa prostopadle do chwilowego kierunku ruchu wahadła; w konsekwencji taka siła reakcji nie wykonuje pracy). Energia mechaniczna wahadła poddanego działaniu tylko tych sił byłaby zachowana i w konsekwencji powodowałaby jego nieustanny ruch[20].
Rzeczywiste układy drgające wprawione w ruch po pewnym czasie zatrzymują się pod wpływem oporów ruchu (np. oporów powietrza), chyba że działa na nie siła wymuszająca ruch, jak to jest w przypadku wahadeł zegarów. Uwzględnienie sił oporów ruchu lub sił wymuszających ruch prowadzi do równań wahadła tłumionego lub wymuszonego (por. Ruch harmoniczny tłumiony oraz Oscylator harmoniczny wymuszony)[21]. Rozwiązanie tych równań w ogólnym przypadku jest niezwykle złożone. Z pomocą przychodzą metody numeryczne, niżej omówione.
Całkowanie numeryczne równań ruchu wahadła matematycznego
Wahadło nietłumione
Równanie ogólne drgań wahadła matematycznego, podane na początku artykułu, jest równaniem różniczkowym zwyczajnym drugiego rzędu. Jest to ponadto równanie nieliniowe – nie da się go rozwiązać analitycznie. Można jednak rozwiązać je efektywnie metodami numerycznymi, np. stosując metody Rungego-Kutty. Metoda polega na tym, że wprowadzając nową zmienną (która de facto ma sens prędkości kątowej wahadła) równanie to sprowadza się do układu dwóch równań różniczkowych pierwszego rzędu
a następnie układ tych równań rozwiązuje się iteracyjnie, znajdując np. zależność kata odchylenia wahadła w zależności od dyskretnych chwil czasu czy też okres drgań (Przykład kodu programu w C++ podano tutaj)
Wahadło tłumione, z siłą wymuszającą
Znalezienie rozwiązania analitycznego równania ruchu złożonych układów fizycznych jest w ogólnym przypadku niemożliwe. Jednak metody numeryczne pozwalają efektywnie rozwiązywać te równania ruchu. Np. równanie ruchu wahadła z tłumieniem i z siłą wymuszającą ma postać
gdzie:
- – współczynnik tłumienia,
- – siła wymuszająca, zależna dowolnie od czasu.
Wprowadzając nową zmienną powyższe równanie sprowadza się do układu dwóch równań różniczkowych pierwszego rzędu
Układ ten rozwiązuje się iteracyjnie, np. metodą Rungego-Kutty, co prowadzi do znalezienia dyskretnych wartości w zależności od dyskretnych chwil czasu w zadanym przedziale całkowania równań. (Przykład kodu programu w języku python, wraz z generowaniem wykresów drgań, podano tutaj).
Wahadło fizyczne
Wahadło fizyczne jest to bryła sztywna zawieszona na stałej osi poziomej w jednorodnym polu grawitacyjnym, która wykonuje drgania polegające na obrotach wokół tej osi raz w jedną, raz w drugą stronę. Na wychylone z położenia równowagi wahadło działa moment siły[1]:
Stąd, korzystając z II zasady dynamiki ruchu obrotowego, otrzymuje się równanie ruchu wahadła fizycznego
gdzie:
- – odległość od punktu zawieszenia do środka ciężkości bryły,
- – przyspieszenie ziemskie,
- – moment bezwładności wahadła względem osi obrotu,
- – masa wahadła.
Podstawiając do powyższego równania tzw. długość zredukowaną
równanie ruchu wahadła fizycznego przyjmie identyczną postać jak równanie ruchu wahadła matematycznego. Oznacza to, że wszystkie wnioski dotyczące ruchu wahadła fizycznego są identyczne z wnioskami dotyczącymi wahadła matematycznego. Przykładowo okres drgań wahadła fizycznego dany jest wzorem[1]:
czyli wzorem, jakie miałoby wahadło matematyczne o długości równej długości zredukowanej wahadła fizycznego.
Wahadło matematyczne jako szczególne wahadło fizyczne
Z drugiej strony wahadło matematyczne, które jest masą punktową zawieszoną na nieważkiej nici, może być traktowane jako szczególna bryła sztywna, ma bowiem ściele określony moment bezwładności i odległość od środka obrotu do środka masy, danych oczywistymi wzorami:[1]
Zastosowania wahadła fizycznego
Wahadło fizyczne stosuje się jako przyrząd do dokładnego pomiaru przyspieszenia ziemskiego[1]. Przykładem wahadła do pomiaru przyspieszenia ziemskiego (oraz jako przyrządu dydaktycznego) jest wahadło rewersyjne.
Inne rodzaje wahadeł
W fizyce rozważa się kilka modeli wahadeł, które nie spełniają założeń wahadła matematycznego lub fizycznego. Przykładami są:
- wahadło sferyczne – ciało na nierozciągliwej nici, ale jego ruch nie jest ograniczony do płaszczyzny,
- wahadło stożkowe – ciało na nierozciągliwej nici, a ciało porusza się po okręgu,
- wahadło podwójne – ciało wahadła jest punktem zawieszenia kolejnego wahadła, może być rozważane jako płaskie i sferyczne, matematyczne i fizyczne,
- wahadło z rozciągliwą nicią,
- tautochrona (wahadło cykloidalne) – wahadło o okresie niezależnym od amplitudy drgań[22][23] (niżej omówione),
- wahadło Foucaulta (niżej omówione)
Inne układy drgające:
Wahadło o okresie niezależnym od amplitudy
Problemem konstruowania dokładnych zegarów wahadłowych, których szybkość chodu nie zależy od amplitudy drgań, zajmował się Christiaan Huygens, wykazał, że niezależność szybkości chodu zapewni zmniejszanie długość nici wahadła wraz z wychyleniem; następnie wykazał, że zrealizuje to wahadło cykloidalne, tj. wahadło, w którym nić lub elastyczny element zawieszenia, będzie owijać się na cykloidzie o poziomej osi i promieniu równym ćwierci długości wahadła. W ten sposób skonstruował wahadło o okresie niezależnym od amplitudy[24].
Problem wahadła o okresie niezależnym od amplitudy sprowadza się do wyznaczenia takiej krzywej, że ciało poruszając pod działaniem stałej siły grawitacji po niej w takim samym czasie przemieści się od punktu ruszenia do jej najniższego punktu. Krzywa zwana jest tautochroną i jest cykloidą[25].
Wahadło Foucaulta
Płaszczyzna drgań wahadła znajdującego się na Ziemi poza równikiem powoli obraca się względem Ziemi. Zjawisko to można wyjaśnić jako efekt działania siły Coriolisa wywołanej ruchem wahadła na obracającej się Ziemi. Wahadło umożliwiające obserwację tego efektu, jest nazywane wahadłem Foucaulta[26].
Okres obrotu płaszczyzny wahadła w dla obserwatora znajdującego się na obracającej się Ziemi opisuje wzór[26]:
gdzie – szerokość geograficzna, na której znajduje się wahadło. Np. dla szerokości geograficznej 52° (okolice Warszawy) okres wahadła Foucaulta wynosi około 30 h 27 min i maleje ze wzrostem szerokości geograficznej. Na biegunach okres ten wynosi 24 h.
Zobacz też
- Przyrządy będące wahadłami
- wahadło balistyczne
- wahadło Newtona
- wahadło Oberbecka
- wahadło radiestezyjne
- wahadło rewersyjne
- wahadło zegarowe
- Oscylatory
- Inne
Uwagi
Przypisy
- ↑ a b c d e f g h Resnick i Halliday 1999 ↓, s. 361–364.
- ↑ Huygens’ Clocks. [dostęp 2015-05-01].
- ↑ a b Victor F. Lenzen, Robert P. Multhauf. Development of Gravity Pendulums in the 19th Century. „On Science and Technology”. papers 34-44. On Science and Technology, Smithsonian Institution.
- ↑ a b c d Królikowski i Rubinowicz 2012 ↓, s. 91.
- ↑ Królikowski i Rubinowicz 2012 ↓, s. 44.
- ↑ Oscillations and Fourier Analysis. [dostęp 2016-09-10].
- ↑ a b Landau i Lifszyc 2011 ↓, s. 41–46.
- ↑ a b c Królikowski i Rubinowicz 2012 ↓, s. 97.
- ↑ a b c G. Białkowski, s. 244
- ↑ Gaetan Kerschen, Douglas Adams, Alex Carrella: Topics in Nonlinear Dynamics, Volume 1. T. 1: Proceedings of the 31st IMAC, A Conference and Exposition on Structural Dynamics. ISBN 978-1-4614-6570-6.
- ↑ Kalkulator okresów drgań wahadła matematycznego dla dowolnej amplitudy. (ang.).
- ↑ a b Wróblewski i Zakrzewski 1976 ↓, s. 340–343.
- ↑ a b c d Kittel, Knight i Ruderman 1993 ↓, s. 256–257.
- ↑ Eric W.E.W. Weisstein Eric W.E.W., Elliptic Integral of the First Kind, [w:] MathWorld, Wolfram Research [dostęp 2015-04-25] (ang.).
- ↑ a b Królikowski i Rubinowicz 2012 ↓, s. 96.
- ↑ Derek F.D.F. Lawden Derek F.D.F., Elliptic Functions and Applications, Springer-Verlag, 1989, ISBN 0-387-96965-9 . Eq. 2.7.9:
- ↑ W.P.W.P. Reinhardt W.P.W.P., P.L.P.L. Walker P.L.P.L., Jacobian Elliptic Functions .
- ↑ Wróblewski i Zakrzewski 1984 ↓, s. 371.
- ↑ a b c Wróblewski i Zakrzewski 1984 ↓, s. 340–342.
- ↑ Królikowski i Rubinowicz 2012 ↓, s. 64, 72.
- ↑ Królikowski i Rubinowicz 2012 ↓, s. 47–58.
- ↑ Alan Emmerson: Things Are Seldom What They Seem – Christiaan Huygens, the Pendulum and the Cycloid. [dostęp 2015-04-25].
- ↑ Eric W.E.W. Weisstein Eric W.E.W., Tautochrone Problem, [w:] MathWorld, Wolfram Research [dostęp 2015-04-25] (ang.).
- ↑ Alan Emmerson: Things Are Seldom What They Seem – Christiaan Huygens, the Pendulum and the Cycloid. [dostęp 2015-04-25].
- ↑ Marek Kordos: Pierwszy nowoczesny zegarmistrz. [dostęp 2015-05-10]. [zarchiwizowane z tego adresu (2007-07-09)].
- ↑ a b Wróblewski i Zakrzewski 1976 ↓, s. 171–172.
Bibliografia
- G. Białkowski, Mechanika klasyczna, Warszawa: PWN, 1975, s.238-245.
- M. Jeżewski, Fizyka, PWN, Warszawa 1966.
- C. Kittel, W.D. Knight, M.A. Ruderman: Mechanika. Warszawa: PWN, 1993.
- W. Królikowski, W. Rubinowicz: Mechanika teoretyczna. Warszawa: PWN, 2012.
- Lew Dawidowicz Landau, E.M. Lifszyc: Mechanika. Warszawa: PWN, 2011.
- Robert Resnick, David Halliday: Fizyka. Wydawnictwo Naukowe PWN, 1999. ISBN 83-01-09323-4.
- A.K. Wróblewski, J.A. Zakrzewski: Wstęp do fizyki. T. 1. Warszawa: PWN, 1976.
Linki zewnętrzne
Zobacz hasło wahadło w Wikisłowniku |
- Ćwiczenie Pomiar natężenia pola grawitacyjnego w Siedlcach przy pomocy modelu wahadła matematycznego. kf.imif.uph.edu.pl. [zarchiwizowane z tego adresu (2015-04-13)]. – w ramach zestawu ćwiczeń do fizyki na Wydziale Nauk Ścisłych Uniwersytetu Przyrodniczo-Humanistycznego w Siedlcach
- Wahadło fizyczne (ang.)
- p
- d
- e
działy |
| ||||||||
---|---|---|---|---|---|---|---|---|---|
formalizmy | |||||||||
wielkości |
| ||||||||
rodzaje ruchu | |||||||||
zjawiska | |||||||||
inne pojęcia | |||||||||
prawa |
| ||||||||
uogólnienia | |||||||||
uczeni według daty narodzin |
|
- SNL: matematisk_pendel