8015x
001660
2020-10-14

Dynamic Analysis of Structures Under Blast Loads

Basics

A structural system must be planned and designed in such a way that it bears up against possible actions and impacts beyond its service life and fulfills the required serviceability. In this regard, actions are classified according to their temporal change as follows:

  • Permanent actions (for example, self-weight)
  • Variable actions (for example, live loads, snow, and wind loads)
  • Exceptional actions (for example, a blast or a vehicle impact)

This technical article deals with the extraordinary action of an explosion. An extraordinary action is of short duration and does not occur with any nameable probability. However, it can have significant consequences on the structure's stability.

"An explosion is a" suddenly occurring, extremely rapid "oxidation or decomposition reaction with a sudden increase in temperature and pressure. This leads to a sudden volume expansion of gases and the release of large amounts of energy in a small space (...). A sudden volume expansion causes a blast that can be described for an ideal explosion (emanating from a single point source) by a blast wave model." [1] In addition to the air blast loading, an explosion comes with more impact due to high temperatures or projection (splinters, debris). In this article, the loading of a remote detonation is represented as a pure air blast loading on a structure, but no other effects of the explosion.

Air Blast Loading of Remote Detonation

The air blast load can be displayed schematically as a pressure-time curve (from [2]).

The free air shock wave hits the structure abruptly with peak overpressure. The curve includes an overpressure period acting on the structure until the time period td is reached, which is reduced by an underpressure period until the ambient air pressure is reached. This exponential approach is often simplified to the overpressure region. In this case, a virtual time t~d (t~d < td) can be calculated that describes the approach linearized with the same amount of momentum, but completely neglects the underpressure period.

The relevant input values for calculating the explosion are the distance to the explosion center R as well as the explosive mass as TNT equivalent MTNT. The following formulas refer to the load model developed in [2]. A scaled distance Z is determined from both input values, R and MTNT.

Z = RMTNT3Z > 0.5 (remote detonation)
Z Scaled distance [m/kg1/3] for Z > 2.8
R Distance to explosion center [m]
MTNT Mass of TNT equivalent [kg]

In the following text, the maximum peak overpressure, the positive specific impulse, and the shape coefficient are calculated. The shape factor has a significant influence on the expression of the underpressure period.

p^10 = p0 · 808 · 1 + Z4.521 + Z0.0482 · 1 + Z0.322 · 1 + Z1.352
p10 Maximum peak pressure of remote explosion (Kinney & Graham) [kPa]
p0 Atmospheric pressure under normal conditions (101.3 [kPa])
Z Scaled distance [m/kg1/3]

i+ = 2.1 · RZ2
i+ Positive specific impulse [kPa ms]
R Distance to explosion center [m]
Z Scaled distance [m/kg1/3] for Z > 2.8

α = 1.5 · Z-0.38
α Shape factor
Z Scaled distance [m/kg1/3] for 0.1 < Z < 30

In the next step, the duration of the positive pressure action td as well as the virtual duration of the positive pressure action t~d can be calculated.

td = i+p^10 · α2α - 1 + e-α
td Duration of the positive pressure action
i+ Positive specific momentum [kPa ms]
p10 Maximum peak overpressure of remote explosion (Kinney & Graham) [kPa]
α Shape factor
e Euler's number

t~d =2 · i+p^10
t~d Virtual duration of the positive pressure action
i+ Positive specific momentum [kPa ms]
p10 Maximum peak overpressure of remote explosion (Kinney & Graham) [kPa]

To determine the reflected pressure-time curve, a reflection factor for the overpressure period cr and a reflection factor for the underpressure period c-r are determined. An infinitely perpendicular reflection face is assumed. For details on the values, see [2].

cr = 8 · p^10 + 14 · p0p^10 + 7 · p0
cr Reflection factor overpressure
p10 Maximum peak overpressure of remote explosion (Kinney & Graham) [kPa]
p0 Atmospheric pressure under normal conditions (101.3 [kPa])

cr- = 1.9 · Z - 0.45Z
cr- Reflection factor underpressure
Z Scaled distance [m/kg1/3] for Z > 0.5

