Oscylator Duffinga: Kiedy sprężyna staje się chaotyczna

Sprężyna, która nie chce być liniowa

Idealna sprężyna przestrzega prawa Hooke'a, czyli siły proporcjonalnej do przemieszczenia, a jej równanie ruchu to prosta funkcja sinusoidalna w nieskończoność. Równanie Duffinga z 1918 roku dodaje człon sztywności sześciennej, współczynnik tłumienia i okresowy napęd zewnętrzny, co prowadzi do pozornie prostego równania, które może generować wszystko od prostych drgań po pełny chaos, w zależności jedynie od amplitudy napędu.

ẍ + δẋ + αx + βx³ = γ·cos(ωt)

δ = damping,  α = linear stiffness,  β = cubic (nonlinear) stiffness
γ = drive amplitude,  ω = drive frequency

Twardnienie, zmiękczanie i złożona krzywa rezonansowa

Gdy β > 0 sprężyna staje się twardsza w miarę jej dalszego rozciągania (twardniejąca sprężyna); gdy β < 0 staje się zmiękczona (zmiękczająca). To samo powoduje odchylenie klasycznego szczytu rezonansowego oscylatora liniowego na boki, w kształcie krzywej pochylonej, a przy krytycznej amplitudzie napędowej ta krzywa zaginana jest i przesuwa się do siebie. W obszarze złożonym matematycznie możliwe są trzy różne amplitudy drgań przy tej samej częstotliwości – system staje się dwustany, skacząc nagle między odpowiedziami o małej i dużej amplitudzie w miarę powolnego przesuwania częstotliwości w górę i w dół, tworząc cykl pętlowy, który sprężyna liniowa nigdy nie byłaby w stanie wygenerować.

Od okresowego do chaotycznego w miarę wzrostu napędów

Gdy γ jest małe, oscylator ustabilizuje się w prostym ruchu okresowym o częstotliwości napędu. Zwiększając γ, odpowiedź może ulec podwójnemu okresowi – ustabilizować się cykl powtarzający się co drugie określenie napędu, następnie co czwarte, a następnie co osmé, co stanowi typową ścieżkę prowadzącą do chaosu w systemach nieliniowych, nie tylko w tym przypadku. Po punkcie gromadzenia się tej kaskady ruch staje się bezokresowy: ograniczony, deterministyczny, ale nigdy się nie powtarza, rysując dziwaczny atraktor w przestrzeni faz z tą samą strukturą rozciągania i składania co system Lorenz, mimo że pochodzi z zupełnie innej konfiguracji fizycznej.

Sekcja Poincaré: zamiana splątania w kropki

Rysowanie położenia względem prędkości bezpośrednio generuje gęsty, trudny do odczytania węzeł, gdy ruch jest chaotyczny. Standardową metodą jest sekcja Poincaré: zamiast rejestrować każdy moment, próbuj stan raz na cykl jazdy (stroboscopically, w tym samym stadium cos(ωt) za każdym razem). Okrągły okres ulega skurczeniu do pojedynczej kropki lub niewielkiej liczby punktów; okrąg podwojony pod względem okresu pokazuje dwa, cztery lub osiem punktów, a chaotyczny ruch ujawnia się jako skomplikowany, samopodobny rozproszenie punktów – fraktalna przekrój ukrytego dziwnego atraktora, widoczny z łatwością zamiast ukryty w skomplikowanej trajektorii.

Dlaczego równanie z 1918 roku nadal ma znaczenie

Równanie Duffinga jest jednym z najprostszych modeli, które zawiera twardość/miękkość w funkcji obciążenia, wymuszone rezonanse, histerezę oraz pełną ścieżkę do chaosu w jednej linii. Dlatego pozostaje standardowym modelem podręcznikowym i testowym dla inżynierów w przypadku drgań nieliniowych – od rezonatorów MEMS i elementów konstrukcyjnych, które uginają się pod obciążeniem, po analogiczne układy elektroniczne historycznie wykorzystywane do badania chaosu przed tym, gdy tanie komputery uczyniły integrację numeryczną łatwą.

Frequently asked questions

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Duffing Oscillator i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.