Effect of Blasting Stress Wave on Dynamic Crack Propagation

Stress waves affect the stress field at the crack tip and dominate the dynamic crack propagation.Therefore,evaluating the influence of blasting stress waves on the crack propagation behavior and the mechanical characteristics of crack propagation is of great significance for engineering blasting.In this study,ANSYS/LS-DYNA was used for blasting numerical simulation,in which the propagation characteristics of blasting stress waves and stress field distribution at the crack tip were closely observed.Moreover,ABAQUS was applied for simulating the crack propagation path and calculating dynamic stress intensity factors(DSIFs).The universal function was calculated by the fractalmethod.The results show that:the compressive wave causes the crack to close and the reflected tensile wave drives the crack to initiate and propagate,and failure mode is mainly tensile failure.The crack propagation velocity varies with time,which increases at first and then decreases,and the crack arrest occurs due to the attenuation of stress waves and dissipation of the blasting energy.In addition,crack arrest toughness is smaller than the crack initiation toughness,applied pressure waveforms(such as the peak pressure,duration,waveforms,wavelengths and loading rates)have a great influence on DSIFs.It is conducive to our deep understanding or the study of blasting stress waves dominated fracture,suggesting a broad reference for the further development of rock blasting in engineering practice.

Blast waveform tailoring using controlled venting in blast simulators and shock tubes

A critical challenge of any blast simulation facility is in producing the widest possible pressure-impulse range for matching against equivalent high-explosive events.Shock tubes and blast simulators are often constrained with the lack of effective ways to control blast wave profiles and as a result have a limited performance range.Some wave shaping techniques employed in some facilities are reviewed but often necessitate extensive geometric modifications,inadvertently cause flow anomalies,and/or are only applicable under very specific configurations.This paper investigates controlled venting as an expedient way for waveforms to be tuned without requiring extensive modifications to the driver or existing geometry and could be widely applied by existing and future blast simulation and shock tube facilities.The use of controlled venting is demonstrated experimentally using the Advanced Blast Simulator(shock tube)at the Australian National Facility of Physical Blast Simulation and via numerical flow simulations with Computational Fluid Dynamics.Controlled venting is determined as an effective method for mitigating the impact of re-reflected waves within the blast simulator.This control method also allows for the adjustment of parameters such as tuning the peak overpressure,the positive phase duration,and modifying the magnitude of the negative phase and the secondary shock of the blast waves.This paper is concluded with an illustration of the potential expanded performance range of the Australian blast simulation facility when controlled venting for blast waveform tailoring as presented in this paper is applied.

Experimental and numerical study on protective effect of RC blast wall against air shock wave

Prototype experiments were carried out on the explosion-proof performance of the RC blast wall.The mass of TNT detonated in the experiments is 5 kg and 20 kg respectively.The shock wave overpressure was tested in different regions.The above experiments were numerically simulated,and the simulated shock wave overpressure waveforms were compared with that tested and given by CONWEP program.The results show that the numerically simulated waveform is slightly different from the test waveform,but similar to CONWEP waveform.Through dimensional analysis and numerical simulation under different working conditions,the equation for the attenuation rate of the diffraction overpressure behind the blast wall was obtained.According to the corresponding standards,the degree of casualties and the damage degree of the brick concrete building at a certain distance behind the wall can be determined when parameters are set.The above results can provide a reference for the design and construction of the reinforced concrete blast wall.

Non-dimensional analysis on blast wave propagation in foam concrete:Minimum thickness to avoid stress enhancement

Foam concrete is a prospective material in defense engineering to protect structures due to its high energy absorption capability resulted from the long plateau stage.However,stress enhancement rather than stress mitigation may happen when foam concrete is used as sacrificial claddings placed in the path of an incoming blast load.To investigate this interesting phenomenon,a one-dimensional difference model for blast wave propagation in foam concrete is firstly proposed and numerically solved by improving the second-order Godunov method.The difference model and numerical algorithm are validated against experimental results including both the stress mitigation and the stress enhancement.The difference model is then used to numerically analyze the blast wave propagation and deformation of material in which the effects of blast loads,stress-strain relation and length of foam concrete are considered.In particular,the concept of minimum thickness of foam concrete to avoid stress enhancement is proposed.Finally,non-dimensional analysis on the minimum thickness is conducted and an empirical formula is proposed by curve-fitting the numerical data,which can provide a reference for the application of foam concrete in defense engineering.

