Streszczenie

Prędkość detonacji, brizancja i zdolność podmuchowa to trzy różne pojęcia, które w popularnych opisach bywają bezsensownie wrzucane do jednego worka jako „moc wybuchu”. Tymczasem VoD opisuje szybkość frontu detonacji w materiale, brizancja dotyczy lokalnego efektu kruszącego i gwałtowności wzrostu ciśnienia, a zdolność podmuchowa jest silniej związana z pracą fali w powietrzu i jej impulsem.1,2

Rozróżnienie jest potrzebne także w technice jądrowej. W soczewkach wybuchowych nie szuka się po prostu „najszybszego” materiału, lecz materiału o takich własnościach, które dadzą właściwą geometrię fali. Artykuł pokazuje, dlaczego liczba m/s nie zastępuje pełnego opisu działania materiału wybuchowego.1,3

Rozszerzenie tematu

Najpierw prędkość detonacji. To szybkość, z jaką front reakcji przemieszcza się przez materiał wybuchowy. Jest parametrem ważnym, bo mówi, jak szybko dany materiał potrafi wygenerować strefę wysokiego ciśnienia i jak będzie zachowywał się w geometrii ładunku. Ale sama wysoka VoD nie oznacza jeszcze, że materiał będzie najlepszy do każdego zadania.1

Co więcej, nawet sama liczba VoD nie jest wielkością absolutnie stałą. Polskie prace z AGH, WITU i KGHM Cuprum przypominają, że na zmierzoną prędkość detonacji wpływają średnica ładunku, gęstość, otoczka, rozdrobnienie kryształów, sposób mieszania składników, średnica otworu strzałowego, temperatura i charakter pobudzacza. To dlatego wartość katalogowa bywa tylko punktem startu, a nie końcem opisu materiału.5,6,7

Drugi termin to brizancja. W praktyce odnosi się on do zdolności materiału do bardzo gwałtownego, lokalnego kruszenia i rozrywania otaczającej materii. Materiał może mieć wysoką brizancję, bo daje stromy front ciśnienia i silny efekt lokalny, a jednocześnie nie być szczególnie dobrym źródłem długiego impulsu podmuchowego w powietrzu. To dlatego amunicja odłamkowo-burząca i amunicja o podwyższonej zdolności podmuchowej nie są tym samym problemem projektowym.2,4

Trzecia sprawa to zdolność podmuchowa. Tu bardziej niż sama gwałtowność liczy się to, jaki przebieg ma fala w powietrzu: jakie jest nadciśnienie, jak długo trwa dodatnia faza i jaki jest jej impuls. Artykuł WAT o improwizowanych ładunkach wybuchowych pokazuje to bardzo dobrze: autorzy badają rzeczywiste charakterystyki fal i na tej podstawie wyciągają wnioski o skutkach, zamiast opierać się tylko na jednej liczbie związanej z materiałem. Ten sam problem wraca w artykule o równoważniku trotylowym, nadciśnieniu i impulsie: energia całkowita jest ważna, ale nie zastępuje przebiegu fali.2

W praktyce oznacza to, że materiał o wyższej prędkości detonacji nie musi być „lepszy” do niszczenia ścian, bunkrów albo lekkich konstrukcji, jeśli krótsza i ostrzejsza fala daje mniejszy impuls całkowity niż dłuższa, wolniejsza fala innego układu. Właśnie z tego powodu amunicja termobaryczna i FAE bywa oceniana jako szczególnie skuteczna w pomieszczeniach: nie dlatego, że ma rekordową VoD, lecz dlatego, że daje dłużej trwającą falę podmuchową o bardzo nieprzyjemnym przebiegu dla ludzi i lekkich konstrukcji.2,4

To rozróżnienie jest kluczowe również dla implozji jądrowej. W Fat Manie i w soczewkach wybuchowych nie chodziło o użycie „najsilniejszego materiału”, tylko o precyzyjne zestawienie materiału szybkiego i wolnego, aby ukształtować falę zbieżną. Z punktu widzenia takiego układu zbyt szybka lub zbyt brizantna mieszanina w złym miejscu może pogorszyć wynik, a nie go poprawić. To zarazem tłumaczy, dlaczego TNT jako materiał odlewany pozostawał użyteczny mimo parametrów gorszych niż czyste wysokoenergetyczne kruszywa.3

Z perspektywy pomiarowej ważne jest jeszcze jedno: inne metody badawcze pokazują nieco inne strony tego samego zjawiska. Artykuł z WITU porównuje metodę czujników zwarciowych z ciągłym pomiarem detonacji i podkreśla, że pomiar ciągły pozwala zobaczyć nie tylko średnią prędkość, ale też rozwijanie się, stabilizowanie i wygaszanie detonacji w kolejnych częściach ładunku. Z kolei tekst z KGHM Cuprum zestawia urządzenia MicroTrap i Explomet-Fo-2000, pokazując, że w praktyce polowej najbardziej interesuje projektanta nie laboratoryjna liczba z normy, lecz wartość eksploatacyjna w rzeczywistym otworze strzałowym.6,7

Najkrótszy wniosek jest taki: VoD, brizancja i podmuch opisują trzy różne twarze wybuchu. Mylenie ich prowadzi do błędnych intuicji, a w technice wojskowej i jądrowej właśnie te niuanse decydują, czy ładunek naprawdę robi to, czego oczekuje projektant.

Prędkość detonacji (VoD) — parametr kinetyczny

