# Calculation Model and Decoupling Coefficient Sensitivity Study of Periphery Hole for Eccentric Decoupled Charge in Highway Tunnels

Authors: Yuanqing Wang; Junfeng Liu
Journal: Shock and Vibration
 (2018)

---

## Abstract

On the basis of the eccentric decoupled charge form on site and borehole concentric decoupled computational theory, the calculation model of eccentric decoupling charge would be established in the present study. The calculation model takes into account the influence factors such as damage to the surrounding rock, air, and water decoupled media. Besides, based on the proposed calculation model, blasting effect of a small-diameter EJ-102 emulsion explosive is calculated by using ANSYS LS-DYNA. The influence of decoupling coefficient, different rock grade, and different RQD value on blasting effect is obtained by sensitivity analysis of the simulated parameters. Then the calculated results are compared with the empirical results. The blasting parameters of fractured rock mass are also optimized by numerical simulation.


---

## Body

## 1. Introduction

Overexcavation is a common phenomenon in the blasting construction of highway tunnel so far. The overexcavation has a serious negative effects on the structural stability of the rock mass around the cavern [1, 2]. Besides, the overexcavation of road tunnel will inevitably affect the field construction of tunnel project [3, 4]. When the smooth blasting layer is thin, the blast split will occur before the smooth perforative cracks between the peripheral holes, and the clef surface formed inside the rock mass at the same time, which will lead to overexcavation [5–7].In order to avoid or alleviate the overexcavation phenomenon, the decoupled charge is always adopted for most tunnel periphery hole charge methods [8, 9]. Since the diameter of stick dynamite is generally smaller than that of the blast-hole, stick dynamite will be eccentric to the hole as the gravity force. The eccentric decoupled charge structure is formed then [1–3]. Wei [10] analyzed the eccentric decouple charge structure theoretically and considered that the crack is firstly produced at the contact point between the explosive and hole wall, and then it played an inhibitory role in the initiation of other points. Liu and Li [11] deduced the calculation formula of eccentric decouple charge, and the result was verified by comparison with the actual engineering application. Besides, Yu et al. [12] studied the explosion characteristics of water medium eccentric decouple charge and finally obtained the optimum formula of charge position. Pu [13] studied the cement mortar model test of eccentric decouple charge and concentric charge and shows that the initial eccentric stress field is formed on the surrounding rock of the borehole by an eccentric decoupling charge. Furthermore, Zhang [14] studied the explosive stress field distribution of the eccentric decoupled charge structure by means of holographic dynamic photoelastic test. The results show that the charge structure has a great influence on the stress field of the borehole. Zhang et al. [15] carried out field tests on the smooth blasting of eccentric decoupling charge. The results show that the initial pressure of the charging hole is 5~10 times away from the charging side. Guan et al. [16] used plexiglass to investigate the crack propagation law of eccentric decoupled charge and they found that when the eccentricity decoupling coefficient is 1.71, the difference of crack reaches the maximum between the coupled and decoupled sides. Shi et al. [17] carried out field and laboratory tests on mortar specimens; the results suggest the initial pressure of the hole at the charging side of the blast-hole is much larger than that of the noncharging side. Ханукаев [18] pointed out that the heterogeneity, isotropy, own defects, joints, and faults of the rock mass have great influence on the performance of rock blasting. Fourney et al. [19] also carried out further research on eccentric decoupling charge. However, these studies assume that the surrounding rock of the borehole is homogeneous and isotropic, without considering the heterogeneity of the surrounding rock, the initial damage of the surrounding rock, and the damage of surrounding rock caused by the blasting. At the same time, the influence of the decoupling medium of the borehole is not considered.In the present study, the calculation model of eccentric decoupling charge is established based on the eccentric decoupled charge form on site and borehole concentric decoupled computational theory. It should be noted that the present calculation model takes into account the influence of the factors including damage to the surrounding rock, air, and water decoupled media. Using ANSYS LS-DYNA, the proposed calculation model is applied to several case studies. Finally, the calculated results are compared with the empirical results and the blasting parameters of fractured rock mass are optimized by numerical simulation.


## 2. Calculation Model of Peripheral Hole Eccentric Decoupled Charge

