# Numerical Simulation of Blast Vibration and Crack Forming Effect of Rock-Anchored Beam Excavation in Deep Underground Caverns

Authors: XinPing Li; JunHong Huang; Yi Luo; Qian Dong; YouHua Li; Yong Wan; TingTing Liu
Journal: Shock and Vibration
 (2017)

---

## Abstract

Aiming at surrounding rock damage induced by dynamic disturbance from blasting excavation of rock-anchored beam in rock mass at moderate or far distance in underground cavern, numerical model of different linear charging density and crustal stress in underground cavern is established by adopting dynamic finite element software based on borehole layout, charging, and rock parameter of the actual situation of a certain hydropower station. Through comparison in vibration velocity, contour surface of rock mass excavation, and the crushing extent of excavated rock mass between calculation result and field monitoring, optimum linear charging density of blast hole is determined. Studies are also conducted on rock mass vibration in moderate or far distance to blasting source, the damage of surrounding rock in near-field to blasting source, and crushing degree of excavated rock mass under various in situ stress conditions. Results indicate that, within certain range of in situ stress, the blasting vibration is independent of in situ stress, while when in situ stress is increasing above certain value, the blasting vibration velocity will be increasing and the damage of surrounding rock and the crushing degree of excavated rock mass will be decreasing.


---

## Body

## 1. Introduction

Jinsha River Basin, which lies on the upstream of Yangtze River in China, has witnessed the development and construction of large numbers of water conservancy and hydropower projects. Due to the fact that most basin area of Jinsha River locates on the towering and steep mountains, the work of layout, construction, and operation of building hub should be conducted by the framework of underground cavern. Meanwhile, a certain amount of underground caverns in a large-scale and ultra-large type appear along with the enlargement of hydropower in scale. The rock-anchored beam is frequently used as the load bearing of crane in caverns in order to reduce the span of underground powerhouses and accelerate the progress of construction.The rock-anchored beam is widely used in hydropower projects as the essential construction in underground caverns which makes the first success in Norway [1, 2]. There is certainly no need for the existence of crane columns or walls in the application of rock-anchored beam since it can take full advantage of the bearing capacity of surrounding walls. The application of rock-anchored beam has several advantages. It can not only reduce the excavation span of underground caverns which benefits the surrounding rock stability and the amount of excavation of tunnels but also make it possible to install the crane in advance and brings convenience for concrete placements of caverns and unit installation. Thus, the construction progress could be accelerated and the investment could be saved, bringing remarkable economical benefits. The application of rock-anchored beam is widespread in many large-scale overseas hydropower stations [3].However, with considerable difficulty, the excavation of rock-anchored beams is in high quality requirements, which directly affects the operation conditions of bridge cranes after construction. In order to make cracks form successfully in adjacent blast holes and reduce the damage of remaining rock mass to minimum, the construction site usually meets the requirements of excavation by means of controlling the charging of the blast hole and blast holes spacing.Considering the noncoupling charge structure applied in driving blasting, it is the most reasonable to adopt the fluid-structure interaction algorithm to simulate the blasting process [4]. In the simulation of explosion process, by adopting Euler algorithm when it comes to the explosive and noncoupling air in blast holes, explosive could mingle with air unit as fluid. Meanwhile, by adopting Lagrange algorithm aiming at mass rock structures, the interaction between fluid and solid is dealt with by means of fluid-structure coupling. The advantages of this kind of algorithm mare that explosive and fluid materials flow in the Euler unit without the problem of element distortion and the interaction between fluid and solid is dealt with by means of fluid-structure coupling [5, 6]. When the explosive blast hole blasts, explosive stress wave and detonation gas take effects on the surrounding rock, which makes the cracks on blast hole walls crack and drives their expansion, therefore, forming the penetrating cracks on lines of centers of blast holes and making detonation gas pressure acting on both sides of the cracks.But the aforesaid algorithms could only simulate the shock response of the near-field of blasting source. There is still a way to go for the simulation of the whole process from the blasting of near-field of blasting source to the blasting vibration spreading of explosion of far-field [7]. The reason is that the explosive material units and meshes around blast holes in the model needed in such algorithm are mapped incredibly small. Furthermore, the calculation model of numerical simulation of explosive vibration field is quite large. All the above factors lead to the heavy amounts of units after the meshing, which would further hamper the progress of calculation. In order to meet the demands of practical application, an equivalent simulation is strongly needed.In recent years, Ma et al. [8] and Wu et al. [9] managed to ascertain the changing course of blast holes pressure through calculation as well as with the assistance of laboratory and field test. However, there is still a way to go for these achievements to resolve the problem of blast holes pressure theoretically. In practical application of programs, the selecting of changing courses of blast holes pressures is still semitheoretical and semiexperimental exponential decay load. Generally speaking, exerting the explosion load on the blast holes walls or exerting the explosion load on the surface where lines of centers of blast holes lie on are two frequently used equivalent simulation methods [10, 11]. According to Saint-Venant principle, the latter method is more reasonable and effective, which adopts triangle pulse waves as equivalent blasting loads to conduct the three-dimensional simulations.The paper is based on the basic mechanical properties and blasting parameters of rock mass in the process of excavation program of underground caverns in Baihetan hydropower station in Yunnan, China. First, the research adopts three-dimensional models in the calculation of vibration speed of the blasting of rock mass in the moderate and far distance and two-dimensional models in the simulation of crushing conditions of rock mass in near-field of blasting sources and explosion. At the same time, the research also includes the comparison of contour surface and vibration speed in the measuring point between the model and the filed. Then, the paper does the study of rock mass vibration in the moderate or far distance of blasting, the damage of surrounding rock in near-field of blasting source, and crushing degree of excavated rock mass under different crustal stress.


## 2. The Calculation Method of Rock Mass Response under Blasting Load

This article uses dynamic finite elements analysis software ANSYS/LS-DYNA, which is a typical dynamic finite element calculation software. It does not open source, so the programming code can not be downloaded. The solving equation used in LS-DYNA for dynamic calculation is motion equation:(1)Mu¨+Cu˙+Ku=f,where M is global mass matrix, C is damping matrix, K is global stiffness matrix, u¨, u˙, and u are nodes’ acceleration vector, velocity vector, and displacement vector, respectively, and f is equivalent nodal force vector [13].ANSYS uses Newmark integration method to solve kinetic equation. It’s basic idea is to change the requirement that any timet should satisfy the displacement vector ut of (1) to only require time discrete points t,t+Δt,… to satisfy the above equation of motion.The solving equation used for static calculation is as follows:(2)Kδ=W,where K is global stiffness matrix, δ is displacement field of rock mass, and W is stress field of rock mass.Comparing (1) with (2), the solving equation used in static calculation does not contain mass matrix, so it can not reflect the influences that inertia force has on surrounding rock deformation in the whole process of tunnel blast excavation. Xie et al. [14] performed further research about the damage situation of tunnel in the process of blasting excavation by LS-DYNA.


### 2.1. Equivalent Load Application Method

Due to the fact that the blast holes studied in this paper are small in size and large in amount, it is inconvenient to build the blast hole model through the method of explosion load exertion in three-dimensional model. According to the equilibrium principle, exerting the peak pressure of shock waves on the blast holes walls is equivalent to the value on lines of centers of blast holes [15]. The disturbances of rock mass under explosion load mainly depend on the peak load produced by blasting and the affecting time of explosion load. The paper simplifies the explosion load curve to triangle load which is exerted on blast holes walls. The explosion load curve is shown in Figure 1.Figure 1
The curved diagram of explosion load.Overseas scholars have come up with C-J theory, which simplifies the detonation wave to strong discontinuity surface contained chemical reactions and gives the conclusion of stable spread of detonation wave under the circumstances of C-J conditions. The simplified load peak valuepb exerting on blast holes is(3)pb=n·ρ0D221+γdcdb2γlelbγ,where ρ0 represents the density of explosive (g/cm3); D represents explosive velocity (m/s); γ represents the isentropic exponent of explosive, usually takes 3; dc represents the charge diameter; db represents borehole diameter; le represents the length of blast holes; lb represents the length of charge column; here smooth blasting holes take le=0.6lb and auxiliary holes take le=0.8lb; n represents pressure increasing factor when the detonation gas expansion collides with hole wall, n= 8–11, and it takes 10 for emulsion explosive.As shown in Figure2, the equivalent peak pressure calculated through force balance principle is(4)pe=2r0apb,where r0 represents blast hole radius; a represents distances between blast holes.Figure 2
The equivalent load exerting on lines of centers of blast holes.


### 2.2. The Principle of Rock Mass Damage under Explosion Action

The principle of rock mass damage depends on the property of rock mass as well as the practical force conditions. The pressure of rock mass, taking the Mises damage rule, forms the crushing area of rock mass blasting, while the cracks area is the result of the damage of tensile force. The damage rule of rock mass is as follows:(5)σVM>σcd,crushs  area,σt>σtd,cracks  area,σVM=32σijσij,where σVM is the von Mises effective stress of any point in rock mass; σij(i,j=1,2,3) is the stress components of rock mass; σt is the tensile stress of explosion load of any point in rock mass; σcd,σtd are known as uniaxial dynamic compressive strength and tensile strength of rock mass, respectively.In order to reduce the damage of surrounding rock on site, grooving has been done on part of blast holes as shown in Figure3. The design of grooved blast holes in models according to the actual grooved blast holes sizes on site is shown in Figure 4. The next Figure 5 is explosion load figure of the groove borehole wall and the circular part of borehole wall is main pressure effect, while the groove tip of borehole wall is mainly tensile effect.Figure 3
Grooved blast holes.Figure 4
Grooved blast holes size.Figure 5
Hole wall of groove hole load distribution.As shown in Figure5, during the detonation of explosives in borhole, tensile stress would be concentrated at groove tip A instantaneously, and crack initial from element A must satisfy the formula(6)T>σtd,T=Prcos⁡α2.The dynamic compressive stress of rock increases with the improvement of loaded strain rate, generally approximated by the following equation [16]:(7)σcd=σcε˙1/3,where σc in the equation refers to the uniaxial static compressive stress of rock.The loaded strain rate of rockε˙ in blasting during programs is within 100–105 s−1, among which the strain rate in crushed zones could be ε˙=102–104s-1 and in cracked zones could be ε˙=100–103s-1.For the lack of corresponding analytical data of experiments and theories, the value of dynamic tensile strength approximates(8)σtd=σtε˙1/3,where σt in the equation refers to the uniaxial static tensile strength of rock mass.The paper applies the fluid-structure interaction algorithm on the calculation of explosive, rock, and air through two-dimensional model, taking the tensile strength of rock as the criteria of rock mass fracture, in order to observe the damage of surrounding rock and the crushing condition of excavated rock mass under the circumstances of noncoupling charge. The damage and crushing rate of rock mass under the effect of blast by explosive have long been researched. Both domestic and overseas scholars have drawn quite a few remarkable conclusions through experimental research and theoretical analysis. However, these conclusions differ from one another due to the differences between experimental conditions. The paper will research on the damage of rock mass and crushing rate under different crustal stress with the reference of actual operation situation on site.


### 2.3. Plastic Kinematic Strengthened Constitutive Model

In rock medium, the analysis of static problems is based on isotropic strengthened constitutive model, while the analysis of cyclic loading and dynamic issue is based on two constitutive models, namely, the kinematic strengthened constitutive and the mixed strengthened constitutive models. Particularly, during the blasting process, under the conditions of comparable large strain of rock mass in near-field, the effect of strain rate is quite apparent. Considering the aforesaid factors, it is more appropriate to build the model of plastic kinematic strengthened constitutive in the analysis of dynamic problems of underground blasting programs, for the model conforms with the conditions in real project. The strain rate is considered in Cowper-Symonds model; the yield stress is expressed in factors related to the strain rate [17]; and the criteria of material which yields in the model of plastic kinematic strengthened constitutive are as follows:(9)ϕ=σi2-σy2=0,σy=1+ε˙C1/Pσ0+βEpεpeff,Ep=E0EtanE0-Etan.In the formula, ε˙ is strain rate and C and P are relative parameters in Cowper-Symonds, which are generally determined by features of material strain rate. The paper refers to the conclusion in reference [17], which indicates that C is 2.5 s−1, P is 4, σ0 represents initial yield stress of rock, and β is the hardening coefficient. When β=0, β represents plastic kinematic hardening; when β=1, it represents isotropic hardening; when 0<β<1, it represents mixed hardening. The paper adopts the plastic following hardening in calculation; therefore, β=0 is taken. Ep, E0, and Etan represent plastic strain-hardening modulus, Young’s modulus, and tangent modulus separately.When the material is in hardening stage, the efficient plastic strain rateεeffp corresponded could be expressed as follows:(10)εeffp=∫0t23ε˙ijpε˙ijp1/2dt,ε˙effp=ε˙ij-ε˙ije,where ε˙ijp and ε˙ije represent plastic strain rate and elastic strain rate separately.


## 2.1. Equivalent Load Application Method

Due to the fact that the blast holes studied in this paper are small in size and large in amount, it is inconvenient to build the blast hole model through the method of explosion load exertion in three-dimensional model. According to the equilibrium principle, exerting the peak pressure of shock waves on the blast holes walls is equivalent to the value on lines of centers of blast holes [15]. The disturbances of rock mass under explosion load mainly depend on the peak load produced by blasting and the affecting time of explosion load. The paper simplifies the explosion load curve to triangle load which is exerted on blast holes walls. The explosion load curve is shown in Figure 1.Figure 1
The curved diagram of explosion load.Overseas scholars have come up with C-J theory, which simplifies the detonation wave to strong discontinuity surface contained chemical reactions and gives the conclusion of stable spread of detonation wave under the circumstances of C-J conditions. The simplified load peak valuepb exerting on blast holes is(3)pb=n·ρ0D221+γdcdb2γlelbγ,where ρ0 represents the density of explosive (g/cm3); D represents explosive velocity (m/s); γ represents the isentropic exponent of explosive, usually takes 3; dc represents the charge diameter; db represents borehole diameter; le represents the length of blast holes; lb represents the length of charge column; here smooth blasting holes take le=0.6lb and auxiliary holes take le=0.8lb; n represents pressure increasing factor when the detonation gas expansion collides with hole wall, n= 8–11, and it takes 10 for emulsion explosive.As shown in Figure2, the equivalent peak pressure calculated through force balance principle is(4)pe=2r0apb,where r0 represents blast hole radius; a represents distances between blast holes.Figure 2
The equivalent load exerting on lines of centers of blast holes.


## 2.2. The Principle of Rock Mass Damage under Explosion Action

The principle of rock mass damage depends on the property of rock mass as well as the practical force conditions. The pressure of rock mass, taking the Mises damage rule, forms the crushing area of rock mass blasting, while the cracks area is the result of the damage of tensile force. The damage rule of rock mass is as follows:(5)σVM>σcd,crushs  area,σt>σtd,cracks  area,σVM=32σijσij,where σVM is the von Mises effective stress of any point in rock mass; σij(i,j=1,2,3) is the stress components of rock mass; σt is the tensile stress of explosion load of any point in rock mass; σcd,σtd are known as uniaxial dynamic compressive strength and tensile strength of rock mass, respectively.In order to reduce the damage of surrounding rock on site, grooving has been done on part of blast holes as shown in Figure3. The design of grooved blast holes in models according to the actual grooved blast holes sizes on site is shown in Figure 4. The next Figure 5 is explosion load figure of the groove borehole wall and the circular part of borehole wall is main pressure effect, while the groove tip of borehole wall is mainly tensile effect.Figure 3
Grooved blast holes.Figure 4
Grooved blast holes size.Figure 5
Hole wall of groove hole load distribution.As shown in Figure5, during the detonation of explosives in borhole, tensile stress would be concentrated at groove tip A instantaneously, and crack initial from element A must satisfy the formula(6)T>σtd,T=Prcos⁡α2.The dynamic compressive stress of rock increases with the improvement of loaded strain rate, generally approximated by the following equation [16]:(7)σcd=σcε˙1/3,where σc in the equation refers to the uniaxial static compressive stress of rock.The loaded strain rate of rockε˙ in blasting during programs is within 100–105 s−1, among which the strain rate in crushed zones could be ε˙=102–104s-1 and in cracked zones could be ε˙=100–103s-1.For the lack of corresponding analytical data of experiments and theories, the value of dynamic tensile strength approximates(8)σtd=σtε˙1/3,where σt in the equation refers to the uniaxial static tensile strength of rock mass.The paper applies the fluid-structure interaction algorithm on the calculation of explosive, rock, and air through two-dimensional model, taking the tensile strength of rock as the criteria of rock mass fracture, in order to observe the damage of surrounding rock and the crushing condition of excavated rock mass under the circumstances of noncoupling charge. The damage and crushing rate of rock mass under the effect of blast by explosive have long been researched. Both domestic and overseas scholars have drawn quite a few remarkable conclusions through experimental research and theoretical analysis. However, these conclusions differ from one another due to the differences between experimental conditions. The paper will research on the damage of rock mass and crushing rate under different crustal stress with the reference of actual operation situation on site.


## 2.3. Plastic Kinematic Strengthened Constitutive Model

In rock medium, the analysis of static problems is based on isotropic strengthened constitutive model, while the analysis of cyclic loading and dynamic issue is based on two constitutive models, namely, the kinematic strengthened constitutive and the mixed strengthened constitutive models. Particularly, during the blasting process, under the conditions of comparable large strain of rock mass in near-field, the effect of strain rate is quite apparent. Considering the aforesaid factors, it is more appropriate to build the model of plastic kinematic strengthened constitutive in the analysis of dynamic problems of underground blasting programs, for the model conforms with the conditions in real project. The strain rate is considered in Cowper-Symonds model; the yield stress is expressed in factors related to the strain rate [17]; and the criteria of material which yields in the model of plastic kinematic strengthened constitutive are as follows:(9)ϕ=σi2-σy2=0,σy=1+ε˙C1/Pσ0+βEpεpeff,Ep=E0EtanE0-Etan.In the formula, ε˙ is strain rate and C and P are relative parameters in Cowper-Symonds, which are generally determined by features of material strain rate. The paper refers to the conclusion in reference [17], which indicates that C is 2.5 s−1, P is 4, σ0 represents initial yield stress of rock, and β is the hardening coefficient. When β=0, β represents plastic kinematic hardening; when β=1, it represents isotropic hardening; when 0<β<1, it represents mixed hardening. The paper adopts the plastic following hardening in calculation; therefore, β=0 is taken. Ep, E0, and Etan represent plastic strain-hardening modulus, Young’s modulus, and tangent modulus separately.When the material is in hardening stage, the efficient plastic strain rateεeffp corresponded could be expressed as follows:(10)εeffp=∫0t23ε˙ijpε˙ijp1/2dt,ε˙effp=ε˙ij-ε˙ije,where ε˙ijp and ε˙ije represent plastic strain rate and elastic strain rate separately.


## 3. The Numerical Model of Blasting Excavation of Rock-Anchored Beam

Baihetan hydropower station, located in the upstream of Jinsha River in southwest China, has a capacity of 2060 billion cubic meters in reservoir and the underground powerhouses install 16 hydroelectric generating sets with the capacity per set of 1000 MW. The installed capacity drafted in the beginning is 16 million KW, and the average annual electric energy production is 6024 billion KWH. The station will become the second largest hydropower station which is next only to Three Gorges after completion, which makes the high demand of quality of rock-anchored beam of rock mass in underground plants. Thus, there is strong need for the fine control of blasting excavation of the rock-anchored beam in order to keep the damage of surrounding rock caused by blasting in a small extent.The charging parameter of blast holes in single blasting excavation on site is shown in Table1. Electronic detonator is used for initiation delay to reduce the blasting vibration and improve the effect of rock fracture in many engineering experiments [18, 19], and at present the minimum delay has already reached 1 ms, so I used MS1 and MS9 electronic detonators to achieve initiation delay for auxiliary blasting holes and smooth blasting holes. The plane figure of blast holes detonation network and the profile map of blast holes are shown in Figure 6.Table 1
Different blast hole layout and charge parameter.
Name
Hole diameter
Hole distance
Hole depth
Hole amount
Cartridge diameter
Charging density
① Vertical holes
Φ 42 mm
30 cm
248 cm
68
Φ 25 mm
65/70/85 g/m
② Oblique holes
Φ 42 mm
30 cm
260 cm
68
Φ 25 mm
65/70/85 g/m
③ Auxiliary holes
Φ 42 mm
90 cm
242 cm
22
Φ 25 mm
186/206/250 g/mFigure 6
The plane figure of blast holes detonation network and the profile map of blast holes (unit: cm).


### 3.1. The Three-Dimensional Model of the Blasting Excavation of Rock-Anchored Beam

Taking the blasting excavation of the protective layers of rock-anchored beams in the underground main power house of the right bank of the dam as an example, the blasting excavation of protective layers of rock-anchored beams drills the vertical smooth blasting holes, oblique smooth blasting holes, and auxiliary blast holes manually in the rate of 20.4 m excavation length each time. Meanwhile, the blasting vibration monitoring is conducted on the sidewalls near excavated rock-anchored beams. The measuring point vibration monitoring of near-field rock during the process of blasting excavation is conducted in order to make sure the disturbances of blasting excavation to surrounding rocks in a safety range. Meanwhile, for ascertaining the safety of monitoring equipment, the nearest measuring points are arranged 10 m away from the boundary of the excavation area, and then one measuring point is arranged every 5 m, and a total of five measuring points are arranged. It is as shown in Figure7.Figure 7
The schematic diagram of measuring point layout of blast vibration velocity (unit: m).The elevation of rock-anchored beams excavated on site ranges from EL. 602.40 m to EL. 606.90 m. Part of blast holes are grooved to ensure the quality of excavation and observe the protection effect on surrounding rock compared with circular blast holes. Figure8 is the front view and side view of excavation area and monitoring point area of rock-anchored beam blasting. The vertical crustal stress of rock mass in excavation area is 12 MPa and horizontal crustal stress is 21 MPa. The blasting load is directly exerted on the lines of centers of blast holes, so the model of protecting layer of rock-anchored beam only builds the part between lines of centers of auxiliary blast holes to those of main blast holes. Exert vertical crustal stress on the above and right plane of model, vertical restraints on bottom, left and back plane, and nonreflected border condition on each plane besides free surface.Figure 8
The front view and side view of blasting test of rock-anchored beam (unit: m).The rock parameter [20] of rock-anchored beam of underground caverns based on the field survey dates is shown in Table 2.Table 2
Rock mass physical and mechanical material.
ρ
1/(g/cm3)
E
0/GPa
μ
σ
t/MPa
σ
c/MPa
E
t
a
n/GPa
2.7
50.0
0.22
6.0
78.0
8.0
Among them,ρ1 represents the density of rock mass; E0 represents Young’s modulus of rock mass; μ represents Poisson’s ratio; σt represents the tensile strength of rock mass. σc represents the compress strength of rock mass. Etan represents the tangent modulus of rock mass.According to the construction design of the main underground powerhouse, the dynamic finite element program is used to establish the calculation model shown in Figure9 with the reference of Figure 8. The model is 14.5 m in height, 55 m in width, and 5 m in thickness. The blasting load is exerted on lines of centers of auxiliary holes first and then on lines of centers of smooth blasting holes.Figure 9
The front view and side view of calculation model of blasting excavation of rock-anchored beam (unit: m).


### 3.2. The Two-Dimensional Model of Blasting Excavation of Rock-Anchored Beam

Due to the fact that the V-shape grooving has been done on part of smooth blasting holes in the blasting excavation process of rock-anchored beam on site, two-dimensional crack forming model of blast holes is added in order to compare the effect of crack forming of grooved blast holes with that of circular blast holes. In the model, the spacing and array pitch of blast holes are shown in Figure10. The vertical symmetry constraints are exerted on front and back planes; the crustal stress is exerted on left plane; the symmetry constraints are exerted on right plane; the nonreflecting boundary and constraints are exerted on upper plane and the bottom is free surface. The two smooth blasting holes on left plane are grooved blast holes and the two on the right plane are circular blast holes.Figure 10
The calculation model of two-dimensional plane (unit: cm).
(a)
Size diagram of explosive and air
(b)
Size diagram of explosive and rockThe state equation of detonation products applies the JWL equation of state [21] specifically set for describing blasting pressure characteristics of energetic material such as explosive; its expression is as follows:(11)P=A1-wR1Ve-R1V+B1-wR2Ve-R2V+wE2V,where V represents relative volume of detonation gas; A, B, R1, R2, w, and E2 all represent material parameters of the state equation.The specific parameters of rock emulsion explosive in model and the JWL state equation are shown in Table3.Table 3
Explosive and JWL state equation parameters [12].
D/(m/s)
A/Gpa
B/Gpa
R
1
R
2
w
E
2/(J/m)
3200
214
0.18
4.15
0.95
0.15
4.19 × 109
Among them,D represents the velocity of detonation wave.


## 3.1. The Three-Dimensional Model of the Blasting Excavation of Rock-Anchored Beam

Taking the blasting excavation of the protective layers of rock-anchored beams in the underground main power house of the right bank of the dam as an example, the blasting excavation of protective layers of rock-anchored beams drills the vertical smooth blasting holes, oblique smooth blasting holes, and auxiliary blast holes manually in the rate of 20.4 m excavation length each time. Meanwhile, the blasting vibration monitoring is conducted on the sidewalls near excavated rock-anchored beams. The measuring point vibration monitoring of near-field rock during the process of blasting excavation is conducted in order to make sure the disturbances of blasting excavation to surrounding rocks in a safety range. Meanwhile, for ascertaining the safety of monitoring equipment, the nearest measuring points are arranged 10 m away from the boundary of the excavation area, and then one measuring point is arranged every 5 m, and a total of five measuring points are arranged. It is as shown in Figure7.Figure 7
The schematic diagram of measuring point layout of blast vibration velocity (unit: m).The elevation of rock-anchored beams excavated on site ranges from EL. 602.40 m to EL. 606.90 m. Part of blast holes are grooved to ensure the quality of excavation and observe the protection effect on surrounding rock compared with circular blast holes. Figure8 is the front view and side view of excavation area and monitoring point area of rock-anchored beam blasting. The vertical crustal stress of rock mass in excavation area is 12 MPa and horizontal crustal stress is 21 MPa. The blasting load is directly exerted on the lines of centers of blast holes, so the model of protecting layer of rock-anchored beam only builds the part between lines of centers of auxiliary blast holes to those of main blast holes. Exert vertical crustal stress on the above and right plane of model, vertical restraints on bottom, left and back plane, and nonreflected border condition on each plane besides free surface.Figure 8
The front view and side view of blasting test of rock-anchored beam (unit: m).The rock parameter [20] of rock-anchored beam of underground caverns based on the field survey dates is shown in Table 2.Table 2
Rock mass physical and mechanical material.
ρ
1/(g/cm3)
E
0/GPa
μ
σ
t/MPa
σ
c/MPa
E
t
a
n/GPa
2.7
50.0
0.22
6.0
78.0
8.0
Among them,ρ1 represents the density of rock mass; E0 represents Young’s modulus of rock mass; μ represents Poisson’s ratio; σt represents the tensile strength of rock mass. σc represents the compress strength of rock mass. Etan represents the tangent modulus of rock mass.According to the construction design of the main underground powerhouse, the dynamic finite element program is used to establish the calculation model shown in Figure9 with the reference of Figure 8. The model is 14.5 m in height, 55 m in width, and 5 m in thickness. The blasting load is exerted on lines of centers of auxiliary holes first and then on lines of centers of smooth blasting holes.Figure 9
The front view and side view of calculation model of blasting excavation of rock-anchored beam (unit: m).


## 3.2. The Two-Dimensional Model of Blasting Excavation of Rock-Anchored Beam

Due to the fact that the V-shape grooving has been done on part of smooth blasting holes in the blasting excavation process of rock-anchored beam on site, two-dimensional crack forming model of blast holes is added in order to compare the effect of crack forming of grooved blast holes with that of circular blast holes. In the model, the spacing and array pitch of blast holes are shown in Figure10. The vertical symmetry constraints are exerted on front and back planes; the crustal stress is exerted on left plane; the symmetry constraints are exerted on right plane; the nonreflecting boundary and constraints are exerted on upper plane and the bottom is free surface. The two smooth blasting holes on left plane are grooved blast holes and the two on the right plane are circular blast holes.Figure 10
The calculation model of two-dimensional plane (unit: cm).
(a)
Size diagram of explosive and air
(b)
Size diagram of explosive and rockThe state equation of detonation products applies the JWL equation of state [21] specifically set for describing blasting pressure characteristics of energetic material such as explosive; its expression is as follows:(11)P=A1-wR1Ve-R1V+B1-wR2Ve-R2V+wE2V,where V represents relative volume of detonation gas; A, B, R1, R2, w, and E2 all represent material parameters of the state equation.The specific parameters of rock emulsion explosive in model and the JWL state equation are shown in Table3.Table 3
Explosive and JWL state equation parameters [12].
D/(m/s)
A/Gpa
B/Gpa
R
1
R
2
w
E
2/(J/m)
3200
214
0.18
4.15
0.95
0.15
4.19 × 109
Among them,D represents the velocity of detonation wave.


## 4. The Comparative Analysis of Computed Results and Monitoring Result

### 4.1. The Comparative Analysis of Numerical Calculation and Measured Results Blasting Vibration Velocity

Conducting three groups of numerical calculations according to the linear charging density on site, calculations of three groups of equivalent load are based on formulas (1) and (2) and the parameters in Table 1. The results of smooth blasting holes are expressed as (a) 17.9 MPa, (b) 19.3 MPa, and (c) 23.4 MPa while those of auxiliary blasting holes are (a) 14.1 MPa, (b) 15.2 MPa, and (c) 18.5 MPa. The three groups of blasting load are exerted on lines of centers of smooth blasting holes as well as those of auxiliary blast holes, respectively.Through extracting the vibration velocity on measuring points shown in Figure8 through postprocessor and making a comparison with that on same points on site, as shown in Figure 11. It will take a lot of space to compare every measuring points’ wave velocity curves, so here comparison and analysis of only the first 0.5 s of MP1 are done, during which the vibration velocity is large, and the equivalent load of blast hole is 17.9 Mpa–14.1 Mpa. From Figures 11(a), 11(b), and 11(c) it can be seen that there will be a peak vibration velocity of measuring points after auxiliary blast holes detonate first and a bigger peak value after smooth blast holes detonate. Then the comparative analysis of peak vibration velocity on measuring points by numerical simulation and site measurement is conducted in three directions. It is shown on the Figures 11(d), 11(e), and 11(f) that the vibration variation trends in different directions are similar on conditions of the same linear charging density. Furthermore, the vibration velocity on measuring points of numerical simulation approximates that of site measuring. With the increase of linear charging density, vibration velocity of same measuring points will increase correspondingly. Moreover, because the increased value of the vibration velocity on measuring points away from blasting source is relatively small and the vibration velocity of points in a rather near distance is much more affected by linear charging density, in the process of blasting excavation, whether the explosive load of blast holes corresponds the safety standards mainly depends on the values of vibration rate of rock mass, which is close to the blasting source. The vibration velocity of the measuring point at 10 m away from the blasting sources should be less than 10 cm/s [22]. We could conclude from the figure that charging and layout of blast holes on site accord with the regulation.Figure 11
The comparison of vibration velocity of measuring points in numerical simulation and field measurement.
(a)
MP1 vibration curve in theX direction
(b)
MP1 vibration curve in theY direction
(c)
MP1 vibration curve in theZ direction
(d)
Equivalent load 17.9 MPa–14.1 MPa
(e)
Equivalent load 19.3 MPa–15.2 MPa
(f)
Equivalent load 23.4 MPa–18.5 MPaSelecting one group unit from the rock-anchored beam region of rock mass shown in Figure5, exerting the average crustal stress suffered from three directions is exerted on two-dimensional planform as the original value of crustal stress. As is shown in Table 4, the average unit stress in X direction which is perpendicular to tunnel axis is 0.28 MPa; the average unit stress in Y direction is 20.32 MPa, and the average unit stress in Z direction which parallels tunnel axis is 20.32 MPa. Due to the fact that two-dimensional plane graph is the XOZ planform of three-dimensional diagram, there is no need for the exertion of crustal stress on Y direction. While the stress on X direction is incredible small compared to that on Z direction, therefore, it could be ignored for the convenience of model calculation, only exerting the stress of Z direction on two-dimensional plane graph.Table 4
Stress values of five measuring points of rock-anchored beam.
Direction
Measuring points
Average value
1#
2#
3#
4#
5#
X/MPa
0.45
0.33
0.25
0.20
0.17
0.28
Y/MPa
11.55
15.92
13.81
13.28
13.22
13.56
Z/MPa
20.42
20.46
20.08
20.22
20.43
20.32


### 4.2. The Comparison between Numerical Simulated Crack Forming Figure and Practical Situation on Site

A comparison is made among the two-dimensional plane figure of cracks and the figures of contour surface and rock mass crushing on site in three values of different linear charging density. As shown in Figure12, when, in the same linear charging density, the modeling figure of crack forming consists with the excavation figure on site.(1)
When the linear charging densities of the smooth lasting holes and auxiliary holes are 65 g/m and 206 g/m, respectively, it is shown in the figure that the cracks between grooved smooth blasting holes are quite flat and the damage to surrounding rock is rather small, while there are no penetrating cracks formed between circular smooth blasting holes, which indicates that such charging is too small to form cracks between circular blast holes. However, the grooved blast holes take effect on energy gathering and guiding so as to utilize the more blasting power on crack forming of rock between blast holes. Furthermore, the crushing degree of excavated rock mass is rather small and the boulder frequency is large, so the demand of crushing degree could not be reached.(2)
When the linear charging densities of smooth blast holes and auxiliary holes are 70 g/m and 206 g/m, respectively, penetrating cracks could be formed between blast holes after the blasting. However, cracks between grooved blast holes are more flat than those of circular blast holes, which indicate that penetrating cracks could be formed among blast holes in different shapes under such charging. Moreover, the crushing degree of excavated rock mass is large so it is convenient for the shipment of crushing rock mass.(3)
When the linear charging densities of smooth blast holes and auxiliary holes are 85 g/m and 250 g/m, respectively, besides the cracks forming among smooth blast holes of different shapes, serious damage of surrounding rocks also takes place and the crushing degree of rock mass is rather large.Figure 12
The comparison of excavated contour and crushing rate of rock mass in numerical simulation and field measurement under different charge amount.
(a)
The linear charging density is 65 g/m
(b)
The linear charging density is 70 g/m
(c)
The linear charging density is 85 g/mAccording to the contour surface of cracks forming on site, the damage reaches the largest extent where vertical smooth blast holes intersect with the oblique ones. The fact consists with the research conclusion of Dong et al. [23]. It indicates that the charging amount of blast holes must be in fine control in order to make sure the oriented cracks forming of grooved blast holes could take place, ensuring the flat and smooth of cracks forming among blast holes and meanwhile reducing the charge amount as much as possible.Based on the aforesaid analysis, when the linear charging density is 70 g/m, the penetrating cracks could be formed among smooth blasting holes and the damage to surrounding rocks is in low extent. Furthermore, when grooved blast hole is applied, the penetrating cracks could be much easier to form among blast holes and the damage to surrounding rocks is rather small.


### 4.3. The Analysis of Blast Vibration and Crack Forming under Different Crustal Stress

According to the analysis of Sections4.1 and 4.2, when linear charging density is 70 g/m, whether in numerical simulation or on-site test, the blast vibration velocity in near-field of blasting source could be assured to meet the safety requirements. Moreover, the cracks forming among smooth blasting holes could be guaranteed and the damage to surrounding rock is lowered in certain extent when the density is in such value. Supposing the value of linear charging density is 70 g/m, the proportion of the horizontal crustal stress and vertical crustal stress is setting 1 : 1 in value of 0 MPa, 5 MPa, 10 MPa, 20 MPa, 40 MPa, and 80 MPa, respectively. Taking into account the blasting vibration velocity of the same measuring points under different crustal stress and the cracks forming between the smooth blasting holes, the blast vibration velocity of measuring points could be worked out by numerical calculation as shown in Figure 13.Figure 13
The comparison of vibration velocity in the same measuring point under different crustal stress.
(a)
Perpendicular to tunnel axis directionX
(b)
Upright directionY
(c)
Parallels to tunnel axis directionZAccording to the vibration velocity in three directions shown in figures, the vibration velocity of measuring points in directionX perpendicular to tunnel axis is in the largest value while that in upright direction Y is the smallest. When the crustal stress is 0–40 MPa, the changing rule of vibration velocity of the same measuring point in the same directions consists with each other, which indicates that there is no apparent relation between vibration velocity of surrounding rock in near-field of blasting and magnitude of crustal stress under the identical charging amount and charging structures. When the crustal stress is 80 MPa, the vibration velocities of measuring points in three directions all increase, which indicates that when the crustal stress reaches a certain value, a large amount of potential energy in underground rock mass is released. When explosive loads exert force disturbances on rock mass, larger vibration has been produced under the effect of blasting kinetic energy and crustal stress potential energy on the rock mass nearby blasting source.According to the calculation and analysis in Section4.2, only exert stress in Z direction parallel to tunnel axis on two-dimensional plane model. Thus, the stress of Z direction on the same measuring point is selected to exert on two-dimensional plane model as the original crustal stress. Table 5 shows the value of stress in Z direction of same measuring points but different crustal stress.Table 5
Stress values ofZ direction of the same measuring points under different crustal stress.
Crustal stress/MPa
Measuring points
Average value/MPa
1#
2#
3#
4#
5#
0
0
0
0
0
0
0
5
4.65
4.57
4.60
4.64
4.64
4.62
10
9.30
9.13
9.18
9.28
9.28
9.23
20
18.60
18.26
18.38
18.56
18.55
18.47
40
37.18
36.50
36.73
37.09
37.07
36.91
80
74.90
73.70
75.20
76.73
77.80
75.67Keeping the charging of blast holes in a certain amount and exerting the average value of stress inZ direction that is worked out in Table 5 on two-dimensional plane model, the crack forming figure is worked out under different stresses as shown in Figure 14. The damage depth and area of surrounding rocks and the crushing rate of excavated rock mass are summarized as shown in Table 6.Table 6
The damage depth, area of surrounding rocks, and the crushing rate of excavated rock mass.
Crustal stress/MPa
0
5
10
20
40
80
Damage depth of surrounding rocks/cm
>50
>50
≈20
≈10
<10
≈10
Damage area of surrounding rocks/%
>50
>50
≈40
≈20
<20
≈6.7
Crushing rate of excavated rock mass/%
>75
>70
>70
>50
<30
<10Figure 14
The crack forming diagram of rock mass under different crustal stress.
(a)
0 MPa
(b)
5 MPa
(c)
10 MPa
(d)
20 MPa
(e)
40 MPa
(f)
80 MPaAs shown in Figure14, the smaller the crustal stress is, the larger the crushing degree of rock mass and the damage to surrounding rocks are. When crustal stress increases and although the damage to surrounding rocks reduces, the crushing degree of rock mass lowers as well, so the required crushing rate of excavated rock mass on site could not be reached. With such charging structure and charging amount, larger crushing degree and low damage can be achieved, only when the exerted crustal stress is 20 MPa, known as the most approximate value to crustal stress on site, and the effect of larger crushing degree of rock mass and the lower damage to surrounding rocks near blast holes can be achieved, which indicates that charging amount must be adjusted correspondingly according to different crustal stress when the underground cavern is excavated by blasting.According to the reference of Figure13 and Table 6, when the crustal stress is 0–40 MPa, the charging amount of blast holes should be increased followed by the increase of crustal stress so as to meet the requirement of blasting excavation. When the crustal stress continues to increase, the vibration velocity of same measuring points on moderate and far area of blasting under the same charging amount will be increased, which could need the reduction of charge amount. However, the lower crushing degree of rock mass indicates the increase of charge amount. Due to the above analysis, in the blasting excavation conducted in underground caverns where the crustal stress is rather large, it is not simply the matter of increase or reduction of the charge amount but rather the new layout of blast holes as well as the adjustment of charge portion of smooth blasting holes and auxiliary holes.


## 4.1. The Comparative Analysis of Numerical Calculation and Measured Results Blasting Vibration Velocity

Conducting three groups of numerical calculations according to the linear charging density on site, calculations of three groups of equivalent load are based on formulas (1) and (2) and the parameters in Table 1. The results of smooth blasting holes are expressed as (a) 17.9 MPa, (b) 19.3 MPa, and (c) 23.4 MPa while those of auxiliary blasting holes are (a) 14.1 MPa, (b) 15.2 MPa, and (c) 18.5 MPa. The three groups of blasting load are exerted on lines of centers of smooth blasting holes as well as those of auxiliary blast holes, respectively.Through extracting the vibration velocity on measuring points shown in Figure8 through postprocessor and making a comparison with that on same points on site, as shown in Figure 11. It will take a lot of space to compare every measuring points’ wave velocity curves, so here comparison and analysis of only the first 0.5 s of MP1 are done, during which the vibration velocity is large, and the equivalent load of blast hole is 17.9 Mpa–14.1 Mpa. From Figures 11(a), 11(b), and 11(c) it can be seen that there will be a peak vibration velocity of measuring points after auxiliary blast holes detonate first and a bigger peak value after smooth blast holes detonate. Then the comparative analysis of peak vibration velocity on measuring points by numerical simulation and site measurement is conducted in three directions. It is shown on the Figures 11(d), 11(e), and 11(f) that the vibration variation trends in different directions are similar on conditions of the same linear charging density. Furthermore, the vibration velocity on measuring points of numerical simulation approximates that of site measuring. With the increase of linear charging density, vibration velocity of same measuring points will increase correspondingly. Moreover, because the increased value of the vibration velocity on measuring points away from blasting source is relatively small and the vibration velocity of points in a rather near distance is much more affected by linear charging density, in the process of blasting excavation, whether the explosive load of blast holes corresponds the safety standards mainly depends on the values of vibration rate of rock mass, which is close to the blasting source. The vibration velocity of the measuring point at 10 m away from the blasting sources should be less than 10 cm/s [22]. We could conclude from the figure that charging and layout of blast holes on site accord with the regulation.Figure 11
The comparison of vibration velocity of measuring points in numerical simulation and field measurement.
(a)
MP1 vibration curve in theX direction
(b)
MP1 vibration curve in theY direction
(c)
MP1 vibration curve in theZ direction
(d)
Equivalent load 17.9 MPa–14.1 MPa
(e)
Equivalent load 19.3 MPa–15.2 MPa
(f)
Equivalent load 23.4 MPa–18.5 MPaSelecting one group unit from the rock-anchored beam region of rock mass shown in Figure5, exerting the average crustal stress suffered from three directions is exerted on two-dimensional planform as the original value of crustal stress. As is shown in Table 4, the average unit stress in X direction which is perpendicular to tunnel axis is 0.28 MPa; the average unit stress in Y direction is 20.32 MPa, and the average unit stress in Z direction which parallels tunnel axis is 20.32 MPa. Due to the fact that two-dimensional plane graph is the XOZ planform of three-dimensional diagram, there is no need for the exertion of crustal stress on Y direction. While the stress on X direction is incredible small compared to that on Z direction, therefore, it could be ignored for the convenience of model calculation, only exerting the stress of Z direction on two-dimensional plane graph.Table 4
Stress values of five measuring points of rock-anchored beam.
Direction
Measuring points
Average value
1#
2#
3#
4#
5#
X/MPa
0.45
0.33
0.25
0.20
0.17
0.28
Y/MPa
11.55
15.92
13.81
13.28
13.22
13.56
Z/MPa
20.42
20.46
20.08
20.22
20.43
20.32


## 4.2. The Comparison between Numerical Simulated Crack Forming Figure and Practical Situation on Site

A comparison is made among the two-dimensional plane figure of cracks and the figures of contour surface and rock mass crushing on site in three values of different linear charging density. As shown in Figure12, when, in the same linear charging density, the modeling figure of crack forming consists with the excavation figure on site.(1)
When the linear charging densities of the smooth lasting holes and auxiliary holes are 65 g/m and 206 g/m, respectively, it is shown in the figure that the cracks between grooved smooth blasting holes are quite flat and the damage to surrounding rock is rather small, while there are no penetrating cracks formed between circular smooth blasting holes, which indicates that such charging is too small to form cracks between circular blast holes. However, the grooved blast holes take effect on energy gathering and guiding so as to utilize the more blasting power on crack forming of rock between blast holes. Furthermore, the crushing degree of excavated rock mass is rather small and the boulder frequency is large, so the demand of crushing degree could not be reached.(2)
When the linear charging densities of smooth blast holes and auxiliary holes are 70 g/m and 206 g/m, respectively, penetrating cracks could be formed between blast holes after the blasting. However, cracks between grooved blast holes are more flat than those of circular blast holes, which indicate that penetrating cracks could be formed among blast holes in different shapes under such charging. Moreover, the crushing degree of excavated rock mass is large so it is convenient for the shipment of crushing rock mass.(3)
When the linear charging densities of smooth blast holes and auxiliary holes are 85 g/m and 250 g/m, respectively, besides the cracks forming among smooth blast holes of different shapes, serious damage of surrounding rocks also takes place and the crushing degree of rock mass is rather large.Figure 12
The comparison of excavated contour and crushing rate of rock mass in numerical simulation and field measurement under different charge amount.
(a)
The linear charging density is 65 g/m
(b)
The linear charging density is 70 g/m
(c)
The linear charging density is 85 g/mAccording to the contour surface of cracks forming on site, the damage reaches the largest extent where vertical smooth blast holes intersect with the oblique ones. The fact consists with the research conclusion of Dong et al. [23]. It indicates that the charging amount of blast holes must be in fine control in order to make sure the oriented cracks forming of grooved blast holes could take place, ensuring the flat and smooth of cracks forming among blast holes and meanwhile reducing the charge amount as much as possible.Based on the aforesaid analysis, when the linear charging density is 70 g/m, the penetrating cracks could be formed among smooth blasting holes and the damage to surrounding rocks is in low extent. Furthermore, when grooved blast hole is applied, the penetrating cracks could be much easier to form among blast holes and the damage to surrounding rocks is rather small.


## 4.3. The Analysis of Blast Vibration and Crack Forming under Different Crustal Stress

According to the analysis of Sections4.1 and 4.2, when linear charging density is 70 g/m, whether in numerical simulation or on-site test, the blast vibration velocity in near-field of blasting source could be assured to meet the safety requirements. Moreover, the cracks forming among smooth blasting holes could be guaranteed and the damage to surrounding rock is lowered in certain extent when the density is in such value. Supposing the value of linear charging density is 70 g/m, the proportion of the horizontal crustal stress and vertical crustal stress is setting 1 : 1 in value of 0 MPa, 5 MPa, 10 MPa, 20 MPa, 40 MPa, and 80 MPa, respectively. Taking into account the blasting vibration velocity of the same measuring points under different crustal stress and the cracks forming between the smooth blasting holes, the blast vibration velocity of measuring points could be worked out by numerical calculation as shown in Figure 13.Figure 13
The comparison of vibration velocity in the same measuring point under different crustal stress.
(a)
Perpendicular to tunnel axis directionX
(b)
Upright directionY
(c)
Parallels to tunnel axis directionZAccording to the vibration velocity in three directions shown in figures, the vibration velocity of measuring points in directionX perpendicular to tunnel axis is in the largest value while that in upright direction Y is the smallest. When the crustal stress is 0–40 MPa, the changing rule of vibration velocity of the same measuring point in the same directions consists with each other, which indicates that there is no apparent relation between vibration velocity of surrounding rock in near-field of blasting and magnitude of crustal stress under the identical charging amount and charging structures. When the crustal stress is 80 MPa, the vibration velocities of measuring points in three directions all increase, which indicates that when the crustal stress reaches a certain value, a large amount of potential energy in underground rock mass is released. When explosive loads exert force disturbances on rock mass, larger vibration has been produced under the effect of blasting kinetic energy and crustal stress potential energy on the rock mass nearby blasting source.According to the calculation and analysis in Section4.2, only exert stress in Z direction parallel to tunnel axis on two-dimensional plane model. Thus, the stress of Z direction on the same measuring point is selected to exert on two-dimensional plane model as the original crustal stress. Table 5 shows the value of stress in Z direction of same measuring points but different crustal stress.Table 5
Stress values ofZ direction of the same measuring points under different crustal stress.
Crustal stress/MPa
Measuring points
Average value/MPa
1#
2#
3#
4#
5#
0
0
0
0
0
0
0
5
4.65
4.57
4.60
4.64
4.64
4.62
10
9.30
9.13
9.18
9.28
9.28
9.23
20
18.60
18.26
18.38
18.56
18.55
18.47
40
37.18
36.50
36.73
37.09
37.07
36.91
80
74.90
73.70
75.20
76.73
77.80
75.67Keeping the charging of blast holes in a certain amount and exerting the average value of stress inZ direction that is worked out in Table 5 on two-dimensional plane model, the crack forming figure is worked out under different stresses as shown in Figure 14. The damage depth and area of surrounding rocks and the crushing rate of excavated rock mass are summarized as shown in Table 6.Table 6
The damage depth, area of surrounding rocks, and the crushing rate of excavated rock mass.
Crustal stress/MPa
0
5
10
20
40
80
Damage depth of surrounding rocks/cm
>50
>50
≈20
≈10
<10
≈10
Damage area of surrounding rocks/%
>50
>50
≈40
≈20
<20
≈6.7
Crushing rate of excavated rock mass/%
>75
>70
>70
>50
<30
<10Figure 14
The crack forming diagram of rock mass under different crustal stress.
(a)
0 MPa
(b)
5 MPa
(c)
10 MPa
(d)
20 MPa
(e)
40 MPa
(f)
80 MPaAs shown in Figure14, the smaller the crustal stress is, the larger the crushing degree of rock mass and the damage to surrounding rocks are. When crustal stress increases and although the damage to surrounding rocks reduces, the crushing degree of rock mass lowers as well, so the required crushing rate of excavated rock mass on site could not be reached. With such charging structure and charging amount, larger crushing degree and low damage can be achieved, only when the exerted crustal stress is 20 MPa, known as the most approximate value to crustal stress on site, and the effect of larger crushing degree of rock mass and the lower damage to surrounding rocks near blast holes can be achieved, which indicates that charging amount must be adjusted correspondingly according to different crustal stress when the underground cavern is excavated by blasting.According to the reference of Figure13 and Table 6, when the crustal stress is 0–40 MPa, the charging amount of blast holes should be increased followed by the increase of crustal stress so as to meet the requirement of blasting excavation. When the crustal stress continues to increase, the vibration velocity of same measuring points on moderate and far area of blasting under the same charging amount will be increased, which could need the reduction of charge amount. However, the lower crushing degree of rock mass indicates the increase of charge amount. Due to the above analysis, in the blasting excavation conducted in underground caverns where the crustal stress is rather large, it is not simply the matter of increase or reduction of the charge amount but rather the new layout of blast holes as well as the adjustment of charge portion of smooth blasting holes and auxiliary holes.


## 5. Conclusion

The following conclusions could be drawn according to the analytical comparison of dynamic finite element calculation and on-site test:(1)
The vibration velocity test of three directions of surrounding rocks in moderate and far area blasting excavation of rock-anchored beams found out that the vibration velocity of rock mass perpendicular to the blasting excavation surface is the largest which should be used as the reference of vibration velocity control in rock mass excavation.(2)
The numerical simulation and on-site test indicate that when grooved and circular blast holes blast under the same charge amount, the grooved one is proved to carry out fine energy gathering effect as well as reduce the damage to surrounding rocks on a certain extent.(3)
Research on blasting excavation of underground caverns under different crustal stress found that, in identical charge amount, the crushing degree of rock mass and the damage to surrounding rocks decreases with the increase of crustal stress. When crustal stress is 0–40 MPa, the blasting vibration velocities of moderate and far area of rock mass approximate and when it reaches 80 MPa, the blasting vibration velocity on the same measuring point increases, so adjustments of blast holes layout and charge amount need to be done based on different crustal stress.Due to the absence of practical situation of fine blasting project of rock-anchored beam excavation on underground caverns in high crustal stress, the results of numerical simulation could not be proved. Experiments and researches on such problem will be continued.



---

*Source: 1812080-2017-04-23.xml*