Spotlight #61 – Fizyka i inżynieria — Archimedes, wyboczenie, regulacja PID i wahadło 3D

I. Prawo Archimedesa — wyporność i stateczność statku

Odkrycie Archimedesa — że na zanurzone ciało działa siła wyporu równa ciężarowi wypartego płynu — to jeden z najstarszych ilościowych wyników w fizyce, wciąż codziennie wykorzystywany przez konstruktorów okrętów. Skierowana ku górze siła wyporu wynosi:

F_w  =  ρ_płynu · g · V_zanurzonej

Dla ciała pływającego w stanie równowagi statycznej F_w dokładnie równoważy ciężar W = m·g. Interesująca fizyka zaczyna się, gdy jednostka się przechyla. Gdy statek przechyla się o kąt θ, objętość podwodna przesuwa się na jedną stronę: środek wyporu (centroid wypartej objętości) przemieszcza się na zewnątrz. Metacentrum M to punkt, w którym pionowa linia przechodząca przez przesunięty środek wyporu przecina oś symetrii kadłuba. Jeśli M leży powyżej środka ciężkości G, moment przywracający wyprostowuje statek; jeśli M leży poniżej G, moment wywracający powoduje wywrotkę.

Wysokość metacentryczna i moment przywracający

Odległość pionowa GM (wysokość metacentryczna) w pełni charakteryzuje stateczność dla małych kątów. Moment przywracający wynosi:

M_przyw  =  W · GM · sin(θ)

Duża dodatnia GM (statek „sztywny") szybko się wyprostowuje, ale powoduje gwałtowną kołysankę, niekomfortową dla pasażerów i niebezpieczną dla ładunku. Mała GM (statek „miękki") kołysze się ospale i może nie powrócić do równowagi po dużych falach. Liniowce pasażerskie projektuje się w wąskich tolerancjach; w pełni załadowany kontenerowiec ma zupełnie inną krzywą stateczności niż ta sama jednostka w balaście. IMO wymaga, by statki posiadały Broszurę Informacji o Stateczności, opisującą, jak wysokość metacentryczna zmienia się w zależności od warunków załadunku — bezpośrednie zastosowanie prawa Archimedesa, obliczane dla każdego rejsu.

Symulacja numerycznie oblicza wypartą objętość przy każdym kącie przechyłu, całkując geometrię kadłuba poniżej linii wodnej. Dzięki temu uchwycone jest nieliniowe zachowanie, którego nie oddaje liniowa teoria metacentryczna: krzywa GZ (ramię prostujące w funkcji kąta) może mieć maksimum, spaść do zera (kąt zaniku stateczności), a nawet stać się ujemna przy skrajnych kątach, gdy zanurzenie krawędzi pokładu zalewa kadłub.

II. Wyboczenie pręta — obciążenie krytyczne Eulera

Smukły pręt poddany ściskaniu osiowemu nie ulega zniszczeniu przez bezpośrednie zmiażdżenie. Zamiast tego, przy obciążeniu krytycznym nagle ugina się na bok — wybacza się — tracąc niemal natychmiast zdolność do przenoszenia obciążenia. Leonhard Euler wyprowadził wzór na obciążenie krytyczne w 1744 roku:

P_kr  =  π²EI / (KL)²

Tutaj E to moduł Younga, I to moment bezwładności przekroju poprzecznego, L to długość pręta, a K to współczynnik długości efektywnej, zależny od warunków podparcia końców:

• K = 1.0 — oba końce przegubowe (swobodny obrót, brak przemieszczenia)
• K = 0.5 — oba końce utwierdzone (zablokowany obrót)
• K = 0.7 — jeden koniec utwierdzony, jeden przegubowy
• K = 2.0 — jeden koniec utwierdzony, jeden swobodny (maszt flagowy)

Wzór od razu pokazuje, dlaczego wysokie, smukłe pręty są niebezpieczne: P_kr maleje jak L². Podwojenie długości zmniejsza obciążenie krytyczne czterokrotnie. Odpowiednim parametrem geometrycznym jest smukłość KL/r, gdzie r = √(I/A) to promień bezwładności. Dla KL/r poniżej mniej więcej 120 (dla stali) dominuje wyboczenie niesprężyste: materiał zaczyna płynąć, zanim osiągnięte zostanie teoretyczne sprężyste obciążenie wyboczeniowe.

Parabola Johnsona i przejście do uplastycznienia

Parabola Johnsona przybliża przejście między uplastycznieniem a wyboczeniem sprężystym. Dla smukłości poniżej wartości krytycznej C_c = √(2π²E/F_y) naprężenie krytyczne wynosi:

F_kr  =  F_y · [ 1  -  (KL/r)² / (2·C_c²) ]

Powyżej C_c obowiązuje sprężysty wzór Eulera. Ten odcinkowy model stanowi podstawę projektowania stalowych słupów według specyfikacji AISC. Symulacja koloruje pręt według stanu naprężenia — sprężysty (niebieski), niesprężysty (żółty), uplastyczniony (czerwony) — dzięki czemu można obserwować, jak zmienia się mechanizm zniszczenia w miarę zmiany smukłości.

III. Regulator PID — algorytm, który rządzi światem

Regulacja proporcjonalno-całkująco-różniczkująca (PID) to najpowszechniej stosowany algorytm w automatyce przemysłowej. Szacuje się, że ponad 90% pętli regulacji w przemyśle procesowym wykorzystuje PID. Sygnał sterujący ma postać:

u(t)  =  K_p · e(t)  +  K_i · ∫e(τ)dτ  +  K_d · de(t)/dt

gdzie e(t) = r(t) - y(t) to uchyb między wartością zadaną r a wielkością zmierzoną y. Człon proporcjonalny reaguje na bieżący uchyb. Człon całkujący eliminuje uchyb ustalony, sumując uchyb z przeszłości. Człon różniczkujący przewiduje przyszły uchyb, reagując na tempo jego zmian — działa jak hamulec, tłumiąc przeregulowanie.

Strojenie metodą Zieglera-Nicholsa i zabezpieczenie przed przesyceniem całki

Klasyczna metoda Zieglera-Nicholsa stroi nastawy PID eksperymentalnie: przy aktywnym wyłącznie członie proporcjonalnym zwiększaj K_p, aż wyjście układu zacznie oscylować w sposób ciągły, o stałej amplitudzie. Zapisz to wzmocnienie graniczne K_u oraz okres graniczny T_u. Następnie ustaw:

K_p  =  0.6 · K_u
T_i  =  0.5 · T_u       (czyli K_i = K_p / T_i)
T_d  =  0.125 · T_u     (czyli K_d = K_p · T_d)

Daje to odpowiedź układu zamkniętego z przeregulowaniem rzędu 25% i szybkim czasem ustalania. Dla bardziej konserwatywnego strojenia stosuje się metody Cohena-Coona lub IMC. Praktycznym problemem, który zawsze się pojawia, jest przesycenie całki (windup): gdy obiekt się nasyca (zawór nie może otworzyć się bardziej niż na 100%), człon całkujący dalej sumuje uchyb, narastając do dużej wartości, która następnie powoduje długotrwałe przeregulowanie, gdy obiekt wychodzi z nasycenia. Strategie zabezpieczające przed przesyceniem obejmują całkowanie warunkowe (zatrzymanie sumowania, gdy wyjście jest nasycone) oraz metodę odwrotnego przeliczania (odjęcie nasycenia od wartości całki).

Symulacja łączy regulator PID z obiektem drugiego rzędu, modelowanym jako układ masa-sprężyna-tłumik. Przeciągając suwaki wzmocnień, można bezpośrednio obserwować przejście między odpowiedzią nadkrytycznie tłumioną (wolna, bez przeregulowania), podkrytycznie tłumioną (oscylacyjna) i krytycznie tłumioną (najszybsze ustalenie się bez przeregulowania).

IV. Wahadło 3D — tryby Foucaulta, stożkowy i rozetkowy

Wahadło sferyczne ma dwa stopnie swobody — kąt biegunowy θ od pionu i kąt azymutalny φ wokół pionu — powiązane więzem, że masa punktowa pozostaje na sferze o promieniu L. Sformułowanie Lagrange'a, śledzące różnicę energii kinetycznej i potencjalnej, daje dwa sprzężone równania różniczkowe drugiego rzędu. Dla małych amplitud ruch rozkłada się na dwa niezależne oscylatory harmoniczne o tej samej częstości, co oznacza, że dowolny warunek początkowy daje w płaszczyźnie poziomej figurę przypominającą krzywą Lissajous: linię prostą (wahadło płaskie), okrąg (wahadło stożkowe) lub elipsę (przypadek ogólny).

Precesja Foucaulta i obrót Ziemi

Gdy obrót Ziemi zostanie uwzględniony jako człon Coriolisa, pozorna płaszczyzna wahań wahadła płaskiego powoli precesuje. Szybkość precesji zależy od szerokości geograficznej φ:

ω_F  =  Ω · sin(φ)

gdzie Ω = 7.27 × 10⁻⁵ rad/s to prędkość kątowa obrotu Ziemi. Na biegunie północnym (φ = 90°) płaszczyzna wahań wykonuje pełny obrót w ciągu jednej doby gwiazdowej (23 godz. 56 min). Na szerokości Londynu (51,5°N) szybkość precesji wynosi Ω · sin(51,5°) ≈ 5,68 × 10⁻⁵ rad/s, co daje okres pełnego obrotu wynoszący około 30,5 godziny. Jean-Bernard-Léon Foucault zawiesił 67-metrowe wahadło pod kopułą paryskiego Panteonu w 1851 roku i zaprosił publiczność do obserwowania, jak piasek kreśli na podłodze obracający się wzór — pierwszy bezpośredni, wizualny dowód na to, że Ziemia się obraca.

Symulacja kompresuje w czasie obrót Ziemi nawet do 1000×, dzięki czemu precesja staje się widoczna w ciągu kilku sekund. Na średnich szerokościach geograficznych tor kreśli charakterystyczny wzór rozetkowy: kolejne wahnięcia są przesunięte o kąt Foucaulta, tworząc gwiazdopodobną figurę, która stopniowo wypełnia całą płaszczyznę.

V. Soczewka akustyczna — ogniskowanie fal metodą FDTD

Metoda różnic skończonych w dziedzinie czasu (FDTD) rozwiązuje równanie falowe poprzez dyskretyzację zarówno przestrzeni, jak i czasu na siatce. Dla dwuwymiarowego pola ciśnienia akustycznego p(x,y,t) w ośrodku o lokalnej prędkości fali c(x,y):

∂²p/∂t²  =  c(x,y)² · ∇²p

Schemat FDTD wykorzystuje dyskretyzację typu „żabi skok" (leapfrog): ciśnienie i prędkość są przechowywane na przesuniętych siatkach, oddalonych o pół kroku czasowego. Ten schemat o dokładności drugiego rzędu jest stabilny, gdy spełniony jest warunek Couranta-Friedrichsa-Lewy'ego (CFL): c·Δt/Δx ≤ 1/√2 w 2D. Naruszenie warunku CFL powoduje niestabilność numeryczną — symulacja rozbiega się wykładniczo.

Absorbujące warunki brzegowe Mura i projektowanie soczewki

Warunki brzegowe mają kluczowe znaczenie w FDTD. Odbijające (twarde) ściany po prostu wymuszają ∂p/∂n = 0 na brzegu. Jednak odbicia od granicy obszaru obliczeniowego zniekształcają symulacje promieniowania. Absorbujący warunek brzegowy Mura (1981) przybliża falę wychodzącą jako równanie jednokierunkowe i stosuje je na krawędzi siatki, absorbując 95–99% padającej energii przy minimalnym koszcie obliczeniowym — nowocześniejsza alternatywa, doskonale dopasowana warstwa (PML), osiąga niemal idealną absorpcję dzięki wprowadzeniu gradientu dyssypatywnego w komórkach brzegowych.

Akustyczna soczewka z gradientem współczynnika załamania działa poprzez zmianę c(x,y) na całej aperturze: wolniejsza prędkość fali w centrum, szybsza na krawędziach (lub odwrotnie, zależnie od typu soczewki). Wolniej propagujące się promienie centralne pozostają w tyle za promieniami zewnętrznymi, wyginając front falowy w kształt zbieżny. Ogniskowa podąża za akustycznym odpowiednikiem równania soczewkowego. Płytki strefowe Fresnela osiągają podobne ogniskowanie poprzez selektywne blokowanie, a nie załamanie: naprzemienne pierścieniowe strefy otwarte i zablokowane powodują interferencję konstruktywną w punkcie ogniskowym.