# Numerical Study on the Relationship between the Dominant Frequency of Blasting Vibration and the Development of Plastic Zone in Cylindrical Charge

Authors: Da Liu; Xiaohua Zhao; Jianglin Gao; Binghan Xue; Xiaobin Wu; Songtao Hu; Fang Chen
Journal: Geofluids
 (2022)

---

## Abstract

Blasting vibration is a widely studied harmful effect of rock blasting excavations. Many factors affect its dominant frequency, which makes analyzing, evaluating, and predicting it difficult. This study explored the factors influencing the dominant frequency of blasting vibrations in the case of a spherical charge. Based on symmetry, a theoretical analysis in terms of a spherical explosion source is generally sufficient to describe a cavity excited by a spherical explosion charge. The elastic cavity radius and the dominant frequency of vibration induced by the spherical blasting source are closely related. However, there is a lack of relevant research on cylindrical charges. Therefore, a calculation model for a single-hole cylindrical charge was established. There is a relationship between the corresponding dominant frequency and the range of the plastic zone. The results indicate that the dominant frequency of the blasting vibration for both cylindrical and spherical charges is closely related to the range of the plastic zone formed by the rock blasting. As the elastic cavity radius increases, both the zero-crossing and the dominant Fourier frequencies decrease, and the amplitude spectrum shifts to lower frequencies. However, increasing the cylindrical charge diameter causes more changes in the plastic zone in the direction perpendicular to the cylindrical explosive axis. Moreover, increasing the charge length causes more changes in the plastic zone along this axis. It is therefore difficult to identify a unique dimensional parameter that characterizes the range of the plastic zone formed by the blasting. Because the plastic zone around a cylindrical charge is less regular than around a spherical charge, the charge weightQ is a more favorable parameter than the elastic cavity radius a when used as the main influencing factor and in an attenuation analysis for the dominant frequency of cylindrical charge blasting.


---

## Body

## 1. Introduction

As an efficient excavation method, borehole blasting is easy to induce harmful vibration effects, which has always been the focus of attention [1–3]. The dominant frequency induced by cylindrical charge is an important index to characterize blasting vibration [4, 5], but the theoretical research is not sufficient compared with spherical charge. Therefore, to fully understand the vibration characteristics of cylindrical charge qualitatively and quantitatively, it is often necessary to learn from the relevant spherical charge theory [6]. Previous studies have demonstrated that there are two stages of blasting vibration effect induced by spherical charge. The first stage is the initial state of explosive expansion, where the energy released by the explosive expands symmetrically in all directions, and P wave is generated, causing compression deformation in the surrounding rock mass. The second stage is the roof-like uplift caused by the expansion of soil within a certain epicenter distance due to the residual pressure of explosive gas in the spherical cavity, which is considered as the elastic vibration in the near region. For the spherical explosion source, due to the spherical symmetry of the explosion load, the problem can be simplified to a one-dimensional problem, and the solution of the problem can be obtained through analytical calculation [7] so that the relationship between the vibration spectrum and the radius of the plastic zone can be analyzed and obtained. However, for cylindrical charge, the analytical solution of its dominant frequency is difficult to obtain due to the influence of many factors such as charge structure, explosion energy release form, and free surface [8, 9].According to the degree of damage to the rock mass around the borehole blasting, it can be simply divided into crushed zone, cracked zone, and elastic formation zone [10]. In fact, the blasting damage of rock mass is a complex dynamic accumulation process. When discriminating the range of rock damage, most of the current methods to evaluate rock blasting damage are based on numerical modeling [11]. Grady [12] proposed an isotropic damage model for rock blasting, which uses a scalar to describe the deterioration of rock stiffness and assumes that the number of cracks obeys the two parameter Weibull distribution. Yang et al. [13] believe that the crack can propagate only when the bulk strain is greater than a certain value and proposed the concept of fracture probability, which is introduced in the definition of damage variable. Ma [14] introduced the Johnson-Holmquist (J-h) material model into LS-DYNA to simulate the rock fracture damage process caused by blasting.To obtain the range of rock blasting damage area, most methods are trying to calculate the peak particle velocity (PPV) produced by the detonating charge compared against a PPV that is known or adopted (most of the time based on site-specific field tests) that will produce some damage in the rock or rock mass [15]. Holmberg-Perssion [16] assumes that the cylindrical charge is regarded as a small section connected in series, and the PPV distribution in the rock mass can be obtained by combining each small section of charge. Based on the cylindrical wave theory and the wavelet theory of spherical wave and long cylindrical charge, Hustrulid and Lu [17] deduced the theoretical formula for the attenuation of peak vibration velocity of blasting vibration particles by considering the detonation velocity of explosives from the viewpoint proposed by Starfield [18]. Most of the studies combine the damage range of rock with PPV, and generally, the larger the PPV induced by rock blasting, the larger the range of cracked zone caused by blasting. In fact, similar to PPV, the dominant frequency of blasting vibration is also an important parameter to characterize blasting vibration [19], which has two most commonly used definitions, including zero-crossing frequency fz (the frequency corresponding to the peak vibration period) and Fourier frequency fF (the frequency corresponding to the maximum amplitude in the spectrum), shown in Figure 1. And the dominant frequency is also closely related to the damage range formed by rock blasting [20]. However, there is much less research on the relationship between the dominant frequency of blasting vibration and the range of cracked zones caused by blasting. Therefore, through the research, it can provide some meaningful perspectives for exploring the propagation characteristics and attenuation characteristics of the dominant frequency of blasting vibration.Figure 1
Typical waveform and spectrum of blasting vibration.
(a)
Velocity–time histories of blasting(b)
Amplitude spectrumOn the basis of the stress wave excited by spherical cavity explosion in an elastic rock mass, the theoretical solution of spherical wave propagation and the velocity spectrum are analyzed. By establishing the dynamic finite element model of single borehole cylindrical charge, the distribution characteristics of rock blasting cracked zone are analyzed, and the numerical relationship between the cracked zone of blasting and the dominant frequency is explored, which is helpful to understand the factors affecting the dominant frequency of blasting vibration by cylindrical charge.


