Od nadciśnienia do impulsu: praktyczne modele Kingery-Bulmash, Brode i CONWEP

Streszczenie

Model fali podmuchowej nie odpowiada tylko na pytanie, „jak wysokie było ciśnienie”. Dla konstrukcji, czujnika, szyby, ściany albo płyty stalowej równie ważne jest, jak długo ciśnienie działało, czy fala była padająca czy odbita, jaki był impuls dodatniej fazy i czy punkt leży w zakresie ważności użytego modelu. Modele Kingery-Bulmash, Brode i CONWEP są praktycznymi narzędziami przejścia od energii wybuchu, geometrii i odległości do takich wielkości jak Ps, Pr, Is, Ir, czas przybycia i przebieg P(t).^1,2,3

W serwisie ten temat pełni rolę mostu między prostym równoważnikiem trotylowym, kalkulatorem Podmuch i późniejszymi artykułami o pressure-impulse, SDOF, wysokości wybuchu oraz symulacjach konstrukcyjnych. Artykuł nie jest instrukcją projektowania wybuchu ani narzędziem do optymalizacji szkód. Pokazuje, jak czytać modele, jak nie mylić piku nadciśnienia z impulsem i gdzie kończy się sensowność prostych kalkulatorów.

Rozszerzenie tematu

Najprostszy wykres fali podmuchowej wygląda niewinnie: ciśnienie skacze nagle do wartości maksymalnej, potem spada, a po chwili może przejść w słabszą fazę podciśnienia. W praktyce ten jeden wykres zawiera kilka różnych informacji. Szczytowe nadciśnienie mówi, jak gwałtowny był front. Impuls mówi, jak dużo „pchnięcia” przekazano w czasie. Czas przybycia mówi, kiedy front dotarł do punktu obserwacji. Czas dodatniej fazy mówi, czy obciążenie było krótkim uderzeniem, czy dłuższym naciskiem.^1,2

To rozróżnienie jest kluczowe. Cienka szyba, lekka ściana, masywna płyta żelbetowa i czujnik ciśnienia nie reagują na tę samą falę identycznie. Dla bardzo lekkiego elementu decydujący może być pik ciśnienia. Dla cięższego elementu o własnym okresie drgań często ważniejszy staje się impuls i stosunek czasu dodatniej fazy do czasu własnego konstrukcji. Dlatego w inżynierii skutków wybuchu szybko przechodzi się od samej wartości Ps do przestrzeni P-I, czyli pressure-impulse.⁴

W artykułach o skutkach broni jądrowej łatwo spotkać wartości progowe typu 3 psi, 5 psi, 20 psi albo 70 psi. Są użyteczne, ale bywają zdradliwe. Progi nadciśnienia dobrze nadają się do szybkiej mapy skutków, lecz nie zastępują czasu działania fali, fali odbitej, geometrii celu i odpowiedzi dynamicznej. Właśnie dlatego lokalny kalkulator Podmuch pokazuje nie tylko nadciśnienie padające, ale też nadciśnienie odbite, impuls, czas trwania dodatniej fazy, czas przybycia i orientacyjny przebieg Friedlandera.⁵

Fizyczna interpretacja parametrów fali blast: skąd biorą się liczby

Zanim przejdziemy do modeli, warto zadać pytanie: dlaczego Ps, Pr, Is i Ir są ważne i co fizycznie oznaczają w kontekście zachowania realnych obiektów?

Nadciśnienie (Ps, Pr) to siła na jednostkę powierzchni — przekraczając próg wytrzymałości materiału, powoduje pęknięcie, zawalenie lub odkształcenie plastyczne. Ciśnienie atmosferyczne to 101 kPa (1 atm). Już kilka kPa (0,5–1 psi) powoduje rozbicie szyb w oknach. Przy ok. 70 kPa (10 psi) typowe budynki murowane ulegają poważnemu zniszczeniu. Przy 350 kPa (50 psi) większość zabudowy ulega całkowitemu zniszczeniu.

Impuls (Is, Ir) to całka ciśnienia po czasie — opisuje, ile łącznego „pchnięcia" otrzymała jednostka powierzchni. Siła pomnożona przez czas równa się zmianie pędu. Dlatego impuls lepiej opisuje odpowiedź masywnych obiektów, które nie zdążyły się poruszyć w czasie krótszym niż ich własny okres drgań. Lekki panel reaguje na szczyt ciśnienia, masywna ściana żelbetonowa reaguje na impuls.

Czas przybycia (Ta) decyduje o czasie ostrzeżenia — dla prędkości dźwięku (ok. 340 m/s) 1 km to ok. 3 sekundy. Fala blast przy silnych eksplozjach porusza się szybciej niż dźwięk: dla Ps = 100 kPa prędkość czoła wynosi ok. 500 m/s.

Czas dodatniej fazy (td, t+) wyznacza trwanie obciążenia. Jeśli td jest dużo krótsze od okresu własnego konstrukcji, element reaguje impulsowo (jak uderzenie młotem). Jeśli td jest porównywalne lub dłuższe od okresu własnego, element reaguje jak na obciążenie quasi-statyczne. Właśnie ta relacja jest tematem artykułu o pressure-impulse i SDOF.²

Odległość skalowana

Wspólnym językiem większości praktycznych modeli jest skalowanie Hopkinsona-Cranza:

Z = R / W^(1/3).

R oznacza odległość od źródła, a W masę lub równoważną masę TNT. Jeżeli dwa przypadki mają podobne Z, to w pierwszym przybliżeniu mogą dawać podobny kształt fali, mimo że jeden dotyczy małego ładunku w laboratorium, a drugi większego zdarzenia w terenie. To nie jest prawo doskonałe, lecz bardzo użyteczna reguła porządkowania danych eksperymentalnych.^1,3

