Pressure-Impulse i model SDOF: dlaczego konstrukcje nie reagują tylko na pik ciśnienia

Streszczenie

Konstrukcja obciążona falą podmuchową nie pyta wyłącznie o maksymalne nadciśnienie. Pyta także, jak długo to nadciśnienie działało, jaki impuls przekazano, jaka była masa i sztywność elementu, jaki jest jego okres własny oraz ile pracy plastycznej może przyjąć bez utraty nośności. Dlatego wykres pressure-impulse i model SDOF są naturalnym następnym krokiem po modelach Kingery-Bulmash, Brode i CONWEP.^1,2

Ten artykuł wyjaśnia, dlaczego dwa obciążenia o tym samym piku ciśnienia mogą dawać różne skutki, jeżeli różnią się impulsem i czasem dodatniej fazy. Pokazuje też, co robi lokalny kalkulator Pressure-Impulse: bierze falę odbitą, zastępuje dodatnią fazę trójkątem o tym samym piku i impulsie, a następnie liczy odpowiedź ściany jako jednej masy, sprężyny i tłumika na jednostkę powierzchni. To jest model dydaktyczny, nie projekt rzeczywistej ściany.

Rozszerzenie tematu

Najłatwiej myśleć o wybuchu przez jedną liczbę: Pmax. Jeżeli nadciśnienie jest małe, szyby pękają rzadziej; jeżeli jest większe, uszkodzenia są poważniejsze. Taka skala ma sens do mapy wstępnej, ale konstrukcja nie jest wskaźnikiem manometru. Ma masę, sztywność, połączenia, tłumienie i własny czas reakcji. Dlatego odpowiedź ściany na obciążenie krótkie i ostre może być inna niż na obciążenie niższe, ale dłuższe.

Właśnie tu pojawia się przestrzeń pressure-impulse, zwykle oznaczana jako P-I. Na jednej osi odkłada się ciśnienie szczytowe (P lub Ps), a na drugiej impuls dodatniej fazy, czyli całkę z nadciśnienia po czasie. Dla danego elementu konstrukcyjnego można narysować krzywe podobnego uszkodzenia: poniżej krzywej element pozostaje w założonym stanie, a powyżej krzywej wymaga większej pracy plastycznej albo ulega zniszczeniu. Jedna liczba Pmax przestaje wystarczać.¹

Ten język dobrze pasuje do całej ścieżki kalkulatorów serwisu. Podmuch liczy parametry fali. Height-of-burst i Mach stem pokazuje, jak geometria może zmienić obciążenie przy gruncie. Pressure-Impulse przenosi falę na przegrodę. Odporność płyty idzie krok dalej i pyta o geometrię oraz materiał konkretnego elementu.

Historia diagramów P-I: od projektowania schronów do antyterroryzmu

Diagramy pressure-impulse rozwinęły się w latach 50.–60. jako narzędzie projektowania schronów jądrowych. USACE (US Army Corps of Engineers) i BRL (Ballistic Research Laboratory) pracowały nad praktyczną metodą oceny odporności budowli na blast bez pełnej analizy dynamicznej.

Kluczowe dokumenty historyczne:

• TM 5-1300 (1969, 1990 revision) — "Structures to Resist the Effects of Accidental Explosions". Pierwszy oficjalny US Army podręcznik projektowania obiektów odpornych na blast. Zawiera krzywe P-I dla różnych typów elementów konstrukcyjnych.
• UFC 3-340-02 (2008) — "Structures to Resist the Effects of Accidental Explosions". Aktualizacja TM5-1300, szeroko stosowana w US DoD. Zawiera mapy P-I, wytyczne projektowania, tabele duktylności, odpowiedź płyt betonowych i stalowych.
• ASCE/SEI 59-11 (2011) — "Blast Protection of Buildings". Norma ASCE dla cywilnych budynków, opisuje procedury analizy i projektowania odporności na blast.
• ISO 16933 i ISO 16934 — Europejskie normy dla szyb odpornych na blast.

W kontekście broni jądrowej P-I diagrams pojawiają się w klasycznym "Effects of Nuclear Weapons" (Glasstone i Dolan 1977) — nie pod tą nazwą, ale jako krzywe prawdopodobieństwa uszkodzenia dla różnych typów budynków w zależności od nadciśnienia. Metoda P-I eksplicytnie jest nowsza, ale intuicja pochodzi właśnie z analiz jądrowych.

Po zamachu na WTC (1993) i Oklahoma City (1995) P-I diagramy znalazły szerokie zastosowanie w ochronie budynków rządowych przed terrorystami (VBIED — vehicle borne IED). Program GSA (General Services Administration) w USA i program CPNI w UK wydały wytyczne do oceny ryzyka blast dla budynków publicznych — wszystkie bazujące na P-I framework.¹

Pik, impuls i czas dodatniej fazy

Impuls dodatniej fazy można zapisać intuicyjnie jako pole pod wykresem nadciśnienia:

I = ∫ P(t) dt.

Jeżeli przyjmiemy prosty trójkątny przebieg dodatniej fazy, to:

I = 0,5 * Pmax * td,

czyli:

td = 2I / Pmax.

Ta prosta zależność jest bardzo pouczająca. Dwa obciążenia mogą mieć ten sam Pmax, ale różny I, a więc różny czas td. Mogą też mieć ten sam impuls, ale różny pik. Dla konstrukcji te przypadki nie są równoważne, bo element może zdążyć ruszyć, wejść w drgania albo praktycznie tylko otrzymać krótki impuls prędkości.