Based on all the determined values, the loading in RF -DYNAM Pro - Forced Vibrations can then be represented by means of the loading model for the complete reflected pressure-time curve

pr0(t) = cr · p^10 · φ(t)t  tdcr- · p^10 · φ(t)t > td
pr0(t) Load model for the fully reflected pressure-time diagram
cr Reflection factor overpressure
p10 Maximum peak overpressure of remote explosion (Kinney & Graham) [kPa]
φ(t) Load function (constant/linear/exponential approach)
td Duration of the positive pressure action
cr- Reflection factor underpressure

and any selected loading functions as time diagrams (functions).

p1(t) = p^r0t  td0t > tdp2(t) = p^r0 · (1 - ttd)t  td0t > tdp3(t) = p^r0 · (1 - tt~d)t  t~d0t > t~dp4(t) = p^r0 · (1 - ttd) · e-α · ttd
p1(t) Load function of the constant momentum
p2(t) Load function of the linear momentum
p3(t) Load function of the linear momentum with the virtual time
p4(t) Load function exponential (approach by Friedlander)
t~d Virtual duration of the positive pressure action
td Duration of the positive pressure action
e Euler's number
α Shape factor
pr0 Fully reflected pressure-time diagram

Input in RF -DYNAM Pro - Forced Vibrations

In the add-on module, the loading functions can be entered as time diagrams. Time diagrams can be defined either transiently, periodically, or directly as a function. They excite the structure in a specific position. The position of the load is defined in static load cases. Almost any load type can be entered here. The static load cases are linked to the time diagrams. This happens in the dynamic load cases. The multiplier k is used to determine the final magnitude of the excitation force.

For the following calculations, a remote explosion of MTNT = 1 kg at a distance of R = 10 m is shown. This results in the following values when using the parameterized input.

In the parameter list stored in the RFEM model file, only the values for R and MTNT have to be adjusted. If these are in the value ranges for the scaled distance of 5 < Z < 30, the calculation model presented in [2] can be used.

With the values calculated in the parameter list, the entries for the four displayed time diagrams in the add-on module have been made as follows. As with many numerical programs, the pressure is not applied directly at t = 0 s, but in our example at t = 0.01 s. The use of nested If-functions is useful here to represent the desired functions.

To compare the four functions in one file, four identical subsystems are analyzed in a dynamic load case. A load case is assigned to each subsystem affecting the front surface by 1 kN/m². Moreover, a different time diagram and thus a different load function is assigned to each subsystem.

Finally, the Rayleigh damping of the subsystems is entered, which can be determined from the two dominant mode shapes of the subsystems in the considered direction.

Results

After calculating and determining the results, you can compare the four loading functions and their effects on the subsystems in the file. This article compares acceleration and displacement in the global X-direction only briefly. You can evaluate the results with the graphical user interface in the Result Navigator. Various result values for the calculated time steps can be displayed here. In addition, after analyzing a dynamic load case, you can access the time history diagram by displaying and comparing the further values of the points. The values in the middle of the front surfaces are considered here.

The application of the constant momentum p1(t) shows the largest values, as expected. The two linearized curves p2(t) and p3(t) are very similar, with the values of p2(t) > p3(t) as expected. In the end, the course of p4(t) shows that considering the underpressure period cannot be neglected and that larger values act on the structure compared to the common, linearized approach of p3(t).

Summary

Representing the real pressure-time curve of a remote detonation by means of time diagrams in RF-DYNAM Pro - Forced Vibrations is an effective way to determine the effects of the over- and underpressure periods on the structure. The model's parametrization allows you to represent and compare different blast scenarios by adjusting R and MTNT.


Author

Mr. Hoffmann is responsible for development in the areas of dynamic analysis, membrane structures, and RWIND. Furthermore, he provides technical support for our customers.

Links
References
  1. Lexikon chemie.de: Explosion
  2. Teich, M.: Berichte aus dem Konstruktiven Ingenieurbau - Interaktionen von Explosionen mit flexiblen Strukturen. Neubiberg: Universität der Bundeswehr München, 2012