Trzeba jednak uważać na jednostki. W lokalnym kalkulatorze fali podmuchowej model Kingery-Bulmash liczy Z w m/kg^(1/3), bo wejściowy uzysk w kilotonach jest przeliczany na kilogramy TNT. W literaturze spotyka się też stopy, funty, metry na kilogram do potęgi 1/3, metry na kilotonę do potęgi 1/3 albo wykresy referencyjne dla 1 kt. Sama postać R/W^(1/3) jest ta sama, lecz liczby w tabelach i współczynnikach zależą od konwencji jednostek.⁵

Skalowanie działa najlepiej, gdy porównujemy zjawiska tego samego rodzaju: falę w powietrzu z falą w powietrzu, wybuch powierzchniowy z wybuchem powierzchniowym, a nie dowolne zdarzenia tylko dlatego, że mają podobną energię. Inaczej zachowuje się wybuch swobodny, inaczej półsferyczny wybuch przy powierzchni, inaczej wybuch z kraterowaniem, inaczej fala odbita od dużej przegrody. Dlatego w kalkulatorze tryb air i surface nie są ozdobą interfejsu, lecz wyborem innego zestawu zależności.⁵

Kingery-Bulmash

Zależności Kingery-Bulmash powstały jako empiryczne ujęcie parametrów fali podmuchowej od sferycznego wybuchu w powietrzu i półsferycznego wybuchu powierzchniowego TNT. Dokument CONWEP opisuje, że równania użyte w programie oparto na raporcie ARBRL-TR-02555 Kingery'ego i Bulmasha z 1984 roku, zawierającym kompilację danych z testów z masami od mniej niż 1 kg do ponad 400 000 kg materiału wybuchowego. To ważny szczegół: K-B nie jest teorią pierwszych zasad, lecz krzywą dopasowaną do dużego zbioru danych.²

W zapisie zestawionym przez Trzcińskiego ogólna postać modelu jest wielomianowa:

Y = C0 + C1 U + C2 U^2 + ... + Cn U^n,

gdzie:

U = K0 + K1 log10(Z).

Y jest logarytmem dziesiętnym szukanego parametru fali, na przykład nadciśnienia padającego, nadciśnienia odbitego, impulsu padającego albo impulsu odbitego. Dla różnych parametrów i zakresów Z stosuje się różne zestawy współczynników. Trzciński podaje je w układzie SI, co ułatwia bezpośrednie użycie w kalkulatorze bez przepisywania anglosaskich jednostek z oryginalnych tabel.¹

W lokalnym kalkulatorze fali podmuchowej te współczynniki są używane jako oddzielne zestawy dla wybuchu swobodnego i powierzchniowego: Ps, Pr, Is oraz Ir. Nadciśnienie wynika bezpośrednio z wielomianu, natomiast impuls skalowany jest następnie przez W^(1/3). To odpowiada intuicji: przy tej samej odległości skalowanej większa eksplozja ma podobny charakter fali, ale działa dłużej, więc impuls rośnie ze skalą liniową zdarzenia.⁵

Zakres ważności także jest częścią modelu. W kodzie zakresy K-B są ograniczone do Z od około 0,0532 dla wybuchu swobodnego i 0,0674 dla powierzchniowego do 40 m/kg^(1/3). Jeżeli punkt wykracza poza ten zakres, kalkulator nie powinien udawać precyzji. Lokalna implementacja przechodzi wtedy awaryjnie na starszy model Kinney-Graham i wyświetla ostrzeżenie, że wybrany K-B nie obejmuje danego punktu.⁵

To podejście jest redakcyjnie ważne. Dydaktyczny kalkulator może pomagać w zrozumieniu zależności, ale nie może sugerować, że każdy wynik z dokładnością do dwóch miejsc po przecinku ma sens fizyczny. Przy bardzo małych Z, blisko źródła, dochodzą lokalne produkty detonacji, niejednorodność źródła i złożona geometria. Przy bardzo dużych odległościach coraz większe znaczenie mają warunki atmosferyczne, teren i tłumienie. Prosty model empiryczny nie znika, ale staje się tylko jedną warstwą analizy.^1,2

CONWEP jako program, nie osobna fizyka

CONWEP, czyli Conventional Weapons Effects Program, bywa przywoływany tak, jakby był oddzielnym modelem fizycznym. W rzeczywistości jest przede wszystkim implementacją zestawu procedur z podręczników projektowania ochronnego, obejmującą między innymi airblast, odłamki, penetrację i kraterowanie. Dla fali podmuchowej kluczowe są zależności Kingery-Bulmash oraz procedury z TM 5-855-1.²

Dokumentacja CONWEP podkreśla kilka założeń, które warto zapamiętać. Obliczenia airblast są sprowadzane do równoważnej masy TNT. Dla wybuchu powierzchniowego zakłada się kształt półsferyczny, a dla wybuchu w powietrzu kształt sferyczny. Wpływ obudowy ładunku traktowany jest ostrożnie, bo nie ma jednej prostej poprawki, która działałaby dla wszystkich przypadków. Faza ujemna fali jest w tej procedurze pomijana jako zwykle mniej istotna dla klasycznego projektowania ochronnego.²

CONWEP nie ogranicza się jednak do samego piku. Dokumentacja opisuje użycie zmodyfikowanego równania Friedlandera do odtworzenia historii ciśnienia w czasie:

P(t) = Pso * [1 - (t - Ta) / To] * exp[-A * (t - Ta) / To].