Prędkość detonacji to fundamentalny parametr kinetyczny opisujący, jak szybko front reakcji detonacyjnej przesuwa się przez materiał wybuchowy. Jego fizycznym odpowiednikiem w teorii Chapman-Jouguet jest prędkość D_CJ — prędkość stacjonarnego frontu detonacyjnego wyznaczona z warunków Rankine-Hugoniot dla fali detonacyjnej przy założeniu, że za frontem strumień gazów porusza się z prędkością lokalnej dźwięku (warunek CJ).

VoD określa kilka ważnych właściwości operacyjnych materiału:

Czas dotarcia fali: W ładunku o długości 0,5 m materiał o VoD = 8 000 m/s zainicjuje detonację całego ładunku w czasie t = 0,5 / 8 000 = 62,5 µs. Ta liczba jest ważna wszędzie tam, gdzie czas synchronizacji detonacji ma znaczenie — np. w soczewkach wybuchowych, gdzie przesunięcia fazowe frontu wybuchu są mierzone w jednostkach nanosekund.

Ciśnienie Chapman-Jouguet: P_CJ ≈ ρ₀ × D² / 4, gdzie ρ₀ jest gęstością materiału wybuchowego. Dla TNT (ρ₀ ≈ 1 630 kg/m³, D ≈ 6 930 m/s): P_CJ ≈ 1 630 × (6 930)² / 4 ≈ 19,5 GPa. Dla RDX (ρ₀ ≈ 1 760 kg/m³, D ≈ 8 700 m/s): P_CJ ≈ 1 760 × (8 700)² / 4 ≈ 33,3 GPa. Wyższe ciśnienie CJ oznacza silniejsze bezpośrednie działanie mechaniczne w pobliżu ładunku.1,3

Tabela VoD wybranych materiałów wybuchowych (przy podanych gęstościach):

Materiał Gęstość [g/cm³] VoD [m/s] P_CJ [GPa]
ANFO (94/6) 0,85 4 200 1,5
Emulsja AN 1,15 5 500 3,4
TNT 1,63 6 930 19,5
Composition B (RDX/TNT) 1,72 7 980 27,4
PETN 1,77 8 300 30,5
RDX 1,76 8 700 33,3
HMX 1,90 9 100 39,4
CL-20 2,04 9 700 48,0

Wzrost VoD i P_CJ idzie w parze z gęstością i zawartością energii na jednostkę objętości. Ale sam VoD nie mówi nic o tym, jak ta energia jest rozdzielona między kruszenie, podmuch i fale sejsmiczne — do tego potrzebne są inne parametry.1,2

Brizancja — parametr efektu kruszącego

Brizancja (od fr. briser — kruszyć) opisuje zdolność materiału wybuchowego do miejscowego, gwałtownego kruszenia i rozerwania otaczającej materii bezpośrednio w strefie kontaktu. Jest ona miarą gwałtowności działania, nie totalnej energii — dwa materiały mogą mieć podobną energię wybuchu, ale zupełnie inną brizancję.

Z punktu widzenia fizyki brizancja wynika ze stromości frontu fali ciśnienia. Im wyższy jest gradient ciśnienia (dP/dt) na czole fali, tym bardziej gwałtowny jest efekt miejscowy. Materiały o wysokiej VoD i dużym P_CJ mają zazwyczaj wyższą brizancję — ale nie jest to ścisła korelacja.

Testy brizancji

Test Trauzla (Trauzl lead block test): Cylindryczny blok ołowiu o masie 200 g (lub 500 g) z wydrążoną jamą jest detonowany przez standardową ilość badanego materiału (10 g). Ekspansja jamy, mierzona w cm³, jest liczbą Trauzla — empiryczną miarą zdolności do kruszenia. Wyższe wartości = wyższa brizancja. Przykłady: TNT ≈ 300–330 cm³; PETN ≈ 450–480 cm³; RDX ≈ 480–500 cm³; HMX ≈ 490–520 cm³; ANFO ≈ 80–100 cm³.

Test BATES (ballistic air-traffic explosive sand) i test piaskowy: Mierzą kompresję piasku lub innego standardowego medium przez falę detonacyjną — stosowane w klasyfikacji materiałów górniczych.

Heusinger von Waldegg test: Mierzy efekt kruszenia płyty stalowej ładunkiem w kontakcie — bardziej zbliżony do warunku bojowego niż test w bloku ołowiu.

Udział brizancji w projektowaniu ładunków:

  • Amunicja kumulacyjna (shaped charge, HEAT): Wysoka brizancja i VoD potrzebna do formowania strumienia kumulacyjnego z wykładziny (liner), który penetruje pancerz.
  • Amunicja burząco-odłamkowa (HE-Frag): Brizancja ważna, ale zbalansowana z wymaganiem dobrego impulsu odłamkowego.
  • Ładunki górnicze: Brizancja często nieistotna lub wręcz niepożądana (zbyt kruszący materiał niszczy nadmiar skały, zwiększa koszty); ważniejsza jest zdolność rozrzutu i podmuch.
  • Ładunki do obróbki metali (explosive forming, explosive cladding): Wymagana jest określona wartość brizancji odpowiednia do sposobu odkształcenia plastycznego blach i rur bez ich rozrywania. Zbyt wysoka brizancja prowadzi do pęknięć zamiast plastycznej deformacji.
  • Ładunki do robót rozbiórkowych (demolition): Często korzysta się z materiałów o pośredniej brizancji, ustawione precyzyjnie na kluczowych elementach konstrukcji. Celem nie jest maksymalne nadciśnienie, lecz selektywne osłabienie struktury nośnej.2,4