Experimental Research on Amplitude Change of Blasting Seismic Wave with Topography

The propagation characteristics of the amplitude of the blasting seismic wave under the conditions of various topographies are approached by means of experiments. Some factors affecting the effects of quake insulation groove, such as the size, the depth and the position of the quake insulation groove, are studied. The amplitudes of the blasting seismic waves under the conditions of the different sizes of the quake insulation groove are measured. According to the experiments, the effects of the quake insulation groove are related to the position, the distance, the energy of the explosion source and the size of the quake insulation groove itself. The farther it is from the explosion source, the smaller the energy is. The lower its position is and the larger its size is, the more remarkable the effects of the quake insulation groove are.

Numerical study on the fracturing mechanism of shock wave interactions between two adjacent blast holes in deep rock blasting

With the application of electronic detonators, millisecond blasting is regarded as a signifi cant promising approach to improve the rock fragmentation in deep rock blasting. Thus, it is necessary to investigate the fracturing mechanisms of short-delay blasting. In this work, a rectangle model with two circle boreholes is modeled as a particles assembly based on the discrete element method to simulate the shock wave interactions induced by millisecond blasting. The rectangle model has a size of 12 × 6 m (L × W) and two blast holes have the same diameter of 12 cm. The shock waves are simplifi ed as time-varying forces applied at the particles of walls of the two boreholes. Among a series of numerical tests in this study, the spacing between two adjacent boreholes and delay time of millisecond blasting are considered as two primary variables, and the decoupling charge with a coeffi cient of 1.5 is taken into account in each case. The results show that stress superposition is not a key factor for improving rock fragmentation (tensile stress interactions rather than compressive stress superposition could aff ect the generation of cracks), whereas collision actions from isolated particles or particles with weakened constraints play a crucial role in creating the fracture network. The delay time has an infl uence on causing cracks in rock blasting, however, whether it works heavily depends on the distance between the two holes.

Mesh Size Effect in Numerical Simulation of Blast Wave Propagation and Interaction with Structures

Numerical method is popular in analysing the blast wave propagation and interaction with structures.However,because of the extremely short duration of blast wave and energy trans-mission between different grids,the numerical results are sensitive to the finite element mesh size.Previous numerical simulations show that a mesh size acceptable to one blast scenario might not be proper for another case,even though the difference between the two scenarios is very small,indicating a simple numerical mesh size convergence test might not be enough to guarantee accu-rate numerical results.Therefore,both coarse mesh and fine mesh were used in different blast scenarios to investigate the mesh size effect on numerical results of blast wave propagation and interaction with structures.Based on the numerical results and their comparison with field test re-sults and the design charts in TM5-1300,a numerical modification method was proposed to correct the influence of the mesh size on the simulated results.It can be easily used to improve the accu-racy of the numerical results of blast wave propagation and blast loads on structures.

Numerical simulation of the blast wave of a multilayer composite charge

The numerical simulation of a blast wave of a multilayer composite charge is investigated.A calculation model of the near-field explosion and far-field propagation of the shock wave of a composite charge is established using the AUTODYN finite element program.Results of the near-field and far-field calculations of the shock wave respectively converge at cell sizes of 0.25-0.5 cm and 1-3 cm.The Euler--fluxcorrected transport solver is found to be suitable for the far-field calculation after mapping.A numerical simulation is conducted to study the formation,propagation,and interaction of the shock wave of the composite charge for different initiation modes.It is found that the initiation mode obviously affects the shock-wave waveform and pressure distribution of the composite charge.Additionally,it is found that the area of the overpressure distribution is greatest for internal and external simultaneous initiation,and the peak pressure of the shock wave exponentially decays,fitting the calculation formula of the peak overpressure attenuation under different initiation modes,which is obtained and verified by experiment.The difference between numerical and experimental results is less than 10%,and the peak overpressure of both internal and external initiation is 56.12% higher than that of central single-point initiation.