W tym zapisie Pso jest szczytowym nadciśnieniem padającym, Ta czasem przybycia, To czasem dodatniej fazy, a A bezwymiarowym współczynnikiem zaniku. Program korzysta z obliczonych wartości piku, impulsu i czasu trwania, aby dobrać kształt krzywej P(t). Dla użytkownika najważniejszy wniosek brzmi: impuls nie jest parametrem dekoracyjnym, lecz całką z tej historii ciśnienia.²

Lokalny kalkulator robi podobny krok dydaktyczny, choć prościej. K-B dostarcza Ps, Pr, Is i Ir, a funkcja friedlanderDurationFromImpulse() odwraca uproszczony kształt Friedlandera, aby z piku i impulsu odzyskać orientacyjny czas dodatniej fazy. Wykres w wynikach kalkulatora nie jest więc „nagraną falą”, lecz rekonstrukcją zgodną z podstawowymi parametrami modelu.⁵

Brode

Metoda Brode'a pojawia się szczególnie często w kontekście jądrowych fal podmuchowych i implementacji w kodach takich jak LS-DYNA. W pracy Kosiuczenki opisano ją jako empiryczny sposób wyznaczania szczytowego nadciśnienia z wykresów dla wybuchu powietrznego 1 kt w standardowych warunkach na poziomie morza, z uwzględnieniem skalowania Hopkinsona-Cranza i poprawek wysokościowych. Ten opis jest dobrym przykładem, jak historyczne wykresy i współczynniki stają się wejściem do współczesnej symulacji konstrukcyjnej.³

W cytowanej pracy porównywane są metody empiryczne, takie jak Brode i CONWEP, z bardziej złożonymi symulacjami ALE. Wniosek nie jest taki, że jedna metoda zawsze zastępuje drugą. Modele empiryczne są szybkie i bardzo użyteczne dla otwartej przestrzeni oraz prostych geometrii. Metody numeryczne lepiej radzą sobie z lokalną geometrią, wielokrotnymi odbiciami i złożonym przepływem, ale kosztują znacznie więcej obliczeniowo i wymagają danych materiałowych.³

Praca Kosiuczenki pokazuje też typowy łańcuch: najpierw dobiera się model fali, potem mapuje się obciążenie na elementy powierzchni, a następnie liczy odpowiedź płyty w LS-DYNA z modelem materiału. Lokalny kalkulator odpornosc-plyty zachowuje bezpieczny, dydaktyczny szkielet tej idei: nie dobiera parametrów źródła, tylko przyjmuje zadany przebieg p(t) i liczy modalną odpowiedź sprężystej płyty prostokątnej. Dzięki temu można zobaczyć, ile wnosi przejście z jednego stopnia swobody do kilku modów, a jednocześnie gdzie nadal zaczyna się potrzeba pełnej symulacji FEM.⁶

Co dokładnie oznaczają Ps, Pr, Is i Ir

Ps to szczytowe nadciśnienie fali padającej w wolnym polu. To wartość, którą zwykle rysuje się jako izobary na mapie skutków. Pr to nadciśnienie odbite od powierzchni ustawionej normalnie do fali albo sprowadzone do takiego przypadku w uproszczonym modelu. Przy silniejszej fali Pr może być znacznie większe niż Ps, dlatego ściana, płyta i czujnik na powierzchni „widzą” inne obciążenie niż punkt w powietrzu.^2,5

Is to impuls dodatniej fazy fali padającej. W najprostszym sensie jest to pole pod wykresem nadciśnienia w czasie. Ir jest analogicznym impulsem dla fali odbitej. Jednostki w kalkulatorach mogą wyglądać różnie: kPa*ms, kPa*s, psi-ms, Pa*s. Zawsze chodzi jednak o ciśnienie pomnożone przez czas, czyli o wielkość bliższą przekazowi pędu niż samej wartości maksymalnej.^2,5

Kiedy obciążenie trafia na konstrukcję, te wielkości nie są równoważne. Wysoki Ps i mały Is może oznaczać gwałtowne lokalne uszkodzenie cienkiego elementu. Niższy Ps, ale większy Is, może mocno rozpędzić większą masę. To dlatego artykuł o pressure-impulse i modelu SDOF jest naturalną kontynuacją tego tekstu.

Kingery-Bulmash, Brode, Kinney-Graham i proste wzory

W źródłach spotyka się wiele zależności: Sadowskiego, Henrycha, Własowa, Kinney-Grahama, Brode'a, Kingery-Bulmasha, procedury NATO i polskie przepisy. Trzciński porównuje część z nich i pokazuje rzecz bardzo praktyczną: dla pewnych zakresów Z wyniki są tego samego rzędu, ale w innych przedziałach mogą rozjechać się bardzo mocno, zwłaszcza dla impulsu. W szczególności ostrzega, że niektóre wzory wyprowadzone z danych jądrowych nie powinny być bezpośrednio stosowane do ładunków TNT w nieodpowiednim zakresie.¹

To jest dobry argument za tym, aby nie traktować modeli jako zamiennych tabel magicznych. Model trzeba dobrać do rodzaju źródła, geometrii, zakresu skalowanej odległości i parametru, który nas interesuje. Jeżeli dwa modele dają podobne Ps, ale bardzo różne Is, to nie jest drobiazg zaokrąglenia. Dla konstrukcji taka różnica może oznaczać zupełnie inną odpowiedź dynamiczną.

W lokalnym kalkulatorze zachowano dwa tryby: Kingery-Bulmash / CONWEP oraz historyczny Kinney-Graham. Pierwszy jest preferowany tam, gdzie mieści się w zakresie danych K-B. Drugi jest przydatny jako tryb porównawczy i awaryjny. Użytkownik powinien patrzeć nie tylko na liczby, ale też na komunikat modelu i ostrzeżenie o zakresie ważności.⁵

