Research Article  Open Access
Xiangfeng Lv, Hongyuan Zhou, Aiwen Wang, Chun Feng, Xiaochun Xiao, "Characteristics of Stress Transfer and Progressive Fracture in Overlying Strata due to MiningInduced Disturbances", Advances in Civil Engineering, vol. 2018, Article ID 8967010, 13 pages, 2018. https://doi.org/10.1155/2018/8967010
Characteristics of Stress Transfer and Progressive Fracture in Overlying Strata due to MiningInduced Disturbances
Abstract
In this study, based on the mining of the 13210 working face in the Yima coal mine of the Gengcun village, China, a simplified mechanical model for the analysis of dynamic destabilization of the overlying strata during underground mining was constructed. The numerical simulation was used to analyze the stress patterns in the advanced abutments of the tunnel face and the characteristics of dynamic failures in the overlying strata. Furthermore, similitude experiments were conducted to study the process of stress release and deformation in the overlying strata, and to analyze the effects of overburden destabilization on the ground surface settlement. The theoretical analysis indicated that if the geometric parameters of a working face are fully determined, a stiffness ratio no greater than 1 is required for dynamic destabilization to occur. The numerical simulation results show that the stress in the overlying strata decreases with a decrease in distance from the tunnel face. The stresses in the advanced abutments initially increase with an increase in distance from the tunnel face, followed by a decrease in stress, and an eventual stabilization of the stress levels; this corresponds to the existence of a “stress buildup zone,” “stress reduction zone,” and “native rock stress zone.” In similitude experiments, it was observed that a “pseudoplastic beam” state arises after the local stresses of the overlying strata have been completely released, and the “trapezoidal” fractures begin to form at stress concentrations. If the excavation of the working face continues to progress, the area of collapse expands upward, thereby increasing the areas of the fracture and densification zones. Owing to the nonuniform settlement of the overlying strata and the continuous development of bedseparating cracks, secondary fractures will be generated on both sides of the working face, which increase the severity of the ground surface settlement.
1. Introduction
In recent years, problems caused by underground mining, such as surface subsidence, rock burst failures in the working face, and overlying strata fractures, have become increasingly prominent [1]. Since the inception of mining subsidence studies in the beginning of the twentieth century, a large body of research works has accumulated on mininginduced ground settlement and overlying strata failures. These research works generally include theoretical analyses, numerical calculations, and experimental studies. The theoretical studies on mining subsidence helped develop the masonry beam theory [2], key strata theory [3, 4], and supporting plate theory [5]. The numerical approaches may be generally divided into two categories. The first category pertains to methods based on continuum mechanics, for example, the finite element method (FEM), finite difference method (FDM), boundary element method (BEM), and elementfree method (EFM). The second category of methods is based on discontinuum mechanics, which includes the discrete element method (DEM) and discontinuous deformation analysis (DDA) [6–8]. In experimental studies, the similitude experiment is the most broadly used experimental tool for studying ground subsidence. These experiments are straightforward to perform and provide an accurate reflection upon mininginduced ground deformations and change patterns in overburden stresses.
In this study, we briefly discuss the most important gaps in the current literature regarding stress transfers and failures in the overlying strata during underground mining processes. In the theoretical aspect, the mechanisms of rock burst failures in the overburden during underground mining remain unclear, and a consistent understanding is yet to be achieved on the overburden’s deformation and failure processes, or its stress distribution patterns. In numerical simulations, the FEM and FDM are highly effective methods for simulating the properties of a material in continuum states. However, these methods are incapable of simulating the process in which a material changes from a continuous to discontinuous state, and they cannot reveal the dynamics of a material in a discontinuous state. The blockdiscrete element method (BDEM) and granular DEM have certain advantages in the simulation of noncontinuum dynamics; however, it is still difficult to simulate the continuous deformation processes of a material using these methods [9]. Although similitude experiments can produce patterns that are consistent with actual engineering works, these experiments lack in precision. Consequently, the similitude experiment results often have to be verified against theoretical analyses or numerical calculations.
In this study, we constructed a mechanical model for the analysis of dynamic failures in a working face, in the presence of disturbances caused by underground mining, based on the effects of working face excavation on the stress states of the overlying strata. Further, this model was used to provide the criterion for the dynamic failures in the overburden of the working face. Based on the numerical methods of the continuumdiscontinuum element method (CDEM), FEM and DEM were coupled to each other to simulate the deformations and dynamics of the overburden in continuous and discontinuous states. This allows the progressive failure of the overlying strata to be analyzed as it transitions from a continuous mass to a discontinuous mass. The similitude experiments were conducted to investigate the stress relief and deformation in the overburden and to analyze the effects of overburden failure on the ground settlement. Finally, the dynamics of stress transfers and failures in the overburden during mining process were obtained via a comprehensive analysis of the results obtained from our theoretical analyses, numerical calculations, and similitude experiments.
2. A Dynamic Model for Underground Resource Mining
During the underground mining processes, the stress state of the overburden changes from the native state to a state with rising axial stresses and decreasing confining stresses, until failure finally occurs. The overburden failure and destabilization is triggered by the superposition of dynamic and static loading; this is already broadly recognized and accepted as a dynamic problem [10–12]. In [13–15], it was observed that the excavationinduced unloading transforms the triaxial stress state of the overburden around the working face into a biaxial stress state, causing a rapid release in the stresses that were concentrated within the overburden, thereby triggering dynamic failures. Based on this analysis, we devised a simplified mechanical model for the analysis of dynamic overburden failures during the excavation of a working face using the 13210 working face in the Yima coal mine of the Gengcun village as an example, as shown in Figure 1.
Supposing that is the pressure acting on the roof, is the stiffness of the roof, is the deformation of the roof, is the load acting on the elastoplastic zone, is the stiffness of the elastoplastic zone, is the load acting on the elastic zone, is the stiffness of the elastic zone, is the total sinkage of the roof, is the deformation of the excavation layer, and the roof may be treated as an elastic body. The conditions for equilibrium are as follows:
In [16], it was observed that strain softening occurs in the elastoplastic zone of rock masses. Therefore, a nonlinear constitutive equation is used here:where is the initial elastic modulus and is the peak strain.
The elastoplastic zone of the excavation strata has a width of and a thickness of . The relationship between load and deformation under unit length conditions is as follows:where is the initial stiffness and is the settlement of the roof when the roof’s strain reached the peak value. The second derivative of (3) was taken to determine the inflection point of the loadversusdeformation curve for the elastoplastic zone; an inflection point occurs at , and the absolute gradient of the curve is . Refer to [13] to determine the energy composition of the overlying strata. The system strain energy consists of the top plate and excavation layer strain energies. The total potential energy is as follows:
The first derivative of (4) was set to 0 to obtain the equilibrium surface equation:
The deformation at the cusp of the equilibrium surface, , is exactly the same as the inflection point of the elastoplastic zone’s constitutive curve. The Taylor expansion was performed on (5) with up to the third term, which gives
The point mutation theory discussed in [17] was used to obtain the standard form of the surface equation:where is the state variable of the system, while and are the control variables:where is a dimensionless parameter that is a function of . The bifurcation set is as follows:
The substitution of (8), (9), (10), and (11) into (12) yields the following:
Equation (13) is the mechanical criterion for the occurrence of rock burst failures in this system. Based on (13), it can be inferred that a sudden change occurs in the system when ; which implies that is a prerequisite for the occurrence of sudden destabilizations in the system. Therefore, one of the conditions for the occurrence of dynamic failures in the overlying strata of this study is a stiffness ratio no greater than 1. Furthermore, a dynamic failure can occur only if soft softening has occurred in the overburden after the geometric parameters were fully determined.
3. ContinuousDiscontinuous Deformation Theory of Overlying Strata during Underground Mining
3.1. Basic Principle
According to the dynamic model for underground resource mining, the continuumbased discrete element theory is adopted is this study [18]. It contains two kinds of elements: blocks and contacts [19] (Figure 2). A discrete block consists of one or more FEM elements that share the same nodes and faces. A contact contains several normal and tangent springs, and each spring owns two nodes that belong to two different blocks. The FEM is used inside a block; while the DEM is adopted for a contact.
3.2. Governing Equation
In the discontinuous deformation theory, the governing equation is the motion equation. The block is subjected to internal and external forces. The internal force includes the damping and deformation forces. The external force includes the external boundary force and the force between springs. As the block body is regarded as a continuous and isotropic linear elastic body, its mechanical properties can be described by the basic differential equations of threedimensional elastodynamics [20, 21], which are as follows:
equilibrium equation:
geometric equation:
physical equation:
Boundary condition: (on the displacement boundary of ) and (on the force boundary of ).
Here, , , , and are the stress, displacement, volume force, and area force, respectively; and denote the rock block region and its boundary, respectively; , and are the Lame constants; and denote the mass density and damping coefficient, respectively; and is the Kronecker delta function.
Based on the elastic variational principle [20], the governing equation is as follows:where , , and are the acceleration, speed, and displacement arrays, respectively. , , , and are the mass matrix, damping matrix, stiffness matrix, and nodal load array, respectively.
When we solve the governing equations, the first step is to loop each deformable block and calculate the corresponding continuous deformation. The second step is to calculate the surface force. When we calculate the surface force, the stiffness matrix and nodal displacement are used to calculate the elastic force; the damping matrix and nodal velocity are used to calculate the damping force. Finally, the equation of motion is solved using the direct integral method. The specific equations are as follows:
elastic force:
damping force:
Solution of the equation of motion using the direct integral method:
The acceleration, velocity, and displacement of the block body nodes can be calculated by (20), where includes the interface and contact forces. It is observed that the interface force is determined by the boundary conditions.
3.3. Model Boundary
Based on the block contact theory proposed in [19], a mechanical model is proposed to describe the interaction between the two blocks. When the deformation of the material is continuous, the force between adjacent blocks is calculated based on the nodal displacement:where is the normal force of the spring, is the normal stiffness of the spring, is the normal displacement vector difference for the springconnected blocks, is the tangential force of the spring, is the tangential stiffness of the spring, and is the tangential displacement vector difference for the springconnected blocks.
When the deformation of the material is discontinuous, the fracture of the spring is used to simulate deformation. The spring has two types of damage: tensile and shear failures. When the normal force of the spring is greater than the tensile strength of the material, the material suffers from tensile failure. The normal and tangential forces between the adjacent blocks are zero, that is,
After the spring suffers from the shear failure, the blocks will slide relative to each other; thus, there is pressure and friction. Among them, the pressure is equal to the product of the normal stiffness of the spring and the difference of the displacement vector, and the friction is equal to the product of the normal force and the friction coefficient. Under these conditions, the direction of friction is opposite to the motion direction of the block. Thus,where is the friction coefficient and is the internal friction angle of the material.
3.4. Calculation Steps
The dynamic relaxation method discussed in [22] is used to solve equations. The static problem is transformed into a dynamic problem by introducing damping. Figure 3 shows the main process of this method.
4. Numerical Simulation and Physical Model Test for Underground Mining
4.1. Calculation Method and Model Parameters
The CDEM was used to perform the numerical calculations in this study. The numerical model shown in Figure 4 is established according to the actual size of the 13210 working face at the YIMA Gengcun (China) coal mine. The model size is set as 300 × 150 × 83 m. The model is divided into four layers, from top to bottom, as conglomerate, roof, coal, and floor; and their thicknesses are 44 m, 15 m, 9 m, and 15 m, respectively. The horizontal and vertical stresses are 20 and 17 MPa, respectively, and the depth of each mining is 4 m. The rectangular roadway section has a length of 5 m and a width of 4 m. The bottom surface of the model is fully constrained, the horizontal level of the model is surrounded, and the roof is a free surface. The working face is shown in Figure 4, and the numerical calculation model is shown in Figure 5. Referring to the arrangement of the stress and displacement monitors (#1, #2, #3, and #4) shown in Figure 5, the monitor points for the numerical simulation were set. In Figure 5, the different colors represent the geologic bodies with different mechanical properties. Groups 1–4 represent the conglomerate, roof, coal, and floor, respectively. Groups 5–17 represent the raw coal with different mechanical properties in front of the working face. The basic mechanical parameters of groups 1–17 are shown in Table 1.
(a)
(b)

4.2. Stress Transfer Law and Distribution Characteristics of Overlying Strata during the Working Face Mining
The stresses of the advanced abutments and the deformational patterns of the overlying strata during mining conditions are shown in Figure 6. Figure 6(a) shows that the deformation of the overburden is relatively small during the initial stages of working face mining (advancement up to 4 m), and the deformation of the overburden increases with a decrease in distance from the tunnel face. The deformation of the overburden decreases and eventually stabilizes with an increase in distance from the tunnel face. In addition, it is observed that the stress values decrease with a decrease in distance from the tunnel face. As this distance increases, the advanced abutments of the tunnel face exhibit an increase in stress, followed by a decrease, and finally the stabilization of stress levels; this is indicative of the existence of a “stress buildup zone,” “stress reduction zone,” and “native rock stress zone.” The stress monitoring data also indicate that the area of effect of the working face’s excavation is approximately 30 m, and the stress peak is located approximately 16 m in front of the tunnel face.
(a)
(b)
The continued excavation of the working face expands the area of the goaf and increases the deformation of the overlying strata. The stresses of the abutments gradually propagate toward the front of the working face, whereas a gradual decrease in stress occurs in the strata behind the excavation surface, accompanied by an increase in vertical displacement. This is because the strata overlying the goaf are beginning to sink under their own weight, which is because the yield strength of the strata is gradually being exceeded [23–25]. The overburden will ultimately collapse when its deformation reaches a critical value. As shown in Figure 6(b), the displacement of the overburden in front of the working face is relatively small when the excavation of the working face has advanced up to 48 m, whereas the displacement of the strata overlying the goaf in the back has increased. At the 40 m point, the displacement of the overburden reaches its maximum at 6 m and the corresponding stress is almost 31 MPa. This indicates that the overburden will completely collapse at this point of the roof; that is, the first weighting of the main roof occurs with a weighting interval of 40 m. The stress data also show that large stresses occur at the 50 m and 60 m points. Therefore, periodic collapses will occur in the overlying strata of the roof as the working face continues to be excavated, with a periodic weighting interval of 10 m.
4.3. Fracture Characteristics and Deformation Law of Overlying Strata during the Working Face Mining
The processes of failure and deformation in the overlying strata during working face mining are shown in Figures 7 and 8. The calculations indicate that local failures occur in the overburden of the coal seam when the working face excavation has progressed up to 4 m; when the working face reaches 8 m, the soft and weak strata of the false roof will collapse as excavation occurs. The immediate roof collapses when the excavation of the working face has advanced up to 20 m. When the excavation of the working face advances up to 48 m, the roof will fracture at the 40 m point; thus, forming the first weighting event. Therefore, the first weighting interval is approximately 40 m. The rock mass near the working face will collapse after the occurrence of bending deformations, and the cavedin rock masses will fill the goaf. The overlying strata that are far from the tunnel face will collapse on a large scale after losing their bearing capacity [26, 27]. When the excavation of the working face has reached 100 m, the rock masses that are far from the working face will collapse continuously. The area of collapse will expand upward, leading to an appearance of fractures and bed separations in the overburden. Further, the rock masses of the separated strata sink under their own weight. This is followed by the appearance of three distinct zones in the overburden of the goaf: the collapse zone, fracture zone, and gradual settlement zone. If the excavation of the working face progresses up to 180 m, the bedseparating cracks will develop upward and the nonuniform settlement of the overburden leads to the generation of secondary fractures. The locally developed secondary fractures will then propagate to the surface, and the area of collapse in the overburden will extend in both directions of the working face. A severe degree of surface settlement ultimately results from these processes.
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
4.4. The Progressive Fracture Characteristics and Deformation Law of the Overlying Strata under the Working Face Mining
A similitude model was constructed according to the methods discussed in [28–30] for studying the progressive failure and deformation of the overlying strata during the mining of the working face. Based on the similarity principles, the geometric, time, and stress similarity ratios were defined as 1/250, 1/2, and 450/1, respectively. The mechanical parameters of the prototype are listed in Table 1. The similarity materials and their mixing ratios were determined according to [31, 32]. The progressive fracture and deformation of the overlying strata in a variety of mining conditions are shown in Figure 9. Referring to the numerical simulation results, the model test results are analyzed, and the similarities are given. In Figure 9(a), it is shown that the overburden remained largely intact when the excavation of the working face advanced up to 4 m, without showing any significant signs of damage. Figure 9(b) shows that the softer and weaker sections of the roof began to cavein when the excavation of the working face reached 8 m; collapses also occurred in local sections of the overburden. In Figure 9(c), it is shown that the immediate roof collapsed when the excavation of the working face reached 20 m. The area of collapse in the overlying strata was expanded and the localized fractures that were formed previously began to propagate to the rest of the overburden. The bedseparating cracks have also formed in the overburden, which indicates that the stability of the overburden has begun to deteriorate. The settlement of the overburden also becomes quite apparent in this stage. Figure 9(d) shows that the entirety of the overburden begins to collapse when the excavation of the working face reaches 40 m, which indicates that the overlying strata exhibits “pseudoplastic beam”like mechanical characteristics upon the relief of local stresses in the overlying strata. This manifests as a downward sag at the center of the overburden and large displacements at the center. Fractures also occurred at stress concentrations on both ends of the working face, which led to the formation of a “trapezoidal” fracture. In addition, the bedseparating cracks continued to develop until they reached the surface of the overburden; thus, these cracks display a tendency to propagate upward along the overburden. Figure 9(e) indicates that in the absence of abutments, the entirety of the overburden will collapse when the excavation of the working face reaches 100 m. In this figure, it is shown that the area of collapse has expanded upward and formed a distinct “threezone” distribution. During the formation of this “threezone” distribution, the constraints imposed by the compression of the strata and the “trapezoidal” fracture formed around the fractures cause the boundaries of the “threezone” distribution to take a trapezoidal form. Thus, a trapezoidal compacted zone was formed in the deformed strata. In addition, the continued development of bedseparating cracks leads to the formation of 35–40 m wide fractured zones on both sides of the compacted zone. The continued excavation of the working face up to 180 m, under these conditions, led to the propagation of the collapsed zone to the ground surface, as shown in Figure 9(f). As the excavation continued, stress was transferred toward both sides of the overburden, which accelerated the development of fractures in the fractured zones. Some of the fractures in the fractured zone developed up to the ground surface. Furthermore, the nonuniform settlement of the overburden induced the generation of secondary fractures in the adjacent regions. Due to the disturbances caused by the secondary fractures and the subsidence of the overlying strata, one may surmise that severe ground surface settlement will occur after rock burst failures in the overburden [33, 34]. Based on [35–37], further excavation of the working face is no longer possible, and the excavation must be stopped to allow for an appropriate level of reinforcement. The results of this similitude experiment are consistent with those of the numerical calculations, thereby indicating that the results of both the experiments are valid. This also indicates that the results of the theoretical derivations are applicable for practical applications. Thus, the findings of this work will provide a useful source of scientific guidance for practical engineering works in the future.
(a)
(b)
(c)
(d)
(e)
(f)
5. Conclusions
In this study, the destabilization mechanisms and the patterns of overburden deformation and failure during underground resource mining were investigated via theoretical analyses, numerical simulations, and model experiments, based on the excavation of the 13210 working face in the Yima coal mine of the Gengcun village. The following conclusions were drawn from this study:(1)The degradation of overburden mechanical properties due to mininginduced disturbances is an important factor for the occurrence of rock burst failures in the overlying strata. If the geometric parameters of the working face are fully determined, a stiffness ratio no greater than 1 is required for the occurrence of rock burst failures in the overburden.(2)During the excavation of a working face, the stress values of the advanced abutments decrease with a decrease in distance from the tunnel face. As the distance from the tunnel face increases, the stresses acting upon the tunnel face’s advanced abutments initially increase before decreasing and stabilizing at some value; this corresponds to the existence of the “stress buildup zone,” “stress reduction zone,” and “native rock stress zone” in the overburden.(3)During the mining processes, the overburden exhibits a “pseudoplastic beam” stress state upon the complete release of local stresses in the overburden. Further, a “trapezoidal” fracture is formed by the cracks at stress concentrations. The continued progression of excavation expands the area of collapse in the overlying strata and increases the areas of fracture and compacted zones. More stress will then be transferred toward both ends of the working face to form highstress zones in these areas, which ultimately leads to the formation of a “threezone” distribution.(4)The progression of the working face excavation is accompanied by upward expansions in the overburden’s area of collapse and propagation of upwardgrowing bedseparating cracks; these are the key factors for the occurrence of ground settlement. Further, the nonuniform settlement of the overburden induces the generation of secondary fractures, which in turn expands the area of collapse of the overburden toward either end of the working face; this is an important factor that causes the exacerbation of surface settlement.(5)The results of our theoretical analyses, numerical calculations, and model experiments were completely consistent with each other. Thus, the findings of this work will be a useful source of scientific guidance for practical engineering works in the future.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Special Funds for Central Guiding Local Science and Technology Development (Project no. Z161100004516014), National Natural Science Foundation of China (Project nos. 51504029 and 51774048), and Beijing Nova Program (Project no. Z161100004916083).
References
 J. Tang, J. A. Wang, and L. Wang, “Dynamic laws and characteristics of surface movement induced by mining under thin alluvium,” Rock and Soil Mechanics, vol. 35, no. 10, pp. 2958–2968, 2014. View at: Google Scholar
 M. G. Qian, P. W. Shi, and J. L. Xv, Mining Pressure and Strata Control, China University of Mining and Technology Press, Xuzhou, China, 2010.
 M. G. Qian, X. X. Miao, and J. L. Xu, “Theoretical study of key stratum in ground control,” Journal of China Coal Society, vol. 21, no. 3, pp. 225–230, 1996. View at: Google Scholar
 X. X. Miao, R. H. Chen, and H. Pu, “Analysis of breakage and collapse of thick key strata around coal face,” Chinese Journal of Rock Mechanics and Engineering, vol. 24, no. 8, pp. 1289–1295, 2005. View at: Google Scholar
 L. X. Wu, J. Z. Wang, and Y. A. Liu, Coal Strip Mining under Buildings: Theory and Practice, China University of Mining and Technology Press, Xuzhou, China, 1994.
 P. A. Cundall, “A computer model for simulating progressive largescale movements in blocky rock systems,” in Proceedings of the ISRM Symposium on Rock Fracture, pp. 11–18, Nancy, France, 1971. View at: Google Scholar
 K. M. O’Connor and C. H. Dowding, “Distinct element modeling and analysis of mining induced subsidence,” Rock Mechanics and Rock Engineering, vol. 25, no. 1, pp. 1–24, 1992. View at: Google Scholar
 S. H. Li, M. H. Zhao, Y. N. Wang, and Y. Rao, “A new computational model of threedimensional DEMblock and particle model,” International Journal of Rock Mechanics and Mining Sciences, vol. 41, no. 3, p. 436, 2004. View at: Publisher Site  Google Scholar
 F. Li, J. A. Wang, P. F. Li, and K. Huang, “Research on movement behavior and failure mechanism of overlying strata caused by mining at mountainous field,” Rock and Soil Mechanics, vol. 37, no. 4, pp. 1089–1095, 2016. View at: Google Scholar
 L. M. Dou, J. He, and A. Cao, “Rock burst prevention methods based on theory of dynamic and static combined load induced in coal mine,” Journal of China Coal Society, vol. 40, no. 7, pp. 1469–1476, 2015. View at: Google Scholar
 X. B. Li, S. M. Wang, L. Weng, L. Q. Huang, T. Zhou, and J. Zhou, “Damage constitutive model of different age concretes under impact load,” Journal of Central South University, vol. 22, no. 2, pp. 693–700, 2015. View at: Publisher Site  Google Scholar
 Z. Q. Yin, X. B. Li, J. F. Jin, X. Q. He, and K. Du, “Failure characteristics of high stress rock induced by impact disturbance under confining pressure unloading,” Transactions of Nonferrous Metals Society of China, vol. 22, no. 1, pp. 175–184, 2012. View at: Publisher Site  Google Scholar
 M. F. Cai, J. A. Wang, and S. H. Wang, “Analysis on energy distribution and prediction of rock burst during deep mining excavation in Linglong gold mine,” Chinese Journal of Rock Mechanics and Engineering, vol. 20, no. 1, pp. 38–42, 2001. View at: Google Scholar
 S. C. Li, W. S. Zhu, and W. Z. Chen, “Application of elastoplastic large displacement finite element method to the study of deformation prediction of soft rock tunnel,” Chinese Journal of Rock Mechanics and Engineering, vol. 21, no. 4, pp. 466–470, 2002. View at: Google Scholar
 H. W. Song, Y. X. Jia, and Y. Y. Duan, “Study on characteristics of rock broken by excavation and roof supporting object,” Journal of China University of Mining and Technology, vol. 35, no. 2, pp. 192–196, 2006. View at: Google Scholar
 F. Q. Gong and G. F. Zhao, “Dynamic indirect tensile strength of sandstone under different loading rates,” Rock Mechanics and Rock Engineering, vol. 47, no. 6, pp. 2271–2278, 2014. View at: Publisher Site  Google Scholar
 Y. P. Wu, Y. G. Yang, and P. S. Xie, “Application of cusp catastrophe theory to research of mine water inrush law from roof,” Coal Science and Technology, vol. 35, no. 3, pp. 37–40, 2007. View at: Google Scholar
 C. Feng and S. H. Li, “2D particle contactbased meshfree method in CDEM and its application in geotechnical problems,” Engineering Computations, vol. 32, no. 4, pp. 1080–1103, 2015. View at: Publisher Site  Google Scholar
 Z. S. Ma, C. Feng, T. P. Liu, and S. H. Li, “A GPU accelerated continuousbased discrete element method for elastodynamics analysis,” Advanced Materials Research, vol. 320, pp. 329–34, 2011. View at: Publisher Site  Google Scholar
 H. C. Hu, Variational Principle of Elasticity and Its Application, Science Press, Beijing, China, 1981.
 W. C. Qian, “Variational principles in elasticity with nonliner stressstrain relation,” Applied Mathematics and Mechanics, vol. 8, no. 7, pp. 567–577, 1987. View at: Google Scholar
 F. Qiu and H. Ding, “Finite element method simulating failure of rock material,” Chinese Journal of Rock Mechanics and Engineering, vol. 26, no. S1, pp. 2663–2668, 2006. View at: Google Scholar
 R. Q. Huang, “Geostress distribution and unloading fracturing mechanism of high rock slopes in western part of China,” Journal of Engineering Geology, vol. 12, pp. 7–13, 2004. View at: Google Scholar
 F. Q. Wu, T. Liu, X. L. Tang et al., “Research on unloading and zonation of rock mass dam foundation excavationa case study of Xiaowan hydropower station,” Chinese Journal of Rock Mechanics and Engineering, vol. 28, no. 6, pp. 1091–1098, 2009. View at: Google Scholar
 D. Huang, Q. Tan, and R. Q. Huang, “Study of micromesoscopic characteristics of marble fracture surface and correlation with unloading rock mass strength under high stress and unloading,” Rock and Soil Mechanics, vol. 33, no. S2, pp. 7–15, 2012. View at: Google Scholar
 X. P. Zhou, Q. H. Qian, and B. H. Zhang, “Zonal disintegration mechanism of deep crackweakened rock masses under dynamic unloading,” Acta Mechanica Solida Sinica, vol. 22, no. 3, pp. 240–250, 2009. View at: Publisher Site  Google Scholar
 H. Q. Xie and C. H. He, “Study of the unloading characteristics of a rock mass using the triaxial test and damage mechanics,” International Journal of Rock Mechanics and Mining Sciences, vol. 41, no. S1, pp. 74–80, 2004. View at: Publisher Site  Google Scholar
 P. L. Gong, Y. Q. Hu, Y. S. Zhao, and D. Yang, “Threedimensional simulation study on law of deformation and breakage of coal floor on mining above aquifer,” Chinese Journal of Rock Mechanics and Engineering, vol. 24, no. 23, pp. 4396–4402, 2005. View at: Google Scholar
 J. W. Wu, H. S. Tong, S. J. Tong, and D. Q. Tang, “Study on similar material for simulation of mining effect of rock mass at fault zone,” Chinese Journal of Rock Mechanics and Engineering, vol. 26, no. S2, pp. 4170–4175, 2007. View at: Google Scholar
 W. D. Song, J. X. Fu, and D. X. Wang, “Study on physical and numerical simulation of failure laws of wall rock due to transformation from openpit to underground mining,” Journal of China Coal Society, vol. 37, no. 2, pp. 186–191, 2012. View at: Google Scholar
 W. M. Cheng, L. L. Sun, G. Wang, and X. C. Huang, “Similar material simulation test of steepinclined extrathick coal seam,” Journal of Mining and Safety Engineering, vol. 33, no. 3, pp. 387–392, 2016. View at: Google Scholar
 J. Xiao, “Selection of similar materials for model test and research on similar material proportioning test,” Beijing Jiaotong University, Beijing, China, 2013, M.S. thesis. View at: Google Scholar
 B. OrleckaSikora, S. Lasocki, G. Lizurek, and L. Rudziński, “Response of seismic activity in mines to the stress changes due to mining induced strong seismic events,” International Journal of Rock Mechanics and Mining Sciences, vol. 53, pp. 151–158, 2012. View at: Publisher Site  Google Scholar
 A. Sainoki and H. S. Mitri, “Dynamic behavior of mininginduced fault slip,” International Journal of Rock Mechanics and Mining Sciences, vol. 66, pp. 19–29, 2014. View at: Publisher Site  Google Scholar
 D. P. Ma, T. X. Wang, and Y. Liu, “An analysis of “spacetime” relationship of gobside entry driving in dynamic pressure area,” Journal of Mining and Safety Engineering, vol. 32, no. 3, pp. 465–470, 2015. View at: Google Scholar
 CUMT, “Study on stress transmission laws of mining floor and its influence on stability of floor roadway,” China University of Mining and Technology, Xuzhou, China, 2011, M.S. thesis. View at: Google Scholar
 H. Wang, Y. Jiang, S. Xue et al., “Influence of fault slip on mininginduced pressure and optimization of roadway support design in faultinfluenced zone,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 8, no. 5, pp. 660–671, 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Xiangfeng Lv 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.