Propagation characteristics of vibration waves induced in surrounding rock by tunneling blasting

The effect of blasting vibration waves on surrounding rock and supporting structures is an important field in underground engineering. In this paper, the separation variable method is used to solve the displacement potential function for the propagation of the blasting vibration waves. In the axis coordinate system, the particle motion and stress change with axial distance, radial distance and time is obtained in surrounding rock. The peak particle velocity law in surrounding rock under different blast loads and surrounding rock parameters is discussed.In addition, the particle vibration characteristics in the surrounding rock are studied using numerical simulations method. The results shows that the peak particle velocity in surrounding rock appears negative exponent attenuation with the increase of axial distance, but it appears positive and negative fluctuations in radial direction. This phenomenon is a new discovery and it has been rarely investigated before. Moreover, the peak particle velocity attenuates more quickly and intensely in the near blasting field,which means that the supporting structure in a shorter distance away from the heading face is vulnerable to the impact of blasting vibration. Theattenuation of blasting vibration velocity is closely related to charge length, blasting load amplitude,attenuation index and rock elastic modulus. The numerical simulation accomplishes the same results and then demonstrates the validity of theoretical results.

Numerical simulation on radiation and energy of blast-induced seismic waves in deep rock masses

With regard to blasting in deep rock masses,it is commonly thought that an increase in the in-situ stress will change the blast-induced rock crack propagation and ultimately affect rock fragmentation.However,little attention has been given to the change in seismic wave radiation when the fractured zone changes with the in-situ stress.In this study,the influences of in-situ stress on blast-induced rock fracture and seismic wave radiation are numerically investigated by a coupled SPH-FEM simulation method.The results show that the change in blast-induced rock fracture with in-situ stress has a considerable effect on the seismic wave energy and composition.As the in-situ stress level increases,the size of the fractured zone is significantly reduced,and more explosion energy is transformed into seismic energy.A reduction in the size of the fractured zone(seismic wave source zone)results in a higher frequency content of the seismic waves.In a nonhydrostatic in-situ stress field,blast-induced cracks are most suppressed in the direction of the minimum in-situ stress,and thus the seismic waves generated in this direction have the highest energy density.In addition to P-waves,Swaves are also generated when a circular explosive is detonated in a nonhydrostatic in-situ stress field.The S-waves result from the asymmetrical release of rock strain energy due to the anisotropic blast-induced fracture pattern.

Seismic effort of blasting wave transmitted in coal-rock mass associated with mining operation

Microseismic effects during the transmission of seismic waves in coal and rock mass associated with mining operation were studied by on-site blasting tests and microseismic monitoring in LW704 of Southern Colliery,Australia,by using spread velocities,amplitudes and frequency contents as the main analysis parameters.The results show that the average P-wave velocity,mean values of combined maximal amplitudes and frequencies of the first arrivals are all reduced significantly along with goaf expanding and intensity weakening of overlying strata during mining process.A full roof fracturing can make the average P-wave velocities,combined maximal amplitudes and frequencies of first arrivals reduce to about 69.8%,92.2% and 60.0%,respectively.The reduction of the above seismic parameters reveals dynamic effects of the variation of strata structure and property to the wave transmission and energy dissipation of blasting wave.The research greatly benefits further study on stability of surrounding rock under the destructive effort by mine tremor,blasting,etc,and provides experimental basis for source relocation and parameter optimization of seismic monitoring as well.

Analysis of Blast Wave Propagation Inside Tunnel