Historia eksperymentów: od trotyle przez Nevada Test Site do baz danych

Modele empiryczne nie biorą się znikąd. Są wynikiem dziesiątek lat eksperymentów z materiałami wybuchowymi, konwencjonalnymi i jądrowymi.

WWII i eksperymenty Kirkwood-Bethe. Pierwsze systematyczne pomiary parametrów fali blast w powietrzu przeprowadzono podczas II wojny światowej, m.in. dla potrzeb projektowania bomb lotniczych. Kirkwood i Bethe (1942) opracowali teorię propagacji sferycznej fali uderzeniowej. Dane te stały się podstawą dla powojennych baz.

Nevada Test Site (NTS), lata 50.–60. Setki testów jądrowych i konwencjonalnych dostarczyły danych w nieosiągalnym laboratoryjnie zakresie energii. Dla każdego testu ustawiano sieci czujników ciśnienia wokół punktu wybuchu, mierzono czas przybycia, szczytowe nadciśnienie, czas dodatniej fazy i impuls. Wyniki kompilowane przez DASA (Defense Atomic Support Agency) i potem przez BRL (Ballistic Research Laboratory, Aberdeen Proving Ground) trafiły do raportów, z których Kingery i Bulmash opracowali swoje tabele wielomianowe w 1984 roku.²

Masa testowych ładunków: kilogramy do setek ton. Tabele K-B obejmują zakres mas od ułamka kilograma TNT (testy laboratoryjne) do ponad 400 ton TNT. Tak szeroki zakres pozwolił na wiarygodne sprawdzenie prawa skalowania Hopkinsona-Cranza w praktyce. Wniosek: prawo sześcienne działa dobrze dla fali w swobodnym powietrzu przez ok. 6 rzędów wielkości energii.

DIPRO, MABS, DDESB. NATO i US Army kontynuują zbieranie danych blast. DIPRO (Data on Internal and Proximate Blast) i MABS (Military Application of Blast Simulation) to konferencje i bazy danych. DDESB (Department of Defense Explosives Safety Board) utrzymuje standardy bezpieczeństwa. Z tych baz korzystają inżynierowie projektujący fortyfikacje, schrony i testujący sprzęt wojskowy na odporność blast.

Faza ujemna: niedoceniany parametr modeli

Klasyczne modele K-B i CONWEP pomijają fazę ujemną lub traktują ją jako drugorzędną. To uproszczenie bywa kosztowne.

Po przejściu czoła fali i pozytywnej fazy ciśnienie spada poniżej atmosferycznego, tworząc fazę ujemną (suction, underpressure). Typowo faza ujemna jest kilkakrotnie dłuższa od pozytywnej, ale ma mniejszy szczyt — rzędu kilku kPa nawet gdy pozytywna sięgała setek kPa.

Dla jakich obiektów faza ujemna jest ważna?

• Lekkie konstrukje dachowe. Dach hali lub magazynu najpierw jest obciążony od góry (pozytywne nadciśnienie), a chwilę później poderwany od dołu (podciśnienie). Kombinacja może spowodować awarię dachownicy przy nadciśnieniu poniżej progu normalnego zniszczenia.
• Szyby i fasady szklane. Szklany panel po pęknięciu w fazie pozytywnej może być „wciągany" na zewnątrz w fazie ujemnej, co zmienia wzór fragmentów.
• Elementy podwieszane. Elementy sufitu, oświetlenia, instalacji podwieszanych mogą być wyrwane w fazie ujemnej, choć pozytywna ich nie naruszyła.
• Rezonans konstrukcji. Jeżeli okres własny drgań elementu jest bliski dwukrotności td pozytywnej, dochodzi do rezonansu między kolejnymi wymuszeniami; uwzględnienie fazy ujemnej jako wymuszenia odwrotnego znacząco zmienia wyliczone maksymalne przemieszczenie.

W modelu Friedlandera pełnym (nie przyciętym do zera) parametr A wyznacza głębokość fazy ujemnej. Niektóre implementacje numeryczne kalibrują A tak, by zerowy impuls łączny był spełniony (warunek zachowania pędu). Inne po prostu przycinają P(t) = 0 dla t > ta + td, co jest przeuproszczeniem.^1,5

Równanie stanu i gaz doskonały: kiedy model się sypie

Wszystkie modele empiryczne zakładają powietrze jako jednorodny gaz doskonały (lub gaz van der Waalsa z poprawkami) i eksplozję jako natychmiastowy wzrost energii wewnętrznej. To działa dla odległości, w których fala jest już oddzielona od produktów detonacji i promieniowania.

Blisko centrum wybuchu (małe Z) te założenia zawodzą:

1. Produkty detonacji mają własne ciśnienie i strukturę chemiczną. Nie są tym samym co powietrze.
2. Kula ognista wybuchu jądrowego emituje promieniowanie cieplne. Ogrzewa powietrze przed falą — powstaje „termicznie wzmocnione" tło. Tego efektu nie ma w modelu TNT.
3. Jonizacja powietrza. Przy bardzo wysokich temperaturach bliskich wybuchowi jądrowemu powietrze jest zjonizowane, a równanie stanu gazu doskonałego przestaje obowiązywać. Stąd modele specjalnie tworzone dla broni jądrowej (np. Brode) różnią się od modeli TNT w pobliżu źródła.
4. Ukształtowanie terenu i zabudowa. Żaden uproszczony model 1D nie widzi wzgórza, kanału, ulicy ani budynku.