Teoria energii detonacji i podział na frakcje

Całkowita energia chemiczna uwolnionej przy detonacji (energia eksplozji Q_ex, wyznaczana kalorymetrycznie) jest dzielona między różne mechanizmy:

  1. Energia fali podmuchowej w powietrzu — przenosi się daleko od ładunku, maleje jak R⁻³ (na bliskie pole) i R⁻¹ (na dalekie pole Mach/Friedlander)
  2. Energia kinetyczna odłamków — przekazana ściance koszulki lub odłamkom płaszcza
  3. Energia fali sejsmicznej — pochłaniana przez grunt (przy wybuchu przyziemnym lub podziemnym)
  4. Energia promieniowania cieplnego — ważna przy wybuchu jądrowym, marginalna dla MW konwencjonalnych
  5. Energia dźwięku i fal sprężystych w zabudowie — straty przez drgania struktur

Dla MW konwencjonalnych podział jest orientacyjnie: ~30–50% energii trafia w falę podmuchową (zależy od geometrii i konfinementu), ~30–50% w kinetykę odłamków lub deformację plastyczną, reszta w fale sejsmiczne i ciepło.

Dla materiałów FAE (fuel-air explosives) i termobarycznych, które detonują w powietrzu z fazowym wybuchem, udział energii w fali podmuchowej jest znacząco wyższy — do 70–80% — kosztem brizancji i kinetyki odłamków. To wyjaśnia, dlaczego amunicja termobaryczna jest oceniana jako szczególnie skuteczna przy ataku na ludzi w okopach, bunkrach i zamkniętych przestrzeniach, gdzie długa fala podmuchowa ma najdotkliwszy efekt.4

Fala podmuchowa: parametry i opis Friedlandera

Fala podmuchowa w powietrzu jest opisywana przez kilka kluczowych parametrów:

  • Nadciśnienie szczytowe P_s [kPa lub PSI]: Chwilowe maksymalne nadciśnienie powyżej wartości atmosferycznej na czole fali
  • Czas trwania fazy dodatniej t_d [ms]: Czas od momentu dotarcia fali do czasu, gdy ciśnienie spada z powrotem do wartości atmosferycznej
  • Impuls fazy dodatniej i_s [kPa·ms]: Całka nadciśnienia po czasie w fazie dodatniej; i_s = ∫P(t)dt dla 0 < t < t_d
  • Faza ujemna (underpressure): Po fazie dodatniej ciśnienie spada chwilowo poniżej atmosferycznego — efekt „zasysania"; może mieć znaczenie dla lekkich konstrukcji wrażliwych na siły odrywające

Model Friedlandera opisuje przybliżony kształt ciśnienia w czasie po detonacji:

P(t) = P_s × (1 − t/t_d) × exp(−α × t/t_d)

gdzie α jest stałą kształtu (decay parameter). Model jest użyteczny do analitycznych obliczeń odpowiedzi konstrukcji, choć realna fala jest bardziej skomplikowana (efekty wielokrotnych odbić, Mach stem itp.).2

Skalowanie Hopkinsona-Cranza: Nadciśnienie, impuls i czas trwania fali są funkcją skalowanej odległości Z:

Z = R / W^(1/3) [m/kg^(1/3)]

gdzie R to odległość od ładunku, W to masa ekwiwalentna TNT. Z tym samym Z uzyska się te same wartości P_s i i_s, niezależnie od rzeczywistych wartości R i W (w granicach skalowalności). To prawo podobieństwa umożliwia przenoszenie wyników testów z małych ładunków na duże.2

Równoważnik trotylowy i jego ograniczenia

Równoważnik trotylowy (TNT equivalent) to powszechnie stosowana, ale głęboko niedoskonała miara porównawcza dla materiałów wybuchowych i ładunków jądrowych.

Definicja: Masa TNT produkująca taki sam efekt (np. to samo nadciśnienie na danej odległości) jak badany ładunek o masie W_x:

TNT_eq = W_x × (Q_x / Q_TNT)

gdzie Q_x i Q_TNT to ciepła wybuchu. Dla masy nitrogliceryny Q_NG ≈ 6 200 J/g, Q_TNT ≈ 4 600 J/g, więc TNT_eq(NG) ≈ 6 200/4 600 ≈ 1,35.

Problemy z TNT_eq:

  1. Różny podział energii: TNT_eq obliczony na podstawie energii chemicznej nie uwzględnia, że różne materiały inaczej dzielą energię między podmuch, odłamki i fale sejsmiczne. Materiał termobaryczny może mieć TNT_eq = 1,5 na podstawie energii chemicznej, ale TNT_eq ≈ 2–3 pod kątem efektu podmuchowego w zamkniętej przestrzeni.

  2. Zależność od odległości: Kształt krzywej P(R) różni się dla różnych materiałów przy małych R (bliska strefa) i dużych R (dalekie pole). TNT_eq może być „2" na 10 m od ładunku i „0,8" na 100 m od tego samego ładunku.

  3. Detonacja nieidealną: Materiały ANFO i emulsje w warunkach polowych mają VoD znacznie poniżej D_CJ, co obniża efektywne TNT_eq w porównaniu z wartościami laboratoryjnymi.

  4. Efekty termiczne: Przy ładunkach jądrowych TNT_eq jest wyrażany zwykle przez energię ogólną, ale efekty termiczne (promieniowanie cieplne) i EMP (electromagnetic pulse) nie mają odpowiednika w konwencjonalnym TNT.