W tomie Safety of structures under terrorist threat rozdział o obciążeniach wybuchem rozdziela metody empiryczne i numeryczne, opisuje skalowanie pierwiastkiem sześciennym, falę odbitą, ciśnienie dynamiczne i generowanie historii ciśnienia w czasie przez narzędzia takie jak CONWEP. To ważne tło: zanim policzymy odpowiedź ściany, musimy wiedzieć, czy wejściem jest fala padająca, odbita, powierzchniowa, powietrzna, czy przebieg złożony przez odbicia w pomieszczeniu.¹

Trzy reżimy odpowiedzi konstrukcji

W przestrzeni P-I można wyróżnić trzy intuicyjne reżimy.

Pierwszy jest impulsowy. Obciążenie trwa znacznie krócej niż okres własny elementu. Konstrukcja nie zdąża się w pełni ugiąć podczas działania ciśnienia, ale otrzymuje prędkość początkową. W tym reżimie ważny jest impuls, a sam pik może być mniej rozstrzygający, jeżeli trwa bardzo krótko.

Drugi jest quasi-statyczny. Obciążenie trwa długo w porównaniu z okresem własnym. Element ma czas podążać za obciążeniem, więc bliżej mu do problemu siły statycznej. Wtedy bardziej rozstrzyga wartość ciśnienia niż sama krótkotrwała całka impulsu.

Trzeci jest dynamiczny pośredni. Czas dodatniej fazy jest porównywalny z okresem własnym elementu. To najtrudniejszy obszar, bo odpowiedź zależy od fazy drgań, tłumienia, kształtu obciążenia i nieliniowości materiału. Właśnie w tym środku wykresu P-I krzywe uszkodzenia są najbardziej informacyjne.

Model SDOF

SDOF oznacza single degree of freedom, czyli model jednego stopnia swobody. Ściana, płyta albo panel zostaje sprowadzony do zastępczej masy, sprężyny i tłumika. Zamiast analizować każdy punkt konstrukcji, liczymy jedno uogólnione przemieszczenie y(t).

Najprostsza postać równania wygląda tak:

m y'' + c y' + k y = p(t).

m to masa na jednostkę powierzchni, c tłumienie, k sztywność na jednostkę powierzchni, a p(t) ciśnienie działające na element. Jeżeli znamy okres własny Tn, to często wygodnie zapisać:

omega = 2π / Tn,

k = m omega^2,

c = 2 zeta m omega,

gdzie zeta jest ułamkiem tłumienia krytycznego.

Taki model jest brutalnym uproszczeniem, ale bardzo przydatnym. Pozwala zobaczyć, czy problem jest raczej impulsowy, quasi-statyczny czy dynamiczny. Pozwala też wyjaśnić, dlaczego cięższa ściana nie zawsze jest automatycznie bezpieczniejsza, jeżeli ma długi okres własny i małą sztywność, oraz dlaczego sztywny lekki element może mieć duże przyspieszenia.

Co robi lokalny kalkulator

Lokalny kalkulator Pressure-Impulse przyjmuje efektywną masę TNT, odległość, powierzchnię ściany, masę powierzchniową, okres własny, tłumienie i umowne ugięcie graniczne. Następnie używa BlastPhysics::compute() w trybie Kingery-Bulmash dla wybuchu powierzchniowego i bierze przede wszystkim falę odbitą: reflected_overpressure oraz reflected_impulse_kPa_ms.²

Potem kalkulator przelicza:

• ciśnienie szczytowe na siłę szczytową F = P * A,
• impuls ciśnienia na impuls całkowity J = I * A,
• impuls i pik na czas dodatniej fazy trójkąta,
• masę, okres i tłumienie na równanie SDOF,
• ugięcie maksymalne, prędkość, dynamic load factor i żądanie duktylności.

Solver w kodzie używa wariantu Newmark-beta ze średnim przyspieszeniem (beta = 0,25, gamma = 0,5). Obciążenie jest trójkątne i ma ten sam pik oraz impuls co wejściowa fala dodatnia. To dobry kompromis dydaktyczny: wynik zależy od P, I, Tn, masy i tłumienia, ale nie wymaga pełnego modelu materiałowego.²

Żądanie duktylności

Kalkulator raportuje także mu, czyli żądanie duktylności:

mu = ymax / yy.

ymax to maksymalne ugięcie dynamiczne, a yy to umowne ugięcie graniczne podane przez użytkownika. Jeżeli mu jest znacznie poniżej 1, model sugeruje odpowiedź sprężystą z zapasem względem tego progu. Jeżeli mu zbliża się do 1, element jest blisko granicy użytkowej. Jeżeli mu przekracza 1, prosty model liniowy zaczyna tylko sygnalizować, że potrzebna byłaby praca plastyczna.²

To nie jest formalne kryterium zniszczenia. W prawdziwej ścianie ważne są zbrojenie, pękanie, odspojenie warstw, połączenia, ścinanie, lokalne przebicie, podparcia i wiele trybów drgań. mu w kalkulatorze jest wskaźnikiem dydaktycznym, który mówi: „ten przypadek wymaga nieliniowej mechaniki, jeżeli mamy traktować wynik poważnie”.

Dynamic Load Factor

Dynamic load factor w kalkulatorze jest stosunkiem maksymalnego ugięcia dynamicznego do ugięcia statycznego od ciśnienia szczytowego:

DLF = ymax / ystat.

Jeżeli obciążenie narasta i zanika w czasie porównywalnym z okresem własnym elementu, DLF może być istotnie różny od 1. Dla bardzo krótkiego impulsu ugięcie statyczne od samego piku może być mylące, bo pik działa zbyt krótko, aby element osiągnął taką deformację. Dla długiego obciążenia wynik zbliża się do statycznej intuicji.