Dlatego dobrą zasadą jest: modele K-B/CONWEP stosować dla Z > 0,2 m/kg^(1/3) dla konwencjonalnych ładunków TNT w otwartym terenie. Dla mniejszych Z lub dla ładunków jądrowych albo specjalnych chemicznie — szukać dedykowanych danych lub symulacji.³

Porównanie modeli: tabela wybranych parametrów

Zestawienie orientacyjne dla Ps (nadciśnienie padające) jako funkcji Z, na podstawie Trzcińskiego i literatury:

Z [m/kg^(1/3)]	K-B (CONWEP)	Kinney-Graham	Brode (jądrowy)
0,2	~12 MPa	~15 MPa	~10 MPa
0,5	~2 MPa	~2,5 MPa	~1,8 MPa
1,0	~500 kPa	~550 kPa	~450 kPa
2,0	~100 kPa	~110 kPa	~90 kPa
5,0	~15 kPa	~18 kPa	~13 kPa
10,0	~4 kPa	~5 kPa	~3,5 kPa
40,0	~0,5 kPa	~0,6 kPa	~0,45 kPa

Wartości orientacyjne — dla konkretnych projeków korzystać z tabel K-B lub kodów CONWEP. Dla impulsu Is rozbieżności między modelami są większe, zwłaszcza przy małych Z. Trzciński ostrzega, że wzory Sadowskiego stosowane w polskich przepisach budowlanych mogą znacząco zaniżać impuls w porównaniu z K-B w pewnych zakresach.¹

Korekcja ekwiwalentu TNT dla różnych materiałów wybuchowych

Tabele K-B opracowano dla TNT. W praktyce stosuje się wiele innych materiałów wybuchowych o różnych właściwościach: ANFO, RDX, PETN, HMX, Comp B, Semtex, amonit. Każdy ma inny ekwiwalent TNT, który zależy od:

• energii właściwej detonacji (Q_det),
• prędkości detonacji (VOD),
• charakteru produktów (H₂O para czy ciecz, CO₂ czy CO itd.),
• trybu porównania (szczytowe nadciśnienie czy impuls).

Ważny niuans: ekwiwalent TNT dla szczytowego nadciśnienia blisko wybuchu (Z < 1) może być inny niż ekwiwalent dla impulsu na dalszych odległościach (Z > 3). RDX ma ekwiwalent ciśnienia ok. 1,60 względem TNT, ale ekwiwalent impulsu w pewnych zakresach jest inny. Dlatego „ekwiwalent TNT" jest zawsze przypisany do konkretnego parametru i odległości, a nie do materiału w ogólności.

W kalkulatorach i obliczeniach cywilnych zwykle stosuje się jeden konserwatywny ekwiwalent TNT — wyższy z dostępnych. To bezpieczne podejście dla ochrony budynków, ale może być mylące przy porównywaniu różnych ładunków w celach badawczych.¹

Modele dla fal wewnątrz budynków (internal blast)

K-B/CONWEP opisują falę w otwartej przestrzeni zewnętrznej. Przy wybuchu wewnątrz pomieszczenia lub tunelu zachowanie jest zupełnie inne:

1. Wielokrotne odbicia. Fala odbija się od ścian, stropu i podłogi, tworząc złożony przebieg ciśnienia z wieloma szczytami.
2. Quasi-statyczne ciśnienie. Produkty detonacji zajmują objętość pomieszczenia i podnoszą średnie ciśnienie powyżej atmosferycznego przez czas dłuższy niż td zewnętrznej fali. To quasi-statyczne ciśnienie jest często groźniejsze dla konstrukcji niż dynamiczny szczyt.
3. Skala efektów. W zamkniętym pomieszczeniu nawet mały ładunek (100 g TNT) może być śmiertelny z powodu wielokrotnych odbić i quasi-statycznego ciśnienia.

Dla modeli internal blast stosuje się inne procedury — np. ESP (Explosion Safety Procedures) DDESB lub podejście van Dorpa (Vapor Cloud Explosion). Lokalne kalkulatory nie omawiają internal blast, bo wymaga on osobnych danych i procedur.

Polska norma i standardy NATO

W Polsce projektowanie obiektów odpornych na blast regulowane jest przez:

• PN-EN 1991-1-7 — Oddziaływania wyjątkowe: obciążenia wybuchem. Podaje ogólne zasady obliczania obciążeń od zewnętrznych eksplozji chemicznych.
• STANAG 2280 i ATP-55 — Standardy NATO dla ochrony personelu i budynków w obiektach wojskowych przed wybuchem.
• Przepisy techniczne WAT/MON — Dla budynków fortyfikacyjnych i schronów wojskowych Ministerstwo Obrony Narodowej stosuje własne przepisy techniczne, bazujące na KingeryBulmash i danych NTS.

W praktyce polskich inżynierów budowlanych norma EN 1991-1-7 jest rzadko stosowana (wybuchy są zdarzeniami wyjątkowymi, dla których trudno zebrać dane), a obliczenia projektowe dla ochrony infrastruktury krytycznej wykonywane są na zlecenie BBN (Biuro Bezpieczeństwa Narodowego) lub ABW (Agencja Bezpieczeństwa Wewnętrznego) przez WAT lub wyspecjalizowane biura projektowe.³

Czego te modele nie widzą

Empiryczny model fali podmuchowej nie widzi miasta. Nie widzi ulicy-korytarza, lasu, wzgórza, okien, narożników budynku, wielokrotnych odbić ani lokalnego skupienia fali między przeszkodami. Nie widzi też jakości wykonania konstrukcji, zbrojenia, połączeń, starzenia betonu, kruchości szkła i sposobu zamocowania płyty. Te wszystkie czynniki mogą zmienić realny skutek, nawet jeżeli wejściowy Ps wygląda rozsądnie.^3,4