## 2. Relationship between the Cracked Zone and the Dominant Frequency

The explosive blasting causes a large transient impact load acting on the rock, and the energy released by the explosive propagates outward in the form of stress wave. The regions encountered by a blasting stress wave propagating away from the blasting center can be approximately classified as the crushed, cracked, and elastic areas, according to the degrees of damage in the rock, and the propagation properties of shock, stress, and seismic waves vary accordingly. In the crushed zone, the rock is crushed because the pressure caused by the blasting load greatly exceeds its compressive strength. In the cracked (fracture) zone, the circumferential tensile stress of the rock mass often exceeds its tensile strength; the rock mass also undergoes an irreversible tensile deformation. In the elastic deformation zone, the load on the rock mass is only sufficient to cause elastic vibration and not plastic deformation. The blasting damage zone of the rock mass is illustrated in Figure2.Figure 2
Damage zones surrounding a spherical detonation.To facilitate the analysis, the nonelastic zones (crushed zone and cracked (fracture) zone) near the blasthole are generally regarded as the equivalent blasting source, in which the elastic-plastic boundaries are applied on the blasting load [21–23]. Therefore, the theoretical solution of an elastic wave excited by a spherical cavity in an elastic medium can be adopted [7].Introducing a potential functionφ to represent the radial displacement
(1)u=∂φr,t∂r.Assuming the load function acting on the inner wall of the spherical cavity ispt,
(2)φr,t=−aρωr∫0sps−te−ητsinωτdτ,(3)s=t−r−aCP,(4)η=1−2ν1−νCPa,(5)ω=1−2ν1−νCPa,where pt is the blasting load acting on the elastic cavity, λ and μ are the constants of lame, ν is Poisson’s ratio, ρ is the density, and Cp is the longitudinal wave velocity. a is the radius of the elastic cavity (cracked zone). ω is the phase of the function φ, which is also the natural frequency of the equivalent system.When the dynamic load of stress waves over time is adopted in the form(6)pt=0t<−τ1,σmax1+t/τ1−τ1≤t≤0,σmax1−t/τ20≤t≤τ2,0t>τ2,where σmax is the peak value of blasting load, τ1 is the rising time of blasting load, and τ2 is the time of blasting load reduction from peak to zero.By means of Fourier transform, the velocity spectrum of blasting vibration in elastic rock mass due to the action ofpt on the elastic boundary can be obtained as follows:
(7)Fω=SσjωaωCpCp2+r2ω24μr2Cp/a4+1−λ+2μ/2μω2Cp/a2+λ+2μ/4μ2ω4,Sσjω=σmaxaebeτω21+ae2+be2+2aebecosωτ−2aecosbeωτ+becosaeωτ1/2,where τ=τ1+τ2, ae=τ1/τ, and be=τ2/τ.The same rock mass parameters are selected: densityρ=2700kg/m3, Poisson’s ratio=0.22, peak value of blasting load σ=50MPa, duration of blasting load τ=10ms, and rising time of blasting load τ1=2ms. The influence of elastic cavity radius on the blasting vibration spectrum is analyzed by substituting the rock parameters into the formula. The velocity amplitude spectra corresponding to different elastic cavity radius (elastic cavity radius a=10m, 20 m, 30 m, 40 m, 50 m) are plotted in Figure 3(a), and the normalized velocity amplitude spectra are plotted in Figure 3(b).Figure 3
Amplitude spectrum variation law of different elastic cavity radius.
(a)
Velocity amplitude spectrum(b)
Normalized to the peak value of the spectrumAs shown in Figure3, increasing the radius of the elastic cavity increases the amplitude of the vibration velocity spectrum and shifts the peak toward lower frequencies. The dominant Fourier frequency decreases with increasing radius of the elastic cavity.


