## Serviços Personalizados

## Journal

## Artigo

## Indicadores

- Citado por SciELO
- Acessos

## Links relacionados

- Citado por Google
- Similares em SciELO
- Similares em Google

## Compartilhar

## Latin American Journal of Solids and Structures

##
*versão On-line* ISSN 1679-7825

### Lat. Am. j. solids struct. vol.11 no.11 Rio de Janeiro 2014

#### https://doi.org/10.1590/S1679-78252014001100006

**Minimum mesh design criteria for blast wave development and structural response - MMALE method**

**Jovan Trajkovski ^{I,}^{*}; Robert Kunc^{I}; Jasenko Perenda^{I,II}; Ivan Prebil^{I}**

^{I}University of Ljubljana, Faculty of Mechanical Engineering, Chair of Modelling in Engineering Sciences and Medicine, Askerceva c. 6, 1000 Ljubljana, Slovenia

^{II}Sistemska Tehnika d. o. o, Koroska cesta 14, 2390 Ravne na koroskem, Slovenia

**ABSTRACT**

The Multi Material Arbitrary Lagrange Euler (MMALE) method is widely used method for numerical investigation of structural response under blast loading. However the method is very demanding for use at the other hand. In this paper are presented the results of the detailed numerical investigation in order to simplify some decisions contributing to the accuracy and efficiency of this model. The influence of mesh properties (particularly mesh size, its biasing and distance of the boundary (DoB) conditions from the deforming structure) on blast wave loading parameters and structural response is investigated in detail and based on the results minimum mesh design criteria is proposed. The results obtained are presented as a function of the scaled distance and additionally related to the radius of the charge. Validation studies were also done successfully.

**Keywords:** Blast wave, blast response, blast loading, Fluid-Structure interaction

** INTRODUCTION**

Examination of blast loading by explosive devices is of great importance in both, design and protection of different military equipment. Beside this, probability of terroristic attack and exposure to blast loading of public structures like bridges, metros, embassies, and other governmental buildings, is dramatically increased especially in the last two decades. Therefore, it is necessary to improve the response of such structures to this high-intensity, short-duration loads. However, preparation and performance of blast loading experiments using high energetic materials is extremely expensive and time consuming even for small scale tests in laboratory conditions, making the numerical analysis the most valuable tool for examination. On the other hand, precise modeling of blast loading of structures is one of the most complex dynamic analyses. However, considerable effort has been made in the last two decades in the field of both, numerical approximation of blast loads and numerical response of structures to this kind of loads. Being two separated branches in the past they are now interactively coupled in the most advanced finite element solvers like Abaqus, AUTODYN, LS-DYNA.

The simplest numerical approximations of blast loads are triangular, rectangular, or bilinear pressure-time curves, more often used in the past but also at present (Biggs (1964), Jama* et al. *(2009), Krauthammer and Ku (1996), Louca* et al.* (1996), Tavakoli and Kiakojouri (2014)), and in analytical methods (Biggs (1964), Jones (1997), Smith and Hetherington (1994)). More precise, build-in numerical representation of blast loads appeared when Randers-Pehrson and Bannister (1997) implemented previous work of Kingery and Bulmash (1984) (ConWep blast model Bruce and Jon (1991)) in LS-DYNA (LSTC (2013a, b)). Shortly after that authors started using it for simulation in different scenarios, comparing their experimental results with the numerical. Xu and Lu (2006) applied ConWep to study the spallation of reinforced concrete, Neuberger* et al. *(2007a), to examine scaling of air blast loaded plates, Børvik* et al.* (2009) investigated the blast response of ISO container, Hussein* et al.* (2011) investigated the response of Retrofitted RC columns. Nearly at the same time, a new Arbitrary Lagrange Euler (ALE) method for modeling explosions became widely available, attacking the reality of the physical process itself. This method allows explosive and air to be represented as separated materials with a possibility to shear the same elements in the model (Alia and Souli (2006), Souli* et al.* (2000)). All the processes, starting with the detonation, blast wave formation and its interaction with the surrounding structures through Fluid Structure Interaction (FSI) algorithm were joined together in this model. This contributed to a more precise approximation of blast loads allowing more complicated real scenarios to be modeled and more optimized design to be achieved. Using ALE method, Zhao* et al.* (2012) investigated the response of reinforced concrete containment of an nuclear power plant under internal blast loading, Neuberger* et al.* (2007a, b) applied it to study armor plate response under air blast and buried charges detonation. Zakrisson et al. (2012), Zakrisson* et al.* (2011) investigated the plate response under blast loading modeling cylindrical explosive geometry in a steel pot according to NATOs standard AEP-55) requirement. Soutis* et al.* (2011) studied the response of GLARE panels to air blast loading using both, ALE and ConWep methods. In other studies ALE method was also used in variety of different applications and purposes: Chafi* et al. *