To właśnie odpowiada na główne pytanie artykułu: konstrukcje nie reagują tylko na pik, bo mają własny czas reakcji. Ciśnienie jest wymuszeniem, a konstrukcja jest układem dynamicznym.

Analityczne asymptoty diagramu P-I

Krzywa P-I dla prostego modelu SDOF z odpowiedzią plastyczną ma dwie analityczne asymptoty, które dobrze ilustrują fizykę.

Asymptota impulsowa (lewa strona diagramu, małe td/Tn):

Gdy obciążenie jest bardzo krótkie (td << Tn), element nie zdąży się ugić podczas obciążenia. Całe ciśnienie przekazuje prędkość początkową v₀ = I/m. Energia kinetyczna E_k = m v₀²/2 = I²/(2m) musi być pochłonięta przez pracę plastyczną. Warunek ugięcia granicznego y_y daje:

I_crit = sqrt(2 m k) * y_y (dla odpowiedzi sprężystej)

lub po uwzględnieniu duktylności μ: krzywa tylko przesuwa się w prawo. Asymptota pionowa w lewej części diagramu ma stałą wartość I_crit niezależną od ciśnienia.

Asymptota quasi-statyczna (dolna część diagramu, duże td/Tn):

Gdy obciążenie jest bardzo długie (td >> Tn), element zdąży dojść do równowagi statycznej pod ciśnieniem bieżącym. Warunek zniszczenia to po prostu P_crit = k * y_y. Asymptota pozioma ma stałą wartość P_crit niezależną od impulsu.

Obszar dynamiczny: Między asymptotami krzywa jest hiperbolą podobną do (P - P∞)(I - I∞) = const. W tym obszarze oba parametry mają znaczenie.

Asymptoty można znaleźć analitycznie dla SDOF i różnych kształtów obciążeń (trójkąt, prostokąt, Friedlander), co jest dobrym ćwiczeniem algebry dla studentów inżynierii. Kalkulator lokalny trafnie pokazuje, że punkt (P, I) leżący blisko jednej z asymptot oznacza, że cały wynik zależy tylko od odpowiadającego parametru.

Różne typy elementów: jak bardzo różnią się krzywe P-I

Krzywe P-I są silnie zależne od rodzaju elementu. Kilka charakterystycznych przykładów:

Szyby ze szkła float. Małe yy (pęknięcie przy ugięciu < 5 mm), mała masa (ok. 20–25 kg/m² dla 10 mm), mała sztywność (zależy od wymiarów tafli). Krzywa P-I leży nisko — szyby pękają już przy nadciśnieniu kilku kPa (0,5–1 psi) i impulsie kilku kPa·ms. Wzmocnione szyby (laminat PVB) mają wyższą krzywą, bo PVB pochłania energię po pęknięciu szkła.

Ściany murowane. Niezbrojone mury są kruche: małe ugięcie graniczne, brak duktylności. Krzywa P-I leży w obszarze dużych ciśnień i małych impulsów (element reaguje na pik, nie na impuls). Po przekroczeniu wytrzymałości mur kruszy się bez absorpcji dużej energii.

Płyty żelbetowe. Lepsza duktylność, możliwość pracy plastycznej. Krzywa P-I leży wyżej i bardziej w prawo — element może wytrzymać większy impuls przy umiarkowanym ciśnieniu. Odpowiedź zależy od zbrojenia, grubości, stopnia zbrojenia i sposobu podparcia.

Panele stalowe. Wysoka duktylność, zdolność do dużych odkształceń plastycznych. Krzywa P-I jest przesunięta w prawo — stalowe panele mogą pochłaniać duże energie impulsowe bez zniszczenia (np. w barierach ochronnych wokół budynków rządowych).

Porównanie tych typów na jednym wykresie P-I dobitnie pokazuje, że „to samo nadciśnienie" ma zupełnie inne skutki dla szyby, muru i stali. Właśnie dlatego progi nadciśnienia z map skutków wybuchu są tylko kategoryzacją grubą, a nie projektem.⁴

UFC 3-340-02 i standardy projektowania ochronnego

UFC 3-340-02 (dawniej TM5-1300) to główny standard projektowania budynków odpornych na blast w USA (US DoD). Dokument jest jawny i dostępny na stronie wbdg.org. Polscy inżynierowie i badacze korzystają z niego jako referencji de facto, mimo że Polska ma własne przepisy (Eurokod EN 1991-1-7).

Kluczowe elementy UFC dla modelu P-I:

1. Tabele właściwości dynamicznych (tzw. Equivalent SDOF parameters): dla płyt betonowych, belek stalowych, paneli murowanych UFC podaje przeliczniki do modelu SDOF — współczynnik masy transformacji KM i współczynnik obciążenia KL, które pozwalają przeliczyć rzeczywistą geometrię na zastępczą masę i sztywność modelu liniowego.

2. Krzywe P-I dla poszczególnych elementów: wykresy zniszczenia i stopnia uszkodzenia dla ścian żelbetowych, murowanych, betonowych bez zbrojenia, paneli stalowych, szyb.

3. Poziomy ochrony (Protection Level): UFC definiuje cztery poziomy — od „Very Low" (brak ochrony ludzi, możliwe zniszczenie budynku) do „Very High" (ochrona niezbędnych funkcji, bez zniszczeń strukturalnych). Każdy poziom odpowiada różnym wymaganiom mu (żądanie duktylności).

