Contact Us

Blast Resistant Building Design With Finite Element Analysis

Modular Blast Resistant Building
Modular Blast Resistant Building
Explosive hazards are prevalent in many industries. The source of these explosions can be accidental or malevolent. Regardless, site personnel need protection from these hazards. Blast-resistant buildings shield workers and critical assets from blast effects. In this post, we highlight Finite Element Analysis, and its application in the design of Blast Resistant Buildings.  In this post we’ll cover:
  • Basic blast loading design concepts
  • Differences between hand calculations and Finite Element Analysis
  • Blast analysis using Finite Element Analysis

Blast Loading Design Concepts

Blast overpressure loads have a large magnitude and short duration. Resolving these loads using traditional elastic design methods results the specification of large, costly structural members.
Instead, blast resistant structures are designed using ductility. A ductile structure yields and deforms before failing. This deformation dissipates the energy of the blast. Its’s like the crumple zone of a car, you want the car to absorb the energy not the occupants. Also like a car, the design of these structures is a balance of economy, strength and ductility.

Determining Blast Loads

Two governing factors determine the blast pressure time history for a given location:
  1. Charge Size, W: The amount of energy released by the explosive blast, typically quantified as an equivalent weight of TNT
  2. Standoff, R: Distance between the blast resistant structure and enter of the blast

Blast loading results in a short, but intense positive phase loading followed by a longer, but much less intense negative phase.  The negative phase is often neglected because the loads are much lower.  The pressure vs. time relationship is nonlinear (shown in red in the figure below) and varies with charge size and standoff distance. A conservative, linear pressure history, shown in green, is often used for simplicity.

Free Field Blast Pressure Time History
Free Field Blast Pressure Time History

Selected parameters for the curve come from tables such as Figure 2-15 of UFC 3-340-02 (shown below). Parameters in blast engineering are based on the scaled distance, Z.  Note that parameters and constants are in lbs, psi, m-sec!

    \[ Z=\frac{R}{W^{1/3}} \]

Positive Phase Shock Wave Parameters for TNT on Surface At Sea Level
Positive Phase Shock Wave Parameters for TNT on Surface at Sea Level

Parameters include:

  1. Peak Static Overpressure, P_S_O : Peak blast pressure in the free field (neglecting clearing effects).
  2. Peak Reflected Overpressure, P_r : Peak reflected blast pressure due to structure
  3. Peak Incident Impulse, i_s :  Area under the pressure-time curve, essentially energy carried by the blast wave
  4. Fictitious Positive Phase Duration, t_{of} : The duration of the positive phase of blast over pressure.

The pressure-time history developed above is for the free field. When the blast wave impacts a structure, the wave is reflected.  This reflection amplifies the P_{SO} and is known as Peak Reflected Overpressure, P_r.

Reflected Blast Pressure Time History
Reflected Blast Pressure Time History

Selected parameters for the curve are calculated based on charge size, standoff and projected frontal area of the structure.  Required parameters include:

  1. Clearing Time, t_c :  Time it takes for the reflected blast wave to reach edge of structure from center face
  2. Positive Reflected Impulse, i_r :  Energy carried by the reflected blast wave
  3. Drag Coefficient, C_D : Drag coefficient, typically C_D=1 for front walls of structures
  4. Peak Dynamic Pressure q : Peak dynamic pressure
  5. Fictitious Reflected Positive Phase Duration, t_{rf} : The duration of the reflected positive phase of blast over pressure.

Two pressure time histories are developed:

  • Peak Reflected Pressure History: Developed from the peak reflected overpressure P_r and Fictitious Reflected Positive Phase Duration, t_{rf}.
  • Peak Dynamic Pressure & Drag Contribution: Developed from the summation of peak static pressure P_S_O and drag contributions C_Dq_o along with the Fictitious Positive Phase Duration, t_o_f

The pressure history with the least impulse (area under the triangular load curve) is used for the analysis of the structure.  This is the minimum of  i_r and 1/2*(C_D q_o)*t_{of}

Structural Response

