Przepływ ruchu na autostradzie: traktowanie drogi jako…

Samochody jako kontynuum

Model Lighthill-Whitham-Richards (LWR), opublikowany niezależnie przez Lighthilla i Whithama oraz Richardsa w latach 1955-56, dokonuje odwagomej uproszczenia: zamiast śledzić każdy samochód osobno, traktuje ruch pojazdów jako ciągłe, spompowalne ciecz opisanej jednym polem — gęstości pojazdów k(x, t) – liczbę pojazdów na jednostkę długości drogi w punkcie x i czasie t. To zignoruje wszystkie informacje dotyczące indywidualnych kierowców i zmiany torów jazdy, ale pozwoliło to złapać duże skala zjawisk, które naprawdę są ważne dla komunikanta: jak formują się zatłoczenia, jak się rozprzestrzeniają i jak szybko porusza się tył tłumu w stosunku do drogi.

demo na żywo · powiązana symulacja● LIVE

Zachowanie pojazdów i podstawowy diagram

Model opiera się na jednym dokładnym prawie fizycznym — pojeźdzy nie są tworzone ani zniszczone na przeciętnej trakcie drogowej — oraz jednej empirycznej relacji zamknięcia:

Δk/δt + Δq/δx = 0 zachowanie pojazdów (zakonczona równość kontynuitetu) q(k) = k · v(k) przepływ = gęstość × prędkość (podstawowy diagram) v(k) = v_max · (1 - k/k_jam) najprostsza relacja zamknięcia: prędkość maleje liniowo z gęstością Podstawowy diagram q(k) jest empirycznym sercem całościowej teorii: przepływ zaczyna się od zera (brak pojazdów, brak przepływu), wzrasta wraz z zwiększeniem gęstości i rozpoczęciem ruchu pojazdów, osiąga maksimum na pewnej krytycznej gęstości, gdzie droga nosi swoje możliwe maksymalne przepustowość, a następnie maleje do zera przy zbliżaniu się gęstości do gęstości zatłoczenia k_jam, w której pojazdy są zapakowane bokiem do boku i nie mogą poruszać się wcale. Podstawienie relacji zamknięcia do zachowania pojazdów przekształca model w jedno równanie różniczkowe cząstkowe pierwszego rzędu nieliniowego znane jako równanie LWR, należące do tej samej rodziny matematycznej (scalarne zasady zachowania hyperboliczne) jak równania opisujące falovate ciepła w dynamiczności gazowej.

∂k/∂t + ∂q/∂x = 0        conservation of vehicles (a continuity equation)

q(k) = k · v(k)            flow = density × speed  (the fundamental diagram)

v(k) = v_max · (1 - k/k_jam)     simplest closure: speed falls linearly with density

Wolne falowe kinematyczne: dlaczego zaburzenia podróżują w tył

Dana wartość gęstości propaguje się po drodze z prędkością fali kinematycznej c(k) = dq/dk, a nie z prędkością ruchu samochodów v(k) — to kluczowa rozbieżność. W małych przepływnościach c(k) jest blisko v_max i jest dodatnia: informacje podróżują w przód, prawie razem z samochodami. Jednak po przejściu gęstości krytycznej, na zatłoczonym odcinku diagramu podstawowego, dq/dk staje się ujemne. Małe zwiększenie gęstości tam powoduje spadek prądu ruchu, a zaburzenie to propaguje się w tył względem drogi — dokładnie takie doświadczenie jest znane z lustrzanych świateł hamulcowych rozprzestrzeniających się przed wami szybciej niż jadę, daleko przed dotarciem do rzeczywistego powodu spowolnienia.

Czujki w punktach zatłoczeń

Dzięki temu, że równanie LWR jest nieliniowe, krzywe charakterystyczne przesyłające różne wartości gęstości mogą się łączyć i przecinać — matematyczna konieczność, gdy szybko ruch light traffic (RLT) napotyka zatłoczony, gęstszy ruch przed nim, dokładnie to, co dzieje się przy zamknięciu kanału, łączeniu z drogi lub zatrutym pojazdzie. Gdy krzywe charakterystyczne przecinają się, spokojna rozwiązanie przestaje być liniowe i musi formować się rzeczywista nieciągłość — czujkę. Szybkość propagacji tej czujki jest zawsze ustalonej niezależnie od lokalnych warunków PDE (równania różniczkowego parabolicznego), ale z argumentu zachowania globalnego, warunku Rankine-Hugoniot:

s = (q_right - q_left) / (k_right - k_left) s = szybkość propagacji czujki (granicznej zatłoczenia) subskrypcje „lewo”/„prawo” = stan wolnego ruchu i zatłoczonych pojazdów po i przed frontem Ujemna wartość s oznacza, że czujka — tył zatłoczenia — porusza się w kierunku przepływu (w górę), dokładnie to, co obserwują kierowcy, gdy zatłoczenie wydaje się rosnąć w tył, nawet gdy tył zatłoczenia, gdzie samochody ostatecznie przebywają przez punkt zatłoczenia, może pozostawać stały lub przemieszczać się wolno w kierunku przodu. Gdzie inaczej wymagałyby bycia krzywe charakterystyczne divergentne — przechodząc od gęstej do rzadkiej sytuacji, jak przy zmianie z czerwonego światła na zielone — poprawna fizyczna rozwiązanie to zamiast tego spokojny fani rarefakcji rozszerzający przejście ciągle w przestrzeni, dlatego samochody zwalniają się wolno z czerwonego światła, a nie opuszczają jak nagle stwarzać ścianę.

s = (q_right - q_left) / (k_right - k_left)

s = shock (jam-boundary) propagation speed
subscripts "left"/"right" = the free-flow and congested states either side of the front

Co model prawidłowo przewiduje, a co go o nim omija

LWR poprawnie przewiduje kwantalną istnienie i odwrotną propagację fal szokowych, położenie oraz zachowanie fal szokowych w punktach zatłoczeń, a także rolę gęstości krytycznej jako ograniczającego wydajność — ciekawe, że to wszystko się zdarza, mimo że model nie zawiera żadnych informacji na temat reakcji pojazdów pojedynczych, odległości między pojazdami ani zachowania zmiany kierunku jednolitej drogi. To, co omija, to dokładnie te szczegółowe mikroskopowe aspekty: nie może powtórzyć samonapływających „przeraźliwych” zatorów, które wydzielają się ze zdeterminowanej reakcji kierowców na jednolitej drodze bez żadnego punktu zatłoczenia. Te zatory wymagają modeli odchętnicowych (car-following) lub automatonów komórkowych (takich jak Nagel-Schreckenberg), które mają jasno określony opóźnienie reakcji, aby je zachować. W praktyce inżynierowie ruchu drogowego używają LWR do planowania makroskopicznych koridorów i autostrad, a przejście do mikroskopowych lub automatycznych modeli następuje wtedy, gdy pytanie dotyczy samonapływających zatorów wynikających bezpośrednio z zachowania pojedynczych kierowców.

Często zadawane pytania

Dlaczego szarża ruchu porusza się w tył, nawet jeśli każda kasa porusza się w przód?

Szarża jest wzorcem gęstości, a nie zestawem samochodów — pojedyncze kasy wpadają do tyłu szarży, przebywają przez nią i wyjazd z jej początku, ale granica zagęszczonego obszaru propaguje się w górę, ponieważ każdy kierowca reaguje na samochód przed nim ze spóźnieniem. Model LWR pokrywa to dokładnie: prędkość fali dq/dk jest ujemna w zagęszczonym ramie diagramu podstawowego, więc zaburzenie porusza się w tył względem drogi, nawet gdy prędkość każdego pojazdu pozostaje dodatnia.

Co to jest diagram podstawowy przepływu ruchu?

To empiryczna zależność q = k · v(k) między przepływem ruchu (pojazdów na godzinę) a gęstością (pojazdów na kilometr). Przepływ wzrasta z zera przy zera gęstości, osiąga maksimum w krytycznej gęstości, gdzie droga nosi maksymalny przepustowość, a następnie spada do zera przy gęstości szarży, gdzie samochody są bujno i nie mogą poruszać się. Model LWR używa tej jednej krzywej, razem z zachowaniem ilości pojazdów, jako swojej całościowej założenia.

Co wywołuje fale szoku w modelu LWR?

Każdy lokalny skok gęstości, gdzie szybko poruszające się ruch napotyka wolniejsze ruch przed nim — zamknięcie kanału, łączenie rampy, wypadek lub prosto nagle zatrzymanie. Ponieważ równanie LWR jest nieliniowym prawem zachowania, cechy przewożące różne wartości gęstości mogą się zbiegać i przecinać; gdy to się dzieje, spokojne rozwiązanie zanika i formuje się niespojna fala szoku, która propaguje się z prędkością ustalona przez warunek skaczący Rankine-Hugoniot.

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz LWR Traffic Flow: Treating the Motorway as a Fluid i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.

▶ Otwórz symulację LWR Traffic Flow: Treating the Motorway as a Fluid

Przepływ ruchu na autostradzie: traktowanie drogi jako cieczy