Tensometry, szybkie kamery i czujniki ciśnienia w badaniach fali uderzeniowej

Streszczenie

Badanie fali uderzeniowej nie polega na „zrobieniu wybuchu i popatrzeniu”. Jest to problem synchronizacji czujników ciśnienia, tensometrów, rejestratorów, kamer i modeli numerycznych, które widzą różne projekcje tego samego zdarzenia. Czujnik ciśnienia mierzy historię obciążenia w punkcie, tensometr pokazuje odpowiedź konstrukcji, a szybka kamera lub obrazowanie schlieren ujawnia geometrię frontu fali.^1,2

Taki artykuł jest potrzebny jako bezpieczne zaplecze metrologiczne dla tekstów o fali podmuchowej, Mach stem, modelach pressure-impulse i odporności konstrukcji. Celem nie jest instrukcja przeprowadzania prób wybuchowych, lecz pokazanie, jak myśli się o aparaturze, kalibracji, czasie, niepewności i zgodności pomiaru z modelem.

Rozszerzenie tematu

Jedno zdarzenie, wiele kanałów

Fala uderzeniowa jest krótkim, gwałtownym zdarzeniem. W jednym miejscu widzimy nagły skok ciśnienia. W innym miejscu konstrukcja zaczyna się uginać. Kamera rejestruje ruch frontu, pyłu albo znacznika. Symulacja numeryczna zwraca pola ciśnienia, przemieszczenia i naprężenia. Żaden z tych obrazów sam nie wystarcza.

Dlatego stanowisko badawcze składa się z wielu kanałów:

• czujników ciśnienia lub mikrofonów impulsowych,
• tensometrów na elementach konstrukcji,
• akcelerometrów i czujników przemieszczenia,
• szybkich kamer, kamer smużeniowych albo obrazowania schlieren,
• rejestratorów o znanym paśmie i czasie próbkowania,
• wspólnego triggera i zegara,
• procedury kalibracji oraz opisu niepewności.

Warto od razu rozdzielić dwie rzeczy. Pomiar fali jako takiej i pomiar odpowiedzi konstrukcji to nie to samo. Czujnik ciśnienia mówi, co „przyszło” do przegrody. Tensometr mówi, co przegroda z tym zrobiła. Szybka kamera pozwala sprawdzić, czy geometria fali lub ruch konstrukcji zgadza się z interpretacją.

Czujnik ciśnienia: historia obciążenia

Najbardziej bezpośrednim pomiarem jest przebieg ciśnienia w czasie. Dla fali podmuchowej interesuje nas nie tylko pik nadciśnienia, ale także dodatnia faza, czas przyjścia, impuls właściwy i ewentualna faza ujemna. Artykuł Xiao, Andrae i Gebbeken pokazuje, że w przypadku wybuchu na wysokości przebieg ciśnienia zależy od tego, czy punkt obserwacji znajduje się poniżej, czy powyżej punktu potrójnego: poniżej może pojawić się pojedynczy pik Mach stem, powyżej dwa kolejne piki fali padającej i odbitej.¹

To jest ważne dydaktycznie. Ten sam eksperyment nie ma jednego „ciśnienia wybuchu”. Ma pole ciśnienia zmienne w przestrzeni i czasie. Czujnik mierzy tylko lokalny wycinek tego pola. Jeśli nie znamy położenia czujnika względem frontu fali, powierzchni odbijających i obiektu, interpretacja piku może być myląca.

Czujnik ciśnienia musi mieć odpowiednie pasmo, czas narastania, zakres amplitudy i mocowanie. Zbyt wolny czujnik wygładzi pik. Zbyt długi przewód może dodać zakłócenia. Źle ustawiony czujnik może mierzyć mieszaninę fali padającej, odbitej i lokalnych zaburzeń. Dlatego w publikacji naukowej sam wykres ciśnienia bez geometrii stanowiska jest niewystarczający.

Tensometr: odpowiedź konstrukcji

Tensometr nie mierzy fali w powietrzu, lecz odkształcenie elementu. Zwykle jest częścią mostka tensometrycznego: niewielka zmiana rezystancji przekłada się na napięcie, które trzeba wzmocnić, odfiltrować, przeliczyć i zarejestrować. W statycznych pomiarach siły i masy taki tor może być stosunkowo wolny. W badaniach dynamicznych trzeba zwrócić uwagę na pasmo całego toru, nie tylko samego tensometru.

Karta miernika CL 363 pokazuje typowy przyrząd do czujników tensometrycznych zasilanych napięciowo. Ma trzy kanały, zasilanie mostków 5 V lub 10 V, zakresy wejściowe w miliwoltach, komunikację szeregową i funkcje minimum/maksimum, tarowania i uśredniania.³ To dobry przykład aparatury przemysłowej i laboratoryjnej do tensometrii, ale jednocześnie przypomnienie: miernik o czasie pomiaru rzędu dziesiątych części sekundy nie jest rejestratorem nanosekundowego frontu fali. Nadaje się do wolniejszych obciążeń, kalibracji, pomiarów quasi-statycznych albo monitoringu, a nie do uchwycenia samego piku szoku.

W artykule edukacyjnym warto więc pytać: czy tensometr ma mierzyć szybki impuls, czy końcowe ugięcie? Czy rejestrator ma wystarczające pasmo? Czy mostek jest skompensowany temperaturowo? Czy przewody są odporne na zakłócenia? Czy odkształcenie lokalne reprezentuje globalną odpowiedź konstrukcji?

Szybka kamera i schlieren: geometria frontu