## 3. Relationship between the Plastic Zone and the Dominant Frequency in Cylindrical Charge

### 3.1. Range of Plastic Zone

Under the condition of cylindrical charge, when the external load of the rock around the explosive reaches a certain value, the stress exceeds the yield strength of the rock, which brings irreversible plastic deformation. Cylindrical charge is different from the spherical charge that the plastic zone is not an approximate spherical zone, and the evolution radius of the plastic zone is not easy to define. In this study, the plastic deformation range of rock around the borehole is compared under different working conditions, and the plastic radius of the middle section of charge length with the bottom initiating cylindrical charge (hereinafter referred to as “the radius of the middle section”) is taken for comparison under all different working conditions. In LS-DYNA program, when the equivalent stress is equal to the yield stress of the material, it is also used as the judgment standard for the material to enter plastic deformation. The yield criterion at coordinatesi,j is as follows:
(8)12sijsij−13σs2=0,where sij=σij−σmδij is the deviatoric stress tensor, σm=1/3σ11+σ23+σ33 is the average stress, σs is the yield stress of rock, when the plastic deformation of rock occurs, and σ11,σ22,σ33 is the first, second, and third principal stresses, respectively.
(9)σ¯=12σ1−σ22+σ2−σ32+σ3−σ12.Under the action of explosion load, the equivalent stress of rockσ¯ is greater than the yield stress of rock σs.


### 3.2. Numerical Model and Calculation Conditions

