Research Article  Open Access
XinPing Li, JunHong Huang, Yi Luo, Qian Dong, YouHua Li, Yong Wan, TingTing Liu, "Numerical Simulation of Blast Vibration and Crack Forming Effect of RockAnchored Beam Excavation in Deep Underground Caverns", Shock and Vibration, vol. 2017, Article ID 1812080, 13 pages, 2017. https://doi.org/10.1155/2017/1812080
Numerical Simulation of Blast Vibration and Crack Forming Effect of RockAnchored Beam Excavation in Deep Underground Caverns
Abstract
Aiming at surrounding rock damage induced by dynamic disturbance from blasting excavation of rockanchored 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 nearfield 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.
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 largescale and ultralarge type appear along with the enlargement of hydropower in scale. The rockanchored 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 rockanchored 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 rockanchored beam since it can take full advantage of the bearing capacity of surrounding walls. The application of rockanchored 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 rockanchored beam is widespread in many largescale overseas hydropower stations [3].
However, with considerable difficulty, the excavation of rockanchored 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 fluidstructure 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 fluidstructure 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 fluidstructure 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 nearfield of blasting source. There is still a way to go for the simulation of the whole process from the blasting of nearfield of blasting source to the blasting vibration spreading of explosion of farfield [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 SaintVenant principle, the latter method is more reasonable and effective, which adopts triangle pulse waves as equivalent blasting loads to conduct the threedimensional 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 threedimensional models in the calculation of vibration speed of the blasting of rock mass in the moderate and far distance and twodimensional models in the simulation of crushing conditions of rock mass in nearfield 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 nearfield 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/LSDYNA, 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 LSDYNA for dynamic calculation is motion equation:where is global mass matrix, is damping matrix, is global stiffness matrix, , , and are nodes’ acceleration vector, velocity vector, and displacement vector, respectively, and 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 time should satisfy the displacement vector of (1) to only require time discrete points to satisfy the above equation of motion.
The solving equation used for static calculation is as follows:where is global stiffness matrix, is displacement field of rock mass, and 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 LSDYNA.
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 threedimensional 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.
Overseas scholars have come up with CJ 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 CJ conditions. The simplified load peak value exerting on blast holes iswhere represents the density of explosive (g/cm^{3}); represents explosive velocity (m/s); represents the isentropic exponent of explosive, usually takes 3; represents the charge diameter; represents borehole diameter; represents the length of blast holes; represents the length of charge column; here smooth blasting holes take and auxiliary holes take ; represents pressure increasing factor when the detonation gas expansion collides with hole wall, = 8–11, and it takes 10 for emulsion explosive.
As shown in Figure 2, the equivalent peak pressure calculated through force balance principle iswhere represents blast hole radius; represents distances between 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:where is the von Mises effective stress of any point in rock mass; is the stress components of rock mass; is the tensile stress of explosion load of any point in rock mass; 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 Figure 3. 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.
As shown in Figure 5, during the detonation of explosives in borhole, tensile stress would be concentrated at groove tip instantaneously, and crack initial from element must satisfy the formula
The dynamic compressive stress of rock increases with the improvement of loaded strain rate, generally approximated by the following equation [16]:where in the equation refers to the uniaxial static compressive stress of rock.
The loaded strain rate of rock in blasting during programs is within 10^{0}–10^{5} s^{−1}, among which the strain rate in crushed zones could be and in cracked zones could be .
For the lack of corresponding analytical data of experiments and theories, the value of dynamic tensile strength approximateswhere in the equation refers to the uniaxial static tensile strength of rock mass.
The paper applies the fluidstructure interaction algorithm on the calculation of explosive, rock, and air through twodimensional 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 nearfield, 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 CowperSymonds 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:In the formula, is strain rate and and are relative parameters in CowperSymonds, which are generally determined by features of material strain rate. The paper refers to the conclusion in reference [17], which indicates that is 2.5 s^{−1}, is 4, represents initial yield stress of rock, and is the hardening coefficient. When , represents plastic kinematic hardening; when , it represents isotropic hardening; when , it represents mixed hardening. The paper adopts the plastic following hardening in calculation; therefore, is taken. , , and represent plastic strainhardening modulus, Young’s modulus, and tangent modulus separately.
When the material is in hardening stage, the efficient plastic strain rate corresponded could be expressed as follows:where and represent plastic strain rate and elastic strain rate separately.
3. The Numerical Model of Blasting Excavation of RockAnchored 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 rockanchored beam of rock mass in underground plants. Thus, there is strong need for the fine control of blasting excavation of the rockanchored 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 Table 1. 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.

3.1. The ThreeDimensional Model of the Blasting Excavation of RockAnchored Beam
Taking the blasting excavation of the protective layers of rockanchored beams in the underground main power house of the right bank of the dam as an example, the blasting excavation of protective layers of rockanchored 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 rockanchored beams. The measuring point vibration monitoring of nearfield 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 Figure 7.
The elevation of rockanchored 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. Figure 8 is the front view and side view of excavation area and monitoring point area of rockanchored 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 rockanchored 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.
The rock parameter [20] of rockanchored beam of underground caverns based on the field survey dates is shown in Table 2.
 
Among them, represents the density of rock mass; represents Young’s modulus of rock mass; represents Poisson’s ratio; represents the tensile strength of rock mass. represents the compress strength of rock mass. 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 Figure 9 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.
3.2. The TwoDimensional Model of Blasting Excavation of RockAnchored Beam
Due to the fact that the Vshape grooving has been done on part of smooth blasting holes in the blasting excavation process of rockanchored beam on site, twodimensional 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 Figure 10. 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.
(a) Size diagram of explosive and air
(b) Size diagram of explosive and rock
The 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:where represents relative volume of detonation gas; , , , , , and 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 Table 3.
 
Among them, 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 Figure 8 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.
(a) MP1 vibration curve in the direction
(b) MP1 vibration curve in the direction
(c) MP1 vibration curve in the 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 MPa
Selecting one group unit from the rockanchored beam region of rock mass shown in Figure 5, exerting the average crustal stress suffered from three directions is exerted on twodimensional planform as the original value of crustal stress. As is shown in Table 4, the average unit stress in direction which is perpendicular to tunnel axis is 0.28 MPa; the average unit stress in direction is 20.32 MPa, and the average unit stress in direction which parallels tunnel axis is 20.32 MPa. Due to the fact that twodimensional plane graph is the planform of threedimensional diagram, there is no need for the exertion of crustal stress on direction. While the stress on direction is incredible small compared to that on direction, therefore, it could be ignored for the convenience of model calculation, only exerting the stress of direction on twodimensional plane graph.

4.2. The Comparison between Numerical Simulated Crack Forming Figure and Practical Situation on Site
A comparison is made among the twodimensional 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 Figure 12, 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.
(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/m
According 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 Sections 4.1 and 4.2, when linear charging density is 70 g/m, whether in numerical simulation or onsite test, the blast vibration velocity in nearfield 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.
(a) Perpendicular to tunnel axis direction
(b) Upright direction
(c) Parallels to tunnel axis direction
According to the vibration velocity in three directions shown in figures, the vibration velocity of measuring points in direction perpendicular to tunnel axis is in the largest value while that in upright direction 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 nearfield 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 Section 4.2, only exert stress in direction parallel to tunnel axis on twodimensional plane model. Thus, the stress of direction on the same measuring point is selected to exert on twodimensional plane model as the original crustal stress. Table 5 shows the value of stress in direction of same measuring points but different crustal stress.

Keeping the charging of blast holes in a certain amount and exerting the average value of stress in direction that is worked out in Table 5 on twodimensional 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.

(a) 0 MPa
(b) 5 MPa
(c) 10 MPa
(d) 20 MPa
(e) 40 MPa
(f) 80 MPa
As shown in Figure 14, 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 Figure 13 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 onsite test:(1)The vibration velocity test of three directions of surrounding rocks in moderate and far area blasting excavation of rockanchored 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 onsite 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 rockanchored 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.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work is supported by the National Natural Science Foundation of China (51274157 and 51309183), PLA University of Science and Technology (no. DPMEIKF201410), the Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences (Grant no. Z015005).
References
 R. Lien, J. Kristiansen, and L. S. Pran, “Suspended support for crane rails in Norwegian underground powerhouses,” International Water Power and Dam Construction, vol. 42, no. 4, pp. 48–49, 1990. View at: Google Scholar
 P. M. Johansen, Model Field Tests on Crane Rail Beams at Skjomen and Kvilldal Power Stations, Norwegian Geotechnical Institute, Oslo, Norway, 1984.
 R. D. Gill and D. S. Lawrenson, “Underground works at Kenya's kiambere project,” International Water Power and Dam Construction, vol. 38, no. 5, pp. 13–16, 1986. View at: Google Scholar
 N. Pelekis, B. Theodoulidis, I. Kopanakis, and Y. Theodoridis, “Literature review of spatiotemporal database models,” Knowledge Engineering Review, vol. 19, no. 3, pp. 235–274, 2004. View at: Publisher Site  Google Scholar
 D. Shi, Y. Li, and S. Zhang, The Explicit Dynamic Analysis on the Basis of ANSYS/LSDYNA 8.1, Tsinghua University Press, Beijing, China, 2004 (Chinese).
 B. Jinze, Theory Basis and Examples Analysis of LSDYNA 3D, Science Press, Beijing, China, 2005 (Chinese).
 G. Ling, Y. Chen, and W. Huang, “A classification method of spatiotemporal data model,” in Proceedings of the Nationwide Cartography and GIS Academic Conferences, pp. 734–740, 2004. View at: Google Scholar
 G. W. Ma, H. Hao, and Y. X. Zhou, “Modeling of wave propagation induced by underground explosion,” Computers and Geotechnics, vol. 22, no. 34, pp. 283–303, 1998. View at: Publisher Site  Google Scholar
 C. Wu, Y. Lu, and H. Hao, “Numerical prediction of blastinduced stress wave from largescale underground explosion,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 28, no. 1, pp. 93–109, 2004. View at: Publisher Site  Google Scholar
 J. Toraño, R. Rodríguez, I. Diego, J. M. Rivas, and M. D. Casal, “FEM models including randomness and its application to the blasting vibrations prediction,” Computers and Geotechnics, vol. 33, no. 1, pp. 15–28, 2006. View at: Publisher Site  Google Scholar
 H. Ding and Z. Zheng, “Source model for blasting vibration,” Science in China, Series E: Technological Sciences, vol. 45, no. 4, pp. 395–407, 2002. View at: Google Scholar
 D. S. Preece and B. J. Thorne, “A study of detonation timing and fragmentation using 3D finite element techniques and a damage constitutive model,” in Rock Fragmentation by Blasting, B. Mohanty, Ed., pp. 147–156, Balkema, Rotterdam, The Netherlands, 1996. View at: Google Scholar
 W. Lu, J. Yang, P. Yan et al., “Dynamic response of rock mass induced by the transient release of insitu stress,” International Journal of Rock Mechanics and Mining Sciences, vol. 53, pp. 129–141, 2012. View at: Publisher Site  Google Scholar
 L. X. Xie, W. B. Lu, Q. B. Zhang, Q. H. Jiang, G. H. Wang, and J. Zhao, “Damage evolution mechanisms of rock in deep tunnels induced by cut blasting,” Tunnelling and Underground Space Technology, vol. 58, pp. 257–270, 2016. View at: Publisher Site  Google Scholar
 H. Xu, W. Lu, and X. Zhou, “An equivalent approach for acting blasting load in dynamic finite element simulation of blasting vibration,” Engineering Journal of Wuhan University, vol. 41, no. 1, pp. 67–71, 2008 (Chinese). View at: Google Scholar
 X. Xia, H.B. Li, J.R. Li, K.Q. Xiao, H. Tang, and G.W. Ma, “Numerical simulation of blastinduced cracks in rock,” Rock and Soil Mechanics, vol. 27, no. 11, pp. 1987–1991, 2006. View at: Google Scholar
 R. Yang, W. F. Bawden, and P. D. Katsabanis, “A new constitutive model for blast damage,” International Journal of Rock Mechanics and Mining Sciences and Geomechanics, vol. 33, no. 3, pp. 245–254, 1996. View at: Publisher Site  Google Scholar
 M. Khandelwal and T. N. Singh, “Prediction of blastinduced ground vibration using artificial neural network,” International Journal of Rock Mechanics and Mining Sciences, vol. 46, no. 7, pp. 1214–1222, 2009. View at: Publisher Site  Google Scholar
 D. Johansson and F. Ouchterlony, “Shock wave interactions in rock blasting: the use of short delays to improve fragmentation in modelscale,” Rock Mechanics and Rock Engineering, vol. 46, no. 1, pp. 1–18, 2013. View at: Publisher Site  Google Scholar
 The Construction Bidding Documents with Civil Engineering and Metal Structure Installation Project on the Right Bank Diversion Generating System of Jinsha River Baihetan Hydropower Station, China's Hydropower Consulting Group East China Survey Design and Research Institute, 2013 (Chinese).
 E. L. Lee, H. C. Hornig, and J. W. Kury, “Lawrence Livermore laboratory report,” Tech. Rep. UCRL50442, 1986. View at: Google Scholar
 W.B. Lu, Y. Luo, M. Chen, and D.Q. Shu, “An introduction to Chinese safety regulations for blasting vibration,” Environmental Earth Sciences, vol. 67, no. 7, pp. 1951–1959, 2012. View at: Publisher Site  Google Scholar
 Q. Dong, X. Li, and Y. Luo, “Numerical simulation of damage characteristics of underground cavern smooth blasting excavation based on the improved KUS model,” Electronic Journal of Geotechnical Engineering, vol. 21, no. 19, pp. 6381–6392, 2016, http://ejge.com. View at: Google Scholar
Copyright
Copyright © 2017 XinPing Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.