Fala uderzeniowa w powietrzu jest niewidoczna gołym okiem, ale zmienia gęstość powietrza, a więc współczynnik załamania światła. Technika schlieren i techniki cieniowe wykorzystują gradient gęstości, aby uwidocznić front fali. Xiao i współautorzy generowali numeryczne obrazy schlieren z gradientów gęstości uzyskanych w symulacjach, a następnie wyznaczali ścieżkę punktu potrójnego.¹

To piękny przykład połączenia pomiaru i modelu. Obraz schlieren mówi, gdzie jest front fali i jak rozwija się Mach stem. Czujnik ciśnienia mówi, co dzieje się w konkretnym punkcie. Model numeryczny pozwala porównać oba obrazy. Dopiero razem tworzą spójną interpretację.

Szybka kamera może też mierzyć ruch elementu konstrukcji: ugięcie płyty, przemieszczenie markera, prędkość odłamka w bezpiecznym eksperymencie materiałowym albo propagację fali w przezroczystym medium. W dydaktyce można używać bezpiecznych analogów: impulsu akustycznego, balonu, sprężonego powietrza o bardzo małej energii, fali na wodzie albo danych syntetycznych. Nie trzeba wykonywać prób z materiałami wybuchowymi, aby nauczyć synchronizacji kamery i czujnika.

Historia aparatury do badań fali uderzeniowej

Dynamiczna aparatura pomiarowa dla badań fal uderzeniowych rozwijała się równolegle z bronią i ochroną przed nią, a jej historia jest ściśle spleciona z historią XX-wiecznych konfliktów zbrojnych.

Pierwsza wojna światowa: Artyleria o doniosłości przekraczającej 100 km generowała fale akustyczne i uderzeniowe, które były rejestrowane przez sieć mikrofonów i czujników sejsmicznych — pierwotnie po to, by triangulować pozycje baterii przeciwnika. To była prawdopodobnie pierwsza masowa aplikacja pomiarów fali ciśnienia w czasie rzeczywistym w celach wojskowych. Technika „sound ranging" była stosowana przez wszystkie strony konfliktu.

Druga wojna światowa i Project Manhattan: Projekt Manhattan wymagał kwantyfikacji efektów fali podmuchowej dla oceny skuteczności broni jądrowej i konwencjonalnej. Pierwsze testy empiryczne przeprowadzono na poligonach z manometrami impulsowymi i szybkimi kamerami. Jednym z przełomów było użycie kamery smużeniowej (streak camera), gdzie obracające się zwierciadło zamieniało obraz na film z rozdzielczością mikrosekund. Luis Alvarez z Berkeley budował układy elektroniczne do synchronizacji detonatorów i rejestracji czasu implozji — to były pierwsze wielokanałowe systemy synchronizowane co do nanosekund.

Lata 50-60: program badań nad falą i schronami. USA i ZSRR prowadziły intensywne badania nad skutkami broni jądrowej na struktury i ludzki organizm. Programy Kirtland Air Force Base, Edgewood Arsenal i Natick Labs produkowały dane o odporności konstrukcji, granicach biologicznych tolerancji i skuteczności schronów. Aparatura obejmowała: czujniki piezoelektryczne (PCB Piezotronics), manometry impulsowe, tensometry przewodowe na stalowych elementach, szybkie kamery Beckman & Whitley i mikrofony kondensatorowe.

Lata 70-80: cyfryzacja. Pojawienie się szybkich ADC i oscyloskopów DSO umożliwiło rejestrację danych cyfrowych zamiast analogowych. Programy takie jak DYNA2D i DYNA3D (Lawrence Livermore) pozwoliły porównywać eksperymentalne przebiegi ciśnienia z symulacjami numerycznymi. Tensometry stały się standardowym narzędziem kwalifikacji schronów wojskowych.

XXI wiek: Cyfrowa korelacja obrazu (DIC), ultra-szybkie kamery CMOS (Photron, Phantom Vision), laserowe mierniki przemieszczenia (LDV) i systemy akwizycji danych z FPGA zastąpiły i uzupełniły tradycyjne metody. Programy obliczeniowe LS-DYNA, AUTODYN i Abaqus pozwalają na symulacje wielofizyczne z jawną detonacją, propagacją fali, kontaktem i zniszczeniem. Walidacja modeli przez dane eksperymentalne jest standardową procedurą publikacyjną.

Typy czujników ciśnienia dla fal dynamicznych

Dla impulsów wysokiego ciśnienia o krótkim czasie trwania (mikro–milisekundy) stosuje się kilka typów czujników:

Czujniki piezoelektryczne (PVDF i ceramiczne PZT): Generują ładunek elektryczny proporcjonalny do naprężenia. Bardzo szybkie (czas narastania < 1 μs), odporne na uszkodzenia mechaniczne, ale dają sygnał proporcjonalny do pochodnej naprężenia (co wymaga całkowania elektronicznego lub korekcji). Producenci: PCB Piezotronics, Kistler, Endevco. Zakresy: 0,01–7000 bar.

Czujniki piezoresistywne (MEMS): Zmiana rezystancji pod wpływem naprężenia. Dają sygnał DC (absolutna wartość ciśnienia), szybkość poniżej 50 μs, dobre liniowość. Producenci: Kulite, Meggitt, Honeywell. Stosowane np. w czujnikach ciśnienia powietrza w samochodach, ale też w wersjach szybkich dla badań udarowych.

