Research Article  Open Access
Hengfeng Liu, Jixiong Zhang, Nan Zhou, Qiang Sun, Meng Li, Zhizhong Cui, "Experimental Investigation into the Bedding Plane Slip Effect on the Overlying Strata Behavior in Longwall Top Coal Caving of Soft Coal Seam", Advances in Civil Engineering, vol. 2019, Article ID 1718751, 14 pages, 2019. https://doi.org/10.1155/2019/1718751
Experimental Investigation into the Bedding Plane Slip Effect on the Overlying Strata Behavior in Longwall Top Coal Caving of Soft Coal Seam
Abstract
Bedding plane shear slip becomes more obvious in rock strata with quite different mechanical properties. A composite beam model considering the behaviors of the main roof and the immediate roof in the “threesoft” coal seam is established based on physical similarity simulation experiments and the slip theory in this paper. The movement and failure of the overlying strata and the mechanism of the bedding plane slip at the 2211 working face are studied by experimental and theoretical analysis. The results suggest that the front abutment stress distribution occurs 50 m ahead of the working face, the initial caving interval of the main roof is 55 m, and the peak stress appears at a distance 20 to 32 m ahead of the working face. The bedding plane slip areas can be divided into the obvious slip area and the slight slip area along the mining direction. The range of the obvious slip area becomes wider and the range of the slight slip area grows to be narrower as the working face advances. The bedding plane slip becomes steady after gradual increase and leads to the subsidence of the overlying strata in the “threesoft” coal seam. The observed initial caving interval of the main roof by field measurement is 51 m, which is consistent with the results of physical similarity simulation experiments and theoretical analysis. The results demonstrate that the beam slip model proposed in this paper is reasonable and able to describe the behaviors of overlying strata and bedding plane slip.
1. Introduction
Complex problems related to soft rocks are often encountered in engineering projects. As a worldwide issue from the 1960s, the soft rock problem has great influence on coal mine construction and mining safety [1], such as resulting in great mining difficulty, high technical cost, and low production efficiency. At present, it is of great significance for sustainable development, and there are many “threesoft” coal seams in China. The “threesoft” coal seam (TSCS) in this paper refers to a coal seam with soft roof, soft coal, and soft floor in coal mining [2]. The roof of TSCS is extremely fragile, and the coal walls often break during the mining process; that is, roof fall, rib spalling, and constraint of support structures occur from time to time. This not only affects the mining safety and efficiency, but threatens the life safety of mine workers. Thus, it is necessary to study the overburden movement and failure during the TSCS mining. A significant amount of research has been carried out on soft rocks, top coal caving, and TSCS mining.
Tsesarsky et al. [3] addressed the stability analysis in two different scales (localdiscontinuous and globalcontinuous) by developing a realistic approach for a threedimensional finite element analysis on the shallow underground openings excavated in soft and discontinuous rock masses. Yadav et al. [4] studied on the deformation evolution of model soft rocks during orthogonal cutting through highresolution imaging of orthogonal cutting experiments under 2D planestrain conditions. Bizjak and Petkovšek [5] suggested that the settlement of vault is most prominent after the excavation of a soft rock tunnel and timely construction of a reverse arch and closed lining can effectively reduce the displacement of the tunnel. Corthésy et al. [6] developed a modified doorstopper technique and affirmed its advantages in soft rock engineering. Shin et al. [7] proposed a preliminary guideline for evaluating the protection zone of the blast vibration and utilized a numerical method for dynamic modeling of a tunnel in soft rocks. Yang et al. [8] analyzed the mechanism of large deformation and control technology of deep soft rock roadways and put forward an improved boltcablereticulated shell supporting mode. Kikumoto et al. [9] established a constitutive model for soft rocks incorporating structural healing and decay in the context of the extended critical state theory. Singh et al. [10] pointed out that large deformation occurs on the premise of weak surrounding rocks under high in situ stress. Qin et al. [11] considered that the support for deep soft rock roadways under dynamic pressure is a major technical challenge, and a support technique was proposed based on the control effect of different support methods. He et al. [11, 12] proposed the technique of constant resistance and large deformation bolt support based on nonlinear large deformation theory of soft rocks. Wang and Pang [13] studied the adaptability evaluation and technical principle of fully mechanized top coal caving mining in extrathick coal seams. Vakili and Hebblewhite [14] introduced a new conceptual model for investigating top coal caving mechanism through analytical, observational and empirical methods. Unver and Yasitli [15] studied the top coal caving mechanism at the M3 longwall panel of Omerler Underground Mine using a finite difference code. Gao et al. [16] put forward an innovative numerical approach to simulate progressive caving of strata above a longwall face. Wu et al. [17] investigated the mechanism of coal wall spall and disastercausing mechanism in fully mechanized top coal caving face with a large dip angle. Wang et al. [18] studied the support technology of key parts of roadways with large dip angle in the soft coal seam. Xu et al. [19] performed an experimental research on the overburden movement and fissure evolution of fully mechanized caving face. Yao et al. [20] studied the applicability of large mining heightcombined sliding support in “threesoft” coal seams. These research studies promote the technique and theory development of TSCS mining. However, the strata movement law of TSCS mining still needs to be further studied [21] and many scholars ignore the influence of bedding plane slip on the overburden movement, especially in the case of huge differences between rock strata. In the process of mininginduced bending deformation, the shear slip among rock layers is more obvious due to incompatible deformation. Therefore, the slip effect is introduced in this work to study the overburden movement failure and bedding plane slip mechanism in the “threesoft” coal seams.
Based on the geological conditions at Linnancang Mine and laboratory tests, physical similarity simulation experiments are carried out and the overburden movement and bedding plane slip mechanism of fully mechanized top coal caving face in the TSCS are analyzed through mechanical analysis and similarity simulation experiments. The research results are expected to provide some reference for the safe mining of TSCSs.
2. Engineering Background
The 11# coal seam studied in this paper is in the No. 2 mining district at Linnancang Mine, China, which has a small change of dip angle and thickness variation. The roof and floor are mainly partially carbonaceous Claystone which are soft and fragile, and the coal body is soft as well. The immediate roof is composed of dark gray Claystone with a thickness of 2.1 m and the immediate floor consists of light gray Claystone with a thickness of 2.5 m. The 2211 working face is located northeast of the coal seam. The elevation of the ground is +5.30 m to +5.97 m, the elevation of the working face is −490.7 to −515.5 m, the strike length is 456.3 m, the slope length is 59.5 m, and the dip angle is 13°. Having a bulk density of 1.43 t/m^{3}, coal seam has an average thickness of 4.2 m and is mainly powdery durian that is soft and easy to fall. The detailed layout is shown in Figure 1.
There are many explanations about the concept of soft rocks, which can be divided into descriptive definition, indexing definition, and engineering definition. Among them, the International Society of Rock Mechanics (ISRM) offers an indexing definition: the soft rock is defined as the rock having uniaxial compressive strength between 0.25 MPa and 25 MPa [22]. Figure 2 shows the characteristics of the “threesoft” coal seam. The pressure behaviors at the working face of the “threesoft” coal seam include severe roof damage, serious rib spalling and roof collapse (gangue inflow), nonobvious first and periodic weighting, uneven supporting force, and low working resistance. The pressure features of the coal walls contain large pressure exerted on the support structures, huge deformation, obvious front abutment stress, which can cause severe rib spalling and floor heave. Hence, the mining pressure behaviors in the fully mechanized top coal caving face of the “threesoft” coal seam can be summarized as follows: (1) the first and periodic weighting is not obvious; (2) the resultant pressure point of the top coal moves forward, causing unreasonable performance of the rear support; (3) the deformed roof is easy to fall and inflow gangue; damage of the coal walls may lead to serious rib spalling; and low strength of the floor can easily induce floor heave. These features make the problem of roof fall more serious and result in more support work, which greatly increases mining and supporting difficulty, leading to heavy economic losses and casualties. Therefore, deformation and bedding plane slip mechanism are studied in this paper based on experimental and theoretical analysis.
3. Model Construction
3.1. Test Procedure
Before conducting the physical similarity simulation experiment, the geomechanical properties of coal and rocks are obtained based on sampling and laboratory tests. Then, the components of the simulation materials and the similarity ratio are determined according to the test results. The specific research process is illustrated as below (Figure 3).
3.2. Mechanical Property Test
The rock and coal samples collected from the mining area at Linnancang Mine were processed to numbered specimens. Physical and mechanical tests were carried out on these specimens. The test design and equipment are shown in Figure 4 [23] and Table 1.