Mimo tych ograniczeń TNT_eq pozostaje powszechnym skrótem konwencjonalnym w literaturze i normach — ale należy go czytać ze świadomością, do jakich celów był wyznaczony.2,3

Brizancja a zdolność podmuchowa: konkretne przykłady

Zestawienie brizancji i impulsu podmuchowego dla wybranych materiałów pozwala zobaczyć, że te dwa parametry nie są ze sobą ściśle skorelowane.

Materiał VoD [m/s] Liczba Trauzla [cm³] Impuls podmuchowy (orientacyjny)
ANFO 4 200 80–100 Wysoki na dalekie pole
TNT 6 930 300–330 Średni (referencja)
RDX 8 700 480–500 Niższy na dalekie pole niż ANFO
PETN 8 300 450–480 Podobny do RDX
FAE (propylene oxide) 1 500–2 000 ~10–30 Bardzo wysoki na dalekie pole
Termobar. TBG-7 1 800–2 200 niska Wyjątkowo wysoki w ograniczonej przestrzeni

Obserwacja: ANFO ma dramatycznie niższy VoD i brizancję niż RDX, ale jego impuls podmuchowy na dalekich odległościach (duże pole) bywa wyższy dzięki długiej fazie dodatniej fali (dłuższy czas trwania wybuchu przy niższym ciśnieniu szczytowym = dłuższa fala). To wyjaśnia, dlaczego ANFO jest standardowym materiałem w górnictwie odkrywkowym (kruszenie skały na duże odległości od otworu strzałowego), podczas gdy HMX jest używany w kumulatywnych głowicach i soczewkach wybuchowych.1,4

Detonacja nieidealna i efekt krytycznej średnicy

Pojęcie detonacji nieidealnej jest kluczowe dla zrozumienia zachowania materiałów wybuchowych poza warunkami laboratoryjnymi. Teoria Chapman-Jouguet zakłada płaski, stacjonarny front detonacyjny o zerowej grubości strefy reakcji. W praktyce strefa reakcji ma skończoną grubość — rzędu 0,1–10 mm dla różnych materiałów — i boczne fale rarefakcyjne nieustannie penetrują tę strefę.

Efekt bocznych fal rarefakcyjnych: Gdy ładunek ma małą średnicę, fale rarefakcji ze ścian docierają do osi symetrii w czasie porównywalnym z czasem przejścia materiału przez strefę reakcji. Skutkuje to lokalnym obniżeniem ciśnienia i temperatury w strefie reakcji — część reakcji nie zdąża się zakończyć zanim nie nastąpi ekspansja. Prowadzi to do obniżenia zmierzonego VoD poniżej D_CJ. Im mniejsza średnica, tym większy jest relatywny wpływ tych bocznych strat energii.

Krytyczna średnica d_crit: Poniżej pewnej granicznej średnicy ładunku detonacja nie może się utrzymać — gasnie lub przechodzi w deflagrację. Ta granica zależy od materiału:

  • TNT: d_crit ≈ 8–10 mm w otwartej przestrzeni
  • PETN: d_crit ≈ 1–2 mm (mały, bo szybka reakcja)
  • ANFO: d_crit ≈ 100–200 mm (duży, materiał niejednorodny, wolna reakcja)
  • Emulsja AN: d_crit ≈ 50–100 mm (zależy od formulacji)

Znajomość d_crit jest niezbędna przy projektowaniu ładunków górniczych (muszą mieć średnicę otworu > d_crit) i ładunków wojskowych (ogranicza miniaturyzację). Materiały o małym d_crit (PETN, RDX, HMX) mogą detonować w bardzo cienkich filmach — co czyni je przydatnymi w detonatorach inicjujących i linach detonujących (detonating cord).1,2

Energia Gurneya i akceleracja metalu

Jednym z najważniejszych parametrów praktycznych materiałów wybuchowych — zwłaszcza w zastosowaniach wojskowych — jest energia Gurneya (Gurney energy, E_G), opisująca zdolność materiału do akceleracji metalowych płyt, koszulek i fragmentów.

Formuła Gurneya dla cylindrycznej koszulki metalowej o masie M i masie MW W wewnątrz:

v_metal = √(2E_G) × [M/W + 1/2]^(−1/2)

gdzie √(2E_G) jest prędkością Gurneya (Gurney velocity) wyrażaną w km/s — empiryczną stałą materiałową.

Przykładowe wartości prędkości Gurneya:

Materiał √(2E_G) [km/s]
ANFO 2,0
TNT 2,44
Composition B 2,70
PETN 2,93
RDX 2,93
HMX 3,05
CL-20 3,2

Prędkość Gurneya koreluje z P_CJ i VoD, ale nie jest od nich liniowo zależna. To oddzielna miara właściwości materiału — przydatna bezpośrednio przy projektowaniu ładunków kumulacyjnych, pocisków odłamkowych i penetratorów kinetycznych.1,2

Odbita fala podmuchowa i efekt Macha

Gdy fala podmuchowa natrafi na sztywną powierzchnię (podłoże, ścianę budynku), następuje odbicie. Zachowanie fali odbitej zależy krytycznie od kąta padania θ fali względem przeszkody.

Normalne odbicie (θ = 0°, fala pada prostopadle na powierzchnię): Ciśnienie odbite P_r jest znacznie wyższe niż nadciśnienie padające P_s:

P_r / P_s = 2 × (7P₀ + 4P_s) / (7P₀ + P_s)