The blasting process of a cylindrical charge in an infinite rock mass was simulated numerically using a finite-element model and numerically analyzed using the LS-DYNA software. The model represents a quarter cylinder of radius 80 m, as shown in Figure4.Figure 4
Quarter cylinder finite-element model.The dynamic impact of explosives is simulated by using the fluid-solid coupling algorithm, and the relationship between pressure and volume during the explosion is calculated by combining with JWL equation of state, in which the calculation parameters of the emulsion explosives used are shown in Table1.Table 1
Calculation parameters of JWL equation of state for explosives.
MaterialDensity (g/m3)V (m/s)CJ pressure (Pa)A (GPa)B (GPa)R1R2ωExplosion1×10336003.24×1092140.184.20.90.15V is the velocity of detonation, A and B are the product JWL coefficient, R1 and R2 are the unreacted JWL coefficient, and CJ pressure is the initial pressure.The rock mass is simplified as an isotropic ideal elastic-plastic material simulated by the bilinear elastic-plastic model in Table2.Table 2
Calculation parameters of rock materials.
MaterialElastic modulus (GPa)Density (g/m3)Poisson’s ratioShear modulus (GPa)Bulk modulus (GPa)Yield stress (MPa)Rock2024000.248.0612.840The development of the elastic cavity radius was analyzed by changing the radius and length of the charge, which are the dimensions commonly used in the context of foundation excavation in hydraulic engineering. The variations of the elastic cavity radius for different charge parameters are listed in Tables3 and 4.Table 3
The radius of the middle section under different charge diameters.
Charge diameterd (mm)Charge lengthl (m)Charge weightQ (kg)The radius of the middle section (m)90850.891.96110876.032.591308106.193.19Table 4
Radius of plastic zone under different charge diameters.
Charge diameterd (mm)Charge lengthl (m)Charge weightQ (kg)The radius of the middle section (m)110657.022.59110876.032.591101095.032.59Under the condition of a constant charge length, we calculated the radii of the middle section formed by blasting, for charge diametersd=90, 110, and 130 mm.Under the condition of a constant charge diameter, we calculated the radii of the middle section formed by blasting, for charge lengthl=6m, 8 m, and 10 m.As shown in Table3, the plastic radius of the middle section of the charge length increases with the charge diameter between 90 and 130 mm. According to Table 4, the range of the plastic zone increases when the charge diameter remains constant at 110 mm while the charge length increases from 6 to 10 m. However, its size increases mainly along the axial direction of the cylindrical charge, and the plastic radius of the middle section for charge lengths a1,a2,a3 shows no significant change.Figures5 and 6 suggest that an increase in the charge diameter and the charge length increases the range of the plastic zone. An increase in the charge diameter tends to increase the range of the plastic zone in the direction perpendicular to the explosive axis, and an increase in the charge length tends to increase the range of the plastic zone along the explosive axis direction. Therefore, the diameter of the plastic range of the cylindrical charge is more difficult to define owing to the shape of the charge in contrast to a spherical charge.Figure 5
Elastic-plastic boundary under different charge diameter.Figure 6
Elastic-plastic boundary under different charge lengths.


### 3.3. Analysis Results under Different Charge Diameters

Measuring points at the elastic-plastic boundary in the middle section were selected for analysis, in the case of different charge diametersd=90, 110, and 130 mm. The waveform and frequency spectrum at these measuring points are compared and analyzed in Figure 7 and Table 5.Figure 7
Waveforms and spectra of the elastic-plastic boundary under different charge diameters.
(a)
Velocity–time histories(b)
Amplitude spectrumTable 5
Radius of plastic zone under different charge diameters.
Charge diameterd (mm)Charge weightQ (kg)The plastic range radius of the middle sectiona (m)Zero-crossing dominant frequency (Hz)Fourier dominant frequency (Hz)9050.891.96322.6371.811076.032.59270.3332.7130106.193.19263.2313.1The plastic radius of the middle section shows little change for charge lengthsl=6, 8, and 10 m. Measuring points on the elastic-plastic boundary of the middle section were selected for analysis, and the corresponding waveforms and spectra are compared and analyzed in Figure 8 and Table 6 for various working conditions.Figure 8
Waveforms and spectra of the elastic-plastic boundary under different charge diameters.
(a)
Velocity–time histories(b)
Amplitude spectrumTable 6
Radius of plastic zone under different charge length.
Charge lengthl (m)Charge weightQ (kg)The plastic range radius of the middle sectiona (m)Zero-crossing dominant frequency (Hz)Fourier dominant frequency (Hz)657.022.59294.1371.8876.032.59270.3332.71095.032.59263.2313.1When increasing either the charge diameter or length, the plastic zone range also increases. However, changing the charge diameter and length causes the plastic zone to change in two different directions relative to the explosive geometry. This is a consequence of the more complex charge structure and load energy-release mechanism pertaining to a cylindrical charge, compared to a spherical charge. Therefore, the diameter of the plastic zone induced by the cylindrical charge is more difficult to define than in the case of a spherical charge.Figures7 and 8 suggest that increasing the charge diameter and length causes an increase in the range of the plastic zone and a decrease in the zero-crossing dominant frequency and in the dominant Fourier frequency. This is generally consistent with the results obtained for a spherical charge. However, the charge diameter and the charge length cause the plastic zone to change in two different directions relative to the explosive geometry. Specifically, an increase in the charge diameter causes an increase in the plastic range in the direction perpendicular to the explosive axis, and an increase in the charge length causes an increase in the plastic range along the explosive axis.According to Equation (5), the dominant frequency of a blasting vibration often has a proportional relationship with Cp/a or Cp/Q. Most existing formulas for this frequency show Cp/a or Cp/Q as the proportionality factor [19, 24, 25]. Under the condition of spherical charge, the charge weight Q and the radius of elastic cavity a are interchangeable according to the relationship expressed as follows:
(10)Q=43πa3q,whereq is the unit explosive consumption of rock.The plastic zone radiusa of a cylindrical charge is challenging to determine in an actual blasting process. (For example, as shown in Table 6, for a constant charge diameter and a charge length varying from 6 to 10 m, the plastic radius of the middle section does not change significantly.) It is therefore often replaced by the charge weight Q. A reasonable correlation is observed between the radius of the plastic zone and the dominant frequency of the cylindrical charge in Figures 9 and 10. However, owing to the complexity of the definition of plastic zone, the charge weight Q is a more favorable parameter than the elastic cavity radius a when used as the main influencing factor and in the attenuation analysis of the dominant frequency of cylindrical charge blasting.Figure 9
Linear fitting relationship between the dominant frequency andCp/a.
(a)
Zero-crossing frequency(b)
Fourier frequencyFigure 10
Linear fitting relationship between the dominant frequency andCp/Q.
(a)
Zero-crossing frequency(b)
Fourier frequency