Czujniki przetwornika ładunku (charge mode vs. ICP): Piezoelektryczne czujniki mogą pracować w trybie ładunkowym (wymagają zewnętrznego wzmacniacza ładunku) lub ICP (Integrated Circuit Piezoelectric), gdzie wzmacniacz jest wbudowany w czujnik. ICP jest wygodniejszy (standardowy kabel BNC), ale nie działa w środowiskach wysokich temperatur.

Manometry zapieczętowane (sealed): Do ciśnień stałych lub quasi-statycznych, nie dla impulsów nanosekund. Użyteczne do pomiaru szczytowego nadciśnienia po fali, gdy pik minął.

Mikrofony impulsowe: Czujniki kondensatorowe kalibrowane w dB re 20 μPa. Stosowane w pomiarach akustycznych pola dalekiego wybuchu, gdzie ciśnienia są niskie (do 200 Pa). W polu bliskim (nadciśnienia > 0,01 MPa) używa się wyżej wymienionych czujników.

Typ czujnika	Zakres ciśnienia	Czas narastania	Zakres częstotliwości	Zastosowanie
Piezoelektryczny (PCB)	0,1–1000 bar	< 1 μs	DC–1 MHz	udary bliskie, detonacja
Piezoresistywny (Kulite)	0,01–700 bar	< 50 μs	DC–100 kHz	lekkie udary, silniki
Mikrofon wolnociśnieniowy	0,1 Pa–200 Pa	1 μs	0,1 Hz–100 kHz	pole dalekie, akustyka
Manometr quasi-statyczny	0–1 bar	wolny	DC	monitoring, ciśnienie tła

Tensometry — fizyka i praktyka

Tensometr rezystancyjny (strain gauge) to cienki drut lub folia układająca się w sposób zwiększający długość ścieżki aktywnej, umieszczone na nośniku poliestrowym lub poliimidowym. Jego działanie opiera się na efekcie Piezoresistivnym: przy odkształceniu ε zmienia się rezystancja:

ΔR/R = GF × ε

gdzie GF (Gauge Factor, współczynnik czułości) typowo wynosi ~2,0 dla folii konstantanowych i ~150–200 dla tensometrów półprzewodnikowych.

Mostki Wheatstone'a: tensometry są łączone w mostki rezystancyjne, co pozwala:

• eliminować temperaturę (tensometr aktywny + tensometr referencyjny w sąsiednim ramieniu)
• zwiększać czułość (full-bridge z czterema aktywnymi tensometrami)
• mierzyć kierunek odkształcenia (rozety trójwymiarowe)

Dla badań dynamicznych ważne są:

• Czas odpowiedzi: tensometr foliowy < 1 μs, ale konduktancja przewodów dodaje RC ≈ 1–10 μs
• Pasmo kondycjonera: miernik tensometryczny CL 363 ma pasmo do ~1 kHz — dla badań fali uderzeniowej o odzwierciedleniu w ms konieczny jest szybszy wzmacniacz (np. do 100 kHz)
• Samonagrzewanie: prąd zasilający powoduje nagrzewanie mostka; przy impulsach krótkich to nie problem, ale przy ciągłym monitoringu może zaburzać pomiar temperatury

Szybkie kamery — technologie i parametry

Nowoczesne szybkie kamery (high-speed cameras) rozróżniamy na kilka kategorii:

Kamery CMOS z globalnym zcięciem (global shutter): Sensor CMOS pobiera cały kadr jednocześnie, bez toczenia linii po linii. Kluczowe dla badań fali, gdzie obiekt szybko się porusza. Phantim Miro (Vision Research), Photron Fastcam, IDT MotionPro. Parametry:

• Fastcam SA-Z: do 2 100 000 fps przy rozdzielczości 256×16
• Phantom v2512: do 1 000 000 fps przy 128×8 px, 25 600 fps przy 1280×800

Kamery smużeniowe (streak cameras): Zamknięta przekształceniak optyczny rozciąga obraz szczeliny na oś czasu. Rozdzielczość czasowa: 0,2–50 ps dla optycznych streak cameras Hamamatsu. Nie dają obrazu 2D w czasie, ale dają czas z rozdzielczością pikosekund wzdłuż jednej osi. Stosowane do diagnostyki laserów, implozji i emisji plazmowych.

Kamery z wzmocniaczem obrazu (ICCD): Intensifier (tube MCP) przed matrycą CCD pozwala na ekspozycje nanosekund i bramkowanie elektryczne (electronic gating). Stosowane w badaniach plazmowych, fluorescencji i emisji lekkich z implozji.

Dla badań fali uderzeniowej typowe zastosowanie: Fastcam SA-Z przy 100 000–500 000 fps i rozdzielczości 512×256 pozwala filmować propagację frontu z rozdzielczością przestrzenną 1 mm i czasową 2–10 μs. Synchronizacja: wejście TTL trigger od układu inicjacyjnego.

Cyfrowa korelacja obrazu (DIC) jako nowoczesna metoda

DIC (Digital Image Correlation) to technika pomiarowa, która wyznacza przemieszczenia i odkształcenia powierzchni poprzez porównanie obrazów przed i po deformacji. Jest szczególnie cenna w badaniach dynamicznych, bo:

• Nie wymaga przyklejania tensometrów — pomiar jest bezkontaktowy
• Daje pole odkształceń (2D lub 3D z dwóch kamer), a nie pomiar punktowy
• Można mierzyć do milionów punktów jednocześnie
• Współpracuje ze szybkimi kamerami CMOS