gdzie P₀ = 101,3 kPa to ciśnienie atmosferyczne. Dla P_s = 100 kPa: P_r/P_s ≈ 2,08 × (7×101,3 + 400)/(7×101,3 + 100) ≈ 2,5. Dla bardzo dużych P_s >> P₀: P_r/P_s → 8. Budynki i konstrukcje trafiane pod kątem bliskim 0° otrzymują odbite ciśnienie wielokrotnie wyższe niż na otwartej przestrzeni.

Mach stem (θ przekroczenie kąta krytycznego): Dla kątów padania powyżej ~45°–50° (zależy od P_s) fala padająca i odbita nie mogą koegzystować jako oddzielne fale i łączą się w jeden Mach front — pionową falę o wyższym ciśnieniu niż każda z fal składowych osobno. Mach stem rośnie w miarę oddalania się od ładunku. W strefie Macha (między podłożem a wysokością, na której Mach stem jeszcze nie sięga) nadciśnienie i impuls są wyższe niż w przypadku padania bezpośredniego. Jest to ważne dla planowania wybuchu powietrznego (airburst) — detonacja powyżej gruntu może dawać silniejsze uszkodzenia naziemne niż detonacja na powierzchni, ze względu na optymalną warstwę Macha.

Efekt Macha jest standardowo uwzględniany w obliczeniach systemem CONWEP i normach UFC 3-340-02 (DoD USA) dla wszystkich scenariuszy wybuchów przyziemnych i powietrznych.2

Dla projektantów budynków odpornych na wybuchy (blast-resistant buildings) najważniejsze jest pytanie: w jaki sposób wybrać jako „wejście" do modelu budynku ciśnienie, impuls lub ich kombinację? Odpowiedź zależy od scenariusza i charakterystyk konstrukcji. Ciężkie masywne ściany betonowe są raczej „wrażliwe na impuls" (impulse-sensitive): reagują na całkę sił, nie na chwilowe maksimum. Lekkie metalowe elewacje i okna są bardziej „ciśnieniowrażliwe" (pressure-sensitive): pękają, gdy nadciśnienie chwilowe przekroczy próg wytrzymałości materiału, niezależnie od długości fazy. Wykresy P-I (Pressure-Impulse diagram) są standardowym narzędziem łączącym obie wrażliwości w jednej płaszczyźnie, pozwalając projektantowi ocenić margines bezpieczeństwa dla danego scenariusza zagrożenia.2,4

Zastosowanie w ładunkach kumulacyjnych i soczewkach

W technice militarnej kluczowym zastosowaniem zrozumienia zależności VoD–brizancja–podmuch są ładunki kumulacyjne (shaped charges) i soczewki wybuchowe (explosive lenses).

Ładunki kumulacyjne: Głównym efektem jest formowanie metalowego strumienia z wykładziny (metallic liner) przez ciśnienie gazów detonacyjnych. Prędkość strumienia koreluje z P_CJ (i pośrednio z VoD): dla RDX prędkości wierzchołka strumienia sięgają 8–10 km/s. W tym zastosowaniu wysoki VoD i P_CJ jest pożądany, brizancja (kruszenie bezpośrednie) jest mniej istotna niż zdolność do przyspieszenia liniowej wykładziny, a podmuch jest efektem ubocznym.

Soczewki wybuchowe w urządzeniach implozyjnych: Jak opisuje artykuł o soczewkach Baratol i Kompozyt B, soczewka wymaga dwóch materiałów o bardzo różnych VoD. Materiał „szybki" (fast explosive, np. Kompozyt B: VoD ≈ 7 980 m/s) i „wolny" (slow explosive, np. Baratol: VoD ≈ 4 900 m/s) są ukształtowane tak, że różnica opóźnień czasowych wynikająca z różnych VoD koryguje geometrię fali detonacyjnej z kulisty-rozbieżny na zbieżny (converging spherical wave). W tym zastosowaniu absolutna wartość VoD jest mniej ważna niż różnica między materiałami — i niż jednorodność VoD w partii produkcyjnej, której tolerancje sięgają ±10–20 m/s.3

Skala Hopkinsona-Cranza: praktyczna analiza

Prawo podobieństwa Hopkinsona-Cranza umożliwia szybkie obliczenie parametrów fali podmuchowej dla dowolnego ładunku TNT na dowolnej odległości, jeśli znana jest krzywa P_s(Z) i i_s(Z).

Standardowe dane dla ładunku powierzchniowego (hemispherical surface burst, HSB) z modelu Kingery-Bulmash używanego w CONWEP:

Z [m/kg^(1/3)] P_s [kPa] i_s [kPa·ms/kg^(1/3)] t_d [ms/kg^(1/3)]
0,5 10 000 800 0,12
1,0 2 400 500 0,22
2,0 500 200 0,50
3,0 200 110 0,85
5,0 70 55 1,50
10,0 17 22 3,20
20,0 4,5 8 7,00

Aby obliczyć absolutne wartości dla konkretnego ładunku W (kg TNT) na odległości R (m), wyznacza się Z = R/W^(1/3), odczytuje P_s(Z) i i_s(Z)/kg^(1/3), a następnie:

i_s [kPa·ms] = i_s(Z)/kg^(1/3) × W^(1/3)

Przykład: W = 100 kg TNT, R = 50 m, Z = 50 / 100^(1/3) = 50 / 4,64 = 10,8 m/kg^(1/3). Z tabeli P_s ≈ 15 kPa, i_s ≈ 20 kPa·ms/kg^(1/3). Absolutny impuls: i_s = 20 × 100^(1/3) = 20 × 4,64 ≈ 93 kPa·ms.

