• Journal of Korean Tunnelling and Underground Space Association
• /
• v.16 no.1
• /
• pp.51-61
• /
• 2014

It is necessary to perform a dynamic analysis to numerically evaluate the effect of blasting on nearby facilities. The blast load time history, which cannot be directly measured, is most often determined from empirical equation. The load has to be adjusted to account for various factors influencing the load and the frequency, but there is not a clear guideline on how to adjust the load. In this study, a series of 2D dynamic numerical analyses that simulates a closely monitored test blasting is performed, from which the blast load that matches the measured vibrations are derived. In the analyses, it is assumed that the hole generated by the blasting is in the form of a circle, and the load was applied normally to the wall of the opening. Special attention was given in selecting the damping ratio for the ground, since it has important influence on the wave propagation and attenuation characteristics of the blast induce waves. The damping ratio was selected such that it matches favorably with the attenuation curve of the measurement. The analyses demonstrate that the empirical blast load widely used in practice highly overstimates the vibration since it does not account for the energy loss due to rock fragmentation. If the empirical load is used without proper adjustment, the numerical analysis may seriously overstimate the predicted vibration, and thus has to be reduced in the analysis.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.5
• /
• pp.287-296
• /
• 2019

In this study, the effect of ground boundary conditions on the evaluation of blast resistance performance of precast arch structures was evaluated by a numerical analysis method. Two types of boundary conditions, namely, fixed boundary conditions and a perfectly matched layer (PML) were applied to numerical models. Blast loads that were much higher than the design load of the target structure were applied to compare the effects of the boundary conditions. The distribution and path of the ground explosion pressure, structural displacement, fracture of concrete, stress of concrete, and reinforcing bars were compared according to the ground boundary condition settings. As a result, the reflecting pressure shock wave at the ground boundaries could be effectively eliminated using PML elements; furthermore, the displacement of the foundation was reduced. However, no distinct difference could be observed in the overall structural behavior including the fracture and stress of the concrete and rebar. Therefore, when blast simulations are performed in the design of protective structures, it is rational to apply the fixed boundary condition on the ground boundaries as conservative design results can be achieved with relatively short computation times.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2000.04a
• /
• pp.63-69
• /
• 2000

Korea Institute of Geology Mining and Materials(KIGAM) cooperating with Southern Methodist University(SMU) has been operating seismo-acoustic array in Chul-Won area to discriminate man-made explosions from natural earthquakes since at the end of July 1999. In order to characterize propagation parameters of detected seismo-acoustic signal and to associate these signals as a blast event accompanying seismic and acoustic signals simultaneously it is necessary to understand infrasound wave propagation in the atmosphere. Two comparable Effective Sound Velocity Structures(ESVS) in atmosphere were constructed by using empirical model (MSISE90 and HWM93) and by aerological observation data of Korea Meteorological Administration (KMA) at O-San area. Infrasound propagation path computed by empirical model resulted in rare arival of refracted waves on ground less than 200km from source region. On the other hand Propagation paths by KMA more realistic data had various arrivals at near source region and well agreement with analyzed seismo-acoustic signals from Chul-Won data. And infrasound propagation in specific direction was very influenced by horizontal wind component in that direction. Linear travel time curve drawn up by 9 days data of the KMA in autumn season showed 335.6m/s apparent sound velocity in near source region. The propagation characteristics will be used to associate seismo-acoustic signals and to calculate propagation parameters of infrasound wave front.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2005.11a
• /
• pp.86-89
• /
• 2005

A computational work has been performed to investigate the aerodynamics of a projectile which is launched from the two-stage light gas gun. A moving coordinate method for a multi-domain technique is employed to simulate unsteady projectile flows with a moving boundary. The effect of a virtual mass is added to the axisymmetric unsteady Euler equation system. The computed results reasonably capture the major flow characteristics which we generated in launching the projectile supersonically, such as the interaction between the shock wave and the blast wave, the interaction between the vortical flow and the barrel shock, and the steady under-expanded jet. The present computational results properly predict the velocity, acceleration, and drag histories of the projectile.

• Tunnel and Underground Space
• /
• v.16 no.6 s.65
• /
• pp.476-485
• /
• 2006

During blasting, both shock wave and gas are generated in detonation process of explosives and the generated wave and gas expansion may create new fractures and damage rock mass. In order to explain and understand completely the fracture mechanism by blasting, we have to consider both effects of the wave and gas expansion simultaneously. In this study, we use a discrete element code, PFC2D and develop an algorithm which is capable of modeling both detonation and gas pressures acting on blasthole wall and visualizing generated cracks within rock mass. Moreover, the gas-pressure modeling method which applies a corresponding external force of gas pressure to parent particles of radial fractures is adopted to simulate a coopting between rock mass and gas penetrating created radial fractures. The developed algorithm is verified by reproducing numerical simulations of a lab-scale test blast successfully.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.10
• /
• pp.649-657
• /
• 2017

A rupture in a high-pressure pipe causes the fluid in the pipe to be discharged in the atmosphere at a high speed resulting in a supersonic jet that generates the compressible flow. This supersonic jet may display complicated and unsteady behavior in general. In this study, Computational Fluid Dynamics (CFD) analysis was performed to investigate the compressible flow generated by a supersonic jet ejected from a high-pressure pipe. A Shear Stress Transport (SST) turbulence model was selected to analyze the unsteady nature of the flow, which depends upon the various gases as well as the diameter of the pipe. In the CFD analysis, the basic boundary conditions were assumed to be as follows: pipe of diameter 10 cm, jet pressure ratio of 5, and an inlet gas temperature of 300 K. During the analysis, the behavior of the shockwave generated by a supersonic jet was observed and it was found that the blast wave was generated indirectly. The pressure wave characteristics of hydrogen gas, which possesses the smallest molecular mass, showed the shortest distance to the safety zone. There were no significant difference observed for nitrogen gas, air, and oxygen gas, which have similar molecular mass. In addition, an increase in the diameter of the pipe resulted in the ejected impact caused by the increased flow rate to become larger and the zone of jet influence to extend further.