Zasada działania: wzorzec losowy (random speckle pattern) jest nakładany na powierzchnię (spray lub naturalna faktura). Kolejne klatki kamery są porównywane przez korelację bloków pikseli. Przemieszczenie bloku pikseli jest wyznaczane z sub-pikselową dokładnością (typowo 0,01–0,05 px = 0,1–0,5 μm przy powiększeniu 10 μm/px).

Dla badań fali uderzeniowej DIC pozwala na:

• Pomiar ugięcia płyty lub belki w każdym punkcie
• Wyznaczanie pola odkształceń naprężeń głównych
• Identyfikację miejsc koncentracji naprężeń
• Porównanie z wynikami LS-DYNA w każdym punkcie

Polska perspektywa — aparatura dynamiczna

Wojskowy Instytut Techniczny Uzbrojenia (WITU) w Zielonce jest jedynym w Polsce certyfikowanym laboratorium do badań materiałów wybuchowych i detonatorów z kompleksową infrastrukturą do pomiarów dynamicznych. Wyposażenie obejmuje szybkie kamery, czujniki ciśnienia PCB/Kistler, tensometry na stanowiskach strzałowych i systemy synchronizacji.

Wojskowa Akademia Techniczna (WAT) w Warszawie — Wydział Mechatroniki — prowadzi badania nad ochroną przed wybuchem, w tym tensometryczne pomiary ugięcia elementów pancernych i numeryczne modelowanie LS-DYNA. WAT jest jedną z głównych uczelni badawczych w obszarze mechaniki wybuchu i ochrony balistycznej.

Politechnika Krakowska — Zakład Budownictwa Obronnego — prowadzi badania nad ochroną budowli przed skutkami wybuchu, z uwzględnieniem walidacji eksperymentalnej. Tensometria i szybkie kamery są używane do pomiarów ugięcia i uszkodzeń betonowych elementów.

Instytut Chemii i Techniki Jądrowej (IChTJ) w Warszawie prowadzi badania w obszarach radiochemii i dozymeetrii, ale ma też kompetencje w zakresie instrumentacji pomiarowej dla systemów jądrowych. Techniki pomiaru dynamicznego pojawiają się przy badaniu inicjatorów i elementów bezpieczeństwa systemów jądrowych.

Numeryczne przykłady dla dydaktyki

Przykład 1: Wymagania pasma czujnika ciśnienia

Typowy impuls Friedlandera dla 1 kg TNT w odległości 3 m:

• Czas narastania frontu: t_r ≈ 50–200 μs
• Czas trwania fazy dodatniej: t_d ≈ 2–5 ms
• Maksymalne nadciśnienie: P_0 ≈ 30–80 kPa

Wymagane pasmo czujnika:
BW = 0,35/t_r = 0,35/(50×10^-6) ≈ 7 kHz

Dla tego impulsu czujnik o paśmie 10 kHz jest wystarczający. Miernik CL363 z pasmem ~1 kHz byłby zbyt wolny — zaokrągliłby pik. Czujnik PCB o paśmie 100 kHz daje duży zapas.

Przy mniejszych ładunkach i mniejszych odległościach czas narastania może spaść do 1–10 μs, wymagając pasma 35–350 kHz. Stąd profejsonalne czujniki piezoelektryczne mają pasmo > 100 kHz.

Przykład 2: Rozdzielczość kamery dla frontu fali

Prędkość propagacji fali uderzeniowej dla nadciśnienia 100 kPa w powietrzu:
M_s ≈ 1,29 (liczba Macha), więc v_s = M_s × c = 1,29 × 343 m/s ≈ 442 m/s

Aby rozdzielić dwa położenia frontu odległe o 1 mm:
Δt = 1 mm / 442 m/s ≈ 2,3 μs

Wymagana częstotliwość klatki kamery:
fps = 1/Δt = 1/(2,3×10^-6) ≈ 430 000 fps

Fantom v2512 przy tej fps ma rozdzielczość 256×256 px — wystarczającą do śledzenia frontu z 1 mm rozdzielczością przy polu 256 mm × 256 mm.

Przykład 3: Tensometr — przeliczenie sygnału

Tensometr foliowy K-AA (Vishay), GF = 2,07, mostek pełny (cztery tensometry):

• Wzmocnienie instrumentacji: A = 100
• Zasilanie mostka: V_ex = 5 V
• Zmierzone napięcie mostka: V_out = 8,5 mV

Odkształcenie obliczone:
ε = (V_out × 4) / (V_ex × A × GF) = (0,0085 × 4) / (5 × 100 × 2,07) = 0,034/1035 ≈ 32,9 μstrain = 32,9 × 10^-6

Naprężenie przy module Younga stali E = 210 GPa:
σ = E × ε = 210×10^9 × 32,9×10^-6 ≈ 6,9 MPa

Przeliczenie jest proste, ale wymaga dobrze zmierzonego V_out bez szumu, kompensacji temperatury i wiedzy o GF z karty tensometru.

Czas jest osobnym kanałem pomiarowym

W pomiarach dynamicznych czas nie jest tłem. Jest mierzoną wielkością. Trzeba znać:

• czas przyjścia fali do czujnika,
• opóźnienie triggera,
• opóźnienie elektroniki,
• różnice długości kabli,
• czas ekspozycji kamery,
• jitter zegara,
• częstotliwość próbkowania,
• filtrację i jej opóźnienie fazowe.

