Wszechnica Wszechwiedzy - Baner

Tarcie w statyce – tarcie statyczne, graniczne i tarcie cięgna o krążek

Tarcie jest oddziaływaniem występującym pomiędzy stykającymi się powierzchniami. Przeciwdziała ich względnemu przesuwaniu albo tendencji do rozpoczęcia takiego ruchu. W statyce szczególne znaczenie ma tarcie statyczne, dzięki któremu ciało może pozostawać w równowadze mimo działania sił próbujących wprawić je w ruch.

Siła tarcia nie zawsze przyjmuje swoją największą możliwą wartość. Dopiero na granicy poślizgu występuje tak zwane tarcie graniczne lub całkowicie rozwinięte. Osobnym zagadnieniem jest tarcie cięgna o nieruchomy krążek, bęben albo słup, opisywane równaniem Eulera – Eytelweina.

Przed rozpoczęciem analizy warto przypomnieć sobie zasady sporządzania diagramu ciała swobodnego, omówione w artykule Czym zajmuje się statyka? Podstawowe pojęcia i zasady. Dodawanie sił oraz rozkładanie ich na składowe przedstawiliśmy natomiast w materiale Rachunek wektorowy w statyce – siły, składowe i momenty.

Czym jest tarcie?

Tarcie jest siłą styczną do powierzchni kontaktu. Powstaje wskutek oddziaływań pomiędzy nierównościami stykających się powierzchni, ich odkształceń oraz zjawisk zachodzących na poziomie mikroskopowym.

W prostych zadaniach ze statyki najczęściej stosuje się model tarcia Coulomba. Zakłada on, że maksymalna wartość tarcia zależy od:

Współczynnik tarcia jest wielkością bezwymiarową. Zależy między innymi od rodzaju materiałów, stanu powierzchni, obecności smaru, wilgotności oraz innych warunków kontaktu.

Jak wyznaczyć kierunek siły tarcia?

Siła tarcia jest skierowana przeciwnie do względnego ruchu stykających się powierzchni albo przeciwnie do tendencji do rozpoczęcia takiego ruchu.

Nie należy automatycznie przyjmować, że tarcie jest zawsze skierowane przeciwnie do pojedynczej siły zewnętrznej. Najpierw trzeba ustalić, jak poruszyłoby się ciało, gdyby powierzchnia była idealnie gładka.

Jeżeli klocek miałby tendencję do przesuwania się w prawo, tarcie działające na klocek będzie skierowane w lewo. Jeżeli natomiast inne siły powodowałyby tendencję ruchu w lewo, zwrot tarcia należałoby odwrócić.

Klocek na poziomej chropowatej powierzchni z siłą P, ciężarem G, reakcją normalną N i siłą tarcia T
Siła tarcia statycznego jest skierowana przeciwnie do przewidywanego kierunku względnego ruchu powierzchni.

Na przedstawionym schemacie działają cztery siły:

Jeżeli przyjęty zwrot tarcia okaże się niewłaściwy, z równań równowagi otrzymamy ujemną wartość \(T\). Nie oznacza to zwykle błędu rachunkowego, lecz informuje, że rzeczywisty zwrot siły jest przeciwny do założonego.

Tarcie statyczne

Tarcie statyczne występuje wtedy, gdy stykające się powierzchnie nie przesuwają się względem siebie. Jego wartość dostosowuje się do warunków równowagi.

\[ |T|\leq \mu_s N, \] 

gdzie:

Nierówność jest niezwykle istotna. W ogólnym stanie równowagi tarcie nie musi być równe \(\mu_sN\). Może przyjmować każdą wartość od zera do wartości granicznej, zależnie od tego, jaka siła jest potrzebna do utrzymania równowagi.

Przykład: jeżeli ciało jest ciągnięte siłą \(P=3\ \mathrm{N}\), a maksymalna możliwa siła tarcia wynosi \(10\ \mathrm{N}\), to w stanie równowagi tarcie ma wartość \(3\ \mathrm{N}\), a nie \(10\ \mathrm{N}\).

Tarcie statyczne może być również równe zeru. Dzieje się tak wtedy, gdy żadna z pozostałych sił nie wywołuje tendencji do względnego przesunięcia powierzchni.

Tarcie graniczne i całkowicie rozwinięte

Gdy ciało znajduje się na granicy poślizgu, tarcie statyczne osiąga największą możliwą wartość. Mówimy wtedy o tarciu granicznym albo tarciu całkowicie rozwiniętym.

\[ |T|=T_{\max}=\mu_sN. \]

Równość tę stosujemy wtedy, gdy z treści zadania wynika, że:

Ważne: w zadaniach ze statyki nie należy od razu podstawiać \(T=\mu_sN\). Najpierw trzeba ustalić, czy rozpatrywany stan rzeczywiście odpowiada granicy poślizgu.

Wykres siły tarcia pokazujący tarcie statyczne, maksymalne tarcie graniczne i mniejsze tarcie kinetyczne
Tarcie statyczne dostosowuje się do wymuszenia aż do wartości granicznej, po której rozpoczęciu poślizgu występuje tarcie kinetyczne.

W zakresie tarcia statycznego wzrost siły próbującej przesunąć ciało powoduje odpowiedni wzrost siły tarcia. Po osiągnięciu wartości \(T_{\max}\) dalsze zwiększanie wymuszenia prowadzi do rozpoczęcia poślizgu.

Tarcie kinetyczne

Tarcie kinetyczne, nazywane również tarciem ślizgowym, występuje podczas względnego przesuwania się powierzchni. W prostym modelu Coulomba jego wartość określa zależność:

\[ |T_k|=\mu_kN, \]

gdzie \(\mu_k\) jest współczynnikiem tarcia kinetycznego.

Zwykle zachodzi:

\[ \mu_k<\mu_s. \]

Oznacza to, że do rozpoczęcia ruchu potrzebna jest zazwyczaj większa siła niż do jego późniejszego podtrzymywania. Ponieważ w statyce analizujemy ciała pozostające w spoczynku, podstawowe znaczenie ma współczynnik \(\mu_s\). Tarcie kinetyczne należy już zasadniczo do zagadnień dynamiki.

Reakcja całkowita i kąt tarcia

Oddziaływanie chropowatego podłoża na ciało można przedstawić za pomocą dwóch składowych:

Wypadkową obu składowych jest reakcja całkowita \(\vec{R}\):

\[ \vec{R}=\vec{N}+\vec{T}. \]

Jej wartość wynosi:

\[ R=\sqrt{N^2+T^2}. \]

Kąt pomiędzy reakcją całkowitą a kierunkiem normalnym nazywamy kątem tarcia \(\varphi\). W stanie granicznym:

\[ \tan\varphi = \frac{T_{\max}}{N} = \mu_s. \]

Stąd:

\[ \varphi=\arctan\mu_s. \]
Reakcja normalna N i siła tarcia T złożone w reakcję całkowitą R odchyloną o kąt tarcia
W stanie granicznym reakcja całkowita podłoża jest odchylona od normalnej o kąt tarcia spełniający zależność \(\tan\varphi=\mu_s\).

W zagadnieniach przestrzennych możliwe kierunki reakcji całkowitej w stanie granicznym tworzą tak zwany stożek tarcia. Jeżeli wypadkowa oddziaływań znajduje się wewnątrz stożka, równowaga jest możliwa. Na jego powierzchni występuje stan graniczny, a poza nim dochodzi do poślizgu.

Tarcie na równi pochyłej

Klasycznym przykładem jest ciało spoczywające na chropowatej równi nachylonej pod kątem \(\alpha\). Ciężar \(G\) rozkładamy na dwie składowe:

\[ G_{\parallel}=G\sin\alpha, \] \[ G_{\perp}=G\cos\alpha. \]

Jeżeli na ciało nie działają inne siły posiadające składową normalną do równi, reakcja podłoża wynosi:

\[ N=G\cos\alpha. \]

Składowa \(G\sin\alpha\) próbuje zsunąć ciało w dół równi. Siła tarcia jest więc skierowana w górę równi.

Klocek na chropowatej równi pochyłej z ciężarem G, reakcją N, tarciem T oraz składowymi ciężaru
Tarcie działa przeciwnie do tendencji zsuwania się klocka, a na granicy poślizgu kąt nachylenia spełnia zależność \(\tan\alpha=\mu_s\).

W stanie równowagi:

\[ T=G\sin\alpha. \]

Równowaga jest możliwa, jeżeli:

\[ G\sin\alpha \leq \mu_sG\cos\alpha. \]

Po skróceniu przez \(G\cos\alpha\) otrzymujemy:

\[ \tan\alpha\leq\mu_s. \]

Na granicy zsuwania:

\[ \tan\alpha=\mu_s. \]

Jeżeli na ciało działa dodatkowa siła próbująca przesunąć je w górę równi, kierunek tarcia może się odwrócić. Tarcie nie jest więc z definicji zawsze skierowane w górę równi — zależy od tendencji ruchu całego układu.

Jak rozwiązywać zadania z tarciem?

W typowym zadaniu można zastosować następującą procedurę:

  1. Wyodrębnij analizowane ciało i sporządź diagram ciała swobodnego.
  2. Narysuj ciężar, reakcję normalną, siły zewnętrzne oraz możliwą siłę tarcia.
  3. Ustal, w którą stronę ciało miałoby tendencję się poruszyć, gdyby tarcia nie było.
  4. Skieruj siłę tarcia przeciwnie do przewidywanego ruchu względnego.
  5. Zapisz równania równowagi sił i momentów.
  6. Oblicz wartość tarcia wymaganą do zachowania równowagi.
  7. Sprawdź, czy spełniony jest warunek \(|T|\leq\mu_sN\).

Jeżeli obliczona wartość spełnia nierówność, tarcie statyczne może utrzymać ciało w równowadze. Jeżeli natomiast:

\[ |T|>\mu_sN, \]

utrzymanie spoczynku nie jest możliwe i rozpocznie się poślizg.

Jeżeli zadanie dotyczy bezpośrednio granicy ruchu, można od początku przyjąć:

\[ |T|=\mu_sN. \]

Reakcja normalna nie zawsze jest równa ciężarowi. Zależność \(N=G\) obowiązuje tylko w prostym przypadku poziomego podłoża, gdy pozostałe siły nie mają składowych pionowych. Siła przyłożona pod kątem może zwiększać albo zmniejszać nacisk, a tym samym zmieniać maksymalną wartość tarcia.

Tarcie cięgna o krążek lub bęben

Jeżeli lina, pas albo inne wiotkie cięgno opasuje chropowaty krążek, bęben, słup lub kabestan, siły naciągu po obu stronach nie muszą być jednakowe.

Tarcie rozłożone na całej długości kontaktu umożliwia utrzymanie znacznie większej siły po stronie napiętej przy stosunkowo niewielkiej sile po stronie luźnej. Zjawisko to wykorzystuje się między innymi w kabestanach, hamulcach taśmowych, urządzeniach linowych i układach transportowych.

Wolne odcinki cięgna odchodzą od bębna stycznie do jego powierzchni. Siły naciągu działają wzdłuż osi tych odcinków, a więc również stycznie do krążka w punktach zejścia liny.

Lina stycznie opasująca nieruchomy bęben z większym i mniejszym naprężeniem oraz zaznaczonym kątem opasania
Dla cięgna znajdującego się na granicy poślizgu obowiązuje równanie Eulera – Eytelweina: \(S_{\mathrm{większa}}=S_{\mathrm{mniejsza}}e^{\mu\beta}\), gdzie \(\mu\) jest współczynnikiem tarcia, a \(\beta\) — kątem opasania wyrażonym w radianach.

Równanie Eulera – Eytelweina

W stanie granicznym, gdy cięgno znajduje się na granicy poślizgu względem nieruchomego bębna, obowiązuje zależność:

\[ \frac{S_{\mathrm{większa}}} {S_{\mathrm{mniejsza}}} = e^{\mu_s\beta}. \]

Można ją również zapisać jako:

\[ S_{\mathrm{większa}} = S_{\mathrm{mniejsza}}e^{\mu_s\beta}, \]

gdzie:

Równanie wynika z analizy bardzo małego fragmentu cięgna obejmującego kąt \(d\theta\). Dla stanu granicznego przyrost siły naciągu spełnia zależność:

\[ dS=\mu_sS\,d\theta. \]

Po rozdzieleniu zmiennych:

\[ \frac{dS}{S}=\mu_s\,d\theta. \]

Całkując od strony luźnej do strony napiętej oraz od kąta \(0\) do \(\beta\), otrzymujemy:

\[ \int_{S_{\mathrm{mniejsza}}}^{S_{\mathrm{większa}}} \frac{dS}{S} = \mu_s\int_0^\beta d\theta, \] \[ \ln \frac{S_{\mathrm{większa}}} {S_{\mathrm{mniejsza}}} = \mu_s\beta, \]

a po podniesieniu liczby \(e\) do obu stron równania:

\[ \frac{S_{\mathrm{większa}}} {S_{\mathrm{mniejsza}}} = e^{\mu_s\beta}. \]

Jeżeli cięgno pozostaje w równowadze, lecz tarcie nie jest jeszcze całkowicie rozwinięte, zachodzi jedynie nierówność:

\[ \frac{S_{\mathrm{większa}}} {S_{\mathrm{mniejsza}}} \leq e^{\mu_s\beta}. \]

Równość stosujemy zatem na granicy poślizgu. Podczas rzeczywistego ślizgania można zastosować analogiczną zależność ze współczynnikiem tarcia kinetycznego \(\mu_k\), o ile przyjęty model tarcia jest odpowiedni.

Kąt opasania musi być podany w radianach

W wykładniku funkcji wykładniczej kąt \(\beta\) musi być wyrażony w radianach:

\[ \beta_{\mathrm{rad}} = \beta_{\mathrm{deg}} \frac{\pi}{180^\circ}. \]

Przykładowo:

Każde dodatkowe owinięcie zwiększa kąt \(\beta\) o \(2\pi\), dlatego nawet niewielka siła po stronie luźnej może równoważyć bardzo dużą siłę po stronie napiętej.

Strona napięta i strona luźna cięgna

Nie można z góry przyjąć, że większa siła znajduje się zawsze po lewej albo po prawej stronie rysunku. Zależy to od układu pozostałych sił oraz przewidywanego kierunku poślizgu cięgna.

Aby wskazać stronę napiętą, należy:

  1. przeanalizować obciążenia działające na oba końce cięgna;
  2. ustalić, w którą stronę lina próbowałaby przesunąć się względem bębna;
  3. wskazać stronę wywołującą ten ruch jako stronę napiętą;
  4. oznaczyć po tej stronie \(S_{\mathrm{większa}}\);
  5. po przeciwnej stronie oznaczyć \(S_{\mathrm{mniejsza}}\).

Tarcie działające na cięgno jest skierowane przeciwnie do przewidywanego poślizgu. Wzdłuż łuku kontaktu siła naciągu zmniejsza się od strony napiętej w kierunku strony luźnej.

Dwa przypadki cięgna opasującego bęben z większą siłą raz po lewej, a raz po prawej stronie
Strona, po której występuje większa siła naciągu cięgna, zależy od przewidywanego kierunku jego poślizgu względem bębna.

Najbezpieczniejszy zapis: najpierw posługujemy się oznaczeniami \(S_{\mathrm{większa}}\) i \(S_{\mathrm{mniejsza}}\), a dopiero po ustaleniu kierunku możliwego ruchu przypisujemy im symbole \(S_1\), \(S_2\), stronę lewą lub prawą.

Krążek gładki a chropowaty bęben

Równania Eulera – Eytelweina nie należy stosować automatycznie do każdego układu zawierającego linę i krążek.

W idealnym modelu gładkiego krążka, nieważkiej liny i pomijalnego oporu osi przyjmuje się:

\[ S_1=S_2. \]

Równanie Eulera – Eytelweina stosuje się natomiast wtedy, gdy:

Typowymi przykładami są lina opasująca nieruchomy słup, bęben hamulcowy, kabestan albo nieruchomy chropowaty krążek.

Jednostki i kontrola poprawności

WielkośćOznaczenieJednostka
Siła tarcia\(T\)\(\mathrm{N}\), \(\mathrm{kN}\)
Reakcja normalna\(N\)\(\mathrm{N}\), \(\mathrm{kN}\)
Reakcja całkowita\(R\)\(\mathrm{N}\), \(\mathrm{kN}\)
Siła naciągu cięgna\(S\)\(\mathrm{N}\), \(\mathrm{kN}\)
Współczynnik tarcia\(\mu\)wielkość bezwymiarowa
Kąt tarcia lub opasania\(\varphi\), \(\beta\)radian lub stopień, przy czym w wykładniku używamy radianów

W zależności:

\[ T_{\max}=\mu_sN \]

współczynnik \(\mu_s\) jest bezwymiarowy, dlatego prawa strona ma jednostkę siły.

W równaniu Eulera – Eytelweina wykładnik również musi być bezwymiarowy:

\[ [\mu\beta]=1. \]

Obie wielkości \(S_{\mathrm{większa}}\) i \(S_{\mathrm{mniejsza}}\) są siłami i muszą być wyrażone w tych samych jednostkach.

Uwaga terminologiczna: symbol \(S\) oznacza tutaj siłę naciągu cięgna wyrażoną w niutonach lub kiloniutonach. Nie jest to naprężenie normalne \(\sigma\), którego jednostką byłby paskal.

Najczęstsze błędy

Podsumowanie

Tarcie statyczne dostosowuje swoją wartość do warunków równowagi i spełnia nierówność:

\[ |T|\leq\mu_sN. \]

Dopiero na granicy poślizgu tarcie jest całkowicie rozwinięte:

\[ |T|=\mu_sN. \]

Po rozpoczęciu względnego ruchu powierzchni występuje tarcie kinetyczne:

\[ |T_k|=\mu_kN. \]

W stanie granicznym reakcja całkowita jest odchylona od kierunku normalnego o kąt tarcia:

\[ \tan\varphi=\mu_s. \]

Dla cięgna opasującego nieruchomy chropowaty bęben obowiązuje na granicy poślizgu równanie Eulera – Eytelweina:

\[ S_{\mathrm{większa}} = S_{\mathrm{mniejsza}}e^{\mu_s\beta}. \]

Najważniejszym etapem rozwiązania zadania jest prawidłowe ustalenie przewidywanego kierunku ruchu. Dopiero na tej podstawie można określić zwrot siły tarcia, wskazać napiętą i luźną stronę cięgna oraz zdecydować, czy należy zastosować nierówność opisującą tarcie statyczne, czy równość odpowiadającą stanowi granicznemu.

Metody naiwne prognozowania – poziom, przyrost i tempo zmian

Metody naiwne prognozowania należą do najprostszych sposobów przewidywania przyszłych wartości szeregu czasowego. Nie wymagają budowania rozbudowanych modeli ani szacowania wielu parametrów. Prognoza powstaje przez przyjęcie, że w przyszłości utrzyma się ostatni znany poziom zjawiska, ostatni przyrost albo ostatnie tempo zmian.

Określenie „naiwne” nie oznacza jednak, że metody te są zawsze nieracjonalne lub bezużyteczne. Dla zjawisk względnie stabilnych prosta prognoza może być dobrze uzasadniona. Metody naiwne są również ważnym punktem odniesienia podczas oceny bardziej zaawansowanych modeli.

W artykule omawiamy trzy warianty często spotykane w polskim nauczaniu akademickim: poziom bez zmian, przyrost bez zmian oraz przyrost procentowy bez zmian. Przedstawiamy także prognozowanie na podstawie średniego tempa zmian oraz pokazujemy, jak poprawnie oceniać prognozy jednookresowe i wielookresowe.


Spis treści

Czym są metody naiwne prognozowania?

Metody naiwne są prostymi metodami prognozowania, w których przyszłe wartości wyznacza się bez budowania złożonego modelu statystycznego. Wykorzystuje się przede wszystkim ostatnie znane obserwacje oraz najprostsze prawidłowości zauważone w szeregu czasowym.

Podstawowe informacje o szeregach czasowych i ich składnikach zostały przedstawione w artykule wprowadzającym do prognozowania.

W metodzie naiwnej można przyjąć między innymi, że:

Metoda naiwna nie próbuje dokładnie wyjaśnić mechanizmu badanego zjawiska. Zakłada jedynie, że najprostsza prawidłowość zaobserwowana w ostatnich okresach będzie kontynuowana w przyszłości.

Metody te mogą być stosowane samodzielnie, ale bardzo często pełnią funkcję modelu odniesienia. Bardziej skomplikowany model powinien tworzyć prognozy wyraźnie dokładniejsze niż odpowiednio dobrana metoda naiwna. Jeżeli tak nie jest, stosowanie bardziej złożonego rozwiązania może nie przynosić rzeczywistej korzyści.

Oznaczenia stosowane we wzorach

Przyjmijmy, że:

Prognoza na jeden okres naprzód jest więc oznaczana jako:

\[ \widehat{y}_{T+1\mid T} \]

Prognoza na \(h\) okresów naprzód ma natomiast postać:

\[ \widehat{y}_{T+h\mid T} \]

Metoda naiwna – poziom bez zmian

Najprostsza metoda naiwna zakłada, że przyszła wartość zmiennej będzie równa ostatniej znanej obserwacji. Jest to wariant określany jako poziom bez zmian.

Prognozę zapisujemy jako:

\[ \widehat{y}_{T+h\mid T}=y_T \]

Oznacza to, że niezależnie od długości horyzontu prognozy wszystkie przyszłe wartości są równe ostatniej obserwacji.

Przykład

Załóżmy, że sprzedaż w ostatnim znanym miesiącu wyniosła 120 sztuk:

\[ y_T=120 \]

Prognozy na trzy kolejne miesiące wynoszą:

\[ \widehat{y}_{T+1\mid T}=120 \]

\[ \widehat{y}_{T+2\mid T}=120 \]

\[ \widehat{y}_{T+3\mid T}=120 \]

Na wykresie ścieżka prognozy jest poziomą linią wychodzącą z ostatniej znanej wartości.

Kiedy metoda może działać dobrze?

Metoda poziomu bez zmian może być użyteczna, gdy:

Metoda może natomiast tworzyć systematycznie zaniżone prognozy w szeregu rosnącym oraz zawyżone prognozy w szeregu malejącym.

Metoda naiwna – przyrost bez zmian

W metodzie przyrost bez zmian zakładamy, że w przyszłości będzie powtarzał się ostatni zaobserwowany przyrost absolutny.

Ostatni przyrost obliczamy jako:

\[ \Delta_T=y_T-y_{T-1} \]

Prognoza na \(h\) okresów naprzód ma postać:

\[ \widehat{y}_{T+h\mid T} = y_T+h\Delta_T \]

Po podstawieniu definicji przyrostu otrzymujemy:

\[ \widehat{y}_{T+h\mid T} = y_T+h\left(y_T-y_{T-1}\right) \]

Metoda przedłuża w przyszłość ostatnią zaobserwowaną różnicę pomiędzy kolejnymi okresami.

Przykład wzrostu

Załóżmy, że sprzedaż w dwóch ostatnich miesiącach wynosiła:

\[ y_{T-1}=100, \qquad y_T=108 \]

Ostatni przyrost wynosi:

\[ \Delta_T=108-100=8 \]

Prognoza na następny okres:

\[ \widehat{y}_{T+1\mid T}=108+8=116 \]

Prognoza na dwa okresy naprzód:

\[ \widehat{y}_{T+2\mid T}=108+2\cdot8=124 \]

Prognoza na trzy okresy naprzód:

\[ \widehat{y}_{T+3\mid T}=108+3\cdot8=132 \]

Przykład spadku

Jeżeli wartość zmniejszyła się ze 100 do 94, ostatni przyrost wynosi:

\[ \Delta_T=94-100=-6 \]

Prognozowane wartości będą więc zmniejszały się o 6 jednostek w każdym kolejnym okresie:

\[ 88,\quad82,\quad76,\ldots \]

Wadą metody jest duża wrażliwość na dwie ostatnie obserwacje. Jeżeli ostatnia zmiana była nietypowa, jej mechaniczne powtarzanie może prowadzić do szybko narastających błędów.

Metoda naiwna – przyrost procentowy bez zmian

Metoda przyrost procentowy bez zmian zakłada, że w przyszłości będzie utrzymywało się ostatnie zaobserwowane tempo zmian.

Ostatnią stopę wzrostu obliczamy jako:

\[ g_T= \frac{y_T-y_{T-1}}{y_{T-1}} \]

Możemy również wyznaczyć ostatni indeks łańcuchowy:

\[ i_T= \frac{y_T}{y_{T-1}} = 1+g_T \]

Prognoza na \(h\) okresów naprzód wynosi:

\[ \widehat{y}_{T+h\mid T} = y_T(1+g_T)^h \]

Równoważnie można ją zapisać przy użyciu indeksu łańcuchowego:

