Příručka pro navrhování ocelových konstrukcí podle AISC 26 – Návrh konstrukcí odolných proti výbuchu [1] a zejména Příklad 2.1 – Předběžné vyhodnocení odolnosti jednopodlažní konstrukce proti výbuchu je ideálním referencí pro inženýry. zjednodušené návrhové zatížení při výbuchu.
Idealizovaný časový průběh tlaku při zatížení výbuchem
Idealizovaný časový průběh tlaku ukazuje, jak se mění tlaková síla v čase po výbuchu.
Některé z nejdůležitějších parametrů jsou v grafu zakresleny přímo:
- Peak overpressure (Pr or Pso) … The instantaneous pressure arriving at the structure above the ambient atmospheric pressure.
- Positive phase duration (td) … The time period for the pressure to return to ambient.
- Positive impulse (I) … The total pressure-time energy applied during the positive duration calculated by the area under the curve.
- Negative phase duration (td-) … The time period following the positive phase where the pressure falls below the atmospheric pressure.
Je třeba poznamenat, že idealizovaný časový průběh tlaku ukazuje dvě různé křivky pro „boční zatížení výbuchem“ (dopadajícím přetlakem) a „zatížení odraženým výbuchem“ (odraženým přetlakem), které jsou vyznačeny čárkovanou resp. plnou čarou. Boční zatížení výbuchem (nazývané také zatížení dopadajícím přetlakem) má v literatuře většinou index „so“ (side-on). Při bočním zatížení postupuje rázová vlna rovnoběžně s povrchem spíše než kolmo na něj. Zatížení se v zásadě přes plochu bez překážek přežene. Příkladem může být boční stěna, která je rovnoběžná se zatížením výbuchem, nebo zadní stěna, která není přímo vystavena výbuchu.
Zatížení odraženým výbuchem označované indexem „r“ nastává tam, kde tlaková vlna naráží na šikmou plochu, která není rovnoběžná. Pro stanovení odrazného tlaku Pr lze použít následující rovnici.
Pr = Cr Pso
Kde Pso je dopadající (boční) přetlak a Cr je odrazný koeficient. Cr je funkcí úhlu dopadu a dopadajícího přetlaku.
Na následujícím obrázku je znázorněno, jak lze vypočítat úhel dopadu se zohledněním počátečního směru rázové vlny a vlny odražené kolmo od plochy.
Once the angle of incidence is determined, Figure 2-193 given in the United Facilities Criteria (UFC) 3-340-02 – Structures to Resist the Effects of Accidental Explosions [2] can be used to provide the Cr value based on the Peak Incident Overpressure value.
Zjednodušený časový průběh tlaku při rázové vlně
Pro posouzení se výše popsané idealizované zobrazení zjednoduší na trojúhelníkové rozdělení s okamžitým vzestupem a lineárním poklesem pod kladnou fázi. Pro zachování hodnot maximálního přetlaku a impulzu (plochy pod křivkou) z idealizovaného průběhu se aproximuje fiktivní doba trvání te jako te = 2(I/P).
Vztahy mezi hmotností nálože, vzdáleností mezi konstrukcí a výbuchem a parametry výbuchu definovanými v časovém průběhu tlaku byly zdokumentovány v rozsáhlých studiích. Technical manuals such as resource [2] include the air blast parameters as a function of the scaled distance in the form of empirical blast parameter curves.
Pro zjednodušení se často u jednoduchých konstrukcí nebere v úvahu záporná fáze, protože má malý dopad na analýzu účinků výbuchu. Záporná fáze se však stává důležitější pro konstrukční prvky slabší ve směru zpětného zatížení nebo pro ty, které mají krátkou základní periodu vůči době trvání zatížení.
Další proměnné, které mohou ovlivnit analýzu účinku výbuchu, jako například odporové síly v důsledku větru nebo dynamického tlaku, stínění přilehlými budovami (redukce zatížení) a odraz (zesílení zatížení) a také vnitřní zatížení vlivem tlakové vlny vstupující do konstrukčních otvorů, nebyly v tomto článku zohledněny.
AISC Design Guide 26 – Example 2.1 in RFEM
AISC Design Guide 26 – Example 2.1 [1] is an ideal reference example to apply the blast load analysis in RFEM which follows the above assumptions. Příkladem konstrukce je jednopodlažní ocelová budova o rozměrech 50 ft (šířka) 70 ft (délka) a 15 ft (výška). Přes šířku budovy jsou v programu RFEM modelovány vyztužené rámy s W-profily válcovanými za tepla, zatímco v podélném směru jsou tuhé rámy také s W-profily. Paždíky a vaznice jsou modelovány pomocí C-profilů válcovaných za tepla. Fasádu budovy tvoří trapézové plechy.
Exploze má hmotnost nálože 500 liber a probíhá mírně nad zemí 50 ft od čela konstrukce. Na základě této informace se vypočítá redukovaná vzdálenost Z podle následující rovnice.
| R | Vzdálenost od prvku k vložce |
| W | Hmotnost ekvivalentní nálože TNT |
Čelní stěna
Using the scaled distance, Figure 2-15 from [2] can be utilized to directly determine the positive blast wave parameters for the reflected and side-on pressure listed below in Table 1.
| Parametry zatížení výbuchem | From Figure 2-15 [2] | vypočítaná hodnota |
|---|---|---|
| Odražený maximální přetlak (+) | Pr = 79.5 psi | - |
| Dopadající maximální přetlak (+) | Pso = 24.9 psi | - |
| Odražený impulz (+) | Ir = 31.0W1/3 | Ir = 246 psi ms |
| Dopadající impulz (+) | Iso = 12.1W1/3 | Iso = 96.0 psi ms |
| Doba dosažení konstrukce | ta = 1.96W1/3 | ta = 15.6 ms |
| Doba exponenciálního zatížení (+) | td = 1.77W1/3 | td = 14.0 ms |
| Rychlost rázové vlny | U = 1.75 ft/ms | - |
Vzhledem k tomu, že čelní strana směřuje přímo k výbuchu, platí pro tuto plochu „odražené“ proměnné z tabulky 1. Pro zjednodušený trojúhelníkový přístup je nutné, aby byla vypočítána ekvivalentní doba trvání tak, aby bylo zajištěno, že impulz (plocha pod křivkou) je pro kladnou fázi zachován.
te,r = 2Ir / Pr = 2(246 psi ms) / 29,5 psi = 6,19 ms
Počáteční časový průběh tlaku je nyní pro čelní stěnu kompletní.
Boční stěny a střecha
Pro zjednodušení se pro výpočet proměnných výbuchu pro boční stěny a střechu budovy použije redukovaná vzdálenost Z vypočítaná pro čelní stěnu. Proto se pro stanovení časového průběhu tlaku pro tyto strany budovy používají "dopadové" hodnoty v tabulce 1 výše. Pro zohlednění redukce rázové vlny v závislosti na vzdálenosti stěny a střechy od výbuchu lze provést podrobnější výpočet.
Ekvivalentní doba te se vypočítá pomocí "dopadových" proměnných.
te,so = 2Iso / Pso = 2(96,0 psi ms) / 24,9 psi = 7,71 ms
Zadní stěna
Redukovaná vzdálenost Z pro zadní stěnu se upraví tak, aby navíc zohledňovala délku budovy. Vzdálenost je nyní 50 ft + 70 ft, což je celkem 120 ft. Z se tak vypočítá následovně.
| R | Vzdálenost od prvku k vložce |
| W | Hmotnost ekvivalentní nálože TNT |
Figure 2-15 from [2] can be utilized again to determine the positive blast wave parameters for the side-on pressure listed below in Table 2.
| Parametry zatížení výbuchem | From Figure 2-15 [1] | vypočítaná hodnota |
|---|---|---|
| Dopadající maximální přetlak (+) | Pso = 4.60 psi | - |
| Dopadající impulz (+) | Iso = 5.54W1/3 | Iso = 44.0 psi ms |
| Doba dosažení konstrukce | ta = 8.32W1/3 | ta = 66.0 ms |
| Doba exponenciálního zatížení (+) | td = 3.11W1/3 | td = 24.7 ms |
| Rychlost rázové vlny | U = 1.26 ft/ms | - |
Ekvivalentní doba trvání pro zadní stěnu te lze vypočítat pomocí příslušných proměnných uvedených výše.
te,so = 2Iso / Pso = 2(44,0 psi ms) / 4,60 psi = 19,1 ms
Vzhledem k tomu, že výška zadní stěny je 15 ft nad úrovní terénu v níž dochází k výbuchu, nedochází k okamžitému nárůstu tlaku. Rather, the velocity of the blast wave, the rear wall height, and time of arrival are used to calculate the time to peak pressure, t².
t2 = L1 / U + ta = 15,0 ft / 1,26 ft/ms + 66,0 ms = 77,9 ms
Nyní je možné stanovit čas do konce zatížení přetlakem tf.
tf = t2 + te,so = 77,9 ms + 19,1 ms = 97,0 ms
Kombinací všech výše vypočítaných proměnných pro zadní stěnu je časový průběh tlaku pro tuto část budovy kompletní.
Shrnutí zatížení výbuchem
Čelní, boční/střešní a zadní stěny lze sestavit dohromady, aby bylo možné zobrazit celkový tlak v čase a ukázat, jak rázová vlna v průběhu času ovlivní různé oblasti konstrukce.
Tyto informace lze nyní zadat do programu RFEM a přídavného modulu RF-DYNAM Pro - Forced Vibrations pro definování časového diagramu.
Použití v programu RFEM
Nyní, když byly definovány časové průběhy tlaku pro různé části budovy, lze tyto informace převzít do přídavného modulu RF-DYNAM Pro - Forced Vibrations v programu RFEM.
Před spuštěním časové analýzy je potřeba v modulu RF-DYNAM Pro - Natural Vibrations stanovit vlastní periody, frekvence a vlastní tvary konstrukce. Tato část analýzy není v tomto článku podrobně rozebrána.
For the time-history analysis, a general area load is applied as three separate load cases in RFEM to emulate the blast load location on the structure including LC1 – Front Wall, LC2 – Side Wall/Roof, and LC3 – Rear Wall. Velikost 1 kip/ft2 se použije pouze jako dočasná hodnota, protože tato hodnota bude později stanovena časovou funkcí.
V modulu RF-DYNAM Pro - Forced Vibrations se časové diagramy definují pro každou oblast konstrukce.
Všimněte si, že každý časový diagram odráží výše stanovené informace, jako je maximální tlak a ekvivalentní doba trvání pro čelní stěnu, boční stěny/střechu a zadní stěnu.
Po zadání časových diagramů se obecná plošná zatížení v programu RFEM přímo spojí s příslušným diagramem.
Před spuštěním analýzy je také třeba v přídavném modulu nastavit další parametry výpočtu, například lineární implicitní Newmarkovu metodu, maximální dobu trvání časové analýzy 0,5 s a časový krok 0,001 s. Kromě toho se v modulu nastaví součinitele Rayleighova útlumu a a β, přičemž se použije úhlová frekvence z obou dominantních vlastních tvarů vypočítaná analýzou vlastních frekvencí spolu s poměrným Lehrovým útlumem 2 %.
Všechny důležité informace pro časovou analýzu výbuchu jsou nyní definovány a lze spustit výpočet v programu RFEM a RF-DYNAM Pro. Pro vyhodnocení odezvy a bezpečnosti konstrukce v průběhu výbuchu lze použít vyhodnocovací nástroje, jako je například sledování časového průběhu v programu RFEM. For a detailed demonstration of AISC Design Guide 26 Example 2.1 [1] in RFEM, refer to the previously recorded webinar Blast Time History Analysis in RFEM.