Jeśli czujnik ciśnienia, kamera i tensometr nie mają wspólnego odniesienia czasu, ich porównanie może być pozorne. W artykule o FPGA jako generatorze sekwencji i akwizytorze danych opisano architekturę, w której jeden układ rozdziela trigger, znaczniki czasu i metadane jakości pomiaru. To właśnie taki sposób myślenia jest potrzebny w badaniach fali.

W starych laboratoriach rolę tę pełniły oscyloskopy, kamery smużeniowe, generatory opóźnień i układy fotograficzne. Współcześnie dochodzą FPGA, szybkie ADC, kamery z wejściem triggera i rejestratory wielokanałowe. Problem pozostaje ten sam: trzeba wiedzieć, czy dwa zdarzenia naprawdę były jednoczesne, czy tylko tak wyglądają na osobnych wykresach.

Systemy rejestracji wielokanałowej i synchronizacja

Badania fali uderzeniowej generują dane z wielu źródeł jednocześnie. System akwizycji danych (DAQ) musi obsłużyć synchronizację i zapis wszystkich kanałów z wspólnym odniesieniem czasu:

Architektura typowego systemu DAQ dla badań fali:

1. Generator triggera: układ elektroniczny (timer z pyrotechnicznym inicjatorem lub generator elektryczny) wyzwala zdarzenie i jednocześnie wysyła sygnał TTL do systemu DAQ.

2. Pretrigger buffer: rejestratory muszą zbierać dane z marginesem przed zdarzeniem. Typowo 10–50 ms przed wyzwoleniem. Pozwala to sprawdzić, czy tło (poziom szumu) przed zdarzeniem był normalny.

3. Wielokanałowy rejestrator o stałym zegarze: wszystkie kanały (czujniki ciśnienia, tensometry, akcelerometry) pobierają próbki z tego samego zegara. Wymóg: skew między kanałami < 1 próbka.

4. Dedykowany kanał triggera: jeden kanał rejestratora jest podłączony do sygnału triggera. Pozwala powiązać timeline danych z chwilą inicjacji.

5. Synchronizacja kamery: szybka kamera ma własny zegar dla klatek, ale jest wyzwalana sygnałem TTL z tego samego systemu. Kamera rejestruje znacznik czasu wyzwolenia (trigger latency), który jest podawany przez producenta lub kalibrowany.

6. Metadane geometryczne: każdy czujnik ma zdefiniowane współrzędne w systemie laboratoryjnym. Bez tego nie można porównać przebiegu ciśnienia z pozycją frontu z kamery.

Dla ćwiczeń dydaktycznych taki system można zbudować w miniaturze: rejestratory USB (np. National Instruments myDAQ, Digilent Analog Discovery) z wejściem triggera, oscyloskop jako drugi rejestrator i smartfon jako szybka kamera (240–960 fps) — dla bezpiecznych eksperymentów akustycznych i mechanicznych.

Skew między kanałami: Dwa kanały ADC w tym samym rejestratorze są zwykle synchroniczne z dokładnością < 1 ns. Ale zewnętrzna kamera z własnym zegarem może mieć skew rzędu mikrosekund lub dziesiątek milisekund, jeśli trigger latency nie jest skalibrowana. Producenci kamer podają trigger latency (np. Phantom Vision: < 1 μs dla external trigger) — wartość tę należy uwzględnić przy porównywaniu timeline klatka-czujnik.

Formatowanie danych: dane z różnych rejestratorów często mają różne formaty (CSV, HDF5, binarny, Matlab). W spójnym systemie DAQ wszystkie dane trafiają do wspólnego pliku z metadanymi (timestamp, konfiguracja, geometria, kalibracja). Oprogramowanie takie jak DIAdem (NI), LabVIEW, lub Python z h5py pozwala unifikować i analizować.

Schlieren i optyczne metody wizualizacji

Fala uderzeniowa jest niewidoczna optycznie, ale zmienia gęstość powietrza, a więc zmienia gradient współczynnika załamania:

n(ρ) = 1 + K_GD × ρ

gdzie K_GD to stała Gladstone-Dale'a (~2,3×10^-4 m^3/kg dla powietrza, w przeliczeniu na n). Zmiana gęstości δρ zmienia n o:

δn ≈ K_GD × δρ

Dla fali uderzeniowej o nadciśnieniu 100 kPa (przyrost gęstości ~40%): δn ≈ 2,3×10^-4 × 0,4 × 1,2 kg/m^3 ≈ 1,1×10^-4

To jest bardzo mała zmiana (11 ppm), ale wystarczająca do detekcji przez układ schlieren lub interferometr.

Układ schlieren z jednym zwierciadłem (Z-schlieren): równoległa wiązka światła przechodzi przez obszar z gęstością gradientu, jest ogniskowana i przecięta nożem schlieren w ognisku. Zmiana kąta ugięcia (proporcjonalna do ∂n/∂y) zmienia kontrast obrazu. Układ wymaga: źródło punktowe (LED lub laser), zwierciadło paraboliczne, nóż schlieren i kamera.

Metoda BSTS (Background Subtracted Schlieren): nowoczesna odmiana, gdzie wiele klatek jest rejestrowanych i odejmowanych od tła. Wynik jest bardzo czysty wizualnie i może być kwantytywny.

Interferometria Macha-Zehndera i Michelson: mierzą absolutną zmianę fazy optycznej, co daje kwantytywne informacje o gęstości. Używana do diagnostyki laserowej, plazmowej i gazowej. Wymagają laserowego źródła i precyzyjnej optyki.