The next step is to determine the structural response.  Traditional elastic analysis is not suitable as blast loads will cause materials to yield.  Blast protection is provided by assuring the structure has the capacity to absorb the blast energy without collapsing.  The peak displacement is checked against the maximum permissible ductility.

First we select an idealized structural system, and determine the ultimate resistance (maximum load) for the system.

Table 3-1 Ultimate Unit Resistances for One-Way Elements
Table 3-1 Ultimate Unit Resistances for One-Way Elements

Next we determine the effective natural period of vibration, T_n the system.  The natural period is determined  as a function of the K_{LM}, the Load Mass Factor and K_E, the equivalent unit stiffness of the system.

    \[ T_n=2\pi\left \{\frac{K_{LM} m}{K_E} \right \}^{1/2} \]

 

Table 3-12 Transformation Factors for One-Way Elements
Table 3-12 Transformation Factors for One-Way Elements

The ductility ratio, \mu, in a SDOF system is defined as the ratio of maximum deflection X_m, to equivalent elastic deflection, X_E.  This can be determined from elasto-plastic response graphs.  One example is Figure 3-54 of UFC 3-340-02. The maximum displacement is a function of  the ratio of applied load to yield load and duration of loading to natural period.  Note that we use the chart for a triangular impulse load to mimic the decay of a blast overpressure.

Maximum Deflection of Elasto-Plastic, One-Degree-of-Freedom System for Triangular Load
Maximum Deflection of Elasto-Plastic, One-Degree-of-Freedom System for Triangular Load

The calculated ductility ratio is then compared to the maximum ductility ratio permitted by the code. These ratios vary by code and material of construction.  For purposes of this example, Table 5-8 of UFC 3-340-02 lists the maximum permitted ductility ratio for different structural elements subjected to blast loads.  If the calculated ductility ratio is less than maximum permitted the structure provides sufficient protection.  If not, the structural response must be re-assessed.

Table 5-8 Summary of Deformation Criteria
Table 5-8 Summary of Deformation Criteria

Hand Calculation Drawbacks

Hand calculations are often used to evaluate blast loads, and they require the simplification of structural systems.  Structural response is determined using Single-Degree of Freedom elasto-plastic analysis, and they’re not sufficient to capture the response of complex structural system. Shortcomings of hand calculations include:
  • Inability to predict localized stresses
  • Failure to predict interaction with adjoining structures
  • Inability to evaluate irregular framing systems.
Finite Element Analysis offers an additional tool when hand calculations aren’t enough.

Finite Element Analysis Advantage

The best way to evaluate the response of a blast resistant structure is with FEA. There are two methods of FEA, implicit and explicit. Explicit FEA models the decay of blast pressure with time. The benefits of explicit FEA blast analysis include:
  • Irregular Geometry: Framing systems around doors and other openings must carry extra load. FEA predicts the stresses on these critical members.
  • Secondary Effects: FEA predicts the resulting accelerations and deformation from blasts. This is used to assess potential internal missile generation following a blast.

Blast Resistant Building Design Example

The example below shows the benefits of using FEA to assess a Blast Resistant Building.

Blast Resistant Building

A modular blast resistant building needs to provide protection for personnel. The postulated blast has a 5msec duration and 5psi peak reflected overpressure.
 
The blast resistant building will be constructed of steel framing with a thin skin. The building is skinned with corrugated steel. Corrugated steel is used instead of flat plate to increase bending moment capacity.
Modular Blast Resistant Building
Modular Blast Resistant Building

Blast Loading

The blast is modeled as an equivalent triangular loading that decays with time. Blast pressures are applied to the exterior surfaces based on orientation.
 

FEA Mesh

Shells and solid elements comprise the FEA mesh of the blast resistant building. The FEA mesh is refined at areas of high stress or geometric irregularity.
Blast Resistant Building with Mesh
Blast Resistant Building with Mesh

Results

