# Simplified Blast Design According to AISC Steel Design Guide 26

### Technical Article

Blast loads from high energy explosives, either accidental or intentional, are rare but may be a structural design requirement. These dynamic loads differ from standard static loads due to their large magnitude and very short duration. A blast scenario can be carried out directly in a FEA program as a time history analysis to minimize loss of life and evaluate varying levels of structural damage.

AISC Steel Design Guide 26 - Design of Blast Resistant Structures [1] and in particular Example 2.1 - Preliminary Evaluation of Blast Resistance of a One-Story Structure is an ideal reference to guide engineers though a simplified blast design load application.

#### Idealized Blast Load Pressure-Time History

An idealized pressure-time history diagram shows how the pressure force varies over time after the explosion takes place.

A few of the most important parameters are listed directly on the diagram including:

- Peak overpressure (P
_{r}or P_{so}) ... The instantaneous pressure arriving at the structure above the ambient atmospheric pressure. - Positive phase duration (t
_{d}) ... 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 (t
_{d}^{-}) ... The time period following the positive phase where pressure falls below atmospheric pressure.

Notice there are two different curves represented in the idealized pressure-time history diagram including the "side-on blast load" and the "reflected blast load" indicated by the dashed line and solid line, respectively. The side-on blast load (also called free-field blast load) includes the subscript "so" used commonly throughout literature. This indicates where blast load travels parallel to a surface rather than perpendicular. Essentially, the load will sweep over the surface with no impeding objects. An example of this includes a side wall parallel to a blast load or a rear wall with no immediate exposure to the blast.

In turn, the reflected blast load, indicated by the subscript "r", is where the blast wave strikes an angled surfaces other than parallel. To determine the reflected pressure, P_{r}, the following equation can be used.

P_{r} = C_{r} P_{so}

Where, P_{so} is the side-on pressure and C_{r} is the reflection coefficient. C_{r} is a function of the angle of incidence and side-on pressure. The image below demonstrates how the angle of incidence can be calculated when considering the initial blast wave direction and the reflected wave normal to the surface.

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 C_{r} value based on the Peak Incident Overpressure value.

#### Simplified Blast Load Pressure-Time History

For design purposes, the idealized plot shown above is simplified to a triangular distribution with an instantaneous rise and linear decay under the positive phase. To maintain peak overpressure from the idealized plot as well as impulse (area under the curve), a fictitious time duration, t_{e}, is approximated as t_{e} = 2(I/P).

Extensive research to determine the relationship between charge weight, the standoff distance (distance from the structure to the explosion), and the blast parameters defined in the pressure-time plot have previously been carried out. 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.

The negative phase is often ignored for simplification with simple structures as there is little impact from the blast analysis. However, the negative phase becomes increasingly important when the structure's elements are weaker in the reverse loading direction or have a short fundamental period with respect to the load duration.

Additional variables which may have influence on the blast analysis for the purposes of this article have not been taken into consideration such as drag forces due to wind or dynamic pressures, adjacent building shielding (load reduction) and reflection (load amplification), and interior loads due to the blast wave entering the structure's openings.

#### 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. The example structure is a one-story steel building with 50 ft (W) ⋅ 70 ft (L) ⋅ 15 ft (H) dimensions. In the structure's short direction, braced frames are modeled in RFEM as hot rolled W-sections, while in the long direction, rigid frames are also modeled with W-sections. The girts and purlins are modeled with hot-rolled C-sections. The building façade includes ribbed metal panels.

The explosion has a charge weight of 500 lbs and occurs 50 ft from the front face of the structure slightly above ground elevation. With this information, the scaled distance, Z, is then calculated according to the following equation.

$\mathrm{Z}=\frac{\mathrm{R}}{\sqrt[3]{\mathrm{W}}}=\frac{50\mathrm{ft}}{\sqrt[3]{500\mathrm{lb}}}=6.3\frac{\mathrm{ft}}{{\mathrm{lb}}^{1/3}}$

R |
Distance from element to charge |

W |
TNT equivalent charge weight |

Front Wall

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.

Blast Loading Parameter | From Figure 2-15 [2] | Calculated Value |
---|---|---|

Reflected peak pressure (+) | P_{r} = 79.5 psi | - |

Side-on peak pressure (+) | P_{so} = 24.9 psi | - |

Reflected impulse (+) | I_{r} = 31.0W^{1/3} | I_{r} = 246 psi ms |

Side-on impulse (+) | I_{so} = 12.1W^{1/3} | I_{so} = 96.0 psi ms |

Time of arrival | t_{a} = 1.96W^{1/3} | t_{a} = 15.6 ms |

Exponential load duration (+) | t_{d} = 1.77W^{1/3} | t_{d} = 14.0 ms |

Shock front velocity | U = 1.75 ft/ms | - |

Because the front wall is directly facing the initial explosion, the "reflected" variables in Table 1 are applicable to this surface. The simplified triangular approach requires the equivalent duration to be calculated to ensure the impulse (area under the curve) is preserved over the positive duration phase.

t_{e,r} = 2I_{r} / P_{r} = 2(246 psi ms) / 29.5 psi = 6.19 ms

The initial pressure-time plot is now complete for the front wall.

Side Walls and Roof

For simplicity, the scaled distance, Z, calculated for the front wall is used to determine the blast variables for the building's side walls and roof. Therefore, the side-on values in Table 1 above are used to define the pressure-time plot for this section of the building. A more detailed calculation could be carried out to consider the blast wave reduction as a function of the side wall and roof distance from the explosion.