Figure1 shows the main influence factors of the eccentric decouple charge, which are the radius of blast-hole r0, the stick dynamite radius rcs, and the type of noncoupling medium.Figure 1
Calculation diagram of eccentric decouple charge, whereo is the geometrical center of the hole, o1 is the geometrical center of the stick dynamite, r0 is the radius of blast-hole, rs is the straight-line distance from the center of the stick dynamite to the wall of blast-hole, and rcs is the radius of stick dynamite.The stress field of surrounding rock with eccentric decouple charge is studied based on the theory of concentric noncoupling charge. The stress field is symmetrized to the stick dynamite center and the adjacent holes. By analysing the stress filed of wave stress, it can be obtained that stress wave field distributes and spreads symmetrized to blast-hole and stick dynamite. As presented in Figure1, the hole pressure is changed with the variation of the roll center hole wall and radial distance, and the pressure of the wall PH for the blast-hole is calculated based on the concentric noncoupling charge blasting by(1)PH=P0×r0rcs-k,where PH is the peak value of the initial incident pressure of surrounding rock, P0=ρeD2/8 is initial average detonation pressure, r0 and rcs are the radius of the blast-hole and the dynamite stick, respectively, D and ρe stand for the bursting speed and density of the explosive material, k is the parameter related to the decoupled medium, and the value of k is taken as 6 and 0.72 for air and water, respectively [20, 21].The pressure value in the hole with the eccentric coupling can be obtained according to the calculation formula of hole wall pressure at the blasting of concentric decoupled charge in this paper. The following can be obtained based on (1):(2)Pcs=P0·rcsrsk=P0·rcsrs0.72waterP0·rcsrs6air,(3)rs=r02+r0-rcs2-2·r0·r0-rcs·cos⁡π-θ.Combining (2) and (3), it can be derived that(4)Pcs=P0·rcsrsk=ρeD22·κ+1·rcsr02+r0-rcs2-2·r0·r0-rcs·cos⁡π-θ0.72waterρeD22·κ+1·rcsr02+r0-rcs2-2·r0·r0-rcs·cos⁡π-θ6air,where Pcs is the initial pressure of the surrounding rock with eccentric decoupling charge. Then based on (4), the distribution of the bore pressure can be obtained (see Figure 2).Figure 2
Initial blast-hole pressure value diagram of eccentric decouple charge structure.It can be found that inside the blast-hole, the pressure value is larger in the side of the stick dynamite. As presented in Figure2, the pressure in B point reaches the maximum value, and the pressure of A point away from the loading side is smaller, which reaches the minimum value.The distribution of the incident pressure of the surrounding rock can be calculated by the following methods [12, 14]:(5)Pb=n·Pcs,where Pb is the initial transmission pressure value of surrounding rock and n is the power augmentation factor, which is generally taken between the values of 8 and 11.The variation law of radial stress amplitude and circumferential stress amplitude of blast stress wave is as follows [22]:(6)σrmax=pd0.5drα+σr0r,(7)σθmax=-pdb0.5drα+σθ0r,(8)b=b^=-1Shockwaveactionzone0.8μs1-0.8μsNonshockwaveactionzone,(9)α=3Shockwaveactionzone1-D00.5ρ0c^λ1λ3+λ2λ3Nonshockwaveaction,(10)D0=1-Ed0E^d,where b is the lateral pressure coefficient of stress wave under the damage condition, d is the diameter of the blast-hole, σr0r is the initial stress of surrounding rock, c^ is the longitudinal wave velocity under lossless condition, Ed0 and E^d are the initial elastic modulus and the undamaged elastic modulus of the rock mass, respectively, and λ1=-4.11×10-8, λ2=2.92, λ3 = 0.5~0.8 are the constants.Meanwhile, when the decoupled medium is water, due to the formation of explosion water vapor, the stress field of surrounding rock should also be considered:(11)pwr=pcs·rsr2=pcs·r02+r0-rcs2-2·r0·r0-rcs·cos⁡π-θr2.The criterion for the radius of rock failure zone is determined by formula (1)~(3). Moreover, bringing (6) and (7) into the criterion, the distribution of crushing zone of surrounding rock with an eccentric decoupled charge structure can be obtained.


## 3. Example Analysis