## 3.1. Range of Plastic Zone

Under the condition of cylindrical charge, when the external load of the rock around the explosive reaches a certain value, the stress exceeds the yield strength of the rock, which brings irreversible plastic deformation. Cylindrical charge is different from the spherical charge that the plastic zone is not an approximate spherical zone, and the evolution radius of the plastic zone is not easy to define. In this study, the plastic deformation range of rock around the borehole is compared under different working conditions, and the plastic radius of the middle section of charge length with the bottom initiating cylindrical charge (hereinafter referred to as “the radius of the middle section”) is taken for comparison under all different working conditions. In LS-DYNA program, when the equivalent stress is equal to the yield stress of the material, it is also used as the judgment standard for the material to enter plastic deformation. The yield criterion at coordinatesi,j is as follows:
(8)12sijsij−13σs2=0,where sij=σij−σmδij is the deviatoric stress tensor, σm=1/3σ11+σ23+σ33 is the average stress, σs is the yield stress of rock, when the plastic deformation of rock occurs, and σ11,σ22,σ33 is the first, second, and third principal stresses, respectively.
(9)σ¯=12σ1−σ22+σ2−σ32+σ3−σ12.Under the action of explosion load, the equivalent stress of rockσ¯ is greater than the yield stress of rock σs.


## 3.2. Numerical Model and Calculation Conditions

