• 한국연소학회:학술대회논문집
• /
• 2012.04a
• /
• pp.283-286
• /
• 2012

The analytical model determining the unknown parameters of reaction rate equation which is necessary to simulate the combustion phenomena of energetic materials is proposed. The relationship between detonation velocity and size effect of energetic materials is derived from simplified JWL++ model. Theoretical model is used to investigate the combustion characteristics of certain energetic materials before running Hydrocode by pre-determination of unknown parameter, b. When b=0.8, the behavior of HANFO gunpowder is in the form of concave-up and ANFO explosives has the concave-down form in case of b=1.5. The analytical model provides efficient and highly accurate results rather than previous method which simulated the unconfined-rate-stick via the numerical means.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.462-468
• /
• 2017

KYP model is a pressure-based chemical kinetics that describes shock to detonation transition of energetic materials. In this research, the parameters of KYP model and JWL EOS for PETN-based explosive, namely PBXN-301, were determined. A series of unconfined rate stick tests and two dimensional hydrodynamic simulation were conducted to obtain the size effect behaviour of the explosive. As a result, it was confirmed that the parameters obtained from KYP modeling have more accuracy to predict the detonation velocities according to the inverse radius of PBXN-301 than the qualitatively obtained LLNL constitutive equations.

• Geotechnical Engineering
• /
• v.13 no.1
• /
• pp.5-18
• /
• 1997

In order to make reliable ground shock predictions in saturated geological media, it is necessary to use multi -phase material models and numerical codes. This paper presents the results of theoretical study of the fundamental behavior of multi-phase porous media subjected to high dynanlic loadings, and deals with the development of numerical code MPDAP with JWL(Jones-Wilkins-Lee) model, which is capable of considering the kinds and characters of explosives. To check the global equilhorium equations of the numerical code, we carried out some verifications. In the cases of the elastic spherical wave propagation in a single phase medium, one-dimensional linear ronsolidation, and one timensional wave propagation in saturated linear elastic soils and rocks, the results calculated by MPDAP show close agreement with closed-form solutions or numerical solutions generated with two phase code.

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

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

• Advances in Computational Design
• /
• v.5 no.4
• /
• pp.417-443
• /
• 2020

Tunnels have been an integral part of human civilization. Due to complexity in its design and structure, the stability of underground structures under extreme loading conditions has utmost importance. Increased terrorism and geo-political conflicts have forced the engineers and researchers to study the response of underground structures, especially tunnels under blast loading. The present study has been carried out to seek the response of tunnel structures under blast load using the finite element technique. The tunnel has been considered in quartzite rock of northern India. The Mohr-Coulomb constitutive model has been adopted for the elastoplastic behaviour of rock. The rock model surrounding the tunnel has dimensions of 30 m x 30 m x 35 m. Both unlined and lined (concrete) tunnel has been studied. Concrete Damage Plasticity model has been considered for the concrete lining. Four different parameters (i.e., tunnel diameter, liners thickness, overburden depth and mass of explosive) have been varied to observe the behaviour under different condition. To carry out blast analysis, Coupled-Eulerian-Lagrangian (CEL) modelling has been adopted for modelling of TNT (Trinitrotoluene) and enclosed air. JWL (Jones-Wilkins-Lee) model has been considered for TNT explosive modelling. The paper concludes that deformations in lined tunnels follow a logarithmic pattern while in unlined tunnels an exponential pattern has been observed. The stability of the tunnel has increased with an increase in overburden depth in both lined and unlined tunnels. Furthermore, the tunnel lining thickness also has a significant effect on the stability of the tunnel, but in smaller diameter tunnel, the increase in tunnel lining thickness has not much significance. The deformations in the rock tunnel have been decreased with an increase in the diameter of the tunnel.

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