\[ \widehat{y}_{T+h\mid T} = y_Ti_T^h \]

Lub bezpośrednio za pomocą dwóch ostatnich obserwacji:

\[ \widehat{y}_{T+h\mid T} = y_T \left( \frac{y_T}{y_{T-1}} \right)^h \]

Przykład

Załóżmy, że wartość zmiennej wzrosła ze 100 do 110:

\[ y_{T-1}=100, \qquad y_T=110 \]

Ostatnie tempo wzrostu wynosi:

\[ g_T= \frac{110-100}{100} = 0{,}10 = 10\% \]

Prognoza na następny okres:

\[ \widehat{y}_{T+1\mid T} = 110\cdot1{,}10 = 121 \]

Prognoza na dwa okresy naprzód:

\[ \widehat{y}_{T+2\mid T} = 110\cdot1{,}10^2 = 133{,}1 \]

W tym przypadku kolejne przyrosty absolutne stają się coraz większe, ponieważ stałe tempo procentowe jest odnoszone do rosnącej podstawy.

Metoda może prowadzić do bardzo wysokich albo bardzo niskich prognoz przy długim horyzoncie. Nawet niewielka, nietypowa zmiana pomiędzy dwoma ostatnimi okresami jest wielokrotnie potęgowana.

Prognoza na podstawie średniego tempa zmian

W zadaniach ze statystyki opisowej i analizy dynamiki spotyka się również prognozę wykorzystującą średnie tempo zmian. Nie opiera się ona wyłącznie na dwóch ostatnich obserwacjach, lecz uwzględnia zmianę zachodzącą w całym analizowanym okresie.

Średni indeks zmian można obliczyć jako średnią geometryczną indeksów łańcuchowych:

\[ \overline{i} = \sqrt[n-1] { i_{2/1}\cdot i_{3/2}\cdot\ldots\cdot i_{n/n-1} } \]

Ponieważ iloczyn indeksów łańcuchowych upraszcza się do ilorazu wartości ostatniej i pierwszej, otrzymujemy:

\[ \overline{i} = \sqrt[n-1] { \frac{y_n}{y_1} } \]

Średnie tempo zmian wynosi:

\[ \overline{g} = \overline{i}-1 \]

Prognoza na jeden okres naprzód może zostać wyznaczona jako:

\[ \widehat{y}_{n+1\mid n} = y_n\overline{i} \]

Dla \(h\) okresów naprzód otrzymujemy:

\[ \widehat{y}_{n+h\mid n} = y_n\overline{i}^{\,h} \]

Metoda ta jest mniej wrażliwa na pojedynczą ostatnią zmianę niż wariant wykorzystujący wyłącznie ostatni indeks łańcuchowy. Nadal zakłada jednak, że średnie tempo zmian z przeszłości będzie utrzymywało się w przyszłości.

Szczegółowe omówienie indeksów indywidualnych, indeksów łańcuchowych, indeksów jednopodstawowych i średniego tempa zmian znajdzie się w osobnym artykule: Indeksy dynamiki i średnie tempo zmian.

Warianty metod naiwnych prognozowania oraz prognoza oparta na średnim tempie zmian
Metody naiwne mogą utrzymywać ostatni poziom, ostatni przyrost lub ostatnie tempo zmian. W dydaktyce spotyka się także prognozę opartą na średnim tempie zmian.

Porównanie metod naiwnych

MetodaZałożeniePostać prognozyPrzebieg prognozy
Poziom bez zmianPrzyszła wartość będzie równa ostatniej obserwacji.\(\widehat{y}_{T+h\mid T}=y_T\)Linia pozioma.
Przyrost bez zmianW przyszłości powtarza się ostatni przyrost absolutny.\(\widehat{y}_{T+h\mid T}=y_T+h(y_T-y_{T-1})\)Linia prosta rosnąca lub malejąca.
Przyrost procentowy bez zmianW przyszłości utrzymuje się ostatnie tempo zmian.\(\widehat{y}_{T+h\mid T}=y_T\left(\frac{y_T}{y_{T-1}}\right)^h\)Ścieżka geometryczna.
Średnie tempo zmianUtrzymuje się przeciętne tempo zmian z całego okresu.\(\widehat{y}_{n+h\mid n}=y_n\overline{i}^{\,h}\)Ścieżka geometryczna oparta na średnim tempie.

Wybór wariantu powinien zależeć od charakteru danych i celu prognozy. Nie należy automatycznie stosować metody przyrostowej tylko dlatego, że dwie ostatnie wartości rosły. Najpierw trzeba sprawdzić, czy ostatnia zmiana nie była przypadkowa lub nietypowa.

Przykład praktyczny – cena cukru

Załóżmy, że przygotowujemy większy model dotyczący kosztów gospodarstwa domowego albo sklepu, a jedną z wykorzystywanych zmiennych jest cena kilograma cukru.

Jeżeli cukier:

to przyjęcie prognozy 5 zł na następny tydzień jest racjonalnym zastosowaniem metody poziomu bez zmian:

\[ \widehat{y}_{T+1\mid T}=y_T=5\text{ zł} \]

Prognoza nie jest pewna. Cena może się zmienić w wyniku promocji, zmiany kosztów hurtowych albo sytuacji rynkowej. Dotychczasowe dane nie wskazują jednak wyraźnego kierunku zmian, dlatego ostatnia wartość stanowi rozsądny punkt odniesienia.

Prostota metody naiwnej nie jest wadą sama w sobie. Wadą byłoby zastosowanie jej bez sprawdzenia, czy przyjęte założenie odpowiada zachowaniu badanego zjawiska.

Prognoza jednookresowa i wielookresowa

Dokładność metody naiwnej może wyglądać zupełnie inaczej w zależności od tego, czy prognozę aktualizujemy po każdej nowej obserwacji, czy wyznaczamy raz całą ścieżkę na kilka okresów naprzód.

Prognozy jednookresowe stale aktualizowane

W podejściu jednookresowym za każdym razem wykorzystujemy najnowszą znaną wartość:

\[ \widehat{y}_{t+1\mid t}=y_t \]

Po poznaniu wartości \(y_{t+1}\) wyznaczamy kolejną prognozę:

\[ \widehat{y}_{t+2\mid t+1}=y_{t+1} \]

Prognoza jest więc ciągle korygowana. W metodzie poziomu bez zmian oznacza to, że model zawsze znajduje się tylko jeden krok za rzeczywistym szeregiem. Przy powolnych zmianach może sprawiać wrażenie bardzo dokładnego.

Prognoza wielookresowa ze stałego momentu

Jeżeli w momencie \(T\) sporządzamy prognozę na kilka okresów naprzód i nie aktualizujemy jej za pomocą nowych obserwacji, dla metody poziomu bez zmian otrzymujemy:

\[ \widehat{y}_{T+1\mid T} = \widehat{y}_{T+2\mid T} = \ldots = \widehat{y}_{T+H\mid T} = y_T \]

Jeżeli rzeczywiste wartości zaczną rosnąć lub spadać, prognoza może szybko oddalać się od przebiegu empirycznego.

Dla metody przyrost bez zmian ścieżka będzie przedłużała ostatnią różnicę:

\[ \widehat{y}_{T+h\mid T} = y_T+h(y_T-y_{T-1}) \]

Jeżeli ostatni przyrost był nietypowy, błąd będzie narastał wraz z horyzontem.

Szersze wyjaśnienie różnicy pomiędzy prognozami aktualizowanymi a prognozami wielookresowymi znajduje się w części Prognozy aktualizowane i prognozy wielookresowe.

Ocena dokładności prognoz naiwnych

Po poznaniu wartości rzeczywistej można obliczyć błąd ex post:

\[ e_{T+h\mid T} = y_{T+h} - \widehat{y}_{T+h\mid T} \]

Dla wielu prognoz można obliczyć między innymi:

Wzory i interpretacje tych miar zostały przedstawione w sekcji Najważniejsze miary błędów prognoz.

Przykład oceny metody poziomu bez zmian

Załóżmy, że wartości szeregu wynoszą:

\[ 100,\quad104,\quad107,\quad111,\quad109 \]

Jednookresowe prognozy naiwne mają postać:

OkresWartość rzeczywista \(y_t\)Prognoza \(\widehat{y}_{t\mid t-1}\)Błąd \(e_t\)\(|e_t|\)
210410044
310710433
411110744
5109111-22

Średni błąd bezwzględny wynosi:

\[ MAE = \frac{4+3+4+2}{4} = 3{,}25 \]

Prognozy różniły się od wartości rzeczywistych przeciętnie o 3,25 jednostki.

Porównuj metody dla tych samych okresów

Jeżeli porównujemy kilka metod, każda z nich powinna być oceniana:

Nie należy porównywać błędu metody aktualizowanej co okres z błędem prognozy wielookresowej bez aktualizacji, nie zaznaczając różnicy pomiędzy tymi procedurami.

Zalety i ograniczenia metod naiwnych

Zalety

Ograniczenia

Podsumowanie

Metody naiwne prognozowania wykorzystują proste założenia dotyczące kontynuacji ostatnio obserwowanych wartości lub zmian.

W metodzie poziom bez zmian przyszła wartość jest równa ostatniej obserwacji:

\[ \widehat{y}_{T+h\mid T}=y_T \]

W metodzie przyrost bez zmian powtarza się ostatni przyrost absolutny:

\[ \widehat{y}_{T+h\mid T} = y_T+h(y_T-y_{T-1}) \]

W metodzie przyrost procentowy bez zmian utrzymuje się ostatnie tempo zmian:

\[ \widehat{y}_{T+h\mid T} = y_T \left( \frac{y_T}{y_{T-1}} \right)^h \]

W zadaniach z indeksów dynamiki stosuje się także wariant oparty na średnim tempie zmian:

\[ \widehat{y}_{n+h\mid n} = y_n\overline{i}^{\,h} \]

Metody naiwne są ważne nie tylko ze względu na prostotę. Pozwalają sprawdzić, czy zastosowanie bardziej skomplikowanego modelu rzeczywiście poprawia dokładność prognoz. Ich wyniki należy jednak oceniać z uwzględnieniem horyzontu oraz sposobu aktualizowania prognozy.

Ruch płaski ciała sztywnego – obliczanie prędkości i przyspieszeń punktów

Ruch płaski ciała sztywnego można traktować jako złożenie ruchu postępowego wybranego punktu oraz obrotu całego ciała względem tego punktu. Taki opis pozwala wyznaczać prędkości i przyspieszenia punktów prętów, tarcz, kół, dźwigni oraz członów mechanizmów.

W artykule przedstawiono dwie podstawowe metody analizy prędkości: zasadę superpozycji oraz metodę chwilowego środka obrotu. Przyspieszenia będą obliczane przede wszystkim za pomocą superpozycji składowych. Omówiony zostanie również chwilowy środek przyspieszeń jako klasyczna, lecz rzadziej stosowana metoda alternatywna.

Materiał rozwija wiadomości z artykułu Kinematyka ciała sztywnego – ruch postępowy, obrotowy i płaski. Przydatne jest również przypomnienie zależności dotyczących ruchu po okręgu.

Na czym polega ruch płaski?

Ruch ciała sztywnego nazywamy płaskim, gdy wszystkie jego punkty poruszają się w płaszczyznach równoległych do jednej ustalonej płaszczyzny. W praktyce analizuje się wtedy przekrój ciała poruszający się w płaszczyźnie xy, a wektory prędkości i przyspieszenia kątowego są skierowane prostopadle do tej płaszczyzny.

\[ \vec\omega=\omega\vec k, \qquad \vec\varepsilon=\varepsilon\vec k. \] 

Położenie ciała w ruchu płaskim można opisać za pomocą dwóch współrzędnych wybranego punktu oraz kąta obrotu:

\[ x_A(t), \qquad y_A(t), \qquad\varphi(t). \]

Oznacza to, że ciało ma trzy stopnie swobody: dwa przesunięcia i jeden obrót.

Wybór punktu odniesienia

Do opisu ruchu wybiera się punkt A, którego prędkość lub przyspieszenie jest znane albo łatwe do wyznaczenia. Może to być środek tarczy, punkt połączenia dwóch członów, suwak poruszający się w prowadnicy albo nieruchomy przegub.

Położenie dowolnego punktu B ciała zapisujemy jako:

\[ \vec r_B=\vec r_A+\vec r_{B/A}. \]

Wektor \(\vec r_{B/A}\) ma stałą długość, ponieważ punkty AB należą do tego samego ciała sztywnego. Może jednak zmieniać kierunek podczas obrotu ciała.

Wybór punktu odniesienia jest dowolny, ale nie obojętny rachunkowo. Najlepiej wybierać punkt o znanym ruchu, ponieważ ogranicza to liczbę niewiadomych w równaniach.

Zasada superpozycji prędkości

Po zróżniczkowaniu równania położenia otrzymujemy podstawową zależność między prędkościami dwóch punktów tego samego ciała:

\[ \boxed{ \vec v_B = \vec v_A + \vec\omega\times\vec r_{B/A} } \]

Równanie to można zapisać krócej:

\[ \vec v_B = \vec v_A+ \vec v_{B/A}, \qquad \vec v_{B/A} = \vec\omega\times\vec r_{B/A}. \]

Składnik \(\vec v_{B/A}\) jest prędkością punktu B względem punktu A, wynikającą z obrotu ciała. W ruchu płaskim:

\[ v_{B/A}=\omega\,AB, \qquad \vec v_{B/A}\perp AB. \] 

Prostopadłość dotyczy prędkości względnej \(\vec v_{B/A}\), a niekoniecznie całkowitej prędkości \(\vec v_B\). Wektor końcowy wyznacza się przez geometryczne dodanie obu składników.

SUPERPOZYCJA PRĘDKOŚCI vB = vA + vB/A A B vA vB/A vB ciało sztywne – odcinek AB vB/A = ω × rB/A Składnik względny jest prostopadły do odcinka AB. Wielobok prędkości vA vB/A vB Wektory dodajemy geometrycznie, a nie przez dodawanie samych wartości.
Superpozycja prędkości punktów ciała sztywnego. Składnik względny \(\vec v_{B/A}\) jest prostopadły do odcinka \(AB\), a prędkość \(\vec v_B\) powstaje przez wektorowe dodanie obu składników. 

Wielobok prędkości

Jeżeli znane są kierunki ruchu punktów, równanie prędkości można rozwiązać graficznie. Z końca znanego wektora \(\vec v_A\) prowadzi się prostą równoległą do kierunku \(\vec v_{B/A}\), czyli prostopadłą do odcinka AB. Z początku wieloboku prowadzi się znany kierunek \(\vec v_B\). Punkt przecięcia zamyka wielobok i pozwala określić poszukiwane wartości.

Wartości prędkości nie zawsze można zwyczajnie dodać. Zależność \(v_B=v_A+v_{B/A}\) jest prawdziwa tylko wtedy, gdy wszystkie trzy wektory mają ten sam kierunek i odpowiednie zwroty. W ogólnym przypadku trzeba wykonywać działania wektorowe.

Chwilowy środek obrotu

Jeżeli prędkość kątowa ciała nie jest równa zeru, jego ruch płaski można w danej chwili interpretować jako obrót wokół punktu I, którego prędkość jest równa zeru:

\[ \vec v_I=0. \]

Punkt I nazywa się chwilowym środkiem obrotu. Dla dowolnego punktu ciała:

\[ \vec v_A = \vec\omega\times\vec r_{A/I}, \qquad v_A=\omega\,IA. \]

Wektor prędkości jest prostopadły do odcinka łączącego rozpatrywany punkt z chwilowym środkiem:

\[ \vec v_A\perp IA, \qquad \vec v_B\perp IB. \] 
WYZNACZANIE CHWILOWEGO ŚRODKA OBROTU Punkt I ma w rozpatrywanej chwili prędkość równą zeru A B I vA vB IA IB vI = 0, vA ⊥ IA, vB ⊥ IB, vA = ω·IA, vB = ω·IB
Chwilowy środek obrotu wyznaczamy jako punkt przecięcia prostych poprowadzonych przez punkty ciała prostopadle do kierunków ich prędkości. 

Jak wyznaczyć chwilowy środek obrotu?

  1. Ustalamy kierunki prędkości dwóch punktów ciała.
  2. Przez pierwszy punkt prowadzimy prostą prostopadłą do jego prędkości.
  3. Przez drugi punkt prowadzimy prostą prostopadłą do jego prędkości.
  4. Punkt przecięcia prostych jest chwilowym środkiem obrotu.

Po wyznaczeniu punktu I można wykorzystać proporcję:

\[ \frac{v_A}{v_B} = \frac{IA}{IB}. \]

W toczeniu bez poślizgu chwilowym środkiem obrotu koła jest punkt styku z nieruchomym podłożem. W ruchu postępowym wszystkie prędkości są równoległe i jednakowe, dlatego chwilowy środek można interpretować jako położony w nieskończoności.

Chwilowy środek obrotu ma zerową prędkość, ale zwykle nie ma zerowego przyspieszenia. Z tego powodu metoda ta jest szczególnie wygodna do obliczania prędkości, lecz nie należy automatycznie przenosić jej na analizę przyspieszeń.

Porównanie metod wyznaczania prędkości

MetodaZaletyOgraniczenia
Superpozycja prędkościUniwersalna, daje pełne równanie wektorowe, dobrze współpracuje z układem współrzędnych.Może wymagać rozpisania składowych i rozwiązania kilku równań.
Chwilowy środek obrotuSzybka i geometrycznie czytelna, szczególnie gdy znane są kierunki prędkości.Bywa niewygodna, gdy punkt I leży daleko poza mechanizmem albo kierunki prędkości są prawie równoległe.

Obie metody opisują ten sam ruch i powinny prowadzić do zgodnych wyników. W praktyce często warto wyznaczyć chwilowy środek w celu szybkiego sprawdzenia kierunków, a obliczenia wykonać równaniem superpozycji.

Zasada superpozycji przyspieszeń

Po zróżniczkowaniu równania prędkości otrzymujemy zależność między przyspieszeniami dwóch punktów tego samego ciała sztywnego:

\[ \boxed{ \vec a_B = \vec a_A + \vec\varepsilon\times\vec r_{B/A} + \vec\omega\times \left( \vec\omega\times\vec r_{B/A} \right) } \]

Równanie można przedstawić jako sumę trzech składników:

\[ \vec a_B = \vec a_A + \vec a_{t,B/A} + \vec a_{n,B/A}. \]

Względna składowa styczna

\[ \vec a_{t,B/A} = \vec\varepsilon\times\vec r_{B/A}, \qquad a_{t,B/A}=\varepsilon\,AB. \]

Składowa styczna jest prostopadła do odcinka AB. Jej zwrot zależy od zwrotu przyspieszenia kątowego.

Względna składowa normalna

\[ \vec a_{n,B/A} = \vec\omega\times \left( \vec\omega\times\vec r_{B/A} \right), \qquad a_{n,B/A}=\omega^2AB. \]

Składowa normalna jest skierowana od punktu B do punktu A. W ruchu płaskim można ją także zapisać jako:

\[ \vec a_{n,B/A} = -\omega^2\vec r_{B/A}. \]
SUPERPOZYCJA PRZYSPIESZEŃ aB = aA + at,B/A + an,B/A A B aA at,B/A an,B/A aB at,B/A = ε·AB,     an,B/A = ω²·AB Składowa styczna jest prostopadła do AB, a normalna jest skierowana od B do A. Wielobok przyspieszeń aA an,B/A at,B/A aB Kolejność składania wektorów jest dowolna, lecz ich kierunki muszą wynikać z geometrii.
Przyspieszenie punktu \(B\) jest sumą przyspieszenia punktu odniesienia \(A\) oraz względnych składowych stycznej i normalnej. 

Wielobok przyspieszeń

Podobnie jak w przypadku prędkości, równanie przyspieszeń można przedstawić za pomocą wieloboku wektorów. Należy jednak uwzględnić dwa względne składniki obrotowe: styczny i normalny. Składowa normalna może być obliczona po wcześniejszym wyznaczeniu \(\omega\), natomiast składowa styczna zawiera zwykle nieznane przyspieszenie kątowe \(\varepsilon\).

Prędkości wyznaczamy przed przyspieszeniami. Wartość \(\omega\) jest potrzebna do obliczenia składowej normalnej \(a_n=\omega^2r\). Próba rozpoczęcia zadania od przyspieszeń prowadzi zwykle do braku danych.

Chwilowy środek przyspieszeń

W klasycznej teorii ruchu płaskiego definiuje się również chwilowy środek przyspieszeń \(Q\), czyli punkt, którego przyspieszenie w rozpatrywanej chwili jest równe zeru:

\[ \vec a_Q=0. \]

Dla dowolnego punktu A można wtedy zapisać:

\[ \vec a_A = \vec\varepsilon\times\vec r_{A/Q} + \vec\omega\times \left( \vec\omega\times\vec r_{A/Q} \right). \]

Przyspieszenie składa się ze składowej stycznej i normalnej względem punktu \(Q\):

\[ a_t=\varepsilon r, \qquad a_n=\omega^2r. \]

Stosunek wartości składowych jest taki sam dla wszystkich punktów ciała:

\[ \frac{a_t}{a_n} = \frac{\varepsilon}{\omega^2}. \]
CHWILOWY ŚRODEK PRZYSPIESZEŃ Metoda alternatywna w analizie ruchu płaskiego A B Q an,A/Q at,A/Q aA an,B/Q at,B/Q aB aQ = 0,     at = εr,     an = ω²r,     tan β = ε/ω² Dla wszystkich punktów ciała kierunek przyspieszenia tworzy taki sam kąt β z kierunkiem do punktu Q. Punkt Q nie musi należeć do fizycznego obszaru ciała. Metoda jest rzadziej stosowana.
Chwilowy środek przyspieszeń \(Q\) jest punktem o zerowym przyspieszeniu. Metoda ma znaczenie teoretyczne i geometryczne, lecz w praktycznych obliczeniach często wygodniejsza jest zasada superpozycji. 

Metoda chwilowego środka przyspieszeń ma znaczenie geometryczne i pojawia się w literaturze z mechaniki technicznej. Jest jednak mniej intuicyjna od chwilowego środka obrotu, ponieważ całkowity wektor przyspieszenia nie jest ani prostopadły, ani równoległy do odcinka łączącego punkt z \(Q\). W praktycznych zadaniach często prostsza i mniej podatna na pomyłki jest zasada superpozycji.

Chwilowy środek obrotu i chwilowy środek przyspieszeń są na ogół różnymi punktami. Punkt o zerowej prędkości nie musi mieć zerowego przyspieszenia, a punkt o zerowym przyspieszeniu nie musi mieć zerowej prędkości.

Prosty przykład z prętem AB

Poziomy pręt ma długość:

\[ AB=0{,}50\ \mathrm{m}. \]

Punkt A porusza się w prawo z prędkością:

\[ \vec v_A=2\vec i\ \mathrm{m/s}. \]

Pręt obraca się przeciwnie do ruchu wskazówek zegara:

\[ \vec\omega=3\vec k\ \mathrm{rad/s}. \]

Wektor położenia punktu B względem A wynosi:

\[ \vec r_{B/A}=0{,}50\vec i\ \mathrm{m}. \]

Prędkość względna:

\[ \vec v_{B/A} = \vec\omega\times\vec r_{B/A} = 3\vec k\times0{,}50\vec i = 1{,}50\vec j\ \mathrm{m/s}. \]

Po zastosowaniu zasady superpozycji:

\[ \vec v_B = 2\vec i+1{,}50\vec j\ \mathrm{m/s}. \]

Wartość prędkości punktu B:

\[ v_B = \sqrt{2^2+1{,}50^2} = 2{,}50\ \mathrm{m/s}. \]
PRZYKŁAD – PRĘDKOŚĆ KOŃCA PRĘTA AB AB = 0,50 m, vA = 2i m/s, ω = 3k rad/s A B 0,50 m vA = 2i vB/A = 1,50j vB ω = 3 rad/s vB/A = ω × rB/A = 1,50j m/s vB = 2i + 1,50j m/s,     |vB| = √(2² + 1,50²) = 2,50 m/s
Prosty przykład zastosowania zasady superpozycji do wyznaczenia prędkości końca poziomego pręta \(AB\).

Krótka kontynuacja dla przyspieszeń

Gdyby dodatkowo znane były \(\vec a_A\) oraz \(\vec\varepsilon\), należałoby obliczyć dwie względne składowe:

\[ \vec a_{t,B/A} = \vec\varepsilon\times\vec r_{B/A}, \qquad \vec a_{n,B/A} = -\omega^2\vec r_{B/A}, \] \[ \vec a_B = \vec a_A + \vec a_{t,B/A} + \vec a_{n,B/A}. \]