4. Ductility ratios: UFC dopuszcza mu od 1,0 (sprężysty) do 10–20 (pełna praca plastyczna, trójniki zbrojeniowe) zależnie od materiału i wymagań.

Przykład liczbowy SDOF z konkretną ścianą

Rozważmy ścianę żelbetową:

• masa powierzchniowa m = 480 kg/m² (20 cm beton, ρ=2400 kg/m³)
• okres własny Tn = 30 ms (mierzone lub z UFC tablic)
• tłumienie ζ = 5%
• ugięcie graniczne yy = 5 mm

Obciążenie z wybuchu (fala odbita, K-B):

• Pr = 150 kPa, Ir = 600 kPa·ms
• Czas dodatniej fazy: td = 2 × 600/150 = 8 ms

Stosunek td/Tn = 8/30 ≈ 0,27 — reżim pośredni, bliski impulsowemu.

Obliczenie SDOF (Newmark-beta, krok dt = 0,1 ms):

• Sztywność: k = m(2π/Tn)² = 480 × (209,4)² ≈ 21,0 MN/m²
• Odpowiedź maksymalna: ymax ≈ 4,3 mm
• DLF: ymax/ystat = 4,3/(150000/21000000) ≈ 0,60
• Żądanie duktylności: mu = 4,3/5,0 = 0,86

Wynik: ściana pozostaje w obszarze sprężystym z zapasem (mu < 1). Gdyby Pr = 300 kPa (podwojenie), wynik ymax ≈ 8,6 mm, mu = 1,72 — ściana wchodzi w pracę plastyczną.

Taki przykład dydaktyczny, przeprowadzony w arkuszu kalkulacyjnym, trwale pokazuje różnicę między statycznym i dynamicznym spojrzeniem na obciążenie.²

Multi-DOF i rozszerzenia modelu

Model SDOF jest pierwszym krokiem. Budynek ma wiele elementów i wiele stopni swobody. Kilka ważnych rozszerzeń:

2-DOF: ściana + strop. Jeżeli ściana i strop są połączone sprężyście, uderzenie w ścianę przenosi siłę na strop. Model 2-DOF pozwala sprawdzić, czy strop nie jest zagrożony przez dynamiczne wzmocnienie od ściany.

Multi-mode response. Płyta dwuprzęsłowa ma kilka modów drgań. Dominuje zwykle mod pierwszy, ale przy blast niedaleko środka rozpiętości mody wyższe mogą być ważne. Odpowiedź sumuje się przez superpozycję modalną.

Geometryczna nieliniowość. Przy dużych ugięciach płyta zaczyna pracować jako membrana (siły podłużne, nie tylko zginanie). Nośność membrany jest wyższa niż nośność zginania, więc płyta może wytrzymać więcej niż sugeruje liniowy SDOF.

Ruch gruntu. Przy wybuchach bliskich powierzchni gruntu trzęsienie gruntu (ground shock) dochodzi do fundamentów i może wzbudzać budynek przez inny kanał niż fala blast. Modele SDOF działające tylko na blast nie widzą tego zjawiska.

W tym serwisie model SDOF jest wystarczający dydaktycznie do wyjaśnienia idei P-I. Pełne modele konstrukcyjne należą do specjalistycznych narzędzi inżynierskich.¹

Polska norma i regulacje blast protection

W Polsce projektowanie ochrony budynków przed wybuchem opiera się na:

PN-EN 1991-1-7 — Eurokod 1, część 1-7: Oddziaływania na budowle. Oddziaływania wyjątkowe. Opisuje wybuchy chemiczne jako zdarzenia wyjątkowe i daje ogólne zasady obliczania obciążeń blast, ale bez szczegółowych krzywych P-I dla poszczególnych elementów.

Przepisy MON i BBN — Obiekty wojskowe i infrastruktura krytyczna mają własne przepisy techniczne MON, bazujące na standardach NATO (STANAG, ATP) i normach US DoD adaptowanych do warunków polskich.

Instrukcja WAT — Wojskowa Akademia Techniczna wydaje wewnętrzne instrukcje obliczeniowe dla fortyfikacji i ochrony aktywnej, częściowo bazujące na UFC i TM5-1300.

Brak normy cywilnej szczegółowej — Polska nie ma dedykowanej normy cywilnej dla blast protection budynków (inaczej niż USA). Budynki ambasad, obiektów NATO i infrastruktura krytyczna (ministerstwa, elektrownie) są projektowane według wytycznych zleceniodawcy (NATO, MON, operatorzy).

Polskie ośrodki badawcze: WAT (Instytut Techniki Uzbrojenia), Politechnika Rzeszowska (badania materiałów), OBRUM Gliwice (pojazdy opancerzone), Instytut Techniki Budowlanej (ITB, odpowiedź konstrukcji).

Dlaczego fala odbita jest wejściem do ściany

Jeżeli fala porusza się w wolnym polu, mówimy o nadciśnieniu padającym. Jeżeli trafia w sztywną przegrodę, ciśnienie na powierzchni może wzrosnąć przez odbicie. Dlatego kalkulator pressure-impulse używa fali odbitej jako obciążenia ściany. To jest zgodne z intuicją z artykułu o fali odbitej i obciążeniu konstrukcji: ściana nie mierzy tego samego, co czujnik w powietrzu obok niej.

Jednocześnie to także uproszczenie. Realna powierzchnia może być pod kątem do fali. Może mieć otwory, wnęki, okna, kolumny, szorstkość, warstwy i elementy podatne. W pomieszczeniu mogą pojawić się odbicia wielokrotne i narastanie ciśnienia. Rozdział o obciążeniach wybuchem w Safety of structures... podkreśla, że przy złożonej geometrii miejskiej albo wewnętrznej może być potrzebny model numeryczny powietrza, a nie tylko punktowy przebieg ciśnienia.¹