Dla celów dydaktycznych schlieren jest najdostępniejszą techniką. Układy DIY z zwierciadłem samochodowego reflektora, LED i aparatem fotograficznym pozwalają zobaczyć termiczne prądy konwekcyjne, sygnał słowny i balonowe impulsy akustyczne. To ta sama fizyka, którą stosują laboratoria badawcze, tylko ze skromniejszą rozdzielczością.

Integracja z modelami numerycznymi

Jednym z najważniejszych celów eksperymentu dynamicznego jest walidacja modelu numerycznego. Schemat walidacji:

1. Zdefiniuj geometrię eksperymentu: wymiary, materiały, odległości, masy.

2. Zbuduj model numeryczny w LS-DYNA, AUTODYN lub Abaqus z tą samą geometrią. Użyj dostępnych modeli materiałowych (Johnson-Cook dla metali, RHT dla betonu, JWL dla materiału wybuchowego).

3. Przeprowadź eksperyment i zarejestruj przebiegi ciśnienia, ugięcia, prędkości w kilku punktach.

4. Porównaj przebiegi: czas przyjścia fali, pik ciśnienia, czas trwania, impuls, maksymalne ugięcie, pozostałe odkształcenie.

5. Oceń rozbieżności: czy model jest zbyt sztywny, zbyt miękki? Czy stała materiałowa musi być skorygowana? Czy geometria siatki jest wystarczająca?

6. Kalibracja modelu: dostosuj jeden parametr (np. dynamiczna granica plastyczności σ_y0) i sprawdź, czy porównanie się poprawia. Nie kalibruj wielu parametrów jednocześnie bez fizycznej motywacji — to prowadzi do overfittingu.

Artykuł Kosiuczenko stosuje tę metodologię dla elementów ochronnych budowli pod wpływem fali wybuchu jądrowego. Model zawiera beton i stal, funkcję Brode'a dla ciśnienia i model Johnson-Cook z krzywą zniszczenia dla stali. Wyniki pokazują mapy ciśnienia, przemieszczenia i naprężeń von Misesa — każda z nich może być teoretycznie porównana z danymi z tensometru i kamery.²

Budżet niepewności pomiaru dynamicznego

Pomiar fali uderzeniowej ma wieloskładnikowy budżet niepewności:

Źródło niepewności	Typowy wkład
Kalibracja czujnika ciśnienia	1–3%
Wzmacniacz i tor analogowy	0,5–2%
Kwantyzacja ADC (12–16-bit)	< 0,01% (zanikowe)
Opóźnienie triggera (jitter triggera)	0,1–1 μs → 0,1–10% przy krótkich impulsach
Skew między kanałami	1–10 ns → mały dla ms-impulsów
Niepewność geometrii (pozycja czujnika)	1–5 mm → 0,5–2% dla 1/r^2
Powtarzalność źródła (ładunku)	2–10% (dla materiałów pyrotechnicznych)
Model Friedlandera (przybliżenie)	5–15% dla dalekiego pola
Suma złożona	~5–15%

Największa niepewność pochodzi zwykle z niepowtarzalności źródła i uproszczenia modelu, a nie z elektroniki pomiarowej. Dlatego przy dobrej aparaturze powtarzalność wynika z fizyki źródła, a nie z jakości ADC.

Pasmo, próbkowanie i aliasing

Szybki impuls wymaga odpowiedniego pasma. Jeżeli czujnik, wzmacniacz albo rejestrator są zbyt wolne, pik zostanie zaniżony i rozciągnięty. Jeżeli próbkowanie jest zbyt rzadkie, czas przyjścia fali i maksymalna wartość będą obarczone dużym błędem. Jeżeli filtracja jest nieopisana, nie wiadomo, czy wygładzenie pochodzi ze zjawiska, czy z toru pomiarowego.

To samo dotyczy kamery. Liczba klatek na sekundę nie wystarcza. Trzeba znać czas ekspozycji, sposób synchronizacji, opóźnienie wyzwolenia i rozdzielczość przestrzenną. Kamera może mieć bardzo dużą liczbę klatek na sekundę, ale jeśli ekspozycja jest za długa, szybko poruszający się front zostanie rozmyty.

W pomiarach studenckich można to pokazać bez żadnego wybuchu: nagrać szybko poruszający się marker, zarejestrować impuls z generatora na oscyloskopie i porównać wynik przy różnym paśmie oraz próbkowaniu. To ta sama lekcja, która później wraca przy czujnikach ciśnienia fali uderzeniowej.

Kalibracja czujników dynamicznych

Kalibracja czujnika dynamicznego różni się od statycznej. Czujnik może być doskonale skalibrowany statycznie (odpowiadać liniowo na ciśnienie stałe), a mimo to zniekształcać impulsy dynamiczne przez własne pasmo i rezonansy.

Kalibracja quasi-statyczna: porównanie z manometrem referencyjnym przy powolnej zmianie ciśnienia. Daje współczynnik przeliczeniowy V/bar, ale nic nie mówi o zachowaniu przy μs-impulsach.

Kalibracja dynamiczna metodą Hopkinsona (SHPB): Pręt Hopkinsona (Split Hopkinson Pressure Bar) generuje kontrolowany impuls o czasie narastania ~10–50 μs przy prędkości odkształcenia 10^3–10^4/s. Używany zarówno do kalibracji czujników jak i do badań właściwości materiałów. Impulsy z SHPB mają dobrze zdefiniowany kształt — porównując sygnał referencyjny z sygnałem z czujnika, można wyznaczić odpowiedź częstotliwościową czujnika.