3.3. Test Model Construction
In order to observe the evolution characteristics of overburden fissures and horizontal slip of bedding planes in the working face, the parameters of coal and rocks for the similarity experiment were determined based on the geological conditions of the 2211 working face and the mechanical property test results (Table 2).
 
Note. m_{1}, weight of sand; m_{2}, weight of calcium carbonate; m_{3}, weight of gyp; m_{4}, weight of water. 
A simulation model of 2.5 m × 2 m × 1.5 m was established based on the theory of similarity (Figure 5). The geometry similarity ratio is 1 : 150, and the height is 225 m. The model mining height was determined as 3.7 m based on field measurement and calculation. 6 sets of pressure boxes with an interval of 30 cm are used to monitor the stress change of the main roof during the mining process. The noncontact field strain measurement system (Vic2D system) was applied [24–26] to monitor the damage process of the roof and overlying strata, the fissure development, and the horizontal slip of bedding planes.
4. Test Results and Analysis
4.1. Mechanical Property Test Results
Through the physical and mechanical tests on the three kinds of specimens, the corresponding compressive, tensile and shear stressstrain curves are obtained (Figures 6–8).
As shown in Figures 6 and 7, the uniaxial compressive strengths of the three kinds of specimens, i.e., roof rocks, coal blocks, and floor rocks, are 22.39 MPa, 11.20 MPa, 22.35 MPa, respectively, and their tensile strengths are 2.03 MPa, 0.45 MPa, and 1.30 MPa, respectively. The normal stresses of the three specimens are positively correlated with their shear stresses (Figure 8). The cohesion (C) of the coal specimens and the rock specimens and their internal friction angles (φ) can be obtained based on the slopes and intercepts of the fitting curves. The test results show that the cohesion of the roof rock is 1.51 MPa, and the internal friction angle is 28°; the floor rock has cohesion of 0.65 MPa and an internal friction angle of 41°; the coal has cohesion of 0.34 MPa and its internal friction angle is 28°. The International Society of Rock Mechanics considers the rocks with the uniaxial compressive strength between 0.5 MPa and 25 MPa as soft rocks. Hence, it can be concluded that the coal, roof, and floor rocks studied in this work have relatively low strength and belong to a typical “threesoft” coal seam (Table 3).