In this section, the blasting field test in Luo Jiang high-speed road, China, is taken as an example to be analyzed. The level of surrounding rock in the blasting field test is level II, and the lithology is sandstone. The blasting parameters are listed as follows: the hole lengthl is 4 m, the emulsion explosive and the detonation velocity of explosive are D=4000 m/s, ρe=1.25×103 kg/m3, charge radius r0 is 0.02 m, density of rock ρc is 3.8×103 kg/m3, compressive strength of surrounding rock σjc is 40 MPa, and tensile strength σl is 2.91 MPa, Poisson’s ratio ν is 0.2, and longitudinal wave velocity of surrounding rock Cp is 4390 m/s. Air density ρa is 1.29 kg/m3, water density ρw is 1.0×103 kg/m3, t1=1 ms [23], and t3=8 ms [23]; supposing that t4=10 ms [24], crack length is taken as 0.001 m. The crack angle β is taken as 45°.


### 3.1. Distribution of Surrounding Rock Pressure

The incident pressure of surrounding rock is calculated and analyzed in this part.As shown in Figure3, the charging hole radius is 40 mm, and the cartridge radius is 20 mm. In addition, the EJ-102 emulsion explosive is adopted. In the form of an eccentric decoupled charge, the water and air are used in the decoupled media, and the calculated results of the incident pressure of the surrounding rock are shown in Figure 3. It can be seen that, for both water decoupling charge and air decoupling charge, the pressure inside the blast-hole increases from the far side to the near side of the charging site. It should be noted that, for water decoupling charge, the hole pressure distribution is more uniform. It can be seen from Figure 3 that the maximum value of the surrounding rock pressure of 2.5 MPa can be achieved at B point. Comparing the water decoupling medium charge and the air decoupling medium charge, we can also find that the pressure of the surrounding rock under the water decoupling medium charge is obviously larger than that of the air decoupling charge. Specifically, Point A, which is in the most distal lateral wall rock, reached a maximum ratio of 377.7 times (1.133/0.003). It can also be seen from Figure 3 that, from the charge side B to the noncharge side A, the surrounding rock pressure for air coupling medium owns faster decay ratio than that for water decoupling charge case. Thus, according to Figure 3, the transmission effect of the water decoupling medium charge on the blasting energy is better than that of the air decoupling medium.Figure 3
Initial blast-hole pressure value diagram of eccentric decouple charge structure when the charge radius is 20 mm, in which the dashed line and the dotted line are the inner pressure distribution of the blast-hole for the decoupling medium of water and air, respectively.Figure4 shows the distribution of the internal pressure when the water medium is eccentrically loaded. It can be seen from the diagram that the initial internal pressure value increases continuously with the increase of charge radius. Because Point B is a coupling charge, the initial entrance pressure value of the hole reaches the maximum value of 2.5 MPa. On the contrary, the incident pressure of the surrounding rock is the minimum. With the decrease of charge diameter, the difference of incident pressure between Point B and Point A is increasing. Figure 5 is the variation of the initial pressure value of the hole for different charge diameters in the case of eccentric charge in the air medium. It can be seen that the smaller the charge diameter is, the faster the initial pressure of the surrounding rock will decay from Point B of the loading side to Point A at the nonloading side. Similarly, because Point B is the coupling charge point, the initial pressure value of the surrounding rock of this point is unchanged, and the maximum value is 2.5 MPa. In addition, the calculating results of the present model are also in agreement with Chen et al.’s research [22].Figure 4
The blast-hole pressure value along with the change of charge radius when the decoupling medium is water, where “d” represents the diameter of the explosive and the dashed line stands for the corresponding inner pressure distribution of the blast-hole.Figure 5
The blast-hole pressure value along with the change of charge radius when the decoupling medium is air, where “d” represents the diameter of the explosive and the dashed line stands for the corresponding inner pressure distribution of the blast-hole.


### 3.2. Distribution of Fractured Zone of Surrounding Rock

A total of four conditions of fractured zone for the surrounding rock are compared and discussed as follows in this section: Conditions 1: the results of this paper for water decoupling charge; Conditions 2: the results of this paper for air decoupling charge; Conditions 3: the calculation results of literature [22] for water decoupling charge; Conditions 4: the calculation results of literature [22] for air decoupling charge.Figure6 is a comparison between the theoretical calculation results of the present study and the results of literature [22] when the charge radius is 20 mm and the eccentric decouple coefficient is 2. By contrast, it can be found that the radius of the fracture zone formed by water medium charge is larger than that of the air medium charge, which is also in agreement with the result of literature [22]. For air decoupling charge, both the result of literature [22] and the numerical results indicate that the fracture zone formed mainly in the region of charge side and account for about 12~15 times of the blast-hole diameter, whereas, for water decoupling charge condition, distribution of the rupture zone around the blast-hole is more uniform, and the fracture zone accounts for about 15~23 times of the blast-hole diameter. The phenomenon above is for the reason that better transferring effect of the blasting energy is produced for water decoupling charge. Meanwhile, it can also be believed that the blast synergistic effect of water decoupling charge is more obvious than that for air decoupling charge.Figure 6
The blast-hole rupture zone when the charge radius is 20 mm, where “1” and “3” stand for the rupture zone around the blast-hole of the present study and Chen et al.’s result when the decoupling medium is water, and “2” and “4” stand for the rupture zone around the blast-hole of the present study and Chen et al.’s result when the decoupling medium is air.


## 3.1. Distribution of Surrounding Rock Pressure

The incident pressure of surrounding rock is calculated and analyzed in this part.As shown in Figure3, the charging hole radius is 40 mm, and the cartridge radius is 20 mm. In addition, the EJ-102 emulsion explosive is adopted. In the form of an eccentric decoupled charge, the water and air are used in the decoupled media, and the calculated results of the incident pressure of the surrounding rock are shown in Figure 3. It can be seen that, for both water decoupling charge and air decoupling charge, the pressure inside the blast-hole increases from the far side to the near side of the charging site. It should be noted that, for water decoupling charge, the hole pressure distribution is more uniform. It can be seen from Figure 3 that the maximum value of the surrounding rock pressure of 2.5 MPa can be achieved at B point. Comparing the water decoupling medium charge and the air decoupling medium charge, we can also find that the pressure of the surrounding rock under the water decoupling medium charge is obviously larger than that of the air decoupling charge. Specifically, Point A, which is in the most distal lateral wall rock, reached a maximum ratio of 377.7 times (1.133/0.003). It can also be seen from Figure 3 that, from the charge side B to the noncharge side A, the surrounding rock pressure for air coupling medium owns faster decay ratio than that for water decoupling charge case. Thus, according to Figure 3, the transmission effect of the water decoupling medium charge on the blasting energy is better than that of the air decoupling medium.Figure 3
Initial blast-hole pressure value diagram of eccentric decouple charge structure when the charge radius is 20 mm, in which the dashed line and the dotted line are the inner pressure distribution of the blast-hole for the decoupling medium of water and air, respectively.Figure4 shows the distribution of the internal pressure when the water medium is eccentrically loaded. It can be seen from the diagram that the initial internal pressure value increases continuously with the increase of charge radius. Because Point B is a coupling charge, the initial entrance pressure value of the hole reaches the maximum value of 2.5 MPa. On the contrary, the incident pressure of the surrounding rock is the minimum. With the decrease of charge diameter, the difference of incident pressure between Point B and Point A is increasing. Figure 5 is the variation of the initial pressure value of the hole for different charge diameters in the case of eccentric charge in the air medium. It can be seen that the smaller the charge diameter is, the faster the initial pressure of the surrounding rock will decay from Point B of the loading side to Point A at the nonloading side. Similarly, because Point B is the coupling charge point, the initial pressure value of the surrounding rock of this point is unchanged, and the maximum value is 2.5 MPa. In addition, the calculating results of the present model are also in agreement with Chen et al.’s research [22].Figure 4
The blast-hole pressure value along with the change of charge radius when the decoupling medium is water, where “d” represents the diameter of the explosive and the dashed line stands for the corresponding inner pressure distribution of the blast-hole.Figure 5
The blast-hole pressure value along with the change of charge radius when the decoupling medium is air, where “d” represents the diameter of the explosive and the dashed line stands for the corresponding inner pressure distribution of the blast-hole.


## 3.2. Distribution of Fractured Zone of Surrounding Rock

A total of four conditions of fractured zone for the surrounding rock are compared and discussed as follows in this section: Conditions 1: the results of this paper for water decoupling charge; Conditions 2: the results of this paper for air decoupling charge; Conditions 3: the calculation results of literature [22] for water decoupling charge; Conditions 4: the calculation results of literature [22] for air decoupling charge.Figure6 is a comparison between the theoretical calculation results of the present study and the results of literature [22] when the charge radius is 20 mm and the eccentric decouple coefficient is 2. By contrast, it can be found that the radius of the fracture zone formed by water medium charge is larger than that of the air medium charge, which is also in agreement with the result of literature [22]. For air decoupling charge, both the result of literature [22] and the numerical results indicate that the fracture zone formed mainly in the region of charge side and account for about 12~15 times of the blast-hole diameter, whereas, for water decoupling charge condition, distribution of the rupture zone around the blast-hole is more uniform, and the fracture zone accounts for about 15~23 times of the blast-hole diameter. The phenomenon above is for the reason that better transferring effect of the blasting energy is produced for water decoupling charge. Meanwhile, it can also be believed that the blast synergistic effect of water decoupling charge is more obvious than that for air decoupling charge.Figure 6
The blast-hole rupture zone when the charge radius is 20 mm, where “1” and “3” stand for the rupture zone around the blast-hole of the present study and Chen et al.’s result when the decoupling medium is water, and “2” and “4” stand for the rupture zone around the blast-hole of the present study and Chen et al.’s result when the decoupling medium is air.


## 4. Numerical Model

This section is based on the calculation model for the three-phase decoupling calculation model of gas, solid, and water. The cm-g-us unit system is adopted in the simulation model. The model size is400cm∗460 cm with 300 cm hole depth and 4 cm hole diameter. The specific schematic diagram of the model is depicted in Figure 7. The simulation scheme adopts the form of reverse detonation. In view of the fact that the surrounding rock mass is semi-infinite, the nonreflecting boundary condition is used for the right and left boundaries and the bottom boundary of the numerical model, and the top is a free surface. The actual blasting surrounding rock is an axisymmetric three-dimensional model, which is simplified by two-dimensional calculation model to facilitate calculations. Referring to the previous scholars’ research, the duration of shock wave, explosion gas, and explosion gas pressure after explosion is set to 300 us. Using ANSYS LS-DYNA and LS-PrePost-4.0, the model is preprocessed and after processed, respectively.Figure 7
The numerical simulation calculation model sketch of eccentric decouple charge (unit/cm).
(a)
Eccentric decoupling charge (b)
Eccentric uncoupling charge under horizontal jointsThe simulation of the fractured zones in the rock mass is introduced as follows. The preprocessing function of ANSYS is employed to simulate the fracture zones such as the joints and the flaws in the rock mass. To begin with, the position of the work plane is changed according to the actual joint orientation and direction, and the joint faces are cut in the analyzed geometric object, which is the necessary condition to produce the joint face. Then take advantage of the multiple loop command of ANSYS to assign different types of physical mechanics parameters to the model. Then the diversity of complexed structures for the established geometric object is reflected. Finally, the types of the units are chosen in the structural model framework formed in ANSYS. Meanwhile, the units of the geometric object are divided according to the joint development characteristics. Then the fractured zones inside the rock mass are formed, which is shown in Figure7.


### 4.1. The Selection of Blasting Materials

#### 4.1.1. Calculation Model Selection of Explosive Materials

The simulated explosive is made of HIGH_EXPLOSIVE_BURE material model and JWL equation of state, and the equation is [25](12)P=A1-ωR1Ve-R1V+B1-ωR2Ve-R2V+ωE0V,where P is the hole pressure value, A=214.4 GPa, B=0.182 GPa, R1=130, R2=0.90, and ω is the Green Elson parameters, which is taken as 0.15. V is the relative volume of detonation products, and E0 is the initial specific energy of the explosive, which is 3.15 GPa.


#### 4.1.2. The Material Model of Water and Air Bag

Grunneisen state equation is used for water and air [25]:(13)P=ρ0C2μ1+1-γ0/2μ-a/2μ21-S1-1μ-S2μ2/μ+1-S3μ3/μ+12+γ0+aμE0,where P is pore pressure and ρ0 is the density of the material, which is 1 g/cm3 for water and 0.001252 g/cm3 for air. γ0 is coefficient of Grunneisen, which is 0.35 for water and 1.4 for air. E0 is internal energy and is taken as 0 for both water and air. C, S1, S2, and S3 are all constants. a is the first-order correction quantity of γ0 and μ, and μ=ρ/ρ0-1.


### 4.2. Calculation Model of Rock Materials

In the present study, the elastoplastic model∗MAT_PLASTIC_KINEMATIC [26], which is shown as (14), is adopted. The rock parameters of the model are introduced as follows: density ρc is 3.8×103 kg/m3, elastic modulus E=25 GPa, and Poisson’s ratio is 0.2. The tangent modulus Et=1.5 G, hardening parameter ς=0.5, and failure strain ε=3.5. The calculating constants P and C are taken as 4.0 and 2.5, respectively.(14)σ=1+εC1/Pσ0+ςEtEE-Etεep.Figure8 is a numerical simulation result of decoupled air medium and water medium charge condition when the eccentric coefficient is 2. It is easy to find that the radius of the rock fracture zone for water medium decoupling charge is larger than that for air medium decoupling charge, suggesting that blast synergistic effect of water decoupling charge is more obvious than that for air decoupling charge.Figure 8
Numerical simulation results of eccentric decouple charge.
(a)
Condition of air medium charge with eccentric decoupling coefficient of 2 (b)
Simulation results of charge by air medium (c)
Condition of water medium charge with eccentric decoupling coefficient of 2 (d)
Simulation results of charge by water mediumTable1 shows the specific data of the rock fracture radius for the numerical simulation results and corresponding theoretical solution. It can be seen that the theoretical and numerical results are in line with each other, and the deviation is less than 15%. Therefore, it can be believed that the proposed eccentric decoupled numerical model is in good agreement with the theoretical study, which may have certain significance in studying the effect of eccentric decoupling charge.Table 1
Comparison of theoretical solution and numerical solution for eccentric decouple charge.
decoupling medium wall rock radius of broken rock zonerccm deviation±% theory solution numerical solution air lateral 2 2 0 medial 52.16 60 15.0 water lateral 76.6 68 12.6 medial 52.16 54 3.5


### 4.3. Study on the Influence of Eccentric Decouple Charge Factor on Overexcavation of Surrounding Rock

Table2 shows the numerical simulation results of the lateral fracture zone for the surrounding rock mass, and Figure 9 is the specific variation law of the fracture zone radius for the surrounding rock as the radial decoupling coefficient kr varies. Obviously, when the charging structure, the surrounding rock grade, and the RQD value are the same, numerical analysis of each group shows that, with the increase of the radial coupling coefficient, the lateral radius of fractured zone takes on a decreasing tendency. Also, when the eccentric decoupling coefficient is larger than or equal to 2, the radius of the lateral rock rupture zone is stabilized around the value of the blast-hole radius. In other words, there is no overexcavation outside the surrounding rock.Table 2
Numerical simulation results of the lateral fracture zone for the surrounding rock mass.
Radial decoupling coefficient radius of broken lateral rock zone (/cm) Surrounding rock grade RQD = 10 y% RQD = 35% RQD = 60% RQD = 85% RQD = 95% 4 -- -- 2 2 2 II -- 2 2 2 -- III 2 2 2 -- -- IV 2 -- -- 2 2 2 II -- 2 2 2 -- III 2 2 2 -- -- IV 1.33 -- -- 33 24 22 II -- 30 31 23 -- III 35 32 32 -- -- IV 1 -- -- 60 52 50 II -- 45 52 52 -- III 40 49 39 -- -- IVFigure 9
Relationship betweenrc of outside surrounding rock and kr.From the above, it can be concluded that under the precondition of no underexcavation, the charging structure with the eccentric decoupling coefficient larger than or equal to 2 should be possibly adopted.


### 4.4. Study about the Relation between the Fracture Zone Radius of Outside Surrounding Rock and the Rock Grade

From Table2 and, in accordance with Figure 10, which is the relationship between the fracture zone radius of outside surrounding rock (rc) and rock grade, it can be obtained from the numerical simulations that when the charge structure, radial eccentric coupling coefficient, and RQD value are in the same condition, the surrounding rock lateral radius of fractured zone decreases or keeps constant as the surrounding rock grade decreases.Figure 10
Relationship between fracture zone radius of outside surrounding rock(rc) and rock grade.


### 4.5. Study on the Effect of Different RQD Values on the Blasting Effect

From the Table2, the variation law between the fracture zone radius of outside surrounding rock (rc) and the RQD value of the surrounding rock can be established, which is depicted as Figure 11. It can be seen that, as the RQD value increases, rc shows a decreasing tendency for surrounding rock of grade II, whereas, for surrounding rock of grades III and IV, rc increases slightly at first and then decrease and finally keeps constant. The results suggest that the grade of the surrounding rock take an obvious effect on the relation between fracture zone radius of outside surrounding rock and RQD value.Figure 11
Relationship between fracture zone radius of outside surrounding rock(rc) and RQD value.
(a)
Surrounding rock of grade II (b)
Surrounding rock of grade III (c)
Surrounding rock of grade IV


## 4.1. The Selection of Blasting Materials

### 4.1.1. Calculation Model Selection of Explosive Materials

The simulated explosive is made of HIGH_EXPLOSIVE_BURE material model and JWL equation of state, and the equation is [25](12)P=A1-ωR1Ve-R1V+B1-ωR2Ve-R2V+ωE0V,where P is the hole pressure value, A=214.4 GPa, B=0.182 GPa, R1=130, R2=0.90, and ω is the Green Elson parameters, which is taken as 0.15. V is the relative volume of detonation products, and E0 is the initial specific energy of the explosive, which is 3.15 GPa.


### 4.1.2. The Material Model of Water and Air Bag

Grunneisen state equation is used for water and air [25]:(13)P=ρ0C2μ1+1-γ0/2μ-a/2μ21-S1-1μ-S2μ2/μ+1-S3μ3/μ+12+γ0+aμE0,where P is pore pressure and ρ0 is the density of the material, which is 1 g/cm3 for water and 0.001252 g/cm3 for air. γ0 is coefficient of Grunneisen, which is 0.35 for water and 1.4 for air. E0 is internal energy and is taken as 0 for both water and air. C, S1, S2, and S3 are all constants. a is the first-order correction quantity of γ0 and μ, and μ=ρ/ρ0-1.


## 4.1.1. Calculation Model Selection of Explosive Materials

The simulated explosive is made of HIGH_EXPLOSIVE_BURE material model and JWL equation of state, and the equation is [25](12)P=A1-ωR1Ve-R1V+B1-ωR2Ve-R2V+ωE0V,where P is the hole pressure value, A=214.4 GPa, B=0.182 GPa, R1=130, R2=0.90, and ω is the Green Elson parameters, which is taken as 0.15. V is the relative volume of detonation products, and E0 is the initial specific energy of the explosive, which is 3.15 GPa.


## 4.1.2. The Material Model of Water and Air Bag

Grunneisen state equation is used for water and air [25]:(13)P=ρ0C2μ1+1-γ0/2μ-a/2μ21-S1-1μ-S2μ2/μ+1-S3μ3/μ+12+γ0+aμE0,where P is pore pressure and ρ0 is the density of the material, which is 1 g/cm3 for water and 0.001252 g/cm3 for air. γ0 is coefficient of Grunneisen, which is 0.35 for water and 1.4 for air. E0 is internal energy and is taken as 0 for both water and air. C, S1, S2, and S3 are all constants. a is the first-order correction quantity of γ0 and μ, and μ=ρ/ρ0-1.


## 4.2. Calculation Model of Rock Materials

In the present study, the elastoplastic model∗MAT_PLASTIC_KINEMATIC [26], which is shown as (14), is adopted. The rock parameters of the model are introduced as follows: density ρc is 3.8×103 kg/m3, elastic modulus E=25 GPa, and Poisson’s ratio is 0.2. The tangent modulus Et=1.5 G, hardening parameter ς=0.5, and failure strain ε=3.5. The calculating constants P and C are taken as 4.0 and 2.5, respectively.(14)σ=1+εC1/Pσ0+ςEtEE-Etεep.Figure8 is a numerical simulation result of decoupled air medium and water medium charge condition when the eccentric coefficient is 2. It is easy to find that the radius of the rock fracture zone for water medium decoupling charge is larger than that for air medium decoupling charge, suggesting that blast synergistic effect of water decoupling charge is more obvious than that for air decoupling charge.Figure 8
Numerical simulation results of eccentric decouple charge.
(a)
Condition of air medium charge with eccentric decoupling coefficient of 2 (b)
Simulation results of charge by air medium (c)
Condition of water medium charge with eccentric decoupling coefficient of 2 (d)
Simulation results of charge by water mediumTable1 shows the specific data of the rock fracture radius for the numerical simulation results and corresponding theoretical solution. It can be seen that the theoretical and numerical results are in line with each other, and the deviation is less than 15%. Therefore, it can be believed that the proposed eccentric decoupled numerical model is in good agreement with the theoretical study, which may have certain significance in studying the effect of eccentric decoupling charge.Table 1
Comparison of theoretical solution and numerical solution for eccentric decouple charge.
decoupling medium wall rock radius of broken rock zonerccm deviation±% theory solution numerical solution air lateral 2 2 0 medial 52.16 60 15.0 water lateral 76.6 68 12.6 medial 52.16 54 3.5


## 4.3. Study on the Influence of Eccentric Decouple Charge Factor on Overexcavation of Surrounding Rock

Table2 shows the numerical simulation results of the lateral fracture zone for the surrounding rock mass, and Figure 9 is the specific variation law of the fracture zone radius for the surrounding rock as the radial decoupling coefficient kr varies. Obviously, when the charging structure, the surrounding rock grade, and the RQD value are the same, numerical analysis of each group shows that, with the increase of the radial coupling coefficient, the lateral radius of fractured zone takes on a decreasing tendency. Also, when the eccentric decoupling coefficient is larger than or equal to 2, the radius of the lateral rock rupture zone is stabilized around the value of the blast-hole radius. In other words, there is no overexcavation outside the surrounding rock.Table 2
Numerical simulation results of the lateral fracture zone for the surrounding rock mass.
Radial decoupling coefficient radius of broken lateral rock zone (/cm) Surrounding rock grade RQD = 10 y% RQD = 35% RQD = 60% RQD = 85% RQD = 95% 4 -- -- 2 2 2 II -- 2 2 2 -- III 2 2 2 -- -- IV 2 -- -- 2 2 2 II -- 2 2 2 -- III 2 2 2 -- -- IV 1.33 -- -- 33 24 22 II -- 30 31 23 -- III 35 32 32 -- -- IV 1 -- -- 60 52 50 II -- 45 52 52 -- III 40 49 39 -- -- IVFigure 9
Relationship betweenrc of outside surrounding rock and kr.From the above, it can be concluded that under the precondition of no underexcavation, the charging structure with the eccentric decoupling coefficient larger than or equal to 2 should be possibly adopted.


## 4.4. Study about the Relation between the Fracture Zone Radius of Outside Surrounding Rock and the Rock Grade

From Table2 and, in accordance with Figure 10, which is the relationship between the fracture zone radius of outside surrounding rock (rc) and rock grade, it can be obtained from the numerical simulations that when the charge structure, radial eccentric coupling coefficient, and RQD value are in the same condition, the surrounding rock lateral radius of fractured zone decreases or keeps constant as the surrounding rock grade decreases.Figure 10
Relationship between fracture zone radius of outside surrounding rock(rc) and rock grade.


## 4.5. Study on the Effect of Different RQD Values on the Blasting Effect

From the Table2, the variation law between the fracture zone radius of outside surrounding rock (rc) and the RQD value of the surrounding rock can be established, which is depicted as Figure 11. It can be seen that, as the RQD value increases, rc shows a decreasing tendency for surrounding rock of grade II, whereas, for surrounding rock of grades III and IV, rc increases slightly at first and then decrease and finally keeps constant. The results suggest that the grade of the surrounding rock take an obvious effect on the relation between fracture zone radius of outside surrounding rock and RQD value.Figure 11
Relationship between fracture zone radius of outside surrounding rock(rc) and RQD value.
(a)
Surrounding rock of grade II (b)
Surrounding rock of grade III (c)
Surrounding rock of grade IV


## 5. Conclusions

In the present study, the analytical calculation model of eccentric decoupling charge is established based on the eccentric decoupled charge form on site and the concentric decoupled computational theory. Moreover, the influence factors including damage to the surrounding rock, air, and water decoupled media are taken into consideration in the present calculation model. Using ANSYS LS-DYNA, the proposed calculation model is applied to a case study. Finally, the following conclusions are obtained:( 1 ) Comparing with the theoretical model, we can learn that the eccentric decoupling charge of air medium has better control of the overbreak-underbreak effect of the periphery hole than the eccentric uncoupling charge of water medium.( 2 ) Under the precondition of no underexcavation, the charging structure with the eccentric decoupling coefficient larger than or equal to 2 should be possibly adopted.( 3 ) Under the same condition of other charging parameters, the radius of the fracture area of surrounding rock decreases with the decrease of the surrounding rock grade, and the change is always decreasing or stable.( 4 ) Under the same condition of other charging parameters, the radius of the fracture area of surrounding rock increases with the increase of RQD value of surrounding rock. The determination of specific charge parameters can be determined by the numerical simulation diagram according to the actual conditions.



---

*Source: 9734529-2018-05-17.xml*