Kalibracja w rurze udarowej (shock tube): Rura z membranną rozdzielającą dwie komory (wysokie i niskie ciśnienie). Pęknięcie membrany generuje falę uderzeniową o kontrolowanym ciśnieniu i czasie narastania poniżej 1 μs. Porównując wynik czujnika z modelem analitycznym rury, kalibruje się dynamicznie.

W Polsce WITU posiada rurę udarową do kalibracji czujników ciśnienia i tensometrów dynamicznych. Jest to rzadki zasób na skalę krajową.

Podsumowanie dydaktyczne

Badania fali uderzeniowej to jeden z nielicznych obszarów inżynierii, gdzie elektronika pomiarowa, fizyka zjawisk szybkich, mechanika strukturalna i modelowanie numeryczne spotykają się w jednym eksperymencie. Dla studenta fizyki jądrowej lub inżynierii wojskowej jest to idealny model złożonego problemu metrologicznego:

• Czujnik mierzy tylko lokalny wycinek zdarzenia — jego interpretacja wymaga geometrii i modelu
• Tor pomiarowy (kabel, wzmacniacz, ADC, software) jest część układu, a nie neutralne medium
• Synchronizacja wielu kanałów jest warunkiem koniecznym sensownego porównania danych
• Walidacja modelu numerycznego przez eksperyment jest procesem iteracyjnym, nie jednorazowym testem
• Budżet niepewności łączy czujnik, geometrię, źródło i model

Wszystkie te lekcje przenoszą się na inne dziedziny: spektrometrię gamma (jeden detektor, jeden aspekt zdarzenia), koincydencje (synchronizacja wielu detektorów), modelowanie reaktorów (walidacja kodu przez dane eksploatacyjne). Metrologia dynamiczna jest zatem szkołą myślenia o pomiarze w szerszym sensie.

Dla celów tego serwisu kluczowe jest, że ten artykuł pozostaje po bezpiecznej stronie granicy: omawia aparaturę, kalibrację, synchronizację i interpretację danych — bez opisu jakichkolwiek prób z materiałami niebezpiecznymi, bez specyfikacji parametrów stanowisk strzałowych, bez ilości ładunków czy geometrii rozmieszczenia. Fizyka fal uderzeniowych i metrologia szybkich zjawisk mogą być nauczane i badane w całości bez wchodzenia w praktykę certyfikowanych laboratoriów strzałowych. Takie podejście — rozumienie fizyki bez wykonywania niebezpiecznych eksperymentów — jest modelowym przykładem odpowiedzialnej edukacji w obszarze wrażliwych technologii. Student może dobrze rozumieć problem synchronizacji czujników, budżet niepewności i relację między eksperymentem a modelem, opierając się na danych literaturowych, syntetycznych przykładach i bezpiecznych modelach dydaktycznych.

Kalibracja nie dotyczy tylko czujnika. Dotyczy całego toru:

zjawisko -> czujnik -> przewód -> kondycjoner -> filtr -> ADC -> oprogramowanie -> wykres

Każdy element ma wzmocnienie, pasmo, szum, opóźnienie i nieliniowość. W tensometrii dochodzi zasilanie mostka, kompensacja temperatury, tarowanie i przeliczanie napięcia na siłę, ciśnienie lub odkształcenie. Karta CL 363 pokazuje, że przemysłowy miernik zawiera funkcje przeliczania wyników, pamięć parametrów i uśrednianie.³ To wygodne, ale w badaniach dynamicznych trzeba wiedzieć, czy te funkcje nie ukrywają szczegółów czasowych.

Dobra dokumentacja pomiaru powinna zawierać:

• typ czujnika i jego zakres,
• pasmo lub czas odpowiedzi,
• kalibrację statyczną i, jeśli potrzebna, dynamiczną,
• długości i typy przewodów,
• ustawienia filtrów,
• częstotliwość próbkowania,
• sposób synchronizacji,
• geometrię stanowiska,
• niepewność przeliczenia i czasu.

Bez tego wykres może być ładny, ale nie musi być dowodem.

Od pomiaru do modelu

Kosiuczenko pokazuje przykład, w którym obciążenie fali podmuchowej jest analizowane numerycznie w LS-DYNA z użyciem funkcji Brode'a i modelu materiałowego Johnson-Cook dla stali.² W pracy pojawiają się mapy ciśnienia, przemieszczeń, naprężeń von Misesa i energie układu. To dobry przykład, że aparatura eksperymentalna i symulacja nie są konkurentami.

Pomiar daje punkty rzeczywistości. Model daje pełne pole, którego nie da się zmierzyć wszędzie. Porównanie obu pozwala sprawdzić, czy model nie jest tylko kolorową animacją. W badaniu fali uderzeniowej szczególnie ważne są:

• przebieg ciśnienia w czasie,
• maksymalne przemieszczenie elementu,
• czas osiągnięcia maksimum,
• rozkład odkształceń,
• energia kinetyczna i wewnętrzna,
• położenie frontów i odbić.

Jeśli czujnik ciśnienia wskazuje inny czas przyjścia niż model, problem może leżeć w geometrii, triggerze, definicji punktu pomiarowego albo samym modelu fali. Jeśli tensometr pokazuje większe odkształcenie niż symulacja, problem może leżeć w modelu materiału, mocowaniu, tłumieniu albo lokalnej koncentracji naprężeń.

Pressure-impulse i SDOF

