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,Aus...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.展开更多
The Rock-soil interface is a common geological interface.Due to mechanical differences between soil and rock,the stress waves generated by underground blasting undergo intense polarization when crossing the rock-soil ...The Rock-soil interface is a common geological interface.Due to mechanical differences between soil and rock,the stress waves generated by underground blasting undergo intense polarization when crossing the rock-soil interface,making propagation laws difficult to predict.Currently,the characteristics of the impact of the rock-soil interface on blasting stress waves remain unclear.Therefore,the vibration field caused by cylindrical charge blasting in elastic rock and partial-saturation poro-viscoelastic soil was solved.A forward algorithm for the underground blasting vibration field in rock-soil sites was proposed,considering medium damping and geometric diffusion effects of stress waves.Further investigation into the influence of rock and soil parameters and blasting source parameters revealed the following conclusions:stress waves in soil exhibit dispersion,causing peak particle velocity(PPV)to display a discrete distribution.Soil parameters affect PPV attenuation only within the soil,while blasting source parameters affect PPV attenuation throughout the entire site.Multi-wave coupling effects induced by the rocksoil interface result in zones of enhanced and attenuated PPV within the site.The size of the enhancement zone is inversely correlated with the distance from the blasting source and positively correlated with the blasting source attenuation rate and burial depth,providing guidance for selecting explosives and blasting positions.Additionally,PPV attenuation rate increases with distance from the rock-soil interface,but an amplification effect occurs near the interface,most noticeable at 0.1 m.Thus,a sufficient safety distance from the rock-soil interface is necessary during underground blasting.展开更多
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 charact...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.展开更多
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 po...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.展开更多
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 constra...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.展开更多
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 miti...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.展开更多
The accurate prediction of peak overpressure of explosion shockwaves is significant in fields such as explosion hazard assessment and structural protection, where explosion shockwaves serve as typical destructive elem...The accurate prediction of peak overpressure of explosion shockwaves is significant in fields such as explosion hazard assessment and structural protection, where explosion shockwaves serve as typical destructive elements. Aiming at the problem of insufficient accuracy of the existing physical models for predicting the peak overpressure of ground reflected waves, two physics-informed machine learning models are constructed. The results demonstrate that the machine learning models, which incorporate physical information by predicting the deviation between the physical model and actual values and adding a physical loss term in the loss function, can accurately predict both the training and out-oftraining dataset. Compared to existing physical models, the average relative error in the predicted training domain is reduced from 17.459%-48.588% to 2%, and the proportion of average relative error less than 20% increased from 0% to 59.4% to more than 99%. In addition, the relative average error outside the prediction training set range is reduced from 14.496%-29.389% to 5%, and the proportion of relative average error less than 20% increased from 0% to 71.39% to more than 99%. The inclusion of a physical loss term enforcing monotonicity in the loss function effectively improves the extrapolation performance of machine learning. The findings of this study provide valuable reference for explosion hazard assessment and anti-explosion structural design in various fields.展开更多
Accurately forecasting the triple point(TP)path is essential for analyzing blast loads and assessing the destructive effectiveness of the height of burst explosion.Empirical models that describe the TP path under norm...Accurately forecasting the triple point(TP)path is essential for analyzing blast loads and assessing the destructive effectiveness of the height of burst explosion.Empirical models that describe the TP path under normal temperature and pressure environments are commonly employed;however,in certain configurations,such as at high-altitudes(HAs),the environment may involve low temperature and pressure conditions.The present study develops a theoretical prediction model for the TP path under reduced pressure and temperature conditions,utilizing the image bursts method,reflected polar analysis,and dimensional analysis.The model's accuracy is evaluated through numerical simulations and experimental data.Results indicate that the prediction model effectively evaluates the TP path under diminished temperature and pressure conditions,with most predictions falling within a±15%deviation.It was found that the TP height increases with altitude.As the altitude rises from 0 m to 10,000 m,the average TP height increases by 61.7%,87.9%,109.0%,and 134.3%for the scaled height of burst of 1.5 m,2.0 m,2.5 m,and 3.0 m,respectively.Moreover,the variation in TP height under HA environments closely mirrors that observed under corresponding reduced pressure conditions.In HA environments,only the effect of low-pressure conditions on the TP path needs to be considered,as the environmental lowtemperature has a minimal effect.展开更多
Most of the existing studies on tunnel blast wave are based on spherical or grouped charges, however,conventional weapons are mostly cylindrical rather than spherical. In order to analyze the impact of cylindrical cha...Most of the existing studies on tunnel blast wave are based on spherical or grouped charges, however,conventional weapons are mostly cylindrical rather than spherical. In order to analyze the impact of cylindrical charges on the tunnel blast wave loads and to develop a quantitative calculation method, this study carried out experimental and numerical research. Initially, external explosion experiments were conducted using both 35 kg spherical charges and cylindrical charges with aspect ratio of 4.8 at two different distances from the tunnel entrance. Comparative analysis of the blast wave parameters in the tunnel revealed that the explosive equivalent of the cylindrical charges was significantly higher than that of the spherical charges. To address this, an equivalent coefficient κ based on the spherical charges was proposed for the cylindrical charges. Subsequently, numerical simulations were conducted for the experimental conditions, and the numerical simulation results match the experiments well. Through numerical calculations, the reliability of the equivalent coefficient κ under the experimental conditions was verified, and comparison analysis indicated that the explosion energy of cylindrical charges spreads more radially, resulting in more explosion energy entering the tunnel, which is the fundamental reason for the increase in tunnel blast wave loads. Additionally, analyzing the explosion energy ratio entering the tunnel is an effective method for calculating the equivalent coefficient κ. Finally, through more than one hundred sets of numerical calculation results, the impact of the proportional distance λ and the ratio of charge mass to the tunnel cross-section dimension φ on the equivalence coefficients κ was investigated. An empirical formula for the equivalence coefficient κ was derived through fitting, and the accuracy of the formula was validated through literature experimental results. The research findings of this paper will provide valuable guidance for the calculation of blast wave loads in tunnel.展开更多
This study designs four types of liquid-filled cylindrical shell structures to investigate their protection characteristics against explosive shock waves and high-speed fragments.Bare charge and charge-driven prefabri...This study designs four types of liquid-filled cylindrical shell structures to investigate their protection characteristics against explosive shock waves and high-speed fragments.Bare charge and charge-driven prefabricated fragments are employed to examine the damage under blast shock waves and combined blast and fragments loading on various liquid-filled cylindrical shell structures.The test results are compared to numerical calculations and theoretical analysis for the structure's deformation,the liquid medium's movement,and the pressure waves'propagation characteristics under different liquid-filling methods.The results showed that the filling method influences the blast protection and the struc-ture's energy absorption performance.The external filling method reduces the structural deformation,and the internal filling method increases the damage effect.The gapped internal filling method improves the structure's energy absorption efficiency.The pressure wave loading on the liquid-filled cylindrical shell structure differs depending on filling methods.Explosive shock waves and high-speed fragments show a damage enhancement effect on the liquid-filled cylindrical shell structure,depending on the thickness of the internal liquid container layer.The specific impulse on the inner surface of the cylindrical shell positively correlates to the radial deformation of the cylindrical shell structure,and the external liquid layer limits the radial structural deformation.展开更多
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 nu...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.展开更多
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 gr...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.展开更多
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 pro...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.展开更多
The present study focuses on the mitigation of shock wave using novel geometric passages in the flow field.The strategy is to produce multiple shock reflections and diffractions in the passage with minimum flow obstru...The present study focuses on the mitigation of shock wave using novel geometric passages in the flow field.The strategy is to produce multiple shock reflections and diffractions in the passage with minimum flow obstruction,which in turn is expected to reduce the shock wave strength at the target location.In the present study the interaction of a plane shock front(generated from a shock tube)with various geometric designs such as,1)zig-zag geometric passage,2)staggered cylindrical obstructions and 3)zigzag passage with cylindrical obstructions have been investigated using computational technique.It is seen from the numerical simulation that,among the various designs,the maximum shock attenuation is produced by the zig-zag passage with cylindrical obstructions which is then followed by zig-zag passage and staggered cylindrical obstructions.A comprehensive investigation on the shock wave reflection and diffraction phenomena happening in the proposed complex passages have also been carried out.In the new zig-zag design,the initial shock wave undergoes shock wave reflection and diffraction process which swaps alternatively as the shock front moves from one turn to the other turn.This cyclic shock reflection and diffraction process helps in diffusing the shock wave energy with practically no obstruction to the flow field.It is found that by combining the shock attenuation ability of zig-zag passage(using shock reflection and diffraction)with the shock attenuation ability of cylindrical blocks(by flow obstruction),a drastic attenuation in shock strength can be achieved with moderate level of flow blocking.展开更多
This paper presents the explosion cratering effects and their propagation laws of blast waves in dry standard sands using a 450 g-t geotechnical centrifuge apparatus.Ten centrifuge model tests were completed with vari...This paper presents the explosion cratering effects and their propagation laws of blast waves in dry standard sands using a 450 g-t geotechnical centrifuge apparatus.Ten centrifuge model tests were completed with various ranges of explosive mass,burial depth and centrifuge accelerations.Eleven accelerometers were installed to record the acceleration response in sand.The dimensions of the explosion craters were measured after the tests.The results demonstrated that the relationship between the dimensionless parameters of cratering efficiency and gravity scaled yield is a power regression function.Three specific function equations were obtained.The results are in general agreement with those obtained by other studies.A scaling law based on the combination of the π terms was used to fit the results of the ten model tests with a correlation coefficient of 0.931.The relationship can be conveniently used to predict the cratering effects in sand.The results also showed that the peak acceleration is a power increasing function of the acceleration level.An empirical exponent relation between the proportional peak acceleration and distance is proposed.The propagation velocity of blast waves is found to be ranged between 200 and 714 m/s.展开更多
Damage assessment of the wing under blast wave is essential to the vulnerability reduction design of aircraft. This paper introduces a critical relative distance prediction method of aircraft wing damage based on the ...Damage assessment of the wing under blast wave is essential to the vulnerability reduction design of aircraft. This paper introduces a critical relative distance prediction method of aircraft wing damage based on the back-propagation artificial neural network(BP-ANN), which is trained by finite element simulation results. Moreover, the finite element method(FEM) for wing blast damage simulation has been validated by ground explosion tests and further used for damage mode determination and damage characteristics analysis. The analysis results indicate that the wing is more likely to be damaged when the root is struck from vertical directions than others for a small charge. With the increase of TNT equivalent charge, the main damage mode of the wing gradually changes from the local skin tearing to overall structural deformation and the overpressure threshold of wing damage decreases rapidly. Compared to the FEM-based damage assessment, the BP-ANN-based method can predict the wing damage under a random blast wave with an average relative error of 4.78%. The proposed method and conclusions can be used as a reference for damage assessment under blast wave and low-vulnerability design of aircraft structures.展开更多
This article investigates the characteristics of shock wave overpressure generated by multi-layer composite charge under different detonation modes.Combining dimensional analysis and the explosion mechanism of the cha...This article investigates the characteristics of shock wave overpressure generated by multi-layer composite charge under different detonation modes.Combining dimensional analysis and the explosion mechanism of the charge,a peak overpressure prediction model for the composite charge under singlepoint detonation and simultaneous detonation was established.The effects of the charge structure and initiation method on the overpressure field characteristics were investigated in AUTODYN simulation.The accuracy of the prediction model and the reliability of the numerical simulation method were subsequently verified in a series of static explosion experiments.The results reveal that the mass of the inner charge was the key factor determining the peak overpressure of the composite charge under single-point detonation.The peak overpressure in the radial direction improved apparently with an increase in the aspect ratio of the charge.The overpressure curves in the axial direction exhibited a multi-peak phenomenon,and the secondary peak overpressure even exceeded the primary peak at distances of 30D and 40D(where D is the charge diameter).The difference in peak overpressure among azimuth angles of 0-90°gradually decreased with an increase in the propagation distance of the shock wave.The coupled effect of the detonation energy of the inner and outer charge under simultaneous detonation improved the overpressure in both radial and axial directions.The difference in peak overpressure obtained from model prediction and experimental measurements was less than 16.4%.展开更多
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 （...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.展开更多
The paper presents modelling and simulation of a local blast wave interaction with a tire of a logistic truck supporting military operations.In the military industry,it is desired to improve effectiveness and strength...The paper presents modelling and simulation of a local blast wave interaction with a tire of a logistic truck supporting military operations.In the military industry,it is desired to improve effectiveness and strength of vehicle components and simultaneously to minimize the risk of crew injuries.Therefore,the main goal of this paper is to present an attempt to improve blast resistance of a tire.Based on the developed,verified and validated finite element model an optimization procedure was conducted to minimize the damage of a tire subjected to a blast wave.The main issue in the performed computations was to estimate an influence of the cords angle in each layer.For this purpose,a pre-processor script was developed in order to easily modify the finite element model so that the generation process was perfo rmed automatically on the basis of optimization variables.Dynamic response of a tire subjected to blast wave in different cases(cords configurations) was analyzed.It was shown that the optimized cords angles configuration reduces tire local damage and increases its blast resistance.展开更多
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...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.展开更多
基金Foundation item: Project(2010CB226805) supported by the National Basic Research Program of ChinaProject(2010QNA30) supported by the Fundamental Research Funds for the Central Universities of China+1 种基金Project supported by the Priority Academic Development Program of Jiangsu Higher Education,ChinaProjects(SZBF2011-6-B35,2012BAK04B06) supported by the National Twelfth Five-year Key Science & Technology Foundation of China
文摘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.
基金supported by the National Natural Science Foundation of China(Grant Nos.41972286 and 42102329).
文摘The Rock-soil interface is a common geological interface.Due to mechanical differences between soil and rock,the stress waves generated by underground blasting undergo intense polarization when crossing the rock-soil interface,making propagation laws difficult to predict.Currently,the characteristics of the impact of the rock-soil interface on blasting stress waves remain unclear.Therefore,the vibration field caused by cylindrical charge blasting in elastic rock and partial-saturation poro-viscoelastic soil was solved.A forward algorithm for the underground blasting vibration field in rock-soil sites was proposed,considering medium damping and geometric diffusion effects of stress waves.Further investigation into the influence of rock and soil parameters and blasting source parameters revealed the following conclusions:stress waves in soil exhibit dispersion,causing peak particle velocity(PPV)to display a discrete distribution.Soil parameters affect PPV attenuation only within the soil,while blasting source parameters affect PPV attenuation throughout the entire site.Multi-wave coupling effects induced by the rocksoil interface result in zones of enhanced and attenuated PPV within the site.The size of the enhancement zone is inversely correlated with the distance from the blasting source and positively correlated with the blasting source attenuation rate and burial depth,providing guidance for selecting explosives and blasting positions.Additionally,PPV attenuation rate increases with distance from the rock-soil interface,but an amplification effect occurs near the interface,most noticeable at 0.1 m.Thus,a sufficient safety distance from the rock-soil interface is necessary during underground blasting.
基金This researchwas supported by the National Natural Science Foundation of China(No.52227805)the Fundamental Research Funds for Central Universities(No.2022JCCXLJ01).Awards were granted to the author Liyun Yang.
文摘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.
基金supported by the National Nature Science Foundation of China(11672112)the Specialized Research Fund for the Doctoral Program of Higher Education of China(20113718110002)+1 种基金the Fund of the State KeyLaboratory of Disaster Prevention&Mitigation of Explosion&Impact(PLA University and Technology)(DPMEIKF201307)Huaqiao University Research Foundation(13BS402)
文摘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.
基金funded partially by the Australian Government through the Australian Research Council’s Linkage Infrastructure,Equipment and Facilities (LIEF)funding scheme (LE130100133)。
文摘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.
基金supported by the National Natural Science Foundation of China (Grant No.52178515)。
文摘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.
文摘The accurate prediction of peak overpressure of explosion shockwaves is significant in fields such as explosion hazard assessment and structural protection, where explosion shockwaves serve as typical destructive elements. Aiming at the problem of insufficient accuracy of the existing physical models for predicting the peak overpressure of ground reflected waves, two physics-informed machine learning models are constructed. The results demonstrate that the machine learning models, which incorporate physical information by predicting the deviation between the physical model and actual values and adding a physical loss term in the loss function, can accurately predict both the training and out-oftraining dataset. Compared to existing physical models, the average relative error in the predicted training domain is reduced from 17.459%-48.588% to 2%, and the proportion of average relative error less than 20% increased from 0% to 59.4% to more than 99%. In addition, the relative average error outside the prediction training set range is reduced from 14.496%-29.389% to 5%, and the proportion of relative average error less than 20% increased from 0% to 71.39% to more than 99%. The inclusion of a physical loss term enforcing monotonicity in the loss function effectively improves the extrapolation performance of machine learning. The findings of this study provide valuable reference for explosion hazard assessment and anti-explosion structural design in various fields.
基金funding from Anhui Engineering Laboratory of Explosive Materials and Technology Foundation(No.AHBP2022B-04)Natural Science Research Project of Anhui Educational Committee(No.2023AH051221)+1 种基金Anhui Provincial Natural Science Foundation(No.2208085QA26)Scientific Research Foundation for High-level Talents of Anhui University of Science and Technology for the project related to this work.
文摘Accurately forecasting the triple point(TP)path is essential for analyzing blast loads and assessing the destructive effectiveness of the height of burst explosion.Empirical models that describe the TP path under normal temperature and pressure environments are commonly employed;however,in certain configurations,such as at high-altitudes(HAs),the environment may involve low temperature and pressure conditions.The present study develops a theoretical prediction model for the TP path under reduced pressure and temperature conditions,utilizing the image bursts method,reflected polar analysis,and dimensional analysis.The model's accuracy is evaluated through numerical simulations and experimental data.Results indicate that the prediction model effectively evaluates the TP path under diminished temperature and pressure conditions,with most predictions falling within a±15%deviation.It was found that the TP height increases with altitude.As the altitude rises from 0 m to 10,000 m,the average TP height increases by 61.7%,87.9%,109.0%,and 134.3%for the scaled height of burst of 1.5 m,2.0 m,2.5 m,and 3.0 m,respectively.Moreover,the variation in TP height under HA environments closely mirrors that observed under corresponding reduced pressure conditions.In HA environments,only the effect of low-pressure conditions on the TP path needs to be considered,as the environmental lowtemperature has a minimal effect.
文摘Most of the existing studies on tunnel blast wave are based on spherical or grouped charges, however,conventional weapons are mostly cylindrical rather than spherical. In order to analyze the impact of cylindrical charges on the tunnel blast wave loads and to develop a quantitative calculation method, this study carried out experimental and numerical research. Initially, external explosion experiments were conducted using both 35 kg spherical charges and cylindrical charges with aspect ratio of 4.8 at two different distances from the tunnel entrance. Comparative analysis of the blast wave parameters in the tunnel revealed that the explosive equivalent of the cylindrical charges was significantly higher than that of the spherical charges. To address this, an equivalent coefficient κ based on the spherical charges was proposed for the cylindrical charges. Subsequently, numerical simulations were conducted for the experimental conditions, and the numerical simulation results match the experiments well. Through numerical calculations, the reliability of the equivalent coefficient κ under the experimental conditions was verified, and comparison analysis indicated that the explosion energy of cylindrical charges spreads more radially, resulting in more explosion energy entering the tunnel, which is the fundamental reason for the increase in tunnel blast wave loads. Additionally, analyzing the explosion energy ratio entering the tunnel is an effective method for calculating the equivalent coefficient κ. Finally, through more than one hundred sets of numerical calculation results, the impact of the proportional distance λ and the ratio of charge mass to the tunnel cross-section dimension φ on the equivalence coefficients κ was investigated. An empirical formula for the equivalence coefficient κ was derived through fitting, and the accuracy of the formula was validated through literature experimental results. The research findings of this paper will provide valuable guidance for the calculation of blast wave loads in tunnel.
基金the National Natural Science Foundation of China(Grant Nos.52371342,52271338,52101378 and 51979277)。
文摘This study designs four types of liquid-filled cylindrical shell structures to investigate their protection characteristics against explosive shock waves and high-speed fragments.Bare charge and charge-driven prefabricated fragments are employed to examine the damage under blast shock waves and combined blast and fragments loading on various liquid-filled cylindrical shell structures.The test results are compared to numerical calculations and theoretical analysis for the structure's deformation,the liquid medium's movement,and the pressure waves'propagation characteristics under different liquid-filling methods.The results showed that the filling method influences the blast protection and the struc-ture's energy absorption performance.The external filling method reduces the structural deformation,and the internal filling method increases the damage effect.The gapped internal filling method improves the structure's energy absorption efficiency.The pressure wave loading on the liquid-filled cylindrical shell structure differs depending on filling methods.Explosive shock waves and high-speed fragments show a damage enhancement effect on the liquid-filled cylindrical shell structure,depending on the thickness of the internal liquid container layer.The specific impulse on the inner surface of the cylindrical shell positively correlates to the radial deformation of the cylindrical shell structure,and the external liquid layer limits the radial structural deformation.
基金Supported by National Natural Science Foundation of China (No.50638030, 50528808)the National Key Technologies R&D Program of China (No.2006BAJ13B02)the Australian Research Council (No.DP0774061).
文摘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.
文摘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.
基金Supported by National Natural Science Foundation of China(No.50678094)
文摘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.
文摘The present study focuses on the mitigation of shock wave using novel geometric passages in the flow field.The strategy is to produce multiple shock reflections and diffractions in the passage with minimum flow obstruction,which in turn is expected to reduce the shock wave strength at the target location.In the present study the interaction of a plane shock front(generated from a shock tube)with various geometric designs such as,1)zig-zag geometric passage,2)staggered cylindrical obstructions and 3)zigzag passage with cylindrical obstructions have been investigated using computational technique.It is seen from the numerical simulation that,among the various designs,the maximum shock attenuation is produced by the zig-zag passage with cylindrical obstructions which is then followed by zig-zag passage and staggered cylindrical obstructions.A comprehensive investigation on the shock wave reflection and diffraction phenomena happening in the proposed complex passages have also been carried out.In the new zig-zag design,the initial shock wave undergoes shock wave reflection and diffraction process which swaps alternatively as the shock front moves from one turn to the other turn.This cyclic shock reflection and diffraction process helps in diffusing the shock wave energy with practically no obstruction to the flow field.It is found that by combining the shock attenuation ability of zig-zag passage(using shock reflection and diffraction)with the shock attenuation ability of cylindrical blocks(by flow obstruction),a drastic attenuation in shock strength can be achieved with moderate level of flow blocking.
基金supported by the National Natural Science Foundation of China (No. 50779073)the Program Foundation of the China Institute of Water Resources and Hydropower Research (IWHR)(No. YANJI ZD0710)
文摘This paper presents the explosion cratering effects and their propagation laws of blast waves in dry standard sands using a 450 g-t geotechnical centrifuge apparatus.Ten centrifuge model tests were completed with various ranges of explosive mass,burial depth and centrifuge accelerations.Eleven accelerometers were installed to record the acceleration response in sand.The dimensions of the explosion craters were measured after the tests.The results demonstrated that the relationship between the dimensionless parameters of cratering efficiency and gravity scaled yield is a power regression function.Three specific function equations were obtained.The results are in general agreement with those obtained by other studies.A scaling law based on the combination of the π terms was used to fit the results of the ten model tests with a correlation coefficient of 0.931.The relationship can be conveniently used to predict the cratering effects in sand.The results also showed that the peak acceleration is a power increasing function of the acceleration level.An empirical exponent relation between the proportional peak acceleration and distance is proposed.The propagation velocity of blast waves is found to be ranged between 200 and 714 m/s.
基金supported by the Natural Science Foundation of Shaanxi Province (Grant No. 2020JQ-122)the Fund support of Science and Technology on Transient Impact Laboratory。
文摘Damage assessment of the wing under blast wave is essential to the vulnerability reduction design of aircraft. This paper introduces a critical relative distance prediction method of aircraft wing damage based on the back-propagation artificial neural network(BP-ANN), which is trained by finite element simulation results. Moreover, the finite element method(FEM) for wing blast damage simulation has been validated by ground explosion tests and further used for damage mode determination and damage characteristics analysis. The analysis results indicate that the wing is more likely to be damaged when the root is struck from vertical directions than others for a small charge. With the increase of TNT equivalent charge, the main damage mode of the wing gradually changes from the local skin tearing to overall structural deformation and the overpressure threshold of wing damage decreases rapidly. Compared to the FEM-based damage assessment, the BP-ANN-based method can predict the wing damage under a random blast wave with an average relative error of 4.78%. The proposed method and conclusions can be used as a reference for damage assessment under blast wave and low-vulnerability design of aircraft structures.
基金funded by the National Natural Science Foundation of China(Grant No.11972018,No.12002336)China Postdoctoral Science Foundation(Grant No.2021M701710)。
文摘This article investigates the characteristics of shock wave overpressure generated by multi-layer composite charge under different detonation modes.Combining dimensional analysis and the explosion mechanism of the charge,a peak overpressure prediction model for the composite charge under singlepoint detonation and simultaneous detonation was established.The effects of the charge structure and initiation method on the overpressure field characteristics were investigated in AUTODYN simulation.The accuracy of the prediction model and the reliability of the numerical simulation method were subsequently verified in a series of static explosion experiments.The results reveal that the mass of the inner charge was the key factor determining the peak overpressure of the composite charge under single-point detonation.The peak overpressure in the radial direction improved apparently with an increase in the aspect ratio of the charge.The overpressure curves in the axial direction exhibited a multi-peak phenomenon,and the secondary peak overpressure even exceeded the primary peak at distances of 30D and 40D(where D is the charge diameter).The difference in peak overpressure among azimuth angles of 0-90°gradually decreased with an increase in the propagation distance of the shock wave.The coupled effect of the detonation energy of the inner and outer charge under simultaneous detonation improved the overpressure in both radial and axial directions.The difference in peak overpressure obtained from model prediction and experimental measurements was less than 16.4%.
文摘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.
基金The research was carried out under a research grant no.RMN 723the support of the Interdisciplinary Centre for Mathematical and Computational Modelling(ICM)University of Warsaw under grant no GB73-19.This support is gratefully acknowledged.
文摘The paper presents modelling and simulation of a local blast wave interaction with a tire of a logistic truck supporting military operations.In the military industry,it is desired to improve effectiveness and strength of vehicle components and simultaneously to minimize the risk of crew injuries.Therefore,the main goal of this paper is to present an attempt to improve blast resistance of a tire.Based on the developed,verified and validated finite element model an optimization procedure was conducted to minimize the damage of a tire subjected to a blast wave.The main issue in the performed computations was to estimate an influence of the cords angle in each layer.For this purpose,a pre-processor script was developed in order to easily modify the finite element model so that the generation process was perfo rmed automatically on the basis of optimization variables.Dynamic response of a tire subjected to blast wave in different cases(cords configurations) was analyzed.It was shown that the optimized cords angles configuration reduces tire local damage and increases its blast resistance.
基金Supported by the National Defense Basic Science Foundation(B1020060357)
文摘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.