(2009), Chung Kim Yuen* et al.* (2012), Fox* et al.* (2011), Jayasinghe* et al.* (2013), Langdon* et al. *(2010), Liu* et al.* (2012), Ma* et al.* (2013), Pi* et al.* (2012), Spranghers* et al.* (2013), Tai* et al. *(2011). In order to efficiently represent FSI on higher scaled distances, other techniques like coupling of the empirical ConWep model with ALE (Todd (2010)) or ALE mapping techniques are used (Luccioni* et al.* (2006), Zakrisson* et al.* (2011)). Some simplifications based on compressibility effect in FSI were also proposed by Kambouchev* et al.* (2006) and applied by Lin* et al. *(2013). Few authors (Barsotti* et al.* (2012), Genevieve and Robert (2008), Wang* et al.* (2005), Xu and Liu (2008)) used Smooth Particle Hydrodynamics (SPH) method to model blast response. The corpuscular approach was also proposed (Olovsson* et al.* (2010)) and successfully applied (Børvik* et al.* (2011)) in blast response scenario. However, MMALE method remains one of the most used methods for blast response analysis. Despite the fact that the accuracy of the numerical results is strongly influenced by the geometrical mesh properties, quite large differences can be observed in the above mentioned studies. These can lead to unreliable results or unnecessary large and nonefficient models on the other hand. As reported in Zukas and Scheffler (2000), given the same task of air blast loaded silo door, four users provided quite different results (reaching even 80% difference among each other) for maximum displacement. Therefore, special care should be given on the most influential numerical parameters in order achieve betters uniqueness and deliver proper loading to the structure while keeping the model cost effective.

In this paper, explosion of 6 kg TNT in free air was simulated using empirical ConWep model, axisymmetric MMALE, and 3D MMALE models using different mesh sizes. The results for incident blast wave parameters corresponding to close range scaled distances of Z=(0.11-0.385) m/kg^{1/3} were compared with the experimental results available in US Army Defense Design Manual UFC-3-340-02 (UFC (2008)) (previously known as TM-5 1300). Numerical asymptotic values of the incident blast wave parameters are presented as a function of number of elements per radius of the charge. After the validation of free air model, the same was upgraded adding a steel plate. This enabled FSI to appear, allowing reflected blast wave parameters to be tracked. The examination was continued further, investigating the influence of DoB conditions from the loaded structure, and the influence of mesh size and its biasing on the reflected blast loading parameters, and consequently on structural displacement. Based on the results, recommendations for more efficient modeling of blast wave phenomena using MMALE method are also given.

**2. BLAST WAVE PARAMETERS **

When explosion occurs in free air, the gaseous products are rapidly expanding out of the initially occupied volume, creating a blast wave. This incident wave has initial velocity close to the detonation velocity of the explosive (6-10 km/s) (Zukas and Walters (2002)). As it is shown on Figure 1, the incident pressure (also called side on pressure) instantly arrives at a certain distance in space at arrival time t_{A}, with its peak value P_{SO}, after which exponentially decreases. The area under pressure-time curve represents the specific impulse from which only the positive one is considered responsible for structural deformation. The negative specific impulse is often neglected due to its small size.

The simplest form of blast wave in free air is described with Friedlander wave equation (1);

in which P_{0} is the atmospheric pressure, P_{SO} is the peak positive pressure, t is current time and t_{d} is time duration of positive pressure.

If the blast wave strikes an object on its way, it reflects and it delivers reflected pressure PR, which could be twice to eight times stronger than the incident one. This appears because the particles at the front of the blast wave are stopped by the structure but they are still forced to move forward by the particles coming from behind.

**3. NUMERICAL MODELS **