Takie obliczenia są podstawą projektowania ochrony budynków, zatwierdzania stref bezpieczeństwa i analizy zagrożeń wybuchami.2

Warto podkreślić ograniczenia tabeli Kingery-Bulmash: dane te zostały wyznaczone empirycznie z setek testów z TNT lub równoważnikami TNT w otwartej atmosferze. Wyniki mogą się znacząco różnić dla: (a) materiałów o innym kształcie P-i-t niż TNT (zwłaszcza materiałów termobarycznych i ANFO), (b) wybuchów w zamkniętych przestrzeniach (wzmocnienie przez wielokrotne odbicia), (c) wybuchów w terenie zabudowanym (canalisation fali między budynkami). Dla tych przypadków zalecane są pełne obliczenia CFD z kodami takimi jak LS-DYNA ALE lub Air3D, kalibrowane odrębnie.

Trzy przykłady numeryczne

Przykład 1: Porównanie brizancji TNT i RDX

Dla TNT: VoD = 6 930 m/s, ρ₀ = 1 630 kg/m³.
P_CJ(TNT) ≈ 1 630 × (6 930)² / 4 ≈ 1 630 × 4,80 × 10⁷ / 4 ≈ 19,6 GPa

Dla RDX: VoD = 8 700 m/s, ρ₀ = 1 760 kg/m³.
P_CJ(RDX) ≈ 1 760 × (8 700)² / 4 ≈ 1 760 × 7,57 × 10⁷ / 4 ≈ 33,3 GPa

Ciśnienie CJ jest 1,70 × wyższe dla RDX niż dla TNT. Odpowiada to przybliżonemu stosunkowi liczb Trauzla: RDX ≈ 490 cm³ / TNT ≈ 315 cm³ = 1,56. Korelacja nie jest idealna ze względu na różne właściwości pola ciśnienia za frontem, ale pokazuje kierunek.1

Przykład 2: Odległość bezpieczna a efekt brizancji

Dla personelu nieosłoniętego przyjmuje się próg uszkodzeń słuchu przy nadciśnieniu P_s ≈ 6 kPa i ryzyko poważnych obrażeń wewnętrznych przy P_s ≈ 100 kPa. Dla ładunku 5 kg RDX (TNT_eq ≈ 6 kg TNT):

Odległość dla P_s = 6 kPa: Z(6 kPa) ≈ 11 m/kg^(1/3) → R = 11 × 6^(1/3) = 11 × 1,82 ≈ 20 m
Odległość dla P_s = 100 kPa: Z(100 kPa) ≈ 4 m/kg^(1/3) → R = 4 × 1,82 ≈ 7,3 m

Wniosek: To działa niezależnie od brizancji. Nawet materiał o bardzo wysokim P_CJ, gdy tylko zdetonuje w otwartej przestrzeni, tworzy falę podmuchową skalowaną tymi samymi prawami Hopkinsona-Cranza.

Przykład 3: Wpływ geometrii ładunku na podmuch vs. kruszenie

Ładunek 1 kg TNT:

  • Geometria kulista w powietrzu (spherical airburst): Maksymalizuje podmuch; energia promieniuje równomiernie, brak strat na podłoże. Efektywny w atakach na obszar.
  • Ładunek w kontakcie ze ścianą betonową: Energia jest częściowo pochłaniana przez beton; lokalne ciśnienia kontaktowe są ogromne (GPa), ale fala w powietrzu jest słabsza. Efektywny w kruszeniu.
  • Ładunek kumulacyjny (skoncentrowany strumień): Cała energia jest skupiona w wąskim strumieniu; brizancja i VoD decydują o głębokości penetracji (płyt pancernych), ale fala podmuchowa jest minimalna.

Te trzy konfiguracje z tym samym 1 kg TNT dają radykalnie różne efekty — co pokazuje, że VoD i masa ładunku to tylko część informacji potrzebnej do oceny zagrożenia.2,4

Historia pojęć: Berthelot, Chapman, Jouguet i Gurney

Współczesna teoria detonacji ma długą historię, a rozróżnienie między VoD, brizancją i podmuchem ewoluowało przez ponad sto lat badań.

Marcelin Berthelot (1827–1907), francuski chemik, jako jeden z pierwszych opisał detonację jako zjawisko fizyko-chemiczne odmienne od spalania. Berthelot i Vieille zmierzyli prędkości fali detonacyjnej w gazowych mieszaninach w latach 80. XIX wieku, a Berthelot opracował wczesne kalorymetryczne metody wyznaczania ciepła wybuchu. Jego praca nad „pracą wybuchu" (travail de l'explosion) była prototypem współczesnego pojęcia energii detonacji.

David Leonard Chapman (1869–1958) i Jacques Charles Émile Jouguet (1871–1943) niezależnie (Chapman 1899, Jouguet 1905–1906) sformułowali warunek pozwalający wyznaczyć prędkość detonacji stacjonarnej. Warunek CJ stanowi, że w stacjonarnej detonacji lokalny strumień gazów na czole fali jest właśnie dźwiękowy (liczba Macha M=1 w układzie związanym z frontem). To warunek termodynamicznie stabilny, który pozwala zamknąć równania bilansu masy, pędu i energii (Rankine-Hugoniot) i wyznaczyć unikalny punkt (p_CJ, v_CJ) na krzywej Hugoniot.

