Hindawi Publishing Corporation e Scientific World Journal Volume 2014, Article ID 379326, 8 pages http://dx.doi.org/10.1155/2014/379326
Research Article Numerical Simulation on Zonal Disintegration in Deep Surrounding Rock Mass Xuguang Chen,1,2 Yuan Wang,1 Yu Mei,1 and Xin Zhang1 1
Institute of Tunnel and Urban Railway Engineering, Key Laboratory of Ministry of Education for Geomechanics and Embankment Engineering, Hohai University, Nanjing 210098, China 2 State Key Laboratory for GeoMechanics and Deep Underground Engineering, China University of Mining & Technology, Xuzhou 221000, China Correspondence should be addressed to Xuguang Chen;
[email protected] Received 13 October 2013; Accepted 21 November 2013; Published 30 January 2014 Academic Editors: B. Lin and B. Uy Copyright © 2014 Xuguang Chen 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. Zonal disintegration have been discovered in many underground tunnels with the increasing of embedded depth. The formation mechanism of such phenomenon is difficult to explain under the framework of traditional rock mechanics, and the fractured shape and forming conditions are unclear. The numerical simulation was carried out to research the generating condition and forming process of zonal disintegration. Via comparing the results with the geomechanical model test, the zonal disintegration phenomenon was confirmed and its mechanism is revealed. It is found to be the result of circular fracture which develops within surrounding rock mass under the high geostress. The fractured shape of zonal disintegration was determined, and the radii of the fractured zones were found to fulfill the relationship of geometric progression. The numerical results were in accordance with the model test findings. The mechanism of the zonal disintegration was revealed by theoretical analysis based on fracture mechanics. The fractured zones are reportedly circular and concentric to the cavern. Each fracture zone ruptured at the elastic-plastic boundary of the surrounding rocks and then coalesced into the circular form. The geometric progression ratio was found to be related to the mechanical parameters and the ground stress of the surrounding rocks.
1. Introduction Recently, many countries have begun to focus on deep resource exploitation. With an increase in embedded depth, the zonal disintegration phenomenon occurs during tunnel excavation. Zonal disintegration refers to “alternating regions of fractured and relatively intact rock masses appearing around or in front of the working stope during the excavation of tunnels in the deep rock mass” [1]. It is a special geological phenomenon that is different from what is observed for shallow embedded tunnels and presents a great hazard to the stability of deep surrounding rocks [2]. Adams and Jager (1980) first observed this phenomenon using a borehole periscope at an embedded depth of 2000–3000 m in the Witwatersrand gold mine, South Africa [3]. Shemyakin et al. [4– 6] explored zonal disintegration by using a resistivity meter in the Taimyrskii deep mine, Russia (Figures 1 and 2). Zonal
disintegration poses a great danger during the excavation of deep tunnels [7]. It is a character of the deep rock mass and has recently been a subject of focus. Many specialists have implemented several types of methods to explain the phenomenon. Sellers and Klerck [8] believed that discontinuous surfaces could be one of the origins of zonal disintegration. Malan and Spottiswoode [9] analyzed the relationship between shock bumps and zonal disintegration. Zhou et al. [10] studied the dynamic excavation of deep tunnels and determined the forming time of fractured zones. Gu et al. [11] carried out compressive tests on cylindrical specimens and believed axial stress to be the important cause of zonal disintegration. Some specialists have used nonequilibrium thermodynamics (Metlov L S, 2002), Hamilton time domain variation (LI SC, 2009), or a non-Euclidean model (GUZEV M A, 2001) to study the forming mechanism.
2
The Scientific World Journal 𝜌( Ω· m
10 0 80 0 60 0 40 0 20 0 0 0
2. Materials and Methods
900
1
5
3
5
1200
𝜌 (Ω·m)
PEZN122 PEZN223
# 2 3#
4#
h (m) 5 3 1
1#
0 200 400 600 800
7
1
3
600
1
3
5 7 h (m)
9
300
𝜌 (Ω·m)
0 30 0 15
)
45
0
) Ω·m 𝜌(
0
PEZN224 PEZN221
Figure 1: Zonal disintegration of Taimyrskii mine [4].
N
40 20
2.1. Introduction to the Model Test and Numerical Simulation. The two methods of physical and numerical simulation have their own merit and are complementary to each other. During the numerical simulation, virtual model is established in computer to research the actual engineering. It can be used to simulate the anisotropism, anisotropic and discontinuity characters of the medium, and the complex boundary condition in engineering. Compared to numerical simulation, model test can reflect failure process in rock mass visually and truly. The analogical model test is an effective reduced-scale method for researching special engineering based on similarity theory. The model is constructed in a manner similar to that in which engineering prototypes and thus the deformation laws can be monitored using precision devices. The data from the model can be converted to that of an engineering prototype using similarity theory to reveal the stress distribution. Hence, realworld problems can be solved using this methodology. The different advantages can be used and complementary to each other and thus the mechanism of the zonal disintegration can be revealed. 2.2. Model Test on Zonal Disintegration. The similarity theory required that the following similarity criteria must be satisfied in mechanical modeling (Fumagalli [15, 16]): 𝐶𝜎 = 𝐶𝛾 𝐶𝐿 ,
Figure 2: Relationship between the damage of support and the distribution of zonal disintegration [4].
𝐶𝛿 = 𝐶𝜀 𝐶𝐿 , 𝐶𝜎 = 𝐶𝜀 𝐶𝐸 ,
(1)
𝐶𝜀 = 𝐶𝑓 = 𝐶𝜙 = 𝐶𝜇 = 1, However, the mechanical behaviour of deep rock mass is nonlinear, very complex, and clearly different from the engineering response observed in shallow embedded tunnel engineering. The zonal structure of fracturing discovered does not fit within the framework of the conventional theoretical models. Indeed, there is no convincing explanation, but some arguments on the forming condition and failure mode have been made. For example, Oparin and Kurlenya [12] believe fractured zone are circular and concentric to the tunnel periphery and developed a formula for calculating the radii of fracture zones. On contrary, Borzykh [13] believes that fractured zones are a group of plastic slippage lines that form after deep rock mass yields plastically. In contrast to the results obtained in the above-mentioned studies, Tang and Zhang’s [14] numerical simulation results show that the fractured zones appear more like spirals. Therefore, it is necessary to research the forming process and fracture pattern and thereby determine the underlying mechanism. Simulation on the excavation process of tunnel in laboratory, using numerical method, is the most effective method. Based on this, the forming process and phenomena of zonal disintegration can be simulated and the mechanism can be revealed accordingly.
where 𝐶𝐿 , 𝐶𝜎 , 𝐶𝜀 , 𝐶𝛾 , 𝐶𝐸 , 𝐶𝜇 , 𝐶𝜌 , 𝐶𝑓 , 𝐶𝑐 , and 𝐶𝜑 represent the similarity constants for geometry, stress, strain, volume weight, deformation modulus, Poisson’s ratio, density, friction coefficient, and cohesion and friction angle, respectively. To carry out the analogical model test, a suitable analogical material is required. The material determines whether the model test can reflect the mechanical response of the engineering prototype. Barite powder, iron powder, and quartz sand are used for the aggregate, and an alcohol solution of rosin is used as the mucilage glue. Through hundreds of sets of proportioning tests, a material referred to as Barites-IronSand cementation analogical material was developed. The proportion of the aggregates and the concentration of the alcohol solution of rosin determine the mechanical behaviour of the material. The proportion of the aggregates and the concentration of the alcohol solution of rosin decide the mechanical behaviour of the material. The similarity ratio of volume weight for the similar material is set to 1 : 1, while the similarity ratio of geostress is 1 : 25. Thus, the mechanical parameters of the mudstone and the corresponding similar material are shown in Table 1. Construction of the model: the model size is limited by the reasonable size of the steel frame and 30 were taken as the optimal similarity coefficient 𝐶𝐿 in the current study.
The Scientific World Journal
3
0. 1 5
2
2
3
4
4
5
6
7
8
9
8 6 5 3
10 9 8 7 6 5 4 3 2 10.5 h (m)
h (m)
1 5 0.
Fracture rock
7
h (m)
Intact rock
9
10
Fracture rock
10 9 8 7 6 5 4 3 2 1 0.5
45∘ 5000 mm
0.51 2 3 4 5 6 7 8 9 10 h (m)
Intact zone Fractured zone
Figure 3: Zonal disintegration of Dingji coal mine of Huainan mine area.
Table 1: Physical-mechanical parameters of the prototype rock and model material. Material type Prototype Volume Weight/KN⋅m−3 2.62 Edef/Mpa 5250 Cohesion/Mpa 10 Friction angle 43 Uniaxial compressive strength/Mpa 29.08 Tensile strength/Mpa 10.22 Passion ratio 0.272
Similar material 2.62 175 0.4 43 1 0.674 0.272
The simulation range of the prototype was 11 m. According to the geometry similarity scale, the dimension of the model was determined to be 0.45 m and the model tunnel was 160 mm. The model was delaminated. The similar material was placed inside the steel drum in layers and then tamped to create the model. The steel drum was 400 mm high, and the wall was 10 mm thick and had an inner diameter of 450 mm. Each layer was 100 mm high, and thus four layers were needed to complete the model. Then, a circular cavern cylinder was preset inside the steel drum to make the tunnel. The cavern diameter was 160 mm and the cavern axis coincided with the central axes of the steel drum. Then the similar material was put inside the steel drum and compressed uniformly to build the model. The measuring components were set up when the material reached the desired height, to monitor the displacement in the surrounding rocks (Figure 3). The model was split and the cracks distribution inside the model was then determined (Figure 4). Circular cracks propagated and coalesced to form a circular pattern. Three to four fracture zones were observed to surround the cavern, in accordance with the fractured shape obtained from the field
Figure 4: Loading procedure of model test.
monitoring. Thus, the occurrence of zonal disintegration was confirmed in the model. As Figure 5 shows, there appears about 3-4 fractured around the tunnel, which is the zonal disintegration. And during the zonal disintegration, it is found that the deformation of the rock mass shows a waved shape. The details could be found in the paper [17, 18]. 2.3. Introduction to XFEM (Extended Finite Element Method). In recent 5 years, many specialists have tried all kinds of numerical methods to simulate the zonal disintegration. Tang and Zhang [14] adapted RFPA and Qian et al. [19] and Wang [20] adapted strain soften model. However, the results are spiral, which is not consistent with the definition of zonal disintegration. So it is needed to find a new numerical method to research it. The elements are needed to subdivide continuously during the crack propagation when using FEM. So it is not
Measuring points
Tunnel 6 5 4 3 2 1 Measuring points
1 2 3 4 5 6
Measuring points
2 1.5 1 0.5 0
Displacement (mm) 2 1.5 1 0.5 0 6 5 4 3 2 1
1 2 3 4 5 6 Measuring points
2 1.5 1 0.5 0
Displacement (mm)
The Scientific World Journal
Displacement (mm)
4
2 1.5 1 0.5 0 Displacement (mm)
(a) Cracks distribution
(b) Displacement distribution surrounding the cavern
Figure 5: Results of the model test.
Applying load
Stiffness matrix by penalty function
Computing stress field
Stress factor of fracture
Extension angle of fracture
Figure 7: The numerical model of crack extension trajectory under uniaxial compression. Decide if the fracture extends
No
Yes The fracture extends a step
Figure 6: Numerical procedure of the extended finite element method (XFEM).
proper to the fracture problem in rock engineering. In 1999, Belytschko and Black [21] developed the XFEM. It is a newly developed means to simulate the fracture propagation for the rock materials. It can be considered as the extension of
regular FEM. It has both the advantages of regular FEM and the special character of itself: (1) the fracture is not need to be considered. The element can be subdivided in any part, including the fracture; (2) the grid is not need to subdivide again, which improves the efficiency of computing highly; (3) the shape function can be changed according to the problem in study, which is very flexible. Many specialists have paid attention to XFEM and applied it to the rock fracture analysis. XFEM is adapted in this study
The Scientific World Journal
(
S, mises : 75%) +4.659e + 03 +4.490e + 03 +4.322e + 03 +4.153e + 03 +3.984e + 03 +3.815e + 03 +3.646e + 03 +3.477e + 03 +3.308e + 03 +3.140e + 03 +2.971e + 03 +2.802e + 03 +2.633e + 03 +2.464e + 03 +2.295e + 03 +2.126e + 03 +1.958e + 03 +1.789e + 03 +1.620e + 03 +1.451e + 03 +1.282e + 03 +1.113e + 03 +9.445e + 02 +7.757e + 02 +6.068e + 02
5
(
S, mises : 75%) +8.857e + 03 +8.545e + 03 +8.233e + 03 +7.921e + 03 +7.609e + 03 +7.297e + 03 +6.984e + 03 +6.672e + 03 +6.360e + 03 +6.048e + 03 +5.736e + 03 +5.424e + 03 +5.112e + 03 +4.799e + 03 +4.487e + 03 +4.175e + 03 +3.863e + 03 +3.551e + 03 +3.239e + 03 +2.927e + 03 +2.615e + 03 +2.302e + 03 +1.990e + 03 +1.678e + 03 +1.366e + 03
(a) S, mises ( : 75%) +5.351e + 04 +5.153e + 04 +4.955e + 04 +4.757e + 04 +4.559e + 04 +4.361e + 04 +4.163e + 04 +3.964e + 04 +3.766e + 04 +3.568e + 04 +3.370e + 04 +3.172e + 04 +2.974e + 04 +2.776e + 04 +2.578e + 04 +2.380e + 04 +2.182e + 04 +1.984e + 04 +1.785e + 04 +1.587e + 04 +1.389e + 04 +1.191e + 04 +9.931e + 03 +7.950e + 03 +5.970e + 03
S, mises : 75%) +2.303e + 04 +2.219e + 04 +2.135e + 04 +2.051e + 04 +1.968e + 04 +1.884e + 04 +1.800e + 04 +1.716e + 04 +1.632e + 04 +1.549e + 04 +1.465e + 04 +1.381e + 04 +1.297e + 04 +1.213e + 04 +1.130e + 04 +1.046e + 04 +9.620e + 03 +8.782e + 03 +7.944e + 03 +7.106e + 03 +6.268e + 03 +5.430e + 03 +4.592e + 03 +3.754e + 03 +2.916e + 03
(
04
(b) (
S, mises : 75%) +1.299e + 05 +1.251e + 05 +1.202e + 05 +1.154e + 05 +1.106e + 05 +1.058e + 05 +1.010e + 05 +9.614e + 04 +9.133e + 04 +8.651e + 04 +8.169e + 04 +7.687e + 04 +7.205e + 04 +6.723e + 04 +6.241e + 04 +5.760e + 04 +5.278e + 04 +4.796e + 04 +4.314e + 04 +3.832e + 04 +3.350e + 04 +2.868e + 04 +2.386e + 04 +1.905e + 04 +1.423e + 04
(c) (
(d)
S, mises : 75%) +7.611e + 09 +7.308e + 09 +7.005e + 09 +6.702e + 09 +6.399e + 09 +6.096e + 09 +5.793e + 09 +5.490e + 09 +5.187e + 09 +4.884e + 09 +4.581e + 09 +4.278e + 09 +3.975e + 09 +3.672e + 09 +3.369e + 09 +3.066e + 09 +2.763e + 09 +2.460e + 09 +2.157e + 09 +1.854e + 09 +1.551e + 09 +1.248e + 09 +9.447e + 08 +6.417e + 08 +3.387e + 08
(e)
(f)
Figure 8: Crack propagated of rock sample under compression simulated by XFEM.
with the criteria proposed to simulate the zonal disintegration. The procedure can be shown (Figure 6).
The angle of crack initiates is
2.4. Simulations on the Crack Development. The simulations on the crack development are carried to examine the efficiency of XFEM, which is shown in Figures 7 and 8.
𝜃0 = arccos
2.5. Simulations on Zonal Disintegration. The numerical tests are carried out on the rock samples containing hole with XFEM, and zonal disintegration is reproduced (Figure 10). The maximum circumference tension stress criterion was taken as the initiate criterion. The critical loading depends on the equation 𝐾𝑒 = cos
𝜃 𝜃 3 [𝐾𝐼 cos2 ( ) − 𝐾𝐼𝐼 sin 𝜃] ≥ 𝐾𝐼𝐶. 2 2 2
2 2 + √𝐾𝐼2 + 8𝐾𝐼2 𝐾𝐼𝐼 3𝐾𝐼𝐼 2 𝐾𝐼2 + 9𝐾𝐼𝐼
(3)
where 𝐾𝐼𝐼 = 0. Substitute it to (2) and (3), then: 𝐾𝐼 = 𝐾𝐼𝐶.
(4)
The stress intensity factor of circular fracture is [17] 𝐾1 (𝑟) =
2𝑃0 √𝜋 √𝑟 sin 𝜃 (3 − cos 𝜃)
(2)
That is, the fracture begins to develop when the equivalent stress intensity factor 𝐾𝑒 is equal to or larger than the fracture toughness of the surrounding rocks 𝐾𝐼𝐶.
,
× {(1 +
𝑟𝑎2 2𝑟𝑎 (3 cos 𝜃 − 1) ) 𝑟2 3𝜋2
− (1 −
𝑟𝑎2 𝜃 ) 𝑟 sin 𝜃 cos } . 2 𝑟 2
(5)
6
The Scientific World Journal
Stress inside the rockmass got from the plastic model
No Decide if the element stress reaches the value to the circular fracture extension by equation (5)
Accomplishment
Yes Circular fracture departures with the rock mass, that is, the excavation, the computing restarts
Yes
Fracture developes the step, direction decided by equation (6)
Circular fracture forms
No Element stress computing after crack extension
Figure 9: Block diagram of the simulating procedure of zonal disintegration.
(
Z
R T
S, mises : 75%) +6.675e + 07 +6.265e + 07 +5.854e + 07 +5.444e + 07 +5.033e + 07 +4.623e + 07 +4.212e + 07 +3.801e + 07 +3.391e + 07 +2.980e + 07 +2.570e + 07 Y +2.159e + 07 X +1.748e + 07
Y X
(a) Numerical model
(b) Forming process of zonal disintegration
Figure 10: Zonal disintegration of circular model obtained with XFEM.
Substitute (5) into (4). The criteria of zonal disintegration can be got by
the equivalent stress intensity factor of fracture; 𝐾𝐼𝐶 is the fracture toughness; 𝜃0 is the angle of crack initiates.
2𝑟𝑎 (3 cos 𝜃 − 1) 𝜃 𝜎1 − 𝑟𝜎2 sin 𝜃 cos 2 3𝜋 2
2.6. Simulation Procedure. The numerical process of zonal disintegration is shown in Figure 9.
𝐾 √𝑟 sin 𝜃 (3 − cos 𝜃) ≥ 1𝐶 . 2√𝜋
(6)
The fracture angle is 𝜃0 = 0, which means that the crack will initiate along the plane where it exists, where 𝐾𝑒 is
2.7. Numerical Simulation of Zonal Disintegration of Deep Cavern. In order to research the forming mechanism deeply, the research on the damage procedure of zonal disintegration is carried on based on the geomechanical model test.
The Scientific World Journal
7
(
T ZR
S, mises : 75%) +9.305e + 07 +8.703e + 07 +8.102e + 07 +7.500e + 07 +6.899e + 07 +6.298e + 07 +5.696e + 07 +5.095e + 07 +4.494e + 07 +3.892e + 07 +3.291e + 07 +2.689e + 07 +2.088e + 07
(a) Model
(
(b) The 21 steps
S, mises : 75%) +9.305e + 07 +8.703e + 07 +8.102e + 07 +7.500e + 07 +6.899e + 07 +6.298e + 07 +5.696e + 07 +5.095e + 07 +4.494e + 07 +3.892e + 07 +3.291e + 07 +2.689e + 07 +2.088e + 07
(c) Computing done
Displacement (cm)
Figure 11: Zonal disintegration phenomenon in surrounding rock mass simulated with XFEM.
Displacement surrounding cavern
6.0 5.5 5.0 4.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0
1 A B C D E
2
3 Spots F G H I
4
5 J K L M
Figure 12: Displacement distribution surrounding the cavern.
The diameter of model is 0.45 m, with the circle cavern of 160 mm in diameter, which is the same as the model test. The model is divided into 30000 elements (four-node). In order to simplify the issue, the hydraulic pressure is applied on the model, which is 0.5 MPa. The numerical model are shown in Figure 8. The result and displacement surrounding cavern are in Figure 9. As Figures 10 and 11 show, there exist 4 circular fractures around the cavern. The result is similar to the model test. So the zonal disintegration is proved by using numerical simulation. And as Figure 12 shows, the deformation law is nonlinear like the result of model test. 2.8. Comparison between Model Test and Numerical Simulation. As shown in the photos of the fracture shapes of model test and numerical simulation, it can be found as follows.
Circular cracks propagated and coalesced to form a circular pattern. Four fracture zones were observed to surround the cavern, in accordance with the fractured shape obtained from the model test. Thus, the occurrence of zonal disintegration was confirmed in the model. Comparing the fracture shape in numerical model (Figure 11) with that in the test (Figure 5), it is found that a system of cracks geometrically similar to the periphery of the tunnel characterised the zonal structure observed both in numerical simulation and the model. Based on the fracture shape of the model, the fracture radius was related to the cavern radius. Suppose that the radii of the fractured zone are 𝑟𝑖 (𝑖 = 0, 1, 2, . . . , 𝑛), and the radius of the cavern is 𝑟0 . Thus, the radius of the fracture zones fulfills the following relation: 𝑟𝑖 = 𝑟0 ⋅ 𝛼𝑖 .
(7)
That is, 𝑟 𝑟1 𝑟2 = = ⋅ ⋅ ⋅ = 𝑛 = 𝛼, 𝑟0 𝑟1 𝑟𝑛−1
(𝑖 = 0, 1, 2, . . . , 𝑛) .
(8)
Some specialists believe that 𝛼 is a constant value; that is, 𝛼 = √2 Oparin and Kurlenya [12]. However, after measuring the disassembled model, it was observed that 𝛼 is not equal to √2; in actuality, 𝛼 ≈ 1.2 in this paper. 2.9. Displacement Law in Surrounding Rocks. Figure 11 shows the displacement surrounding the cavern. The displacement in each measurement line presents an oscillating law in which the peak and valley values alternate with distance. This law is completely different from the law observed for shallow embedded tunnels, in which the displacement decreases monotonically with increasing distance from the tunnel wall. By comparing with the fracture shape of the model, it was observed that the region of peak values showing greater displacement is the fracture zone, whereas the region of minimum values showing less displacement is the intact zone. The additional displacement in the peak value area caused by the circular fracture increases the total displacement. Comparing Figure 11 with Figure 5 reveals that the displacement law obtained from the numerical simulation is analogous to that obtained from model test. They both show the nonmonotonic law.
8
The Scientific World Journal
3. Conclusions (1) Zonal disintegration phenomenon is simulated by XFEM. (2) The criteria of crack propagation are got. And it is appropriate for the simulation of zonal disintegration. (3) The results between the geomechanics model test and numerical simulation are consistent. Both the two methods show that the fracture line is the concentric circles of the cave. There are 4 fractured zones that appear surrounding the cavern. And the radiuses of the fractured zones fulfill the geometric progression.
Conflict of Interests The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments This work was supported by the National Natural Science Foundations of China (Grant no. 51209074); the China Postdoctoral Science Foundation (Grant no. 2012M511189 and 2013T60494); the Fundamental Research Funds for the Central Universities (Grant no. 2012B02714), supported by the State Key Laboratory For Geomechanics and Deep Underground Engineering, China University of Mining and Technology, under Grant SKLGDUEK1206; the Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering, the Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, under Grant no. Z012008, funded by CRSRI Open Research Program CKWV2012306/KY; and the Key Laboratory of Coal-based CO2 Capture and Geological Storage Open Research Program 2012KF08. The authors are deeply grateful for this support.
References [1] E. I. Shemyakin, G. L. Fisenko, M. V. Kurlenya et al., “Zonal disintegration of rocks around underground workings—part 1: data of in situ observations,” Soviet Mining Science, vol. 22, no. 3, pp. 157–168, 1986. [2] Q. H. Qian, “The characteristic scientific phenomena of engineering response to deep rock mass and the implication of deepness,” Journal of East China Institute of Technology, vol. 27, no. 1, pp. 1–5, 2004. [3] G. R. Adams and A. J. Jager, “Etroscopic observations of rock fracturing ahead of the stope faces in deep-level gold mines,” Journal of the South African Institute of Mining and Metallurgy, vol. 21, no. 2, pp. 115–127, 1980. [4] E. I. Shemyakin, G. L. Fisenko, M. V. Kurlenya et al., “Zonal disintegration of rocks around underground workings—part II: rock fracture simulated in equivalent materials,” Soviet Mining Science, vol. 22, no. 4, pp. 223–232, 1986. [5] E. I. Shemyakin, G. L. Fisenko, M. V. Kurlenya et al., “Zonal disintegration of rocks around underground mines—part III: theoretical concepts,” Soviet Mining Science, vol. 23, no. 1, pp. 1–6, 1987.
[6] E. I. Shemyakin, M. V. Kurlenya, V. N. Oparin et al., “Zonal disintegration of rocks around underground workings. IV: practical applications,” Soviet Mining Science, vol. 25, no. 4, pp. 297–302, 1989. [7] B. V. Laptev and R. P. Potekhin, “Burst triggering by zonal disintegration of evaporites,” Soviet Mining Science, vol. 24, no. 3, pp. 238–241, 1988. [8] E. J. Sellers and P. Klerck, “Modelling of the effect of discontinuities on the extent of the fracture zone surrounding deep tunnels,” Tunnelling and Underground Space Technology, vol. 15, no. 4, pp. 463–469, 2000. [9] D. F. Malan and S. M. Spottiswoode, “Time-dependent fracture zone behavior and seismicity surrounding deep level stopping operations,” in Rockbursts and Seismicity in Mines, S. J. Gibowicz and S. Lasocki, Eds., pp. 173–177, A. A. Balkema, Rotterdam, The Netherlands, 1997. [10] X. Zhou, Q. Qian, and B. Zhang, “Zonal disintegration mechanism of deep crack-weakened rock masses under dynamic unloading,” Acta Mechanica Solida Sinica, vol. 22, no. 3, pp. 240– 250, 2009. [11] J. Gu, L. Gu, A. Chen, J. Xu, and W. Chen, “Model test study on mechanism of layered fracture within surrounding rock of tunnels in deep stratum,” Chinese Journal of Rock Mechanics and Engineering, vol. 27, no. 3, pp. 433–438, 2008. [12] V. N. Oparin and M. V. Kurlenya, “Gutenberg velocity section of the earth and its possible geomechanical explanation. I. Zonal disintegration and the hierarchical series of geoblocks,” Journal of Mining Science, vol. 30, no. 2, pp. 97–108, 1994. [13] A. F. Borzykh, “Features of the zonal disintegration of roof rocks and a coal seam around mine workings,” Soviet Mining Science, vol. 26, no. 5, pp. 418–427, 1990. [14] C. Tang and Y. Zhang, “Discussion on mechanism and evolution laws of fracture spacing in rock mass,” Chinese Journal of Rock Mechanics and Engineering, vol. 27, no. 7, pp. 1362–1369, 2008. [15] E. Fumagalli, Statical and Geomechanical Models, Springer, New York, NY, USA, 1973. [16] E. Fumagalli, “Geomechanical models of dam foundation,” in Proceedings of the International Colloquium on Physical and Geomechanical Models, Bergamo, Italy, 1979. [17] C. Xu-Guang, Z. Qiang-Yong, and W. Yuan, “Model test of anchoring effect on zonal disintegration in deep surrounding rock masses,” The Scientific World Journal, vol. 2013, Article ID 935148, 16 pages, 2013. [18] C. Xuguang and W. Yuan, “Energy dissipation and release during rock burst,” Disaster Advances, vol. 9, no. 6, pp. 48–54, 2013. [19] Q. H. Qian, X. P. Zhou, H. Q. Yang, Y. X. Zhang, and X. H. Li, “Zonal disintegration of surrounding rock mass around the diversion tunnels in Jinping II Hydropower Station, Southwestern China,” Theoretical and Applied Fracture Mechanics, vol. 51, no. 2, pp. 129–138, 2009. [20] H.-Y. Wang, “Numerical simulation on disintegration of surrounding rock mass in deep mine roadways,” Journal of the China Coal Society, vol. 35, no. 4, pp. 535–540, 2010. [21] T. Belytschko and T. Black, “Elastic crack growth in finite elements with minimal remeshing,” International Journal for Numerical Methods in Engineering, vol. 45, no. 5, pp. 601–620, 1999.