The blasting process of a cylindrical charge in an infinite rock mass was simulated numerically using a finite-element model and numerically analyzed using the LS-DYNA software. The model represents a quarter cylinder of radius 80 m, as shown in Figure4.Figure 4
Quarter cylinder finite-element model.The dynamic impact of explosives is simulated by using the fluid-solid coupling algorithm, and the relationship between pressure and volume during the explosion is calculated by combining with JWL equation of state, in which the calculation parameters of the emulsion explosives used are shown in Table1.Table 1
Calculation parameters of JWL equation of state for explosives.
MaterialDensity (g/m3)V (m/s)CJ pressure (Pa)A (GPa)B (GPa)R1R2ωExplosion1×10336003.24×1092140.184.20.90.15V is the velocity of detonation, A and B are the product JWL coefficient, R1 and R2 are the unreacted JWL coefficient, and CJ pressure is the initial pressure.The rock mass is simplified as an isotropic ideal elastic-plastic material simulated by the bilinear elastic-plastic model in Table2.Table 2
Calculation parameters of rock materials.
MaterialElastic modulus (GPa)Density (g/m3)Poisson’s ratioShear modulus (GPa)Bulk modulus (GPa)Yield stress (MPa)Rock2024000.248.0612.840The development of the elastic cavity radius was analyzed by changing the radius and length of the charge, which are the dimensions commonly used in the context of foundation excavation in hydraulic engineering. The variations of the elastic cavity radius for different charge parameters are listed in Tables3 and 4.Table 3
The radius of the middle section under different charge diameters.
Charge diameterd (mm)Charge lengthl (m)Charge weightQ (kg)The radius of the middle section (m)90850.891.96110876.032.591308106.193.19Table 4
Radius of plastic zone under different charge diameters.
Charge diameterd (mm)Charge lengthl (m)Charge weightQ (kg)The radius of the middle section (m)110657.022.59110876.032.591101095.032.59Under the condition of a constant charge length, we calculated the radii of the middle section formed by blasting, for charge diametersd=90, 110, and 130 mm.Under the condition of a constant charge diameter, we calculated the radii of the middle section formed by blasting, for charge lengthl=6m, 8 m, and 10 m.As shown in Table3, the plastic radius of the middle section of the charge length increases with the charge diameter between 90 and 130 mm. According to Table 4, the range of the plastic zone increases when the charge diameter remains constant at 110 mm while the charge length increases from 6 to 10 m. However, its size increases mainly along the axial direction of the cylindrical charge, and the plastic radius of the middle section for charge lengths a1,a2,a3 shows no significant change.Figures5 and 6 suggest that an increase in the charge diameter and the charge length increases the range of the plastic zone. An increase in the charge diameter tends to increase the range of the plastic zone in the direction perpendicular to the explosive axis, and an increase in the charge length tends to increase the range of the plastic zone along the explosive axis direction. Therefore, the diameter of the plastic range of the cylindrical charge is more difficult to define owing to the shape of the charge in contrast to a spherical charge.Figure 5
Elastic-plastic boundary under different charge diameter.Figure 6
Elastic-plastic boundary under different charge lengths.


## 3.3. Analysis Results under Different Charge Diameters

Measuring points at the elastic-plastic boundary in the middle section were selected for analysis, in the case of different charge diametersd=90, 110, and 130 mm. The waveform and frequency spectrum at these measuring points are compared and analyzed in Figure 7 and Table 5.Figure 7
Waveforms and spectra of the elastic-plastic boundary under different charge diameters.
(a)
Velocity–time histories(b)
Amplitude spectrumTable 5
Radius of plastic zone under different charge diameters.
Charge diameterd (mm)Charge weightQ (kg)The plastic range radius of the middle sectiona (m)Zero-crossing dominant frequency (Hz)Fourier dominant frequency (Hz)9050.891.96322.6371.811076.032.59270.3332.7130106.193.19263.2313.1The plastic radius of the middle section shows little change for charge lengthsl=6, 8, and 10 m. Measuring points on the elastic-plastic boundary of the middle section were selected for analysis, and the corresponding waveforms and spectra are compared and analyzed in Figure 8 and Table 6 for various working conditions.Figure 8
Waveforms and spectra of the elastic-plastic boundary under different charge diameters.
(a)
Velocity–time histories(b)
Amplitude spectrumTable 6
Radius of plastic zone under different charge length.
Charge lengthl (m)Charge weightQ (kg)The plastic range radius of the middle sectiona (m)Zero-crossing dominant frequency (Hz)Fourier dominant frequency (Hz)657.022.59294.1371.8876.032.59270.3332.71095.032.59263.2313.1When increasing either the charge diameter or length, the plastic zone range also increases. However, changing the charge diameter and length causes the plastic zone to change in two different directions relative to the explosive geometry. This is a consequence of the more complex charge structure and load energy-release mechanism pertaining to a cylindrical charge, compared to a spherical charge. Therefore, the diameter of the plastic zone induced by the cylindrical charge is more difficult to define than in the case of a spherical charge.Figures7 and 8 suggest that increasing the charge diameter and length causes an increase in the range of the plastic zone and a decrease in the zero-crossing dominant frequency and in the dominant Fourier frequency. This is generally consistent with the results obtained for a spherical charge. However, the charge diameter and the charge length cause the plastic zone to change in two different directions relative to the explosive geometry. Specifically, an increase in the charge diameter causes an increase in the plastic range in the direction perpendicular to the explosive axis, and an increase in the charge length causes an increase in the plastic range along the explosive axis.According to Equation (5), the dominant frequency of a blasting vibration often has a proportional relationship with Cp/a or Cp/Q. Most existing formulas for this frequency show Cp/a or Cp/Q as the proportionality factor [19, 24, 25]. Under the condition of spherical charge, the charge weight Q and the radius of elastic cavity a are interchangeable according to the relationship expressed as follows:
(10)Q=43πa3q,whereq is the unit explosive consumption of rock.The plastic zone radiusa of a cylindrical charge is challenging to determine in an actual blasting process. (For example, as shown in Table 6, for a constant charge diameter and a charge length varying from 6 to 10 m, the plastic radius of the middle section does not change significantly.) It is therefore often replaced by the charge weight Q. A reasonable correlation is observed between the radius of the plastic zone and the dominant frequency of the cylindrical charge in Figures 9 and 10. However, owing to the complexity of the definition of plastic zone, the charge weight Q is a more favorable parameter than the elastic cavity radius a when used as the main influencing factor and in the attenuation analysis of the dominant frequency of cylindrical charge blasting.Figure 9
Linear fitting relationship between the dominant frequency andCp/a.
(a)
Zero-crossing frequency(b)
Fourier frequencyFigure 10
Linear fitting relationship between the dominant frequency andCp/Q.
(a)
Zero-crossing frequency(b)
Fourier frequency