Tłumienie w modelu SDOF i jego znaczenie dla blast

Tłumienie w konstrukcjach budowlanych jest złożone: pochodzi z tarcia w złączach, pękania betonu, plastyczności stali, tłumienia hydrodynamicznego (woda w betonie), tłumienia radiatywnego (promieniowanie energii z fundamentów) i wymuszonego (poduszka powietrzna).

Typowe wartości ułamku tłumienia krytycznego ζ dla konstrukcji:

• Beton niezbrojony, kruchy: ζ ≈ 2–4%
• Beton zbrojony: ζ ≈ 5–7%
• Stal: ζ ≈ 2–5%
• Mur: ζ ≈ 7–10%
• Drewno: ζ ≈ 10–15%

Dla krótkich obciążeń blast (td < Tn) tłumienie ma małe znaczenie na maksymalnym ugięciu, bo element nie zdąży wykonać pełnego cyklu. Większe znaczenie ma dla:

• ugięcia resztkowego (po zakończeniu obciążenia),
• odpowiedzi na wielokrotne odbicia (faza negatywna, echo fali),
• drgań po uderzeniu, ważnych dla wytrzymałości zmęczeniowej połączeń.

W kalkulatorze lokalnym ζ = 5% jest wartością domyślną, bezpieczną dla żelbetu. Zmiana do 1% lub do 15% znacząco zmienia ugięcie resztkowe, ale nie maksymalne przy krótkim impulsie blast.

Odpowiedź MDOF przez dekompozycję modalną

Dla budynku lub płyty wieloprzęsłowej zamiast jednego modelu SDOF można użyć superpozycji modalnej:

1. Wyznaczyć mody własne i częstości ωᵢ przez analizę modalną (FEM lub analitycznie).
2. Zrzutować obciążenie p(t) na każdy mod przez całkowanie z kształtem modowym Φᵢ.
3. Rozwiązać niezależne równanie SDOF dla każdego modu i:
   Mᵢ qᵢ'' + Cᵢ qᵢ' + Kᵢ qᵢ = Fᵢ(t)
4. Zsumować odpowiedzi ze wszystkich modów: u(t) = Σ Φᵢ qᵢ(t)

Dla blast pierwotny mod dominuje zazwyczaj, bo kształt obciążenia blast (jednorodny po powierzchni) najlepiej współgra z pierwszym modem kształtowym (zwykle symetrycznym). Wyższe mody mogą być istotne przy niesymetrycznych obciążeniach lub gdy źródło blast leży blisko połowy rozpiętości.

Metoda superpozycji modalnej w postaci SRSS (Square Root Sum of Squares) lub CQC (Complete Quadratic Combination) jest standardem w analizie sejsmicznej, ale dla blast problem jest inny: obciążenie jest deterministyczne i krótkie, więc pełne całkowanie w dziedzinie czasu (time-history analysis) daje dokładniejszy wynik niż SRSS.

Lokalny kalkulator pressure-impulse pomija mody wyższe i używa tylko jednego stopnia swobody. To jest dozwolone jako pierwsze przybliżenie dla elementów jednorozpiętościowych (pola między słupami, jednotrzęsłowe panele). Kalkulator odporności płyty idzie o krok dalej i sumuje mody płyty prostokątnej, ale dla dużych płyt, ram wielokondygnacyjnych, nieliniowych podpór i uplastycznienia nadal konieczna jest pełna analiza dynamiczna FEM.

Metoda całkowania Duhamela

Dla arbitralnego kształtu obciążenia p(t) odpowiedź liniowego SDOF można zapisać przez całkę splotową (całka Duhamela):

y(t) = (1/m·ωₐ) ∫₀ᵗ p(τ) · e^{-ζω(t-τ)} · sin(ωₐ(t-τ)) dτ

gdzie ωₐ = ω√(1-ζ²) to częstość tłumionych drgań własnych.

Fizyczna interpretacja: całka Duhamela dodaje wkłady impulsów elementarnych p(τ)dτ do odpowiedzi. Każdy elementarny impuls wzbudza tłumione drgania własne. Superpozycja daje pełną odpowiedź.

Dla trójkątnego impulsu blast (czas wzrostu 0 do Pmax i liniowy zanik od Pmax do 0 w czasie td) całka Duhamela ma analityczne wyrażenie, ale przy małym tłumieniu łatwiej numerycznie Newmark-beta z małym krokiem. Kalkulator lokalny używa Newmark-beta z krokiem dt = min(Tn/100, 0,05 ms) dla dobrej dokładności.

Szyby jako wrażliwe elementy — szczegółowy przykład

Szyby są kluczowym elementem oceny skutków wybuchu w środowisku miejskim, bo pękają przy bardzo niskim nadciśnieniu (1–3 kPa = ok. 0,15–0,45 psi) i stanowią główne źródło odłamków wtórnych mogących zabić lub poważnie zranić ludzi wewnątrz budynku.

Model SDOF dla szyby float (10 mm, 1 m × 1,5 m):

• Masa: 25 kg/m² × 1,5 m² = 37,5 kg
• Naturalny okres drgań: około 80–100 ms (zależnie od grubości i wymiarów)
• Sztywność (metoda FEM lub normy ISO): kilkaset kN/m
• Ugięcie przy pęknięciu: ok. 3–5 mm dla tafli float

