• Explosives and Blasting
• /
• v.41 no.1
• /
• pp.19-31
• /
• 2023

There are various types of explosives, and each explosive has different characteristics such as water resistance, energy required for detonation, and crushing power, so understanding the characteristics of explosives is important for safe use and performance improvement. Computer simulation is used indirectly along with various experiments to understand the characteristics of explosives, and a state equation is used to express the explosive detonation process through computer simulation. In this study, the explanation of JWL EOS, which is mainly used among the state equations of explosives, and the cylinder expansion experiment to calculate the coefficient of JWL EOS were implemented as ANSYS AUTODYN and compared and analyzed with the actual experimental results. As a result, an error rate of around 20% occurred, and it was found that the overall change pattern of pressure and energy was consistent with the previously published experimental results.

• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.29-35
• /
• 2006

A code is developed to analyze a spherically symmetric underwater explosion. The arbitrary Lagrangian-Eulerian(ALE) Godunov scheme for two-phase flow is used to calculate numerical fluxes through moving control surfaces. For detonation gas of TNT and liquid water, the Jones-Wilkins-Lee(JWL) equation of states and the isentropic Tait relation are used respectively. It is suggested to use the Godunov variable to estimate the velocity of a material interface. The code is validated through comparisons with other results on the gas-water shock tube problem. It is shown that the code can handle generation of discontinuity and recovering of continuity in the normal velocity near the material interface during shock waves interact with the material interface. The developed code is applied to analyze a spherically symmetric underwater explosion. Repeated transmissions of shock waves are clearly captured. The calculated period and maximum radius of detonation gas bubble show good agreements with experimental and other numerical results.

• Geomechanics and Engineering
• /
• v.25 no.2
• /
• pp.111-121
• /
• 2021

Large deformation and rapid pressure propagation take place inside the rock mass under the dynamic loads caused by the explosives, on quarry faces in order to extract aggregate material. The complexity of the science of rock blasting is due to a number of factors that affect the phenomenon. However, blasting engineering computations could be facilitated by innovative software algorithms in order to determine the results of the violent explosion, since field experiments are particularly difficult to be conducted. The present research focuses on the design of a Finite Element Analysis (FEA) code, for investigating in detail the behavior of limestone under the blasting effect of Ammonium Nitrate & Fuel Oil (ANFO). Specifically, the manuscript presents the FEA models and the relevant transient analysis results, simulating the blasting process for three types of limestone, ranging from poor to very good quality. The Finite Element code was developed by applying the Jones-Wilkins-Lee (JWL) equation of state to describe the thermodynamic state of ANFO and the pressure dependent Drucker-Prager failure criterion to define the limestone plasticity behavior, under blasting induced, high rate stress. A progressive damage model was also used in order to define the stiffness degradation and destruction of the material. This paper performs a comparative analysis and quantifies the phenomena regarding pressure, stress distribution and energy balance, for three types of limestone. The ultimate goal of this research is to provide an answer for a number of scientific questions, considering various phenomena taking place during the explosion event, using advanced computational tools.

• Geomechanics and Engineering
• /
• v.23 no.3
• /
• pp.245-259
• /
• 2020

Tunnels have become an indispensable part of metro cities. Blast resistance design of tunnel has attracted the attention of researchers due to numerous implosion event. Present paper deals with the non-linear finite element analysis of rock tunnel having shear zone subjected to internal blast loading. Abaqus Explicit schemes in finite element has been used for the simulation of internal blast event. Structural discontinuity i.e., shear zone has been assumed passing the tunnel cross-section in the vertical direction and consist of Highly Weathered Granite medium surrounding the tunnel. Mohr-Coulomb constitutive material model has been considered for modelling the Highly Weathered Granite and the shear zone material. Concrete Damage Plasticity (CDP), Johnson-Cook (J-C), Jones-Wilkins-Lee (JWL) equation of state models are used for concrete, steel reinforcement and Trinitrotoluene (TNT) simulation respectively. The Coupled-Eulerian-Lagrangian (CEL) method of modelling for TNT explosive and air inside the tunnel has been adopted in this study. The CEL method incorporates the large deformations for which the traditional finite element analysis cannot be used. Shear zone orientations of 0°, 15°, 30°, 45°, 60°, 75° and 90°, with respect to the tunnel axis are considered to see their effect. It has been concluded that 60° orientation of shear zone presents the most critical situation.

• Computers and Concrete
• /
• v.7 no.2
• /
• pp.159-167
• /
• 2010

This paper demonstrates numerical techniques for complex large-scale modeling with microplane constitutive theories for reinforced high strength concrete, which for these applications, is defined to be around the 7000 psi (48 MPa) strength as frequently found in protective structural design. Applications involve highly impulsive loads, such as an explosive detonation or impact-penetration event. These capabilities were implemented into the authors' finite element code, ParaAble and the PRONTO 3D code from Sandia National Laboratories. All materials are explicitly modeled with eight-noded hexahedral elements. The concrete is modeled with a microplane constitutive theory, the reinforcing steel is modeled with the Johnson-Cook model, and the high explosive material is modeled with a JWL equation of state and a programmed burn model. Damage evolution, which can be used for erosion of elements and/or for post-analysis examination of damage, is extracted from the microplane predictions and computed by a modified Holmquist-Johnson-Cook approach that relates damage to levels of inelastic strain increment and pressure. Computation is performed with MPI on parallel processors. Several practical analyses demonstrate that large-scale analyses of this type can be reasonably run on large parallel computing systems.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.6
• /
• pp.521-527
• /
• 2016

This paper is concerned with the numerical analysis of dynamic response of floating offshore wind turbine subject to underwater explosion using an effective non-reflecting technique. An infinite sea water domain was truncated into a finite domain, and the non-reflecting technique called the perfectly matched layer(PML) was applied to the boundary of truncated finite domain to absorb the inherent reflection of out-going impact wave at the boundary. The generalized transport equations that govern the inviscid compressible water flow was split into three PML equations by introducing the direction-wise absorption coefficients and state variables. The fluid-structure interaction problem that is composed of the wind turbine and the sea water flow was solved by the iterative coupled Eulerian FVM and Largangian FEM. And, the explosion-induced hydrodynamic pressure was calculated by JWL(Jones-Wilkins-Lee) equation of state. Through the numerical experiment, the hydrodynamic pressure and the structural dynamic response were investigated. It has been confirmed that the case using PML technique provides more reliable numerical results than the case without using PML technique.