Deformation: The peak displacement in the center of the wall is 24 inches (Node 31515); however, this is an absolute displacement.  The un-anchored blast resistant building also slides 17 inches on the ground (Node 16277).  The relative displacement of the wall is 7″. The building needs to be anchored or have have flexible connections that accommodate this displacement.
Blast Resistant Building: Z Displacement
Blast Resistant Building: Z Displacement
The displacement history at two critical points is shown below.  Node 31515 is located in the center of the wall, while node 16277 is located on the exterior structural frame.  The figure shows that the flexible wall quickly deforms under the blast load and then vibrates following the blast load.  As the structure is loaded by the blast it begins to slide.  This momentum is slowly dissipated by friction as the building slides along the ground.
Blast Resistant Building: Displacement Time History
Blast Resistant Building: Displacement Time History
 
Rotations: Blast analysis procedures limit the rotation at support points. The rotation at the top and bottom of the skin can be calculated, knowing that the maximum relative displacement of the wall is 7″ and is structure is 10′ tall.  Support rotation is calculated to be

    \[ \theta_rotation=\arctan\left \{\frac{7inches}{120inches/2} \right \}=6.7 degrees \]

Depending on material, the calculated support rotation is high in terms of allowables in Table 5-8 of UFC 3-340-02.  The wall materials should be reviewed to determine if it is plate steel or cold formed sections.
 
Stresses: Stresses in the skin surpass the yield point of the steel. Increased loading results in plastic hinging. These hinges develop in the skin and structural framing.
Blast Resistant Building: Von-Mises Stress
Blast Resistant Building: Von-Mises Stress
 
Ductility: Plastic strain accumulates following material yield. This plastic strain is a measure of ductility. This strain occurs at areas of plastic hinging. Plastic strain is less than the minimum tensile plastic strain of 0.20 and is acceptable.
Blast Resistant Building: Effective Plastic Strain
Blast Resistant Building: Effective Plastic Strain
 
Local Acceleration: When subjected to the blast, the wall accelerates at a peak of 400G! Any contents anchored to the wall could become missiles. Wall-mounted fixtures need to take this acceleration into account!
Acceleration Time History of Wall of Blast Resistant Building
Acceleration Time History of Wall of Blast Resistant Building

Conclusion

Finite element analysis is the best method for evaluating blast resistant buildings.  Only FEA can predict the localized strains, stresses and structural accelerations in complex situations.

Read More!

Blast loading concepts described in this blog were taken from UFC 3-340-02  “STRUCTURES TO RESIST THE EFFECTS OF ACCIDENTAL EXPLOSIONS.  This document is used by the U.S. Department of Defense to design their facilities. It is approved for public release with unlimited distribution and can be downloaded here. Recommended reading and worked examples include:

  • UFC 3-340-02 Chapter 2: Blast Fragment & Shock Loads
  • UFC 3-340-02 Appendix 2: Blast Fragment & Shock Loads: Sample Problem 2A-3 Surface Burst

 

Let Us Be Your Trusted Partner

Our Professional Engineers are focused on giving our clients outstanding support.  We’ve been doing blast analysis and physical security design for years, using hand calculations and FEA.  We’ve worked with owners, and fabricators to develop economical, safe designs.
  • Scoping: Our team of engineers use FEA to develop economical blast resistant structures. Our scoping process considers cost-saving measures such as offsite fabrication, minimizing welds and protection afforded by existing buildings.
  • Fabrication: Our FEA software communicates directly with CADD geometry and drawings. This communication reduces design and analysis duration. 
  • Evaluating New Hazards: We use FEA to scope modifications of existing structures for greater hazards. Owners can see substantial savings through small modifications to existing structures.
  • Documentation: We understand Quality Assurance. Our team of engineers prepares reports documenting full compliance with code requirements, and our assessments satisfy regulatory and oversight requirements.

Have Questions?

We’re here to help. If you have any questions regarding Blast Analysis, please call us at 585-340-7277. Our FEA consultants are ready to help 24 hours a day, 7 days per week. You can also email us at info@xceed-eng.com, or use the contact us form at the side of this page.
 

Leave a comment

Your email address will not be published. Required fields are marked *

5 thoughts on “Blast Resistant Building Design With Finite Element Analysis”

  • Jeff Gardiner
    • Jeff Gardiner
  • Jeff Gardiner
    • Jeff Gardiner
  • Jeff Gardiner