Model nie widzi również pełnego zjawiska jądrowego. Wybuch jądrowy daje falę podmuchową, ale równolegle występuje promieniowanie cieplne, prompt radiation, opad, jonizacja atmosfery i efekty wysokości detonacji. Sprowadzenie wszystkiego do TNT equivalent jest użyteczne tylko dla mechanicznej skali podmuchu. Nie opisuje poparzeń, dawki, skażenia ani impulsu elektromagnetycznego.⁷

Nawet w ramach samej mechaniki trzeba rozdzielać cel dydaktyczny od projektowego. Kalkulator może pomóc studentowi zobaczyć, dlaczego przy stałym Z fala ma podobny charakter, dlaczego impuls skaluje się inaczej niż pik i dlaczego fala odbita jest groźna dla płaskiej przegrody. Nie zastępuje badań poligonowych, walidacji czujników, pełnego modelu CFD/FEM ani normowego projektu obiektu ochronnego.

Jak wyznaczyć ekwiwalent TNT z dowolnego źródła energii

Zanim cokolwiek policzymy modelem K-B, musimy zamienić energię wybuchu na masę TNT. Procedura różni się zależnie od rodzaju źródła:

Materiały chemiczne (TNT, RDX, PETN itd.). Dane z kart bezpieczeństwa (SDS) lub tabel wartości energii detonacji. Masę W_TNT = W_ładunku × e_eq, gdzie e_eq to tabelaryczny ekwiwalent. Dla TNT e_eq = 1,0 (definicja). Dla PETN e_eq ≈ 1,27, dla RDX e_eq ≈ 1,60, dla ANFO e_eq ≈ 0,82.

Para skroplonego gazu (LPG, LNG, propan, butan). Energia wybuchu oparów zależy od stopnia odparowania i mieszania z powietrzem. Procedury TNO "Multi-Energy Method" lub Baker-Strehlow-Tang dają lepszy wynik niż prosty ekwiwalent TNT, ale do szybkich szacunków stosuje się e_eq ≈ 0,3–0,5 (uwzględniając, że tylko część gazu bierze udział w deflagracji lub detonacji).

Wybuchy jądrowe. Konwencja: 1 kt TNT = 4,184 × 10¹² J. Dla wybuchu powietrznego ok. 50% energii idzie w falę blast, reszta w promieniowanie (cieplne i jonizujące). Efektywny ekwiwalent blast (TNTY) jest więc mniejszy niż całkowita energia kilograma. W tablicach Glasstone'a-Dolana dla wybuchu na poziomie morza standardowo przyjmuje się W_TNT_blast = 0,5 × Y × 10^6 kg, gdzie Y jest w kilotonach — czyli 50% całkowitej energii trafia do fali podmuchowej.

Katastrofy przemysłowe (BLEVE, UVCE). Metodologia TNO lub API RP 2Z definiuje procedury. Wartości e_eq bardzo różnią się od „idealnego" wybuchu materiału wybuchowego, bo deflagracja ma niższy Pmax i dłuższy td niż detonacja.

Diagram P-I jako wizualne narzędzie analizy skutków

Diagram pressure-impulse (P-I) jest narzędziem łączącym modele K-B z analizą skutków. Na osi X impuls, na osi Y nadciśnienie. Na tym wykresie wykreśla się:

1. Krzywą odpowiedzi obiektu — zestaw par (P, I) powodujących dane kryterium (np. 5% ryzyko śmierci, zniszczenie okna, złamanie belki). Krzywa ta jest asymptotyczna: po lewej stronie leży obszar impulsowy (krótkie obciążenie), po prawej — quasi-statyczny (długie).

2. Punkt z modelu K-B — dla danej odległości R i masy W kalkulator daje Ps i Is. Ten punkt leżący powyżej krzywej odpowiedzi oznacza, że skutek zostanie osiągnięty; poniżej — nie.

3. Izolinie odległości — punkty dla tej samej odległości skalowanej Z ale różnych W tworzą charakterystyczne krzywe na diagramie P-I.

Diagram P-I wyjaśnia, dlaczego nie można oceniać skutków samym nadciśnieniem. Lekkie elementy (witryny, gipsowe płyty sufitowe) reagują impulsowo — leżą w lewym obszarze diagramu. Masywne ściany żelbetowe leżą w prawym obszarze quasi-statycznym. Dla danej eksplozji punkt (Ps, Is) może leżeć „bezpiecznie" dla grubej ściany, ale w obszarze zniszczenia dla okna. To właśnie wyjaśnia, dlaczego progi nadciśnienia podawane dla całych klas budynków są tylko orientacją, a nie dokładną predykcją.

Lokalny kalkulator pressure-impulse implementuje uproszczoną wersję tej idei dla dydaktyki — bez pełnych krzywych uszkodzeń, ale z intuicją o zależności od stosunku td/T_n (czas trwania do okresu własnego).^4,6

Numeryczne modele blast: kiedy K-B nie wystarczy

Empiryczne modele K-B/CONWEP działają doskonale dla symetrii sferycznej lub półsferycznej w otwartej przestrzeni. Gdy geometria staje się złożona, sięga się po metody numeryczne:

ALE (Arbitrary Lagrangian-Eulerian) w AUTODYN lub LS-DYNA. Powietrze i produkty detonacji traktowane jako ciągły fluid. Mesh eulerowski ślizga się przez materiał. Możliwe modelowanie budynku w polu blast, kilku ładunków, ściany ekranującej. Czas obliczeniowy: od minut (2D) do dni (3D pełna geometria).