W bardziej złożonych mechanizmach trzeba jednocześnie uwzględnić więzy wynikające z prowadnic, przegubów i połączeń kilku członów. Wtedy równania zapisuje się kolejno dla każdego elementu mechanizmu.

Schemat postępowania

  1. Narysuj aktualne położenie mechanizmu i zaznacz jego geometrię.
  2. Określ znane kierunki prędkości wynikające z więzów.
  3. Wybierz punkt odniesienia o znanej prędkości.
  4. Zastosuj równanie superpozycji lub wyznacz chwilowy środek obrotu.
  5. Oblicz prędkość kątową członu.
  6. Dopiero następnie przejdź do przyspieszeń.
  7. Rozdziel przyspieszenie względne na składową styczną i normalną.
  8. Sprawdź kierunki, jednostki oraz zgodność wyniku z warunkami ruchu.
SCHEMAT ROZWIĄZYWANIA ZADAŃRuch płaski ciała sztywnego 1 Narysuj geometrię Zaznacz aktualne położenie członów, podpory, prowadnice i punkty. 2 Określ kierunki prędkości Wykorzystaj więzy: prowadnice, osie obrotu i toczenie bez poślizgu. 3 Wybierz metodę dla prędkości Superpozycja wektorów albo chwilowy środek obrotu. 4 Wyznacz prędkość kątową Skorzystaj z v = ωr lub z równania prędkości dwóch punktów. 5 Oblicz przyspieszenia Zastosuj superpozycję oraz składową styczną i normalną. 6 Sprawdź wynik Skontroluj kierunki, znaki, jednostki i zgodność z więzami.
Zalecana kolejność analizy zadania z prędkości i przyspieszeń w ruchu płaskim ciała sztywnego.

Najczęstsze błędy

Podsumowanie

Podstawowym równaniem prędkości w ruchu płaskim ciała sztywnego jest:

\[ \boxed{ \vec v_B = \vec v_A + \vec\omega\times\vec r_{B/A} } \]

Prędkości można również wyznaczać za pomocą chwilowego środka obrotu:

\[ \vec v_I=0, \qquad v_A=\omega\,IA. \]

Dla przyspieszeń podstawowym narzędziem pozostaje zasada superpozycji:

\[ \boxed{ \vec a_B = \vec a_A + \vec\varepsilon\times\vec r_{B/A} + \vec\omega\times \left( \vec\omega\times\vec r_{B/A} \right) } \]

Chwilowy środek przyspieszeń jest poprawną konstrukcją teoretyczną, jednak ze względu na mniej intuicyjną geometrię w wielu zadaniach wygodniej jest bezpośrednio składać przyspieszenie punktu odniesienia ze względnymi składowymi styczną i normalną.

Indeksy dynamiki, przyrosty i średnie tempo zmian

Indeksy dynamiki, przyrosty absolutne i względne oraz średnie tempo zmian służą do opisywania zmian wartości zjawiska w czasie. Pozwalają ustalić nie tylko, czy wartość wzrosła lub zmalała, lecz również o ile jednostek, o jaki procent albo ile razy zmieniła się w porównaniu z wybranym okresem.

Miary dynamiki można klasyfikować według dwóch niezależnych kryteriów. Pierwszym jest rodzaj mierzonej zmiany, a drugim sposób wyboru okresu stanowiącego podstawę porównania. Z tego powodu zarówno przyrosty absolutne, przyrosty względne, jak i indeksy dynamiki mogą być jednopodstawowe albo łańcuchowe.

W artykule przedstawiamy wzory, interpretacje i zależności między poszczególnymi miarami. Pokazujemy również, jak zmieniać podstawę indeksów jednopodstawowych, przeliczać indeksy jednopodstawowe na łańcuchowe i odwrotnie oraz obliczać średnie tempo zmian.


Spis treści

Szereg czasowy i oznaczenia

Miary dynamiki stosuje się przede wszystkim do analizy szeregów czasowych, czyli ciągów wartości tego samego zjawiska uporządkowanych według kolejnych momentów lub okresów.

Przyjmijmy, że dysponujemy szeregiem:

\[ y_0,\ y_1,\ y_2,\ \ldots,\ y_t,\ \ldots,\ y_n \]

gdzie \(y_t\) oznacza wartość badanego zjawiska w okresie \(t\). Okresy mogą oznaczać lata, kwartały, miesiące, dni albo inne jednostki czasu.

W zależności od sposobu porównywania wartości będziemy stosowali:

Klasyfikacja miar dynamiki

Miary dynamiki można uporządkować według dwóch kryteriów.

Kryterium 1: rodzaj miary

Kryterium 2: sposób przyjęcia podstawy

Określenia „jednopodstawowy” i „łańcuchowy” nie oznaczają osobnego rodzaju miary. Informują jedynie o tym, względem której wartości obliczamy przyrost lub indeks.

Klasyfikacja miar dynamiki według rodzaju miary oraz sposobu przyjęcia podstawy
Miary dynamiki można klasyfikować według rodzaju miary oraz sposobu przyjęcia podstawy: stałej albo zmiennej.

Przyrosty absolutne

Przyrost absolutny informuje, o ile jednostek zmieniła się wartość badanego zjawiska. Jest wyrażony w tych samych jednostkach co analizowana zmienna, na przykład w złotych, kilogramach, osobach albo sztukach.

Przyrost absolutny jednopodstawowy

Jeżeli stałym okresem bazowym jest okres \(0\), przyrost absolutny jednopodstawowy ma postać:

\[ \Delta_t^{(0)}=y_t-y_0 \]

Miara pokazuje, o ile jednostek wartość w okresie \(t\) różni się od wartości w okresie bazowym.

Jeżeli sprzedaż wzrosła z 1000 do 1320 sztuk, przyrost jednopodstawowy wynosi:

\[ \Delta_t^{(0)}=1320-1000=320 \]

Sprzedaż jest więc o 320 sztuk większa niż w okresie bazowym.

Przyrost absolutny łańcuchowy

Przyrost absolutny łańcuchowy porównuje wartość z okresem bezpośrednio poprzedzającym:

\[ \Delta_t=y_t-y_{t-1} \]

Jeżeli sprzedaż wzrosła z 1210 do 1320 sztuk, otrzymujemy:

\[ \Delta_t=1320-1210=110 \]

W badanym okresie sprzedaż wzrosła o 110 sztuk względem okresu poprzedniego.

Dodatni przyrost absolutny oznacza wzrost, przyrost ujemny — spadek, a wartość równa zero — brak zmiany.

Przyrosty względne

Przyrost względny informuje, jaką część wartości odniesienia stanowi zaobserwowana zmiana. Najczęściej przedstawia się go w procentach.

Przyrost względny jednopodstawowy

Dla podstawy w okresie \(0\):

\[ \delta_t^{(0)} = \frac{y_t-y_0}{y_0} \]

W ujęciu procentowym:

\[ \delta_t^{(0)}\cdot100\% = \frac{y_t-y_0}{y_0}\cdot100\% \]

Jeżeli wartość wzrosła z 1000 do 1320, przyrost względny wynosi:

\[ \delta_t^{(0)} = \frac{1320-1000}{1000} = 0{,}32 = 32\% \]

Przyrost względny łańcuchowy

Przyrost względny łańcuchowy odnosi zmianę do wartości z okresu poprzedniego:

\[ \delta_t = \frac{y_t-y_{t-1}}{y_{t-1}} \]

W procentach:

\[ \delta_t\cdot100\% = \frac{y_t-y_{t-1}}{y_{t-1}}\cdot100\% \]

Jeżeli wartość wzrosła z 1210 do 1320, przyrost względny wynosi:

\[ \delta_t = \frac{1320-1210}{1210} \approx 0{,}0909 = 9{,}09\% \]

Indeksy dynamiki

Indeks dynamiki jest ilorazem wartości z okresu badanego i wartości stanowiącej podstawę porównania. Informuje, ile razy wartość w okresie badanym jest większa lub mniejsza od wartości odniesienia.

Indeks jednopodstawowy

Jeżeli okresem bazowym jest okres \(0\), indeks jednopodstawowy ma postać:

\[ I_t^{(0)} = \frac{y_t}{y_0} \]

Wskaźnik procentowy otrzymujemy przez pomnożenie indeksu przez \(100\%\):

\[ I_t^{(0)}\cdot100\% = \frac{y_t}{y_0}\cdot100\% \]

Indeks łańcuchowy

Indeks łańcuchowy porównuje wartość z okresem bezpośrednio poprzedzającym:

\[ I_t = \frac{y_t}{y_{t-1}} \]

W ujęciu procentowym:

\[ I_t\cdot100\% = \frac{y_t}{y_{t-1}}\cdot100\% \]

Interpretacja indeksów dynamiki

Indeks może być przedstawiony jako współczynnik, na przykład \(1{,}15\), albo jako wskaźnik procentowy, na przykład \(115\%\).

Przykładowo indeks:

\[ I_t=1{,}25=125\% \]

oznacza, że wartość w okresie badanym stanowi 125% wartości odniesienia, a zatem jest od niej większa o 25%.

Indeks:

\[ I_t=0{,}82=82\% \]

oznacza, że wartość stanowi 82% wartości odniesienia, czyli jest od niej mniejsza o 18%.

Między indeksem i przyrostem względnym zachodzi prosta zależność:

\[ I_t=1+\delta_t \]

oraz:

\[ \delta_t=I_t-1 \]

Indeks równy 120% nie oznacza wzrostu o 120%. Oznacza wartość równą 120% podstawy, czyli wzrost o 20%.

Wybór podstawy porównania

W indeksach jednopodstawowych podstawą jest bardzo często pierwsza obserwacja szeregu. Nie jest to jednak obowiązkowe. Jako podstawę można wybrać dowolny okres, który ma sens z punktu widzenia analizy.

Podstawą może być na przykład:

Jeżeli podstawą jest okres \(k\), indeks jednopodstawowy zapisujemy jako:

\[ I_t^{(k)} = \frac{y_t}{y_k} \]

Dla okresu bazowego zawsze otrzymujemy:

\[ I_k^{(k)}=1=100\% \]

Zmiana podstawy nie zmienia rzeczywistych wartości szeregu ani przebiegu zjawiska. Zmienia jedynie punkt odniesienia i skalę indeksów.

Przeliczanie indeksów

Zmiana podstawy indeksów jednopodstawowych

Załóżmy, że znamy indeksy względem okresu \(0\):

\[ I_t^{(0)} = \frac{y_t}{y_0} \]

Chcemy przeliczyć je tak, aby nową podstawą był okres \(k\). Nowy indeks ma postać:

\[ I_t^{(k)} = \frac{y_t}{y_k} \]

Dzieląc indeks okresu \(t\) przez indeks nowego okresu bazowego, otrzymujemy:

\[ I_t^{(k)} = \frac{I_t^{(0)}}{I_k^{(0)}} \]

Przy zmianie podstawy indeksów jednopodstawowych wszystkie stare indeksy dzielimy przez indeks okresu, który ma zostać nową podstawą.

Z indeksów jednopodstawowych na łańcuchowe

Indeks łańcuchowy jest ilorazem dwóch kolejnych indeksów jednopodstawowych:

\[ I_t = \frac{I_t^{(0)}}{I_{t-1}^{(0)}} \]

Dlaczego dzielimy? Ponieważ:

\[ \frac{I_t^{(0)}}{I_{t-1}^{(0)}} = \frac{\frac{y_t}{y_0}}{\frac{y_{t-1}}{y_0}} = \frac{y_t}{y_{t-1}} = I_t \]

Z indeksów łańcuchowych na jednopodstawowe

Aby uzyskać indeks jednopodstawowy, mnożymy kolejne indeksy łańcuchowe:

\[ I_t^{(0)} = I_1\cdot I_2\cdot\ldots\cdot I_t \]

W zapisie iloczynowym:

\[ I_t^{(0)} = \prod_{j=1}^{t}I_j \]

Iloczyn upraszcza się, ponieważ kolejne wartości w licznikach i mianownikach skracają się:

\[ \frac{y_1}{y_0} \cdot \frac{y_2}{y_1} \cdot \frac{y_3}{y_2} \cdot\ldots\cdot \frac{y_t}{y_{t-1}} = \frac{y_t}{y_0} \]

Z indeksów jednopodstawowych na łańcuchowe przechodzimy przez dzielenie kolejnych indeksów. Z indeksów łańcuchowych na jednopodstawowe przechodzimy przez mnożenie.

Przeliczanie przyrostów

Przyrosty absolutne łańcuchowe a jednopodstawowe

Jednopodstawowy przyrost absolutny jest sumą kolejnych przyrostów łańcuchowych:

\[ \Delta_t^{(0)} = \Delta_1+\Delta_2+\ldots+\Delta_t \]

Czyli:

\[ \Delta_t^{(0)} = \sum_{j=1}^{t}\Delta_j \]

Wynika to z zależności:

\[ (y_1-y_0)+(y_2-y_1)+\ldots+(y_t-y_{t-1}) = y_t-y_0 \]

Z kolei przyrost łańcuchowy można otrzymać jako różnicę dwóch kolejnych przyrostów jednopodstawowych:

\[ \Delta_t = \Delta_t^{(0)} - \Delta_{t-1}^{(0)} \]

Dlaczego nie dodajemy przyrostów względnych?

Przyrostów względnych i procentowych zmian z kolejnych okresów nie należy zwykle dodawać, ponieważ każda zmiana jest liczona od innej podstawy.

Jeżeli wartość wzrosła najpierw o 10%, a następnie ponownie o 10%, łączny wzrost nie wynosi 20%, lecz:

\[ 1{,}10\cdot1{,}10 = 1{,}21 \]

czyli 21%. W przypadku zmian względnych korzystamy więc z indeksów i mnożenia, a nie z prostego sumowania procentów.

Średnie tempo zmian

Średnie tempo zmian opisuje przeciętną zmianę wartości zjawiska przypadającą na jeden okres. Ponieważ kolejne indeksy dynamiki łączą się przez mnożenie, do obliczeń wykorzystuje się średnią geometryczną, a nie średnią arytmetyczną.

Jeżeli obserwujemy wartości:

\[ y_1,\ y_2,\ \ldots,\ y_n \]

to występuje \(n-1\) zmian pomiędzy kolejnymi obserwacjami. Średni indeks zmian wynosi:

\[ \overline{I} = \sqrt[n-1] { I_2\cdot I_3\cdot\ldots\cdot I_n } \]

Ponieważ iloczyn indeksów łańcuchowych jest równy ilorazowi wartości ostatniej i pierwszej, otrzymujemy:

\[ \overline{I} = \sqrt[n-1] { \frac{y_n}{y_1} } \]

Średni przyrost względny, nazywany także średnim tempem zmian, wynosi:

\[ \overline{\delta} = \overline{I}-1 \]

W procentach:

\[ \overline{\delta}\cdot100\% = (\overline{I}-1)\cdot100\% \]

Jeżeli \(\overline{I}=1{,}08\), oznacza to przeciętny wzrost o 8% w jednym okresie. Jeżeli \(\overline{I}=0{,}97\), oznacza to przeciętny spadek o 3% na okres.

Stopień pierwiastka odpowiada liczbie zmian, a nie liczbie obserwacji. Dla pięciu kolejnych wartości występują cztery zmiany, dlatego stosujemy pierwiastek czwartego stopnia.

Średnie tempo zmian a prognozowanie

Średni indeks dynamiki może zostać wykorzystany do sporządzenia prostej prognozy. Zakładamy wtedy, że przeciętne tempo zmian obserwowane w przeszłości utrzyma się również w przyszłości.

Prognoza na jeden okres naprzód ma postać:

\[ \widehat{y}_{n+1\mid n} = y_n\overline{I} \]

Dla \(h\) okresów naprzód:

\[ \widehat{y}_{n+h\mid n} = y_n\overline{I}^{\,h} \]

Metoda ta została również omówiona w artykule Metody naiwne prognozowania – poziom, przyrost i tempo zmian.

Należy pamiętać, że utrzymanie średniego tempa z przeszłości jest tylko założeniem. Prognoza może być nietrafna, jeżeli zmienią się warunki funkcjonowania badanego zjawiska.

Przykład obliczeniowy

Załóżmy, że sprzedaż produktu w kolejnych latach wynosiła:

RokSprzedaż \(y_t\)
20211000
20221100
20231210
20241331

Przyrosty absolutne łańcuchowe

\[ \Delta_{2022}=1100-1000=100 \]

\[ \Delta_{2023}=1210-1100=110 \]

\[ \Delta_{2024}=1331-1210=121 \]

Indeksy łańcuchowe

\[ I_{2022} = \frac{1100}{1000} = 1{,}10 = 110\% \]

\[ I_{2023} = \frac{1210}{1100} = 1{,}10 = 110\% \]

\[ I_{2024} = \frac{1331}{1210} = 1{,}10 = 110\% \]

W każdym roku sprzedaż rosła o 10% względem roku poprzedniego.

Indeksy jednopodstawowe względem 2021 roku

\[ I_{2021}^{(2021)}=1=100\% \]

\[ I_{2022}^{(2021)} = \frac{1100}{1000} = 1{,}10 = 110\% \]

\[ I_{2023}^{(2021)} = \frac{1210}{1000} = 1{,}21 = 121\% \]

\[ I_{2024}^{(2021)} = \frac{1331}{1000} = 1{,}331 = 133{,}1\% \]

Sprzedaż w 2024 roku stanowiła 133,1% sprzedaży z 2021 roku, czyli była od niej większa o 33,1%.

Średnie tempo zmian

W szeregu występują trzy zmiany, dlatego:

\[ \overline{I} = \sqrt[3] { \frac{1331}{1000} } = \sqrt[3]{1{,}331} = 1{,}10 \]

Średnie tempo zmian wynosi:

\[ \overline{\delta} = 1{,}10-1 = 0{,}10 = 10\% \]

Sprzedaż rosła przeciętnie o 10% rocznie.

Prognoza na 2025 rok

Przy założeniu utrzymania średniego tempa zmian:

\[ \widehat{y}_{2025\mid2024} = 1331\cdot1{,}10 = 1464{,}1 \]

Prognozowana sprzedaż wynosi około 1464 sztuki.

Najczęstsze błędy

Podsumowanie

Przyrosty absolutne pokazują zmianę w jednostkach badanego zjawiska, przyrosty względne — zmianę w stosunku do podstawy, a indeksy dynamiki — relację wartości z dwóch okresów.

Każda z tych miar może być:

Przy zmianie podstawy indeksów jednopodstawowych dzielimy każdy indeks przez indeks nowego okresu bazowego:

\[ I_t^{(k)} = \frac{I_t^{(0)}}{I_k^{(0)}} \]

Indeksy łańcuchowe uzyskujemy przez dzielenie kolejnych indeksów jednopodstawowych:

\[ I_t = \frac{I_t^{(0)}}{I_{t-1}^{(0)}} \]

Indeksy jednopodstawowe odtwarzamy natomiast przez mnożenie indeksów łańcuchowych:

\[ I_t^{(0)} = \prod_{j=1}^{t}I_j \]

Średnie tempo zmian opiera się na średniej geometrycznej i opisuje przeciętną zmianę przypadającą na jeden okres:

\[ \overline{\delta} = \sqrt[n-1] { \frac{y_n}{y_1} } -1 \]

Poprawna interpretacja wymaga zawsze wskazania rodzaju miary, sposobu przyjęcia podstawy oraz okresów, których dotyczy porównanie.

Prognozowanie metodą średniej ruchomej – średnia prosta i ważona

Średnia ruchoma jest jedną z najprostszych metod analizy i prognozowania szeregów czasowych. Jej wartość oblicza się na podstawie określonej liczby ostatnich obserwacji, tworzących tak zwane okno średniej ruchomej. Po pojawieniu się nowej obserwacji okno przesuwa się o jeden okres, najstarsza wartość zostaje pominięta, a do obliczeń włączana jest wartość najnowsza.

Średnie ruchome mogą służyć zarówno do wygładzania szeregu czasowego, jak i do sporządzania krótkookresowych prognoz. Pozwalają ograniczyć wpływ przypadkowych wahań, ułatwiają rozpoznanie tendencji rozwojowej oraz mogą być wykorzystywane do wyznaczania prognoz na kolejny okres.

W artykule omawiamy prostą średnią ruchomą, ważoną średnią ruchomą, wpływ długości okna na wyniki, zastosowanie średnich do wygładzania sezonowości oraz sposoby wyznaczania prognoz wielookresowych po zakończeniu zbioru obserwacji empirycznych.


Spis treści

Czym jest średnia ruchoma?

Średnia ruchoma jest średnią obliczaną dla kolejnych, przesuwających się fragmentów szeregu czasowego. Każdy taki fragment obejmuje ustaloną liczbę obserwacji, nazywaną długością okna.

Jeżeli okno obejmuje \(k\) obserwacji, w okresie \(t\) wykorzystujemy wartości:

\[ y_t,\ y_{t-1},\ y_{t-2},\ldots,y_{t-k+1} \]

Po przejściu do następnego okresu najstarsza obserwacja wypada z okna, a do obliczeń włączana jest nowa wartość. Okno „porusza się” więc wraz z upływem czasu.

Podstawowe informacje o szeregach czasowych oraz ich budowie zostały przedstawione w artykule poświęconym szeregom statystycznym.

Średnia ruchoma nie wykorzystuje jednego stałego zestawu danych. Jej okno jest regularnie przesuwane, dlatego każda kolejna średnia uwzględnia nowszy fragment szeregu.

Prosta średnia ruchoma

W prostej średniej ruchomej każda obserwacja znajdująca się w oknie ma jednakowe znaczenie. Dla okna długości \(k\) średnia ruchoma w okresie \(t\) wynosi:

\[ MA_t(k) = \frac{1}{k} \sum_{j=0}^{k-1}y_{t-j} \]

W zapisie rozwiniętym:

\[ MA_t(k) = \frac{ y_t+y_{t-1}+y_{t-2}+\ldots+y_{t-k+1} }{k} \]

Jest to, po prostu, średnia arytmetyczna \(k\) obserwacji bezpośrednio poprzedzających okres, dla którego liczymy średnią ruchomą.

Jeżeli stosujemy średnią czterookresową, otrzymujemy:

\[ MA_t(4) = \frac{ y_t+y_{t-1}+y_{t-2}+y_{t-3} }{4} \]

Przykładowo dla wartości:

\[ 98,\quad102,\quad105,\quad103 \]

średnia wynosi:

\[ MA_t(4) = \frac{98+102+105+103}{4} = 102 \]

Jeżeli średnią wykorzystujemy do prognozowania, wartość \(102\) może stanowić prognozę na kolejny okres.

Związek z metodą naiwną

Prosta metoda naiwna poziomu bez zmian jest szczególnym przypadkiem średniej ruchomej. Jeżeli długość okna zmniejszymy do jednej obserwacji, otrzymujemy:

\[ MA_t(1)=y_t \]

Prognoza na następny okres wynosi wtedy:

\[ \widehat{y}_{t+1\mid t}=y_t \]

Jest to dokładnie metoda naiwna poziomu bez zmian.

Metoda naiwna poziomu bez zmian odpowiada średniej ruchomej z oknem \(k=1\).

Jak dobrać długość okna?

Długość okna jest jednym z najważniejszych parametrów metody. Decyduje o tym, jak szybko średnia reaguje na nowe informacje i jak silnie wygładza przypadkowe wahania.

Krótkie okno

Średnia ruchoma z krótkim oknem:

Średnia trzyokresowa \(MA(3)\) zwykle szybciej dostosowuje się do zmiany poziomu szeregu niż średnia pięcio- lub siedmiookresowa.

Długie okno

Średnia ruchoma z dłuższym oknem:

Porównanie danych rzeczywistych, średniej ruchomej MA(3), średniej ruchomej MA(5) oraz ważonej średniej ruchomej
Krótsze okno średniej ruchomej szybciej reaguje na zmiany, natomiast dłuższe silniej wygładza szereg. Ważona średnia ruchoma nadaje większe znaczenie nowszym obserwacjom.

W szeregu rosnącym średnia o długim oknie może być systematycznie niższa od bieżących obserwacji. W szeregu malejącym może być od nich wyższa. Zjawisko to jest określane jako opóźnienie średniej ruchomej.

Dlaczego dla pierwszych okresów nie ma średniej?

Do obliczenia średniej ruchomej potrzebna jest pełna liczba obserwacji określona przez długość okna. Jeżeli stosujemy średnią \(k\)-okresową, pierwszą wartość można wyznaczyć dopiero po zgromadzeniu \(k\) obserwacji.

Dla średniej \(MA(3)\) pierwszą wartość wyznaczamy z obserwacji:

\[ y_1,\quad y_2,\quad y_3 \]

Jeżeli średnią przypisujemy do końca okna, pierwsza średnia pojawia się w okresie \(3\). Dla okresów \(1\) i \(2\) brakuje dostatecznej liczby wcześniejszych danych.