4.2. Failure Mechanism of Overburden Movement
Mining operations can cause bending subsidence and fracture of overlying strata, and the fissure development is characterized into two kinds, namely, separation fissures and vertical fissures. With the advancing of the working face, the fissures will experience the development of initiation, propagation, and closure.
As shown in Figure 9, the overburden begins to undergo bending subsidence with the advancement of the 2211 working face. When the working face is advanced to 12.5 m, the first fracture occurs in the immediate roof, and the separation between the main roof and the upper rock layers appears. When the working face advances to 55 m, the initial fracture occurs in the main roof. Hence, it is considered that the initial caving interval of the main roof is 55 m. In addition to the main roof fractured, the overlying strata begin to subside, the separation grows larger, and the fissures continue to develop upward. The main roof breaks many times when the working face advances to 117 m, and there are many separations in the overburden. The vertical fissures penetrate upward and connect with each other and the phenomenon of bending subsidence becomes more apparent. As the working face advances to 210 m, bending subsidence develops to the upper part of the model, the separations gradually close, and the vertical fissures continue to propagate upward. When the process advances to 242 m, the upper surface of the model is severely damaged and the maximum fissure height develops up to 123 m. The overburden subsides completely, the separations reclose, and the vertical fissures have a full development when the working face is advanced to 315 m (to the mining stoppage line).
From Figures 10 and 11(a), it can be noted that the main roof begins to subside, and the subsidence amount increases gradually, and the subsidence peak value moves forward with the advancing of the working face. When the working face advances to 107 m, the maximum subsidence value of the main roof is 3.0 m. As the working face advances to 242 m, a large area of collapse occurs, and the maximum subsidence value is 3.37 m. When the working face is advanced to the stoppage line (315 m), the main roof having a maximum subsidence value of 3.5 m is no longer in a stable state with the passage of time.
(a)
(b)
The overburden subsidence is accompanied by the redistribution of stress. The distribution and evolution characteristics of the front abutment stress are obtained based on four groups of effective data on the vertical stress of the main roof collected in the similarity experiment, presented in Figure 11(b). As the working face gradually approaches the pressure boxes, the stress (in situ rock stress) at each measuring point shows the trend of being stable first, then increasing, and finally decreasing. In the range of 20 m to 32 m ahead of the working face, the stress reaches a peak value; then it slowly drops to 0 MPa, and the influence range of the front abutment stress is about 50 m. The front abutment stress of the measuring point 4 reaches a maximum value of 26.6 MPa, and its average stress concentration factor is 1.52. These indicate that the influence range of the front abutment stress of the 2211 working face is large and the position of the peak stress ahead the working face moves forward accordingly.
4.3. Analysis of Bedding Plane Slip Mechanism
According to the contour plots of the horizontal slip of bedding planes with different advance distances (Figure 12), bedding plane slip along the mining direction has obvious zoning area, which can be divided into the obvious slip area and the slight slip area. The obvious slip area corresponds to the deformation and failure area of overlying strata above the goaf. The obvious slip area shows a positivenegative layerline distribution in the vertical direction. Each slip layer line presents a positive and negative symmetry distribution with respect to the midpoint of the goaf along the horizontal direction. The slip mainly occurs in the indirect contact position of the rock layers with obvious horizontal fissures. The other area belongs to the slight slip area.
(a)
(b)
(c)
(d)
With the advancing of the working face, the obvious slip area increases gradually along the horizontal and vertical directions. The horizontal length of the slip is the same as that of the overburden bedding plane fissures. The vertical length of the slip is shorter than the height of the vertical fissures. According to field measurement, the horizontal fissures are shorter than the vertical fissures. There is no obvious bending subsidence in the strata, which means the bedding plane slip is not obvious. The overlying strata in the slight slip area in the working face have no obvious bending deformation, indicating that there is no obvious horizontal slip in this area. The horizontal slip value oscillates near zero. Therefore, the range of the obvious slip area mainly depends on the size of the mining affected areas.
Figure 13 gives the horizontal slip curves of overburden when the 2211 working face is advanced to 78 m, 117 m, 160 m, and 315 m, respectively. The slip curves in Figure 13 are obtained from the valuetaking curve in the physical similarity simulation experiment (Figure 12). Figure 13(a) is the bedding plane slip curve between the immediate roof and the main roof when the working face advanced to 78 m. It can be seen that the positive and negative symmetrical distribution of the bedding plane slip curves in the obvious slip zone is consistent with the results obtained from the contour plot. Meanwhile, the distribution of the slip curve of the obvious slip area is approximately an odd function relating to the center of the goaf. The maximum and minimum values of bedding plane slip are 0.06 m and −0.08 m, respectively, when the mining distance is 78 m, appearing at the positions of 50.3 m and 20.4 m ahead of the working face, respectively. The maximum and minimum values when the working face advances to 117 m, 160 m, and 315 m and their positions relative to the working face are given in Figure 14.
(a)
(b)
(c)
(d)
(a)
(b)
According to the contour plots, the maximum and minimum values of horizontal slip and their positions and distances away the working face can be obtained (Figure 14). It can be seen that the maximum (minimum) slip value in the obvious slip area gradually increases (decreases) and then tends to be stable, mainly because with the gradual advance of the working face, the overburden movement is fast at the beginning and then gradually becomes stable, that is, bedding plane slip in the overlying strata tends to be stable after the gradual increase. Besides, the positions of the maximum and minimum slip are moving forward, but the moving speed is lower than the advancing speed of the working face. Hence, the positions are gradually away from the work face.
5. Discussion
The mechanical models of the main roof with or without slip effect are established based on theoretical analysis. The effect of bedding plane slip on overburden movement and failure are obtained through comparative analysis. In the process of fully mechanized top coal caving mining, with the advancement of the working face, the top coal will pass through the top coal intact zone, the top coal failure and development zone, the fissure development zone, and the top coal caving zone until the coal is completely broken and released by the coal mouth [27]. The immediate roof (claystone) and the main roof (fine sandstone) of the fully mechanized caving face in the “threesoft” coal seam are combined as a doublelayer composite beam with two sides being fixed. The simplified slip mechanical model is established (Figure 15). Without considering the bedding plane slip, the mechanical model of the main roof can be directly established.
When the shear stress at the interface of the rock beam is higher than the allowable shear stress, bedding plane slip occurs in the doublelayer rock beam. Bending subsidence occurs in the doublelayer rock beam due to the external loading, and bending occurs in both the upper and lower strata. The differential equation of the horizontal slip between the rock layers is obtained [28, 29]:where is the slip displacement, h_{1} is the thickness of rock, and is the deflection equation. According to the structural characteristics of overlying strata, the firstorder expression of deflection of the rock beam satisfying boundary conditions is constructed by the Rayleigh method [30, 31]:where l is the advance distance of the working face. The function of slip displacement between the strata can be obtained according to equation (1):
The Rayleigh–Ritz method is used to solve the deflection equation, and the deformation potential energy of the doublelayer rock beam iswhere E_{s} and E_{x} are elastic moduli of the upper and lower strata of the doublelayer rock beam, respectively, ξ_{s} is the total strain of the upper stratum, and ξ_{x} is the total strain of the lower stratum, which can be expressed as follows:where S_{s} is the crosssectional area of the upper stratum and S_{x} is the crosssectional area of the lower stratum. The potential energy (V_{k}) of the doublelayer rock beam induced by the external loading of the overlying strata can be derived by
Slip potential energy produced by relative slip of upper and lower strata of doublelayer rock beam is obtained bywhere k_{s} is the equivalent average slip stiffness of the doublelayer rock beam. The total potential energy of the doublelayer rock beam (Π_{y} = W + V_{k} + U_{h}) can be achieved, namely,
In equation (8),where d_{1} and d_{2} are the distances from the neutral axes of the upper and lower strata to the neutral axes of the cross section of the doublelayer rock beam. The deformation potential energy of a thin plate is given in equation (10) if the slip effect is ignored:
Potential energy V generated by external loading is derived:
Total potential energy of the rock beam is Π_{n} and Π_{n} = W + V:
The total potential energy of the rock beam can be obtained as the following if the slip effect is ignored:
According to the principle of minimum potential energy, the variation of the total potential energy can be taken as zero [32]. The values can be obtained when the slip effect is considered or not. The stationary condition is
Based on the elastic theory, the relationship among deflection, bending moment, and shear force can be obtained. Then the bending moment equation of the doublelayer rock beam can be derived.
In order to verify the accuracy of the theoretical calculation, the theoretically calculated slip curve and the slip curve obtained by the similarity simulation experiment are compared. The parameters in the theoretical calculation are taken as follows: h_{1} is 10.6 m, h_{2} is 2.1 m, E_{1} is 12 GPa, E_{2} is 9.6 GPa, the overburden load q is 13.25 × 106 N/m, the equivalent average slip stiffness k_{s} is 3 MPa, and the ultimate bending moment of the upper rock beam is 3.879 × 107 N·m. These parameters are substituted above equations to obtain the theoretical calculation curves.
Figure 16 compares the theoretical calculation curves and the physical simulation curves. The maximum and average values in the simulation curve in Figure 16(a) (working face advanced to 32 m) are 0.096 m and 0.052 m, respectively. The maximum and average values of the subsidence curve calculated with considering the slip effect are 0.097 m and 0.048 m, respectively. The maximum and average values of the subsidence curve calculated without considering the slip effect are 0.078 m and 0.038 m, respectively. Therefore, the calculated subsidence curve of the main roof coincides well with the physical simulation results if the slip effect is considered while the calculated subsidence of the main roof is relatively small if the slip effect is ignored. It is concluded that the slip effect enhances the subsidence of the main roof. In addition, the average value and the maximum value of the fitting curve (R^{2} = 0.98) obtained by the physical similarity simulation are 0.13 m and 0.40 m, respectively, if the range of the working face is focused on 0 to 178.80 m (in the strike direction), and these values become −0.14 m and −0.34 m if the range is focused on 178.80 to 375 m. In the theoretical calculation, the average and maximum slip values are 0.16 m and 0.37 m, respectively, when the working face ranges from 0 to 186.80 m, while these values are −0.16 m and −0.37 m when the range is focused on 186.80 to 375 m. The fitting curves obtained by theoretical calculation are in good agreement with that derived from the physical simulation results. Hence, it is concluded that the theoretical calculation is reliable and is able to reflect the behaviors of the overlying strata and bedding plan slip.
(a)
(b)
The initial caving intervals of the main roof are theoretically calculated when the slip theory is considered or not (Figure 17). The calculated initial caving intervals are 52 m and 60 m, respectively. The initial caving intervals obtained by field measurement and physical simulation are 51 m and 55 m, respectively. The contour plot of the positive stress in the main roof beam is given in Figure 18. If the initial caving interval is assumed to be 52 m, three peak compressive stresses and three peck tensile stresses can be found in Figure 18. Peak tensile stresses appear at the top of the front part and the back part of the rock beam, as well as the bottom of the middle part. Two peak compressive stresses appear at the bottom of the front part and the back part of the rock beam with another locating at the top of the middle part. The value of the peak tensile stress at the bottom of the main roof beam is 2.03 MPa. The result suggests that the bottom of the middle part of the rock beam tends to fracture in the early stage.
6. Conclusion
The results obtained lead to the following conclusions:(1)Results of the similarity simulation experiment show that the front abutment stress distributes 50 m ahead of the 2211 working face at Linnancang Mine and the peak stress appears at a distance 20 to 32 m ahead of the working face. The initial caving interval of the main roof is 55 m, and the periodic weighting is not obvious. The bedding plane slip in the mining direction has obvious zoning phenomenon, and the slip area is divided into the obvious slip area and the slight slip area.(2)Based on the slip theory and the characteristics of overburden movement in the “threesoft” coal seam, a mechanical model of a doublelayer rock beam is established, and the calculation equations for the deflection of the beam, horizontal slip value, and bending moment with or without slip effect are deduced. The initial caving interval under the influence of slip is calculated to be 52 m, that is, 1 m different from the 51 m obtained by field measurement. The initial caving interval is 60 m without considering the slip effect, which is quite different from the result of field measurement (51 m). In addition, the normal stress distribution in the main roof is achieved, and the first breaking position is determined according to the calculation result.(3)Physical similarity simulation experiment, theoretical analysis, and field measurement are adopted to study the overburden movement and bedding plane slip mechanism in the “threesoft” coal seam, and the results by these different methods are compared. The comparison demonstrates that the doublelayer rock beam model proposed in this paper is reasonable and able to describe the behaviors of overlying strata and bedding plane slip, and the slip effect accelerates the overburden movement in the “threesoft” coal seam; hence, attention should be paid to the dangers caused by the slip in the “threesoft” coal seam mining.
Notations
:  Slip displacement 
h_{1}:  Thickness of rock 
:  Deflection equation 
l:  Advance distance of the working face 
E_{s}:  Elastic modulus of the upper strata of the doublelayer rock beam 
E_{x}:  Elastic modulus of the lower strata of the doublelayer rock beam 
ξ_{s}:  Total strain of the upper stratum 
ξ_{x}:  Total strain of the lower stratum 
S_{s}:  Crosssectional area of the upper stratum 
S_{x}:  Crosssectional area of the lower stratum 
V_{k}:  Potential energy 
k_{s}:  Equivalent average slip stiffness of the doublelayer rock beam 
Π_{y}:  Total potential energy of the doublelayer rock beam with slip effect 
d_{1}:  Distances from the neutral axes of the upper strata to the neutral axes of the cross section of the doublelayer rock beam 
d_{2}:  Distances from the neutral axes of the lower strata to the neutral axes of the cross section of the doublelayer rock beam 
Π_{n}:  Total potential energy of the doublelayer rock beam with no slip effect. 
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Disclosure
The authors are solely responsible for the content.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
This research was funded by the National Key R&D Program of China (grant no. 2018YFC0604704), the National Science Fund for Distinguished Young Scholars (grant no. 51725403), and the National Natural Science Foundation of China (grant no. 51874287).
References
 M. C. He, H. H. Jing, and X. M. Sun, Soft Rock Engineering Mechanics, Science Press, Beijing, China, 2002.
 B. T. Zhao and B. Q. Lin, Gas Emission and Prevention Technology in Three Soft and Unstable Coal Seam with Low Permeability, China University of Mining and Technology, Xuzhou, China, 2007.
 M. Tsesarsky, E. Gal, and E. Machlav, “3D globallocal finite element analysis of shallow underground caverns in soft sedimentary rock,” International Journal of Rock Mechanics and Mining Sciences, vol. 57, pp. 89–99, 2013. View at: Publisher Site  Google Scholar
 S. Yadav, C. Saldana, and T. G. Murthy, “Experimental investigations on deformation of soft rock during cutting,” International Journal of Rock Mechanics and Mining Sciences, vol. 105, pp. 123–132, 2018. View at: Publisher Site  Google Scholar
 K. F. Bizjak and B. Petkovšek, “Displacement analysis of tunnel support in soft rock around a shallow highway tunnel at Golovec,” Engineering Geology, vol. 75, no. 1, pp. 89–106, 2004. View at: Publisher Site  Google Scholar
 R. Corthésy, M. H. Leite, D. E. Gill, and B. Gaudin, “Stress measurements in soft rocks,” Engineering Geology, vol. 69, no. 34, pp. 381–397, 2003. View at: Publisher Site  Google Scholar
 J.H. Shin, H.G. Moon, and S.E. Chae, “Effect of blastinduced vibration on existing tunnels in soft rocks,” Tunnelling and Underground Space Technology, vol. 26, no. 1, pp. 51–61, 2011. View at: Publisher Site  Google Scholar
 X. Yang, E. Wang, Y. Wang, Y. Gao, and P. Wang, “A study of the large deformation mechanism and control techniques for deep soft rock roadways,” Sustainability, vol. 10, no. 4, pp. 1100–1120, 2018. View at: Publisher Site  Google Scholar
 M. Kikumoto, V. P. Q. Nguyen, H. Yasuhara, and K. Kishida, “Constitutive model for soft rocks considering structural healing and decay,” Computers and Geotechnics, vol. 91, pp. 93–103, 2017. View at: Publisher Site  Google Scholar
 B. Singh, R. K. Goel, J. L. Jethwa, and A. K. Dube, “Support pressure assessment in arched underground openings through poor rock masses,” Engineering Geology, vol. 48, no. 12, pp. 59–81, 1997. View at: Publisher Site  Google Scholar
 D. Qin, X. Wang, D. Zhang, and X. Chen, “Study on surrounding rockbearing structure and associated control mechanism of deep soft rock roadway under dynamic pressure,” Sustainability, vol. 11, no. 7, pp. 1892–1907, 2019. View at: Publisher Site  Google Scholar
 M. He, W. Gong, J. Wang et al., “Development of a novel energyabsorbing bolt with extraordinarily large elongation and constant resistance,” International Journal of Rock Mechanics and Mining Sciences, vol. 67, pp. 29–42, 2014. View at: Publisher Site  Google Scholar
 G. F. Wang and Y. H. Pang, “Fullmechanized coal mining and caving mining method evaluation and key technology for thick coal seam,” Journal of China Coal Society, vol. 43, no. 1, pp. 33–42, 2018. View at: Google Scholar
 A. Vakili and B. K. Hebblewhite, “A new cavability assessment criterion for longwall top coal caving,” International Journal of Rock Mechanics and Mining Sciences, vol. 47, no. 8, pp. 1317–1329, 2010. View at: Publisher Site  Google Scholar
 B. Unver and N. E. Yasitli, “Modelling of strata movement with a special reference to caving mechanism in thick seam coal mining,” International Journal of Coal Geology, vol. 66, no. 4, pp. 227–252, 2006. View at: Publisher Site  Google Scholar
 F. Gao, D. Stead, and J. Coggan, “Evaluation of coal longwall caving characteristics using an innovative UDEC trigon approach,” Computers and Geotechnics, vol. 55, pp. 448–460, 2014. View at: Publisher Site  Google Scholar
 Y. P. Wu, D. Lang, and P. S. Xie, “Mechanism of disaster due to rib spalling at fullymechanized top coal caving face in soft steeply dipping seam,” Journal of China Coal Society, vol. 41, no. 8, pp. 1878–1884, 2016. View at: Google Scholar
 X. F. Wang, Y. Wang, and D. S. Zhang, “Enhanced support technology for key area of the roadway in large inclined angle “threesoft” coal seam,” Journal of Mining and Safety Engineering, vol. 34, no. 2, pp. 208–213, 2017. View at: Google Scholar
 M. G. Xu, P. Wei, S. G. Li, and J.J. Wang, “Experimental study on overburden migration and fracture evolution law of “three soft” coal seam fully mechanized workingface,” Journal of China Coal Society, vol. 42, no. 1, pp. 122–127, 2017. View at: Google Scholar
 Q. L. Yao, S. G. Cao, F. H. Wang, and Z.G. Jiang, “Applicability of combined slip support of large mining height in threesoft coal seam of mediumthickness,” Journal of Mining and Safety Engineering, vol. 27, no. 2, pp. 185–189, 2010. View at: Google Scholar
 J. C. Wang, “Engineering practice and theoretical progress of topcoal caving mining technology in China,” Journal of China Coal Society, vol. 43, no. 1, pp. 43–51, 2018. View at: Google Scholar
 M. He, “Latest progress of soft rock mechanics and engineering in China,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 6, no. 3, pp. 165–179, 2014. View at: Publisher Site  Google Scholar
 H. Liu, J. Zhang, W. Zhang, F. Gao, H. Yan, and T. L. An, “Experimental investigation of perceptual characteristics of functional cemented backfilling materials in coal mines,” Minerals, vol. 9, no. 1, pp. 55–70, 2019. View at: Publisher Site  Google Scholar
 Q. Sun, J. Zhang, Q. Zhang, and X. Zhao, “Analysis and prevention of geoenvironmental hazards with highintensive coal mining: a case study in China’s western ecoenvironment frangible area,” Energies, vol. 10, no. 6, pp. 786–801, 2017. View at: Publisher Site  Google Scholar
 Y. Huang, J. Zhang, W. Yin, and Q. Sun, “Analysis of overlying strata movement and behaviors in caving and solid backfilling mixed coal mining,” Energies, vol. 10, no. 7, pp. 1057–1073, 2017. View at: Publisher Site  Google Scholar
 Q. Sun, J. Zhang, and N. Zhou, “Study and discussion of short strip coal pillar recovery with cemented paste backfill,” International Journal of Rock Mechanics and Mining Sciences, vol. 104, pp. 147–155, 2018. View at: Publisher Site  Google Scholar
 J. H. Wang, Z. Z. Huang, and L. Yu, “Three plies in top coal” theory and its application in top coal caving mining for ultrathick coal seams,” Journal of China Coal Society, vol. 42, no. 4, pp. 809–816, 2017. View at: Google Scholar
 H. T. Sun, “Study on the deformation fracture mechanism of the surface borehole with the excavation disturbance,” Chongqing University, Chongqing, China, 2008, Ph.D. thesis. View at: Google Scholar
 D. Y. Liu, H. T. Sun, and Y. Zhang, “A model of shear slipping of overlying strata under mining disturbance,” Rock and Soil Mechanics, vol. 31, no. 2, pp. 609–614, 2010. View at: Google Scholar
 Y. Pan, S. T. Gu, and Y. S. Qi, “Analytic solution of tight roof’s bending moment deflection and shear force under advanced supercharger load and supporting resistance before first weighting,” Chinese Journal of Rock Mechanics and Engineering, vol. 32, no. 8, pp. 1544–1553, 2013. View at: Google Scholar
 Z. L. Xu, Concise Tutorial on Elastic Mechanics, Advanced Education, Beijing, China, 2013.
 H. A. F. A. Santos and C. I. Almeida Paulo, “On a pure complementary energy principle and a forcebased finite element formulation for nonlinear elastic cables,” International Journal of NonLinear Mechanics, vol. 46, no. 2, pp. 395–406, 2011. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Hengfeng Liu 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.