SPH w AUTODYN lub OpenFOAM. Bezsiatkowe cząstki śledzą ruch gazu. Dobre dla dużych odkształceń i fragmentacji. Problem: trudniejsza kontrola rozdzielczości i dokładności.

Sprzężone ALE+FEM. Powietrze modelowane ALE, konstrukcja modelowana FEM (belki, płyty, powłoki). Sprzężenie przez warunki brzegowe ciśnienia na powierzchni. Standard w projektowaniu ochronnym budynków rządowych i bunkrów.

OpenFOAM + compressibleInterFoam. Open-source, możliwa symulacja blast w miejskiej zabudowie. Stosowany przez uczelnie europejskie do badań naukowych. Mniej „certyfikowany" niż AUTODYN dla zastosowań wojskowych.

Kryteria wyboru narzędzia:

• Prosty zasięg nadciśnienia w otwartym terenie → K-B/CONWEP
• Zasięg i czas przybycia w terenie z jedną przeszkodą → K-B + korekta geometryczna lub 2D ALE
• Pełna odpowiedź budynku → ALE+FEM
• Planowanie ewakuacji lub wstępna ocena ryzyka → K-B + progi z GlasstoneDolan

Błędy i pułapki przy korzystaniu z K-B/CONWEP

Błąd 1: Pomijanie jednostek. K-B w oryginalnym raporcie 1984 używa angloamerykańskich jednostek (stopy, funty, psi, ms·psi). Implementacje metryczne wymagają przeliczenia współczynników. Lokalny kalkulator używa jednostek SI (metry, kilogramy, kPa), ale cytowane tabele K-B z innych źródeł mogą być w różnych układach.

Błąd 2: Traktowanie Ps jako obciążenia budynku. Ps to nadciśnienie w wolnym polu (free-field). Ściana prostopadła do fali dostaje Pr — kilkakrotnie wyższe. Płaski dach dostaje coś pośredniego. Czujnik „na wprost" eksplozji dostaje Pr. Używanie tylko Ps systematycznie zaniża obciążenie ścian.

Błąd 3: Użycie K-B poza zakresem ważności. Dla małych Z < 0,1 m/kg^(1/3) (blisko wybuchu) błędy mogą sięgać rzędu wielkości. Dla dużych Z > 40 m/kg^(1/3) (daleko, słabe fale) modele K-B też mają gorszą dokładność, bo fale zbliżają się do fal akustycznych i wrażliwość na atmosferę rośnie.

Błąd 4: Ignorowanie warunków atmosferycznych. K-B dotyczy standardowej atmosfery na poziomie morza (P₀ = 101 kPa, T₀ = 15°C). Dla terenów górskich (np. 2000 m n.p.m.) ciśnienie jest ok. 80 kPa — fala blast propaguje się inaczej. Korekta gęstości powietrza jest potrzebna, choć uproszczona.

Błąd 5: Zamiana modeli bez sprawdzenia zakresu. Kinney-Graham i K-B dają podobne Ps dla Z = 1–10, ale dla bardzo małych Z KG może zawyżać, a K-B niedoszacowywać. Najgorsza decyzja to wybranie wyników z dwóch modeli wyłącznie dlatego, że są korzystniejsze w różnych zakresach — to wymaga walidacji.

Jak czytać wynik kalkulatora Podmuch

Najpierw sprawdza się model i tryb geometrii. Air oznacza wybuch swobodny w powietrzu, surface oznacza model powierzchniowy. Potem patrzy się na Z, bo ono mówi, gdzie leży punkt względem zakresu danych. Następnie warto porównać Ps z Pr: duża różnica oznacza, że sama mapa nadciśnienia w wolnym polu może zaniżać obciążenie przegrody ustawionej prostopadle do fali.⁵

Dopiero po tym sensownie czyta się progi skutków. Próg 35 kPa dla zabudowy mieszkalnej, 70 kPa dla cięższych uszkodzeń przemysłowych czy niższe progi szyb są orientacyjne. Nie mówią, że każdy budynek w tej izobarze zachowa się identycznie. Mówią raczej, że w tym przedziale zaczyna się określona klasa uszkodzeń w typowych warunkach. Dla rzeczywistej konstrukcji potrzebna jest jej odpowiedź dynamiczna, nie tylko etykieta progu.⁵

Ostatnim krokiem jest spojrzenie na impuls i czas. Jeżeli wynik ma być użyty w ćwiczeniu z SDOF, to Pr i Ir są zwykle ważniejsze niż samo Ps, bo element konstrukcyjny jest obciążany falą odbitą. Jeżeli celem jest porównanie kilku modeli, trzeba zestawić nie tylko piki, ale też impulsy i czasy dodatniej fazy.

Walidacja: jak sprawdzić, czy kalkulator działa poprawnie

Każda implementacja K-B/CONWEP powinna dawać wyniki zgodne z tabelami referencyjnymi dla kilku punktów kontrolnych. Wartości referencyjne dla Z = 1,0 m/kg^(1/3) (wybuch swobodny, powietrze):

• Ps ≈ 490–510 kPa (ok. 70–74 psi)
• Pr ≈ 2000–2100 kPa (ok. 290–300 psi)
• Is ≈ 200–220 kPa·ms
• Ta (czas przybycia przy W = 1 kg) ≈ 0,5 ms

Dla Z = 5,0 m/kg^(1/3) (wybuch swobodny):

• Ps ≈ 14–16 kPa (ok. 2 psi)
• Pr ≈ 30–35 kPa
• Is ≈ 35–40 kPa·ms