Dla średniej \(MA(5)\) pierwsza wartość pojawia się dopiero w okresie \(5\). Oznacza to, że dla czterech pierwszych okresów średnia nie może zostać wyznaczona.

Ogólnie dla średniej o oknie \(k\) pomijamy na początku:

\[ k-1 \]

obserwacji, jeżeli średnia jest wyrównana do końca okna.

Ma to znaczenie przy porównywaniu błędów różnych metod. Średnia o dłuższym oknie tworzy mniej historycznych prognoz jednookresowych, ponieważ część początkowych obserwacji służy jedynie do zbudowania pierwszego pełnego okna.

Dla porównania:

Przy porównywaniu dokładności modeli należy wykorzystywać wspólny zakres okresów, dla których każda z ocenianych metod może wyznaczyć prognozę.

Ważona średnia ruchoma

W prostej średniej ruchomej wszystkie obserwacje mają jednakowy wpływ na wynik. W ważonej średniej ruchomej poszczególnym obserwacjom przypisuje się różne wagi.

Dla okna długości \(k\) ważona średnia ruchoma ma postać:

\[ WMA_t = \sum_{j=0}^{k-1}w_jy_{t-j} \]

W zapisie rozwiniętym:

\[ WMA_t = w_0y_t+w_1y_{t-1}+\ldots+w_{k-1}y_{t-k+1} \]

Wagi powinny spełniać warunek:

\[ \sum_{j=0}^{k-1}w_j=1 \]

Najczęściej przyjmuje się również:

\[ w_j\geq0 \]

Dzięki temu ważona średnia pozostaje kombinacją wartości znajdujących się w oknie.

Przykład

Załóżmy, że wykorzystujemy trzy ostatnie obserwacje, uporządkowane od najstarszej do najnowszej:

\[ 104,\quad112,\quad118 \]

Przypisujemy im wagi:

\[ 0{,}2,\quad0{,}3,\quad0{,}5 \]

Suma wag wynosi:

\[ 0{,}2+0{,}3+0{,}5=1 \]

Ważona średnia ruchoma wynosi:

\[ WMA_t = 0{,}2\cdot104 + 0{,}3\cdot112 + 0{,}5\cdot118 \]

\[ WMA_t = 20{,}8+33{,}6+59 = 113{,}4 \]

Najnowsza obserwacja ma największą wagę, dlatego wynik jest bardziej zbliżony do wartości \(118\) niż w przypadku zwykłej średniej arytmetycznej.

Zasada postarzania informacji

W prognozowaniu często zakłada się, że najnowsze obserwacje zawierają więcej informacji o najbliższej przyszłości niż obserwacje starsze. Z tego powodu w ważonej średniej ruchomej stosuje się zasadę postarzania informacji.

Wagi są wtedy uporządkowane następująco:

\[ w_0>w_1>w_2>\ldots>w_{k-1} \]

gdzie \(w_0\) jest wagą obserwacji najnowszej, a \(w_{k-1}\) wagą obserwacji najstarszej.

Przykładowy układ wag dla średniej czterookresowej może mieć postać:

\[ 0{,}4,\quad0{,}3,\quad0{,}2,\quad0{,}1 \]

Najświeższa obserwacja otrzymuje wagę \(0{,}4\), a najstarsza wagę \(0{,}1\).

Malejące wagi są często uzasadnione, ale nie stanowią bezwzględnego wymogu matematycznego. Dobór wag powinien odpowiadać charakterowi zjawiska. W niektórych zastosowaniach większe znaczenie mogą mieć obserwacje pochodzące z analogicznych okresów sezonowych.

Średnia ruchoma a sezonowość

Średnie ruchome są często wykorzystywane do wygładzania wahań sezonowych. Jeżeli długość okna odpowiada długości pełnego cyklu sezonowego, wartości wysokie i niskie występujące w różnych częściach cyklu mogą się wzajemnie równoważyć.

Typowe długości okna to:

Średnia czterokwartalna może ograniczyć regularne różnice pomiędzy poszczególnymi kwartałami i ułatwić rozpoznanie długookresowego trendu.

W tym zastosowaniu średnia ruchoma jest przede wszystkim narzędziem wygładzania i dekompozycji szeregu. Nie każda średnia służąca do wygładzania jest jednocześnie bezpośrednią prognozą przyszłej wartości.

Centrowana średnia ruchoma

Jeżeli długość okna jest nieparzysta, środek okna przypada na konkretny okres. Dla średniej trzyokresowej:

\[ \frac{y_{t-1}+y_t+y_{t+1}}{3} \]

wartość można przypisać do okresu \(t\).

Przy oknie parzystym, na przykład czterookresowym, środek okna znajduje się pomiędzy dwoma okresami. W analizie szeregu stosuje się wtedy centrowanie średniej ruchomej.

Najpierw oblicza się kolejne średnie czterookresowe:

\[ MA_{t-\frac12}(4), \qquad MA_{t+\frac12}(4) \]

Następnie oblicza się ich średnią:

\[ CMA_t(4) = \frac{ MA_{t-\frac12}(4)+MA_{t+\frac12}(4) }{2} \]

Centrowana średnia ruchoma jest szczególnie przydatna podczas wyodrębniania trendu i składnika sezonowego. Nie jest jednak klasyczną prognozą na końcu szeregu, ponieważ do jej obliczenia potrzebne są obserwacje leżące zarówno przed, jak i po analizowanym okresie.

Prognoza na jeden okres naprzód

Najbardziej naturalnym zastosowaniem średniej ruchomej jest prognoza na jeden okres naprzód. Dla prostej średniej o oknie \(k\):

\[ \widehat{y}_{T+1\mid T} = \frac{1}{k} \sum_{j=0}^{k-1}y_{T-j} \]

Dla ważonej średniej ruchomej:

\[ \widehat{y}_{T+1\mid T} = \sum_{j=0}^{k-1}w_jy_{T-j} \]

Po poznaniu rzeczywistej wartości \(y_{T+1}\) przesuwamy okno o jeden okres i wyznaczamy prognozę na okres \(T+2\).

Jest to przykład prognozowania stale aktualizowanego. Każda nowa prognoza korzysta z najnowszej dostępnej obserwacji.

Różnicę pomiędzy prognozami stale aktualizowanymi a prognozami rozwijanymi ze stałego momentu przedstawiono w artykule Prognozowanie – podstawowe pojęcia, metody i błędy prognoz.

Prognoza na kilka okresów naprzód

Po zakończeniu szeregu obserwacji pojawia się problem: aby przesunąć okno dalej, potrzebne są wartości, których jeszcze nie znamy. Można wtedy zastosować kilka procedur.

Stała prognoza dla wszystkich horyzontów

Najprostszym rozwiązaniem jest utrzymanie pierwszej wyznaczonej prognozy dla wszystkich kolejnych okresów:

\[ \widehat{y}_{T+h\mid T} = \widehat{y}_{T+1\mid T} \]

Ścieżka prognozy jest wtedy pozioma. Podejście jest proste, ale nie wykorzystuje konstrukcji przesuwającego się okna w kolejnych horyzontach.

Prognozowanie rekurencyjne

Drugim podejściem jest zastępowanie brakujących przyszłych obserwacji wcześniej wyznaczonymi prognozami. Procedurę tę można określić jako prognozowanie rekurencyjne albo iteracyjne.

Dla średniej trzyokresowej pierwsza prognoza wynosi:

\[ \widehat{y}_{T+1\mid T} = \frac{ y_T+y_{T-1}+y_{T-2} }{3} \]

W kolejnym kroku nie znamy wartości \(y_{T+1}\), dlatego zastępujemy ją prognozą:

\[ \widehat{y}_{T+2\mid T} = \frac{ \widehat{y}_{T+1\mid T} + y_T + y_{T-1} }{3} \]

Następnie:

\[ \widehat{y}_{T+3\mid T} = \frac{ \widehat{y}_{T+2\mid T} + \widehat{y}_{T+1\mid T} + y_T }{3} \]

W kolejnych krokach rzeczywiste obserwacje stopniowo wypadają z okna, a ich miejsce zajmują wartości prognozowane.

Rekurencyjne prognozy średniej ruchomej często stopniowo stabilizują się w pobliżu pewnego poziomu. Wraz z horyzontem maleje wpływ rzeczywistych obserwacji, a rośnie udział wcześniejszych prognoz.

Przy prognozowaniu wielookresowym trzeba zawsze wskazać, czy kolejne prognozy są utrzymywane na stałym poziomie, czy też są obliczane rekurencyjnie przez podstawianie wcześniejszych prognoz.

Przykład obliczeniowy

Załóżmy, że wartości szeregu wynoszą:

\[ 100,\quad106,\quad104,\quad112,\quad118 \]

Prosta średnia ruchoma \(MA(3)\)

Prognoza na okres \(6\) opiera się na trzech ostatnich obserwacjach:

\[ \widehat{y}_{6\mid5} = \frac{104+112+118}{3} = \frac{334}{3} \approx 111{,}33 \]

W prognozowaniu rekurencyjnym wartość ta zastępuje nieznaną obserwację z okresu \(6\).

Prognoza na okres \(7\):

\[ \widehat{y}_{7\mid5} = \frac{ 112+118+111{,}33 }{3} \approx 113{,}78 \]

Prognoza na okres \(8\):

\[ \widehat{y}_{8\mid5} = \frac{ 118+111{,}33+113{,}78 }{3} \approx 114{,}37 \]

Ważona średnia ruchoma

Przypiszmy trzem ostatnim obserwacjom wagi:

\[ 0{,}2,\quad0{,}3,\quad0{,}5 \]

Największa waga dotyczy wartości najnowszej. Prognoza wynosi:

\[ \widehat{y}_{6\mid5} = 0{,}2\cdot104 + 0{,}3\cdot112 + 0{,}5\cdot118 \]

\[ \widehat{y}_{6\mid5} = 20{,}8+33{,}6+59 = 113{,}4 \]

Ważona średnia daje wyższą prognozę niż prosta średnia \(MA(3)\), ponieważ najnowsza i najwyższa obserwacja otrzymała największą wagę.

Ocena dokładności prognoz

Dokładność średnich ruchomych można oceniać za pomocą błędów ex post i typowych miar, takich jak:

Wzory i interpretacje tych miar znajdują się w części Najważniejsze miary błędów prognoz.

Porównując kilka długości okna, należy używać tych samych okresów testowych. Jeżeli \(MA(5)\) nie może wyznaczyć prognoz dla pierwszych czterech obserwacji, nie należy porównywać jej błędu z \(MA(3)\) na innym, dłuższym zakresie.

Dobór okna może zostać przeprowadzony przez obliczenie błędów dla kilku wariantów, na przykład:

Następnie wybiera się wariant tworzący najmniejsze błędy na danych, które nie były wykorzystywane do arbitralnego dopasowania długości okna.

Zalety i ograniczenia

Zalety średnich ruchomych

Ograniczenia średnich ruchomych

Ograniczenia średnich ruchomych prowadzą do zastosowania metod adaptacyjnych, w których wpływ starszych obserwacji zmniejsza się w bardziej systematyczny sposób. Należą do nich między innymi proste wygładzanie wykładnicze, metoda Browna, metoda Holta oraz metoda Holta – Wintersa.

Podsumowanie

Prosta średnia ruchoma o oknie \(k\) jest średnią arytmetyczną z \(k\) ostatnich obserwacji:

\[ MA_t(k) = \frac{1}{k} \sum_{j=0}^{k-1}y_{t-j} \]

Metoda naiwna poziomu bez zmian jest jej szczególnym przypadkiem dla \(k=1\).

Krótkie okno szybciej reaguje na nowe informacje, lecz słabiej wygładza przypadkowe wahania. Długie okno daje bardziej stabilny przebieg, ale powoduje większe opóźnienie.

Ważona średnia ruchoma pozwala różnicować znaczenie obserwacji:

\[ WMA_t = \sum_{j=0}^{k-1}w_jy_{t-j}, \qquad \sum_{j=0}^{k-1}w_j=1 \]

W prognozowaniu często nadaje się większe wagi obserwacjom nowszym, zgodnie z zasadą postarzania informacji.

Po zakończeniu szeregu prognozę można utrzymać na stałym poziomie albo rozwijać rekurencyjnie, zastępując nieznane obserwacje wcześniejszymi prognozami. Wybrana procedura powinna być jasno wskazana, ponieważ wpływa na przebieg i błędy prognozy wielookresowej.

Średnie ruchome są prostymi i użytecznymi narzędziami, ale ich wyniki zależą od długości okna, sposobu doboru wag oraz charakteru badanego szeregu czasowego.

Wygładzanie wykładnicze – metody Browna, Holta i Wintersa

Wygładzanie wykładnicze obejmuje grupę adaptacyjnych metod prognozowania szeregów czasowych. Ich podstawową cechą jest regularne aktualizowanie modelu po pojawieniu się każdej nowej obserwacji. Najnowsze dane otrzymują zwykle największe znaczenie, natomiast wpływ obserwacji starszych stopniowo maleje.

Do najczęściej omawianych metod należą metoda Browna, metoda Holta oraz metoda Wintersa, nazywana również metodą Holta – Wintersa. Metodę Browna stosuje się przede wszystkim do szeregów o względnie stałym poziomie, metodę Holta do szeregów z trendem, a metodę Wintersa do szeregów zawierających jednocześnie trend i sezonowość.

W artykule przyjmujemy konsekwentną konwencję, w której:

Osobno zapisujemy szeregi pomocnicze modelu i osobno wynikającą z nich prognozę. Takie podejście jest bardziej czytelne niż bezpośrednie włączanie prognozy do równania aktualizującego model.


Spis treści

Idea wygładzania wykładniczego

Wygładzanie wykładnicze należy do metod adaptacyjnych. Oznacza to, że wartości opisujące model są aktualizowane za każdym razem, gdy pojawia się nowa obserwacja.

W odróżnieniu od średniej ruchomej model nie ogranicza się do jednego, sztywnego okna ostatnich \(k\) obserwacji. Pośrednio uwzględnia całą historię szeregu, ale znaczenie starszych danych maleje wraz z ich wiekiem.

Wartość wygładzona jest połączeniem:

Najprostsze równanie ma postać:

\[ F_t = \alpha y_t + (1-\alpha)F_{t-1} \]

gdzie:

Wykładniczo malejące wagi

Nazwa metody wynika ze sposobu, w jaki maleje wpływ kolejnych obserwacji. Po wielokrotnym podstawieniu wcześniejszych wartości wygładzonych otrzymujemy:

\[ F_t = \alpha y_t + \alpha(1-\alpha)y_{t-1} + \alpha(1-\alpha)^2y_{t-2} +\ldots \]

Najnowsza obserwacja otrzymuje wagę:

\[ \alpha \]

obserwacja poprzednia:

\[ \alpha(1-\alpha) \]

a obserwacja sprzed dwóch okresów:

\[ \alpha(1-\alpha)^2 \]

Wagi tworzą więc ciąg geometryczny. Starsze informacje nadal wpływają na wynik, ale ich znaczenie stopniowo maleje.

Wygładzanie wykładnicze nie odrzuca nagle starszych obserwacji, jak średnia ruchoma po przesunięciu okna. Ich wpływ zanika stopniowo.

Brown, Holt i Winters – porównanie

Wybór metody zależy przede wszystkim od budowy szeregu czasowego.

MetodaPoziomTrendSezonowośćParametry
Browntaknienie\(\alpha\)
Holttaktaknie\(\alpha,\beta\)
Winterstaktaktak\(\alpha,\beta,\gamma\)

Można zatem przedstawić rozwój modeli w następujący sposób:

\[ \text{Brown} \rightarrow \text{Holt} \rightarrow \text{Winters} \]

Każdy kolejny model dodaje nowy element struktury szeregu. 

Porównanie metod Browna, Holta i Wintersa w wygładzaniu wykładniczym
Metoda Browna służy do prognozowania szeregu bez trendu i sezonowości, metoda Holta uwzględnia trend, a metoda Wintersa pozwala modelować jednocześnie trend i sezonowość.

Metoda Browna

Najprostsza metoda wygładzania wykładniczego jest w wielu polskojęzycznych materiałach nazywana metodą Browna. Stosuje się ją do szeregów, w których nie występuje wyraźny trend ani sezonowość, a wartości oscylują wokół względnie stałego poziomu.

W literaturze można spotkać różne konwencje nazewnicze. Określenie „metoda Browna” bywa również odnoszone do podwójnego wygładzania wykładniczego. Dlatego w tym artykule przez metodę Browna rozumiemy proste wygładzanie wykładnicze szeregu bez trendu i sezonowości.

Wartość wygładzona

Szereg pomocniczy wartości wygładzonych obliczamy według wzoru:

\[ F_t = \alpha y_t + (1-\alpha)F_{t-1} \]

Prognoza

Wartość wygładzona z okresu \(t\) staje się prognozą na okres następny:

\[ \widehat{y}_{t+1\mid t}=F_t \]

Prognoza wyznaczona w poprzednim okresie i przypisana do okresu \(t\) wynosi zatem:

\[ \widehat{y}_{t\mid t-1}=F_{t-1} \]

W tabeli obliczeniowej warto rozdzielić kolumny:

\[ t,\qquad y_t,\qquad F_t,\qquad \widehat{y}_{t\mid t-1},\qquad e_t \]

Dzięki temu łatwo zauważyć, że wartość \(F_t\) nie jest prognozą dla tego samego okresu, lecz prognozą na okres kolejny.

Wartość początkowa w metodzie Browna

Do rozpoczęcia obliczeń potrzebna jest wartość początkowa \(F_1\). Najprostsza konwencja polega na przyjęciu:

\[ F_1=y_1 \]

Można również przyjąć średnią arytmetyczną z pierwszych kilku obserwacji:

\[ F_1 = \frac{y_1+y_2+\ldots+y_r}{r} \]

Drugie rozwiązanie może ograniczyć wpływ nietypowej pierwszej obserwacji, ale powoduje, że właściwe obliczenia rozpoczynają się później.

Możliwe jest również potraktowanie wartości początkowej jako dodatkowego parametru wyznaczanego numerycznie. W przykładzie dydaktycznym przyjmiemy jednak najprostszą konwencję:

\[ F_1=y_1 \]

Metoda Holta

Metoda Holta rozszerza proste wygładzanie wykładnicze o składnik trendu. Tworzy dwa szeregi pomocnicze:

Aktualizacja poziomu

\[ F_t = \alpha y_t + (1-\alpha)(F_{t-1}+S_{t-1}) \]

Nowy poziom jest kombinacją najnowszej obserwacji oraz poziomu z poprzedniego okresu powiększonego o dotychczasowy trend.

Aktualizacja trendu

\[ S_t = \beta(F_t-F_{t-1}) + (1-\beta)S_{t-1} \]

Nowy składnik trendu jest kombinacją najnowszej zmiany poziomu i wcześniejszego trendu.

Prognoza

Prognoza na jeden okres naprzód wynosi:

\[ \widehat{y}_{t+1\mid t} = F_t+S_t \]

Prognoza na \(h\) okresów naprzód:

\[ \widehat{y}_{t+h\mid t} = F_t+hS_t \]

W tabeli dla okresu \(t\) obliczamy \(F_t\) i \(S_t\), natomiast suma \(F_t+S_t\) jest prognozą dotyczącą okresu \(t+1\).

Przy zapisie prognozy w wierszu okresu, którego dotyczy, otrzymujemy:

\[ \widehat{y}_{t\mid t-1} = F_{t-1}+S_{t-1} \]

Wartości początkowe w metodzie Holta

Metoda Holta wymaga ustalenia początkowego poziomu i trendu. W literaturze i programach komputerowych stosuje się różne sposoby inicjalizacji.

Najprostsza konwencja

Poziom początkowy:

\[ F_1=y_1 \]

Trend początkowy:

\[ S_1=y_2-y_1 \]

Jest to konwencja najłatwiejsza do zastosowania w przykładach rachunkowych.

Średni początkowy przyrost

Aby ograniczyć wpływ przypadkowej różnicy pomiędzy dwiema pierwszymi obserwacjami, można wykorzystać średnią z kilku początkowych przyrostów:

\[ S_1 = \frac{ (y_2-y_1)+(y_3-y_2)+\ldots+(y_r-y_{r-1}) }{r-1} \]

Po skróceniu:

\[ S_1 = \frac{y_r-y_1}{r-1} \]

Trend wyznaczony za pomocą regresji

Na początkowym fragmencie szeregu można oszacować model ekonometryczny trendu liniowego:

\[ y_t=a+bt+\varepsilon_t \]

Współczynnik kierunkowy \(b\) może zostać przyjęty jako wartość początkowa trendu:

\[ S_1=b \]

Takie rozwiązanie jest często stabilniejsze od pojedynczej różnicy, ale wymaga dodatkowych obliczeń.

Metoda Wintersa

Metoda Wintersa, nazywana również metodą Holta – Wintersa, rozszerza model Holta o składnik sezonowy.

Tworzymy trzy szeregi pomocnicze:

Załóżmy, że pełny cykl sezonowy obejmuje \(m\) okresów. Dla danych kwartalnych z cyklem rocznym \(m=4\), a dla danych miesięcznych \(m=12\).

Sezonowość addytywna i multiplikatywna

Metoda Wintersa występuje w dwóch podstawowych wariantach. Wybór zależy od tego, jak zmienia się amplituda wahań sezonowych.

Wariant addytywny

Model addytywny stosujemy, gdy wahania sezonowe mają zbliżoną amplitudę niezależnie od poziomu szeregu. Sezonowość jest wyrażana w tych samych jednostkach co analizowana zmienna.

Ogólny zapis szeregu można przedstawić jako:

\[ y_t = \text{poziom} + \text{trend} + \text{sezonowość} + \text{składnik losowy} \]

Równania modelu addytywnego:

\[ F_t = \alpha(y_t-C_{t-m}) + (1-\alpha)(F_{t-1}+S_{t-1}) \]

\[ S_t = \beta(F_t-F_{t-1}) + (1-\beta)S_{t-1} \]

\[ C_t = \gamma(y_t-F_t) + (1-\gamma)C_{t-m} \]

Prognoza na \(h\) okresów naprzód:

\[ \widehat{y}_{t+h\mid t} = F_t+hS_t+C_{t-m+r} \]

gdzie \(r\) wskazuje składnik sezonowy właściwy dla prognozowanego okresu. W praktyce indeks sezonowy dobiera się cyklicznie, zgodnie z położeniem prognozy w sezonie.

Wariant multiplikatywny

Model multiplikatywny stosujemy, gdy amplituda wahań sezonowych rośnie lub maleje wraz z poziomem szeregu. Sezonowość jest wtedy wyrażana w postaci mnożnika.

Ogólny zapis:

\[ y_t = (\text{poziom}+\text{trend}) \cdot \text{sezonowość} \cdot \text{składnik losowy} \]

Równania modelu multiplikatywnego:

\[ F_t = \alpha\frac{y_t}{C_{t-m}} + (1-\alpha)(F_{t-1}+S_{t-1}) \]

\[ S_t = \beta(F_t-F_{t-1}) + (1-\beta)S_{t-1} \]

\[ C_t = \gamma\frac{y_t}{F_t} + (1-\gamma)C_{t-m} \]

Prognoza:

\[
\widehat{y}_{t+h\mid t}
=
(F_t+hS_t)C_{t-m+r}
\]

Porównanie sezonowości addytywnej i multiplikatywnej na wykresach szeregu czasowego
W sezonowości addytywnej amplituda wahań pozostaje zbliżona, natomiast w sezonowości multiplikatywnej zmienia się wraz z poziomem szeregu.

Stała amplituda wahań wskazuje zwykle na wariant addytywny. Amplituda zwiększająca się wraz z poziomem szeregu wskazuje na wariant multiplikatywny.

Wartości początkowe w metodzie Wintersa

Inicjalizacja metody Wintersa jest bardziej złożona niż w metodach Browna i Holta. Potrzebujemy początkowego poziomu, trendu oraz zestawu \(m\) początkowych składników sezonowych.

Stosowane są różne konwencje. Programy statystyczne mogą estymować wartości początkowe numerycznie, korzystać z procedur heurystycznych albo przyjmować wartości wskazane przez użytkownika.

Początkowy poziom

Najprostszą wartością początkową jest średnia arytmetyczna z pierwszego pełnego cyklu:

\[ F_m = \frac{1}{m} \sum_{t=1}^{m}y_t \]

Początkowy trend

Jeżeli dysponujemy co najmniej dwoma pełnymi cyklami, możemy obliczyć średnie:

\[ \overline{y}_1 = \frac{1}{m} \sum_{t=1}^{m}y_t \]

\[ \overline{y}_2 = \frac{1}{m} \sum_{t=m+1}^{2m}y_t \]

Trend przypadający na jeden okres można przyjąć jako:

\[ S_m = \frac{\overline{y}_2-\overline{y}_1}{m} \]

Możliwe jest również oszacowanie trendu liniowego na początkowym fragmencie szeregu i wykorzystanie jego współczynnika kierunkowego.

Początkowe składniki sezonowe – wariant addytywny

Dla pierwszego cyklu:

\[ C_t=y_t-F_m, \qquad t=1,\ldots,m \]

Składniki sezonowe powinny przeciętnie dawać zero:

\[ \sum_{t=1}^{m}C_t=0 \]

Początkowe składniki sezonowe – wariant multiplikatywny

W wariancie multiplikatywnym:

\[ C_t = \frac{y_t}{F_m}, \qquad t=1,\ldots,m \]

Ich średnia powinna wynosić jeden:

\[ \frac{1}{m} \sum_{t=1}^{m}C_t=1 \]

czyli:

\[ \sum_{t=1}^{m}C_t=m \]

W artykule i przykładach możemy przyjmować tę prostą konwencję, zaznaczając, że inne podręczniki i programy mogą wykorzystywać odmienne procedury inicjalizacji.

Znaczenie parametrów alfa, beta i gamma

Parametry wygładzania spełniają warunki:

\[ 0\leq\alpha\leq1, \qquad 0\leq\beta\leq1, \qquad 0\leq\gamma\leq1 \]

Parametr \(\alpha\)

Parametr \(\alpha\) steruje aktualizacją poziomu.

Parametr \(\beta\)

Parametr \(\beta\) określa szybkość aktualizacji trendu.

Parametr \(\gamma\)

Parametr \(\gamma\) steruje aktualizacją składników sezonowych.

Graniczne wartości parametrów i związki z prostszymi metodami

Parametry wygładzania określają, w jakim stopniu model reaguje na nowe informacje. Dla niektórych wartości granicznych metody wygładzania wykładniczego stają się równoważne prostszym metodom prognozowania albo przestają aktualizować wybrane składniki modelu.

Metoda Browna dla \(\alpha=1\)

W metodzie Browna wartość wygładzona wynosi:

\[ F_t = \alpha y_t + (1-\alpha)F_{t-1} \]

Jeżeli przyjmiemy:

\[ \alpha=1 \]

otrzymujemy:

\[ F_t=y_t \]

Ponieważ prognoza na kolejny okres jest równa:

\[ \widehat{y}_{t+1\mid t}=F_t \]

to dla \(\alpha=1\):

\[ \widehat{y}_{t+1\mid t}=y_t \]

Metoda Browna staje się wtedy dokładnie metodą naiwną poziomu bez zmian. Model całkowicie pomija wcześniejszą wartość wygładzoną i jako prognozę przyjmuje ostatnią obserwację.

Metoda Browna dla \(\alpha=0\)

Jeżeli:

\[ \alpha=0 \]

to:

\[ F_t=F_{t-1} \]

Wartość wygładzona nie reaguje wówczas na nowe obserwacje. Wszystkie kolejne prognozy pozostają równe przyjętej wartości początkowej:

\[ \widehat{y}_{t+1\mid t}=F_1 \]

Parametr równy zero oznacza więc całkowite zamrożenie poziomu modelu.

Metoda Holta a metoda naiwna przyrostu bez zmian

W metodzie Holta:

\[ F_t = \alpha y_t + (1-\alpha)(F_{t-1}+S_{t-1}) \]

\[ S_t = \beta(F_t-F_{t-1}) + (1-\beta)S_{t-1} \]

Jeżeli przyjmiemy:

\[ \alpha=1, \qquad \beta=1 \]

to poziom jest równy bieżącej obserwacji:

\[ F_t=y_t \]

a trend staje się ostatnim przyrostem:

\[ S_t = F_t-F_{t-1} = y_t-y_{t-1} \]

Prognoza wynosi wtedy:

\[ \widehat{y}_{t+1\mid t} = F_t+S_t = y_t+(y_t-y_{t-1}) \]

Otrzymujemy zatem metodę naiwną przyrostu bez zmian. Ostatni zaobserwowany przyrost jest przenoszony na kolejny okres.

Metoda Holta dla \(\beta=0\)

Jeżeli:

\[ \beta=0 \]

to równanie trendu przyjmuje postać:

\[ S_t=S_{t-1} \]

Trend nie jest aktualizowany i pozostaje równy wartości początkowej. Model prognozuje wówczas z wykorzystaniem stałego przyrostu ustalonego podczas inicjalizacji.

Od metody Wintersa do metody Holta

Metoda Wintersa przechodzi w metodę Holta, jeżeli składnik sezonowy jest neutralny.

W wariancie addytywnym neutralna sezonowość oznacza:

\[ C_t=0 \]

W wariancie multiplikatywnym:

\[ C_t=1 \]

Po usunięciu wpływu sezonowości prognoza opiera się wyłącznie na poziomie i trendzie:

\[ \widehat{y}_{t+h\mid t} = F_t+hS_t \]

czyli przyjmuje postać prognozy metody Holta.

Znaczenie granicznych wartości \(\gamma\)

W metodzie Wintersa parametr \(\gamma\) steruje aktualizacją sezonowości.

Metody Browna, Holta i Wintersa nie są całkowicie odrębnymi konstrukcjami. Dla odpowiednich wartości parametrów i składników modelu przechodzą w prostsze metody, w tym w metody naiwne.

Infografika wyjaśniająca znaczenie parametrów alfa, beta i gamma w wygładzaniu wykładniczym
Parametr alfa odpowiada za poziom, beta za trend, a gamma za sezonowość. Ich dobór wpływa na szybkość reakcji modelu na nowe obserwacje.

Dobór parametrów za pomocą Solvera

Parametry można przyjąć arbitralnie, sprawdzić kilka wariantów albo wyznaczyć za pomocą optymalizacji numerycznej. W arkuszach kalkulacyjnych często wykorzystuje się do tego Solver.

Dla każdego okresu, dla którego znamy prognozę jednookresową i wartość empiryczną, obliczamy błąd:

\[ e_t = y_t-\widehat{y}_{t\mid t-1} \]

Następnie obliczamy sumę kwadratów błędów:

\[ SSE = \sum_{t=t_0}^{T}e_t^2 \]

Solver powinien minimalizować \(SSE\), zmieniając odpowiednie parametry:

Do komórek zmienianych dodajemy ograniczenia:

\[ 0\leq\alpha,\beta,\gamma\leq1 \]

Dlaczego wystarczy minimalizować SSE?

Średni błąd kwadratowy ma postać:

\[ MSE = \frac{SSE}{n} \]

Pierwiastek średniego błędu kwadratowego:

\[ RMSE = \sqrt{\frac{SSE}{n}} \]

Liczba \(n\) jest podczas optymalizacji stała, a pierwiastek kwadratowy jest funkcją rosnącą. Dlatego:

\[ \arg\min SSE = \arg\min MSE = \arg\min RMSE \]

Minimalizacja samej sumy kwadratów jest prostsza obliczeniowo i prowadzi do tych samych optymalnych parametrów.

W Solverze najlepiej minimalizować bezpośrednio sumę kwadratów błędów. Dzielenie przez liczbę obserwacji i wyciąganie pierwiastka nie zmienia położenia minimum.

Próba ucząca, testowa i walidacja krocząca

Minimalizacja błędu na całym dostępnym szeregu pokazuje przede wszystkim jakość dopasowania wewnątrz próby. Model może bardzo dobrze odtwarzać dane historyczne, ale gorzej prognozować nowe obserwacje.

Dlatego warto podzielić szereg na dwie części.

Próba ucząca

Na próbie uczącej:

Próba testowa

Na próbie testowej obliczamy prognozy bez ponownego dopasowywania parametrów. Dzięki temu sprawdzamy, jak model radzi sobie z danymi, których nie wykorzystano podczas strojenia.

Walidacja krocząca

W walidacji kroczącej prognozujemy kolejno każdy następny okres, korzystając wyłącznie z informacji dostępnych przed jego rozpoczęciem.

Procedura może wyglądać następująco:

  1. wyznacz prognozę na okres \(t+1\) na podstawie danych do okresu \(t\),
  2. porównaj prognozę z rzeczywistą wartością \(y_{t+1}\),
  3. dołącz nową obserwację do zbioru informacji,
  4. wyznacz prognozę na kolejny okres.

Wszystkie porównywane modele powinny być oceniane na tym samym zakresie, przy tym samym horyzoncie oraz za pomocą tej samej miary błędu.

Przykład metody Browna

Załóżmy, że wartości szeregu wynoszą:

\[ 100,\quad106,\quad103,\quad108,\quad110 \]

Przyjmujemy:

\[ \alpha=0{,}4, \qquad F_1=y_1=100 \]

Dla okresu drugiego:

\[ F_2 = 0{,}4\cdot106 + 0{,}6\cdot100 = 102{,}4 \]

Dla okresu trzeciego:

\[ F_3 = 0{,}4\cdot103 + 0{,}6\cdot102{,}4 = 102{,}64 \]

Dla okresu czwartego:

\[ F_4 = 0{,}4\cdot108 + 0{,}6\cdot102{,}64 = 104{,}784 \]

Dla okresu piątego:

\[ F_5 = 0{,}4\cdot110 + 0{,}6\cdot104{,}784 = 106{,}8704 \]

Prognoza na okres szósty:

\[ \widehat{y}_{6\mid5} = F_5 = 106{,}8704 \]

Przykład metody Holta

Załóżmy, że obserwacje wynoszą:

\[ 100,\quad104,\quad109,\quad113 \]

Przyjmujemy:

\[ \alpha=0{,}5, \qquad \beta=0{,}4 \]

oraz wartości początkowe:

\[ F_1=100, \qquad S_1=104-100=4 \]

Dla okresu drugiego:

\[ F_2 = 0{,}5\cdot104 + 0{,}5(100+4) = 104 \]

\[ S_2 = 0{,}4(104-100) + 0{,}6\cdot4 = 4 \]

Prognoza na okres trzeci:

\[ \widehat{y}_{3\mid2} = F_2+S_2 = 108 \]

Po poznaniu wartości \(y_3=109\) aktualizujemy model:

\[ F_3 = 0{,}5\cdot109 + 0{,}5(104+4) = 108{,}5 \]

\[ S_3 = 0{,}4(108{,}5-104) + 0{,}6\cdot4 = 4{,}2 \]

Prognoza na okres czwarty:

\[ \widehat{y}_{4\mid3} = 108{,}5+4{,}2 = 112{,}7 \]

Prognoza na dwa okresy naprzód, wyznaczona po okresie trzecim:

\[ \widehat{y}_{5\mid3} = F_3+2S_3 = 108{,}5+2\cdot4{,}2 = 116{,}9 \]

Zalety i ograniczenia

Zalety

Ograniczenia

Podsumowanie

Wygładzanie wykładnicze jest grupą adaptacyjnych metod prognozowania, w których wpływ starszych obserwacji maleje wykładniczo.

Metoda Browna wykorzystuje jeden szereg wartości wygładzonych:

\[ F_t = \alpha y_t + (1-\alpha)F_{t-1} \]

a prognoza na kolejny okres wynosi:

\[ \widehat{y}_{t+1\mid t}=F_t \]

Metoda Holta uwzględnia poziom i trend:

\[ \widehat{y}_{t+h\mid t} = F_t+hS_t \]

Metoda Wintersa dodaje składnik sezonowy \(C_t\) i występuje w wariancie addytywnym oraz multiplikatywnym.

Parametry \(\alpha,\beta,\gamma\) można wyznaczać przez minimalizację sumy kwadratów błędów prognoz jednookresowych:

\[ SSE = \sum e_t^2 \]

Rzetelna ocena modelu powinna obejmować nie tylko dopasowanie wewnątrz próby, lecz również próbę testową albo walidację kroczącą.

Ze względu na różnorodność konwencji przy każdej analizie należy jasno określić sposób inicjalizacji, oznaczenia szeregów pomocniczych, funkcję celu oraz procedurę wyznaczania prognoz.

Rachunek wektorowy w statyce – siły, składowe i momenty

Rachunek wektorowy jest jednym z podstawowych narzędzi statyki. Pozwala zapisywać siły, momenty i położenia punktów za pomocą składowych, dodawać obciążenia oraz obliczać moment siły jako iloczyn wektorowy.

W artykule Czym zajmuje się statyka? Podstawowe pojęcia i zasady pokazaliśmy, że równowaga ciała zależy od sumy działających na nie sił oraz momentów. Teraz rozwiniemy te zagadnienia od strony rachunkowej, ze szczególnym uwzględnieniem zapisu za pomocą wersorów \(\vec{i}\), \(\vec{j}\) i \(\vec{k}\).

Dlaczego w statyce używa się wektorów?

Nie każdą wielkość fizyczną można opisać jedną liczbą. Masa, temperatura i czas są wielkościami skalarnymi, ponieważ do ich określenia wystarcza wartość liczbowa wraz z jednostką. Inaczej jest w przypadku siły.

Siła jest określona nie tylko przez swoją wartość, ale również przez:

Z tego powodu siła jest wielkością wektorową. Wielkością wektorową jest także moment siły, którego kierunek i zwrot wynikają z położenia siły oraz reguły prawej dłoni.

Rachunek wektorowy umożliwia między innymi:

Wektor w układzie współrzędnych

W statyce najczęściej posługujemy się prostokątnym układem współrzędnych. W przestrzeni wykorzystuje się trzy osie: \(x\), \(y\) i \(z\). Z każdą z nich związany jest odpowiedni wektor jednostkowy, nazywany również wersorem:

Każdy z wersorów ma wartość równą \(1\). W prostokątnym układzie współrzędnych wersory są wzajemnie prostopadłe.

Dowolny wektor siły można zapisać jako sumę składowych:

\[ \vec{F}=F_x\vec{i}+F_y\vec{j}+F_z\vec{k}. \]

Równoważny jest zapis kolumnowy:

\[ \vec{F}= \begin{bmatrix} F_x\\ F_y\\ F_z \end{bmatrix}. \]

Na przykład wektor:

\[ \vec{F}=3\vec{i}-2\vec{j}+4\vec{k}\ \mathrm{N} \]

można również zapisać jako:

\[ \vec{F}= \begin{bmatrix} 3\\ -2\\ 4 \end{bmatrix} \mathrm{N}. \]

W płaskim układzie sił wszystkie wektory leżą w jednej płaszczyźnie, na przykład \(xy\). Ich składowa wzdłuż osi \(z\) jest wtedy równa zeru:

\[ \vec{F}=F_x\vec{i}+F_y\vec{j}. \]
Wektor siły w przestrzennym układzie współrzędnych oraz jego składowe Fx, Fy i Fz
Wektor w przestrzeni można zapisać jako sumę składowych skierowanych wzdłuż osi układu współrzędnych.

Wartość wektora

Jeżeli znamy składowe wektora, jego wartość obliczamy ze wzoru:

\[ |\vec{F}|=\sqrt{F_x^2+F_y^2+F_z^2}. \]

W przypadku płaskim:

\[ |\vec{F}|=\sqrt{F_x^2+F_y^2}. \]

Rozważmy wektor:

\[ \vec{F}=3\vec{i}-4\vec{j}\ \mathrm{N}. \]

Jego wartość wynosi:

\[ |\vec{F}| = \sqrt{3^2+(-4)^2} = 5\ \mathrm{N}. \]

Ważne rozróżnienie: zapis \(\vec{F}\) oznacza cały wektor, natomiast \(F=|\vec{F}|\) oznacza wyłącznie jego wartość.

Wektor wyznaczony przez dwa punkty

W statyce często wyznacza się wektor łączący dwa punkty. Jest on potrzebny między innymi do określania kierunku pręta lub cięgna oraz do obliczania momentu siły.

Jeżeli dane są punkty:

\[ A(x_A,y_A,z_A), \qquad B(x_B,y_B,z_B), \]

to wektor skierowany od punktu \(A\) do punktu \(B\) ma postać:

\[ \overrightarrow{AB} = (x_B-x_A)\vec{i} + (y_B-y_A)\vec{j} + (z_B-z_A)\vec{k}. \]

Dla punktów:

\[ A(1,2,0), \qquad B(4,3,2) \]

otrzymujemy:

\[ \overrightarrow{AB} = (4-1)\vec{i} + (3-2)\vec{j} + (2-0)\vec{k} = 3\vec{i}+\vec{j}+2\vec{k}. \]

Kolejność punktów ma znaczenie:

\[ \overrightarrow{BA}=-\overrightarrow{AB}. \]

Wersor wybranego kierunku

Jeżeli znamy kierunek działania siły, możemy opisać go za pomocą wersora. Dla niezerowego wektora kierunkowego \(\vec{r}\) odpowiedni wektor jednostkowy wynosi:

\[ \vec{e}_r=\frac{\vec{r}}{|\vec{r}|}. \]

Siłę o wartości \(F\), działającą wzdłuż tego kierunku, można zapisać jako:

\[ \vec{F}=F\vec{e}_r. \]

Zapis ten jest szczególnie użyteczny w zadaniach dotyczących lin, cięgien i prętów dwusiłowych, ponieważ występująca w nich siła działa wzdłuż osi elementu.

Dodawanie wektorów

Wektory można dodawać zarówno geometrycznie, jak i algebraicznie.

Metoda trójkąta

Początek drugiego wektora umieszcza się w końcu pierwszego. Wektor poprowadzony od początku pierwszego do końca drugiego jest ich sumą:

\[ \vec{R}=\vec{F}_1+\vec{F}_2. \]

Metoda równoległoboku

Oba wektory zaczepia się we wspólnym punkcie, a następnie buduje równoległobok. Jego przekątna wychodząca ze wspólnego początku jest sumą wektorów.

Dodawanie składowych

Jeżeli:

\[ \vec{F}_1 = F_{1x}\vec{i} + F_{1y}\vec{j} + F_{1z}\vec{k}, \] \[ \vec{F}_2 = F_{2x}\vec{i} + F_{2y}\vec{j} + F_{2z}\vec{k}, \]

to ich suma wynosi:

\[ \vec{R} = (F_{1x}+F_{2x})\vec{i} + (F_{1y}+F_{2y})\vec{j} + (F_{1z}+F_{2z})\vec{k}. \]

Rozważmy siły:

\[ \vec{F}_1=3\vec{i}+2\vec{j}\ \mathrm{N}, \qquad \vec{F}_2=-\vec{i}+4\vec{j}\ \mathrm{N}. \]

Ich suma jest równa:

\[ \vec{R} = (3-1)\vec{i} + (2+4)\vec{j} = 2\vec{i}+6\vec{j}\ \mathrm{N}. \]

Wartość wypadkowej wynosi:

\[ R = \sqrt{2^2+6^2} = 2\sqrt{10}\ \mathrm{N}. \]

Odejmowanie wektorów

Odejmowanie wektorów sprowadza się do dodania wektora przeciwnego:

\[ \vec{F}_1-\vec{F}_2 = \vec{F}_1+(-\vec{F}_2). \]

Wektor przeciwny ma tę samą wartość i kierunek, lecz przeciwny zwrot:

\[ -\vec{F}_2 = -F_{2x}\vec{i} -F_{2y}\vec{j} -F_{2z}\vec{k}. \]

Na przykład dla:

\[ \vec{F}_1=5\vec{i}+3\vec{j}, \qquad \vec{F}_2=2\vec{i}-\vec{j} \]

otrzymujemy:

\[ \vec{F}_1-\vec{F}_2 = (5-2)\vec{i} + \bigl(3-(-1)\bigr)\vec{j} = 3\vec{i}+4\vec{j}. \]
Dodawanie wektorów metodą trójkąta i równoległoboku oraz odejmowanie przez dodanie wektora przeciwnego
Wektory można dodawać metodą trójkąta lub równoległoboku, a odejmowanie sprowadza się do dodania wektora przeciwnego.

Rozkład siły na składowe

W statyce często trzeba przejść od wartości i kierunku siły do jej składowych w przyjętym układzie współrzędnych.

Jeżeli siła \(F\) działająca w płaszczyźnie tworzy z dodatnim kierunkiem osi \(x\) kąt \(\alpha\), to:

\[ F_x=F\cos\alpha, \qquad F_y=F\sin\alpha. \]

Znaki składowych zależą jednak od rzeczywistego zwrotu siły. Dla siły skierowanej w lewo i do góry może być na przykład:

\[ F_x=-F\cos\alpha, \qquad F_y=F\sin\alpha. \]

W przestrzeni można posłużyć się cosinusami kierunkowymi:

\[ \vec{F} = F\cos\alpha\,\vec{i} + F\cos\beta\,\vec{j} + F\cos\gamma\,\vec{k}, \]

gdzie \(\alpha\), \(\beta\) i \(\gamma\) są kątami między wektorem siły a osiami \(x\), \(y\) i \(z\). Zachodzi przy tym zależność:

\[ \cos^2\alpha+\cos^2\beta+\cos^2\gamma=1. \]

Iloczyn skalarny

Iloczyn skalarny dwóch wektorów \(\vec{a}\) i \(\vec{b}\) jest liczbą zdefiniowaną wzorem:

\[ \vec{a}\cdot\vec{b} = |\vec{a}|\,|\vec{b}|\cos\theta, \]

gdzie \(\theta\) oznacza kąt między wektorami.

W zapisie składowym:

\[ \vec{a}\cdot\vec{b} = a_xb_x+a_yb_y+a_zb_z. \]

Iloczyn skalarny można wykorzystywać do:

Jeżeli \(\vec{e}_u\) jest wersorem kierunku \(u\), rzut skalarny siły \(\vec{F}\) na ten kierunek wynosi:

\[ F_u=\vec{F}\cdot\vec{e}_u. \]

Po uwzględnieniu kąta między wektorami:

\[ F_u=|\vec{F}|\cos\theta. \]

Rzut może być dodatni, ujemny albo równy zeru. Wynik zerowy oznacza, że wektor siły jest prostopadły do danego kierunku.

Rzut wektora siły na wybrany kierunek obliczany za pomocą iloczynu skalarnego
Iloczyn skalarny wektora i wersora danego kierunku wyznacza rzut skalarny wektora na ten kierunek.

Iloczyn wektorowy

W statyce szczególne znaczenie ma iloczyn wektorowy, ponieważ właśnie za jego pomocą oblicza się moment siły.

Iloczyn wektorowy:

\[ \vec{a}\times\vec{b} \]

jest wektorem:

Wartość iloczynu wektorowego wynosi:

\[ |\vec{a}\times\vec{b}| = |\vec{a}|\,|\vec{b}|\sin\theta. \]

Iloczyny wersorów bazowych

W prawoskrętnym, prostokątnym układzie współrzędnych zachodzą zależności:

\[ \vec{i}\times\vec{j}=\vec{k}, \qquad \vec{j}\times\vec{k}=\vec{i}, \qquad \vec{k}\times\vec{i}=\vec{j}. \]

Po zamianie kolejności wektorów zmienia się znak wyniku:

\[ \vec{j}\times\vec{i}=-\vec{k}, \qquad \vec{k}\times\vec{j}=-\vec{i}, \qquad \vec{i}\times\vec{k}=-\vec{j}. \]

Iloczyn wektora przez wektor równoległy do niego jest równy zeru. W szczególności:

\[ \vec{i}\times\vec{i} = \vec{j}\times\vec{j} = \vec{k}\times\vec{k} = \vec{0}. \]

Iloczyn wektorowy nie jest przemienny:

\[ \vec{a}\times\vec{b} = -\vec{b}\times\vec{a}. \]

Iloczyn wektorowy w zapisie wyznacznikowym

Iloczyn wektorowy można obliczać za pomocą formalnego zapisu przypominającego wyznacznik:

\[ \vec{a}\times\vec{b} = \begin{vmatrix} \vec{i} & \vec{j} & \vec{k}\\ a_x & a_y & a_z\\ b_x & b_y & b_z \end{vmatrix}. \]

Po rozwinięciu otrzymujemy:

\[ \vec{a}\times\vec{b} = (a_yb_z-a_zb_y)\vec{i} + (a_zb_x-a_xb_z)\vec{j} + (a_xb_y-a_yb_x)\vec{k}. \]

Ten sposób obliczeń jest szczególnie wygodny przy wyznaczaniu momentów sił w przestrzeni. Zasady obliczania i rozwijania wyznaczników zostały szerzej omówione w artykule Wyznacznik macierzy kwadratowej.

Obliczanie iloczynu wektorowego za pomocą wyznacznika i zapis jego składowych
Iloczyn wektorowy można obliczyć przez rozwinięcie formalnego wyznacznika zawierającego wersory i składowe obu wektorów.

Moment siły jako iloczyn wektorowy

Najważniejszym zastosowaniem iloczynu wektorowego w statyce jest obliczanie momentu siły.

Moment siły \(\vec{F}\) względem punktu \(O\) jest równy:

\[ \vec{M}_O=\vec{r}\times\vec{F}, \]