The explosion inside tunnel would generate blast wave which transmits through the longitudinal tunnel. Because of the close-in effects of the tunnel and the reflection by the confining tunnel structure, blast wave propagation inside tunnel is distinguished from that in air. When the explosion happens inside tunnel, the overpressure peak is higher than that of explosion happening in air. The continuance time of the blast wave also becomes longer. With the help of the numerical simulation finite element software LS-DYNA, a three-dimensional nonlinear dynamic simulation analysis for an explosion experiment inside tunnel was carried out. LS-DYNA is a fully integrated analysis program specifically designed for nonlinear dynamics and large strain problems. Compared with the experimental results, the simulation results have made the material parameters of numerical simulation model available. By using the model and the same material parameters, many results were adopted by calculating the model under different TNT explosion dynamites. Then the method of dimensional analysis was used for the simulation results. As overpressures of the explosion blast wave are the governing factor in the tunnel responses, a formula for the explosion blast wave over-pressure at a certain distance from the detonation center point inside the tunnel was derived by using the dimensional analysis theory. By comparing the results computed by the formula with experimental results which were obtained before, the formula was proved to be very applicable at some instance. The research may be helpful to estimate rapidly the effect of internal explosion of tunnel on the structure.

Blasting cumulative damage effects of underground engineering rock mass based on sonic wave measurement

The principle of sonic wave measurement was introduced, and cumulative damage effects of underground engineering rock mass under blasting load were studied by in situ test, using RSM-SY5 intelligent sonic wave apparatus. The blasting test was carried out for ten times at some tunnels of Changba Lead-Zinc Mine. The damage depth of surrounding rock caused by old blasting excavation （0.8-1.2 m） was confirmed. The relation between the cumulative damage degree and blast times was obtained. The results show that the sonic velocity decreases gradually with increasing blast times, hut the damage degree （D） increases. The damage cumulative law is non-linear. The damage degree caused by blast decreases with increasing distance, and damage effects become indistinct. The blasting damage of rock mass is anisotropic. The damage degree of rock mass within charging range is maximal. And the more the charge is, the more severe the damage degree of rock mass is. The test results provide references for researches of mechanical parameters of rock mass and dynamic stability analysis of underground chambers.

Effects of Blast Wave-induced Biomechanical Changes on Lung Injury in Rats

Objective To observe the dynamic impacts of shock waves on the severity of lung injury in rats with different injury distances.Methods Simulate open-field shock waves;detect the biomechanical effects of explosion sources at distances of 40,44,and 48 cm from rats;and examine the changes in the gross anatomy of the lungs,lung wet/dry weight ratio,hemoglobin concentration,blood gas analysis,and pathology.Results Biomechanical parameters such as the overpressure peak and impulse were gradually attenuated with an increase in the injury distance.The lung tissue hemorrhage,edema,oxygenation index,and pathology changed more significantly for the 40 cm group than for the 44 and 48 cm groups.The overpressure peak and impulse were significantly higher for the 40 cm group than for the 44 and 48 cm groups(P<0.05 or P<0.01).The animal mortality was significantly higher for the 40 cm group than for the other two groups(41.2%vs.17.8%and 10.0%,P<0.05).The healing time of injured lung tissues for the 40 cm group was longer than those for the 44 and 48 cm groups.Conclusions The effects of simulated open-field shock waves on the severity of lung injuries in rats were correlated with the injury distances,the peak overpressure,and the overpressure impulse.

Mathematical Model to Locate Interference of Blast Waves from Multi-Hole Blasting Rounds

Maximum charge per delay in a blasting round is universally accepted as the influencing parameter to quantify magni-tude of vibration for any distance of concern. However, for any blasting round experimental data reveals that for same charge per delay magnitude of vibration varies with total charge. Considering linear transmission of blast waves, the paper firstly investigates into the influence of explosive weight, blast design parameters and geology of strata on magnitude and characteristics of vibration parameters and thereafter communicates that possibly interference of blast waves generated from same and different holes of a blasting round result into variation in vibration magnitude. The paper lastly developed a mathematical model to evaluate points of interference of blast waves generated from single- and multi-hole blasting round.

Numerical simulation of the parabolic shell surface under blast wave loading

Deformation of parabolic shell surface under explosion shock waves is a complex dynamic problem. Because of reflection and interference of blast wave, it＇s hard to analytically delineate the dynamic responds of radar parabolic shell surface on blast wave. To gain the characteristics of thin shell deformation under impulsive loading of blast wave, numerical simulation methods for blast load on the shell structure was studied and analyzed. Euler-Lagrange numerical simulation was implemen- ted by AUTODYN code to simulate the problem. Through analysis on deflection feature of radial po- sition under different explosive mass and detonation height, an equation was founded by fitting the deflection results from numerical simulation results of shockwave loading. Experiments were ar- ranged to confirm the validity of the formula. The results gained by simulation are consistent with experiments, and the formula can be used to delineate the deflection of aluminum alloy parabolic shell under blast loading.