• Journal of computational fluids engineering
• /
• v.19 no.4
• /
• pp.68-74
• /
• 2014

This computational study features the transient compressible and inviscid flow analysis on a metallic plasma discharge from the opposing composite electrodes which is subjected to pulsed electric power. The computations have been performed using the flux corrected transport algorithm on the axisymmetric two-dimensional domain of electrode gap and outer space along with the calculation of plasma compositions and thermophysical properties such as plasma electrical conductivity. The mass ablation from aluminum electrode surfaces are modeled with radiative flux from plasma column experiencing intense Joule heating. The computational results shows the highly ionized and highly under-expanded supersonic plasma discharge with strong shock structure of Mach disk and blast wave propagation, which is very similar to muzzle blast or axial plasma jet flows. Also, the geometrical effects of composite electrodes are investigated to compare the amount of mass ablation and penetration depth of plasma discharge.

• Journal of the Korean Geotechnical Society
• /
• v.16 no.6
• /
• pp.59-68
• /
• 2000

현재까지 진동이 양생중인 콘크리트에 미치는 영향을 알아보기 위해 수행된 대부분의 연구에서는 콘크리트 공시체나 콘크리트 블록에 대해 100Hz 미만의 주주파수를 가지는 충격진동이나 진동테이블을 이용한 진동을 가한 후 콘크리트의 강도 변화를 평가하는 방법이 사용되었다. 이 연구에서는 발파진동이 양생중인 라이닝 콘크리트에 미치는 영향을 알아보기 우해 실험실 충격진동 시험과 터널 현장에서의 발파진동 시험을 수행하였다. 터널발파진동과 유사한 100~300Hz의 주주파수를 가지 충격진동을 각각 재령 3, 7, 12 시간에 다한 실험실 시험결과 2cm/sec의 진동속도는 모르타르 라이닝의 P파속도를 증가시키지만, 5 cm/sec, 10cm/sec의진동의 모르타르 라이닝의 P파 속도를 감소시킬 수 있는 것으로 나타났다. 양생기간동안 2.5 cm/sec 이하의 발파진동이 가해진 양생중인 라이닝 콘크리트는 진동을 가하지 않고 양생시킨 콘크리트 공시체에 비해 압축강도가 더 큰 값을 나타내었다. 재령 5시간에 콘크리트 라이닝애 대한 소규모 시험발파로 발파진동을 가한 콘크리트 시료와 진동을 가하지 않고 터널 내에서 양생시킨 공시체에 대해 압축강도를 비교한 결과 콘크리트의 강도 및 탄성파 속도를 저하시킬 수 있는 진동수준은 3~4cm/sec 인 것 으로 나타났다.로 나타났다.

• Explosives and Blasting
• /
• v.23 no.1
• /
• pp.47-54
• /
• 2005

An explosion modeling technique was developed by using the spherical discrete element code, $PFC^{3D}$, which can be used to model the dynamic stress wave propagation phenomenon. The modeling technique is simply based on an idea that the explosion pressure should be applied to a $PFC^{3D}$ particle assembly not in the form of an external force (body force), but in the form of a contact force (surface force). A test blast was conducted for a RC column, whose dimension was $600\times300\times1800$ in millimeters. The initial velocities of the surface movements were measured to be in the range of $14\~18\;m/s$ with the initiation times of $1.5\~2.0m$. Then the blasting procedure was simulated by using the modeling technique. The particle assembly representing the concrete was made of cement mortar and coarse aggregates, whose mirco-properties were obtained from the calibration processes. As a result, the modeling technique developed in this study made it possible for the burden to move with the velocity of $17\~24\;m/s$, which are slightly higher values compared to those of the test blast.

• Journal of the Korean Geotechnical Society
• /
• v.31 no.4
• /
• pp.45-55
• /
• 2015

The stability of the adjcent structures or slopes under blasting is typically evaluated using an empirical vibration attenuation curve or dynamic numerical analysis. To perform a dynamic analysis, it is necessary to determine the blast load and the damping ratio of rock mass. Various empirical methods have been proposed for the blast load. However, a study on representative values of damping ratio of a rock mass has not yet been performed. Therefore, the damping ratio was either ignored or selected without a clear basis in performing a blast analysis. Selection of the dampring ratio for the rock mass is very difficult because the vibration propagation is influenced by the layout and properties of the rock joints. Besides, the vibration induced by blasting is propagated spherically, whereas plane waves are generated by an earthquake. Since the geometrical spreading causes additional attenuation, the damping ratio should be adjusted in the case of a 2D plane strain analysis. In this study, we proposed equivalent damping ratios for use in continuum 2D plane strain analyses. To this end, we performed 2D dynamic analyses for a wide range of rock stiffness and investigated the characteristics of blast vibration propagation. Based on numerical simulations, a correlation between the attenuation equation, shear wave velocity, and equivalent damping ratio of rock mass is presented. This novel approach is the first attempt to select the damping ratio from an attenuation relationship. The proposed chart is easy to be used and can be applied in practice.