gdzie:

Wartość momentu wynosi:

\[ |\vec{M}_O| = |\vec{r}|\,|\vec{F}|\sin\theta, \]

gdzie \(\theta\) jest kątem między wektorem położenia i wektorem siły.

Ponieważ:

\[ d=|\vec{r}|\sin\theta, \]

gdzie \(d\) jest prostopadłą odległością punktu \(O\) od linii działania siły, otrzymujemy również znany wzór:

\[ M_O=Fd. \]

Interpretacja przestrzenna

Jeżeli wektor \(\vec{r}\) i siła \(\vec{F}\) leżą w płaszczyźnie \(xy\), to ich iloczyn wektorowy jest do niej prostopadły. Wektor momentu jest więc skierowany wzdłuż osi \(z\) albo przeciwnie do niej.

To właśnie dlatego w zadaniach płaskich moment zapisuje się zwykle jako dodatni albo ujemny. Znak określa zwrot wektora momentu wzdłuż osi prostopadłej do płaszczyzny rysunku.

Jeżeli patrzymy od dodatniego końca osi \(z\), moment skierowany w stronę \(+z\) odpowiada obrotowi przeciwnemu do ruchu wskazówek zegara, natomiast moment skierowany w stronę \(-z\) odpowiada obrotowi zgodnemu z ruchem wskazówek zegara.

Wektor położenia i siła w płaszczyźnie xy oraz prostopadły do niej wektor momentu
Moment siły jest iloczynem wektorowym wektora położenia i siły, dlatego jest prostopadły do wyznaczonej przez nie płaszczyzny.

Przykład obliczenia momentu siły

Załóżmy, że wektor położenia punktu przyłożenia siły oraz wektor siły wynoszą:

\[ \vec{r}=2\vec{i}+3\vec{j}\ \mathrm{m}, \qquad \vec{F}=4\vec{i}-2\vec{j}\ \mathrm{N}. \]

Dopisujemy zerowe składowe wzdłuż osi \(z\):

\[ \vec{r}=2\vec{i}+3\vec{j}+0\vec{k}, \] \[ \vec{F}=4\vec{i}-2\vec{j}+0\vec{k}. \]

Moment siły obliczamy jako iloczyn wektorowy:

\[ \vec{M}_O = \begin{vmatrix} \vec{i} & \vec{j} & \vec{k}\\ 2 & 3 & 0\\ 4 & -2 & 0 \end{vmatrix}. \]

Składowe wzdłuż osi \(x\) i \(y\) są równe zeru. Pozostaje:

\[ \vec{M}_O = \bigl(2\cdot(-2)-3\cdot4\bigr)\vec{k}. \]
\[ \vec{M}_O = (-4-12)\vec{k} = -16\vec{k}\ \mathrm{N\,m}. \]

Wartość momentu wynosi zatem:

\[ |\vec{M}_O|=16\ \mathrm{N\,m}. \]

Znak ujemny oznacza, że wektor momentu jest skierowany przeciwnie do dodatniego zwrotu osi \(z\). Przy standardowej orientacji osi odpowiada to obrotowi zgodnemu z ruchem wskazówek zegara, gdy patrzymy od dodatniego końca osi \(z\).

Moment względem osi

W przestrzennych zagadnieniach statyki może nas interesować nie cały wektor momentu, lecz jego składowa względem wybranej osi.

Jeżeli \(\vec{e}_u\) jest wersorem osi \(u\), to moment względem tej osi wynosi:

\[ M_u=\vec{e}_u\cdot\vec{M}_O. \]

Po uwzględnieniu definicji momentu:

\[ M_u = \vec{e}_u\cdot \left(\vec{r}\times\vec{F}\right). \]

Jest to iloczyn mieszany trzech wektorów. Jego wynik jest liczbą oznaczającą rzut wektora momentu na rozpatrywaną oś.

Sumowanie wielu sił i momentów

Jeżeli na ciało działa wiele sił, główny wektor układu otrzymujemy przez ich dodanie:

\[ \vec{R}=\sum_i\vec{F}_i. \]

Główny moment układu względem punktu \(O\) wynosi:

\[ \vec{M}_O = \sum_i\vec{r}_i\times\vec{F}_i. \]

Jeżeli na ciało działają także przyłożone momenty skupione albo pary sił, należy dodać je bezpośrednio:

\[ \vec{M}_O = \sum_i\vec{r}_i\times\vec{F}_i + \sum_j\vec{M}_j. \]

W ten sposób dowolny układ sił redukuje się w wybranym punkcie do głównego wektora sił i głównego momentu.

Równania równowagi w zapisie wektorowym

Warunki równowagi ciała sztywnego można zapisać w zwartej postaci wektorowej:

\[ \sum_i\vec{F}_i=\vec{0}, \] \[ \sum_i\vec{r}_i\times\vec{F}_i + \sum_j\vec{M}_j = \vec{0}. \]

Pierwsze równanie oznacza, że główny wektor sił jest równy zeru. Drugie wymaga, aby zerowy był również główny moment układu.

Po rozpisaniu tych równań na składowe w przestrzeni otrzymujemy sześć równań skalarnych:

\[ \sum F_x=0, \qquad \sum F_y=0, \qquad \sum F_z=0, \] \[ \sum M_x=0, \qquad \sum M_y=0, \qquad \sum M_z=0. \]

W układzie płaskim pozostają trzy niezależne warunki równowagi: dwa dotyczące sił oraz jeden dotyczący momentów.

Moment siły w anaglifie 3D

Wzajemne położenie wektora \(\vec{r}\), siły \(\vec{F}\) i momentu \(\vec{M}_O\) ma charakter przestrzenny. Można je obejrzeć również na ilustracji anaglifowej przeznaczonej do czerwono-cyjanowych okularów 3D.

Czerwony filtr należy umieścić przed lewym okiem, a cyjanowy przed prawym. Najlepszy efekt uzyskuje się, patrząc na obraz prostopadle i dobierając taką odległość od ekranu, przy której linie obu kolorów łączą się w jeden przestrzenny rysunek.

Anaglif czerwono-cyjanowy pokazujący wektor położenia, siłę i prostopadły wektor momentu
Anaglif umożliwia przestrzenne obejrzenie wzajemnego położenia wektora \(\vec{r}\), siły \(\vec{F}\) i momentu \(\vec{M}_O\). Kliknięcie otwiera obraz w pełnym rozmiarze.

Najczęstsze błędy

Kontrola wyniku: jednostką siły jest niuton, natomiast jednostką momentu niutonometr. Jeżeli wynik iloczynu \(\vec{r}\times\vec{F}\) ma jednostkę niutona zamiast niutonometra, oznacza to, że w obliczeniach pominięto jednostkę długości.

Podsumowanie

Rachunek wektorowy jest podstawowym narzędziem statyki, ponieważ pozwala w jednoznaczny sposób opisywać siły, kierunki ich działania, momenty oraz warunki równowagi.

Najważniejsze zależności omówione w artykule to:

\[ \vec{F} = F_x\vec{i}+F_y\vec{j}+F_z\vec{k}, \] \[ |\vec{F}| = \sqrt{F_x^2+F_y^2+F_z^2}, \] \[ \vec{a}\cdot\vec{b} = a_xb_x+a_yb_y+a_zb_z, \] \[ \vec{a}\times\vec{b} = (a_yb_z-a_zb_y)\vec{i} + (a_zb_x-a_xb_z)\vec{j} + (a_xb_y-a_yb_x)\vec{k}, \] \[ \vec{M}_O=\vec{r}\times\vec{F}, \] \[ \sum_i\vec{F}_i=\vec{0}, \qquad \sum_i\vec{r}_i\times\vec{F}_i+\sum_j\vec{M}_j=\vec{0}. \]

Opanowanie tych działań znacznie ułatwia redukcję układów sił, obliczanie reakcji podporowych oraz analizowanie płaskich i przestrzennych układów obciążeń.

Czym zajmuje się statyka? Podstawowe pojęcia i zasady

Statyka jest działem mechaniki zajmującym się równowagą ciał oraz układów materialnych poddanych działaniu sił. Pozwala ustalić, czy ciało pozostanie w spoczynku, a także wyznaczyć siły i momenty występujące w podporach, połączeniach oraz elementach konstrukcji.

Z metod statyki korzysta się podczas analizowania belek, ram, kratownic, maszyn, mostów i budynków. Zanim jednak przejdziemy do rozwiązywania konkretnych zadań, trzeba poznać podstawowe pojęcia: siłę, moment siły, układ sił, reakcje więzów oraz warunki równowagi.

Czym jest statyka?

Mechanika klasyczna opisuje ruch ciał oraz oddziaływania, które mogą ten ruch powodować. Najczęściej dzieli się ją na trzy podstawowe części:

W typowych zadaniach statyki rozpatruje się ciała pozostające w spoczynku względem przyjętego układu odniesienia. Nie oznacza to, że nie działają na nie żadne siły. Sił może być wiele, ale ich łączny skutek nie może powodować ani przesuwania, ani obracania ciała.

Statyka nie zajmuje się jeszcze bezpośrednio wytrzymałością elementu. Najpierw pozwala wyznaczyć działające siły, reakcje i momenty. Dopiero później, na przykład w wytrzymałości materiałów, na ich podstawie oblicza się naprężenia, odkształcenia i wymagane wymiary konstrukcji.

Co oznacza równowaga ciała?

Ciało znajduje się w równowadze, jeżeli nie zmienia swojego stanu ruchu. W statyce interesuje nas przede wszystkim przypadek, w którym ciało pozostaje w spoczynku.

W przypadku ciała sztywnego trzeba rozpatrzyć dwa możliwe skutki działania obciążeń:

  1. przesunięcie ciała;
  2. obrót ciała.

Brak przesunięcia wymaga, aby suma wszystkich sił była równa zeru:

\[ \sum \vec{F}=\vec{0}. \]

Brak obrotu wymaga natomiast, aby suma momentów sił względem dowolnie wybranego punktu \(O\) była równa zeru:

\[ \sum \vec{M}_O=\vec{0}. \]

Dopiero jednoczesne spełnienie obu warunków oznacza pełną równowagę ciała sztywnego. Zerowa suma sił nie wystarcza, ponieważ układ może mieć zerową wypadkową, a mimo to powodować obrót.

Punkt materialny i ciało sztywne

Rzeczywiste obiekty są zwykle zbyt złożone, aby uwzględniać wszystkie ich cechy. W mechanice zastępuje się je więc uproszczonymi modelami.

Punkt materialny

Punkt materialny jest modelem ciała, którego wymiary można w danym zagadnieniu pominąć. Całą jego masę traktuje się jako skupioną w jednym punkcie.

Punkt materialny może się przemieszczać, ale nie analizujemy jego obrotu wokół własnej osi. Warunek jego równowagi sprowadza się zatem do zerowania sumy działających sił:

\[ \sum \vec{F}=\vec{0}. \]

Ciało sztywne

Ciało sztywne jest modelem obiektu, którego odkształcenia pomijamy. Zakładamy, że odległości między dowolnymi punktami ciała pozostają stałe.

Ciało sztywne może zarówno przesuwać się, jak i obracać. Z tego powodu przy jego analizie trzeba uwzględniać nie tylko siły, lecz również ich momenty. Model ten stosuje się między innymi do belek, ram, dźwigni, tarcz oraz elementów maszyn.

Zasady Newtona jako podstawa statyki

Podstawy statyki wynikają z trzech zasad dynamiki Newtona. Pierwsza zasada wyjaśnia, czym jest stan równowagi, druga pozwala zapisać jego matematyczne warunki, natomiast trzecia opisuje wzajemne oddziaływania między ciałami.

Pierwsza zasada Newtona

Pierwsza zasada Newtona, nazywana także zasadą bezwładności, mówi, że jeżeli na ciało nie działa żadna siła wypadkowa, pozostaje ono w spoczynku albo porusza się ruchem jednostajnym prostoliniowym.

\[ \sum \vec{F}=\vec{0} \quad\Longrightarrow\quad \vec{v}=\mathrm{const}. \]

Brak siły wypadkowej nie musi więc zawsze oznaczać spoczynku. Może również oznaczać ruch ze stałą prędkością. W typowych zadaniach statyki rozpatruje się jednak ciało pozostające w spoczynku względem inercjalnego układu odniesienia.

Druga zasada Newtona

Druga zasada Newtona wiąże wypadkową siłę działającą na ciało z jego przyspieszeniem:

\[ \sum \vec{F}=m\vec{a}. \]

Jeżeli ciało pozostaje w równowadze, jego przyspieszenie jest równe zeru, a więc:

\[ \vec{a}=\vec{0} \quad\Longrightarrow\quad \sum \vec{F}=\vec{0}. \]

Dla ciała sztywnego trzeba dodatkowo uwzględnić możliwość obrotu. Brak przyspieszenia kątowego prowadzi do warunku zerowania sumy momentów:

\[ \sum \vec{M}_O=\vec{0}. \]

Równania statyki można więc traktować jako szczególny przypadek równań dynamiki, w którym przyspieszenie liniowe i kątowe są równe zeru.

Trzecia zasada Newtona

Trzecia zasada Newtona, nazywana zasadą akcji i reakcji, mówi, że jeżeli ciało \(A\) działa na ciało \(B\) pewną siłą, to ciało \(B\) działa na ciało \(A\) siłą o tej samej wartości i przeciwnym zwrocie:

\[ \vec{F}_{AB}=-\vec{F}_{BA}. \]

Siły akcji i reakcji mają wspólny kierunek oraz jednakowe wartości, ale działają na różne ciała. Z tego powodu nie równoważą się na schemacie jednego, odizolowanego elementu.

Trzecia zasada Newtona ma podstawowe znaczenie podczas wyznaczania reakcji podpór, sił w przegubach, nacisków, sił tarcia oraz sił występujących w prętach i cięgnach.

Infografika przedstawiająca pierwszą, drugą i trzecią zasadę Newtona
Pierwsza zasada Newtona opisuje bezwładność, druga wiąże siłę z przyspieszeniem, a trzecia wyjaśnia wzajemne oddziaływania ciał.

Siła jako wielkość wektorowa

Siła jest miarą oddziaływania jednego ciała na drugie. Może zmieniać stan ruchu ciała albo powodować jego odkształcenie.

Siła jest wielkością wektorową. Do jej pełnego określenia potrzebne są:

Jednostką siły w układzie SI jest niuton:

\[ 1\ \mathrm{N} = 1\ \mathrm{kg}\cdot\frac{\mathrm{m}}{\mathrm{s}^2}. \]

Graficznie siłę przedstawia się za pomocą skierowanego odcinka. Długość strzałki odpowiada wartości siły, natomiast jej ustawienie określa kierunek i zwrot.

Linia działania siły

Prosta przechodząca przez punkt przyłożenia siły i zgodna z jej kierunkiem nazywa się linią działania siły.

W przypadku ciała sztywnego siłę można przesuwać wzdłuż jej linii działania bez zmiany jej zewnętrznego skutku mechanicznego. Nie można natomiast przenieść jej na inną, równoległą prostą bez dodania odpowiedniego momentu.

Rozkład siły na składowe

Siłę można rozłożyć na składowe zgodne z osiami przyjętego układu współrzędnych. W układzie płaskim:

\[ \vec{F}=F_x\vec{i}+F_y\vec{j}. \]

Jeżeli siła \(F\) tworzy z dodatnim kierunkiem osi \(x\) kąt \(\alpha\), jej składowe wynoszą:

\[ F_x=F\cos\alpha, \qquad F_y=F\sin\alpha. \]

Znaki składowych zależą od rzeczywistego zwrotu siły oraz od przyjętych dodatnich zwrotów osi.

Moment siły

Siła może powodować nie tylko przesunięcie, ale również obrót ciała. Miarą zdolności siły do wywoływania obrotu względem wybranego punktu jest moment siły.

W zapisie wektorowym moment siły \(\vec{F}\) względem punktu \(O\) ma postać:

\[ \vec{M}_O=\vec{r}\times\vec{F}, \]

gdzie \(\vec{r}\) jest wektorem poprowadzonym od punktu \(O\) do dowolnego punktu leżącego na linii działania siły.

W układzie płaskim wartość momentu można obliczyć ze wzoru:

\[ M_O=F d, \] 

gdzie \(d\) jest prostopadłą odległością punktu \(O\) od linii działania siły, nazywaną ramieniem siły.

Jednostką momentu siły jest niutonometr:

\[ [M_O]=\mathrm{N}\cdot\mathrm{m}. \]

Wektor momentu jest prostopadły do płaszczyzny wyznaczonej przez wektor położenia i siłę. Jego zwrot określa się za pomocą reguły prawej dłoni. Oś działania wektora momentu przechodzi przez punkt \(O\) i jest prostopadła do płaszczyzny obrotu. Nie musi więc być pionowa — jej pozorne ustawienie zależy od sposobu przedstawienia układu na rysunku.

Jeżeli linia działania siły przechodzi przez punkt \(O\), ramię \(d\) jest równe zeru. Wtedy również moment tej siły względem punktu \(O\) wynosi zero.

Siła jako wektor oraz wektor momentu prostopadły do płaszczyzny obrotu
Siła działa wzdłuż swojej linii działania, natomiast wektor momentu jest skierowany prostopadle do płaszczyzny obrotu.

Para sił

Parą sił nazywa się układ dwóch sił, które:

Wypadkowa siła pary jest równa zeru, jednak para wywołuje obrót ciała. Jej moment ma wartość:

\[ M=F a, \] 

gdzie \(a\) oznacza odległość między liniami działania obu sił.

Moment pary sił nie zależy od wyboru punktu, względem którego jest obliczany. Jest więc tak zwanym wektorem swobodnym i może być przenoszony w obrębie ciała sztywnego bez zmiany zewnętrznego skutku działania pary.

Układ sił i jego wypadkowa

Na ciało zazwyczaj działa jednocześnie wiele sił. Zbiór wszystkich sił działających na rozpatrywany obiekt nazywamy układem sił.

Układy sił można klasyfikować między innymi jako:

W płaskim układzie sił wszystkie siły działają w jednej płaszczyźnie. W układzie przestrzennym ich linie działania mogą być rozmieszczone w różnych płaszczyznach.

Jeżeli linie działania wszystkich sił przecinają się w jednym punkcie, mamy do czynienia z układem zbieżnym. Moment każdej z tych sił względem punktu zbieżności jest równy zeru.

Wypadkowa sił

Wypadkowa jest pojedynczą siłą wywołującą taki sam zewnętrzny skutek jak rozpatrywany układ sił — o ile zastąpienie całego układu jedną siłą jest możliwe.

\[ \vec{R}=\sum_i \vec{F}_i. \]

W układzie płaskim składowe głównego wektora sił wynoszą:

\[ R_x=\sum_i F_{ix}, \qquad R_y=\sum_i F_{iy}. \]

Jego wartość jest równa:

\[ R=\sqrt{R_x^2+R_y^2}. \]

Sama znajomość wartości i kierunku wektora \(\vec{R}\) nie zawsze wystarcza. Trzeba również ustalić jego linię działania, ponieważ od niej zależy moment wywierany na ciało.

Redukcja układu sił

Redukcja układu sił polega na zastąpieniu wielu sił prostszym układem równoważnym.

Dowolny układ sił działających na ciało sztywne można zredukować w wybranym punkcie \(O\) do:

\[ \vec{R}=\sum_i\vec{F}_i, \qquad \vec{M}_O=\sum_i\vec{M}_O(\vec{F}_i). \]

Cały układ można więc przedstawić jako siłę \(\vec{R}\) przyłożoną w punkcie \(O\) oraz dodatkowy moment \(\vec{M}_O\).

W zależności od wartości głównego wektora i głównego momentu mogą wystąpić następujące przypadki:

Ważne jest rozróżnienie między głównym wektorem sił a siłą wypadkową. Główny wektor można obliczyć dla każdego układu. Pojedyncza siła wypadkowa istnieje natomiast tylko wtedy, gdy może samodzielnie zastąpić cały układ, zachowując zarówno jego skutek postępowy, jak i obrotowy.

Układ trzech sił zredukowany w punkcie O do siły wypadkowej i momentu
Dowolny układ sił można w wybranym punkcie zredukować do głównego wektora sił i głównego momentu.

Więzy, podpory i reakcje

Ciało może mieć ograniczoną możliwość przemieszczania się lub obracania. Ograniczenia te nazywamy więzami. Mogą nimi być podpory, przeguby, utwierdzenia, prowadnice, cięgna albo połączenia z innymi elementami.

Więzy przeciwdziałają ruchowi za pomocą sił i momentów nazywanych reakcjami więzów. Podczas wykonywania schematu ciała swobodnego podporę usuwa się myślowo, a jej działanie zastępuje odpowiednią reakcją lub zestawem reakcji.

W płaskich zadaniach statyki najczęściej spotyka się następujące rodzaje więzów:

Rodzaj więzuOgraniczony ruchReakcje w układzie płaskim
Podpora przegubowa nieprzesuwnaPrzesunięcie w kierunku osi \(x\) i \(y\)\(R_{Ax}\), \(R_{Ay}\)
Podpora przegubowa przesuwnaPrzesunięcie prostopadłe do powierzchni podparciaJedna reakcja \(R_B\)
Utwierdzenie sztywneOba przesunięcia i obrót\(R_{Ax}\), \(R_{Ay}\), \(M_A\)
Prowadnica lub utwierdzenie przesuwnePrzesunięcie prostopadłe do prowadnicy oraz obrót\(R_B\), \(M_B\)
Pręt lub cięgno dwusiłoweRuch wzdłuż osi prętaSiła osiowa \(S\)

Podpora przegubowa nieprzesuwna blokuje dwa przesunięcia, ale umożliwia obrót. Dlatego w układzie płaskim występują w niej dwie składowe reakcji: \(R_{Ax}\) oraz \(R_{Ay}\).

Podpora przegubowa przesuwna blokuje ruch tylko w kierunku prostopadłym do powierzchni podparcia. Wywołuje więc jedną reakcję \(R_B\), której kierunek wynika z ustawienia powierzchni.

Utwierdzenie sztywne nie pozwala ani na przesunięcie, ani na obrót elementu. W układzie płaskim występują w nim dwie składowe siły reakcji oraz moment utwierdzenia.

Prowadnica, nazywana także utwierdzeniem przesuwnym, umożliwia przesuwanie elementu wzdłuż prowadnicy, ale blokuje ruch prostopadły i obrót. Powoduje więc powstanie jednej siły reakcji i jednego momentu.

Idealny pręt dwusiłowy przenosi jedynie siłę \(S\) działającą wzdłuż jego osi. Może być ona siłą rozciągającą albo ściskającą.

Uwaga: przedstawione reakcje dotyczą układu płaskiego. W przestrzeni ciało może przemieszczać się wzdłuż trzech osi oraz obracać wokół trzech osi. Pełne utwierdzenie przestrzenne może więc wywoływać trzy składowe siły i trzy składowe momentu.

Rodzaje podpór w układzie płaskim i odpowiadające im siły oraz momenty reakcji
Rodzaj więzu decyduje o liczbie i kierunkach reakcji występujących w układzie płaskim.

Schemat ciała swobodnego

Jednym z najważniejszych etapów rozwiązywania zadania ze statyki jest wykonanie schematu ciała swobodnego, określanego również skrótem SCS.

Polega on na odizolowaniu analizowanego ciała od otoczenia i zaznaczeniu wszystkich działających na nie sił oraz momentów. Usunięte podpory i połączenia zastępuje się odpowiadającymi im reakcjami.

Na schemacie należy umieścić:

Zwroty nieznanych reakcji można początkowo przyjąć dowolnie. Jeżeli obliczona wartość okaże się ujemna, oznacza to, że rzeczywisty zwrot reakcji jest przeciwny do przyjętego na schemacie.

Poprawny schemat ciała swobodnego jest często najważniejszą częścią rozwiązania. Pominięcie jednej siły albo niewłaściwe zastąpienie podpory prowadzi do błędnych równań równowagi.

Belka na podporze przegubowej i przesuwnej oraz jej schemat ciała swobodnego
Na schemacie ciała swobodnego podpory zastępuje się odpowiadającymi im reakcjami.

Warunki równowagi

Równowaga w układzie płaskim

Dla dowolnego płaskiego układu sił działającego na ciało sztywne trzeba zapisać trzy niezależne równania równowagi. Najczęściej wykorzystuje się dwa równania rzutów sił oraz jedno równanie momentów:

\[ \sum F_x=0, \qquad \sum F_y=0, \qquad \sum M_O=0. \]

Osie nie muszą jednak być oznaczone jako \(x\) i \(y\), nie muszą też być do siebie prostopadłe. Wystarczy wybrać dwa nierównoległe kierunki \(u\) i \(v\):

\[ \sum F_u=0, \qquad \sum F_v=0, \qquad \sum M_O=0, \qquad u\not\parallel v. \] 