W wielu analizach konstrukcyjnych sam pik ciśnienia jest niewystarczający. Ważny jest impuls, czyli całka z nadciśnienia po czasie. Dwa obciążenia o tym samym maksimum mogą dać inne skutki, jeśli jedno trwa krócej, a drugie dłużej. Właśnie dlatego istnieją wykresy pressure-impulse i modele SDOF, opisane w osobnym artykule o pressure-impulse i modelu SDOF.

Aparatura pomiarowa musi więc dostarczyć nie tylko punktu maksymalnego, ale całego przebiegu. Jeśli tor pomiarowy obetnie pik, źle policzymy nadciśnienie. Jeśli źle odtworzy ogon dodatniej fazy, źle policzymy impuls. Jeśli nie znamy czasu, nie porównamy odpowiedzi konstrukcji z obciążeniem.

Tensometry i czujniki ciśnienia spotykają się właśnie w tym miejscu: czujnik opisuje wejście do modelu, a tensometr i kamera opisują odpowiedź. Dobrze zrobiony eksperyment pozwala zamknąć bilans: czy obciążenie o takim ciśnieniu i impulsie dało takie ugięcie, jakie przewiduje model?

Aparatura historyczna i współczesna

W Projekcie Manhattan diagnostyka dynamiczna była równie ważna jak teoria. Artykuł o klasycznej elektronice impulsowej ery atomowej przypomina, że konstruktorzy potrzebowali oscyloskopów, kamer smużeniowych, fotografii błyskowej i czujników czasu przyjścia fali. Bez tego układ implozyjny byłby tylko hipotezą.

Współczesna aparatura jest wygodniejsza, ale nie zwalnia z myślenia. FPGA może rozdzielać trigger. Oscyloskop cyfrowy może zapisać miliony próbek. Kamera może rejestrować tysiące lub miliony klatek na sekundę. Miernik tensometryczny może komunikować się przez MODBUS. Ale źle podłączony czujnik, źle dobrane pasmo albo brak synchronizacji nadal psują wynik.

Postęp polega więc nie tylko na lepszym sprzęcie, ale na lepszej rachunkowości pomiaru: każdy kanał ma znane opóźnienie, pasmo, kalibrację i miejsce w geometrii zdarzenia.

Granice bezpieczeństwa

Publiczny artykuł o aparaturze fali uderzeniowej nie powinien opisywać, jak przygotować ładunek, jak dobrać odległości prób poligonowych ani jak przeprowadzić test z materiałem wybuchowym. Takie informacje należą do procedur certyfikowanych laboratoriów i reżimu bezpieczeństwa materiałów niebezpiecznych.

Można natomiast bezpiecznie omawiać:

• strukturę toru pomiarowego,
• wymagania pasma i synchronizacji,
• interpretację przebiegów,
• kalibrację czujników,
• modele numeryczne,
• dane syntetyczne,
• eksperymenty niskonapięciowe i niskoenergetyczne,
• analizę opublikowanych wyników.

To wystarcza do bardzo dobrego artykułu dydaktycznego. Student uczy się, jak działa metrologia zjawisk szybkich, bez wchodzenia w praktykę niebezpiecznych prób. Umiejętność krytycznej oceny budżetu niepewności, doboru pasma i synchronizacji torów pomiarowych to kompetencja transferowalna — przydatna zarówno w metrologii jądrowej, jak i w szerokim spektrum badań inżynierskich.

Dodatkowe materiały multimedialne

Najkrótsze podsumowanie: aparatura do badań fali uderzeniowej jest systemem czasu, geometrii i kalibracji. Czujnik ciśnienia, tensometr i kamera nie konkurują ze sobą, lecz opisują różne aspekty tego samego krótkiego zdarzenia.

Ćwiczenia praktyczne

Pierwsze ćwiczenie katalogowe: student dostaje kartę miernika tensometrycznego i wypisuje, które parametry są istotne dla pomiaru statycznego, a które dla dynamicznego. Ma wskazać, dlaczego rozdzielczość 24 bit nie zastępuje pasma, a czas pomiaru jednego kanału decyduje o tym, czy przyrząd nadaje się do szybkiego impulsu.

Drugie ćwiczenie syntetyczne: wygenerować w arkuszu lub skrypcie prosty impuls ciśnienia Friedlandera albo trójkątny impuls dodatni. Następnie przefiltrować go filtrem dolnoprzepustowym o trzech różnych pasmach i porównać pik, czas narastania oraz impuls. Ćwiczenie pokazuje, jak tor pomiarowy zmienia wynik bez zmiany zjawiska.

Trzecie ćwiczenie niskoenergetyczne: zsynchronizować kamerę telefonu lub szybką kamerę szkolną z rejestratorem prostego czujnika nacisku w bezpiecznym eksperymencie mechanicznym, np. uderzeniu wahadła w belkę sprężystą. Celem jest nie fala uderzeniowa, lecz nauka wspólnego triggera, opóźnienia i porównania obrazu z sygnałem.

Czwarte ćwiczenie modelowe: wziąć opublikowany wykres ciśnienia lub dane syntetyczne i użyć go jako wejścia do prostego modelu SDOF. Następnie porównać przewidywane ugięcie z hipotetycznym odczytem tensometru. Wynikiem ma być lista źródeł niepewności: czujnik, filtr, model materiału, mocowanie i synchronizacja.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Równoważnik trotylowy, nadciśnienie i impuls
• Od nadciśnienia do impulsu: praktyczne modele Kingery-Bulmash, Brode i CONWEP
• Od nadciśnienia do ugięcia płyty: numeryczne modelowanie skutków fali podmuchowej
• Hydrodynamika fal uderzeniowych w broni jądrowej