The equivalent duration, t_{e}, is calculated using the side-on variables.

t_{e,so} = 2I_{so} / P_{so} = 2(96.0 psi ms) / 24.9 psi = 7.71 ms

Rear Wall

The scaled distance, Z, for the rear wall is modified to consider the additional length of the building. The distance is now 50 ft + 70 ft for a total of 120 ft. Therefore, Z is calculated as the following.

$\mathrm{Z}=\frac{\mathrm{R}}{\sqrt[3]{\mathrm{W}}}=\frac{120\mathrm{ft}}{\sqrt[3]{500\mathrm{lb}}}=15.1\frac{\mathrm{ft}}{{\mathrm{lb}}^{1/3}}$

R |
Distance from element to charge |

W |
TNT equivalent charge weight |

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.

Blast Loading Parameter | From Figure 2-15 [2] | Calculated Value |
---|---|---|

Side-on peak pressure (+) | P_{so} = 4.60 psi | - |

Side-on impulse (+) | I_{so} = 5.54W^{1/3} | I_{so} = 44.0 psi ms |

Time of arrival | t_{a} = 8.32W^{1/3} | t_{a} = 66.0 ms |

Exponential load duration (+) | t_{d} = 3.11W^{1/3} | t_{d} = 24.7 ms |

Shock front velocity | U = 1.26 ft/ms | - |

The rear wall equivalent duration, t_{e}, can be calculated with the relevant variables above.

t_{e,so} = 2I_{so} / P_{so} = 2(44.0 psi ms) / 4.60 psi = 19.1 ms

Because the rear wall height is 15 ft above the ground elevation where the blast is taking place, there is not an instantaneous rise in pressure. 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_{2}.

t_{2} = L_{1} / U + t_{a} = 15.0 ft / 1.26 ft/ms + 66.0 ms = 77.9 ms

The time to the end of the blast load, t_{f}, can now be determined.

t_{f} = t_{2} + t_{e,so} = 77.9 ms + 19.1 ms = 97.0 ms

Combining all rear wall variables calculated above, the pressure-time plot for this section of the building is complete.

Blast Load Summary

The front, side/roof, and rear walls can be compiled together to display the total pressure vs. time and illustrate how the blast wave will impact the different areas of the structure over time.

This information can now be taken into the RFEM and the RF-DYNAM Pro-Forced Vibrations add-on modules for the time diagram definitions.

#### Application in RFEM

Now that the pressure-time diagrams have been defined for the various sections of the building, this information can be taken into the RF-DYNAM Pro-Forced Vibrations add-on module within RFEM.

RF-DYNAM Pro-Natural Vibrations to determine the structures natural periods, frequencies, and mode shapes is required before running the time-history analysis. This portion of the analysis is not discussed in detail for the purposes of this article.

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. A magnitude of 1 kip/ft2 is used only as a placeholder as this value will later be dependent on the time-history function.

In RF-DYNAM Pro-Forced Vibrations, the time diagrams are defined for each region of the structure.

Notice, each time diagram reflects the information determined above such as the peak pressure and equivalent duration for the front wall, side walls/roof, and rear wall.

Once the time diagrams are defined, the general area loads in RFEM are directly linked to the relevant diagram.

Additional variables must also be set in the add-on module before running the analysis such as the linear implicit Newmark analysis solver, a maximum time of 0.5 seconds for the time-history analysis duration, and a time step of 0.001 seconds to be used in the calculation. Additionally, utilizing the angular frequency from the two dominant modes calculated with the natural frequency analysis along with a Lehr's damping ratio of 2 %, the Rayleigh damping coefficients a and β are also set in the module.

All relevant information is now defined for the blast time-history analysis and the RFEM and RF-DYNAM Pro calculation can be run. Evaluation tools such as the time-course monitor in RFEM can be used to assess the structure's response and safety over the course of the blast explosion. 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.

#### Author

#### Amy Heilig, PE

CEO - USA Office

Sales & Technical Support Engineer

Amy Heilig is the CEO of the USA office located in Philadelphia, PA. In addition, she provides sales and technical support and continues to aid in the development of Dlubal Software programs for the North American market.

#### Keywords

Dynamic Time History Blast Explosion FEA AISC DG 26 UFC 3-340-02

#### Reference

#### Links

Write Comment...

Write Comment...

Contact Us

Do you have further questions or need advice? Contact us via phone, email, chat, or forum, or find suggested solutions and useful tips on our FAQ page, available 24/7.

New

Features of the Modal Analysis Add-on for RFEM 6

The dynamic analysis in RFEM 6 and RSTAB 9 is divided into several add-ons. This article will show you some of the main features of the Modal Analysis Add-on for RFEM 6 / RSTAB 9.

- How do I include the redundancy factor(s), ρ in the ASCE 7 load combinations?
- How can I create a response spectrum according to ASCE 7‑22?
- I need many mode shapes to achieve a sum of 90% effective modal mass. Can I limit the number of mode shapes for generating equivalent loads?
- Is the German National Annex 2021‑07 to Eurocode 8 available in RFEM 5 / RSTAB 8?
- Can I export a response spectrum from RFEM 6 and use it in RFEM 5, for example?

- Where can I adjust the standardization of a mode shape in RFEM 6 / RSTAB 9?
- How can I neglect masses in my modal analysis in RFEM 6 / RSTAB 9?
- How can I analyze the failure of an object, such as a column, in my modal analysis?
- Why do the results in a modal analysis differ between the initial prestress and the surface load?
- How can I display mode shapes in RFEM 6 and RSTAB 9?

Associated Products