Dwa pierwsze równania zapewniają, że główny wektor sił jest równy zeru. Równanie momentów wyklucza natomiast możliwość pozostania niezrównoważonej pary sił. Punkt \(O\), względem którego obliczamy momenty, może być wybrany dowolnie.

Dwa równania momentów i jedno równanie rzutów

Równoważnym zestawem mogą być również dwa równania momentów względem różnych punktów \(A\) i \(B\) oraz jedno równanie rzutów sił:

\[ \sum M_A=0, \qquad \sum M_B=0, \qquad \sum F_u=0. \]

Punkty \(A\) i \(B\) muszą być różne, a kierunek \(u\), na który rzutujemy siły, nie może być prostopadły do prostej \(AB\).

Zerowanie momentów względem punktów \(A\) i \(B\) pozostawia bowiem jeszcze możliwość istnienia niezerowej siły wypadkowej działającej wzdłuż prostej \(AB\). Równanie rzutu musi wykluczyć także tę możliwość. Gdyby kierunek rzutowania był prostopadły do \(AB\), rzut takiej siły byłby automatycznie równy zeru, a równania nie byłyby niezależne.

Trzy równania momentów

Warunki równowagi można również zapisać wyłącznie za pomocą trzech równań momentów:

\[ \sum M_A=0, \qquad \sum M_B=0, \qquad \sum M_C=0. \]

W tym przypadku punkty \(A\), \(B\) i \(C\) nie mogą leżeć na jednej prostej. Warunek ten zapewnia niezależność trzech równań i wyklucza zarówno niezrównoważoną siłę, jak i niezrównoważony moment.

Dobór równań: wszystkie przedstawione zestawy są równoważnymi warunkami równowagi płaskiego układu sił, o ile spełnione są podane warunki ich niezależności. W praktyce wybiera się taki zestaw osi i punktów, który upraszcza obliczenia. Szczególnie korzystne jest obliczanie momentów względem punktów, przez które przechodzą linie działania kilku nieznanych reakcji, ponieważ ich momenty są wtedy równe zeru.

Najczęściej stosowany zestaw dwóch równań rzutów i jednego równania momentów jest zwykle najbardziej przejrzysty, ale w wielu zadaniach odpowiedni dobór dwóch albo trzech równań momentów pozwala znacznie skrócić obliczenia.

Przedstawione wcześniej równoważne zestawy trzech równań dotyczą płaskiego układu sił działającego na ciało sztywne. Dla punktu materialnego oraz układu przestrzennego liczba i postać potrzebnych warunków równowagi są inne.

Równowaga punktu materialnego

W przypadku punktu materialnego nie rozpatruje się jego obrotu. W układzie płaskim wystarczają więc dwa równania:

\[ \sum F_x=0, \qquad \sum F_y=0. \]

Równowaga w przestrzeni

W przestrzeni ciało może przemieszczać się w trzech niezależnych kierunkach i obracać wokół trzech osi. Warunki równowagi obejmują więc sześć równań:

\[ \sum F_x=0, \qquad \sum F_y=0, \qquad \sum F_z=0, \] \[ \sum M_x=0, \qquad \sum M_y=0, \qquad \sum M_z=0. \]

Obciążenia skupione i rozłożone

Siła skupiona jest idealizacją obciążenia działającego w jednym punkcie. W rzeczywistości każde obciążenie jest rozłożone na pewnej powierzchni, ale jeżeli obszar jego działania jest niewielki w porównaniu z wymiarami ciała, można zastąpić je pojedynczą siłą.

Obciążenie rozłożone działa na określonej długości, powierzchni albo objętości. Przykładem może być ciężar materiału rozmieszczonego na belce, obciążenie śniegiem działające na dach albo parcie cieczy.

Dla obciążenia ciągłego o intensywności \(q(x)\) wypadkowa ma wartość:

\[ R_q=\int_a^b q(x)\,dx. \]

Jej linia działania przechodzi przez środek ciężkości pola znajdującego się pod wykresem intensywności obciążenia. Dla obciążenia równomiernego o wartości \(q\), działającego na długości \(l\), otrzymujemy:

\[ R_q=q\,l, \]

a siła wypadkowa jest przyłożona w połowie obciążonego odcinka.

Statyczna wyznaczalność układu

Równania równowagi pozwalają wyznaczać nieznane reakcje oraz siły. Nie zawsze jednak ich liczba jest wystarczająca.

Jeżeli wszystkie niewiadome można obliczyć wyłącznie za pomocą równań statyki, układ nazywamy statycznie wyznaczalnym.

Jeżeli liczba niewiadomych jest większa od liczby niezależnych równań równowagi, układ jest statycznie niewyznaczalny. Do jego rozwiązania potrzebne są dodatkowe zależności dotyczące odkształceń, przemieszczeń i właściwości materiału.

Samo porównanie liczby niewiadomych z liczbą równań nie zawsze wystarcza. Trzeba również sprawdzić, czy więzy są właściwie rozmieszczone i czy równania są od siebie niezależne. Układ z niewystarczającą albo nieprawidłowo rozmieszczoną liczbą więzów może być geometrycznie zmienny i zachowywać się jak mechanizm.

Jak rozwiązywać zadania ze statyki?

Rozwiązanie większości podstawowych zadań można uporządkować w następujących etapach:

  1. Wybrać ciało albo fragment konstrukcji przeznaczony do analizy.
  2. Odizolować go od otoczenia.
  3. Zastąpić więzy odpowiednimi reakcjami.
  4. Zaznaczyć wszystkie siły, momenty, wymiary i kąty.
  5. Przyjąć układ współrzędnych i dodatni zwrot momentów.
  6. Rozłożyć siły ukośne na składowe.
  7. Zapisać równania równowagi.
  8. Rozwiązać otrzymany układ równań.
  9. Sprawdzić jednostki, znaki i fizyczny sens wyników.

W bardziej złożonych konstrukcjach trzeba kolejno odcinać poszczególne elementy i dla każdego z nich wykonywać osobny schemat ciała swobodnego.

Najczęstsze błędy

Zastosowania statyki

Statyka jest wykorzystywana wszędzie tam, gdzie trzeba ocenić równowagę i obciążenie elementów. Jej metody stosuje się między innymi do:

Podsumowanie

Statyka zajmuje się równowagą ciał poddanych działaniu sił. Jej podstawą są trzy zasady dynamiki Newtona: zasada bezwładności, zasada wiążąca siłę z przyspieszeniem oraz zasada akcji i reakcji.

Do najważniejszych pojęć statyki należą:

W płaskim układzie sił ciało sztywne pozostaje w równowadze, jeżeli spełnione są trzy równania:

\[ \sum F_x=0, \qquad \sum F_y=0, \qquad \sum M_O=0. \]

Poprawne rozpoznanie wszystkich sił, wykonanie schematu ciała swobodnego i właściwe obliczanie momentów stanowią podstawę rozwiązywania niemal każdego zadania ze statyki.

Rzut ukośny – tor ruchu, czas lotu, zasięg i maksymalna wysokość

Rzut ukośny jest ruchem ciała, któremu nadano prędkość początkową skierowaną pod pewnym kątem do poziomu. Przy pominięciu oporu powietrza ruch można rozłożyć na jednostajny ruch poziomy oraz ruch pionowy ze stałym przyspieszeniem ziemskim. Połączenie obu składowych daje paraboliczny tor ruchu.

Rzut ukośny jest rozwinięciem zagadnień omówionych w artykułach Rzut poziomy – tor ruchu, czas lotu, zasięg i prędkość oraz Swobodne spadanie i rzut pionowy – prędkość, czas lotu i maksymalna wysokość.

Czym jest rzut ukośny?

Rzut ukośny zachodzi wtedy, gdy ciało otrzymuje prędkość początkową \(\vec v_0\) skierowaną pod kątem \(\alpha\) do poziomu. Po rozpoczęciu ruchu, przy pominięciu oporu powietrza, na ciało działa wyłącznie siła grawitacji.

Prędkość początkowa ma dwie składowe:

W kierunku poziomym ciało porusza się ruchem jednostajnym. W kierunku pionowym jego ruch jest ruchem jednostajnie zmiennym z przyspieszeniem skierowanym w dół.

Założenia modelu

W podstawowym modelu rzutu ukośnego przyjmujemy, że:

Wartość przyspieszenia ziemskiego wynosi w przybliżeniu:

\[ g\approx 9{,}81\ \mathrm{m/s^2}. \]

W prostszych zadaniach często przyjmuje się:

\[ g\approx 10\ \mathrm{m/s^2}. \]

Układ współrzędnych

Przyjmujemy początek układu współrzędnych w punkcie wyrzutu. Oś \(x\) jest skierowana poziomo w prawo, a oś \(y\) pionowo w górę.

Składowe przyspieszenia wynoszą wtedy:

\[ a_x=0, \qquad a_y=-g. \]

Wektor przyspieszenia można zapisać jako:

\[ \vec a=-g\vec e_y. \] 

Znak minus oznacza, że przyspieszenie jest skierowane przeciwnie do dodatniego zwrotu osi \(y\).

Rozkład prędkości początkowej

Wektor prędkości początkowej \(\vec v_0\) tworzy z poziomem kąt \(\alpha\). Rozkładamy go na dwie wzajemnie prostopadłe składowe:

\[ \vec v_0 = \vec v_{0x}+\vec v_{0y}. \]

Składowa pozioma wynosi:

\[ \boxed{ v_{0x}=v_0\cos\alpha } \]

Składowa pionowa:

\[ \boxed{ v_{0y}=v_0\sin\alpha } \]

Ponieważ składowe są prostopadłe, wartość prędkości początkowej spełnia zależność:

\[ v_0 = \sqrt{v_{0x}^2+v_{0y}^2}. \]
Rozkład prędkości początkowej w rzucie ukośnym na składową poziomą i pionową
Prędkość początkowa \(\vec v_0\) ma składową poziomą \(v_{0x}=v_0\cos\alpha\) oraz pionową \(v_{0y}=v_0\sin\alpha\).

Równania ruchu

Ruch w kierunku poziomym

W kierunku poziomym przyspieszenie jest równe zeru:

\[ a_x=0. \]

Pozioma składowa prędkości pozostaje stała:

\[ v_x=v_0\cos\alpha. \]

Położenie poziome:

\[ \boxed{ x(t)=v_0\cos\alpha\cdot t } \] 

Ruch w kierunku pionowym

W kierunku pionowym przyspieszenie jest stałe:

\[ a_y=-g. \]

Pionowa składowa prędkości zmienia się zgodnie z równaniem:

\[ \boxed{ v_y(t)=v_0\sin\alpha-gt } \]

Położenie pionowe, gdy punkt wyrzutu przyjęto jako \(y=0\), wynosi:

\[ \boxed{ y(t) = v_0\sin\alpha\cdot t-\frac{gt^2}{2} } \]

W zapisie wektorowym:

\[ \vec r(t) = v_0\cos\alpha\cdot t\,\vec e_x + \left( v_0\sin\alpha\cdot t-\frac{gt^2}{2} \right)\vec e_y, \] \[ \vec v(t) = v_0\cos\alpha\,\vec e_x + \left( v_0\sin\alpha-gt \right)\vec e_y, \] \[ \vec a(t)=-g\vec e_y. \] 

Równanie toru rzutu ukośnego

Równanie toru otrzymujemy przez wyeliminowanie czasu z równań położenia.

Z równania:

\[ x=v_0\cos\alpha\cdot t \] 

wyznaczamy czas:

\[ t = \frac{x}{v_0\cos\alpha}. \]

Po podstawieniu do równania pionowego:

\[ y = v_0\sin\alpha \frac{x}{v_0\cos\alpha} - \frac{g}{2} \left( \frac{x}{v_0\cos\alpha} \right)^2. \]

Otrzymujemy:

\[ \boxed{ y(x) = x\tan\alpha - \frac{gx^2}{2v_0^2\cos^2\alpha} } \]

Jest to równanie paraboli skierowanej ramionami w dół. Przy braku oporu powietrza tor rzutu ukośnego jest więc paraboliczny.

Składowe prędkości podczas lotu

Pozioma składowa prędkości przez cały czas pozostaje stała:

\[ v_x=v_0\cos\alpha=\mathrm{const}. \]

Pionowa składowa prędkości maleje liniowo:

\[ v_y=v_0\sin\alpha-gt. \]

Podczas wznoszenia:

\[ v_y>0. \]

W najwyższym punkcie:

\[ v_y=0. \]

Podczas opadania:

\[ v_y<0. \]

Wartość prędkości w dowolnej chwili wynosi:

\[ v = \sqrt{v_x^2+v_y^2}, \] \[ \boxed{ v = \sqrt{ v_0^2\cos^2\alpha + \left(v_0\sin\alpha-gt\right)^2 } } \]
Składowe prędkości podczas wznoszenia, w najwyższym punkcie i podczas opadania w rzucie ukośnym
Składowa pozioma \(v_x\) pozostaje stała. Podczas wznoszenia \(v_y>0\), w najwyższym punkcie \(v_y=0\), a podczas opadania \(v_y<0\).

Najwyższy punkt toru

W najwyższym punkcie pionowa składowa prędkości jest równa zeru:

\[ v_y=0. \]

Nie oznacza to jednak zatrzymania ciała. Pozioma składowa prędkości nadal wynosi:

\[ v_x=v_0\cos\alpha. \]

Wartość prędkości w najwyższym punkcie jest więc równa:

\[ \boxed{ v_{\mathrm{szczyt}} = v_0\cos\alpha } \]

W najwyższym punkcie nie znikają ani prędkość, ani przyspieszenie. Zerowa jest jedynie pionowa składowa prędkości. Ciało nadal porusza się poziomo, a przyspieszenie pozostaje równe:

\[ a_y=-g. \]

Czas wznoszenia

Czas wznoszenia wyznaczamy z warunku:

\[ v_y=0. \]

Z równania:

\[ 0 = v_0\sin\alpha-gt_{\mathrm{w}} \]

otrzymujemy:

\[ \boxed{ t_{\mathrm{w}} = \frac{v_0\sin\alpha}{g} } \]

Maksymalna wysokość

Maksymalną wysokość ponad poziom wyrzutu można wyznaczyć przez podstawienie czasu wznoszenia do równania położenia pionowego:

\[ h_{\max} = v_0\sin\alpha\cdot t_{\mathrm{w}} - \frac{gt_{\mathrm{w}}^2}{2}. \]

Po podstawieniu:

\[ \boxed{ h_{\max} = \frac{v_0^2\sin^2\alpha}{2g} } \]

Wzór ten określa wysokość ponad poziom punktu wyrzutu. Jeżeli punkt wyrzutu znajduje się na wysokości \(y_0\), najwyższe położenie wynosi:

\[ y_{\max} = y_0+\frac{v_0^2\sin^2\alpha}{2g}. \]

Całkowity czas lotu

Jeżeli ciało wraca na ten sam poziom, z którego zostało wyrzucone, w chwili końcowej:

\[ y=0. \]

Z równania:

\[ 0 = v_0\sin\alpha\cdot t-\frac{gt^2}{2} \]

otrzymujemy:

\[ t \left( v_0\sin\alpha-\frac{gt}{2} \right) = 0. \] 

Rozwiązanie \(t=0\) odpowiada chwili wyrzutu. Drugie rozwiązanie jest całkowitym czasem lotu:

\[ \boxed{ t_{\mathrm{lotu}} = \frac{2v_0\sin\alpha}{g} } \]

Przy jednakowym poziomie początkowym i końcowym czas opadania jest równy czasowi wznoszenia:

\[ t_{\mathrm{opadania}} = t_{\mathrm{w}}, \] \[ t_{\mathrm{lotu}} = 2t_{\mathrm{w}}. \]

Zasięg rzutu ukośnego

Zasięg \(L\) jest poziomą odległością pomiędzy punktem wyrzutu i punktem upadku.

W kierunku poziomym ciało porusza się ze stałą prędkością:

\[ v_x=v_0\cos\alpha. \]

Zatem:

\[ L = v_0\cos\alpha\cdot t_{\mathrm{lotu}}. \]

Po podstawieniu czasu lotu:

\[ L = v_0\cos\alpha \frac{2v_0\sin\alpha}{g}. \]

Korzystamy z tożsamości:

\[ 2\sin\alpha\cos\alpha = \sin2\alpha. \]

Ostatecznie:

\[ \boxed{ L = \frac{v_0^2\sin2\alpha}{g} } \]

Wzór ten obowiązuje wyłącznie wtedy, gdy punkt wyrzutu i punkt upadku znajdują się na tej samej wysokości.

Symetria ruchu i kąt uderzenia

Jeżeli ciało wraca na poziom wyrzutu i pomijamy opór powietrza, tor jest symetryczny względem pionowej prostej przechodzącej przez najwyższy punkt.

W chwili wyrzutu:

\[ v_{x0}=v_0\cos\alpha, \qquad v_{y0}=v_0\sin\alpha. \]

W chwili powrotu na poziom początkowy:

\[ v_x=v_0\cos\alpha, \] \[ v_y=-v_0\sin\alpha. \]

Wartość prędkości podczas uderzenia wynosi więc:

\[ v_{\mathrm{uderzenia}} = \sqrt{ v_0^2\cos^2\alpha + v_0^2\sin^2\alpha }, \] \[ \boxed{ v_{\mathrm{uderzenia}}=v_0 } \]

Jeżeli \(\beta\) oznacza wartość kąta uderzenia mierzonego poniżej poziomu, to:

\[ \tan\beta = \frac{|v_y|}{v_x} = \frac{v_0\sin\alpha}{v_0\cos\alpha} = \tan\alpha. \]

Zatem:

\[ \boxed{ \beta=\alpha } \]

Kąt uderzenia ma więc taką samą wartość jak kąt wyrzutu, lecz jest skierowany poniżej poziomu.

Symetryczny tor rzutu ukośnego z czasem lotu, maksymalną wysokością, zasięgiem oraz jednakowym kątem wyrzutu i uderzenia
Przy braku oporu powietrza i jednakowym poziomie wyrzutu oraz upadku tor jest symetryczny. Czas wznoszenia jest równy czasowi opadania, a kąt uderzenia ma taką samą wartość jak kąt wyrzutu.

Wpływ kąta rzutu na zasięg

Dla ustalonej wartości prędkości początkowej i jednakowych poziomów wyrzutu oraz upadku zasięg zależy od wyrażenia:

\[ \sin2\alpha. \]

Największa możliwa wartość sinusa wynosi:

\[ \sin2\alpha=1. \]

Warunek ten jest spełniony, gdy:

\[ 2\alpha=90^\circ, \] \[ \boxed{ \alpha=45^\circ } \]

W idealnym modelu największy zasięg występuje więc dla kąta \(45^\circ\).

Kąty dopełniające się

Kąty \(\alpha\) oraz \(90^\circ-\alpha\) dają ten sam zasięg, ponieważ:

\[ \sin\left[2(90^\circ-\alpha)\right] = \sin(180^\circ-2\alpha) = \sin2\alpha. \]

Przykładowo:

\[ L_{30^\circ}=L_{60^\circ}. \]

Tor dla \(60^\circ\) osiąga większą wysokość i wiąże się z dłuższym czasem lotu. Tor dla \(30^\circ\) jest niższy i bardziej płaski, ale oba rzuty mają taki sam zasięg.

Porównanie torów rzutu ukośnego dla kątów 30, 45 i 60 stopni
Dla tej samej prędkości początkowej i jednakowych poziomów startu oraz lądowania kąty \(30^\circ\) i \(60^\circ\) dają ten sam zasięg, a największy zasięg występuje dla \(45^\circ\).

Kąt \(45^\circ\) nie zawsze daje największy rzeczywisty zasięg. Wynik ten obowiązuje dla braku oporu powietrza oraz jednakowych wysokości punktu wyrzutu i punktu upadku. Przy oporze powietrza lub różnych wysokościach optymalny kąt może być inny.

Wykresy ruchu

Położenie poziome \(x(t)\)

\[ x(t)=v_0\cos\alpha\cdot t. \] 

Wykres jest linią prostą. Jego nachylenie jest równe stałej składowej poziomej prędkości.

Położenie pionowe \(y(t)\)

\[ y(t) = v_0\sin\alpha\cdot t-\frac{gt^2}{2}. \]

Wykres jest parabolą skierowaną ramionami w dół. Położenie początkowo rośnie, osiąga maksimum, a następnie maleje.

Prędkość pozioma \(v_x(t)\)

\[ v_x(t)=v_0\cos\alpha. \]

Wykres jest poziomą linią, ponieważ składowa pozioma prędkości jest stała.

Prędkość pionowa \(v_y(t)\)

\[ v_y(t)=v_0\sin\alpha-gt. \]

Wykres jest linią prostą o nachyleniu \(-g\). W najwyższym punkcie przecina oś czasu, ponieważ wtedy \(v_y=0\).

Wykresy położenia poziomego i pionowego oraz składowych prędkości w rzucie ukośnym
Położenie \(x(t)\) rośnie liniowo, \(y(t)\) zmienia się parabolicznie, \(v_x(t)\) jest stałe, a \(v_y(t)=v_0\sin\alpha-gt\) maleje liniowo.

Rzut ukośny z wysokości

Jeżeli punkt upadku znajduje się na innej wysokości niż punkt wyrzutu, ruch nie jest symetryczny względem najwyższego punktu. Nie można wtedy stosować wzorów:

\[ t_{\mathrm{lotu}} = \frac{2v_0\sin\alpha}{g}, \qquad L = \frac{v_0^2\sin2\alpha}{g}. \]

Należy skorzystać z ogólnego równania położenia:

\[ y(t) = y_0+v_0\sin\alpha\cdot t-\frac{gt^2}{2}. \]

Jeżeli ciało zostało wyrzucone z wysokości \(h\) nad poziomem upadku, a początek osi umieszczono w punkcie wyrzutu, poziom końcowy ma współrzędną:

\[ y=-h. \]

Warunek uderzenia ma postać:

\[ -h = v_0\sin\alpha\cdot t-\frac{gt^2}{2}. \]

Dodatnie rozwiązanie równania kwadratowego wynosi:

\[ \boxed{ t_{\mathrm{lotu}} = \frac{ v_0\sin\alpha + \sqrt{ v_0^2\sin^2\alpha+2gh } }{g} } \]

Zasięg obliczamy następnie ze wzoru:

\[ L = v_0\cos\alpha\cdot t_{\mathrm{lotu}}. \]

W takim przypadku czas opadania jest dłuższy od czasu wznoszenia, kąt uderzenia nie musi być równy kątowi wyrzutu, a kąt \(45^\circ\) nie musi dawać największego zasięgu.

Wpływ oporu powietrza

W rzeczywistym ruchu opór powietrza powoduje zmniejszanie zarówno poziomej, jak i pionowej składowej prędkości. W konsekwencji:

Znaczenie oporu powietrza zależy między innymi od kształtu i rozmiaru ciała, jego masy, prędkości oraz gęstości powietrza.

Przykłady obliczeniowe

Przykład 1. Rzut pod kątem \(30^\circ\)

Ciało wyrzucono z poziomu podłoża z prędkością:

\[ v_0=20\ \mathrm{m/s} \]

pod kątem:

\[ \alpha=30^\circ. \]

Składowe prędkości początkowej:

\[ v_{0x} = 20\cos30^\circ \approx 17{,}32\ \mathrm{m/s}, \] \[ v_{0y} = 20\sin30^\circ = 10\ \mathrm{m/s}. \]

Czas wznoszenia:

\[ t_{\mathrm{w}} = \frac{10}{9{,}81} \approx 1{,}02\ \mathrm{s}. \]

Maksymalna wysokość:

\[ h_{\max} = \frac{10^2}{2\cdot9{,}81} \approx 5{,}10\ \mathrm{m}. \]

Całkowity czas lotu:

\[ t_{\mathrm{lotu}} = 2t_{\mathrm{w}} \approx 2{,}04\ \mathrm{s}. \]

Zasięg:

\[ L = \frac{20^2\sin60^\circ}{9{,}81} \approx 35{,}3\ \mathrm{m}. \]

W chwili uderzenia ciało ma szybkość \(20\ \mathrm{m/s}\), a jego prędkość jest skierowana pod kątem \(30^\circ\) poniżej poziomu.

Przykład 2. Porównanie kątów \(30^\circ\) i \(60^\circ\)

Dwa ciała wyrzucono z taką samą prędkością \(v_0=20\ \mathrm{m/s}\), odpowiednio pod kątem \(30^\circ\) i \(60^\circ\).

Zasięg dla kąta \(30^\circ\):

\[ L_{30^\circ} = \frac{20^2\sin60^\circ}{9{,}81} \approx 35{,}3\ \mathrm{m}. \]

Zasięg dla kąta \(60^\circ\):