Fala blast o Ps = 5 kPa i Is = 30 kPa·ms (odległość ~50 m od wybuchu 100 kg TNT) spowoduje:

• czas trwania: td = 2 × 30/5 = 12 ms
• td/Tn ≈ 12/100 = 0,12 — silnie impulsowy reżim
• DLF ≈ 0,12 (obciążenie impulsowe — mało skuteczne)
• Ale impuls przenosi prędkość v₀ = I/m = 30/37,5 × 1,5 ≈ 0,8 mm/ms = 0,8 m/s

Taka prędkość może wystarczyć do przekroczenia ugięcia granicznego szyby. Właśnie dlatego szyby pękają przy stosunkowo małych nadciśnieniach — bo reagują na impuls, nie na pik.

Laminowane szyby z PVB (Polyvinyl Butyral) nie pękają przy tych samych impulsach, bo PVB pochłania energię i utrzymuje odłamki. Normy ISO 16933 klasyfikują szyby ochronne od klasy 1 (minimalna ochrona) do klasy 4 (wysokie zagrożenie), gdzie klasy wyższe odpowiadają wyższej asymptoty impulsowej na diagramie P-I.

Praktyczne tabele: progi uszkodzeń dla różnych elementów

Tabela orientacyjna progów nadciśnienia padającego Ps dla typowych elementów przy standardowych geometriach eksplozji:

Element	Próg lekkiego uszkodzenia [kPa]	Próg zniszczenia [kPa]	Dominujący parametr
Szyba float 4 mm	0,5–1	2–5	Impuls
Szyba float 10 mm	1–3	5–15	Impuls
Drzwi drewniane	3–5	15–30	Ciśnienie + impuls
Ściana murowana	10–20	30–50	Ciśnienie
Lekki panel stalowy	5–10	30–70	Impuls + duktylność
Ściana RC 15 cm	30–50	100–200	Ciśnienie + impuls
Ściana RC 30 cm	70–100	200–500	Ciśnienie

Wartości bardzo orientacyjne — zależą od wymiarów, zbrojenia, podparcia i sposobu obciążenia. Tabele z UFC 3-340-02 są bardziej szczegółowe i walidowane testami.

Progi te pokazują, dlaczego w miejskim środowisku po wybuchu (np. bomba samochodowa) okna pękają w promieniu kilkuset metrów, a ściany są uszkodzone tylko w pobliżu epicentrum.⁴

Efekty dynamiczne a odpowiedź quasi-statyczna: DLF w praktyce

Dynamic Load Factor (DLF) jest stosunkiem maksymalnej dynamicznej odpowiedzi do statycznej odpowiedzi pod tym samym ciśnieniem szczytowym. Dla trójkątnego impulsu:

• td/Tn << 1 (impulsowy): DLF → td/Tn × π — małe DLF, element reaguje na impuls, nie na pik
• td/Tn ≈ 1 (rezonans): DLF → ~1,8 (dla ζ=0) — wzmocnienie dynamiczne
• td/Tn >> 1 (quasi-statyczny): DLF → 2 dla prostokąta, 1 dla wolno narastającego

Dla typowego impulsu blast (td/Tn między 0,1 a 0,5):

• DLF jest poniżej 1 w reżimie impulsowym — pik ciśnienia jest za krótki na pełne ugięcie statyczne
• ale samo ugięcie nie jest małe — zależy od impulsu

To wyjaśnia, dlaczego:

1. Szyba reaguje na impuls (małe Tn, duże td/Tn) — DLF bliskie jedności lub większe
2. Masywna ściana z długim Tn ma małe DLF przy tym samym blast — dostaje mało ugięcia, ale za to długo
3. Dla elementu rezonansowego (td ≈ Tn) DLF może przekroczyć 1,5–2 — efekt dynamiczny szkodliwy

Dydaktycznie DLF jest kluczowym pojęciem łączącym statykę z dynamiką. Student który rozumie DLF rozumie, dlaczego konstrukcja nie może być oceniana tylko ciężarem statycznym obciążenia.²

Czego model SDOF nie widzi

Model jednego stopnia swobody nie widzi lokalnego zniszczenia. Może wskazać duże ugięcie, ale nie powie, czy pękło połączenie, czy odpadła warstwa elewacyjna, czy zbrojenie weszło w uplastycznienie, czy nastąpiło przebicie lokalne. Nie rozróżnia trybów drgań płyty, pracy membranowej, propagacji rysy ani efektów szybkości odkształcenia.

Nie widzi też pełnego rozkładu przestrzennego obciążenia. Zakłada, że rozpatrywany element można obciążyć jednym przebiegiem p(t) na jednostkę powierzchni. Jeżeli fala dochodzi ukośnie, jeżeli front ma duży gradient po powierzchni albo jeżeli element jest duży względem długości charakterystycznej fali, takie założenie może być zbyt ubogie.

Właśnie dlatego artykuł o od nadciśnienia do ugięcia płyty jest ważnym uzupełnieniem. Pełniejsza analiza LS-DYNA z modelem Brode'a i materiałem Johnson-Cook pokazuje, czego model SDOF nie obejmuje: siatki elementów, lokalnych naprężeń, warunków brzegowych, szybkości odkształcenia, uszkodzeń i pracy materiałowej.³

Jak czytać wykres P-I

Wykres P-I najlepiej czytać jako mapę pytań, a nie gotową odpowiedź. Lewa część, przy małym impulsie, mówi o krótkich obciążeniach. Dolna część, przy małym ciśnieniu, mówi o słabych obciążeniach. Krzywa uszkodzenia oddziela kombinacje, które dla danego elementu są jeszcze akceptowalne, od kombinacji wymagających większej deformacji.