Explosion of 6 kg spherical in shape trinitrotoluene (TNT) in free air was modeled using the explicit code LS-DYNA, with three different numerical models, pure Lagrangian model (ConWep), axisymmetric MMALE model, and 3D MMALE model. Blast wave parameters (pressure and specific impulse) in all models were tracked at each 0.1 m distance starting from 0.2 to 0.7 m. This corresponds to close range scaled distance 0.11 m/kg^{1/3} < Z < 0.385 m/kg^{1/3}, where Z=R/m^{1/3}, "R" is the distance from the center of the charge to the target (stand-off distance) and "m" is the equivalent mass of TNT charge.

**3.1. Pure Lagrangian model (ConWep) **

The simplest numerical way of modeling and investigating design of components under blast loading is pure Lagrangian model. With this model only the component under investigation is modeled with Lagrangian element formulation and the ConWep function is called to apply the blast pressure over the elements. The ConWep function (2) is based on blast loading equations proposed by Kingery and Bulmash (1984), which are nothing else but empirical approximation of experimental results.

In this equation *θ* is the angle of incidence of the blast wave relative to the reflecting surface. This is numerically the cheapest method for evaluation of blast response of structures. Main drawback is that the method is unable to represent the interaction between the blast wave and the structure when complicated geometries are investigated. This means that the blast wave cannot be focused due to geometry and also the main structure cannot be "hidden" by the sacrificial one, if it is not physically connected to it. Another limitation of this method is that the empirical equations underlying the spherical air blast are valid for the range of scaled distance 0.147 m/kg^{1/3} < Z < 40 m/kg^{1/3} (LSTC (2013a)).

**3.2. Arbitrary Language Euler (ALE) models **

Lagrangian formulation becomes almost useless when very large deformation processes are simulated (Souli* et al.* (2000)), leading to a very small time steps and possible solution failure. For such cases an alternative ALE method was developed in which the mesh does not necessary follows the material. In this method Lagrangian motion of the nodes is computed at every time step, followed by a possible advection stage in which the mesh is either not advected (pure Lagrangian formulation), advected to the original shape (pure Eulerian formulation), or advected to some more advantageous shape (somewhat between Lagrangian and Eulerian). In the current model the explosive charge and the air were modeled with MMALE formulation in which one element can possible contain both ALE materials. Advection between the steps in all MMALE models was controlled using the modified Van Leer advection algorithm. Structural material in the model was represented with Lagrangian formulation. The coupling between ALE fluids and Lagrangian plate was defined with penalty-based FSI algorithm which conserves energy.

**3.3. Geometry and mesh **

In order to study the influence of mesh size on the incident blast loading parameters the square box model with size of 0.9 m was represented with axisymmetric 2D and 3D models (Figure 2) with different mesh sizes of: 25, 50, 100, 200, 400, 800, and 1600 elements per side (corresponding to element lengths of: 36, 18, 9, 4.5, 2.25, 1.125, 0.5625 mm) for 2D axisymmetric model and 25, 50, and 100 elements per side for 3D model. The charge radius of 96.58 mm (6 kg of TNT) was represented with approximately 3, 6, 11, 22, 43, 86, and 172 elements with the finest regular mesh, leading to a massive models of over 2.5 million elements. Despite the fact that spherical mesh better represents the blast wave in free air (Kakogiannis* et al.* (2010)), rectangular mesh was used for all the models because it is better in solving the leakage problem when FSI is involved (Olovsson and Souli ), Swee* et al.* (2012)) (at least in the particular case of flat plate, placed in parallel with the mesh of the fluids). Spherical shape of the explosive was represented using the initial volume fraction geometry card available in LS-DYNA. The fluids in all 3D models were represented with eight-node brick elements with one integration point and default viscous form of hourglass control, while the structural plate was modeled with eight-node fully integrated solid elements. In the axisymmetric MMALE model axisymmetric area weighted element formulation was used. Boundary conditions as depicted on Figure 2 were applied in all MMALE models.

The reflected blast wave parameters were investigated for stand-off distance of 0.3 m (scaled distance Z=0.165 m/kg^{1/3}) using one-fourth symmetry 3D models. To examine in more detail the influence of DoB conditions on structural deformation, and to relate the minimum required distance at which the boundary conditions can be placed with some physical property of the model, two set of simulations were performed. First set was performed at each 0.1 m stand-off distance from 0.1 m to 0.7 m while keeping the mass of the charge (m = 6 kg) constant. In the second set of simulations the stand-off distance was kept constant (R = 0.3 m) for different mass values of the explosive charge (m = 0.4, 1, 8 and 24 kg). An example representing the series of models for case of R = 0.3 m and m = 8 kg is shown on Figure 3.