## 4. Discussion

Explosive blasting in rock forms a rock damage area. For cylindrical charge, there have been many studies on the definition and range of rock damage, which are different according to the rock damage criteria. For instance, Xayyraed [26] shows that when blasting a single hole in a semi-infinite rock medium with a known blasting load, the crushed zone radius r1 and fractured (cracked) zone radius r2 take the form
(11)r1=ρrCP25σc1/2Pσ1/4a,(12)r2=μP1−μσt1/βa,where μ is the Poisson ratio, Cp is the P wave velocity of the rock mass, σc and σt are the unconfined compressive strength and tensile strength of the rock, σ is the compressive strength of the rock under multiaxial stress, P is the radial blasting load, and β is the attenuation exponent of the stress wave.In Equations (10) and (12), the range of rock blasting damage area is comprehensively affected by rock mechanical parameters, blasting load, and charge weight (diameter). It can be seen that the range of rock plastic area is a more comprehensive index for the influencing factors of blasting vibration dominant frequency, compared with charge weight. However, it is not easy to characterize the range of rock damage zone for cylindrical charge.In the process of excitation and propagation of blasting seismic waves, affected by the form of explosion source, the types and dominant components of blasting-induced seismic waves will be different. Different types of seismic waves will inevitably lead to the differences of blasting vibration spectrum and dominant frequency due to their differences in propagation speed and attenuation characteristics. According to the different propagation paths of the wave, it can be divided into body wave and surface wave. The propagation of body wave in rock mass can be divided into compressive wave (P wave) and shear wave (S wave), as shown in Figure11. In terms of the existence of the surface free surface, the surface wave (S wave) will also appear on the free surface of the semi-infinite space, which further increases the complexity of analyzing the dominant frequency of cylindrical charge blasting.Figure 11
Schematic diagram of single-hole blasting.


## 5. Conclusion

In this study, the numerical relationship between the elastic cavity radius and the dominant frequency induced by cylindrical charge has been explored. It has been found that(1)
with the increase of the elastic cavity radius, the zero-crossing and the Fourier dominant frequency decreases, and the amplitude spectrum shifts to the low frequency part, which is similar to the regulation of spherical charge obtained by theoretical analysis(2)
both the increase of charge diameter and charge length will increase the plastic zone induced by cylindrical charge, but as a result of the more complex charge structure and load energy release mechanism of cylindrical charge compared with spherical charge, the charge diameter and charge length will cause the plastic zone to change in the different directions of the explosive(3)
the range of plastic zone of cylindrical charge also has an approximate logarithmic relationship with the dominant frequency. However, since the plastic zone of cylindrical charge is less regular than that of spherical charge, the charge weightQ is a more favorable parameter compared with the elastic cavity radius a when used as the main influencing factor and attenuation analysis on the dominant frequency of cylindrical charge blasting



---

*Source: 9017792-2022-08-05.xml*