Dla lekkiego panelu krzywa może leżeć zupełnie inaczej niż dla masywnej ściany. Dla ściany kruchej inaczej niż dla ściany zdolnej do pracy plastycznej. Dla elementu dobrze zamocowanego inaczej niż dla elementu z podatnymi połączeniami. Dlatego nie ma jednej uniwersalnej krzywej P-I dla „budynku”. Są krzywe dla konkretnych elementów, założeń materiałowych i stanów granicznych.

Lokalny kalkulator nie rysuje jeszcze pełnych krzywych izouszkodzenia. Zamiast tego rozwiązuje jeden punkt: dane P, I, m, Tn, zeta i yy dają jedną odpowiedź dynamiczną. To dobry etap pośredni. Pełna krzywa powstałaby przez wielokrotne przeliczenie siatki punktów P-I i zapisanie granicy, przy której mu osiąga wybrany poziom.

Przykład myślowy

Wyobraźmy sobie dwie fale działające na tę samą ścianę. Pierwsza ma wysoki pik, ale trwa bardzo krótko. Druga ma niższy pik, ale działa kilka razy dłużej. Jeżeli patrzymy tylko na Pmax, pierwsza wydaje się groźniejsza. Jeżeli patrzymy na impuls i okres własny ściany, odpowiedź może się odwrócić.

Dla bardzo sztywnej, lekkiej przegrody wysoki pik może być decydujący lokalnie. Dla masywnej ściany o dłuższym okresie własnym ważniejsze może być, ile pędu i energii przekazano w czasie. Dla elementu z ograniczoną duktylnością najważniejsze może być, czy odpowiedź przekracza próg pracy plastycznej i ile razy.

Ta intuicja ma znaczenie także poza bronią jądrową. Dotyczy ochrony infrastruktury, ścian osłonowych, paneli, szyb, bram, drzwi, tuneli i elementów wyposażenia. W serwisie jest jednak szczególnie użyteczna dydaktycznie, bo pozwala przejść od bardzo ogólnego języka skutków wybuchu do mechaniki konkretnego elementu.

Weryfikacja modelu SDOF przez porównanie z testami

Jak sprawdzić, czy model SDOF daje wiarygodne wyniki? Standardowa metoda to porównanie z testami blastowymi na elementach w skali rzeczywistej lub zredukowanej.

Ośrodki prowadzące testy blast konstrukcji:

• US Army ERDC (Engineer Research and Development Center, Vicksburg, MS) — testy płyt, ścian, paneli na polach blastowych. Wyniki w raportach ERDC/GSL TR.
• DRDC (Defence Research and Development Canada, Suffield, AB) — testy schronów, kontenerów, szyb, drzwi.
• UK DSTL (Defence Science and Technology Laboratory, Porton Down) — testy odporności blast budynków (rządowych, dyplomatycznych).
• TNO (Holandia) — testy przemysłowe, BLEVE, wybuchy chemiczne.

Typowy wynik walidacji: model SDOF z K-B i modelem trójkątnym P(t) przewiduje ugięcie maksymalne z dokładnością ±20–30% dla prostych płyt jednorozpiętościowych. Dokładność spada przy złożonych podparciach, wyraźnej pracy membranowej i dużych ugięciach plastycznych.

Polska: WAT prowadzi testy elementów konstrukcji na własnym polu testowym. Wyniki są częściowo publikowane w Biuletynie WAT i materiałach konferencji MEPS (Military Engineering Problem Solving).⁴

Jak projektować ćwiczenia bezpiecznie

Ćwiczenia do tego artykułu powinny być całkowicie obliczeniowe albo wykonywane na danych syntetycznych. Nie ma potrzeby generowania wybuchu, fali uderzeniowej ani dużego impulsu ciśnienia. Najlepszym materiałem dydaktycznym jest tabela P, I, Tn, m/A, zeta i yy, na której można zobaczyć zmianę odpowiedzi.

W laboratorium można zbudować prosty analog mechaniczny: mała masa na sprężynie i tłumiku, pobudzana krótkim impulsem elektromagnesu albo stuknięciem o kontrolowanej sile. Nie jest to model blastu, lecz model układu dynamicznego. Pokazuje, że ten sam maksymalny nacisk przy innym czasie trwania daje inną odpowiedź. To wystarczy do zrozumienia sensu SDOF bez niebezpiecznej aparatury.

Drugi bezpieczny wariant to symulacja w arkuszu. Studenci implementują równanie m y'' + c y' + k y = p(t), używają trójkątnego p(t) i porównują odpowiedź dla różnych td/Tn. Potem zestawiają wyniki z lokalnym kalkulatorem i sprawdzają, czy rozumieją różnice jednostek.

Odpowiedź konstrukcji na blast jądrowy vs. konwencjonalny

Dla wybuchu jądrowego parametry fali blast są jakościowo takie same jak dla konwencjonalnych materiałów wybuchowych — tym samym można stosować modele P-I i SDOF. Różnice wynikają z:

1. Skali. Wybuch 1 kt TNT tworzy strefę poważnych uszkodzeń o promieniu ~300 m. Wybuch 1 Mt TNT — ~3 km. Impuls rośnie z W^(1/3), więc przy tej samej odległości skalowanej Z impuls na metr kwadratowy jest identyczny, ale w jednostkach bezwzględnych jest ogromny.

