Simplified Blast Design According to AISC Steel Design Guide 26

Technical Article on the Topic Structural Analysis Using Dlubal Software

  • Knowledge Base

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 (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 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, Pr, the following equation can be used.

Pr = Cr Pso

Where, Pso is the side-on pressure and Cr is the reflection coefficient. Cr 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 Cr 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, te, is approximated as te = 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.

Front Wall Scaled Distance

Z = RW3 = 50 ft500 lb3 = 6.3 ftlb1/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 ParameterFrom Figure 2-15 [2]Calculated Value
Reflected peak pressure (+)Pr = 79.5 psi-
Side-on peak pressure (+)Pso = 24.9 psi-
Reflected impulse (+)Ir = 31.0W1/3Ir = 246 psi ms
Side-on impulse (+)Iso = 12.1W1/3Iso = 96.0 psi ms
Time of arrivalta = 1.96W1/3ta = 15.6 ms
Exponential load duration (+)td = 1.77W1/3td = 14.0 ms
Shock front velocityU = 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.

te,r = 2Ir / Pr = 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, te, is calculated using the side-on variables.

te,so = 2Iso / Pso = 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.

Rear Wall Scaled Distance

Z = RW3 = 120 ft500 lb3 = 15.1 ftlb1/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 ParameterFrom Figure 2-15 [2]Calculated Value
Side-on peak pressure (+)Pso = 4.60 psi-
Side-on impulse (+)Iso = 5.54W1/3Iso = 44.0 psi ms
Time of arrivalta = 8.32W1/3ta = 66.0 ms
Exponential load duration (+)td = 3.11W1/3td = 24.7 ms
Shock front velocityU = 1.26 ft/ms-

The rear wall equivalent duration, te, can be calculated with the relevant variables above.

te,so = 2Iso / Pso = 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, t2.

t2 = L1 / U + ta = 15.0 ft / 1.26 ft/ms + 66.0 ms = 77.9 ms

The time to the end of the blast load, tf, can now be determined.

tf = t2 + te,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

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...

  • Views 669x
  • Updated 27 October 2021

Contact us

Contact Dlubal

Do you have questions or need advice?
Contact our free e-mail, chat, or forum support or find various suggested solutions and useful tips on our FAQ page.

+49 9673 9203 0

info@dlubal.com

Event Invitation

2022 NASCC: The Steel Conference

Conference 23 March 2022 - 25 March 2022

Event Invitation

International Mass Timber Conference

Conference 12 April 2022 - 14 April 2022

Event Invitation

Structures Congress 2022

Conference 21 April 2022 - 22 April 2022

Effective BIM Workflows Between RSTAB & RFEM and IDEA StatiCa

Effective BIM Workflows Between RSTAB & RFEM and IDEA StatiCa

Webinar 5 August 2021 11:00 AM - 12:00 PM CEST

Glass Design with Dlubal Software

Glass Design with Dlubal Software

Webinar 8 June 2021 2:00 PM - 2:45 PM CEST

Blast Time History Analysis in RFEM

Blast Time History Analysis in RFEM

Webinar 13 May 2021 2:00 PM - 3:00 PM EDT

Timber Structures | Part 2: Design

Timber Beam and Surface Structures | Part 2: Design

Webinar 11 May 2021 2:00 PM - 3:00 PM CEST

Plate and Shell Buckling Utilizing Dlubal Software

Plate and Shell Buckling Utilizing Dlubal Software

Webinar 30 March 2021 2:00 PM - 2:45 PM CEST

CSA S16: 19 Steel design in RFEM

CSA S16:19 Steel Design in RFEM

Webinar 10 March 2021 2:00 PM - 3:00 PM EDT

The Most Common User Errors With RFEM and RSTAB

The Most Common User Errors With RFEM and RSTAB

Webinar 4 February 2021 2:00 PM - 3:00 PM BST

Member design according to ADM 2020 in RFEM

ADM 2020 Member Design in RFEM

Webinar 19 January 2021 2:00 PM - 3:00 PM EDT

Dlubal Info Day

Dlubal Info Day Online | 15 December 2020

Webinar 15 December 2020 9:00 AM - 4:00 PM BST

Stability Design in Steel Construction with RFEM and RSTAB

Stability Design in Steel Construction with RFEM and RSTAB

Webinar 1 December 2020 2:00 PM - 2:45 PM BST

FEM - Troubleshooting and Optimization in RFEM

FEA Troubleshooting and Optimization in RFEM

Webinar 11 November 2020 2:00 PM - 3:00 PM EDT

Soil-Structure Interaction in RFEM

Soil-Structure Interaction in RFEM

Webinar 27 October 2020 2:00 PM - 2:45 PM BST

NBC 2015 Modal Response Spectrum Analysis in RFEM

Webinar 30 September 2020 2:00 PM - 3:00 PM EDT

Documenting Results in the RFEM Printout Report

Documenting Results in the RFEM Printout Report

Webinar 25 August 2020 2:00 PM - 2:45 PM CEST

ACI 318-19 Concrete Design in RFEM

Webinar 20 August 2020 2:00 PM - 3:00 PM EDT

How to Be More Productive Using RFEM

How to Be More Productive Using RFEM

Webinar 7 July 2020 3:00 PM - 4:00 PM CEST

RFEM 5
RFEM

Main Program

Structural engineering software for finite element analysis (FEA) of planar and spatial structural systems consisting of plates, walls, shells, members (beams), solids and contact elements

Price of First License
3,540.00 USD
RFEM 5

Add-on Module

Dynamic analysis of natural frequencies and mode shapes of member, surface, and solid models

Price of First License
1,030.00 USD
RFEM 5
RF -DYNAM Pro - Nonlinear Time History

Add-on Module

Nonlinear dynamic analysis to external excitations

Price of First License
1,120.00 USD
RFEM 5
RF -DYNAM Pro - Forced Vibrations

Add-on Module

Dynamic and seismic analysis including time history analysis and multi-modal response spectrum analysis

Price of First License
1,120.00 USD