Experimental Investigation on Shock Wave Characteristics of Aluminized Explosives in Air Blast

To investigate the shock wave characteristics of RDX-based aluminized explosives,air blast tests were conducted for measuring the parameters of 10 kg aluminized explosives which contained 0-40% aluminum.The results showed that with the increasing of aluminum content,the overpressures and impulses increase at first and then decrease within 7 m or 5 m,which reached the maximum when aluminum content was 20% or 30%.Power exponential formulas are used to fit the shock wave parameters vs scaled distance,where an equal weight of TNT is used to calculate the scaled distance.The overpressures of HL0 and TNT in tested locations not only conform to the similar law,but also conform to the same attenuation law after gaining the scaled distances of equal TNT mass.The pre-exponential factors of overpressure and impulse,kp and kI,decrease along with the increasing of Al content and keep the same pace as the calculated PCJ）.The attenuation coefficients a_P and aIincrease at first and decrease later with the increasing of aluminum content,and they reached the maximal values with30% Al containing,which keeps the same pace as the calculated QV.

Experimental and numerical study of the blast wave decrease using sandwich panel by granular materials core

Among the intrinsic properties of some materials,e.g.,foams,porous materials,and granular materials,are their ability to mitigate shock waves.This paper investigated shock wave mitigation by a sandwich panel with a granular core.Numerical simulations and experimental tests were performed using Autodyn hydro-code software and a shock tube,respectively.The smoothed particle hydrodynamics(SPH)method was used to model granular materials.Sawdust and pumice,whose properties were determined by several compression tests,were used as granular materials in the sandwich panel core.These granular materials possess many mechanisms,including compacting(e.g.,sawdust)and crushing(e.g.,pumice)that mitigate shock/blast wave.The results indicated the ineffectiveness of using a core with low thickness,yet it was demonstrated to be effective with high thickness.Low-thickness pumice yielded better results for wave mitigation.The use of these materials with a core with appropriate core reduces up to 88%of the shock wave.The results of the experiments and numerical simulations were compared,suggesting a good agreement between the two.This indicates the accuracy of simulation and the ability of the SPH method to modeling granular material under shock loading.The effects of grain size and the coefficient of friction between grains have also been investigated using simulation,implying that increasing the grain size and coefficient of friction between grains both reduce overpressure.

DYNAMIC LOAD ANALYSIS OF UNDERGROUND STRUCTURE UNDER EFFECT OF BLAST WAVE

A semi-analytical method of solving the problem of dynamic stress concentration of arbitrary underground structure under the effect of blast waves was introduced. Using the Fourier transform theory, the shock waves （in the forms of SH-waves） can be converted into frequency bands. After employing complex functions and conformal mapping, the admittance functions of various underground structures were obtained. Then, the problem of the time domain dynamic stress response of underground structure can be easily solved through the Fourier inverse transform. At last, the results and curves of the dynamic stress for the square, triangle and horseshoe cavity were presented.

Biomechanical modeling for the response of human thorax to blast waves

A simplified finite element model of a human thorax had been developed for probing into the mechani- cal response in simple and complex blast environments. The human thorax model was first created by CT images with blast loading applied via a coupled arbitrary Lagrangian- Eulerian method, allowing for a variety of loads to be considered. The goal is to analyze the maximum stress distri- butions of lung tissue and peak inward thorax wall velocity and to know the possible regions and levels of lung injury. In parallel, a mathematical model has been modified from the Lobdell model to investigate the detailed percentage of lung injury at each level. The blast loadings around the human tho- rax were obtained from the finite element model, and were then applied in the mathematical model as the boundary con- ditions to predict the normalized work of the human thorax lung. The present results are found in agreement with the modified Bowen curves and the results predicted by Axels- son＇s model.