\[ L_{60^\circ} = \frac{20^2\sin120^\circ}{9{,}81} \approx 35{,}3\ \mathrm{m}. \]

Zasięgi są takie same, ale maksymalne wysokości różnią się.

\[ h_{30^\circ} \approx 5{,}10\ \mathrm{m}, \] \[ h_{60^\circ} = \frac{20^2\sin^260^\circ}{2\cdot9{,}81} \approx 15{,}3\ \mathrm{m}. \]

Rzut pod kątem \(60^\circ\) trwa dłużej i osiąga większą wysokość, ale ma taki sam zasięg jak rzut pod kątem \(30^\circ\).

Przykład 3. Rzut ukośny z wysokości

Ciało wyrzucono z wysokości \(10\ \mathrm{m}\) z prędkością:

\[ v_0=15\ \mathrm{m/s} \]

pod kątem:

\[ \alpha=40^\circ. \]

Czas lotu:

\[ t_{\mathrm{lotu}} = \frac{ 15\sin40^\circ + \sqrt{ 15^2\sin^240^\circ+2\cdot9{,}81\cdot10 } }{9{,}81} \approx 2{,}72\ \mathrm{s}. \]

Składowa pozioma prędkości:

\[ v_x = 15\cos40^\circ \approx 11{,}49\ \mathrm{m/s}. \]

Zasięg:

\[ L = v_xt_{\mathrm{lotu}} \approx 11{,}49\cdot2{,}72 \approx 31{,}2\ \mathrm{m}. \]

Ponieważ punkt upadku znajduje się niżej niż punkt wyrzutu, czas opadania jest dłuższy od czasu wznoszenia, a kąt uderzenia jest większy od kąta wyrzutu.

Jednostki i kontrola wymiarów

WielkośćOznaczenieJednostka SI
Położenie i wysokość\(x\), \(y\), \(h\)\(\mathrm{m}\)
Zasięg\(L\)\(\mathrm{m}\)
Czas\(t\)\(\mathrm{s}\)
Prędkość\(v_0\), \(v_x\), \(v_y\)\(\mathrm{m/s}\)
Przyspieszenie ziemskie\(g\)\(\mathrm{m/s^2}\)
Kąt\(\alpha\)\(\mathrm{rad}\) lub \(^\circ\)

Kontrola wymiarów wzoru na czas lotu:

\[ \left[ \frac{2v_0\sin\alpha}{g} \right] = \frac{\mathrm{m/s}}{\mathrm{m/s^2}} = \mathrm{s}. \]

Kontrola wzoru na maksymalną wysokość:

\[ \left[ \frac{v_0^2\sin^2\alpha}{2g} \right] = \frac{\mathrm{m^2/s^2}}{\mathrm{m/s^2}} = \mathrm{m}. \]

Kontrola wzoru na zasięg:

\[ \left[ \frac{v_0^2\sin2\alpha}{g} \right] = \frac{\mathrm{m^2/s^2}}{\mathrm{m/s^2}} = \mathrm{m}. \]

Najczęstsze błędy

Przed wykonaniem obliczeń należy sprawdzić tryb kalkulatora. Jeżeli kąt podano w stopniach, kalkulator musi pracować w trybie stopniowym, zwykle oznaczonym jako DEG.

Podsumowanie

Prędkość początkową w rzucie ukośnym rozkładamy na składowe:

\[ v_{0x}=v_0\cos\alpha, \qquad v_{0y}=v_0\sin\alpha. \]

Równania ruchu:

\[ x(t)=v_0\cos\alpha\cdot t, \] \[ y(t) = v_0\sin\alpha\cdot t-\frac{gt^2}{2}. \]

Składowe prędkości:

\[ v_x=v_0\cos\alpha, \] \[ v_y=v_0\sin\alpha-gt. \]

Tor ruchu jest parabolą:

\[ y(x) = x\tan\alpha - \frac{gx^2}{2v_0^2\cos^2\alpha}. \]

Przy jednakowych poziomach wyrzutu i upadku:

\[ t_{\mathrm{w}} = \frac{v_0\sin\alpha}{g}, \] \[ h_{\max} = \frac{v_0^2\sin^2\alpha}{2g}, \] \[ t_{\mathrm{lotu}} = \frac{2v_0\sin\alpha}{g}, \] \[ L = \frac{v_0^2\sin2\alpha}{g}. \]

W idealnym modelu największy zasięg występuje dla kąta \(45^\circ\), a kąty dopełniające się do \(90^\circ\) dają taki sam zasięg. Tor jest symetryczny, szybkość w chwili uderzenia jest równa szybkości początkowej, a kąt uderzenia ma taką samą wartość jak kąt wyrzutu.

Swobodne spadanie i rzut pionowy – prędkość, czas lotu i maksymalna wysokość

Swobodne spadanie i rzut pionowy są szczególnymi przypadkami ruchu prostoliniowego ze stałym przyspieszeniem grawitacyjnym. Ich opis zależy przede wszystkim od przyjętego zwrotu osi pionowej oraz od poprawnego rozróżnienia prędkości algebraicznej i szybkości, czyli wartości bezwzględnej prędkości.

Przy pominięciu oporu powietrza wszystkie ciała znajdujące się w pobliżu powierzchni Ziemi poruszają się z takim samym przyspieszeniem skierowanym pionowo w dół. W rzucie pionowym w górę prędkość początkowo maleje, w najwyższym punkcie jest chwilowo równa zeru, a następnie zmienia znak i ciało zaczyna spadać.

Ruch pionowy korzysta z tych samych zależności, które omówiliśmy w artykule Ruch prostoliniowy – położenie, prędkość, przyspieszenie i równania ruchu. Tutaj zastosujemy je do ruchu pod wpływem przyspieszenia ziemskiego.

Czym jest ruch pionowy?

Ruchem pionowym nazywamy ruch ciała wzdłuż pionowej prostej. Jeżeli pomijamy opór powietrza, jedynym rozpatrywanym przyspieszeniem jest przyspieszenie grawitacyjne skierowane ku powierzchni Ziemi.

Do podstawowych przypadków ruchu pionowego należą:

We wszystkich trzech przypadkach ruch odbywa się wzdłuż tej samej pionowej osi. Różne są natomiast prędkość początkowa oraz jej zwrot.

Założenia modelu

W podstawowych zadaniach dotyczących swobodnego spadania i rzutu pionowego przyjmujemy, że:

W rzeczywistym ruchu opór powietrza zależy między innymi od prędkości, kształtu ciała i gęstości powietrza. W tym artykule rozpatrujemy jednak wyidealizowany model bez oporu.

Przyspieszenie ziemskie

W pobliżu powierzchni Ziemi wartość przyspieszenia grawitacyjnego wynosi w przybliżeniu:

\[ g\approx 9{,}81\ \mathrm{m/s^2}. \]

W prostszych zadaniach szkolnych często stosuje się przybliżenie:

\[ g\approx 10\ \mathrm{m/s^2}. \]

Symbol \(g\) oznacza dodatnią wartość przyspieszenia ziemskiego. Znak algebraiczny przyspieszenia zależy natomiast od przyjętego zwrotu osi.

Wybór zwrotu osi pionowej

W ruchu pionowym można wybrać oś skierowaną dodatnio w górę albo dodatnio w dół. Oba sposoby są poprawne, jeżeli konsekwentnie stosujemy wynikające z nich znaki.

Oś skierowana dodatnio w górę

Jeżeli dodatni zwrot osi \(y\) jest skierowany w górę, przyspieszenie ziemskie ma wartość algebraiczną:

\[ a_y=-g. \]

Prędkość skierowana w górę jest dodatnia, a skierowana w dół ujemna.

Oś skierowana dodatnio w dół

Jeżeli dodatni zwrot osi jest skierowany w dół, przyspieszenie ma znak dodatni:

\[ a=+g. \]

W takim układzie prędkość spadającego ciała także jest dodatnia.

Znak minus nie jest częścią wartości \(g\). Zapis \(a=-g\) oznacza, że przyspieszenie jest skierowane przeciwnie do dodatniego zwrotu osi. Sama wartość \(g\) pozostaje dodatnia.

Prędkość a szybkość

W ruchu jednowymiarowym często posługujemy się algebraiczną składową prędkości \(v_y\). Jej znak wskazuje zwrot ruchu:

Szybkość jest natomiast wartością bezwzględną prędkości i nigdy nie jest ujemna:

\[ v_{\mathrm{szybkość}} = |\vec v| = |v_y|. \]

Przy osi skierowanej w górę podczas swobodnego spadania:

\[ v_y=-gt, \]

ale szybkość ciała wynosi:

\[ |\vec v| = |v_y| = gt. \] 

W wielu zadaniach szkolnych słowo „prędkość” jest używane potocznie w znaczeniu szybkości. Wtedy poprawną odpowiedzią może być dodatnia wartość, na przykład \(20\ \mathrm{m/s}\), wraz z informacją, że ciało porusza się w dół.

Najbardziej jednoznaczna odpowiedź zawiera wartość i kierunek. Zamiast pisać wyłącznie \(v=-20\ \mathrm{m/s}\), można podać: „prędkość ma składową \(v_y=-20\ \mathrm{m/s}\), czyli ciało porusza się w dół z szybkością \(20\ \mathrm{m/s}\)”.

Ogólne równania ruchu pionowego

Przyjmijmy oś \(y\) skierowaną dodatnio w górę. Wtedy:

\[ a_y=-g. \]

Jeżeli w chwili \(t=0\) ciało znajduje się na wysokości \(y_0\) i ma algebraiczną prędkość początkową \(v_{0y}\), równania ruchu przyjmują postać:

\[ v_y(t) = v_{0y}-gt, \] \[ y(t) = y_0+v_{0y}t-\frac{gt^2}{2}, \] \[ v_y^2 = v_{0y}^2-2g(y-y_0). \]

Wielkość \(v_{0y}\) jest tutaj składową algebraiczną:

Porównanie swobodnego spadania, rzutu pionowego w dół i rzutu pionowego w górę wraz z równaniami ruchu
Przy osi \(y\) skierowanej dodatnio w górę we wszystkich przypadkach \(a_y=-g\). Różne są natomiast wartość i zwrot prędkości początkowej.

Swobodne spadanie

Swobodne spadanie zachodzi wtedy, gdy ciało zostaje puszczone bez prędkości początkowej:

\[ v_{0y}=0. \]

Przy osi skierowanej dodatnio w górę:

\[ v_y=-gt, \] \[ y=y_0-\frac{gt^2}{2}. \]

Jeżeli przez \(h\) oznaczymy dodatnią odległość przebytą w dół:

\[ h=y_0-y, \]

to otrzymujemy wygodne zależności dotyczące wartości ruchu:

\[ h=\frac{gt^2}{2}, \] \[ |\vec v|=gt, \] \[ |\vec v|^2=2gh. \]

W tych wzorach szybkość i droga są dodatnie. Kierunek ruchu określamy osobno jako pionowo w dół.

Rzut pionowy w dół

W rzucie pionowym w dół ciało otrzymuje początkową prędkość skierowaną zgodnie z przyspieszeniem ziemskim.

Niech \(v_0>0\) oznacza dodatnią szybkość początkową. Przy osi skierowanej w górę algebraiczna prędkość początkowa wynosi:

\[ v_{0y}=-v_0. \]

Równania ruchu:

\[ v_y=-v_0-gt, \] \[ y=y_0-v_0t-\frac{gt^2}{2}. \]

Szybkość ciała rośnie zgodnie z zależnością:

\[ |\vec v| = v_0+gt. \]

Dodatnią odległość przebytą w dół można zapisać jako:

\[ h = v_0t+\frac{gt^2}{2}. \]

Po wyeliminowaniu czasu:

\[ |\vec v|^2 = v_0^2+2gh. \]

Tak zapisane wzory operują dodatnią szybkością oraz dodatnią drogą mierzoną w dół. Są dlatego często spotykane w zadaniach praktycznych.

Rzut pionowy w górę

W rzucie pionowym w górę prędkość początkowa jest skierowana przeciwnie do przyspieszenia ziemskiego. Przy osi dodatniej w górę:

\[ v_{0y}=v_0>0. \]

Prędkość zmienia się zgodnie z równaniem:

\[ v_y=v_0-gt. \]

Położenie ciała:

\[ y=y_0+v_0t-\frac{gt^2}{2}. \]

Zależność bez czasu:

\[ v_y^2 = v_0^2-2g(y-y_0). \]

Podczas wznoszenia:

\[ v_y>0, \qquad a_y=-g. \]

Prędkość i przyspieszenie mają przeciwne znaki, dlatego szybkość maleje.

Podczas spadania:

\[ v_y<0, \qquad a_y=-g. \]

Prędkość i przyspieszenie mają ten sam znak, dlatego szybkość rośnie.

Wznoszenie, najwyższy punkt i spadanie z zaznaczonymi znakami prędkości oraz przyspieszenia
Przy osi skierowanej w górę podczas wznoszenia \(v_y>0\), w najwyższym punkcie \(v_y=0\), a podczas spadania \(v_y<0\). W całym ruchu \(a_y=-g\).

Najwyższy punkt ruchu

W najwyższym punkcie rzutu pionowego w górę prędkość chwilowa jest równa zeru:

\[ v_y=0. \]

Z równania:

\[ v_y=v_0-gt \]

otrzymujemy czas wznoszenia:

\[ t_{\mathrm{w}} = \frac{v_0}{g}. \]

W najwyższym punkcie przyspieszenie nie jest równe zeru. Nadal wynosi:

\[ a_y=-g. \]

Zerowa jest tylko chwilowa prędkość. To właśnie niezerowe przyspieszenie powoduje, że po osiągnięciu najwyższego punktu ciało rozpoczyna ruch w dół.

Maksymalna wysokość

Maksymalny przyrost wysokości ponad poziom wyrzutu wyznaczamy z równania:

\[ v_y^2 = v_0^2-2g(y-y_0). \]

W najwyższym punkcie \(v_y=0\), więc:

\[ 0 = v_0^2-2gh_{\max}. \]

Stąd:

\[ h_{\max} = \frac{v_0^2}{2g}. \]

Bezwzględna wysokość najwyższego punktu wynosi:

\[ y_{\max} = y_0+h_{\max} = y_0+\frac{v_0^2}{2g}. \]

Czas lotu i powrót na poziom wyrzutu

Jeżeli ciało wraca na ten sam poziom, z którego zostało wyrzucone, spełniony jest warunek:

\[ y=y_0. \]

Z równania położenia:

\[ y_0 = y_0+v_0t-\frac{gt^2}{2} \]

otrzymujemy:

\[ t\left(v_0-\frac{gt}{2}\right)=0. \]

Pierwsze rozwiązanie \(t=0\) odpowiada chwili wyrzutu. Drugie daje całkowity czas lotu:

\[ t_{\mathrm{lotu}} = \frac{2v_0}{g}. \]

Czas lotu jest dwukrotnie większy od czasu wznoszenia:

\[ t_{\mathrm{lotu}} = 2t_{\mathrm{w}}. \]

Po powrocie na poziom wyrzutu algebraiczna prędkość wynosi:

\[ v_y=-v_0. \]

Szybkość ma więc taką samą wartość jak w chwili wyrzutu:

\[ |\vec v|=v_0, \]

ale prędkość jest skierowana w dół.

Schemat rzutu pionowego w górę z zaznaczoną maksymalną wysokością, czasem wznoszenia i czasem całego lotu
Jeżeli ciało wraca na poziom wyrzutu, \(t_{\mathrm{w}}=\frac{v_0}{g}\), \(h_{\max}=\frac{v_0^2}{2g}\), a \(t_{\mathrm{lotu}}=\frac{2v_0}{g}\).

Rzut pionowy z wysokości

Jeżeli ciało nie wraca na poziom wyrzutu, nie można automatycznie stosować wzoru:

\[ t_{\mathrm{lotu}}=\frac{2v_0}{g}. \]

Trzeba wtedy skorzystać z pełnego równania położenia:

\[ y = y_0+v_{0y}t-\frac{gt^2}{2}. \]

Dla chwili uderzenia o ziemię, gdy poziom ziemi przyjęto jako \(y=0\):

\[ 0 = y_0+v_{0y}t-\frac{gt^2}{2}. \]

Jest to równanie kwadratowe względem czasu. Z jego rozwiązań wybieramy czas spełniający warunki fizyczne, najczęściej dodatni.

Prędkość końcową można następnie wyznaczyć z zależności:

\[ v_y = v_{0y}-gt, \]

albo bez obliczania czasu:

\[ v_y^2 = v_{0y}^2-2g(y-y_0). \]

Wykresy położenia, prędkości i przyspieszenia

Wykres położenia \(y(t)\)

W rzucie pionowym w górę:

\[ y(t) = y_0+v_0t-\frac{gt^2}{2}. \]

Wykres jest parabolą skierowaną ramionami w dół. Jego nachylenie w każdej chwili odpowiada prędkości:

\[ v_y(t) = \frac{dy}{dt}. \]

W najwyższym punkcie styczna do wykresu jest pozioma, ponieważ:

\[ v_y=0. \]

Wykres prędkości \(v_y(t)\)

Prędkość zmienia się liniowo:

\[ v_y(t) = v_0-gt. \] 

Nachylenie wykresu wynosi:

\[ \frac{dv_y}{dt} = -g. \] 

Pole algebraiczne pod wykresem prędkości jest równe przemieszczeniu:

\[ \Delta y = \int_{t_1}^{t_2}v_y(t)\,dt. \]

Wykres przyspieszenia \(a_y(t)\)

Przyspieszenie jest stałe:

\[ a_y(t)=-g. \]

Wykres jest poziomą linią znajdującą się poniżej osi czasu. Nie przecina osi nawet w najwyższym punkcie ruchu.

Wykresy położenia, prędkości i przyspieszenia w funkcji czasu dla rzutu pionowego w górę
Dla osi skierowanej w górę \(y(t)=y_0+v_0t-\frac{gt^2}{2}\), \(v_y(t)=v_0-gt\), a \(a_y(t)=-g\).

Symetria rzutu pionowego

Jeżeli ciało powraca na poziom wyrzutu i pomijamy opór powietrza, ruch jest symetryczny względem najwyższego punktu.

Jeżeli na pewnej wysokości podczas wznoszenia:

\[ v_y=+v_1, \]

to podczas spadania na tej samej wysokości:

\[ v_y=-v_1. \]

Szybkość w obu chwilach wynosi jednak:

\[ |\vec v|=v_1. \]

Przykłady obliczeniowe

Przykład 1. Swobodne spadanie z wysokości \(20\ \mathrm{m}\)

Ciało zostaje puszczone bez prędkości początkowej z wysokości:

\[ h=20\ \mathrm{m}. \]

Czas spadania wyznaczamy z równania:

\[ h=\frac{gt^2}{2}. \]

Stąd:

\[ t = \sqrt{\frac{2h}{g}} = \sqrt{\frac{2\cdot20}{9{,}81}} \approx 2{,}02\ \mathrm{s}. \]

Szybkość przed uderzeniem:

\[ |\vec v| = gt \approx 9{,}81\cdot2{,}02 \approx 19{,}8\ \mathrm{m/s}. \]

Przy osi skierowanej w górę algebraiczna prędkość wynosi:

\[ v_y\approx-19{,}8\ \mathrm{m/s}. \]

Ciało uderza w ziemię z szybkością około \(19{,}8\ \mathrm{m/s}\), poruszając się pionowo w dół.

Przykład 2. Rzut pionowy w górę

Piłkę wyrzucono pionowo w górę z prędkością:

\[ v_0=20\ \mathrm{m/s}. \]

Czas wznoszenia:

\[ t_{\mathrm{w}} = \frac{v_0}{g} = \frac{20}{9{,}81} \approx 2{,}04\ \mathrm{s}. \]

Maksymalna wysokość ponad poziom wyrzutu:

\[ h_{\max} = \frac{v_0^2}{2g} = \frac{20^2}{2\cdot9{,}81} \approx 20{,}4\ \mathrm{m}. \]

Całkowity czas lotu przy powrocie na poziom początkowy:

\[ t_{\mathrm{lotu}} = \frac{2v_0}{g} \approx 4{,}08\ \mathrm{s}. \]

Przykład 3. Rzut pionowy w dół

Ciało rzucono pionowo w dół z wysokości:

\[ h=45\ \mathrm{m} \]

z szybkością początkową:

\[ v_0=5\ \mathrm{m/s}. \]

Stosujemy równanie dodatniej drogi mierzonej w dół:

\[ h=v_0t+\frac{gt^2}{2}. \]

Po podstawieniu:

\[ 45 = 5t+\frac{9{,}81t^2}{2}. \]

Dodatnie rozwiązanie równania wynosi w przybliżeniu:

\[ t\approx2{,}56\ \mathrm{s}. \]

Szybkość końcowa:

\[ |\vec v| = v_0+gt \approx 5+9{,}81\cdot2{,}56 \approx 30{,}1\ \mathrm{m/s}. \]

Przy osi skierowanej w górę algebraiczna prędkość końcowa wynosi około:

\[ v_y\approx-30{,}1\ \mathrm{m/s}. \]

Przykład 4. Rzut w górę z dachu

Piłkę wyrzucono pionowo w górę z dachu znajdującego się na wysokości:

\[ y_0=15\ \mathrm{m} \]

z prędkością początkową:

\[ v_0=10\ \mathrm{m/s}. \]

Chwilę uderzenia o ziemię wyznaczamy z warunku \(y=0\):

\[ 0 = 15+10t-\frac{9{,}81t^2}{2}. \]

Po rozwiązaniu równania wybieramy dodatni pierwiastek:

\[ t\approx3{,}22\ \mathrm{s}. \]

Prędkość w chwili uderzenia:

\[ v_y = 10-9{,}81\cdot3{,}22 \approx -21{,}6\ \mathrm{m/s}. \]

Piłka uderza w ziemię z szybkością około \(21{,}6\ \mathrm{m/s}\), skierowaną pionowo w dół.

Jednostki i kontrola wymiarów

WielkośćOznaczenieJednostka SI
Wysokość i położenie\(h\), \(y\)\(\mathrm{m}\)
Czas\(t\)\(\mathrm{s}\)
Prędkość\(v_y\)\(\mathrm{m/s}\)
Szybkość\(|\vec v|\)\(\mathrm{m/s}\)
Przyspieszenie\(a_y\), \(g\)\(\mathrm{m/s^2}\)

Kontrola wymiarów równania położenia:

\[ [y] = [y_0]+[v_0t]+\left[\frac{gt^2}{2}\right], \] \[ \mathrm{m} = \mathrm{m} + \frac{\mathrm{m}}{\mathrm{s}}\mathrm{s} + \frac{\mathrm{m}}{\mathrm{s^2}}\mathrm{s^2}. \]

Wszystkie składniki mają jednostkę długości.

Kontrola wzoru na maksymalną wysokość:

\[ \left[\frac{v_0^2}{2g}\right] = \frac{\mathrm{m^2/s^2}}{\mathrm{m/s^2}} = \mathrm{m}. \]

Najczęstsze błędy

Po powrocie na poziom wyrzutu przemieszczenie jest równe zeru, ale droga nie. Jeżeli ciało osiągnęło maksymalną wysokość \(h_{\max}\) ponad poziom wyrzutu, całkowita droga wynosi:

\[ s=2h_{\max}. \]

Podsumowanie

Przy osi pionowej skierowanej dodatnio w górę przyspieszenie ziemskie ma postać:

\[ a_y=-g. \]

Ogólne równania ruchu pionowego są następujące:

\[ v_y=v_{0y}-gt, \] \[ y=y_0+v_{0y}t-\frac{gt^2}{2}, \] \[ v_y^2=v_{0y}^2-2g(y-y_0). \]

W swobodnym spadaniu z bezruchu:

\[ v_y=-gt, \qquad |\vec v|=gt, \] \[ h=\frac{gt^2}{2}. \]

W rzucie pionowym w dół, gdy \(v_0\) oznacza dodatnią szybkość początkową:

\[ v_y=-v_0-gt, \qquad |\vec v|=v_0+gt. \]

W rzucie pionowym w górę:

\[ v_y=v_0-gt, \] \[ y=y_0+v_0t-\frac{gt^2}{2}. \]

Czas wznoszenia i maksymalna wysokość ponad poziom wyrzutu wynoszą:

\[ t_{\mathrm{w}}=\frac{v_0}{g}, \qquad h_{\max}=\frac{v_0^2}{2g}. \]

Jeżeli ciało wraca na poziom wyrzutu:

\[ t_{\mathrm{lotu}}=\frac{2v_0}{g}. \]

W najwyższym punkcie prędkość jest chwilowo równa zeru, ale przyspieszenie nadal wynosi \(-g\). Podczas spadania algebraiczna prędkość może być ujemna, natomiast szybkość jest jej dodatnią wartością bezwzględną.