Jeżeli kalkulator podaje wyniki mocno odbiegające od tych wartości, to sygnał do sprawdzenia jednostek współczynników, zakresu Z i trybu geometrii (swobodny vs. półsferyczny powierzchniowy). Lokalna implementacja ma zaszyte automatyczne testy jednostkowe weryfikujące te punkty — student może je uruchomić bezpośrednio z interfejsu trybu debug kalkulatora Podmuch.

Dobra praktyka: po ustawieniu nowej odległości lub masy zawsze sprawdzić, czy kalkulator wyświetla komunikat o zakresie ważności. Jeżeli Z leży na granicy zakresu K-B, wynik może być interpolowany lub wsparty modelem KG — co powinno być wyraźnie zaznaczone.⁵

Minimalna procedura bezpiecznego ćwiczenia obliczeniowego

Dobre ćwiczenie dla studentów nie wymaga żadnych materiałów wybuchowych. Wystarczy tabela syntetycznych punktów albo lokalny kalkulator. Dla kilku odległości należy policzyć Z, Ps, Pr, Is i Ir, a potem sprawdzić, jak zmienia się stosunek Pr/Ps i jak rośnie lub maleje impuls. Taki zestaw od razu pokazuje, dlaczego ten sam próg nadciśnienia nie wystarcza do opisu obciążenia konstrukcji.

Drugi wariant ćwiczenia polega na porównaniu dwóch uzysków przy tej samej odległości skalowanej. Studenci mogą dobrać odległości tak, aby Z było podobne, a następnie sprawdzić, że piki są zbliżone, lecz czas dodatniej fazy i impuls nie są identyczne w zwykłej skali czasu. To intuicyjnie tłumaczy, czemu większe eksplozje są nie tylko „mocniejsze”, ale też dłużej obciążają duże konstrukcje.

Trzeci wariant dotyczy walidacji. Należy wziąć trzy modele albo źródła: K-B/CONWEP, Kinney-Graham i jedną aproksymację z literatury przeglądowej. Dla tych samych wartości Z porównuje się Ps i Is, a następnie zaznacza przedziały, w których modele są blisko siebie i przedziały, w których różnica staje się zbyt duża na beztroskie użycie. Wynikiem ćwiczenia powinna być tabela niepewności modelowej, a nie „jedna prawdziwa liczba”.

Czwarty wariant, dla zaawansowanych, polega na zakłóceniu skalowania: student wybiera dwa ładunki o masach 10 kg i 10 000 kg, wyznacza odległości dające to samo Z = 3,0 m/kg^(1/3) i sprawdza, ile wynosi bezwzględna odległość R, bezwzględny czas przybycia i bezwzględna wartość impulsu dla obu przypadków. To ćwiczenie trwale zapada w pamięć, bo widać, jak ogromne odległości kryją się za skalowanymi parametrami dla dużych ładunków, i dlaczego skutki duży eksplozji są tak rozległe.

Dodatkowe materiały multimedialne

Powiązane kalkulatory i narzędzia

• Podmuch - liczy Ps, Pr, impuls, czas przybycia i orientacyjny przebieg P(t) z modeli K-B/CONWEP oraz Kinney-Graham.
• Równoważnik TNT - pokazuje, jak energia i sprawność sprzężenia przechodzą w efektywną masę TNT i odległość skalowaną.
• Pressure-Impulse - przenosi Pr i impuls na dydaktyczny model odpowiedzi ściany.
• Height-of-burst i Mach stem - pokazuje wpływ geometrii wybuchu powietrznego, odbicia i frontu Mach stem.
• Odporność płyty - liczy modalną odpowiedź prostokątnej płyty na zadany przebieg p(t).

Ćwiczenia praktyczne

Pierwsze ćwiczenie powinno polegać na zbudowaniu arkusza porównującego Ps i Is dla kilku wartości Z. Należy:

1. wybrać pięć wartości Z w zakresie ważności K-B,
2. odczytać albo policzyć Ps, Pr, Is i Ir,
3. obliczyć stosunki Pr/Ps oraz Ir/Is,
4. narysować wykres w skali logarytmicznej,
5. opisać, w których zakresach fala odbita najbardziej zmienia obraz obciążenia.

Drugie ćwiczenie powinno dotyczyć kształtu przebiegu ciśnienia w czasie. Należy:

1. przyjąć dane syntetyczne Ps, Is i td,
2. narysować prosty impuls trójkątny,
3. narysować zmodyfikowany przebieg Friedlandera,
4. porównać pola pod wykresami,
5. wyjaśnić, dlaczego dwa przebiegi o tym samym piku mogą dawać różną odpowiedź konstrukcji.

Trzecie ćwiczenie powinno być audytem modelu, a nie obliczeniem „zasięgu”. Należy:

1. dla jednego punktu policzyć wynik K-B/CONWEP i Kinney-Graham,
2. zapisać Z, zakres ważności i ewentualne ostrzeżenia,
3. sprawdzić osobno nadciśnienie i impuls,
4. wskazać, które różnice mają znaczenie dla szkła, a które dla ciężkiej płyty,
5. przygotować krótką notatkę: „czy ten wynik nadaje się do ćwiczenia dydaktycznego, do wstępnego projektu, czy tylko do jakościowej ilustracji?”.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Równoważnik trotylowy, nadciśnienie i impuls
• Fala odbita i obciążenie konstrukcji
• Pressure-Impulse i model SDOF: dlaczego konstrukcje nie reagują tylko na pik ciśnienia
• Mach stem i wysokość wybuchu: pomiary mikrofonowe, szybka kamera i symulacja
• Wysokość detonacji a fala podmuchowa