Geoffrey I. Taylor (1886–1975, Wielka Brytania) i jego teoria fali odprężeniowej (Taylor wave) za frontem CJ była kluczowa dla zrozumienia, jak ciśnienie spada za frontem detonacyjnym. Taylor współpracował nad Projektem Manhattan i jego obliczenia dotyczące dynamiki wybuchów konwencjonalnych były bezpośrednio używane przy projektowaniu ładunków implozyjnych.

Ronald Wilfred Gurney (1898–1953, Los Alamos) opracował w 1943 roku proste formuły opisujące akcelerację metalowych koszulek przez materiały wybuchowe, bez konieczności pełnego numerycznego rozwiązania równań hydrodynamiki. Formuły Gurneya do dziś są używane jako szybkie narzędzie szacunkowe w inżynierii wojskowej.1,3

Symulacje numeryczne: LS-DYNA i Autodyn

Współczesnym narzędziem do symulacji wszystkich trzech zjawisk — propagacji detonacji, efektów brizancji i fali podmuchowej — są kody FEM/FVM klasy LS-DYNA, ANSYS Autodyn i CTH (hydrocode SANL).

LS-DYNA jest komercyjnym kodem Explicit Dynamic FEM, szeroko stosowanym w środowisku naukowym i przemysłowym. Dla symulacji detonacji i efektów wybuchów kluczowe są:

  • Model materiału: karta *MAT_HIGH_EXPLOSIVE_BURN + równanie stanu *EOS_JWL z parametrami A, B, R₁, R₂, ω, E₀ dla każdego materiału wybuchowego
  • Inicjacja detonacji: karta *INITIAL_DETONATION lub *DETONATION_POINT definiująca miejsce i czas inicjacji
  • Modelowanie powietrza: Euler lub ALE (Arbitrary Lagrangian-Eulerian) z modelem powietrza jako gazu doskonałego

Schemat jawny (explicit) LS-DYNA wymaga kroku czasowego Δt < L_min / c, gdzie L_min to minimalny wymiar elementu i c prędkość dźwięku w materiale. Dla elementu 1 mm w stali (c ≈ 5 100 m/s) Δt < 0,2 µs — co dla symulacji kilku milisekund wybuchu daje miliony kroków czasowych.

Autodyn (ANSYS) jest popularną alternatywą, szczególnie dla modelowania zderzeń, penetracji i efektów wybuchowych w materiałach polidomenowych. Posiada wbudowane biblioteki parametrów JWL dla ponad 50 powszechnie stosowanych materiałów wybuchowych.

CTH (Sandia National Laboratories) jest kodem rządowym US (niekomercyjnym), specjalnie opracowanym do modelowania zjawisk podmuchowych w dużych domenach obliczeniowych, w tym atmosferze. Stosowany do symulacji nuklearnych testów powietrznych i zjawisk podmuchowych przy wybuchach przyziemnych.

W Polsce WAT dysponuje licencjami LS-DYNA i Autodyn, a wyniki symulacji tych kodów są porównywane z wynikami pomiarów terenowych w publikacjach biuletynu WAT.4,5

Znaczenie pojęć dla nieproliferacji i techniki jądrowej

Dla doktoranta kluczowe jest też powiązanie tych trzech parametrów z techniką jądrową:

VoD i soczewki implozyjne: Jak wyjaśniono wcześniej, różnica VoD między szybkim i wolnym materiałem w soczewce wybuchowej jest podstawowym parametrem projektowym. Producenci soczewek stosują materiały o ściśle kontrolowanym VoD ±5–20 m/s od specyfikacji, bo każde odchylenie geometryczne fali detonacyjnej przekłada się na asymetrię implozji i obniżenie wydajności ładunku.

Brizancja a inicjacja materiałów rozszczepialnych: Inicjator implosyjny musi dostarczyć wystarczający impuls do homogenicznego skompresowania rdzenia rozszczepialnego. Zbyt mała brizancja → niepełna kompresja → brak wymaganej gęstości krytycznej dla samonapędzającej reakcji łańcuchowej. Zbyt wysoka brizancja w złym miejscu → niejednorodne ciśnienie → predetonacja z niską wydajnością.

Podmuch jądrowy vs. konwencjonalny: Fala podmuchowa wybuchu jądrowego jest opisywana tymi samymi parametrami (P_s, i_s, t_d) i skaluje się wg tych samych praw Hopkinsona-Cranza, ale skala energii jest radykalnie inna. Dla wybuchu 15 kt TNT_eq (bomba Little Boy): na odległości R = 1 km Z = 1 000 / 15 000^(1/3) ≈ 1 000 / 24,7 ≈ 40,5 m/kg^(1/3). Z tabeli P_s ≈ 2,5 kPa — ale dla eksplozji jądrowej dochodzą czynniki termiczne (promieniowanie cieplne) i radiologiczne (wstępne promieniowanie jądrowe), które przy tej odległości są dominującymi czynnikami szkodliwymi.3

Perspektywa polska: zastosowania akademickie i przemysłowe

W Polsce materiałami wybuchowymi i ich parametrami zajmuje się kilka ośrodków:

WAT (Wojskowa Akademia Techniczna): Prowadzi badania nad charakterystykami fal podmuchowych, narzędziami numerycznymi do symulacji (LS-DYNA, ANSYS Autodyn) i oceną skutków wybuchów IED. Artykuł Steckiewicza i Trzcińskiego o falach podmuchowych z IED (Biuletyn WAT 2009) jest przykładem łączenia pomiaru VoD z analizą rzeczywistych fal podmuchowych.