2. Czasu trwania fali. Dla dużych ładunków czas dodatniej fazy td jest proporcjonalny do W^(1/3). Dla 1 Mt przy Z = 3 m/kg^(1/3) czas td wynosi setki milisekund. To oznacza, że praktycznie każda budowla reaguje quasi-statycznie — DLF zbliża się do 2, impuls ma mniejsze znaczenie od ciśnienia.

3. Kombinacji efektów. W wybuchu jądrowym jednocześnie z blastm pojawia się promieniowanie cieplne, prompt radiation, EMP i opad. Modele P-I opisują tylko blast. Łączny efekt wybuchu jądrowego na konstrukcję i na ludzi wymaga oddzielnych analiz dla każdego efektu.

4. Skażenia i ewakuacji. Ochrona konstrukcyjna przed blastem nie gwarantuje ochrony przed odpadem promieniotwórczym, który może wnikać przez zniszczone przegrody lub systemy wentylacji.

Dlatego w literaturze o ochronie cywilnej przed bronią jądrową (Glasstone i Dolan, Raport Sołtan dla PL) rozdziały o schronach łączą wymagania blast-ochronne z wymaganiami filtracji powietrza, grubości zasłony gamma i autonomiczności zaopatrzenia. Model SDOF jest tylko jednym składnikiem tej złożonej analizy.

Prognozy strat ludności — od P-I do maps casualty

Wyniki modeli P-I dla poszczególnych elementów można agregować do prognoz strat ludności w strefach. Uproszczone podejście:

1. Dla każdej odległości od epicentrum obliczyć Ps i Is (K-B/CONWEP).
2. Korzystając z krzywych P-I dla typowych budynków, określić procent budynków w stanie lekkiego, średniego i ciężkiego uszkodzenia.
3. Dla każdego stanu uszkodzenia przyjąć procent ofiar wewnątrz.
4. Zsumować po odległościach i typach zabudowy.

Taka metodologia, opisana w NSO/ENA STANAG 7196 i stosowana przez NATO Nuclear Planning Group, używa krzywych P-I jako pośredniego kroku między blast physics a demografią strat. Kalkulatory skutków nuklearnych (np. NukeMap by Alex Wellerstein) używają podobnej logiki, choć uproszczonej.

W Polsce prognozy strat ludności na wypadek ataku jądrowego były opracowywane w latach zimnej wojny przez Urząd Obrony Cywilnej. Dokumenty te są klasyfikowane lub wewnętrzne, ale metoda P-I/SDOF jest standardowa i znana w środowisku badaczy bezpieczeństwa.¹

Minimalny raport z obliczenia

Każdy wynik P-I/SDOF powinien podawać przynajmniej:

• źródło parametrów fali: K-B/CONWEP, Brode, dane syntetyczne albo pomiar,
• czy użyto fali padającej czy odbitej,
• Pmax, I, td i jednostki,
• masę powierzchniową, okres własny i tłumienie,
• założony próg ugięcia granicznego,
• ymax, DLF, mu i interpretację,
• zastrzeżenie, czego model nie obejmuje.

Bez takiej metryki wynik jest trudny do powtórzenia. Szczególnie łatwo pomylić kPa*ms z kPa*s, falę padającą z odbitą albo masę całej ściany z masą powierzchniową. Te drobne błędy potrafią zmienić wynik bardziej niż elegancki wybór metody całkowania.

Dodatkowe materiały multimedialne

Powiązane kalkulatory i narzędzia

• Pressure-Impulse dla ściany - lokalny solver SDOF dla fali odbitej i trójkątnego obciążenia.
• Podmuch - źródło P, I, czasu dodatniej fazy i fali odbitej.
• Odporność płyty - kolejny poziom, z geometrią i materiałem płyty.
• Równoważnik TNT - przejście od energii do efektywnej masy TNT.
• Height-of-burst i Mach stem - wpływ geometrii wybuchu powietrznego na falę przy gruncie.

Ćwiczenia praktyczne

Pierwsze ćwiczenie powinno polegać na porównaniu dwóch obciążeń o tym samym piku. Należy:

1. przyjąć Pmax = 100 kPa,
2. policzyć odpowiedź SDOF dla td = 5 ms, 20 ms i 80 ms,
3. zachować tę samą masę powierzchniową, okres własny i tłumienie,
4. porównać I, ymax, DLF i mu,
5. wyjaśnić, dlaczego sam Pmax nie wystarczył do przewidzenia wyniku.

Drugie ćwiczenie powinno dotyczyć stosunku td/Tn. Należy:

1. wybrać jeden impuls trójkątny,
2. policzyć odpowiedź dla Tn = 10 ms, 40 ms, 100 ms i 250 ms,
3. wskazać reżim impulsowy, pośredni i quasi-statyczny,
4. sprawdzić, gdzie DLF najbardziej odbiega od 1,
5. zapisać, kiedy wartość mu przekracza przyjęty próg.

Trzecie ćwiczenie powinno być audytem lokalnego kalkulatora. Należy:

1. otworzyć Pressure-Impulse,
2. zmienić tylko masę powierzchniową ściany,
3. następnie zmienić tylko okres własny,
4. następnie zmienić tylko tłumienie,
5. opisać, który parametr najmocniej zmienił ugięcie i dlaczego.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Od nadciśnienia do impulsu: praktyczne modele Kingery-Bulmash, Brode i CONWEP
• Mach stem i wysokość wybuchu: pomiary mikrofonowe, szybka kamera i symulacja
• Fala odbita i obciążenie konstrukcji
• Od nadciśnienia do ugięcia płyty: numeryczne modelowanie skutków fali podmuchowej
• Równoważnik trotylowy, nadciśnienie i impuls