After the selection of the proper DoB conditions, simulations were performed with different mesh densities of: 20, 25, 50, and 100 elements per side (≈ 3, 4, 8, and 16 elements per radius of the charge) to study the mesh size influence. In each model, plate sides were represented with two Lagrangian elements per one Eulerian, while five elements through the thickness (t = 10 mm) were kept for each model. Mesh biasing influence was examined using 3D model of 25 elements per side and three different mesh biasing factors of 2, 4, and 6. Finally, the most optimal values for DoB conditions, mesh size and mesh biasing were used in structural response validation study.

**4. MATERIAL MODELS **

**4.1. Explosive **

The explosive charge in all ALE models, was represented with

*MAT_HIGH_EXPLOSIVE_BURN in combination with Jones-Wilkins-Lee (JWL) equation of state (EOS);

which calculates the blast pressure as a function of relative volume *ν* = *ρ*_{0}/*ρ*, and internal energy E, for an explosive element. In this equation A, B, R_{1}, R_{2}, and ω are parameters related to the explosive material and can be found in most of the explosive textbooks. They were taken from reference (Zukas and Walters (2002)) for TNT high explosive and are given in Table 1.

**4.2. Air **

The air is best approximated as an ideal gas for which *MAT_NULL was used in combination with linear polynomial equation of state;

in which *Υ* =** 1.4** is the ratio of the specific heats, *ρ* is current density, ρ0 =1.29 kg/m^{3} is initial density and E=250 kJ represents initial internal energy of the air at atmospheric pressure of 1 bar.

**4.3. Structural material **

Structural plate (diameter D = 0.5 m, plate thickness t = 10 mm) in all the models is represented with bilinear material model and parameters representing RHA steel were taken from Neuberger *et al.* (2007a).

**5. RESULTS **

**5.1. Incident (Side on) blast wave parameters **

The change of incident blast wave parameters (pressure and specific impulse) with change of mesh density, are presented on Figure 4, for four different stand-off distances using axisymmetric MM-ALE model. For better comparison, corresponding experimental values taken from US Army defense design manual (UFC (2008)) and empirical ConWep data from the simulations are also added. It can be noticed that numerical values are approaching the experimental ones as the mesh is refined. It can also be noticed that the empirical model in LS-DYNA overestimates blast loading parameters for scaled distances Z < 0.147 m/kg^{1/3}.

Pressure contours for three different mesh sizes at the same computational time are shown on Figure 5. It can be observed that as the mesh is refined the blast wave is getting its spherical shape despite the rectangular shape of the mesh. It is also visible on the figures with coarser mesh (Figure 5a and 5b) that pressure contours have maximum values along the diagonals and minimum values along the axes. As the mesh is refined, the pressure reaches nearly the same value in all directions (Figure 5c). This is more precisely shown on Figure 6.

Comparison of the axisymmetric MM ALE model and 3D MM ALE model with the experimental results is shown on Figure 7. The results are here presented as a function of number of elements per radius of the charge. It is clear from Figure 7 that we need considerably refined mesh in order to reach experimentally obtained peak pressures while the values of specific impulses can be approached using coarser meshes (See Figure 7). However, in impulsive loading regime the specific impulse is responsible for structural deformation and not the peak pressure (Westine* et al. *(1975)), which allows us more efficient representation of even complicated 3D models on personal computers. Number of elements per radius of the charge needed to achieve a value of 90% of the exact numerical solution for both, incident pressure and impulse is presented on Figure 8 as a function of scaled distance. The values of the exact numerical solutions were estimated using the finest three solutions according to Grid Convergence Index (GCI) method given by Roache (1998). It can be observed from the same figure that depending on the scaled distance, various mesh densities are needed to achieve the 90% of the exact numerical solution. For the specific impulse 65 elements per charge radii were needed near the surface of the charge, while the number of elements decreased to 4 at the scale distance of 0.165 m/kg^{1/3} and then stabilized at nearly 17 elements per charge radii for higher scaled distances.

**5.2. Reflected blast wave parameters **

Beside the material parameters and advection algorithm scheme, mainly the mesh size has the greatest influence on incident blast wave parameters. However, when some structure is investigated under blast loading, additional parameters like DoB conditions of ALE domain, FSI algorithm, and also the ratio between Lagrangian and ALE mesh sizes can have considerable influence on the reflected loading parameters, and consequently on the structural deformation.

**- Influence of the DoB conditions **

on Figure 9 are shown the results for reflected pressure and reflected specific impulse profiles for the most loaded element located at the corner of the one-fourth symmetry models. It is clearly visible that the DoB conditions have great influence on blast loading parameters especially when they are too close to the object of interest (cases 0.325 and 0.4 m). As we move the boundary conditions away from the plate they are losing the influence and after 0.6 m box case we are getting the same results for pressure and specific impulse (See Figure 9). This appears because ambient pressure applied at the boundaries, influences the FSI and makes it unstable, when boundaries are close to the FSI location.

**- Mesh size influence**

With defined DoB (0.6 m box model) simulations were performed with different mesh sizes of: 20, 25, 50, and 100 elements per side (≈ 3, 4, 8, and 16 elements per radius of the charge) to study the mesh size influence on reflected blast wave parameters. The simulation results are presented on Figure 10. It is clear that reflected specific impulse approaches the experimental value as the mesh is refined. However, it should be noted here, that the curves are not wholly comparable, because the area of the corner element considerably differs among the models.

**- Bias influence **

On Figure 11 mesh biasing influence on the blast loading parameters is presented. The results for three different mesh biasing values of 2, 4, and 6 are presented for 3D model with 25 elements per side (4 el./radius). Additionally the results from 25 and 50 el./side models without biasing are given for better comparison. It can be observed that biasing has also great influence on blast loading parameters and can considerably lower the size of the numerical model. However, biasing values higher than four can decrease the time step. Additionally it can lead to negative influence on FSI algorithm causing instabilities and even simulation breakdown.

**5.3. Structural deformation**

Numerical factors that have influence on blast loading parameters will consequently influence the deformation of the loaded object. On Figure 12a, time-displacement curves are presented for the models with different DoB conditions from the plate. Similar to blast loading parameters, timedisplacement curves are changing until 0.6 m box case while no difference is visible with their further moving away (See Figure 12a). On Figure 12b, the influence of mesh density on structural displacement of the plate is presented. It can be observed here, that structural displacement curves are not in consistence with the blast loading parameters in this case. This is due to area differences of the centrally loaded elements. Maximum displacement is growing as the mesh is refined from 3 to 8 el./radius of the charge while no significant difference is visible with further mesh refinement. On Figure 12c) mesh biasing influence on structural deformation is presented.

The results for minimum required distance from both sets of simulations (described in part 3.3.) are related to the radius of the charge and plotted on Figure 13 as a function of the scaled distance. It can be observed that the results from both sets are quite different, which confirms the complexity of numerous parameters that affect the minimum value of the DoB conditions. However, it can be concluded that DoB conditions higher than four times the radius of the charge, can assure a model in which boundary conditions have no influence on the structural response.

**5.4. Structural response validation**

Beside validation of the blast loading parameters, validation study of the structural response was also performed. For that purpose numerical simulations were done and the results compared with the experimental tests results published by Neuberger* et al.* (2007a). Three experimental cases with different close range stand-off distances (see Table 2) were simulated, and the results for maximum normalized deflection compared. Numerical models were prepared as explained in part 3.3, applying the findings of the study; MMALE materials were represented with a box of size R+4R_{TNT} , meshed with 10 elements per radius of the charge, additionally refined toward the center of the charge with a biasing factor of 4. On Figure 14, numerical normalized deflections for all three cases are plotted versus time and corresponding experimental values are added for better comparison.

**6. APPLICATION TO OTHER SCENARIOS **

In reality, more complicated structures like sandwich plates, military vehicles and civil buildings are often analyzed under blast loading. So the decision where to put the boundary conditions and is it necessary to cover all the structure within MMALE can be challenging. To answer that question two different structures representing a sandwich plate (Figure 15a) and simplified vehicle chassis (Figure 16a) were analyzed under blast loading of 8kg TNT charge and constant stand-off distance 0.3m, while changing the DoB conditions. The results are presented on Figure 15b) and on Figure 16b). It can be observed that in both cases the minimum required DoB conditions is 0.4m. The difference is that the sandwich plate has to be all covered with MMALE materials, thus minimum DoB conditions should be measured from the upper plate (Figure 15a) while only a part of the structure can be covered in the case of simplified vehicle (or building structure) and the DoB conditions should be measured from the bottom plate of the vehicle chassis. This is due to the mutual proximity and also connectivity between different parts of particular structure. In the case of sandwich plate, both plates are locally deformed. In the case of blast response of simplified vehicle chassis, the bottom plate is locally deformed while the roof experiences only the global movement of the vehicle.

**7. SUMMARY AND CONCLUSION**

In this paper the influence of the numerical parameters on blast wave loading and structural deformation was investigated in detail. Numerical simulations with different mesh densities were performed for spherical TNT charge in free air. As a result the asymptotic values of incident pressure and specific impulse were obtained for different scaled distances and compared with the experimental results. From these results the optimal number of elements per charge radius was also determined for different scaled distances. The investigation was continued examining the influence of DoB conditions from the object of interest, mesh size, and mesh biasing on the reflected blast loading parameters and structural displacement. It was find out that combination of these parameters has great influence on the results as well as on the size and efficiency of the numerical model. The proper combination of these parameters can greatly increase the accuracy and decrease the computational time. According to the results obtained minimum mesh design criteria is proposed. Models with 10 el./radius of the charge, boundary conditions placed away from the deforming structure at a minimum distance of four charge radii (4R_{TNT} ) additionally refined at the place of the charge with mesh biasing factor of 3 to 4 are recommended. This should yield reasonably accurate results in a relatively short computational time.

**Acknowledgements**

The authors would like to acknowledge SISTEMSKA TEHNIKA d.o.o. for the financial support and their collaboration during the project.

**References **

AEP-55. Procedures for evaluating the protection level of logistic and light armored vehicles (1st ed.), vol. 2 NATO, 2006. [ Links ]

Alia, A., Souli, M. (2006). High explosive simulation using multi-material formulations, Applied Thermal Engineering, 26(10),1032-1042. [ Links ]

Barsotti, M.A., Puryear, J.M.H., Stevens, D.J., Alberson, R.M., McMahon, P. (2012). Modeling Mine Blast with SPH, 12th International LS-DYNA User Conference, Detroit, USA. [ Links ]

Biggs, J.M. (1964). Introduction to structural dynamics, McGraw-Hill. [ Links ]

Børvik, T., Hanssen, A.G., Langseth, M., Olovsson, L. (2009). Response of structures to planar blast loads - A finite element engineering approach, Computers & Structures, 87(9-10),507-520. [ Links ]

Børvik, T., Olovsson, L., Hanssen, A.G., Dharmasena, K.P., Hansson, H., Wadley, H.N.G. (2011). A discrete particle approach to simulate the combined effect of blast and sand impact loading of steel plates, Journal of the Mechanics and Physics of Solids, 59(5),940-958. [ Links ]

Bruce, R.P., Jon, E.W. (1991). CONWEP - Conventional Weapons Effects Prediction Evaluation Test Series, United States. Army. Corps of Engineers, [ Links ].

Chafi, M.S., Karami, G., Ziejewski, M. (2009). Numerical analysis of blast-induced wave propagation using FSI and ALEmulti-material formulations, International Journal of Impact Engineering, 36(10-11),1269-1275. [ Links ]

Chung Kim Yuen, S., Langdon, G.S., Nurick, G.N., Pickering, E.G., Balden, V.H. (2012). Response of V-shape plates to localised blast load: Experiments and numerical simulation, International Journal of Impact Engineering, 46(0),97-109. [ Links ]

Fox, D.M., Huang, X., Jung, D., Fourney, W.L., Leiste, U., Lee, J.S. (2011). The response of small scale rigid targets to shallow buried explosive detonations, International Journal of Impact Engineering, 38(11),882-891. [ Links ]

Genevieve, T., Robert, D. (2008). Finite element simulation using SPH particles as loading on typical Light Armoured Vehicles, 10th International LS-DYNA users conference. [ Links ]

Hussein, M.E., Tarek, H.A., Husain, A., Yousef, A.A.-S., Saleh, H.A. (2011). Effect of blast loading on CFRP- Retrofitted RC columns - A Numerical Study, Latin American Journal of Solids and Structures, 8(1),55-81. [ Links ]

Jama, H.H., Bambach, M.R., Nurick, G.N., Grzebieta, R.H., Zhao, X.L. (2009). Numerical modelling of square tubular steel beams subjected to transverse blast loads, Thin-Walled Structures, 47(12),1523-1534. [ Links ]

Jayasinghe, L.B., Thambiratnam, D.P., Perera, N., Jayasooriya, J.H.A.R. (2013). Computer simulation of underground blast response of pile in saturated soil, Computers & Structures, 120(0),86-95. [ Links ]

Jones, N. (1997). Structural Impact, Cambridge University Press. [ Links ]

Kakogiannis, D., Van Hemelrijck, D., Wastiels, J., Palanivelu, S., Van Paepegem, W., Vantomme, J., Kotzakolios, A., Kostopoulos, V. (2010). Assessment of pressure waves generated by explosive loading, CMES - Computer Modeling in Engineering and Sciences, 65(1),75-93. [ Links ]

Kambouchev, N., Noels, L., Radovitzky, R. (2006). Nonlinear compressibility effects in fluid-structure interaction and their implications on the air-blast loading of structures, Journal of Applied Physics, 100(6), [ Links ] -.

Kingery, C.N., Bulmash, G. Airblast parameters from TNT spherical air burst and hemispherical surface burst. US Army, Armament Research adn Development Center, Ballistic Research Laboratory, Aberdeen Proving Ground, USA; 1984. [ Links ]

Krauthammer, T., Ku, C.K. (1996). A hybrid computational approach for the analysis of blast resistant connections, Computers & Structures, 61(5),831-843. [ Links ]

Langdon, G.S., Rossiter, I.B., Balden, V.H., Nurick, G.N. (2010). Performance of mild steel perforated plates as a blast wave mitigation technique: Experimental and numerical investigation, International Journal of Impact Engineering, 37(10),1021-1036. [ Links ]

Lin, J., Fei, Y., Zhihua, W., Longmao, Z. (2013). A numerical simulation of metallic cylindrical sandwich shell subjected to air blast loading, Latin American Journal of Solids and Structures, 10(3),631-645. [ Links ]

Liu, X., Tian, X., Lu, T.J., Zhou, D., Liang, B. (2012). Blast resistance of sandwich-walled hollow cylinders with graded metallic foam cores, Composite Structures, 94(8),2485-2493. [ Links ]

Louca, L.A., Punjani, M., Harding, J.E. (1996). Non-linear analysis of blast walls and stiffened panels subjected to hydrocarbon explosions, Journal of Constructional Steel Research, 37(2), 93-113. [ Links ]

LSTC. (2013a). LS-DYNA Keyword User's Manual, Volume I, [ Links ]

LSTC. (2013b). LS-DYNA Keyword User's Manual, Volume II, material models, [ Links ]

Luccioni, B., Ambrosini, D., Danesi, R. (2006). Blast load assessment using hydrocodes, Engineering Structures, 28(12), 1736-1744. [ Links ]

Ma, L., Xin, J., Hu, Y., Zheng, J. (2013). Ductile and brittle failure assessment of containment vessels subjected to internal blast loading, International Journal of Impact Engineering, 52(0),28-36. [ Links ]

Neuberger, A., Peles, S., Rittel, D. (2007a). Scaling the response of circular plates subjected to large and close-range spherical explosions. Part I: Air-blast loading, International Journal of Impact Engineering, 34(5),859-873. [ Links ]

Neuberger, A., Peles, S., Rittel, D. (2007b). Scaling the response of circular plates subjected to large and close-range spherical explosions. Part II: Buried charges, International Journal of Impact Engineering, 34(5),874-882. [ Links ]

Olovsson, L., Hanssen, A.G., Børvik, T., Langseth, M. (2010). A particle-based approach to close-range blast loading, European Journal of Mechanics - A/Solids, 29(1),1-6. [ Links ]

Olovsson, L., Souli, M. ALE and Fluid-Structure Interaction Capabilities in LS-DYNA, Pi, S.J., Cheng, D.S., Cheng, H.L., Li, W.C., Hung, C.W. (2012). Fluid-structure-interaction for a steel plate subjected to non-contact explosion, [ Links ] Theoretical and Applied Fracture Mechanics, 59(1),1-7. [ Links ]

Randers-Pehrson and Bannister, K.A. (1997). Airblast Loading Model for DYNA2D and DYNA3D, Defense Technical Information Center. [ Links ]

Roache, P.J. (1998). Verification and Validation in Computational Science and Engineering, Hermosa Publishers. [ Links ]

Smith, P.D., Hetherington, J.G. (1994). Blast and ballistic loading of structures, Butterworth-Heinemann. [ Links ]

Souli, M., Ouahsine, A., Lewin, L. (2000). ALE formulation for fluid-structure interaction problems, Computer Methods in Applied Mechanics and Engineering, 190(5-7),659-675. [ Links ]

Soutis, C., Mohamed, G., Hodzic, A. (2011). Modelling the structural response of GLARE panels to blast load, Composite Structures, 94(1),267-276. [ Links ]

Spranghers, K., Vasilakos, I., Lecompte, D., Sol, H., Vantomme, J. (2013). Numerical simulation and experimental validation of the dynamic response of aluminum plates under free air explosions, International Journal of Impact Engineering, 54(0),83-95. [ Links ]

Swee, H.T., Jiing, K.P., Roger, C., David, C. (2012). Retrofitting of Reinforced Concrete BEam-Column via Steel Jackets against Close-in Detonation, 12 International LS-DYNA Users Conference 2012. [ Links ]

Tai, Y.S., Chu, T.L., Hu, H.T., Wu, J.Y. (2011). Dynamic response of a reinforced concrete slab subjected to air blast load, Theoretical and Applied Fracture Mechanics, 56(3),140-147. [ Links ]

Tavakoli, H.R., Kiakojouri, F. (2014). Numerical dynamic analysis of stiffened plates under blast loading, Latin American Journal of Solids and Structures, 11(2),185-199. [ Links ]

Todd, P.S. (2010). A coupling of empirical explosive blast loads to ALE air domains in LS-DYNA^{®}, IOP Conference Series: Materials Science and Engineering, 10(1),012146. [ Links ]

UFC. Unified Facilities Criteria UFC 3-340-02, Structures to resist the effects of accidental explosions. U.S. Army Armament Research, Development and Engineering Centre, Armament Engineering Directorate; 2008. [ Links ]

Wang, Z., Lu, Y., Hao, H., Chong, K. (2005). A full coupled numerical analysis approach for buried structures subjected to subsurface blast, Computers & Structures, 83(4-5),339-356. [ Links ]

Westine, P.S., Baker, W.E., TX., S.R.I.S.A. (1975). Energy Solutions for Predicting Deformations in Blast-Loaded Structures, Defense Technical Information Center. [ Links ]

Xu, J.-x., Liu, X.-l. (2008). Analysis of structural response under blast loads using the coupled SPH-FEM approach, J. Zhejiang Univ. Sci. A, 9(9),1184-1192. [ Links ]

Xu, K., Lu, Y. (2006). Numerical simulation study of spallation in reinforced concrete plates subjected to blast loading, Computers & Structures, 84(5-6),431-438. [ Links ]

Zakrisson, B., Häggblad, H.-A., Jonsen, P. (2012). Modelling and simulation of explosions in soil interacting with deformable structures, cent.eur.j.eng, 2(4),532-550. [ Links ]

Zakrisson, B., Wikman, B., Häggblad, H.-A. (2011). Numerical simulations of blast loads and structural deformation from near-field explosions in air, International Journal of Impact Engineering, 38(7),597-612. [ Links ]

Zhao, C.F., Chen, J.Y., Wang, Y., Lu, S.J. (2012). Damage mechanism and response of reinforced concrete containment structure under internal blast loading, Theoretical and Applied Fracture Mechanics, 61(0),12-20. [ Links ]

Zukas, J.A., Scheffler, D.R. (2000). Practical aspects of numerical simulations of dynamic events: effects of meshing, International Journal of Impact Engineering, 24(9),925-945. [ Links ]

Zukas, J.A., Walters, W.P. (2002). Explosive Effects and Applications, Springer London, Limited. [ Links ]

* Author e-mail: jovan.trajkovski@fs.uni-lj.si