WITU (Wojskowy Instytut Techniczny Uzbrojenia): Badania nad materiałami wybuchowymi stosowanymi w amunicji wojskowej, pomiary VoD metodami elektrycznymi i optycznymi, testy brizancji.

AGH Kraków i IGO: Skupiają się na przemysłowych materiałach strzałowych (emulsje, ANFO, dynamity), pomiaru VoD in situ i optymalizacji efektywności strzelania w górnictwie odkrywkowym i podziemnym.

KGHM Cuprum: W kopalniach KGHM Polska Miedź corocznie detonowane są dziesiątki tysięcy ładunków. Optymalizacja VoD i wzorca detonacji jest tu bezpośrednio związana z ekonomiką wydobycia i bezpieczeństwem robotników.5,6,7

Pytania otwarte

Czy prędkość detonacji nieidealnych materiałów da się przewidywać z pierwszych zasad? Dla jednolitych materiałów krystalicznych (RDX, HMX) obliczenia termochemiczne (CHEETAH, EXPLO5) są dość dokładne. Dla heterogenicznych emulsji i ANFO modele homo-geniczne zawodzą: efekt średnicy, hot spots i impedancja składników wymagają podejść mezoskopowych.

Jak nowoczesne kruszące materiały przyszłości (CL-20, ONC, HNIW) zmienią optymalne geometrie ładunków kumulacyjnych? CL-20 (HNIW) ma VoD ≈ 9 700 m/s i P_CJ ≈ 48 GPa — parametry znacznie wyższe niż HMX. Ale jego wrażliwość mechaniczna ogranicza zastosowania. ONC (octanitrocubane, 1,3,5,7-tetranitro-cubane) jest teoretycznie najsilniejszym dostępnym materiałem — ale synteza w skali kilogramowej pozostaje wyzwaniem.

Jak brizancja koreluje z efektami biologicznymi? Dla pełnej oceny zagrożeń militarnych i planowania ochrony ważna jest biofizyka oddziaływania fali podmuchowej z organizmem ludzkim. To aktywna dziedzina badań medycyny wojskowej.

Podsumowanie dydaktyczne

Dla doktoranta pracującego w obszarze materiałoznawstwa energetycznego, inżynierii wybuchów lub oceny zagrożeń kluczowe punkty z tego artykułu są następujące:

  1. VoD opisuje kinetykę frontu detonacyjnego, a pośrednio — ciśnienie CJ i czas synchronizacji. Jest niezbędny do projektowania geometrii ładunku, ale nie wystarcza do pełnego opisu efektów zewnętrznych.
  2. Brizancja opisuje lokalny efekt kruszący wynikający z gwałtowności frontu ciśnienia. Mierzona pośrednio przez testy Trauzla lub BATES, koreluje z P_CJ, ale nie z impulsem fali podmuchowej.
  3. Zdolność podmuchowa zależy od impulsu fali podmuchowej i czasu jej trwania — nie od samego nadciśnienia szczytowego. Materiały o niskiej VoD i brizancji (ANFO, FAE) mogą mieć wyższą podmuchowość na dalekim polu niż materiały kruszące o wysokim VoD.
  4. Skalowanie Hopkinsona-Cranza pozwala przeliczać wyniki testów małoskalowych na realne scenariusze — podstawowe narzędzie w inżynierii bezpieczeństwa wybuchowego.
  5. Równoważnik TNT jest przybliżeniem, którego dokładność zależy od odległości, geometrii i mechanizmu odpowiedzi — nie należy go traktować jako precyzyjnej miary zdolności materiałowej.
  6. Soczewki wybuchowe w ładunkach implozyjnych są doskonałym przykładem, w którym absolutna wartość VoD jest mniej ważna niż różnica VoD między materiałami i jednorodność produkcyjna.1,2,3
  7. Wykresy P-I (Pressure-Impulse) integrują wrażliwość ciśnieniową i impulsową konstrukcji w jednej płaszczyźnie — standardowe narzędzie inżynierii ochrony budynków przed wybuchami, uzupełniające skalowanie Hopkinsona-Cranza.

Dodatkowe materiały multimedialne

Przy kolejnej redakcji warto dodać tabelę porównawczą: VoD / brizancja / impuls fali / typowe zastosowanie.

Powiązane kalkulatory i narzędzia

  • Podmuch — wyznacza nadciśnienie, impuls i falę odbitą z modeli Kingery-Bulmash/CONWEP.

Ćwiczenia praktyczne

Pierwsze ćwiczenie powinno polegać na rozdzieleniu trzech parametrów na przykładach. Należy:

  1. wybrać trzy hipotetyczne materiały wybuchowe,
  2. przypisać im różne wartości VoD, brizancji i impulsu podmuchu,
  3. wskazać, który lepiej nadaje się do kruszenia, który do podmuchu, a który do sterowania falą,
  4. odnieść wynik do soczewki wybuchowej,
  5. sformułować wniosek, dlaczego jedna „miara mocy” nie wystarcza.

Drugie ćwiczenie powinno dotyczyć błędnych intuicji. Należy:

  1. założyć materiał o bardzo wysokiej prędkości detonacji,
  2. sprawdzić, czy automatycznie daje on najlepszy efekt podmuchowy,
  3. porównać jego działanie z ładunkiem o dłuższej, mniej gwałtownej fali,
  4. powiązać to z FAE i amunicją termobaryczną,
  5. wyciągnąć wniosek, co naprawdę powinien mierzyć projektant zależnie od celu bojowego.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły